{"id":"02b11d0a-4f2a-4b65-9f0a-10d7d53d3e79","arxiv_id":"2506.05618","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Post-common-envelope helium white dwarfs with thin hydrogen envelopes cool in roughly 0.3 to 2 billion years, while thicker envelopes sustain hydrogen burning and keep the star bright for billions of years.","lead":"This paper computes new cooling histories for helium-core white dwarfs formed through common envelope binary evolution, showing that the amount of leftover hydrogen decides whether the star cools within hundreds of millions of years or keeps burning for billions. The work matters because these models are the tool astronomers use to convert observed temperatures and gravities into masses and ages for extremely low-mass white dwarfs in close binaries.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The young-age branch depends on the unmodeled choice of the CE bifurcation point: if entropy or binding energy sets MH, Table 2 ages shift from ~10 Myr to hundreds of Myr or Gyr, so the 5–130 Myr ages are conditional, not robust.","rationale":"After reading the paper in good faith, I find the strongest-claim conclusion internally consistent but explicitly conditional on the post-CE bifurcation point. The reader's weakest-assumption identification is correct and is the load-bearing hinge: Table 2 demonstrates that the same stars have ages spanning an order of magnitude between minimum and maximum non-flashing MH, and the paper's own conclusion disclaims CE energetics. I do not see a separate fatal flaw. The mcp criterion is motivated but not derived from a CE simulation, and the XH=0.1 'minimum' is even further below mcp, so the headline ages sit at the extreme end of the uncertainty. The entropy/binding-energy criteria would move remnants into the flashing/H-burning regime with much longer ages. This does not invalidate the paper as a scenario tool: the sequences, tables, and max-MH variants are useful and the public cooling tracks permit testing. Since the reader already conditioned acceptance on resolving the bifurcation-point question, no verdict adjustment is needed.","tokens_in":16271,"tokens_out":7401,"duration_ms":77823,"concrete_test":"Run a 3D hydrodynamic CE simulation of a ~1 M_sun RGB star with a ~0.27 M_sun He core and a low-mass companion, and directly measure the mass coordinate at which ejection stops and the surviving MH; compare with mcp and XH=0.1. If no simulation is available, recompute the Table 2 ages with LPCODE using initial MH set by the entropy-profile and binding-energy bifurcation points shown in Figs. 1–2; if the resulting ages at Teff = 12,000–27,000 K exceed ~300 Myr for a substantial fraction of the sample, the young-age branch is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that CE He WDs cool within 5–130 Myr at 12,000–27,000 K—is not a direct output of CE physics. It follows from choosing the maximum-compression point mcp (and, for the minimum-MH tables, the ten-times-smaller XH=0.1 coordinate) as the post-CE remnant boundary. The paper's own caveat is explicit: Section 4 states this 'is a structural analysis and does not account for the energetics required to reach such bifurcation points during CE ejection.' The sensitivity is large. In Table 2, the same observed objects yield ages of 10.5 vs 159 Myr (J0651+2844), 111 vs 430 Myr (J0822+3048), and 134 vs 708 Myr (J1738+2927) for minimum versus maximum non-flashing MH. The entropy-profile and binding-energy criteria shown in Figs. 1–2 place the bifurcation point at larger MH, in the flashing or H-burning regime, where the models themselves predict ages of several hundred Myr to >1 Gyr. Hence the young-age branch is a scenario prediction contingent on the low-MH end of a bifurcation-point uncertainty that the paper does not resolve.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents new LPCODE evolutionary sequences for helium-core white dwarfs (He WDs) formed through the common envelope (CE) channel, covering WD masses from 0.20 to 0.4352 Msun. The key variable is the residual hydrogen envelope mass MH, set by the assumed CE bifurcation point. Two branches are identified: non-flashing sequences, where MH is small and the WD cools rapidly, and flashing sequences, where hydrogen shell flashes delay the pre-WD phase. For minimal MH, the models predict cooling ages of 5-130 Myr for Teff between 12,000 and 27,000 K, much younger than standard SRLOF tracks. For larger MH, residual H burning can extend ages to several Gyr. The sequences are applied to a sample of observed He WDs with orbital periods below 0.1 day to infer masses, ages, and merger properties, and the cooling tracks are made publicly available.","tokens_in":16581,"tokens_out":3636,"duration_ms":37024,"significance":"If the young-age branch is correct, the paper has substantial implications for the interpretation of ELM WDs in compact binaries: masses and ages inferred with SRLOF models could be systematically biased, and the merger temperatures of CE-formed systems would be significantly lower than previously assumed. The paper provides a useful grid of post-CE tracks and explicit tables for observers, and it is a strength that the sequences are publicly available and that the dependence on MH is explored systematically. However, the central young-age result is not a direct output of CE physics; it is contingent on the adopted choice of the bifurcation point. The paper itself acknowledges this limitation, but the abstract and conclusions present the young ages as the main finding without quantifying the sensitivity to the unmodeled CE energetics.","major_comments":[{"comment":"The central claim that CE He WDs cool within 5-130 Myr at Teff = 12,000-27,000 K is a scenario prediction, not a robust result. It follows from adopting the maximum-compression point mcp and the even smaller XH=0.1 coordinate as the post-CE remnant boundary. The paper explicitly states in Sect. 4 that this 'is a structural analysis and does not account for the energetics required to reach such bifurcation points during CE ejection.' The entropy and binding-energy criteria shown in Figs. 1-2 place the bifurcation point at larger MH, in the flashing or H-burning regime. Table 2 quantifies the impact: the same observed objects yield ages differing by factors of 5-15 between the minimum and maximum non-flashing MH (e.g., J0651+2844: 10.5 vs 159 Myr; J0822+3048: 111 vs 430 Myr; J1738+2927: 134 vs 708 Myr). The manuscript does not provide a quantitative estimate of the energetic cost of stripping the envelope down to mcp, nor a defensible argument that the entropy/binding-energy criteria are excluded. Without such a sensitivity analysis, the young-age branch must be presented as conditional on the low-MH end of the bifurcation-point uncertainty, not as a definitive prediction.","section":"4. Conclusions / 3.1"},{"comment":"The upper bound on MH imposed by the condition that nuclear expansion does not exceed 1 Rsun (dashed red line in Fig. 3) is used to exclude models with larger MH. However, this 1 Rsun limit is motivated by orbital separations of P < 0.05 day systems, while Table 2 includes binaries with periods up to 0.0995 days. For a 0.0995-day orbit, the Roche lobe radius is substantially larger than 1 Rsun, so the adopted cutoff may be too restrictive for the long-period end of the sample. The paper should either derive the radius cutoff consistently with each system's orbital period or state explicitly that the 1 Rsun limit is a conservative simplification and discuss how relaxing it would affect the maximum-MH ages in Table 2.","section":"3.1 / Fig. 3"},{"comment":"The treatment of the adjustment time in the flashing sequences is unclear. Section 3.3 reports maximum pre-WD times from the end of CE of up to 230 Myr for the 0.2026 Msun flashing sequence, and the footnote to Fig. 4 states that in flashing sequences 'tau_adjust does not include the time spent in the first cooling track prior to the occurrence of the H flash.' The paper then states in Sect. 4 that the tables 'are also appropriate for He WD flashing sequences; in these cases, it is essential to consider pre-WD ages.' It is not clear whether the ages listed in Table 2 for the minimum-MH and maximum-MH columns include the full pre-WD phase for flashing cases. If they do not, then ages inferred for objects that actually followed the flashing branch would be underestimated by up to several hundred Myr. The manuscript should state explicitly which times are included in the tabulated ages and, if the flashing pre-WD phase is excluded, provide corrected ages or a prescription for adding tau_adjust.","section":"3.3 / Table 2"}],"minor_comments":[{"comment":"The abstract and introduction state that the sequences cover 0.20 to 0.42 Msun, but Table 1 includes a sequence at 0.4352 Msun. Please harmonize the stated mass range with the actual grid.","section":"Abstract / Table 1"},{"comment":"The binding energy definition as written yields a negative quantity for a bound envelope; the sign convention is not stated. Please define the sign explicitly or use the absolute value so that the 'steep increase' description is unambiguous.","section":"2.1, Eq. (1)"},{"comment":"Several SRLOF ages have lower uncertainties larger than the central value, e.g., J0935+4411 with 2.13 +/- 48 Myr, implying negative ages. This is presumably a consequence of the multi-solution nature of SRLOF tracks, but the presentation should be clarified, e.g., by quoting asymmetric ranges or stating that the uncertainty covers multiple solutions.","section":"Table 2"},{"comment":"The caption refers to the 'solid black line' as 'predictions from SRLOF mass transfer,' but the line presumably shows SRLOF cooling tracks rather than mass-transfer rates. Please reword to avoid ambiguity.","section":"Fig. 4"},{"comment":"The sentence 'these objects... reach the cooling track almost instantaneously after CE (within a century; see Table 1)' is contradicted by Table 1 for the 0.2026 and 0.2390 Msun sequences, which have tau_adjust values of 0.16-0.23 Myr and 0.07 Myr. The statement should be qualified to apply only to the lowest-MH sequences or to the specific masses shown.","section":"3.1"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid scenario-grid calculation with useful public tracks, but the headline result is presented more strongly than the underlying assumption warrants. The authors explicitly disclaim modeling of CE energetics, yet the abstract and conclusions do not carry the corresponding caveat into the central age predictions. I would advise the editor that the revision should either include a quantitative sensitivity analysis of the bifurcation-point choice (e.g., using entropy or binding-energy criteria to define MH) or substantially soften the language so that the young ages are clearly labeled as conditional on the low-MH end of the uncertainty. The paper is not circular; it is an exploratory evolution calculation, but the physical motivation for preferring mcp over the other criteria needs strengthening beyond a citation to Ivanova (2011)."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a well-executed scenario paper, not a measurement. It gives the first systematic grid of post-CE He WD cooling tracks and inference tables, and the tracks are public. But the headline young ages (5–130 Myr for the 12,000–27,000 K sample) depend directly on choosing the maximum-compression point as the post-CE remnant boundary, which the paper itself flags as a structural assumption. The tables are still worth having; the ages are conditional.\n\nWhat's genuinely new: LPCODE sequences for 0.20–0.42 Msun with three envelope-mass choices (minimum, mcp, maximum non-flashing), inference tables that map Teff and log g to mass and age for observed CE ELM WDs, and merger positions for those systems. The two-branch behavior (non-flashing vs. flashing) was already in Scherbak & Fuller 2023, but this grid is more systematic and the application to the Brown et al. 2020 sample is useful. The paper is transparent about the central limitation—Section 4 explicitly says it “does not account for the energetics required to reach such bifurcation points during CE ejection.” That honesty is a real strength.\n\nSoft spots: The young ages are not robust to the bifurcation-point uncertainty. As the stress-test note correctly says, entropy and binding-energy criteria put the bifurcation point at larger MH, in the flashing/H-burning regime, where the models themselves give ages of hundreds of Myr to >1 Gyr. The paper attempts to exclude those larger MH values by requiring that the post-CE remnant not expand beyond ~1 R⊙ (the orbital separation), but that expansion limit is also an assumption, not a CE calculation. For individual objects in Table 2, the inferred ages differ by factors of 5–10 between minimum and maximum non-flashing MH (e.g., J0651+2844: 10.5 vs 159 Myr; J1738+2927: 134 vs 708 Myr). So the central claim is conditional, and the paper should say “if the bifurcation point is mcp, then…” more prominently. That said, the paper does not overclaim—it clearly labels the analysis as structural—and the tables are explicitly offered as scenario tools. The mass shifts between CE and SRLOF sequences are smaller and more robust; those will survive even if the young ages don't. Citation pattern looks fair: Scherbak & Fuller and Vigna-Gómez et al. are cited where they're relevant, and the self-citations point to the group's own SRLOF tracks that serve as benchmarks.\n\nWho it's for: stellar modelers and observers working on ELM WDs in short-period binaries. The inference tables and public tracks will get cited. The paper deserves peer review; a referee should push for a tighter framing of the assumption-dependence, but the work is sound enough to publish after revision.\n\nRecommendation: send to review. Conditional acceptance at most, not desk rejection.","headline":"A solid, honest scenario grid for post-CE He WDs—but the young ages hinge on the unmodeled choice of the CE bifurcation point, so treat the tables as conditional tools.","tokens_in":17090,"tokens_out":3472,"would_cite":true,"duration_ms":31981,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["97.20.Rp","97.10.Cv","97.80.-d"],"model":"deepseek-v4-flash","headline":"Helium-core white dwarfs born in common-envelope ejection may cool much faster than standard Roche-lobe-overflow models predict: with a thin residual hydrogen skin they reach 12,000–27,000 K in 5–130 Myr, while a thicker skin delays…","keywords":["helium-core white dwarfs","extremely low-mass white dwarfs","common envelope evolution","close binaries","white dwarf cooling","hydrogen envelope mass","bifurcation point","LPCODE stellar evolution"],"falsifier":"Measure the hydrogen-layer mass of a short-period (P < 0.1 day) He WD with a known dynamical mass, for example through asteroseismic fits to an ELM pulsator or through the combination of cooling rate and gravitational-wave orbital decay for an eclipsing system such as J0651+2844. If $M_H$ is found to exceed a few times $10^{-4}\\,M_\\odot$ with significant residual burning, the predicted 5–130 Myr ages in the 12,000–27,000 K range would be ruled out, and the predicted low-temperature merges for P $\\gtrsim$ 0.07-day systems would not occur.","tokens_in":2256,"feed_emoji":"⭐","tokens_out":3028,"duration_ms":121482,"temperature":0.7,"pith_summary":"This paper argues that helium-core white dwarfs (He WDs) emerging from common-envelope (CE) binary evolution cool along paths very different from those born by stable Roche-lobe overflow (SRLOF), with the mass of the residual hydrogen envelope as the deciding factor. Using evolutionary sequences for 0.20–0.42 M_sun remnants computed with the LPCODE stellar evolution code, the authors find two regimes: below a threshold envelope mass, residual hydrogen burning is negligible and the white dwarf enters the cooling track almost immediately, while above it, burning resumes and stretches the cooling age from millions to billions of years. The location of the bifurcation point, the mass coordinate where the pre-CE red giant's envelope is cut, sets this envelope mass; at the maximum-compression point it leaves so little hydrogen that CE He WDs cool within 5–130 Myr across 12,000–27,000 K and reach about 300 Myr below 10,000 K. If correct, ages and masses inferred for short-period extremely low-mass white dwarfs from SRLOF-based models are systematically shifted, and the merger temperatures and times of these binaries change accordingly.","feed_headline":"Common-envelope white dwarfs cool in millions, not billions, of years","feed_subtitle":"Helium-core WDs stripped by common envelopes reach 12,000–27,000 K in 5–130 Myr, far younger than Roche-lobe models say.","key_machinery":"The load-bearing object is the bifurcation point, specifically the maximum-compression point $m_{\\rm cp}$, the innermost mass coordinate where $P/\\rho$ peaks in the H-burning shell before the CE ejection; it sets $M_H$ and therefore decides whether the remnant cools with or without nuclear support. The paper additionally maps three threshold lines in the $M_H$–$M_{\\rm WD}$ plane, the line where H burning dominates the luminosity, the line above which H-shell flashes occur, and the line above which nuclear expansion pushes the radius past roughly 1 $R_\\odot$, and uses them to delimit the allowed post-CE envelope masses.","core_discovery":"The paper claims that post-CE He WDs of a given mass follow two distinct cooling channels fixed by the initial residual hydrogen mass $M_H$. With the bifurcation point placed at the maximum-compression point $m_{\\rm cp}$, the local maximum of $P/\\rho$ inside the H-burning shell, $M_H$ is so small, a few times $10^{-4}\\,M_\\odot$, that residual burning supplies less than half the luminosity; these non-flashing sequences cool in 5–130 Myr across the observed $T_{\\rm eff}$ range 12,000–27,000 K and reach roughly 300 Myr below 10,000 K, much younger than SRLOF tracks. With larger $M_H$, still below the roughly 1 $R_\\odot$ envelope-expansion cap set by the tight orbits, residual H burning powers the envelope and ages grow to several Gyr; beyond a higher threshold, an early H-shell flash reshapes the envelope and extends the pre-WD phase to tens or hundreds of Myr. The same $M_H$ shifts the inferred mass at fixed $\\log g$ and $T_{\\rm eff}$, since a thinner envelope makes the WD more compact.","pith_inferences":["If the young-age branch is right, the observable lifetime of short-period He WDs at a given formation rate shrinks to roughly 10–100 Myr, so the systems we detect must be very recently formed, and the implied formation rate of these binaries is higher than SRLOF-based ages suggest.","A three-dimensional or energy-consistent CE simulation could place the bifurcation point deeper in the envelope than $m_{\\rm cp}$; in that case more hydrogen would survive, residual burning would resume, and ages would drift back toward SRLOF values, so the 5–130 Myr numbers are best read as a speed limit set by the most aggressive stripping scenario.","Asteroseismic or spectral inference of the hydrogen-layer mass of an ELM WD in a short-period binary would discriminate the scenarios directly: a measured $M_H$ above a few times $10^{-4}\\,M_\\odot$ would rule out the minimal-envelope tracks and their rapid cooling ages."],"forward_implications":["For the observed P < 0.1-day extremely low-mass white dwarfs, the minimal-envelope CE sequences put most ages below 200 Myr, many below 5 Myr, a large downward shift from SRLOF-based estimates.","At fixed $T_{\\rm eff}$ and $\\log g$, minimal-$M_H$ CE sequences give smaller masses than SRLOF sequences, while maximum-$M_H$ sequences give larger masses, introducing systematic offsets in inferred stellar parameters.","Systems with orbital periods $P \\gtrsim 0.07$ days are expected to merge below 8000 K after up to roughly 2 Gyr when the envelope is thin, whereas hydrogen-richer remnants merge at higher temperatures and younger ages.","Flashing sequences lengthen the pre-WD phase, up to about 260 Myr for a 0.20 $M_\\odot$ remnant, but remain shorter than SRLOF evolution; after a flash, the resulting WD is observationally indistinguishable from one born with a thin envelope.","The allowed range of $M_H$ is capped at $1.2\\times10^{-3}$ to $8\\times10^{-3}\\,M_\\odot$ because a thicker envelope would expand the star past the binary separation of about 1 $R_\\odot$ before cooling begins."],"supporting_citations":[{"why":"Supplies the maximum-compression-point criterion $m_{\\rm cp}$ that defines the minimal post-CE envelope mass and underlies the non-flashing branch.","marker":"Ivanova 2011"},{"why":"Prior work exploring how envelope mass alters post-CE He WD evolution; this paper extends it to physically motivated CE scenarios and a wider mass grid.","marker":"Scherbak & Fuller 2023"},{"why":"Provides the SRLOF He WD cooling sequences used as the comparison baseline for ages and masses.","marker":"Althaus et al. 2013"},{"why":"Supplies the observed sample of P < 0.1-day He + CO WD binaries that the paper treats as CE products and uses to infer masses, ages, and merger properties.","marker":"Brown et al. 2020"},{"why":"Basis for the CE energetics argument that sets the lower mass limit near 0.20 M_sun for CE-formed ELM WDs.","marker":"Li et al. 2019"},{"why":"Reviews the competing bifurcation-point criteria, including the $X_H = 0.1$ cut the paper adopts as a lower envelope-mass reference.","marker":"Tauris & Dewi 2001"},{"why":"Three-dimensional hydrodynamic simulations supporting near-complete envelope removal during CE ejection of low-mass giants.","marker":"Sand et al. 2020"},{"why":"Alternative SRLOF sequences, including rotational effects, used for comparison and as the motivation for future work on rotation in CE remnants.","marker":"Istrate et al. 2016"}],"fun_headline_variants":["Hydrogen envelope mass splits helium white dwarfs into fast and slow coolers","Common-envelope helium WDs: hydrogen mass sets age from Myr to Gyr","Flash or freeze: hydrogen envelope decides helium WD cooling","Two cooling paths for common-envelope helium white dwarfs","Helium-core WDs: residual hydrogen mass sets cooling timescale"],"cache_read_input_tokens":19200,"weakest_assumption_plain":"The rapid-cooling branch assumes common-envelope ejection strips the red giant down to the maximum-compression point $m_{\\rm cp}$, leaving only a very thin hydrogen envelope of order $10^{-4}\\,M_\\odot$, and treats the ejection as instantaneous structural removal without modeling the energetics required to reach such a bifurcation point; the paper explicitly states that it is a structural analysis that does not account for ejection energetics.","fun_headline_variants_meta":{"raw":{"variants":["Hydrogen envelope mass splits helium white dwarfs into fast and slow coolers","Common-envelope helium WDs: hydrogen mass sets age from Myr to Gyr","Flash or freeze: hydrogen envelope decides helium WD cooling","Two cooling paths for common-envelope helium white dwarfs","Helium-core WDs: residual hydrogen mass sets cooling timescale"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001041,"raw_usage":{"total_tokens":4475,"prompt_tokens":1137,"completion_tokens":3338,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":753,"completion_tokens_details":{"reasoning_tokens":3260}},"tokens_in":753,"tokens_out":3338,"duration_ms":28770,"temperature":1.0,"reasoning_tokens":3260,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T10:13:32.996361+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the hydrogen-layer mass of a short-period (P < 0.1 day) He WD with a known dynamical mass, for example through asteroseismic fits to an ELM pulsator or through the combination of cooling rate and gravitational-wave orbital decay for an eclipsing system such as J0651+2844. If $M_H$ is found to exceed a few times $10^{-4}\\,M_\\odot$ with significant residual burning, the predicted 5–130 Myr ages in the 12,000–27,000 K range would be ruled out, and the predicted low-temperature merges for P $\\gtrsim$ 0.07-day systems would not occur.","supporting_citations":[{"cited_title":"& Fuller, J","cited_arxiv_id":null,"evidence_quote":"Prior work exploring how envelope mass alters post-CE He WD evolution; this paper extends it to physically motivated CE scenarios and a wider mass grid."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reviews the competing bifurcation-point criteria, including the $X_H = 0.1$ cut the paper adopts as a lower envelope-mass reference."},{"cited_title":"T., Schneider, F","cited_arxiv_id":null,"evidence_quote":"Three-dimensional hydrodynamic simulations supporting near-complete envelope removal during CE ejection of low-mass giants."}],"review_version":1}