Letters https://doi.org/10.1038/s41550-018-0561-6

The of the young planet b through the astrometric motion of its host

I. A. G. Snellen* and A. G. A. Brown

The young massive discovered with high-contrast Beta Pic is included in the second data release16 (DR2) of Gaia imaging1–4 provide a unique opportunity to study the formation and covers 32 observations for this star, between 1 October 2014 and early evolution of gas giant planets. A key question is to and 19 April 2016, of which 30 were used for the astrometric solu- what extent gravitational energy from accreted gas contributes tion. The individual measurements are not provided in DR2. For to the internal energy of a newly formed planet. This has led to most observed by Hipparcos, Gaia is expected to eventually a range of formation scenarios from ‘cold’ to ‘hot’ start mod- achieve a precision that is 100 times better17. However, β​ Pic is els5–8. For a planet of a given mass, these initial conditions gov- such a bright star (apparent visual magnitude, V =​ 3.86) that its ern its subsequent evolution in luminosity and radius. Except images in the Gaia data are saturated, substantially degrading the for upper limits from radial velocity studies9,10, disk modelling11 achieved positional accuracy, resulting in DR2 uncertainties of and dynamical instability arguments12, no mass measurements 0.3 mas in each positional direction and parallax, and 0.7 mas yr−1 of young planets are yet available to distinguish between in proper motion. However, since the Hipparcos and Gaia mis- these different models. Here, we report on the detection of sions span a baseline of 24 yr, the combined positional data pro- the astrometric motion of Beta Pictoris, the ~21-Myr-old host vide a long-term proper-motion measurement16 with a precision star of an archetypical directly imaged gas giant planet, around of 0.02 mas yr−1 (Table 1). the system’s centre of mass. Subtracting the highly accurate Figure 2 shows the heliocentric Hipparcos positional data after Hipparcos13,14 and Gaia15,16 proper motion from the internal 3 yr subtraction of the 24 yr Hipparcos–Gaia baseline. The residual Hipparcos astrometric data reveals the reflex motion of the proper motion is governed by the reflex motion of β​ Pic around star, giving a model-independent planet mass of 11 ±​ 2 the centre of mass of the system. The planet β ​Pic b is in a nearly . This is consistent with scenarios in which the planet is edge-on orbit at a position angle of 212° (ref. 18), consistent with formed in a high-entropy state as assumed by hot start mod- the observed motion of the star in the Hipparcos data. The planet’s els. The ongoing data collection by Gaia will soon lead to mass discovery was announced in 2009 from high-contrast images taken measurements of other young gas giants and form a great in 20032,3. Only after the planet reappeared on the other side (south- asset to further constrain early-evolution scenarios. west) of the star, it started to be intensely monitored, in particu- The European Space Agency’s Hipparcos13 (1989–1993) and lar with the new generation of high-contrast imaging instruments Gaia15 (2013–present) astrometric space observatories use the same at the and SPHERE measurement principles. Both contain two telescopes with lines (Spectro-Polarimetric High-contrast REsearch) on the of sight perpendicular to the spin axis of a continuously rotating Very Large Telescope18–20. It means that while about half of the spacecraft, which slowly precesses while maintaining a constant orbit is well known, the is still poorly constrained to angle with respect to the spacecraft–Sun direction. The telescopes be between 20.2 yr and 26.3 yr (ref. 18). Since the reference epochs precisely time the crossing of stars in their field of view, providing of Hipparcos (1991.25) and Gaia DR2 (2015.5) are separated by accurate one-dimensional astrometric positions, which result in a 24.25 yr, the observations of the former mission fall in the well- rigid system of reference by combining the relative positions of the determined part of the planet’s orbit, but the precise orbital phase observed stars. is uncertain. Depending on the orbital period, the Hipparcos–Gaia Hipparcos observed Beta Pictoris (β​ Pic) 111 times between baseline also contains a small part of the stellar reflex motion. For 1990.0 and 1993.1 (Julian epochs are used throughout this Letter), example, for an orbital period of 22 yr and a planet mass of 10 distributed over 35 spacecraft orbits. Each set of ~3–5 observations Jupiter masses (MJ), the two reference epochs correspond to orbital taken during one orbit has similar telescope scanning directions. phases that differ by 0.1, resulting in a planet-induced position dif- While the initial data release delivered one position per orbit13, ference of 1.0 mas and an effect on the Hipparcos–Gaia baseline of with typical uncertainties in the range of 1.5–2.5 mas, the subse- 1.0 mas ÷​ 24.25 yr =​ 0.04 mas yr−1. This effect is taken into account in quent reprocessing of these data by van Leeuwen14 maintained the our calculations below. individual data points and delivered uncertainties in the range of Since the projected orientation of the orbit is well constrained18, 0.8–1.0 mas (see Supplementary Table 1), improving the precision so is the expected direction of stellar reflex motion. Hence, we con- by a factor of ~4. The latter analysis provides for β​ Pic a parallax of verted the positional data from Hipparcos into one-dimensional 51.44 ±​ 0.12 mas, corresponding to a distance of 19.44 ±​ 0.05 pc, and position measurements in this direction as a function of time, a proper motion of 4.65 ±​ 0.11 mas yr−1 and 83.10 ±​ 0.15 mas yr−1 shown in the right panel of Fig. 2. The solid curves show the best- in the right ascension (RA) and the declination (dec.) direction, fit stellar motion for several trial orbital periods, indicating that respectively (Fig. 1; Table 1). The star is the largest member of an short periods are not consistent with the Hipparcos data since these association of young stars, the β​ Pic moving group (BPMG), sharing would result in substantial acceleration in the stellar proper motion the same origin and motion through space. in the period 1990–1993, which is not observed. Figure 3 shows the

Leiden Observatory, Leiden University, Leiden, Netherlands. *e-mail: [email protected]

Nature Astronomy | www.nature.com/natureastronomy Letters Nature Astronomy

200 β Pic reflex motion 2 4

1 1990.3 2 24 yr (mas) 23 yr 1990.9 HG 22 yr 0 1991.4 0 1991.9 21 yr Δ Pos (mas) 1992.6 100 –1 –2 Δ dec. – PM

–2 –4 2 1 0 –1 –2 1989 1990 1991 1992 1993

ΔRA – PMHG (mas) of observation

Fig. 2 | The heliocentric Hipparcos positional data after subtraction of the 24 yr Hipparcos–Gaia baseline. Left: the change in position of β​ Pic in 0 the Hipparcos data relative to the 24 yr Hipparcos–Gaia baseline, where

PMHG is the proper motion as derived from the Hipparcos–Gaia baseline. Δ dec. (mas) The one-dimensional positions taken during 35 Hipparcos spacecraft orbits were combined into five two-dimensional astrometric positions. The shaded areas indicate the 1σ​ uncertainty areas (s.e.m.) with their mean epochs. The arrow shows the expected direction for the movement of the star around the centre of mass of the system. Right: the change in position,

–100 Δ​Pos, of the star relative to the Hipparcos–Gaia baseline in the direction of the expected reflex motion direction of the star at 212°, known from the orbital mapping of the planet. Positions averaged over a spacecraft orbit are shown, which strongly vary in accuracy (grey scaling, with the darkest greys showing the smallest errors) depending on the relative orientation of the one-dimensional astrometric measurement. Five measurements are omitted because their measurement vectors are nearly perpendicular to the reflex motion movement. The coloured lines show the best-fit velocity –200 150 100 50 0 –50 curves for orbital periods of 21–24 yr. Planet orbital periods shorter than ΔRA (mas) 22 yr are excluded since they would show substantially more curvature than the Hipparcos data allow. Fig. 1 | The astrometric motion of the young star β Pic as shown by the Hipparcos data (1990–1993). The black curve is the best-fit model to the parallax (caused by the ’s yearly motion around the Sun) constraints to the planet’s mass and orbital period. The green shaded area indicates the 1σ​ uncertainty interval for the orbital period from and the star’s proper motion, with the points indicating the modelled 18 positions for the one-dimensional Hipparcos positions (the dark blue the high-contrast imaging monitoring , and the contours indicate lines are perpendicular to the scanning direction), which are averaged the 1σ​, 2σ​ and 3σ ​limits on orbital period and mass from the study presented here. The planet’s mass is constrained to 11 ±​ 2 MJ (1σ​), per satellite orbit for clarity and have typical error bars of <​1 mas (s.e.m.). The grey circle and line indicate the separate parallax and proper-motion and the orbital period is likely to be >​22.2 yr (2σ ​limit). Previous upper limits to the planet’s mass of <​20–30 MJ from radial veloc- components, respectively. The light blue line shows the proper motion as 9,10 determined from the 24 yr Hipparcos–Gaia baseline, with the red arrow ity monitoring are consistent with the astrometric measurement showing the difference from the internal Hipparcos proper motion caused presented here. by the reflex motion of the star around the system’s barycentre. This We compared our measurement with the mass estimates obtained difference is amplified by a factor of ten to make the effect visible. from gas giant cooling models. The wide range of estimates in the lit- erature are driven by different assumptions on the age of the planet and the internal energy of the planet at time of formation. In partic- ular, it is uncertain what fraction of the released gravitational energy Table 1 | The Hipparcos and Gaia DR2 astrometric data on β Pic of newly accreted material during formation is lost in the accretion (HIP 27321; Gaia DR2 4792774797545105664) shock on the planet’s surface or circumplanetary disk and not incor- porated in the planetary structure, resulting in a range of hot start, Hipparcos (International Celestial Gaia DR2 (International Celestial cold start and intermediate ‘warm’ start models5–8. In addition, dif- Reference System; epoch Reference System; epoch 2015.5) ferent age-determination methods give different stellar ages, and it 1991.25) is also not certain whether the planet was formed at exactly the same RA =​ 86.82118054° ±​ 0.10 mas RA =​ 86.82123366° ±​ 0.314 mas time as the star or possibly a few million later. Isochronal fits dec. =​ −​51.06671341° ±​ 0.11 mas dec. =​ −​51.06614803° ±​ 0.342 mas to low-mass BPMG members provide21 an age of 8–20 Myr, while plx =​ 51.44 ±​ 0.12 mas plx =​ 50.62 ±​ 0.33 mas dynamical ages of the BPMG provide age estimates from 10–12 Myr 22 23 1 1 (ref. ) to 20 Myr (ref. ). The lithium depletion boundary of the pm-RA =​ 4.65 ±​ 0.11 mas yr− pm-RA =​ 2.49 ±​ 0.68 mas yr− BPMG points24 to an age of 21 ​ 4 Myr. The older the system and the −1 −1 ± pm-dec. =​ 83.10 ±​ 0.15 mas yr pm-dec. =​ 82.58 ±​ 0.68 mas yr lower the initial energy, the higher is the theoretically inferred mass Hipparcos Gaia (1991.25–2015.5) 8 + of the planet . This has led to mass estimates for β​ Pic b of 4–11 MJ 1 25 3 26 24 pm-RA =​ 4.94 ±​ 0.02 mas yr− (ref. ), 6–12 MJ (ref. ), 8–13 MJ (ref. ) and 10–17 MJ (ref. ), but 27 1 also, recently, 12.9 ​ 0.2 M (ref. ). The most recent mass estimates pm-dec. =​ 83.93 ±​ 0.02 mas yr− ± J are in line with the model-independent measurement presented The RA, dec., parallax (plx) and proper motion in the RA and dec. direction (pm-RA and pm-dec., here, implying that the underlying assumptions in the hot start mod- respectively) are shown. els and the age estimates of the system are largely correct.

Nature Astronomy | www.nature.com/natureastronomy Nature Astronomy Letters

β Pic b mass 21 yr, the Hipparcos observations fall around the epoch of longest elongation, 40 implying a strong acceleration/deceleration of its angular motion, which is excluded by the Hipparcos data. We subsequently performed a chi-squared analysis on the data presented in the 30 right panel of Fig. 2 over a fine grid of orbital period and planet mass (assuming a ) J stellar mass31 of 1.75 solar masses), the results of which are shown in Fig. 3, with the 1σ​, 2σ ​and 3σ ​uncertainty intervals indicated by the contours. The green area is 20 the 1σ ​uncertainty interval for the orbital period as determined from high-contrast imaging of the planet18. We checked to what extent the deviation in the Gaia parallax from that Planet mass ( M 10 of Hipparcos of −​0.82 mas can influence the Hipparcos–Gaia proper-motion measurement. Gaia DR2 quotes a negligible RA/parallax correlation, and a dec./parallax correlation of +​0.248, implying a deviation in position of ​0.82 mas ​ 0.248 ​ ​0.2 mas. This results in a change in the Hipparcos–Gaia 0 − × =− 18 20 22 24 26 28 proper motion, which has a baseline of 24.25 yr, of only −​0.2 mas ÷​ 24.25 yr =​ −1 Orbital period (yr) −​0.008 mas yr .

Data availability. The data that support the plots within this paper and other Fig. 3 | Constraints on the mass and orbital period of the young exoplanet findings are available in Supplementary Table 1 and the Gaia data archive https:// β Pic b. The contours show the single-parameter 1σ​, 2σ​ and 3σ​ uncertainty gea.esac.esa.int/archive/. areas from the astrometric analysis presented here, with the blue area highlighting the 1σ ​region. Constraints on the orbital period from Code availability. All analyses were performed with custom-built IDL scripts. high-contrast imaging of the planet are shown in green. The planet is determined to have a mass of 11 ​ 2 M (1 ),​ and the orbital period is Received: 5 July 2018; Accepted: 23 July 2018; ± J σ Published: xx xx xxxx likely to be >​22.2 yr (2σ ​limit). Longer orbital periods are permitted by the astrometric data, except for ≫​30 yr when substantial acceleration/ References deceleration in the motion of the star would be expected during the 1. Marois, C. et al. Direct imaging of multiple planets orbiting the star HR 8799. time of the Hipparcos observations. Science 322, 1348–1352 (2008). 2. Lagrange, A.-M. et al. A probable giant planet imaged in the β​ Pictoris disk. VLT/NaCo deep L′​-band imaging. Astron. Astrophys. 493, L21–L25 (2009). The ongoing data collection by Gaia will, in the near future, lead 3. Lagrange, A.-M. et al. A giant planet imaged in the disk of the young star β​ to a further improvement on the mass measurement of β ​Pic b by up Pictoris. Science 329, 57–59 (2010). to one order of magnitude, and the masses of other known young 4. Macintosh, B. et al. Discovery and spectroscopy of the young jovian planet 51 17,28 Eri b with the Gemini Planet Imager. Science 350, 64–67 (2015). planets will be in reach . This means that Gaia will become an 5. Chabrier, G., Barafe, I., Allard, F. & Hauschildt, P. Evolutionary models for important asset to further constrain early-evolution scenarios. In very low-mass stars and brown dwarfs with dusty atmospheres. Astrophys. J. addition, Gaia is poised29 to find many new extrasolar planets. 542, 464–472 (2000). 6. Barafe, I., Chabrier, G., Barman, T. S., Allard, F. & Hauschildt, P. H. Evolutionary models for cool brown dwarfs and extrasolar giant planets. Methods Te case of HD 209458. Astron. Astrophys. 402, 701–712 (2003). We collected the Hipparcos intermediate data products from the van Leeuwen 2007 7. Marley, M. S., Fortney, J. J., Hubickyj, O., Bodenheimer, P. & Lissauer, J. J. reduction14 from a compact disc, which contain for each observation the spacecraft On the luminosity of young Jupiters. Astrophys. J. 655, 541–549 (2007). orbit and epoch, the cosine and sine of the scan orientation, the residual to the best- 8. Spiegel, D. S. & Burrows, A. Spectral and photometric diagnostics of giant fit solution and the formal error (see Supplementary Table 1). All calculations were planet formation scenarios. Astrophys. J. 745, 174 (2012). performed with custom-built IDL scripts. Figure 1 was produced using the make_ 9. Lagrange, A.-M. et al. Constraints on planets around β​ Pic with Harps radial parallax_coords.pro procedure30 to convert the Hipparcos barycentric position of velocity data. Astron. Astrophys. 542, A18 (2012). β​ Pic (epoch 1991.25), the intermediate data products, and the best-fit solution to 10. Lagrange, A.-M. et al. Full exploration of the giant planet population around the parallax and proper motion into a model of the on-sky geocentric movement of β​ Pictoris. Astron. Astrophys. 612, A108 (2018). the star and the accompanying one-dimensional astrometric measurements. The 11. Mouillet, D., Larwood, J. D., Papaloizou, J. C. B. & Lagrange, A. M. A planet Gaia (epoch 2015.5) barycentric position was extracted from the DR2 archive and on an inclined orbit as an explanation of the warp in the β​ Pictoris disc. combined with that of Hipparcos to determine the average proper motion of the Mon. Not. R. Astron. Soc. 292, 896–904 (1997). star over the intermediate 24.25 yr. 12. Fabrycky, D. C. & Murray-Clay, R. A. Stability of the directly imaged The 111 one-dimensional astrometric measurements of Hipparcos, taken multiplanet system HR 8799: resonance and masses. Astrophys. J. 710, during 35 spacecraft orbits, were combined into five two-dimensional positions 1408–1421 (2010). (Fig. 2, left panel), combining observations with seven different scanning 13. Perryman, M. A. C. et al. Te HIPPARCOS catalogue. Astron. Astrophys. 323, directions per position. This was done by combining the individual measurements L49–L52 (1997). per spacecraft orbit to produce a (Δ​RA, Δ​dec.) grid of Δ​χ2 values, assuming a 14. van Leeuwen, F. Validation of the new Hipparcos reduction. Astron. Astrophys. Gaussian noise distribution along the scanning direction. These Δχ​ 2 grid maps 474, 653–664 (2007). were subsequently combined into Δ​χ2 =​ 1 contour for each of the five epochs, as 15. Gaia Collaboration. Te Gaia mission. Astron. Astrophys. 595, A1 (2016). shown in the left panel of Fig. 2. It shows that the movement of the star is in the 16. Gaia Collaboration. Gaia Data Release 2. Summary of the contents and expected direction for the reflex motion induced by planet β ​Pic b, as indicated survey properties. Astron. Astrophys. (in the press). by the arrow. 17. Science Performance (European Space Agency, 2018); https://www.cosmos.esa. Since the direction of the reflex motion of the host star is well determined18 int/web/gaia/science-performance (orientation angle, φ =​ 212°), we converted each orbit-averaged observation into a 18. Wang, J. J. et al. Te orbit and prospects for β​ Pictoris b constrained position measurement along this direction and adjusted its associated uncertainty, with one milliarcsecond astrometry. Astron. J. 152, 97 (2016). which depends on 1/cosine of the angle relative to the scanning direction. This 19. Chauvin, G. et al. Orbital characterization of the β​ Pictoris b giant planet. is shown in the right panel of Fig. 2, where the data points have grey scales, Astron. Astrophys. 542, A41 (2012). depending on their uncertainty, to emphasize the data with the smallest errors. 20. Nielsen, E. L. et al. Te Gemini NICI Planet-Finding Campaign: the orbit of Five measurements were omitted because their scanning directions were nearly the young exoplanet β​ Pictoris b. Astrophys. J. 794, 158 (2014). perpendicular, leading to very large errors of >​10 mas. The shape of the orbit of 21. Zuckerman, B., Song, I., Bessell, M. S. & Webb, R. A. Te β​ Pictoris moving the planet is well determined for the Hipparcos 1990–1993 observing interval group. Astrophys. J. 562, L87–L90 (2001). because it is approximately one planet orbit (20–26 yr) before the current era 22. Ortega, V. G., de la Reza, R., Jilinski, E. & Bazzanella, B. Te origin of the (2010–present), during which the planet’s position has been frequently monitored β​ Pictoris moving group. Astrophys. J. 575, L75–L78 (2002). (see Supplementary Fig. 1). Therefore, although the orbital period is not well 23. Mamajek, E. E. & Bell, C. P. M. On the age of the β​ Pictoris moving group. constrained, the best-fit18 solution to its orbit can be used to fit the reflex motion Mon. Not. R. Astron. Soc. 445, 2169–2180 (2014). of the star, modulo a scaling depending on the planet/star mass ratio and a phase 24. Binks, A. S. & Jefries, R. D. A lithium depletion boundary age of shift governed by the assumed orbital period. This is shown for trial orbital periods 21 Myr for the Beta Pictoris moving group. Mon. Not. R. Astron. Soc. 438, between 21 yr and 24 yr in the right panel of Fig. 2. For a short orbital period of L11–L15 (2014).

Nature Astronomy | www.nature.com/natureastronomy Letters Nature Astronomy

25. Currie, T. et al. A combined and Gemini study of https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has the atmosphere of the directly imaged planet, β​ Pictoris b. Astrophys. J. 776, been provided by national institutions, in particular the institutions participating in 15 (2013). the Gaia Multilateral Agreement. I.A.G.S. acknowledges funding from the European 26. Bonnefoy, M. et al. Te near- spectral energy distribution of Research Council under the European Union’s Horizon 2020 research and innovation β​ Pictoris b. Astron. Astrophys. 555, A107 (2013). programme under grant agreement No. 694513. A.G.A.B. acknowledges funding from 27. Chilcote, J. et al. 1–2.4 μ​m near-IR spectrum of the giant planet β​ Pictoris b the Netherlands Research School for Astronomy (NOVA). obtained with the Gemini Planet Imager. Astron. J. 153, 182 (2017). 28. Michalik, D., Lindegren, L., Hobbs, D. & Lammers, U. Joint astrometric solution of HIPPARCOS and Gaia. A recipe for the Hundred Tousand Author contributions Proper Motions project. Astron. Astrophys. 571, A85 (2014). I.A.G.S. devised the general idea, conducted the main analysis and wrote a first version 29. Perryman, M., Hartman, J., Bakos, G. Á. & Lindegren, L. Astrometric of the manuscript. A.G.A.B. contributed to the implementation and optimization of the exoplanet detection with Gaia. Astrophys. J. 797, 14 (2014). analysis, including important details on the Gaia and Hipparcos data, and commented on 30. make_parallax_coords (NASA/IPAC, accessed 1 March 2018); https://irsa. the manuscript. ipac.caltech.edu/data/SPITZER/docs/dataanalysistools/tools/contributed/ general/make_parallax_coords/ Competing interests 31. Crifo, F., Vidal-Madjar, A., Lallement, R., Ferlet, R. & Gerbaldi, M. The authors declare no competing interests β​ Pictoris revisited by Hipparcos. Star properties. Astron. Astrophys. 320, L29–L32 (1997). Additional information Acknowledgements Supplementary information is available for this paper at https://doi.org/10.1038/ This work was only possible because numerous scientists and engineers have devoted s41550-018-0561-6. large parts of their careers to the design, construction and successful operation of the Reprints and permissions information is available at www.nature.com/reprints. Hipparcos and Gaia missions and the analysis of their data. We are grateful. We thank F. van Leeuwen also for discussions on his analysis of the Hipparcos data. This work has Correspondence and requests for materials should be addressed to I.A.G.S. made use of data from the European Space Agency mission Gaia (https://www.cosmos. Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, published maps and institutional affiliations.

Nature Astronomy | www.nature.com/natureastronomy