REVIEW 2 major objections 5 minor 34 references
Planet engulfment can explain TOI-5882's lithium, but only under uncommon timing and orbits.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · grok-4.5
2026-07-11 05:42 UTC pith:XQGXNT76
load-bearing objection Solid system-specific N-body+MESA test: engulfment is common overall but rare (~5%) inside the Li window under a frozen BD orbit, so the Kotten scenario stays viable but non-generic. the 2 major comments →
Dynamical constraints on planet engulfment as the origin of lithium enhancement in TOI-5882
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Across 300 simulations that span the mass and period ranges suggested for a pre-engulfment planet, planet engulfment is the most probable outcome overall, yet only about 5 percent of the runs produce engulfment inside the short lithium-detectability window (roughly 4.23–4.61 Gyr). Thus recent engulfment remains a viable but uncommon explanation for the observed lithium enhancement under the simplified architecture considered.
What carries the argument
MESA stellar tracks interpolated into REBOUND N-body integrations (with REBOUNDx tides and time-dependent stellar mass, radius and Love number) that evolve an inner planet plus the observed brown dwarf while the star leaves the main sequence; outcomes are timed against the narrow convective-envelope window in which accreted lithium would remain detectable.
Load-bearing premise
The brown dwarf is fixed at its currently observed short-period orbit and wind drag is omitted, so the reported 5 percent fraction is conditional on that frozen architecture and force model.
What would settle it
Additional chemical tracers of planetary material, or tighter constraints on the brown dwarf's earlier orbit and the system's primordial architecture, that either place a recent engulfment inside the lithium window or show that no such planet could have survived long enough to be swallowed then.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper tests whether the lithium enhancement of the subgiant TOI-5882 can be explained by recent engulfment of a super-Earth to Neptune-mass planet, as suggested by Kotten et al. (2026). Combining a MESA track for a 1.3 M⊙, [Fe/H]=0.33 star with REBOUND/REBOUNDx N-body integrations (IAS15, equilibrium-tide constant-time-lag model with envelope-dependent Love number), the authors run 300 simulations of a system containing the observed 22 m_Jup brown dwarf (P=7.1 d, e=0.0339) plus a randomly drawn interior planet (mass 9–95 m⊕, period 0.5–3.5 d). They define a lithium-detectability window from the stellar models (4.23–4.61 Gyr) during which accreted material can produce an observable surface A(Li) excess. Planet engulfment is the most common outcome overall (~65%), but only ~5% of runs produce engulfment inside that window; the authors therefore conclude that recent engulfment remains viable but non-generic under the simplified single-interior-planet architecture explored.
Significance. The work supplies a concrete dynamical test of an externally proposed chemical scenario for a specific system, rather than a generic population statement. Strengths include a carefully defined MESA-based detectability window, self-consistent coupling of time-dependent stellar mass/radius/luminosity and tides into the N-body runs, and an explicit scoping of the ~5% figure as a fraction of explored initial conditions rather than a posterior probability for TOI-5882. The result is useful for the growing literature on planet engulfment and Li-rich subgiants: it shows that timing, not merely the occurrence of engulfment, is the binding constraint, and it flags the need for better primordial-architecture constraints and additional forces (e.g., wind drag). The caveats on frozen BD orbit and omitted physics are already stated in the Discussion, which keeps the claim appropriately limited.
major comments (2)
- Discussion (and Methods setup): The integrations initialize the brown dwarf at its currently observed post-MS orbit and omit both secular main-sequence tidal decay and wind-induced drag. The paper itself cites Narayan et al. (2026) for continued orbital decay on ~25–125 Myr timescales and notes that wind drag was not included. Because the Li-window fraction is the central quantitative result, the manuscript should either (i) add a small set of exploratory runs with a modestly wider/closer initial BD period (or a simple drag term) to show how the ~5% figure moves, or (ii) elevate the existing qualitative caveat into a quantitative bound (e.g., order-of-magnitude estimate of how much the window fraction could change). Without that, the load-bearing claim remains conditional on an architecture that the authors acknowledge is transient.
- Results / Fig. 3 and Appendix B: The reported ~5% (and ~65%) fractions are raw counts over a uniform random draw in mass and period that deliberately spans both Hill-unstable and Hill-stable regimes. The paper correctly states that this is not a probability for TOI-5882, yet the abstract and conclusions still present “only 5%” as the headline number. A short reweighting or stratified reporting (e.g., fraction among runs that remain stable until the subgiant branch, or among runs with Δ_H ≳ 5) would make the non-generic character of successful cases clearer and less sensitive to the arbitrary prior volume. This is a presentation-of-result issue that affects how the central claim is read, not a flaw in the integrations themselves.
minor comments (5)
- Fig. 1 bottom panel and Appendix A: The post-engulfment A(Li) track is described as a post-processing instantaneous-mixing model. A one-sentence statement of the assumed accreted Li mass (or reference to the Kotten et al. budget used) would make the green detectability cutoff fully reproducible.
- Fig. 3 caption: The blue and green dashed stability lines are defined in Appendix B; a brief parenthetical in the caption (Δ_H ~ 2√3 and Δ_H = 5) would help readers who do not immediately turn to the appendix.
- Methods: The stellar tidal Love number is said to follow Becker & Batygin (2013) and to vary with the convective envelope. A short note on whether the planet and BD Love numbers/tidal quality factors are held fixed would complete the force-model description.
- Throughout: Occasional typographical inconsistencies (e.g., “offthe”, “m Jup” spacing, “V owell”) should be cleaned in production.
- Conclusions: The suggestion to measure additional chemical tracers is valuable; naming one or two concrete species (e.g., Be, refractory-to-volatile ratios) would strengthen the observational outlook paragraph.
Circularity Check
No significant circularity: the ~5% Li-window engulfment fraction is an independent output of the N-body ensemble, not forced by definition, fit, or self-citation.
full rationale
The paper tests an external hypothesis (Kotten et al. 2026 Li budget and suggested super-Earth–Neptune engulfment) by combining MESA stellar tracks with 300 REBOUND N-body runs that sample planet mass/period/eccentricity/inclination ranges motivated by that hypothesis. The detectability window (4.23–4.61 Gyr) is defined from the models’ convective-envelope onset and A(Li) dilution/burning curves; the reported ~65% overall engulfment and <5% in-window fractions are direct counts of simulation outcomes under those fixed initial conditions. No parameter is fitted to the target Li signature and then re-presented as a prediction; no uniqueness theorem or ansatz is imported from the authors’ prior work to force the result; self-citations (Aguilera-Gómez et al. 2016, 2020) supply only background Li-abundance context and are not load-bearing for the dynamical fractions. The Discussion explicitly scopes the claim as conditional on the simplified single-interior-planet architecture and the frozen post-MS BD orbit, confirming the derivation is self-contained numerical exploration rather than circular redefinition.
Axiom & Free-Parameter Ledger
free parameters (5)
- Inner-planet mass sampling range =
[9, 95] m_earth
- Inner-planet period sampling range =
[0.5, 3.5] days
- Planet eccentricity and inclination priors =
e~1e-4–1e-3; |i|≤1°
- MESA mixing-length, overshooting, Reimers efficiency =
α_MLT=1.73; f=0.016; η=0.1
- Stellar tidal Love number prescription =
Becker & Batygin (2013) envelope-dependent k
axioms (5)
- domain assumption Equilibrium-tide constant-time-lag model (Hut 1981) adequately captures post-MS tidal evolution for both planet and BD.
- domain assumption Li delivered by an engulfed 9–95 m⊕ planet is instantaneously mixed in the surface convective envelope and then destroyed under the same burning that sets the baseline A(Li) track.
- ad hoc to paper The currently observed BD orbit can be used as the initial post-MS condition for the ensemble (neglecting prior MS tidal decay and wind drag).
- domain assumption Mutual Hill-radius separations (Gladman/Chambers/Obertas-type criteria) usefully organize stable vs unstable initial placements for this two-companion system.
- ad hoc to paper A single interior planet (not multi-planet or exterior architectures) is a sufficient test of dynamical viability for the Li scenario.
invented entities (1)
-
Hypothetical pre-engulfment super-Earth to Neptune-mass planet
no independent evidence
read the original abstract
As stars evolve off the main sequence, changes in stellar structure can alter the dynamical architecture of planetary systems and, in some cases, lead to planet engulfment events capable of producing observable chemical signatures such as lithium enhancement. TOI-5882 is a lithium-rich subgiant hosting a 22 m_Jup brown dwarf on a 7.1-day orbit, where the enrichment could plausibly result from the recent engulfment of a super-Earth to Neptune-mass planet. We assess the dynamical viability of a planet engulfment scenario in the TOI-5882 system that could explain its observed Li enrichment. We combine stellar evolution models with N-body simulations that incorporate time-dependent stellar properties as the host star leaves the main sequence, exploring a broad range of pre-engulfment planetary masses and orbital configurations. Planet engulfment is the most probable outcome under the explored configurations. However, only 5% of the simulations produce engulfment within the short detectability window of the lithium-enrichment signature. Engulfment therefore remains a viable explanation for the observed Li enhancement in TOI-5882, but only for a relatively uncommon subset of the initial conditions considered here. Under our simplified architecture, recent engulfment cannot be ruled out as the origin of the observed lithium enhancement. However, because successful cases represent only a small fraction of the initial conditions explored here, and because additional processes such as wind-induced drag forces were not included, additional constraints on the system's primordial architecture are needed to assess how likely this scenario is for TOI-5882.
Figures
Reference graph
Works this paper leans on
-
[1]
The resonance overlap criterion and the onset of stochastic behavior in the restricted three-body problem. , keywords =. doi:10.1086/112778 , adsurl =
-
[2]
, keywords =
Tidal evolution in close binary systems. , keywords =
-
[3]
Dynamics of Systems of Two Close Planets. , keywords =. doi:10.1006/icar.1993.1169 , adsurl =
-
[4]
, year = 1996, month = feb, volume =
The Stability of Multi-Planet Systems. , year = 1996, month = feb, volume =. doi:10.1006/icar.1996.0019 , adsurl =
-
[5]
Are There Unstable Planetary Systems around White Dwarfs?. , keywords =. doi:10.1086/340291 , archivePrefix =. astro-ph/0202273 , primaryClass =
-
[6]
Can Planets Survive Stellar Evolution?. , keywords =. doi:10.1086/516746 , archivePrefix =. astro-ph/0702724 , primaryClass =
-
[7]
Dynamical Outcomes of Planet-Planet Scattering. , keywords =. doi:10.1086/590227 , archivePrefix =. astro-ph/0703166 , primaryClass =
-
[8]
The Orbital Evolution of Gas Giant Planets Around Giant Stars. , keywords =. doi:10.1088/0004-637X/705/1/L81 , archivePrefix =. 0910.2396 , primaryClass =
-
[9]
Planet Engulfment by ~1.5-3 Solar-Mass Red Giants
Planet Engulfment by -0.5ex 1.5-3 M _ sun Red Giants. , keywords =. doi:10.1088/0004-637X/737/2/66 , archivePrefix =. 1106.2251 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1088/0004-637x/737/2/66
-
[10]
Modules for Experiments in Stellar Astrophysics (MESA). , keywords =. doi:10.1088/0067-0049/192/1/3 , archivePrefix =. 1009.1622 , primaryClass =
-
[11]
Foretellings of Ragnar \"o k: World-engulfing Asymptotic Giants and the Inheritance of White Dwarfs. , keywords =. doi:10.1088/0004-637X/761/2/121 , archivePrefix =. 1210.0328 , primaryClass =
-
[12]
REBOUND: an open-source multi-purpose N-body code for collisional dynamics. , keywords =. doi:10.1051/0004-6361/201118085 , archivePrefix =. 1110.4876 , primaryClass =
-
[13]
Dynamical Measurements of the Interior Structure of Exoplanets
Dynamical Measurements of the Interior Structure of Exoplanets. , keywords =. doi:10.1088/0004-637X/778/2/100 , archivePrefix =. 1309.5363 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1088/0004-637x/778/2/100
-
[14]
First order resonance overlap and the stability of close two planet systems
First-order Resonance Overlap and the Stability of Close Two-planet Systems. , keywords =. doi:10.1088/0004-637X/774/2/129 , archivePrefix =. 1307.8119 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1088/0004-637x/774/2/129
-
[15]
The Resonance Overlap and Hill Stability Criteria Revisited
The resonance overlap and Hill stability criteria revisited. Celestial Mechanics and Dynamical Astronomy , keywords =. doi:10.1007/s10569-015-9646-z , archivePrefix =. 1509.03607 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1007/s10569-015-9646-z
-
[16]
IAS15: a fast, adaptive, high-order integrator for gravitational dynamics, accurate to machine precision over a billion orbits. , keywords =. doi:10.1093/mnras/stu2164 , archivePrefix =. 1409.4779 , primaryClass =
-
[17]
WHFAST: a fast and unbiased implementation of a symplectic Wisdom-Holman integrator for long-term gravitational simulations. , keywords =. doi:10.1093/mnras/stv1257 , archivePrefix =. 1506.01084 , primaryClass =
-
[18]
Post-main-sequence planetary system evolution
Post-main-sequence planetary system evolution. Royal Society Open Science , keywords =. doi:10.1098/rsos.150571 , archivePrefix =. 1601.05419 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1098/rsos.150571
-
[19]
On Lithium-6 as diagnostic of the lithium-enrichment mechanism in red giants
On Lithium-6 as a Diagnostic of the Lithium-enrichment Mechanism in Red Giants. , keywords =. doi:10.3847/2041-8213/ab9d26 , archivePrefix =. 2005.14209 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/2041-8213/ab9d26 2041
-
[20]
How Jupiters save or destroy inner Neptunes around evolved stars
How Jupiters Save or Destroy Inner Neptunes around Evolved Stars. , keywords =. doi:10.3847/2041-8213/aba35f , archivePrefix =. 2007.04337 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/2041-8213/aba35f 2041
-
[21]
The GALAH survey: Data release 4. , keywords =. doi:10.1017/pasa.2025.26 , archivePrefix =. 2409.19858 , primaryClass =
-
[22]
Eleven New Transiting Brown Dwarfs and Very-low-mass Stars from TESS. , keywords =. doi:10.3847/1538-3881/addd17 , archivePrefix =. 2501.09795 , primaryClass =
-
[23]
Observable Signatures of Planet Accretion in Red Giant Stars. I. Rapid Rotation and Light Element Replenishment. , keywords =. doi:10.1088/0004-637X/757/2/109 , archivePrefix =. 1208.1775 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1088/0004-637x/757/2/109
-
[24]
Hot Jupiters and Cool Stars. , keywords =. doi:10.1088/0004-637X/794/1/3 , archivePrefix =. 1407.7879 , primaryClass =
-
[25]
Star-planet interactions. II. Is planet engulfment the origin of fast rotating red giants?. , keywords =. doi:10.1051/0004-6361/201628758 , archivePrefix =. 1606.08027 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1051/0004-6361/201628758
-
[26]
On Lithium-rich Red Giants. I. Engulfment of Substellar Companions. , keywords =. doi:10.3847/0004-637X/829/2/127 , archivePrefix =. 1605.05332 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/0004-637x/829/2/127
-
[27]
The stability of tightly-packed, evenly-spaced systems of Earth-mass planets orbiting a Sun-like star. , keywords =. doi:10.1016/j.icarus.2017.04.010 , archivePrefix =. 1703.08426 , primaryClass =
-
[28]
REBOUNDx: a library for adding conservative and dissipative forces to otherwise symplectic N-body integrations. , keywords =. doi:10.1093/mnras/stz2870 , archivePrefix =. 1908.05634 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1093/mnras/stz2870 1908
-
[29]
Lithium Enrichment Signatures of Planetary Engulfment Events in Evolved Stars
Lithium Enrichment Signatures of Planetary Engulfment Events in Evolved Stars. , keywords =. doi:10.3847/1538-3881/ac273c , archivePrefix =. 2002.05275 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/1538-3881/ac273c 2002
-
[30]
Stellar Evolution and Tidal Dissipation in REBOUNDx
Stellar evolution and tidal dissipation in REBOUNDx. , keywords =. doi:10.1093/mnras/stac043 , archivePrefix =. 2101.12277 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1093/mnras/stac043
-
[31]
Planet Engulfment Detections are Rare According to Observations and Stellar Modeling
Planet engulfment detections are rare according to observations and stellar modelling. , keywords =. doi:10.1093/mnras/stad745 , archivePrefix =. 2210.12121 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1093/mnras/stad745
-
[32]
Planet Engulfment Signatures in Twin Stars
Planet engulfment signatures in twin stars. , keywords =. doi:10.1093/mnras/stac3435 , archivePrefix =. 2210.11679 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1093/mnras/stac3435
-
[33]
arXiv e-prints , keywords =
Lithium Enrichment in a Subgiant Star with a Brown Dwarf Companion: A Planetary Engulfment Candidate. arXiv e-prints , keywords =
-
[34]
The Tale of a Hungry Subgiant and Its Brown Dwarf: Interior Radiative Damping Dominates the Tidal Evolution of TOI-5882. arXiv e-prints , keywords =. doi:10.48550/arXiv.2605.04141 , archivePrefix =. 2605.04141 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2605.04141
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.