Open Access
Issue
A&A
Volume 621, January 2019
Article Number L8
Number of page(s) 7
Section Letters to the Editor
DOI https://doi.org/10.1051/0004-6361/201834302
Published online 11 January 2019

© ESO 2019

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

With its imaged debris disk of dust (see Smith & Terrile 1984, for the discovery image), its falling, evaporating exocomets (Kiefer et al. 2014, and references therein), and an imaged giant planet (Lagrange et al. 2010), the ∼20 Myr old β Pictoris is a unique proxy for the study of the early stages of planetary system formation and evolution, when giant planets are formed, Earth-mass planets may still be forming, and most of the protoplanetary gas has disappeared from the disk. Its proximity to Earth (van Leeuwen 2007, distance = 19.454 ± 0.05 pc) and the relatively short (≃ two decades) orbital period of β Pictoris b enable detailed studies of its orbit and its physical and atmospheric properties. The system also allows us to study the interaction between planet(s) and disks. β Pictoris b can explain, for example, several (but not all) of the dust disk morphologies, in particular its inner warp, and some outer asymmetries (Lagrange et al. 2010). β Pictoris b could also be the trigger for the infall and evaporation of cometary bodies (exocomets) onto the star, if it has a non-zero eccentricity (Beust & Morbidelli 1996). Last, β Pictoris b was suggested to be responsible for the photometric variations observed in 1981 (Lecavelier Des Etangs et al. 1995) and has tentatively been attributed to a planet transit (Lecavelier Des Etangs et al. 1997).

Careful monitoring of the position of the planet relative to the star (referred to as astrometric measurements) with the Nasmyth Adaptive Optics System (NAOS) Near-Infrared Imager and Spectrograph (CONICA) (NaCo) at the Very Large Telescope (VLT) constrained its orbital properties (Bonnefoy et al. 2014; Chauvin et al. 2012). By combining data from NaCo and the Gemini Planetary Imager (GPI), Wang et al. (2016) found a similar bet slightly different orbit. Assuming that the peculiar photometric event observed in 1981 could be due to the transit of β Pictoris b in front of the star, Lecavelier des Etangs & Vidal-Madjar (2016) identified a second family of orbital solutions: a semi-major axis a = 13 au instead of 9 au, that is, a period of 34 yr instead of about 20 yr, and e = 0.3 instead of less than 0.1. Recently, Snellen & Brown (2018) used Gaia and HIPPARCOS measurements to constrain the planet period to ≥22 yr, and reported a mass of 11 ± 2 MJup, which is compatible with the constraints derived from radial velocity data alone (Lagrange et al. 2012a) or by the combined analysis of direct imaging and radial velocity data (Bonnefoy et al. 2014).

β Pictoris b was discovered northeast of the star in data obtained in 2003 with NaCo. All available images in addition to this discovery image were obtained in 2009 and later, as the planet orbited southwest of the star, after it passed behind the star (Lagrange et al. 2010). The orbital plane is close to equatorial, and the disk of β Pictoris is seen almost edge-on. This geometrical configuration prevented following the planet in direct imaging when it was projected too close (typically less than 120 mas) to the star. Altogether, only about 50% of the planetary orbit had been monitored thus far, including the 2003 epoch, which suffers from relatively large uncertainties and thus limits the precision on the orbital parameters of the planet. We here present a homogeneous set of planet observations obtained with the Spectro-Polarimetric High-Contrast Exoplanet Research instrument (SPHERE), and in particular, the recent recovery of the planet in September 2018. The observations are described in Sect. 2, and the results are shown and discussed in Sect. 3.

2. Observations

High-contrast coronagraphic SPHERE (Beuzit et al. 2008) observations were obtained between December 2014 and November 2016, and in September 2018, using the IRDIFS mode in the context of the SpHere INfrared survey for Exoplanets (SHINE; Chauvin et al. 2017a). In this setup, the IRDIS (Dohlen et al. 2008) and IFS (Claudi et al. 2008) instruments operate simultaneously. The data were obtained under various atmospheric conditions (see Table 1) with the H2 (λc = 1.593 μm; Δλ = 0.055 μm) and H3 (λc = 1.667 μm; Δλ = 0.056 μm) narrow bands of IRDIS, except in December 2014, when we used the K1 (λc = 2.1025 μm; Δλ = 0.102 μm) and K2 (λc = 2.255 μm; Δλ = 0.109 μm) narrow bands. The IRDIS images have a field of view (FoV) of ∼10″ × 11″, and a pixel size of approximately 12.25 mas. IFS data were also recorded, but they are not analyzed here. We used apodized Lyot coronagraphs that include either a 185 mas diameter focal mask (N_ALC_YJH_S) or, when the planet was closer to the star, a smaller (145 mas) mask (N_ALC_YJ_S), combined to an apodizer as well as a pupil stop (Carbillet et al. 2011). All coronagraphic data were recorded in stabilized pupil mode so as to apply angular differential imaging (ADI) post-processing techniques to remove the stellar halo, as described in Marois et al. (2006). Most coronagraphic data were also recorded while four satellite footprints of the point spread function (PSF) had been created by the deformable mirror of the instrument and were used for fine monitoring of the frame centering and for photometric purposes. The FoV rotation during the coronagraphic observations varied between 16° and 54° (see Table 1).

Table 1.

Observing log.

Each observing sequence was obtained with the following pattern: PSF – coronagraphic observations – PSF – sky. The PSF data correspond to non-saturated exposures of the star placed out of the coronagraphic mask and obtained using a neutral density filter. They are used for relative photometric reference and to estimate the image quality at the beginning and end of the observations. The sky data were recorded at the end of the coronagraphic sequence to estimate the background level and hot pixels in the science images. Finally, an astrometric field, either Orion or 47 Tuc, was observed with IRDIS for each run (see Maire et al. 2016). In the case of Orion, the sub-field we used was chosen to be part of the one considered with NaCo since 2008 to allow the best match between the astrometric calibrations of NaCo and SPHERE. The pixel scales and north positions are provided in Table 1.

The data were reduced as described in Chauvin et al. (2017b) and using the SpeCal tool developed for SPHERE (Galicher et al. 2018). β Pictoris b is clearly detected in all images taken between 2014 and 2016, orbiting SW of the star at signal-to-noise ratio (S/N) higher than 9 (Table 2). Figure 1 shows images of the planet at various dates, and Table 2 provides the relative position of the planet with respect to the star1. The S/N is relatively poor in the last observation of November 2016, as the projected separation of the planet from the star is less than 125 mas (≃1.5 au only) and the contrast is about 9.5 mag, which leads to larger error bars on its astrometry. To our knowledge, neither β Pictoris b nor any other planet has ever been imaged at such a close projected separation to the star. The 2018 data clearly reveal the planet at 139 mas NE (PA about 30°) from the star.

thumbnail Fig. 1.

Images of β Pictoris b with SPHERE IRDIS from December 2014 to September 2018. Each panel displays a FoV of 1″ × 1″. North is up and east is to the left. The intensity scale is adapted at each epoch according to the intensity peak of the planet.

Table 2.

Relative astrometry of β Pictoris b.

3. Orbital properties of β Pictoris b

With these additional SPHERE data, a large part of the orbit is now sampled, as seen in Fig. 2. We used the positions of the planet relative to the star in these data as well as in the previously published NaCo data (Bonnefoy et al. 2014) together with the Markov chain Monte Carlo (MCMC) Bayesian analysis technique described in Chauvin et al. (2012) to derive the probabilistic distribution of the orbital solutions. The results are shown in Fig. A.1. For comparison, we show in Fig. A.2 the parameters that we deduce with and without the September 2018 data to illustrate the importance of this post-conjunction detection for constraining the orbital properties of the planet.

thumbnail Fig. 2.

NaCo (blue) and SPHERE (magenta) astrometric data points of β Pictoris b shown together with 200 probable solutions of the MCMC analysis (gray) and the best-fit Levenberg–Marquardt solution (black).

From Fig. A.2, we derive the following 1σ confidence intervals for the major orbital parameters: au for the semi-major axis, yr for the orbital period, for the eccentricity, and for the inclination. We note that the orbital period distribution is compatible with a period greater than 22 yr found by Snellen & Brown (2018) within only 2σ error bars.

The sma found is compatible with the value derived by Bonnefoy et al. (2014), and definitely excludes the long-period solution proposed by Lecavelier des Etangs & Vidal-Madjar (2016) that assumed a transit of the planet in 1981. Snellen & Brown (2018) predicted a period longer than 22 yr, while the period we derive here is rather 20 yr. We note, however, that the result from the Snellen & Brown (2018) relies on the assumption that there is only one planet around β Pictoris, which is not necessarily correct. Finally, based on our values, the 2017 conjunction took place in 2017.72 ± 0.04, and the next conjunction will occur in 2038.06.

The orbital solution favors a very low eccentricity but still remains compatible with zero. We note that the peak of the distribution is now slightly off zero, in contrast with previous determinations (Bonnefoy et al. 2014). Previous studies (Beust & Morbidelli 1996, 2000; Thébault & Beust 2001) of the falling evaporating bodies scenario showed that this phenomenon could be explained by the perturbing action of a giant planet orbiting at ∼10 au from the star, if it has a low orbital eccentricity. β Pictoris b might of course be assimilated to this planet, and the fact that its eccentricity might be non-zero is a key point for this scenario.

Finally, the large FoV of IRDIS enables detecting the disk and the planet in the same image. The SPHERE data confirm the conclusions reached in Lagrange et al. (2012b): the planet projection is between the main disk and the warp.

4. Concluding remarks and perspectives

The sensitivity of SPHERE allowed us to follow β Pictoris b down to 125 mas from the star in projected separation. The latest measurements reveal the planet on the NE side of the star, for the first time since its discovery. It was last detected in November 2016 and was redetected in September 2018. Based on the observed data, the semi-major axis of the planet is well constrained to au (1σ), its eccentricity to , and its inclination to deg. The data do not support that the planet was responsible for the photometric event in 1981. More data obtained in the NE part of the disk will allow further refining the orbital properties of the planet. Further work will include the combination of these new, crucial astrometric data points with our latest radial velocity measurements and, if possible, a combination with Gaia and HIPPARCOS data to constrain the dynamical mass of the planet. From an instrumental point of view, these data demonstrate that if they exist, SPHERE can detect young and massive Jupiters as close as 1.6 au from a star located at 20 pc.


1

The planet was hardly or even not at all detected in November 2016 with the TLOCI algorithm (Marois et al. 2014), while it was detected using a principal component analysis (PCA; Soummer et al. 2012) algorithm. For consistency, all position measurements were obtained on PCA images.

Acknowledgments

We acknowledge financial support from the Programme National de Planétologie (PNP) and the Programme National de Physique Stellaire (PNPS) of CNRS-INSU. This work has also been supported by a grant from the French Labex OSUG@2020 (Investissements d’avenir ANR10 LABX56). The project is supported by CNRS, by the Agence Nationale de la Recherche (ANR-14-CE33-0018). This work has made use of the the SPHERE Data Centre, jointly operated by OSUG/IPAG (Grenoble), PYTHEAS/LAM/CESAM (Marseille), OCA/Lagrange (Nice) and Observatoire de Paris/LESIA (Paris). AML thanks Julian Mejia for his help during the September 2018 DVM observations, Nadège Meunier, Eric Lagadec, and Philippe Delorme for their assistance in the DC process, and Pascal Rubini for his help in the scheduling of the observations. E.S., R.G., D.M., S.D. and R.U.C. acknowledge support from the “Progetti Premiali” funding scheme of the Italian Ministry of Education, University, and Research. SPHERE is an instrument designed and built by a consortium consisting of IPAG (Grenoble,France), MPIA (Heidelberg, Germany), LAM (Marseille, France), LESIA (Paris, France), Laboratoire Lagrange (Nice, France), INAF Osservatorio Astronomico di Padova (Italy), Observatoire de Genève (Switzerland), ETH Zurich (Switzerland), NOVA (Netherlands), ONERA (France) and ASTRON (Netherlands) in collaboration with ESO. SPHERE was funded by ESO, with additional contributions from CNRS (France), MPIA (Germany), INAF (Italy), FINES (Switzerland) and NOVA (Netherlands). SPHERE also received funding from the European Commission Sixth and Seventh Framework Programmes as part of the Optical Infrared Coordination Network for Astronomy (OPTICON) under grant number RII3-Ct-2004-001566 for FP6 (2004–2008), grant number 226604 for FP7 (2009–2012) and grant number 312430 for FP7 (2013–2016).

References

  1. Beust, H., & Morbidelli, A. 1996, Icarus, 120, 358 [NASA ADS] [CrossRef] [Google Scholar]
  2. Beust, H., & Morbidelli, A. 2000, Icarus, 143, 170 [NASA ADS] [CrossRef] [Google Scholar]
  3. Beuzit, J.-L., Feldt, M., Dohlen, K., et al. 2008, Proc. SPIE, 7014, 701418 [Google Scholar]
  4. Bonnefoy, M., Marleau, G.-D., Galicher, R., et al. 2014, A&A, 567, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Bonnefoy, M., Marleau, G.-D., Galicher, R., et al. 2014b, A&A, 567, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Carbillet, M., Bendjoya, P., Abe, L., et al. 2011, Exp. Astron., 30, 39 [NASA ADS] [CrossRef] [Google Scholar]
  7. Chauvin, G., Lagrange, A.-M., Beust, H., et al. 2012, A&A, 542, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Chauvin, G., Desidera, S., Lagrange, A. M., et al. 2017a, in SF2A-2017: Proceedings of the Annual Meeting of the French Society of Astronomy and Astrophysics, eds. C. Reylé, P. Di Matteo, F. Herpin, et al., 331 [Google Scholar]
  9. Chauvin, G., Desidera, S., Lagrange, A.-M., et al. 2017b, A&A, 605, L9 [CrossRef] [EDP Sciences] [Google Scholar]
  10. Claudi, R. U., Turatto, M., Gratton, R. G., et al. 2008, Proc. SPIE, 7014, 70143E [Google Scholar]
  11. Dohlen, K., Langlois, M., Saisse, M., et al. 2008, Proc. SPIE, 7014, 70143L [Google Scholar]
  12. Galicher, R., Boccaletti, A., Mesa, D., et al. 2018, A&A, 615, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Kiefer, F., Lecavelier des Etangs, A., Boissier, J., et al. 2014, Nature, 514, 462 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Lagrange, A.-M., Bonnefoy, M., Chauvin, G., et al. 2010, Science, 329, 57 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  15. Lagrange, A.-M., De Bondt, K., Meunier, N., et al. 2012a, A&A, 542, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Lagrange, A.-M., Boccaletti, A., Milli, J., et al. 2012b, A&A, 542, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Lecavelier des Etangs, A., & Vidal-Madjar, A. 2016, A&A, 588, A60 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Lecavelier Des Etangs, A., Deleuil, M., Vidal-Madjar, A., et al. 1995, A&A, 299, 557 [NASA ADS] [Google Scholar]
  19. Lecavelier Des Etangs, A., Vidal-Madjar, A., Burki, G., et al. 1997, A&A, 328, 311 [NASA ADS] [Google Scholar]
  20. Maire, A.-L., Langlois, M., Dohlen, K., et al. 2016, in Ground-based and Airborne Instrumentation VI, Proc. SPIE, 9908, 990834 [Google Scholar]
  21. Marois, C., Lafrenière, D., Doyon, R., Macintosh, B., & Nadeau, D. 2006, ApJ, 641, 556 [Google Scholar]
  22. Marois, C., Correia, C., Véran, J.-P., & Currie, T. 2014, IAU Symp., 299, 48 [NASA ADS] [Google Scholar]
  23. Smith, B. A., & Terrile, R. J. 1984, Science, 226, 1421 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  24. Snellen, I. A. G., & Brown, A. G. A. 2018, Nat. Astron., 2, 883 [NASA ADS] [CrossRef] [Google Scholar]
  25. Soummer, R., Pueyo, L., & Larkin, J. 2012, ApJ, 755, L28 [NASA ADS] [CrossRef] [Google Scholar]
  26. Thébault, P., & Beust, H. 2001, A&A, 376, 621 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. van Leeuwen, F. 2007, A&A, 474, 653 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Wang, J. J., Graham, J. R., Pueyo, L., et al. 2016, AJ, 152, 97 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Additional figures

thumbnail Fig. A.1.

Results of the MCMC fit of the NaCo + SPHERE astrometric data points. The star and the red dotted line correspond to the best-fit solution (best χ2) obtained with a Levenberg–Marquardt fit.

thumbnail Fig. A.2.

Comparison between the orbital parameters obtained with (black) or without (red) the recovery point of September 2018.

All Tables

Table 1.

Observing log.

Table 2.

Relative astrometry of β Pictoris b.

All Figures

thumbnail Fig. 1.

Images of β Pictoris b with SPHERE IRDIS from December 2014 to September 2018. Each panel displays a FoV of 1″ × 1″. North is up and east is to the left. The intensity scale is adapted at each epoch according to the intensity peak of the planet.

In the text
thumbnail Fig. 2.

NaCo (blue) and SPHERE (magenta) astrometric data points of β Pictoris b shown together with 200 probable solutions of the MCMC analysis (gray) and the best-fit Levenberg–Marquardt solution (black).

In the text
thumbnail Fig. A.1.

Results of the MCMC fit of the NaCo + SPHERE astrometric data points. The star and the red dotted line correspond to the best-fit solution (best χ2) obtained with a Levenberg–Marquardt fit.

In the text
thumbnail Fig. A.2.

Comparison between the orbital parameters obtained with (black) or without (red) the recovery point of September 2018.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.