Some stars eat its planets. The chemical paradox of HD 81809

In a binary star system, two stars born from the same molecular cloud carry an implicit contract: they should be chemically identical. Same metallicity, same abundance ratios, same history written in their elemental signatures. When that contract is violated — when one star looks measurably different from its companion — something must have happened after they formed. A new study targeting the binary system HD 81809 confronts exactly that violation, and the most compelling explanation on the table is that one of the stars swallowed its own planets.

The preprint, by Nuno Moedas and Maria Pia Di Mauro and uploaded to arXiv on May 29, 2026 (arXiv:2605.31060), builds on a detailed asteroseismic and spectroscopic characterization of the system published earlier this year in The Astrophysical Journal (Di Mauro et al. 2026, ApJ 1000, 92). Together, these two studies portray HD 81809 as one of the most chemically peculiar and scientifically productive binary systems accessible to current observations.

Two stars, one age, two very different compositions

HD 81809 is an old system with an age of roughly 10 Gyr. The primary, HD 81809A, is a subgiant with mass approximately 0.87 solar masses and radius approximately 1.96 solar radii. The secondary, HD 81809B, remains on the main sequence with mass around 0.85 solar masses and radius around 1.10 solar radii. Both are solar-like G-type stars. Both, in principle, should have emerged from the same reservoir of interstellar material a decade of gigayears ago.

Yet the HD 81809 system presents a peculiar chemical composition, with a large metallicity difference between its two components: the primary has low metallicity at [Fe/H] = -0.57 dex, while the secondary sits at approximately solar metallicity with [Fe/H] = 0.0 dex. In practical terms, a difference of 0.57 dex on the logarithmic [Fe/H] scale means HD 81809B carries roughly 3.7 times the iron-to-hydrogen ratio of its companion. For two stars that formed together as part of the galactic thick disk — a population characteristically metal-poor — this discrepancy is not a rounding error. It is a puzzle that demands a physical mechanism.

A thorough baseline from TESS and decades of monitoring

Before the new accretion study could be done, the system needed a rigorous physical characterization. Using multi-sector TESS photometry, solar-like oscillations were detected in both components, yielding global asteroseismic parameters for each star separately. Combined with high-resolution spectroscopy, Gaia astrometry, radial velocity measurements, and X-ray monitoring accumulated over decades, these data established the fundamental parameters of the system with a precision rarely achieved for binary pairs. The HD 81809 binary has an orbital period of approximately 34.5 years and also displays high-amplitude X-ray emission and a well-defined 8-year solar-like magnetic cycle. The subgiant primary is considered the likely driver of that reactivated dynamo cycle.

This detailed physical portrait is what gives the chemical modelling a firm footing. Without precise masses, radii, ages, and evolutionary states for both stars, simulating accretion events and comparing their outcomes to observations would yield ambiguous results.

The planetary engulfment scenario and what MESA says

Moedas and Di Mauro used the MESA stellar evolution code to model HD 81809B by introducing accretion events with rocky material whose composition was based on CI chondrites — primitive meteorites considered a close proxy for the bulk composition of rocky planetary bodies. Their question was direct: can planetary engulfment explain how HD 81809B ended up so much more metal-rich than its companion?

To reach the observed metallicity, the star must engulf between 25 and 75 Earth masses of metals near its current age. Accretion brings the effective temperature of the stellar models closer to the observed value. Both of these results point in the same direction: an engulfment event is physically consistent with what is seen. As a planetary system evolves over gigayears, gravitational interactions can destabilize orbits and send rocky bodies spiraling inward. The accreted material would settle into the outer convective layers of the star and raise its apparent metallicity in the photosphere, while the companion — having presumably lost its planets by other means or not formed them — retains the original low-metallicity signature.

Where the story complicates itself

Here the paradox deepens. Lithium is a fragile element destroyed inside stellar interiors at relatively modest temperatures, making its photospheric abundance a sensitive probe of mixing history and accretion. The lithium results in the accretion models come out over-enriched. In order to reproduce the observed lithium abundance of HD 81809B, the star must have accreted less than 6 Earth masses.

This creates a genuine tension: to match the metallicity, the accreted mass must be 25 to 75 Earth masses; to match the lithium, it must be below 6. No single value satisfies both constraints simultaneously if the accreted material has the CI chondrite composition assumed in the models.

The possible explanations include planet engulfment, suppressed depletion due to low rotational mixing, or an unusual internal structure. One natural resolution is that the actual composition of the accreted material was not identical to CI chondrites. Rocky material depleted in volatile elements — including lithium — but enriched in refractories like iron, silicon, and magnesium could simultaneously account for the metallicity enhancement without over-producing lithium. In other words, the simple story of swallowing generic rocky bodies may not work, but a more differentiated type of planetary material might.

There is another layer to the lithium story. HD 81809B’s observed lithium abundance shows the star is still relatively lithium-rich, which is unusual. For a low-mass main sequence star, lithium should have been depleted during the main sequence. Planet engulfment, if it delivered even a modest amount of lithium-bearing material at a recent epoch, is one of the few mechanisms that could explain this retention.

Looking for the forensic signature

The authors identify a concrete observational path forward. Detecting rotation and magnetic activity on HD 81809B, which could hold telltale signs of a planet engulfment event, would provide additional evidence. A recent engulfment is expected to spin up the accreting star through the transfer of orbital angular momentum from the swallowed body, and to reactivate its magnetic dynamo — leaving independent traces on the stellar surface that complement the chemical signature.

This study emphasizes the importance of considering such external events when modeling chemically anomalous binary systems. The chemical composition of the accreted material may differ from the assumptions adopted here. That caveat is a call for more refined observations: high-resolution spectroscopy sensitive to beryllium and neutron-capture element abundances, continued activity monitoring, and asteroseismic constraints on interior mixing profiles would all help close the remaining gaps.

HD 81809 is not giving up its answer easily. That is precisely what makes it worth studying.

Source: Moedas, N. & Di Mauro, M. P. (2026). Chemical paradox in a binary system: Exploring metal enrichment in HD 81809B. arXiv:2605.31060 [astro-ph.SR]. DOI: 10.48550/arxiv.2605.31060