{"id":"b242c661-3825-4a9c-8e0e-13db3d198847","arxiv_id":"2412.11712","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A new analytic prescription for hydrogen-rich shells in rapid binary population synthesis improves the match between simulated and observed subdwarf B stars, and a 162-configuration parameter study shows all tested binary physics choices strongly affect the population.","lead":"Astronomers built 162 synthetic binary star populations to study how subdwarf B stars are born, and added a new mathematical recipe for the thin hydrogen layer such stars can keep. The recipe visibly improves how simulated stars match the measured temperatures and surface gravities of real subdwarf B stars.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (15) as printed cannot reproduce the lifetimes in Fig. B.1: the above-MHeF branch has a negative base under a fractional exponent, and the below-MHeF branch gives about 1 Myr instead of roughly 150 Myr. The H-shell prescription's central validation is therefore not reproducible.","rationale":"I read the paper in good faith: the parameter study is extensive, the limitations are acknowledged, and the idea of adding hydrogen-rich envelopes to rapid BPS output is physically motivated and worth pursuing. However, the most load-bearing condition for the central claim is that the analytic prescription actually reproduces the detailed Bauer & Kupfer (2021) models on which it is built. The manuscript asserts this in Section 2.2.1 and visualizes it in Fig. B.1, but direct substitution into the printed Equation (15) gives either an undefined negative base for the above-MHeF branch or a lifetime of about 1 Myr for the below-MHeF branch. I cannot reconcile these numbers with the plotted lifetimes, and the paper provides no code or data to resolve the discrepancy. If Equation (15) is simply mis-transcribed, the corrected version might restore the central claim, which is why my concrete test asks for the corrected equation or fitting script and a re-check of Fig. 12. But as the manuscript stands, the central numerical result is not reproducible, so a conditional acceptance would be inappropriate. This is a different and more basic issue than the reader's weakest assumption about metallicity extrapolation and uniform MH sampling, although those would matter even after Equation (15) is fixed.","tokens_in":33820,"tokens_out":9523,"duration_ms":84597,"concrete_test":"Re-implement Equation (15) with Table B.1 for the Fig. B.1 cases, especially MZAMS = 4.0, MHeMS = 0.479, MH = 0 and MZAMS = 1.2, MHeMS = 0.462, MH = 0, and compare with the plotted lifetimes of about 139 and 150 Myr. If the printed formula does not match, request the corrected age equation or fitting script from the authors; if the corrected version reproduces Fig. B.1, then re-run the Fig. 12 comparison to confirm that the claimed improvement in Kiel-diagram coverage persists.","verdict_should_be":"REJECT","load_bearing_attack":"The central claim rests on the H-shell prescription in Section 2.2.1, and specifically on Equation (15) with the coefficients in Table B.1. As printed, this equation is internally inconsistent. For the 'above MHeF' branch, substituting the Fig. B.1 case (MZAMS = 4.0, MHeMS = 0.479, MH = 0) gives base = 0.05161968 / (0.479 - 0.25380777) - 0.09981282 - 0.479 = -0.3496, raised to 1 / (0.00252778 * 0.479 + 424.26395881) = 0.00236. A negative base to a fractional power is undefined in real arithmetic. For the 'below MHeF' branch, taking MHeMS = 0.462 and MH = 0 gives base = 7.185 and exponent = 1 / (0.07264813 * 0.462 - 61.88868775) = -0.01617, so the computed lifetime is about 0.97 Myr, not the roughly 150 Myr shown in Fig. B.1 and discussed in the text. Thus the printed formula cannot be the function used to generate the model tracks in Fig. B.1 or the improvement in Fig. 12. The random uniform MH draw and the solar-metallicity extrapolation are additional legitimate concerns, but the failure of Equation (15) is more fundamental: without a working prescription, the paper's headline reconciliation with observed effective temperature and surface gravity is unsupported in the submitted manuscript.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper uses the rapid binary population synthesis code COMPAS to generate sdB populations across 162 configurations, varying common-envelope efficiency, metallicity, angular momentum loss, mass accretion efficiency, mass-transfer stability, and the helium-ignition threshold. The authors introduce an analytic post-processing prescription, fit to Bauer & Kupfer (2021) MESA models, for helium-main-sequence stars with hydrogen-rich shells, and apply it to convert COMPAS naked-HeMS tracks into sdB candidates in the Kiel diagram. They report that all studied parameters strongly affect sdB yields and distributions, identify formation channels (stable RLOF, one/two CE episodes, mergers), and claim that the H-shell prescription improves agreement with observed sdB samples. The paper explicitly states several limitations, including solar-metallicity-only fits, no extrapolation testing, and post-processing-only implementation.","tokens_in":34197,"tokens_out":7628,"duration_ms":67905,"significance":"If the H-shell prescription were correct and validated, the paper would be a useful contribution: it extends COMPAS to low-mass non-compact remnants for the first time, provides a systematic sensitivity map over 162 binary population configurations, and offers a fast analytic fitting scheme that other rapid BPS codes could adopt. The authors are transparent about limitations such as solar metallicity, no overshooting in the underlying MESA models, and the post-processing nature of the prescription. However, the reproducibility failure of Eq. (15) and the ad hoc random sampling of the hydrogen-shell mass mean that the paper's headline claim is not currently supported by the submitted manuscript.","major_comments":[{"comment":"Equation (15) as printed cannot reproduce the lifetimes shown in Fig. B.1. For the above-MHeF branch, taking the example MZAMS=4.0, MHeMS=0.479, MH=0 and interpreting the first term as (A1-MH)/(M+MH-A2) yields base = 0.05161968/(0.479-0.25380777) - 0.09981282 - 0.479 = -0.3496, which is negative and cannot be raised to the fractional exponent 1/(A4M+A5) in real arithmetic. For the below-MHeF branch with M=0.462 and MH=0, the same reading gives base = 7.185 and exponent = -0.01617, so tau_He is approximately 0.97 Myr, whereas Fig. B.1 and the text report roughly 150 Myr. As printed, the central prescription is therefore internally inconsistent and cannot be the function that generated the model tracks or the improved Kiel-diagram coverage in Fig. 12. The authors should provide a corrected, unambiguous equation, verified by a code snippet or a table of evaluated lifetimes, and re-run the affected post-processing.","section":"2.2.1, Eq. (15), Table B.1, Fig. B.1"},{"comment":"The claimed reconciliation with observations in Fig. 12 is obtained by drawing the hydrogen-rich shell mass MH uniformly from [0, 3e-3] Msun with no physical motivation, as stated in Section 3.1.6. Because the headline claim is that the H-shell prescription improves agreement with observed Teff and log g, the result is conditional on this ad hoc distribution. A different MH distribution, or one correlated with formation channel or metallicity, would change the spread and location of the synthetic candidates. The authors should quantify the sensitivity of the Kiel-diagram agreement to the assumed MH distribution, or justify the uniform draw from formation physics.","section":"3.1.6, Fig. 12"},{"comment":"The H-shell prescription is fit to Bauer & Kupfer (2021) models computed at solar metallicity with no overshooting, and the authors explicitly state that it has not been tested for extrapolation. Nevertheless, it is applied to all 162 runs, including Z=0.0012 and Z=0.03, and it enters the sdB-candidate selection through the Kiel-diagram cuts used throughout Section 3.1. The metallicity trends reported in Section 3.1.2 could therefore be partly influenced by extrapolating the H-shell fits beyond their validity range. Please either restrict the H-shell prescription to solar-metallicity runs, validate it against MESA models at sub- and super-solar Z, or demonstrate that the candidate selection is insensitive to the extrapolation.","section":"2.2.1, 3.1.2"},{"comment":"The central validation of the H-shell prescription is a visual comparison in Fig. 12; no quantitative metric is given for the claimed improvement, such as the fraction of synthetic candidates inside the observed box, a two-dimensional KS test, or chi-squared statistics on log g and Teff. Given the reproducibility problem in Eq. (15), this visual claim is currently unsupported. A quantitative comparison, ideally performed with the corrected prescription, should be added.","section":"3.1.6, Fig. 12"}],"minor_comments":[{"comment":"The expression for tau_He is ambiguous because the fraction is not typeset unambiguously in the text; please use explicit parentheses and define the order of operations.","section":"2.2.1, Eq. (15)"},{"comment":"The relation [Fe/H] approximately equal to log(Z/Zsun) ignores the dependence on hydrogen mass fraction and alpha enhancement; please state it as an approximation and give the resulting Z values explicitly.","section":"2.1.3, Eq. (7)"},{"comment":"The column header 'FA' is not defined in the table or caption; it appears to denote the mass accretion efficiency beta. Please rename and define it.","section":"Table A.1"},{"comment":"The word 'ellapsed' should be 'elapsed'.","section":"2.2.1"},{"comment":"The manuscript says that all data are available upon reasonable request, but no fitting code is released. For a prescription with dozens of fitted coefficients in Tables B.1-B.3, providing the fitting code or a machine-readable table of evaluated fits would greatly aid reproducibility.","section":"Data Availability"}],"recommendation":"major_revision","confidential_remarks":"The main technical issue is the printed Eq. (15), which must be corrected before the H-shell claim can be evaluated. If the authors can supply a corrected, reproducible prescription and quantitative validation, the paper could be suitable for publication in PASA. I also note that the random MH sampling and the solar-metallicity extrapolation are acknowledged in the text, but they still limit the strength of the headline claim. Because the fitting tables are complex and no code is released, I would strongly encourage the authors to make the fitting code and validation data available alongside the revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the parameter study: 162 COMPAS configurations across common-envelope efficiency, metallicity, angular momentum loss, accretion efficiency, and mass-transfer stability is a genuinely new mapping for sdB populations. The broad finding that all these inputs move the yields, and that the helium ignition threshold matters most, is plausible and internally consistent. The paper is also unusually honest about its limits: solar-metallicity MESA fits, no extrapolation, post-processing only.\n\nNow the soft spot, and it is big. Equation (15) as printed cannot reproduce the lifetimes that the paper itself shows in Fig. B.1. For the above-MHeF case (MZAMS=4.0, MHeMS=0.479, MH=0) the base is negative and the exponent is fractional, so the expression is undefined in real arithmetic. For the below-MHeF case (MHeMS=0.462) the same formula gives about 1 Myr, not the roughly 150 Myr shown. The H-shell prescription is the centerpiece of the paper's headline claim that theoretical log g and Teff now agree better with observations (Fig. 12). If the printed equation is wrong, readers cannot reproduce or evaluate that claim. The authors likely used a working fit, but the manuscript does not provide it, and no code or data are released.\n\nSmaller issues remain: shell masses are drawn from a flat distribution between 0 and 3e-3 Msun with no link to formation channel, and the prescription is applied outside the metallicity range it was fitted on. These are acknowledged, but they still weaken the reconciliation.\n\nBottom line: this is a paper for BPS practitioners working on hot subdwarfs. The parameter study deserves a referee, and the error in Eq. (15) is likely fixable with a corrected formula or a supplementary script. But as submitted, the central claim is unsupported. A revision that corrects the equation, provides the fitting code or the generated tracks, and adds a quantitative comparison metric (e.g., a KS test in the Kiel diagram) would be worth accepting.","headline":"A useful 162-run parameter study undermined by an unreproducible H-shell prescription: Equation (15) as printed contradicts the paper's own Fig. B.1.","tokens_in":34732,"tokens_out":6755,"would_cite":false,"duration_ms":56206,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A new analytic prescription for hydrogen-rich shells makes synthetic subdwarf B stars match observed surface gravity and temperature.","keywords":["hot subdwarf B stars","binary population synthesis","hydrogen-rich shells","common envelope evolution","mass transfer","helium ignition","Kiel diagram","horizontal branch stars"],"falsifier":"Measure or infer the hydrogen shell masses of a large sample of field subdwarf B stars, for example through asteroseismology or through the depth of the hydrogen Balmer lines, and compare the distribution to the uniform zero-to-$3\\times10^{-3}\\,M_\\odot$ assumption used here; a clear metallicity or channel dependence would break the prescription. Alternatively, rerun the same synthesis with metallicity-dependent shell masses and check whether the improved Kiel diagram agreement in Figure 12 survives.","tokens_in":33607,"feed_emoji":"🌟","tokens_out":5487,"duration_ms":50341,"temperature":0.7,"pith_summary":"This paper argues that the mismatch between synthetic and observed hot subdwarf B stars comes largely from a missing ingredient: the thin hydrogen-rich shell that survives on the surface of the stripped star. It presents an analytic prescription, fitted to detailed stellar models, that gives radius, luminosity, and helium-burning lifetime as functions of the helium-star mass, the shell mass, and whether helium ignited in a flash or smoothly. Applied to 162 binary population synthesis realizations, the prescription shifts the synthetic population's surface gravity and effective temperature so it overlaps the observed distribution in the Kiel diagram. The paper also reports that every tested binary-physics parameter, including common envelope efficiency, metallicity, angular momentum loss, mass transfer efficiency and stability, and the helium ignition threshold, changes the yields and properties of the resulting population, with the ignition threshold having the largest effect on total numbers.","feed_headline":"Thin hydrogen shells reconcile synthetic and observed hot subdwarf B stars","feed_subtitle":"A post-processing recipe for binary population synthesis lets rapid codes match observed surface gravity and temperature data.","key_machinery":"The central object is the analytic hydrogen-shell prescription: three fitted formulas, $\\tau_{\\rm He}(M, M_H)$, $R(t_r, M, M_H)$, and $L(t_r, M, M_H)$, with coefficients parameterized by the zero-age helium main sequence mass, the zero-age main sequence mass of the progenitor, and the hydrogen shell mass $M_H$. It carries the argument by converting naked helium stars, as classified by the Hurley et al. (2000) stellar types, into cooler and larger objects whose positions in the Kiel diagram can be compared to observations. The second load-bearing mechanism is the helium ignition threshold: stars that would be classified as helium white dwarfs are reclassified as subdwarf B candidates when their mass lies within 3 or 5 percent of the expected core mass at the tip of the red giant branch, which strongly boosts the canonical $\\sim0.47\\,M_\\odot$ peak.","core_discovery":"On the paper's own terms, the discovery is a workable shortcut: rapid binary population synthesis codes treat subdwarf B progenitors as naked helium stars, and that assumption places the synthetic population at hotter temperatures and higher surface gravity than observed. The paper derives a post-processing prescription, based on the MESA models of Bauer & Kupfer (2021), that adds a hydrogen-rich outer shell of up to $3\\times10^{-3}\\,M_\\odot$ and recomputes radius, luminosity, and lifetime through three fitted relations, with separate coefficient sets for progenitors that ignite helium in a flash and those that ignite it smoothly. When this prescription is applied to the helium main sequence stars produced by COMPAS, the synthetic sample in the Kiel diagram spreads out to cover the observed subdwarf B box. The study's claim is therefore that hydrogen shells, not just the choice of binary parameters, are what reconcile population synthesis predictions with the observed $\\log g$ and $T_{\\rm eff}$ distributions, and that the prescription can be ported to any rapid binary population synthesis output.","pith_inferences":["If real hydrogen shell masses are correlated with formation channel or progenitor metallicity, the uniform random sampling between 0 and $3\\times10^{-3}\\,M_\\odot$ used here would bias the inferred population, and the Kiel coverage would change.","The lack of coverage at low surface gravity in Figure 12 suggests that envelopes heavier than $3\\times10^{-3}\\,M_\\odot$ exist in nature; a targeted asteroseismic or spectroscopic survey of subdwarf B shell masses would test this directly.","The same fitting strategy could be transferred to subdwarf O stars or other stripped stars, where hydrogen shells also regulate temperature and radius.","A direct test is to compare the predicted fraction of subdwarf B plus neutron star systems from the early common envelope channel against the accumulating sample of wide subdwarf B binaries with neutron star companions."],"forward_implications":["If the prescription is correct, any rapid binary population synthesis code can post-process its helium stars to predict subdwarf B surface gravities and effective temperatures without running detailed stellar models.","The parameter study maps which uncertain binary physics matters most: common envelope efficiency and the helium ignition threshold dominate yields, while mass-loss geometry shifts period and mass distributions.","Predicted subdwarf B populations can now be compared to observed catalogs in the Kiel diagram, enabling future current-day Galactic population synthesis with a chosen configuration.","The early common envelope plus stable mass transfer channel is predicted to produce most subdwarf B plus neutron star systems with periods longer than one day, a testable prediction.","Hydrogen shells change candidate counts by only about 1 percent, but they change the observable properties of the population substantially, so counts and distributions should be treated separately."],"supporting_citations":[{"why":"Provides the MESA models with and without hydrogen-rich shells that the new prescription is fitted to.","marker":"Bauer & Kupfer (2021)"},{"why":"Defines the stellar types and naked helium star evolution scheme that COMPAS inherits and that the prescription extends.","marker":"Hurley et al. (2000)"},{"why":"Establishes the binary formation channels for subdwarf B stars and the 5 percent helium ignition threshold used in the population comparison.","marker":"Han et al. (2002)"},{"why":"Documents the COMPAS code and its default binary physics that all 162 configurations build on.","marker":"Riley et al. (2022)"},{"why":"Supplies the observational sample whose Kiel diagram envelope defines the candidate selection box.","marker":"Culpan et al. (2022)"},{"why":"Provides the alternative critical mass ratio prescription tested against the default stability criterion.","marker":"Ge et al. (2020)"},{"why":"Gives the detailed helium ignition mass ranges used to justify the 3 and 5 percent cutoffs.","marker":"Arancibia-Rojas et al. (2024)"}],"fun_headline_variants":["Hydrogen shells align synthetic hot subdwarfs with data","New recipe reconciles subdwarf models and observations","Missing hydrogen shells were the subdwarf model's blind spot","Hot subdwarf simulations get a hydrogen-shell fix","Prescription adds hydrogen shells to subdwarf models for better fits"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The prescription assumes that detailed models built at solar metallicity, without overshooting, for a restricted range of masses can be applied to binary population synthesis helium stars of any metallicity and history, and that the hydrogen shell mass is uniformly distributed between zero and $3\\times10^{-3}\\,M_\\odot$.","fun_headline_variants_meta":{"raw":{"variants":["Hydrogen shells align synthetic hot subdwarfs with data","New recipe reconciles subdwarf models and observations","Missing hydrogen shells were the subdwarf model's blind spot","Hot subdwarf simulations get a hydrogen-shell fix","Prescription adds hydrogen shells to subdwarf models for better fits"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000991,"raw_usage":{"total_tokens":4218,"prompt_tokens":983,"completion_tokens":3235,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":599,"completion_tokens_details":{"reasoning_tokens":3152}},"tokens_in":599,"tokens_out":3235,"duration_ms":21979,"temperature":1.0,"reasoning_tokens":3152,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T14:38:52.913608+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure or infer the hydrogen shell masses of a large sample of field subdwarf B stars, for example through asteroseismology or through the depth of the hydrogen Balmer lines, and compare the distribution to the uniform zero-to-$3\\times10^{-3}\\,M_\\odot$ assumption used here; a clear metallicity or channel dependence would break the prescription. Alternatively, rerun the same synthesis with metallicity-dependent shell masses and check whether the improved Kiel diagram agreement in Figure 12 survives.","supporting_citations":[],"review_version":1}