Are Stars Really Ingesting Their Planets? Examining an Alternative Explanation
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
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-21Published online