Arxiv:2105.12820V1 [Astro-Ph.EP] 26 May 2021 Known with Exquisite Precision
Total Page:16
File Type:pdf, Size:1020Kb
Draft version May 28, 2021 Typeset using LATEX twocolumn style in AASTeX63 The First Dynamical Mass Measurement in the HR 8799 System 1, 1 2 3 G. Mirek Brandt , ∗ Timothy D. Brandt , Trent J. Dupuy , Daniel Michalik , and Gabriel-Dominique Marleau 4, 5, 6 1Department of Physics, University of California, Santa Barbara, Santa Barbara, CA 93106, USA 2Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh, EH9 3HJ, UK 3European Space Agency (ESA), European Space Research and Technology Centre (ESTEC), Keplerlaan 1, 2201 AZ Noordwijk, The Netherlands 4Institut f¨urAstronomie und Astrophysik, Universit¨atT¨ubingen,Auf der Morgenstelle 10, 72076 T¨ubingen,Germany 5Physikalisches Institut, Universit¨atBern, Gesellschaftsstr. 6, 3012 Bern, Switzerland 6Max-Planck-Institut f¨urAstronomie, K¨onigstuhl17, 69117 Heidelberg, Germany (Received April 7, 2021; Revised May 14, 2021; Accepted May 26, 2021) Submitted to ApJ Letters ABSTRACT HR 8799 hosts four directly imaged giant planets, but none has a mass measured from first principles. We present the first dynamical mass measurement in this planetary system, finding that the innermost +1:9 planet HR 8799 e has a mass of 9:6−1:8 MJup. This mass results from combining the well-characterized orbits of all four planets with a new astrometric acceleration detection (5σ) from the Gaia EDR3 version of the Hipparcos-Gaia Catalog of Accelerations. We find with 95% confidence that HR 8799 e +24 is below 13 MJup, the deuterium-fusing mass limit. We derive a hot-start cooling age of 42−16 Myr for HR 8799 e that agrees well with its hypothesized membership in the Columba association but is also consistent with an alternative suggested membership in the β Pictoris moving group. We exclude the presence of any additional &5-MJup planets interior to HR 8799 e with semi-major axes between 3{16 au. We provide proper motion anomalies and a matrix equation to solve for the mass of any of ≈ the planets of HR 8799 using only mass ratios between the planets. Keywords: | 1. INTRODUCTION ics, mass measurements for any of the planets in the HR 8799 is the only star that has been detected with HR 8799 system. four directly imaged exoplanets (Marois et al. 2008; The HR 8799 planets induce accelerations on their Marois et al. 2010; Currie et al. 2011). Each planet host star that are measurable by sufficiently precise ab- orbits counter-clockwise on the sky, in a gently eccen- solute astrometry, enabling direct, dynamical measure- tric orbit, with a nearly face-on configuration (Sudol & ments of the planets' masses. The Hipparcos-Gaia cat- Haghighipour 2012; Pueyo et al. 2015; Konopacky et al. alog of accelerations (HGCA, Brandt 2018, 2021) has 2016; Currie 2016). The orbits of all four planets are cross-calibrated the positions and proper motions so that these accelerations can be used for inference. The arXiv:2105.12820v1 [astro-ph.EP] 26 May 2021 known with exquisite precision. Wang et al.(2018) stud- ied the system in detail and restricted orbital parameter acceleration of HR 8799 A was not significant in the space based on configurations that are dynamically sta- Gaia DR2 (Lindegren et al. 2018) version of the cata- ble. There have been many mass constraints based on log (Brandt 2018). However, Gaia's precision on bright the long term dynamical stability of the system (e.g., stars has increased significantly in EDR3 (Gaia Collab- oration et al. 2020; Lindegren et al. 2020; Brandt 2021). all planets . 5MJup by Esposito et al. 2013; inner three In the Gaia EDR3 version of the HGCA, HR 8799 A planets < 13MJup by Pueyo et al. 2015). However, there are no robust dynamical, i.e., from Newtonian mechan- is accelerating at nearly 5σ. In this paper we infer the mass, and thereafter the age, of HR 8799 e from this acceleration. ∗ NSF Graduate Research Fellow 2 Brandt et al. The HR 8799 system is most commonly thought Wang et al.(2018) used the Gaia DR1 parallax of to be an 40-Myr-old member of the Columba as- 24:76 0:64 mas (Gaia Collaboration et al. 2016)1. ≈ ± sociation (Zuckerman et al. 2011), but Lee & Song Gaia EDR3 measures a much more precise value of (2019) recently suggested membership with the younger 24:462 0:046 mas (Lindegren et al. 2020). We scale ± β Pic Moving Group (BPMG). Given a 40-Myr age, all 1000 MCMC samples to the Gaia EDR3 parallax, ≈ hot-start evolutionary models predict masses of 7 removing its contribution to the uncertainty of the or- ± 2 MJup for the innermost three planets, well below the bital fits. Defining r to be the ratio of a given chain's deuterium-fusion mass boundary ( 12-13 M ; Spiegel parallax to the Gaia EDR3 value, we multiply the semi- ≈ Jup et al. 2011). Lower masses (.5{7 MJup) improve the major axes by r (to preserve the relative astrometry) system's dynamical stability (Fabrycky & Murray-Clay and multiply the system mass by r3 (to preserve the 2010; Sudol & Haghighipour 2012), but resonant locking orbital periods). We keep all other orbital parameters could render the system stable at higher planet masses unchanged. Updating the parallax to its Gaia EDR3 (G¨otberg et al. 2016; Wang et al. 2018; Go´zdziewski& value ultimately improves the fractional error on the fi- Migaszewski 2018, 2020). Higher masses would imply nal mass estimate from 20% to 19.5%. either an older age for the system or entropy loss during We use the open-source tool htof (Brandt et al. formation (a colder start), although there is a limit to 2020, submitted; Brandt & Michalik 2020; Brandt et al. the entropy that can be lost at formation (Marleau & 2021b) to model Hipparcos and Gaia observations. In Cumming 2014). Conversely, lower masses would sug- brief, htof uses the Hipparcos intermediate astromet- gest a younger age and therefore favor membership with ric data (both ESA(1997) and van Leeuwen(2007) re- the BPMG. Dynamical mass measurements can conclu- ductions) and predicted scan angles and observational sively test such hypotheses. The inferred atmospheric epochs of Gaia (via GOST2), to generate synthetic Hip- properties and chemical abundances of the planets also parcos and Gaia astrometry for any orbit. depend on their assumed masses (e.g., Wang et al. 2020; The astrometric measurement that we use is the Molli`ereet al. 2020). proper motion anomaly. This is the difference between We detail our method in Section2. In Section3 we a nearly instantaneous proper motion from Gaia EDR3 consider planet mass ratios estimated from their relative and a long-term proper motion (the difference in posi- luminosities, show that the mass of HR 8799 e is insensi- tion between the Hipparcos and Gaia astrometry mis- tive to these mass ratios, and obtain a robust dynamical sions divided by the time between the measurements). mass measurement. We derive a substellar cooling age We denote these proper motion anomalies as, e.g., for HR 8799 e in Section4 and conclude in Section5. α Gaia α Hip ∆µα∗ = µGaia ∗ − ∗ (1) − tGaia tHip − where α = α cos δ. Parameters with subscript Hip ∗ refers to the average of those parameters from Hippar- cos 2. METHODOLOGY 2007 and 1997, weighted 60/40 as they are in the HGCA (Brandt 2018). A set of orbital parameters then The acceleration that HR 8799 A experiences is the gives predicted values for ∆µα∗ and ∆µδ as functions of sum of the acceleration due to each of its four planetary the masses of the planets, mb; mc; md, and me. Because companions. Because the planet masses are small com- HR 8799 A's motion closely follows a linear combina- +0:11 pared to the mass of HR 8799 A (1:47−0:08 M ; Wang tion of Keplerian orbits, we can represent its predicted et al. 2018), the star's motion is approximately given proper motion anomaly as the Jacobian by a linear combination of the orbits of the planets, 2 3 weighted by their masses. We can then optimize the mb " # " @∆µα∗ @∆µα∗ @∆µα∗ @∆µα∗ # planet masses until the modelled host star acceleration ∆µ ∗ 6m 7 α = @mb @mc @md @me 6 c7 : matches the observed value. @∆µδ @∆µδ @∆µδ @∆µδ 6 7 ∆µδ 4md5 We use all 1000 samples from the orbital posteriors model @mb @mc @md @me m from Wang et al.(2018), published on the github repos- e (2) itory for the resource whereistheplanet.com(Wang et al. 2021). The orbits correspond to the coplanar, dynami- cally stable case of the HR 8799 system (see Table 4 of 1 +u In all cases where we quote a posterior by listing m l or m ± σ: Wang et al. 2018). Each set of orbital parameters, to- m denotes the median with l and u (or singularly,−σ) denoting gether with the masses of the four planets and the star, the 16% and 84% confidence intervals, respectively. 2 predicts the motion of HR 8799 A. https://gaia.esac.esa.int/gost/ The First Dynamical Mass Measurement in the HR 8799 System 3 We compute the partial derivatives by using only the We compute the partial derivatives of Equation (3) orbit of a given planet and assigning that planet unit against ∆µα∗, and ∆µδ. Setting these partials to zero mass. We do this for all 1000 orbital draws. We use gives an under-constrained system of equations. The REBOUND and the ias15 scheme (Rein & Liu 2012; Rein two components of the proper motion anomaly measure & Spiegel 2015) to integrate the orbits in time. a combination of the masses of the four planets, but not When we compute the model partial derivatives, we the four masses individually.