{"id":"c396f19f-2597-44df-a727-0113838be1a6","arxiv_id":"2506.20360","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Applying internal-gravity-wave mixing to intermediate-mass AGB stars raises their predicted lithium yields, giving these stars about 10 percent of Galactic lithium.","lead":"This paper models how internal gravity waves, stirred up by convection inside medium-mass dying stars, can carry beryllium-7 into cooler outer layers where it becomes lithium-7, boosting the lithium these stars release. It estimates that these stars produce about 15 solar masses of lithium, roughly twice previous AGB estimates and about 10 percent of the Galaxy's lithium.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. 6 sets L_IGW = M L_conv, but the cited IGW literature (Lecoanet & Quataert 2013; Schwab 2020) uses L_IGW ∝ M^3 L_conv; if so, the He-shell IGW luminosity falls below the 10^4 L_sun mixing threshold and the extra-mixing mechanism cannot operate.","rationale":"The reader's weakest_assumption correctly flags the IGW prescriptions, but treats them as 'possibly overestimated.' My stress-test goes one step further: Eq. (6) appears to misquote the wave-luminosity scaling of the very references it cites. The central claim is not that IGW mixing is somewhat weaker; it is that the wave power budget may be off by orders of magnitude, in which case the radiative gap is never mixed and the Li doubling vanishes. This is the single most load-bearing point because every quantitative downstream result—the 7Be transport, the yield grid in Table A.1, the 15 M_sun total yield, and the 10% Galactic contribution—inherits it. Other concerns (arithmetic inconsistencies in Table A.1, missing data, no error bars) are important for reproducibility but do not by themselves invalidate the physics. A mis-stated Eq. (6) would. The authors do not report the convective Mach number M or the individual terms of Eq. (4), so the test cannot be done from the manuscript alone; that is precisely why the answer should be CONDITIONAL. I agree partially with the reader: same region of vulnerability, but a sharper and more specific version. I keep the verdict UNCHANGED because the manuscript already received CONDITIONAL, and this concern tightens the condition without requiring a different verdict. If the test fails, the paper should be rejected; if it passes, the extra-mixing mechanism remains plausible.","tokens_in":15523,"tokens_out":9102,"duration_ms":91698,"concrete_test":"Re-derive Eq. (6) from Lecoanet & Quataert (2013) and Schwab (2020). If the correct scaling is L_IGW = M^3 L_conv, extract the convective Mach number M in the He-shell flash convection zone of the 4 M_sun, Z=0.014 model at peak L_He, compute L_IGW,c = M^3 L_conv, and compare with L_mix ≈ 10^4 L_sun from Eq. (4). If L_IGW,c < L_mix, the activation condition fails and the central mechanism does not operate. Then rerun one fiducial model with the corrected L_IGW and check whether 7Be is still transported into the convective envelope and whether the Li yield changes.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim—IGW mixing during each He-shell flash transports 7Be into the convective envelope, doubling the Li yield—rests entirely on L_IGW being large enough to sustain mixing across the radiative gap. Eq. (6) states L_IGW = M L_conv, citing Lecoanet & Quataert (2013) and Schwab (2020). In those references the gravity-wave luminosity generated by convection is L_IGW ~ M^3 L_conv, not M L_conv. The M^3 scaling reflects that wave excitation depends on the convective velocity relative to the sound speed to a higher power, and Schwab (2020) uses it for the RGB He-flash. If Eq. (6) is a mis-transcription, the values shown in Fig. 2 ('10^4-10^5 L_sun') are too high by a factor ~M^-2. For He-shell-flash convection M is typically ≤0.1, so the corrected L_IGW is likely at or below L_mix ≈ 10^4 L_sun from Eq. (4). The L_He > 10^4 L_sun activation condition does not repair this: that condition concerns He luminosity, not wave luminosity; the comparison that matters is L_IGW versus L_mix. If the corrected wave luminosity is sub-threshold, no 7Be transport channel opens and the yield enhancement disappears. This is an apparent factor-10^2-10^4 error in the wave power budget, not a calibration nuance. Eq. (7) with η=0.1 is secondary if the power budget fails.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper uses MESA stellar models to study internal gravity wave (IGW) mixing in intermediate-mass TP-AGB stars and its effect on lithium yields. The authors propose that IGWs excited by He-shell flash convection mix the radiative zone between the convective thermal pulse and the convective envelope, transporting 7Be into the envelope where it decays to 7Li. They construct grids with initial masses 3.5-7.5 Msun and four metallicities, compute Li yields with and without IGW mixing, and use population synthesis with a Kroupa IMF to estimate a total Galactic Li yield of about 15 Msun, twice the no-IGW value and roughly 10% of the present-day Li content of the Milky Way.","tokens_in":15935,"tokens_out":14601,"duration_ms":172782,"significance":"If the central mechanism is correct, the paper offers a concrete physical explanation for extra-mixing in intermediate-mass AGB stars and identifies IGW mixing as a potentially non-negligible source of Li in the Galaxy. The work uses a standard stellar evolution code, provides a yield grid for multiple masses and metallicities in Table A.1, and compares the models with Li-rich AGB stars in the Milky Way and the Magellanic Clouds. The main limitation is that the IGW luminosity, activation threshold, and mixing coefficient are imported without AGB-specific calibration or sensitivity analysis, and the Galactic yield estimate is based on a solar-metallicity-only population synthesis. The quantitative conclusions therefore rest on prescriptions that are not yet established for this evolutionary phase.","major_comments":[{"comment":"The wave luminosity is written as L_IGW = M L_conv, where M is the convective Mach number, and this is attributed to Lecoanet & Quataert (2013) and Schwab (2020). In those references the scaling is instead L_IGW ~ M^3 L_conv, which differs by a factor M^2 relative to the expression used here. For the He-shell convection in these models, M is likely at most ~0.1, so the IGW luminosity used in Fig. 2 may be overestimated by one to four orders of magnitude, placing it at or below the L_mix ~ 10^4 L_sun threshold derived from Eq. (4). Because the entire extra-mixing mechanism and the factor-of-two Li yield enhancement depend on this power budget, the authors must verify the formula against the cited papers, report the correct scaling, and recompute the models if needed.","section":"§2.2.1, Eq. (6)"},{"comment":"The IGW mixing is activated by the condition L_He > 10^4 L_sun, but the physical criterion established in Section 2.2.1 is L_IGW > L_mix, where L_IGW depends on the convective luminosity and Mach number rather than directly on the helium-burning luminosity. The relationship between L_He and L_conv during a He-shell flash is not established in the text, so the adopted threshold is an ad hoc proxy. The model should either switch mixing based on L_IGW versus L_mix or justify why L_He is an equivalent proxy; otherwise the activation of extra mixing may be overestimated.","section":"§2.2.2"},{"comment":"The mixing coefficient D_IGW ~ eta (curl u)^2 K/N^2 is adopted from Herwig (2023) with eta=0.1 determined by Garaud & Kulenthirarajah (2016). This coefficient is not calibrated for the He-shell flash regime in AGB stars, where the radiative gap is transient, N^2 is large, and the geometry may differ from the simulations used to calibrate eta. No sensitivity study is presented: varying eta by even a factor of a few could change whether 7Be is transported to the envelope before decaying. A robustness test of the Li yields against eta and against the activation threshold should be included, given that the factor-of-two yield enhancement is the central claim.","section":"§2.2.2, Eq. (7)"},{"comment":"The total Galactic Li yield of about 15 Msun is obtained by linear interpolation in mass only over the Z=0.014 grid in Table A.1, with no metallicity distribution integrated over the star formation history. This is not a Galactic yield estimate because Table A.1 shows strong metallicity dependence, including cases where IGW mixing reduces the yield (for example Z=0.00014, 6.0-7.5 Msun). The authors should either integrate over a metallicity distribution appropriate for the Milky Way or explicitly argue that solar-metallicity stars dominate the Li production. As written, the 10% Galactic contribution claim is not supported by the calculation shown.","section":"§3.3, population synthesis"},{"comment":"The mechanism requires 7Be to be transported across the radiative gap and into the convective envelope in less than the 53-day decay half-life. The paper does not compare the IGW diffusion timescale with the 7Be lifetime. For the lowest D_IGW values quoted as ~10^10 cm^2/s and a gap thickness that can be inferred from Fig. 3, the transport time may be months or longer, which would mean much of the 7Be decays before reaching the envelope. The authors should report the relevant diffusion timescales across the mixing zone to show that the 7Be actually survives the transport.","section":"§3.1, Fig. 3"}],"minor_comments":[{"comment":"There are numerous typographical and naming errors, including 'Pranztos' for Prantzos, 'Unclear network' for 'Nuclear network', 'Venture' or 'Venture 2001' for Ventura, 'tempature' for temperature, 'enhacement' for enhancement, 'Ruckya' for Rukeya, and 'the the' in the text. These should be corrected.","section":"Throughout"},{"comment":"The quantity L in the scaling |curl u| proportional to L^{1/3} is not defined; the text should state whether this is the local luminosity, convective luminosity, or another quantity, and in which region the scaling is applied.","section":"Eq. (8)"},{"comment":"The statement that both the IGW and no-IGW models cover the observational samples equally well weakens the case for IGW mixing; the discussion should either identify a regime where observations can discriminate between the two sets of models or state explicitly that the comparison is not a discriminating test.","section":"Fig. 7"},{"comment":"Several entries in Table A.1 show that models without IGW mixing have higher yields than those with IGW mixing at low metallicity and high mass. This caveat appears in the text but is absent from the abstract and conclusions, where the effect of IGW is described only as positive; the abstract should be qualified accordingly.","section":"Table A.1"}],"recommendation":"major_revision","confidential_remarks":"For the editor: the key risk in this paper is the IGW luminosity scaling in Eq. (6). If the cited references do indeed use a M^3 scaling, the central mechanism may not operate and the paper would require a fundamental reanalysis. I did not see this explicitly resolved in the manuscript, so the authors should be asked to verify the formula and rerun the models if necessary. The population-synthesis estimate is also too simplistic to support the 10% Galactic contribution claim. These issues are load-bearing but appear fixable in principle, hence major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this paper. It applies the IGW extra-mixing picture from Schwab's RGB He-flash work to intermediate-mass TP-AGB stars, computes Li yields in MESA across a mass/metallicity grid, and ends with a population-synthesis number: about 15 Msun of Li, roughly 10% of the Galactic inventory, double the no-IGW case. That quantitative target is new, and the observational comparison with the Milky Way and Magellanic Cloud Li-rich AGB stars is a genuine plus: they match the A(Li)-T_eff samples without invoking the artificially high mass-loss rates that earlier yield calculations needed. The paper deserved a careful read.\n\nThe soft spots are real, and one is load-bearing. Eq. (6) states L_IGW = M L_conv, citing Lecoanet & Quataert (2013) and Schwab (2020). The cited works use L_IGW proportional to M^3 L_conv. That is not a calibration nicety: it is two powers of M. For He-shell flash convection M is typically at or below 0.1, so the wave luminosity plotted in Fig. 2 (10^4 to 10^5 Lsun) drops by about 10^2 to 10^4, below the L_mix near 10^4 Lsun threshold from Eq. (4). If that is a mis-transcription, the 7Be transport channel never opens and the doubling of the Li yield disappears. The L_He > 10^4 Lsun activation condition does not save it, because the relevant comparison is L_IGW versus L_mix, not the helium luminosity.\n\nThere are also smaller but annoying issues. A few rows in Table A.1 do not subtract correctly: for Z=0.014, M=4.5, the last column says 1.2e-9 while 1.9e-8 minus 6.6e-9 is 1.2e-8; M=6.5 is off too. No inlists or data are released, and there are no error bars on any of the yield numbers, which frustrates a grid that is the basis of the Galactic claim. The non-IGW models already give 8 Msun, ten times the older estimates, so part of the 10% contribution comes from the updated Simonucci 7Be e-capture rate, not from IGW; the IGW-specific increment is the 8 to 15 Msun gap.\n\nI do not think this is a circular or unserious paper. The reaction network, mixing coefficient framework, and population synthesis are honest, and the authors cite the relevant work. But the central quantitative claim currently rests on an equation that appears inconsistent with its own citations. That is exactly what a referee should be able to settle. Send it to review; ask the authors to rederive Eq. (6), rerun the models if needed, and fix the table. If the corrected power budget is sub-threshold, the paper becomes a null result. If the M^3 scaling still clears L_mix for part of the grid, the conclusion survives and this is a useful contribution to Galactic Li evolution.","headline":"The paper's new numbers deserve refereeing, but Eq. (6) misstates the wave-luminosity scaling from its own citations, and if that holds up the central mixing mechanism loses its power budget.","tokens_in":16473,"tokens_out":6038,"would_cite":false,"duration_ms":57874,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Internal gravity waves excited by helium-shell flashes mix 7Be into AGB convective envelopes, doubling the predicted lithium yield from intermediate-mass AGB stars to roughly 15 solar masses.","keywords":["internal gravity waves","asymptotic giant branch stars","lithium production","hot bottom burning","extra mixing","stellar evolution","galactic lithium abundance","MESA stellar models"],"falsifier":"A three-dimensional hydrodynamic simulation of a helium-shell flash in a roughly 4 solar mass AGB star could directly measure the wave luminosity and the turbulent diffusivity in the radiative gap; if the computed $L_{\\rm IGW}$ stays below the approximately $10^4\\,L_\\odot$ needed to sustain mixing, or the effective $D_{\\rm IGW}$ falls below $10^{10}\\,\\mathrm{cm^2\\,s^{-1}}$ during the flash, the claimed 15 solar mass Galactic lithium yield would not follow.","tokens_in":15322,"feed_emoji":"🔭","tokens_out":6243,"duration_ms":60529,"temperature":0.7,"pith_summary":"The paper argues that internal gravity waves (IGW) excited during each helium-shell flash of intermediate-mass asymptotic giant branch (AGB) stars mix the radiative zone between the thermal pulse and the convective envelope, carrying 7Be into the envelope where it decays to 7Li. With this extra mixing included, the surface lithium abundance can exceed A(Li)=5, and most stars of 3.5 to 7.5 solar masses end up as net lithium producers. Population synthesis gives a total Galactic lithium yield from these stars of about 15 solar masses, twice the yield without IGW mixing and roughly 10% of the lithium in the Galaxy's interstellar medium. The paper therefore claims that AGB stars are a non-negligible lithium source under normal mass-loss rates, removing the need for artificially high mass loss that earlier models required.","feed_headline":"Mixing by internal gravity waves doubles lithium yield of AGB stars","feed_subtitle":"Simulations put AGB stars' contribution to Galactic lithium at about 10%, twice previous estimates.","key_machinery":"The load-bearing object is the IGW-driven diffusion coefficient together with its activation threshold. The paper uses $D_{\\rm IGW}\\approx \\eta\\,(\\nabla\\times u)^2 K/N^2$ with $\\eta=0.1$ taken from three-dimensional hydrodynamic simulations, and activates the mixing only when the helium-shell luminosity exceeds $10^4\\,L_\\odot$. The wave luminosity $L_{\\rm IGW}=M\\,L_{\\rm conv}$ is compared against the mixing-maintenance power $L_{\\rm mix}=\\int N^2 K\\,dm$; this comparison is what allows the radiative gap between the thermal pulse and the convective envelope to be treated as mixed, creating the 7Be transport channel that roughly doubles the predicted lithium yield.","core_discovery":"The central claim is that internal gravity waves are the missing extra-mixing agent in thermally pulsing AGB stars. During helium-shell ignition, the wave luminosity from convective motions exceeds the luminosity needed to sustain mixing across the radiative gap ($L_{\\rm IGW}=M\\,L_{\\rm conv} > L_{\\rm mix}\\sim10^4\\,L_\\odot$), so the diffusion coefficient $D_{\\rm IGW}\\approx \\eta\\,(\\nabla\\times u)^2 K/N^2$ stays above $10^{10}\\,\\mathrm{cm^2\\,s^{-1}}$ in that gap. This transports 7Be produced near $4\\!-\\!8\\times10^7\\,\\mathrm{K}$ into the cooler convective envelope before it is destroyed, realizing the Cameron\\,-\\,Fowler mechanism. The simulations yield peak A(Li)>5, with the largest relative lithium gain in stars just massive enough for hot bottom burning, and total yields of about 15 solar masses versus about 8 solar masses without IGW mixing.","pith_inferences":["If the same IGW prescription operates during the first helium flash in low-mass red giants, then lithium-rich giants and lithium-rich AGB stars may share one underlying transport mechanism; the paper notes the RGB analogue but does not quantify its interstellar-medium contribution.","A direct testable extension would be a three-dimensional hydrodynamic simulation of a helium-shell flash that computes the wave luminosity from first principles; if it falls below the $\\sim10^4\\,L_\\odot$ threshold, the claimed yield enhancement would not occur.","The yield grid treats only single stars; binarity or rotation could modify the mixing and change the 15 solar mass figure, an effect this paper does not address."],"forward_implications":["IGW mixing makes intermediate-mass AGB stars net lithium producers across most of the 3.5 to 7.5 solar mass range, with per-star yields up to about $10^{-7}\\,M_\\odot$.","The total Galactic lithium contribution from these stars rises to about 15 solar masses, roughly 10% of the interstellar lithium inventory, an order of magnitude above earlier AGB estimates.","Li-rich, O-rich AGB stars with A(Li) up to roughly 4.3 and mass-loss rates below $10^{-6}\\,M_\\odot\\,\\mathrm{yr^{-1}}$ can be reproduced without invoking anomalous mass loss.","The relative benefit of IGW mixing is largest near the minimum mass for hot bottom burning and declines toward higher masses, so lithium-rich AGB stars should preferentially be found at the low-mass end of the hot-bottom-burning range."],"supporting_citations":[{"why":"Supplies the wave-luminosity relation $L_{\\rm IGW}=M\\,L_{\\rm conv}$ that powers the IGW mixing.","marker":"Lecoanet & Quataert (2013)"},{"why":"Applies IGW mixing to the first helium flash in low-mass stars and sets the $L_{\\rm He}>10^4\\,L_\\odot$ activation condition adopted here.","marker":"Schwab (2020)"},{"why":"Provides the $L_{\\rm mix}=\\int N^2 K\\,dm$ criterion for the power needed to sustain a mixed region.","marker":"Press (1981)"},{"why":"Provides the $D_{\\rm IGW}$ mixing-coefficient framework used to model IGW-induced extra mixing.","marker":"Herwig et al. (2023)"},{"why":"Calibrates the dimensionless efficiency $\\eta=0.1$ in the $D_{\\rm IGW}$ formula from three-dimensional hydrodynamic simulations.","marker":"Garaud & Kulenthirarajah (2016)"},{"why":"Supplies earlier AGB lithium models with a fixed mixed mass; the authors argue IGW mixing is the physical origin of that assumed extra mixing.","marker":"Karakas & Lugaro (2016)"},{"why":"One of the previous AGB lithium yield calculations that the present 15 solar mass result exceeds by an order of magnitude.","marker":"Romano et al. (2001)"},{"why":"The other previous AGB lithium yield calculation used for comparison in the yield grid.","marker":"Ventura et al. (2000)"}],"fun_headline_variants":["Gravity-wave mixing doubles lithium yield of AGB stars","Internal gravity waves lift AGB lithium output twofold","Wave-driven mixing doubles lithium from intermediate AGB stars","AGB stars double lithium production via internal waves","Gravity waves boost AGB lithium yield by a factor of two"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that helium-shell flash convection generates enough internal gravity wave power, and the adopted wave-mixing coefficient is strong enough, to mix the entire radiative gap between the thermal pulse and the convective envelope; if the waves are weaker or the coefficient is smaller, the extra lithium production mostly disappears.","fun_headline_variants_meta":{"raw":{"variants":["Gravity-wave mixing doubles lithium yield of AGB stars","Internal gravity waves lift AGB lithium output twofold","Wave-driven mixing doubles lithium from intermediate AGB stars","AGB stars double lithium production via internal waves","Gravity waves boost AGB lithium yield by a factor of two"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000526,"raw_usage":{"total_tokens":2586,"prompt_tokens":1040,"completion_tokens":1546,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":656,"completion_tokens_details":{"reasoning_tokens":1468}},"tokens_in":656,"tokens_out":1546,"duration_ms":13581,"temperature":1.0,"reasoning_tokens":1468,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T22:50:18.832308+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A three-dimensional hydrodynamic simulation of a helium-shell flash in a roughly 4 solar mass AGB star could directly measure the wave luminosity and the turbulent diffusivity in the radiative gap; if the computed $L_{\\rm IGW}$ stays below the approximately $10^4\\,L_\\odot$ needed to sustain mixing, or the effective $D_{\\rm IGW}$ falls below $10^{10}\\,\\mathrm{cm^2\\,s^{-1}}$ during the flash, the claimed 15 solar mass Galactic lithium yield would not follow.","supporting_citations":[{"cited_title":"2020, ApJ, 901, L18","cited_arxiv_id":null,"evidence_quote":"Applies IGW mixing to the first helium flash in low-mass stars and sets the $L_{\\rm He}>10^4\\,L_\\odot$ activation condition adopted here."},{"cited_title":"R., Mao, H., et al","cited_arxiv_id":null,"evidence_quote":"Provides the $D_{\\rm IGW}$ mixing-coefficient framework used to model IGW-induced extra mixing."},{"cited_title":"& Kulenthirarajah, L","cited_arxiv_id":null,"evidence_quote":"Calibrates the dimensionless efficiency $\\eta=0.1$ in the $D_{\\rm IGW}$ formula from three-dimensional hydrodynamic simulations."},{"cited_title":"2001, A&A, 374, 646","cited_arxiv_id":null,"evidence_quote":"One of the previous AGB lithium yield calculations that the present 15 solar mass result exceeds by an order of magnitude."},{"cited_title":"2000, A&A, 363, 605","cited_arxiv_id":null,"evidence_quote":"The other previous AGB lithium yield calculation used for comparison in the yield grid."}],"review_version":1}