REVIEW 3 major objections 4 minor 32 references
Early Post Asymptotic Giant Branch Instability: Does it Affect White Dwarf Hydrogen Envelope Mass?
T0 review · 3 major / 4 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read 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.
desk verdict 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. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
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.
What would settle it
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.
Extended reading notes
Core claim
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.
Load-bearing premise
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.
Editorial extensions
If this is right
- 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.
Reading between the lines
- 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.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (3)
- [§2.1, §2.2, §2.4, Tables 1–2, Fig. 7] 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.
- [§3, Figs. 2 and 9] 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.
- [§2.3, Table 2, Fig. 12] 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.
minor comments (4)
- [Fig. 7 and text] 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.
- [§2.2 and Fig. 3] 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'.
- [§2.2] The phrase 'the the EPAGBI' contains a duplicated article, and 'radiative driven winds' should be 'radiatively driven winds'.
- [§2.1 and §2.4] 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.
Circularity Check
No significant circularity: MH values, HRD loop timescales, and 7Be abundances are model outputs, not fitted inputs; the acknowledged hydrostatic limitation is a validity caveat, not a circular step.
full rationale
The paper's derivation chain is self-contained and does not reduce to its inputs. Stellar models are evolved from the pre-main sequence with the DEUCES code using standard, externally sourced physics (e.g., Schröder et al. convective overshoot, Vink et al. mass-loss rates). The central outputs—MH at the end of the evolution, HRD loop timescales, and photospheric 7Be abundances—are computed from the models, not imposed as inputs. The comparison run with the EPAGBI suppressed is a control calculated with the same code and physics, so the conclusion that EPAGBI does not cause significant additional mass loss is a difference between two model outputs, not a tautology. The final MH values are compared with asteroseismic determinations as an external benchmark, but no parameter is fitted to those determinations. The self-citation to Lawlor & MacDonald (2023) is for the DEUCES code itself, which is the computational tool rather than a load-bearing cited theorem or uniqueness claim; using one's own code is not circularity. The continuation procedure in which outer layers with T < 10^6 K are forced into thermal equilibrium is an acknowledged limitation, and the paper explicitly cautions that hydrodynamic behavior is not included and that all hydrogen could be removed. This weakens the quantitative reliability of the MH results but is a stated validity caveat, not a circular reduction: the final MH is not assumed by construction, and the same caveat applies to both the EPAGBI and suppressed-EPAGBI runs. No equation in the paper is identical by construction to an input quantity, and no fitted parameter is renamed as a prediction. Therefore, the paper exhibits no significant circularity.
Assumptions & free parameters
free parameters (5)
- Initial stellar mass =
1 M_sun and 2 M_sun
- Initial metallicity =
Z = 0.01661
- Convective overshoot parameter =
on (Schroder et al. 1997) for 2 M_sun; off for 1 M_sun baseline
- Time step floor for EPAGBI suppression =
10 yr (minimum) until Teff = 9700 K
- Mass loss prescriptions =
Reimers (eta_r = 0.477), Vink et al. (1999), radiatively driven wind
assumptions (4)
- domain assumption Hydrostatic equilibrium holds for the stellar structure during almost all of the evolution
- domain assumption Instantaneous convective adjustment with mixing length theory is valid
- domain assumption Low temperature opacity tables are valid for the modeled compositions
- domain assumption Mass loss is described by literature prescriptions
Cite this review
Pith. "Pith review of Early Post Asymptotic Giant Branch Instability: Does it Affect White Dwarf Hydrogen Envelope Mass?." pith.science (2026). https://pith.science/paper/LI4B5RTF
@misc{pith2026250513313,
author = {Pith},
title = {Pith review of: Early Post Asymptotic Giant Branch Instability: Does it Affect White Dwarf Hydrogen Envelope Mass?},
year = {2026},
howpublished = {\url{https://pith.science/paper/LI4B5RTF}},
note = {Machine review of arXiv:2505.13313}
}
read the original abstract
Although most white dwarf stars have hydrogen-dominated atmospheres, a significant fraction have atmospheres in which hydrogen is spectroscopically absent, with the fraction of hydrogen-free atmospheres varying with effective temperature. Estimates of the total mass of hydrogen, MH, in the stellar envelope from either asteroseismology or spectral evolution are at odds with predicted values from theoretical stellar evolution modeling. Recent work has found that models in the early post Asymptotic Giant Branch (AGB) phase of evolution can exhibit thermally and dynamical unstable behavior. Here we investigate whether this Early Post AGB Instability (EPAGBI) can help resolve the conflict in MH values determined from white dwarf spectral evolution, analysis of DAV pulsations and canonical stellar evolution modeling, by evolving models of mass 1 and 2Msun through the AGB phases and to the white dwarf cooling track. The MH values at the end of the calculations are in the range consistent with asteroseismic determinations. The major impact of EPAGBIs is that they cause loops in the HRD, which are absent when the EPAGBI is suppressed. Such loops might be detectable in a long-term monitoring program, or by their imprint on planetary nebula morphology imparted by the cyclically varying mass loss rate. Since the characteristic timescale of the looping in the HRD depends on the stellar mass, it could provide a way to determine the stellar mass just after AGB departure. Another EPAGBI signature is the production of Li by the Cameron-Fowler process. During the EPAGBI phase the photospheric temperature is always too high for the Li I resonance line to be detected. However, 7Be is convected to the photosphere in significant amounts (up to 400 times the solar photospheric mass fraction) at various times in the EPAGBI phase, which may be detectable by observing the Be II resonance doublet.
Figures
Figures from the paper (6 more)
Reference graph
Works this paper leans on
-
[1]
1998, A&A, 330, 641, doi: 10.48550/arXiv.astro-ph/9710158 B´ edard, A
Asplund, M. 1998, A&A, 330, 641, doi: 10.48550/arXiv.astro-ph/9710158 B´ edard, A. 2024, Ap&SS, 369, 43, doi: 10.1007/s10509-024-04307-5 B´ edard, A., Bergeron, P., & Brassard, P. 2023, ApJ, 946, 24, doi: 10.3847/1538-4357/acbb62 B´ edard, A., Bergeron, P., Brassard, P., & Fontaine, G. 2020, ApJ, 901, 93, doi: 10.3847/1538-4357/abafbe
-
[2]
Bergeron, P., Fontaine, G., Bill` eres, M., Boudreault, S., & Green, E. M. 2004, ApJ, 600, 404, doi: 10.1086/379808 Bl¨ ocker, T. 2001, Ap&SS, 275, 1, doi:10.1023/A:1002777931450
doi:10.1086/379808 2004
-
[3]
Boesgaard, A. M., Deliyannis, C. P., Lum, M. G., & Chontos, A. 2022, ApJ, 941, 21, doi: 10.3847/1538-4357/ac9625
-
[4]
2018, MNRAS, 478, 4513, doi: 10.1093/mnras/sty1281
Buder, S., Asplund, M., Duong, L., et al. 2018, MNRAS, 478, 4513, doi: 10.1093/mnras/sty1281
-
[5]
Cameron, A. G. W., & Fowler, W. A. 1971, ApJ, 164, 111, doi: 10.1086/150821
doi:10.1086/150821 1971
-
[6]
Campbell, S. W., & Lattanzio, J. C. 2008, A&A, 490, 769, doi: 10.1051/0004-6361:200809597
-
[7]
Castanheira, B. G., & Kepler, S. O. 2009, MNRAS, 396, 1709, doi: 10.1111/j.1365-2966.2009.14855.x
arXiv 2009
-
[8]
1993, ApJ, 407, 649, doi: 10.1086/172547
Castellani, M., & Castellani, V. 1993, ApJ, 407, 649, doi: 10.1086/172547
doi:10.1086/172547 1993
Show all 32 references
-
[9]
D., Kalirai, J
Cummings, J. D., Kalirai, J. S., Tremblay, P. E., Ramirez-Ruiz, E., & Choi, J. 2018, ApJ, 866, 21, doi: 10.3847/1538-4357/aadfd6 D’Cruz, N. L., Dorman, B., Rood, R. T., & O’Connell, R. W. 1996, ApJ, 466, 359, doi: 10.1086/177515
2018 doi
-
[10]
1987, in IAU Colloq
Fontaine, G., & Wesemael, F. 1987, in IAU Colloq. 95: Second Conference on Faint Blue Stars, ed. A. G. D. Philip, D. S. Hayes, & J. W. Liebert, 319–326
1987
-
[11]
Fujimoto, M. Y. 1977, PASJ, 29, 331
1977
- [12]
-
[13]
Gianninas, A., Bergeron, P., Dupuis, J., & Ruiz, M. T. 2010, ApJ, 720, 581, doi: 10.1088/0004-637X/720/1/581
2010 doi
-
[14]
G., & Bischoff-Kim, A
Hall, W., Castanheira, B. G., & Bischoff-Kim, A. 2023, ApJ, 948, 74, doi: 10.3847/1538-4357/acc52c
2023 doi
-
[15]
1965, ApJ, 141, 993, doi: 10.1086/148193 —
Iben, I. 1965, ApJ, 141, 993, doi: 10.1086/148193 —. 1984, ApJ, 277, 333, doi: 10.1086/161700
1965 doi
-
[16]
B., Truran, J
Iben, I., Kaler, J. B., Truran, J. W., & Renzini, A. 1983, ApJ, 264, 605, doi: 10.1086/160631
1983 doi
-
[17]
1985, ApJ, 296, 540, doi: 10.1086/163473 —
Iben, I., & MacDonald, J. 1985, ApJ, 296, 540, doi: 10.1086/163473 —. 1995, in White Dwarfs, ed. D. Koester & K. Werner, Vol. 443, 48
1985 doi
-
[18]
Iben, I., & Tutukov, A. V. 1984, ApJ, 282, 615, doi: 10.1086/162241
1984 doi
-
[19]
C., Salpeter, E
Joss, P. C., Salpeter, E. E., & Ostriker, J. P. 1973, ApJ, 181, 429, doi: 10.1086/152060
1973 doi
-
[20]
M., & MacDonald, J
Lawlor, T. M., & MacDonald, J. 2003, ApJ, 583, 913, doi: 10.1086/345411 —. 2006, MNRAS, 371, 263, doi: 10.1111/j.1365-2966.2006.10641.x —. 2023, MNRAS, 525, 4700, doi: 10.1093/mnras/stad2582
2003
-
[21]
J., et al
Liebert, J., Wesemael, F., Hansen, C. J., et al. 1986, ApJ, 309, 241, doi: 10.1086/164595 L´ opez-Sanjuan, C., Tremblay, P. E., Ederoclite, A., et al. 2022, A&A, 658, A79, doi: 10.1051/0004-6361/202141746
1986 doi
-
[22]
1991, ApJ, 371, 719, doi: 10.1086/169937
MacDonald, J., & Vennes, S. 1991, ApJ, 371, 719, doi: 10.1086/169937
1991 doi
- [23]
- [24]
-
[25]
2018, ApJ, 857, 56, doi: 10.3847/1538-4357/aab713
Rolland, B., Bergeron, P., & Fontaine, G. 2018, ApJ, 857, 56, doi: 10.3847/1538-4357/aab713
2018 doi
-
[26]
D., C´ orsico, A
Romero, A. D., C´ orsico, A. H., Althaus, L. G., et al. 2012, MNRAS, 420, 1462, doi: 10.1111/j.1365-2966.2011.20134.x
2012
-
[27]
1979, A&A, 79, 108 —
Schoenberner, D. 1979, A&A, 79, 108 —. 1981, A&A, 103, 119 Schr¨ oder, K.-P., Pols, O. R., & Eggleton, P. P. 1997, MNRAS, 285, 696, doi: 10.1093/mnras/285.4.696
1979 doi
-
[28]
2000, A&A, 359, 1042
Unglaub, K., & Bues, I. 2000, A&A, 359, 1042
2000
- [29]
-
[30]
1994, A&A, 290, 807
Wagenhuber, J., & Weiss, A. 1994, A&A, 290, 807
1994
-
[31]
A., Mari, M
Weidmann, W. A., Mari, M. B., Schmidt, E. O., et al. 2020, A&A, 640, A10, doi: 10.1051/0004-6361/202037998
2020 doi
-
[32]
F., & Liebert, J
Wesemael, F., Green, R. F., & Liebert, J. 1985, ApJS, 58, 379, doi: 10.1086/191046 This paper was built using the Open Journal of As- trophysics LATEX template. The OJA is a journal which provides fast and easy peer review for new papers in the astro-ph section of the arXiv, m...
1985 doi
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.