{"id":"dc212eb8-f100-4d5c-8c36-5bb6028e31d5","arxiv_id":"2505.13313","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"The early post-AGB instability produces HRD loops and photospheric 7Be enhancement but does not significantly alter the final white dwarf hydrogen envelope mass.","lead":"This paper simulates the late evolution of 1 and 2 solar mass stars and finds that a recently discovered instability does not change how much hydrogen ends up on the white dwarf, but it does make the star loop repeatedly across the brightness and temperature diagram and push beryllium to its surface. If real, those loops could be used to weigh stars right after they leave the asymptotic giant branch.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Most load-bearing concern: the final MH values in Tables 1–2 and Fig. 7 are obtained by forcing dynamically unstable outer layers (T<10^6 K) into thermal equilibrium, so the title question is not settled until a genuine hydrodynamic calculation is made.","rationale":"The reader's weakest assumption is exactly the hydrostatic continuation through dynamically unstable phases, and I agree that it is the most load-bearing point. The paper is careful and self-aware, explicitly warning that all hydrogen could be removed in a full hydrodynamic treatment, and the conditional verdict is appropriate. The concern is not merely a caveat: the central claim about MH is produced by a prescription whose validity is the very question at issue, and the same prescription is used in both the EPAGBI and EPAGBI-suppressed branches, so the comparison cannot isolate the physical effect of the instability. The suggested hydrodynamic rerun is expensive but would settle whether the tabulated MH values are even upper limits. I would keep the paper as CONDITIONAL rather than rejecting it, because the characterization of the EPAGBI phase, the HRD loops, and the exploratory 7Be predictions are honestly framed model results that merit further study; however, the title question remains open until the hydrodynamics is actually computed.","tokens_in":13009,"tokens_out":7746,"duration_ms":79510,"concrete_test":"As a decisive check, take the 1 and 2 M_sun post-AGB models at the onset of the density-inversion/super-Eddington envelope state (the red circle in Fig. 2 and the corresponding stage in §2.3) and evolve them with a 1D or 2D radiation-hydrodynamics code that includes the outer envelope, using the same composition, mass-loss prescription, and updated opacity tables (including high-Z cases), for at least several dynamical timescales. If the hydrogen-rich layer is ejected, MH at the start of the WD cooling track is not the tabulated value and the paper's MH conclusion fails; if the layer survives with MH within a factor of 2 of the tabulated values, the forced-equilibrium continuation is validated and the concern is retired.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The quantitative answer to the title question rests on the continuation procedure described in §2.1, §2.2, and §2.4: when the envelope becomes dynamically unstable (luminosity exceeding the local Eddington limit and a density inversion), the hydrostatic DEUCES code is advanced by forcing the outer layers with T<10^6 K to be in thermal equilibrium. This is an ad hoc numerical constraint, not a model of envelope ejection, and the author explicitly cautions in the Abstract and §3 that 'because hydrodynamic behavior is not included... it is possible that all hydrogen would be removed.' The comparison in Fig. 7 and §2.2 ('EPAGBI does not lead to significant additional mass loss') is therefore between two calculations that both use the same forced equilibrium at the VLTP/hydrodynamic phases, so it cannot establish whether EPAGBI actually affects MH. Tables 1 and 2 final values (e.g., 5.52e-6 and 3.51e-7 M_sun) should be read as upper estimates from an extrapolated continuation, not as physically validated predictions. The HRD loop timescale–mass relation is also based on only two models, one with and one without convective overshoot, but the primary unsupported step is the hydrodynamic fate of the hydrogen envelope.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper investigates whether the Early Post-AGB Instability (EPAGBI) affects the final hydrogen envelope mass of white dwarfs. Using the DEUCES stellar evolution code, the author evolves 1 and 2 solar-mass, solar-metallicity models from the pre-main sequence through the AGB and post-AGB phases to the white dwarf cooling track, computing cases with the EPAGBI followed in detail and with the EPAGBI suppressed by enforcing large time steps. The reported hydrogen masses at the start of the white dwarf cooling phase are in the range inferred from asteroseismology, but the author cautions that hydrodynamic behavior is not modeled and that all hydrogen might be removed. The paper argues that the main observable impact of EPAGBIs is the production of loops in the Hertzsprung-Russell diagram, with loop timescales that decrease with mass, and the convective transport of 7Be to the photosphere in amounts up to about 400 times the solar photospheric value.","tokens_in":13250,"tokens_out":4586,"duration_ms":42513,"significance":"The paper is an honest and clearly written exploration of a recently identified instability, with explicit statements of its limitations and a control calculation in which the EPAGBI is suppressed. The proposed HRD loop timescale as a potential stellar mass indicator is an interesting, falsifiable prediction, and the predicted 7Be enhancement is a concrete observable signature. The use of a stellar evolution code without free parameters tuned to the target quantities, and the explicit comparison between EPAGBI and EPAGBI-suppressed runs, are strengths. However, the central quantitative conclusion about hydrogen envelope mass is undermined by the need to force the outer layers into thermal equilibrium when the envelope becomes dynamically unstable, so the paper's main claim about MH is not yet fully supported.","major_comments":[{"comment":"The central quantitative claim that the EPAGBI does not significantly affect the final hydrogen envelope mass is not established by the present calculations. In both the EPAGBI and EPAGBI-suppressed runs, the evolution is continued past the hydrodynamic instability by forcing the outer layers at T < 10^6 K to be in thermal equilibrium, which is a numerical continuation procedure rather than a physical treatment of envelope ejection. Because this same forced equilibrium is applied in both branches, the comparison in Fig. 7 measures differences in the timing and outcome of the imposed continuation, not the dynamical effect of the EPAGBI on the hydrogen envelope. The abstract and §3 correctly state that all hydrogen might be removed hydrodynamically; this caveat should be carried through to the wording of the conclusion in §2.2 that 'EPAGBI does not lead to significant additional mass loss,' which currently overstates what a hydrostatic calculation can establish.","section":"§2.1, §2.2, §2.4, Tables 1–2, Fig. 7"},{"comment":"The loop timescale–mass relation is a central positive result, but it rests on only two computed models with AGB-departure masses of 0.567 and 0.642 M_sun and one extrapolated estimate near 0.72 M_sun. No uncertainty estimate is given, nor is the sensitivity to input assumptions such as convective overshoot, mass-loss prescription, or metallicity explored. The statement that the timescale 'could provide a way to determine the stellar mass just after AGB departure' is therefore plausible but not demonstrated; the paper should either add intermediate-mass models or present the relation as a tentative suggestion with explicit caveats about the sparse calibration.","section":"§3, Figs. 2 and 9"},{"comment":"For the 2 M_sun model, the photospheric heavy-element abundance reaches Z = 0.48, and the author notes that the low-temperature opacity tables do not cover such high Z values. This means that the later evolutionary phases, including the final MH value in Table 2 and the second blue loop in Fig. 12, are computed with opacity data outside their stated domain. The 2 M_sun results should be explicitly labeled as exploratory, and the mass-dependence argument in §3 should not rely heavily on them without additional justification or a test of the opacity extrapolation.","section":"§2.3, Table 2, Fig. 12"}],"minor_comments":[{"comment":"The sentence 'The red vertical marks correspond to evolutionary stages given in table 2' should refer to Table 1, since Fig. 7 displays the 1 M_sun models and the stages listed in the text are those of Table 1.","section":"Fig. 7 and text"},{"comment":"There are typographical inconsistencies in the abbreviation: 'EAGBI' appears instead of 'EPAGBI' in §2.2, and the axis label in Fig. 3 reads 'Time from EAGPI' rather than 'Time from EPAGBI onset'.","section":"§2.2 and Fig. 3"},{"comment":"The phrase 'the the EPAGBI' contains a duplicated article, and 'radiative driven winds' should be 'radiatively driven winds'.","section":"§2.2"},{"comment":"Since the forced thermal equilibrium in the outer layers with T < 10^6 K is central to the limitation of the models, the paper would benefit from a brief technical description of how this constraint is imposed (for example, which variables are adjusted and over what timescale), so that readers can assess the effect of the procedure.","section":"§2.1 and §2.4"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is honest about its central limitation, and the referee report focuses on the mismatch between the stated conclusion and the forced-equilibrium continuation. The HRD-loop and 7Be parts are the more novel and better-supported contributions, while the hydrogen-envelope-mass answer is conditional on a future hydrodynamic treatment. The two-model calibration of the loop timescale relation is too thin for the strength of the claim in §3."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe paper is worth reading for its new quantitative predictions on the early post-AGB instability (EPAGBI): HRD loops on timescales of ~100 yr for a 0.567 Msun remnant, ~10 yr for 0.642 Msun, and an extrapolated ~1 yr for 0.72 Msun, plus surface 7Be enhancement up to ~400 times solar. These go beyond Gautschy 2023 and offer testable signatures of the instability. The author is careful and honest about the modeling: the code is hydrostatic, the dynamically unstable phases are handled by forcing outer layers into thermal equilibrium, and the 2 Msun model exceeds the low-T opacity tables. No parameter is fitted to MH, loop times, or 7Be abundances; they are outputs. The suppressed-EPAGBI run is a sensible control.\n\nThe soft spot is exactly the title question. The final MH values in Tables 1-2 and Fig. 7 come from a continuation procedure that is not a model for ejection. Both the EPAGBI and suppressed runs use the same forced equilibrium when the envelope becomes dynamically unstable, so the comparison cannot establish whether the instability changes MH. The author says this himself: all hydrogen could be removed in a hydrodynamic treatment. These numbers are upper estimates, not predictions. The loop-mass scaling also rests on only two computed masses, so it is provisional.\n\nWho should read this? Stellar evolution and white dwarf folks, and observers looking for variability in post-AGB stars. It is a solid, honest modeling study that deserves a serious referee. I would send it to peer review with a request that the author frame the MH result as an upper limit and discuss what a hydrodynamic calculation might change. Not a desk reject.\n\nRegards.","headline":"New predictions for EPAGBI loops and 7Be are worth a look, but the MH answer isn't settled because the hydrostatic continuation limits what the models can claim.","tokens_in":13830,"tokens_out":3451,"would_cite":false,"duration_ms":32201,"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":"A thermally driven instability after the AGB produces HRD loops whose period falls with stellar mass, offering a new mass diagnostic and sending 7Be to the surface.","keywords":["early post-AGB instability","white dwarf hydrogen envelopes","Hertzsprung-Russell diagram loops","7Be convective dredge-up","DAV asteroseismology","asymptotic giant branch evolution","planetary nebula central stars","stellar mass determination"],"falsifier":"A long-term photometric monitoring campaign of post-AGB stars and planetary nebula central stars that finds no HRD loops on the predicted 1 to 100 yr timescales would falsify the loop signature; likewise, a detection or non-detection of the Be II 313 nm doublet at the predicted EPAGBI phase would test the $^7$Be transport. A hydrodynamic simulation showing complete ejection of the hydrogen layer would falsify the paper's conclusion that the EPAGBI can leave thin but nonzero hydrogen envelopes.","tokens_in":2183,"feed_emoji":"🌟","tokens_out":4925,"duration_ms":98915,"temperature":0.7,"pith_summary":"The paper asks whether the Early Post-AGB Instability (EPAGBI), a recently identified thermal and dynamical instability in stars just after they leave the asymptotic giant branch, can explain why white dwarfs have much thinner hydrogen envelopes than standard evolution predicts. Evolving 1 and 2 solar-mass models through the AGB and onto the white-dwarf cooling track, the author finds that the final hydrogen envelope masses fall in the range inferred from asteroseismology, not the canonical $10^{-4}\\,M_\\odot$. However, because the hydrostatic code cannot follow the dynamical phase and the outer layers are forced into thermal equilibrium, these masses are upper limits and it remains possible that all hydrogen is ejected. The paper's central new result is that the instability produces characteristic loops in the Hertzsprung-Russell diagram whose period shrinks with stellar mass, about 100 years at $0.567\\,M_\\odot$, 10 years at $0.642\\,M_\\odot$, and an estimated 1 year near $0.72\\,M_\\odot$, making the loop timescale a potential stellar-mass diagnostic. It also finds that $^7$Be, the lithium precursor, is convected to the photosphere at up to about 400 times the solar photospheric mass fraction, offering a spectroscopic signature via the Be II doublet.","feed_headline":"Post-AGB instability leaves loops that reveal a star's mass","feed_subtitle":"Modeled 1 and 2 solar-mass stars show 100- and 10-year HRD loops, plus surface 7Be that may be detectable.","key_machinery":"The load-bearing object is the Early Post-AGB Instability itself: a thermally driven radial pulsation of the hydrogen-burning shell that develops after a model leaves the AGB, with driving in the partial-ionization zones of hydrogen and helium. In the hydrostatic stellar evolution code used here, the instability appears as an exponentially growing oscillation of hydrogen-burning luminosity; when the envelope exceeds the local Eddington luminosity and develops a density inversion, the calculation is continued only by forcing the outer layers ($T<10^6$ K) into thermal equilibrium. This machinery produces the HRD loops, the episodic dredge-up that mixes helium-burning products into the photosphere, and the convective transport of $^7$Be to the surface. The comparison calculation suppresses the instability by forcing time steps larger than the growth time, showing that the loops disappear and that the final hydrogen mass differs only modestly.","core_discovery":"The central claim is that the EPAGBI does not by itself settle the hydrogen-envelope-mass conflict, because the computed $M_{\\mathrm{H}}$ values, though consistent with DAV asteroseismology, depend on a hydrostatic treatment that is known to break down; the author states directly that all hydrogen could be removed dynamically. What the instability does produce is a distinctive observable: repeated loops in the HRD caused by thermal pulsations of the hydrogen-burning shell, with a loop period that decreases steeply with the mass at AGB departure. In the 1 $M_\\odot$ model (departure mass $0.567\\,M_\\odot$) the loop period is about 100 yr with a period-doubling route to slightly chaotic behavior before a final thermal pulse; in the 2 $M_\\odot$ model (departure mass $0.642\\,M_\\odot$) the period is about 10 yr and no thermal pulse develops during the loop phase. The loop timescale-mass relation, extrapolated to about 1 yr at $0.72\\,M_\\odot$, is proposed as a way to measure the mass of a star just after AGB departure if the loops are detected. The other predicted signature is surface $^7$Be enhancement through the Cameron-Fowler process, reaching approximately 400 times the solar photospheric $^7$Be fraction, which should be sought in the Be II 313 nm resonance doublet.","pith_inferences":["If the HRD loops are real, photometric monitoring campaigns of bright post-AGB stars over decades could catch the roughly 10 yr loop for $0.64\\,M_\\odot$ stars, turning a theoretical instability into a mass measurement independent of atmospheric modeling.","A hydrodynamic simulation of the same models would settle whether the hydrogen envelope survives; if it is fully ejected, the EPAGBI would point toward the origin of hydrogen-deficient central stars rather than of thin-hydrogen white dwarfs.","The $^7$Be surface enhancement is a clean diagnostic of deep convective mixing during the EPAGBI, so surveying post-AGB stars for the Be II doublet would test both the instability and the mixing physics.","The loop timescale-mass relation, if calibrated with more masses, might extend the empirical initial-final mass relation to the post-AGB boundary."],"forward_implications":["HRD loops with periods of roughly 100 yr at $0.567\\,M_\\odot$ and 10 yr at $0.642\\,M_\\odot$ should be absent when the EPAGBI is suppressed, so their presence is a direct test of the instability.","If the loop timescale-mass relation holds, measuring the period of HRD looping in a post-AGB star or planetary nebula central star gives the stellar mass just after AGB departure, especially for masses above about $0.72\\,M_\\odot$ where the period is about 1 yr.","Final hydrogen envelope masses in both the 1 and 2 $M_\\odot$ models are much lower than the canonical $10^{-4}\\,M_\\odot$ and fall in the range inferred from DAV asteroseismology, but should be read as upper limits because hydrodynamics is not included.","$^7$Be is brought to the photosphere at up to about 400 times the solar photospheric mass fraction during the EPAGBI phase and may be detectable through the Be II 313.0 and 313.1 nm resonance doublet, whereas Li I detection is not expected.","EPAGBI-induced cyclic mass loss may leave an imprint on planetary nebula morphology."],"supporting_citations":[{"why":"Identifies the Early Post-AGB Instability and its driving in hydrogen and helium partial-ionization zones, which this paper builds on.","marker":"Gautschy 2023"},{"why":"Provides the DEUCES evolution code used to compute the models.","marker":"Lawlor & MacDonald 2023"},{"why":"Supplies the canonical stellar evolution prediction of a hydrogen envelope mass near $10^{-4}\\,M_\\odot$ that the paper's models are compared against.","marker":"Schoenberner 1981"},{"why":"Gives the asteroseismic range of $\\log(M_{\\mathrm{H}}/M_*)$ for DAV stars that the computed envelope masses are checked against.","marker":"Castanheira & Kepler 2009"},{"why":"Asteroseismic analysis of ZZ Ceti stars showing that most DAVs have thinner hydrogen envelopes than standard evolution predicts.","marker":"Romero et al. 2012"},{"why":"Kepler and K2 asteroseismic analysis providing the modern distribution of hydrogen envelope masses in DAVs.","marker":"Hall et al. 2023"},{"why":"Explains the density inversion mechanism invoked when the envelope exceeds the local Eddington luminosity.","marker":"Joss et al. 1973"},{"why":"Describes the process by which $^7$Be is produced and later decays to $^7$Li, the basis of the predicted beryllium signature.","marker":"Cameron & Fowler 1971"}],"fun_headline_variants":["AGB instability loops reveal stellar mass after departure","7Be spikes in post-AGB stars may be detectable","Loops in HRD track mass of post-AGB stars","Post-AGB instability leaves loops that reveal star mass"],"cache_read_input_tokens":15872,"weakest_assumption_plain":"The calculation relies on a hydrostatic code and continues through dynamic instability only by forcing outer layers with $T<10^6$ K into thermal equilibrium, so the computed hydrogen masses are upper limits and all hydrogen could in reality be ejected.","fun_headline_variants_meta":{"raw":{"variants":["AGB instability loops reveal stellar mass after departure","7Be spikes in post-AGB stars may be detectable","Loops in HRD track mass of post-AGB stars","Post-AGB instability leaves loops that reveal star mass"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000556,"raw_usage":{"total_tokens":2778,"prompt_tokens":1207,"completion_tokens":1571,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":823,"completion_tokens_details":{"reasoning_tokens":1503}},"tokens_in":823,"tokens_out":1571,"duration_ms":11182,"temperature":1.0,"reasoning_tokens":1503,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T20:15:15.874730+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A long-term photometric monitoring campaign of post-AGB stars and planetary nebula central stars that finds no HRD loops on the predicted 1 to 100 yr timescales would falsify the loop signature; likewise, a detection or non-detection of the Be II 313 nm doublet at the predicted EPAGBI phase would test the $^7$Be transport. A hydrodynamic simulation showing complete ejection of the hydrogen layer would falsify the paper's conclusion that the EPAGBI can leave thin but nonzero hydrogen envelopes.","supporting_citations":[{"cited_title":"On an Early - Post-AGB Instability","cited_arxiv_id":"2303.11374","evidence_quote":"Identifies the Early Post-AGB Instability and its driving in hydrogen and helium partial-ionization zones, which this paper builds on."}],"review_version":1}