Welcome to the new version of CaltechAUTHORS. Login is currently restricted to library staff. If you notice any issues, please email coda@library.caltech.edu
Published March 21, 2024 | Published
Journal Article Open

At least one in a dozen stars shows evidence of planetary ingestion

Abstract

Stellar chemical compositions can be altered by ingestion of planetary material1,2 and/or planet formation, which removes refractory material from the protostellar disk3,4. These ‘planet signatures’ appear as correlations between elemental abundance differences and the dust condensation temperature3,5,6. Detecting these planet signatures, however, is challenging owing to unknown occurrence rates, small amplitudes and heterogeneous star samples with large differences in stellar ages7,8. Therefore, stars born together (that is, co-natal) with identical compositions can facilitate the detection of planet signatures. Although previous spectroscopic studies have been limited to a small number of binary stars9,10,11,12,13, the Gaia satellite14 provides opportunities for detecting stellar chemical signatures of planets among co-moving pairs of stars confirmed to be co-natal15,16. Here we report high-precision chemical abundances for a homogeneous sample of ninety-one co-natal pairs of stars with a well defined selection function and identify at least seven instances of planetary ingestion, corresponding to an occurrence rate of eight per cent. An independent Bayesian indicator is deployed, which can effectively disentangle the planet signatures from other factors, such as random abundance variation and atomic diffusion17. Our study provides evidence of planet signatures and facilitates a deeper understanding of the star–planet–chemistry connection by providing observational constraints on the mechanisms of planet engulfment, formation and evolution.

Copyright and License

© The Author(s), under exclusive licence to Springer Nature Limited 2024.

Acknowledgement

This paper includes data gathered with the 6.5 metre Magellan Telescope located at Las Campanas Observatory, Chile. Some of the data presented herein were obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California and the National Aeronautics and Space Administration. The observatory was made possible by the generous financial support of the W. M. Keck Foundation. Based on observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere under ESO programme 108.22EC.001. This research were supported by the Australian Research Council Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), through project number CE170100013. Y.-S.T. acknowledges financial support from the Australian Research Council through DECRA Fellowship DE220101520. M.T.M. acknowledges the support of the Australian Research Council through Future Fellowship grant FT180100194. B.B. thanks the European Research Council (ERC Starting Grant 757448-PAMDORA) for their financial support. M.J. gratefully acknowledges funding of MATISSE: Measuring Ages Through Isochrones, Seismology, and Stellar Evolution, awarded through the European Commission’s Widening Fellowship. This project has received funding from the European Union’s Horizon 2020 research and innovation programme. We thank A. Ji for offering advice on data collection and preparation; and S. Campbell, A. Mustill and Q. Sun for discussions. The C3PO programme is made possible through the Carnegie Observatories’ support and allocation of observation time on the Magellan Telescope. We recognize and acknowledge the very significant cultural role and reverence that the summit of Maunakea has always had within the Indigenous Hawaiian community. We are most fortunate to have the opportunity to conduct observations from this mountain. This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement.

Contributions

F.L. led and played a part in all aspects of the observations and data analysis for this study, and wrote and developed the paper. Y.-S.T. initiated the C3PO program, carried out the Magellan observations and contributed to the statistical analysis of the research. D.Y. carried out the observations and part of spectroscopic analysis of the Magellan and the VLT data, and contributed to designing this study. B.B., M.J., and A.D. contributed to the theoretical interpretations of the observational results. A.K., M.T.M. and F.D. contributed to the development and writing of the paper. All authors read, commented and agreed on the paper.

Data Availability

The spectral data underlying this article are available in Keck Observatory Archive (https://koa.ipac.caltech.edu/cgi-bin/KOA/nph-KOAlogin) and the European Southern Observatory Science Archive Facility (http://archive.eso.org/eso/eso_archive_main.html). They can be accessed with Keck Program ID W244Hr (semester: 2021B, principal investigator F.L.) and European Southern Observatory Programme ID 108.22EC.001 (principal investigator D.Y.), respectively. The spectra from the Magellan Telescope can be shared upon request to the corresponding author. The rest of data underlying this article are available in the article and its Supplementary Information.

Extended Data Fig. 1 A portion of the reduced spectra of an example co-moving pair of stars (HD 185726/185689; Pair 124).

Extended Data Fig. 2 Comparison of abundance results.

Extended Data Fig. 3 The distributions of differences in Bayesian evidence Δln(Z) between the planetary ingestion and the flat models.

Extended Data Fig. 4 Differences in Bayesian evidence Δln(Z)atom between the planetary ingestion and the atomic diffusion models as a function of spatial separations Δs.

Extended Data Fig. 5 Total accreted mass of Earth material (from the best-fitting model) as a function of Δln(Z) between the planetary ingestion and the flat models for the 91 close, co-natal co-moving pairs.

Extended Data Fig. 6 Abundance differences (Δ[X/H]) in our candidate pairs.

Extended Data Fig. 7 Abundance differences (Δ[X/H]) in our candidate pairs.

Extended Data Fig. 8 Abundance differences (Δ[X/H]) in our candidate pairs.

Extended Data Table 1 Fitting results of seven co-moving pairs of stars exhibiting distinct chemical signatures of ingestion of planetary ingestion

Code Availability

The stellar line analysis program MOOG is available at https://www.as.utexas.edu/~chris/moog.html. The stellar model atmospheres are available at http://kurucz.harvard.edu/grids.html. The code for equivalent width measurements is very similar to REVIEW, which is provided in https://github.com/madeleine-mckenzie/REvIEW. The Bayesian modelling program DYNESTY is available at https://github.com/joshspeagle/dynesty.

Conflict of Interest

The authors declare no competing interests.

Files

41586_2024_7091_MOESM2_ESM.csv
Files (366.3 kB)
Name Size Download all
md5:8faff493c2f6fa352286575b3bf837f8
315.7 kB Preview Download
md5:c5f46a4d0d11eb16114d09ab275ae5f1
50.6 kB Preview Download

Additional details

Created:
March 26, 2024
Modified:
March 26, 2024