Published February 1, 2025 | Version Published
Journal Article Open

Are Stars Really Ingesting Their Planets? Examining an Alternative Explanation

  • 1. ROR icon California Institute of Technology

Abstract

Numerous stars exhibit surprisingly large variations in their refractory element abundances, often interpreted as signatures of planetary ingestion events. In this study, we propose that differences in the dust-to-gas ratio near stars during their formation can produce similar observational signals. We investigate this hypothesis using a suite of radiation-dust-magnetohydrodynamic STAR FORmation in Gaseous Environments (or STARFORGE) simulations of star formation. Our results show that the distribution of refractory abundance variations (Δ[X/H]) has extended tails, with about 10%–30% of all stars displaying variations around ∼0.1 dex. These variations are comparable to the accretion of 2–5 M of planetary material into the convective zones of Sun-like stars. The width of the distributions increases with the incorporation of more detailed dust physics, such as radiation pressure and back-reaction forces, as well as with larger dust grain sizes and finer resolutions. Furthermore, our simulations reveal no correlation between Δ[X/H] and stellar separations, suggesting that dust-to-gas fluctuations likely occur on scales smaller than those of wide binaries. These findings highlight the importance of considering dust dynamics as a potential source of the observed chemical enrichment in stars.

Copyright and License

© 2025. The Author(s). Published by the American Astronomical Society. Original content from this work may be used under the terms of the Creative Commons Attribution 4.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.

Acknowledgement

We thank Caleb R. Choban for his valuable feedback on the clarity and organization of this work. We also thank Carlos Saffe for his insightful suggestions and observations, which helped us refine our research. Support for for N.S. and P.F.H. was provided by NSF Research Grants 1911233, 20009234, and 2108318; an NSF CAREER grant 1455342; and NASA grant Nos. 80NSSC18K0562, and HST-AR-15800. Numerical calculations were run on the TACC compute cluster “Frontera,” allocations AST21010 and AST20016, supported by the NSF and TACC, and NASA HEC SMD-16-7592.

Software References

matplotlib (J. D. Hunter 2007), NumPy (C. R. Harris et al. 2020), SciPy (P. Virtanen et al. 2020).

Data Availability

Data necessary to reproduce the results of this study, including simulation outputs and processed data used for figures, are available upon reasonable request to the corresponding author.

Files

Soliman_2025_ApJ_979_98.pdf

Files (745.8 kB)

Name Size Download all
md5:a9eddbe89e10478fcdf2a2d7b0aa4cfb
745.8 kB Preview Download

Additional details

Related works

Is new version of
Discussion Paper: arXiv:2408.15326 (arXiv)

Funding

National Science Foundation
1911233
National Science Foundation
20009234
National Science Foundation
2108318
National Science Foundation
1455342
National Science Foundation
1455342
National Aeronautics and Space Administration
80NSSC18K0562
National Aeronautics and Space Administration
HST-AR-15800
National Aeronautics and Space Administration
HEC SMD-16-7592

Dates

Submitted
2024-08-19
Accepted
2024-12-15
Available
2025-01-21
Published online

Caltech Custom Metadata

Caltech groups
Astronomy Department, TAPIR, Walter Burke Institute for Theoretical Physics, Division of Physics, Mathematics and Astronomy (PMA)
Publication Status
Published