{"id":"5f07b7f1-8762-46ff-ac61-fd41da069a84","arxiv_id":"2505.22201","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A new grid of 2,700 accretion models shows that weak Ba, CH, and CEMP-s stars are well reproduced with AGB donors of 2-3 solar masses, while strong Ba stars require accreted masses beyond the grid's 0.5 solar mass limit, which fails to match their observed mass distribution.","lead":"This paper computes 2,700 binary accretion models to infer the masses of the AGB donor stars and the amounts of accreted material that produced observed barium, CH, and carbon-enhanced metal-poor s-process stars. Generalist readers should care because the results constrain how heavy elements made by slow neutron capture are transferred in binary star systems.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Strong-Ba conclusion extrapolates beyond the grid: with accreted mass capped at 0.50 Msun, the 'must accrete more' claim is not established, and untested higher-accretion models could also resolve the mass-distribution mismatch.","rationale":"Good-faith reading: the paper delivers a useful, clearly described grid of accretion models and applies it to several related stellar classes. The qualitative results for weak Ba, CH, and CEMP-s stars (2-3 Msun AGB donors, modest accretion) are plausible and consistent with prior literature; those parts do not hinge on the grid boundary. The place where the argument is least secure is the headline strong-Ba conclusion. The data shown in Figures 4-5 and Section 4 indicate that the maximum-likelihood models cluster at the largest accretion mass available in the grid; the abstract then states that strong Ba stars 'must' accrete more than 0.50 Msun. From a grid with a hard cap at 0.50, the correct inference is that the fit is truncated, not that the true value lies strictly above 0.50. The paper's own Section 6 sentence about higher accretion masses is effectively a limitation statement, and it counts as evidence against the central claim. Thermohaline mixing is a legitimate concern for CEMP-s and CH abundance dilution, as the reader notes, but it is less central: the strong-Ba conclusion would be vulnerable even if thermohaline mixing were included, because the tested accreted-mass range ends at 0.50. The proposed check is direct: recompute with accreted masses up to 1.0 Msun. If the best fits stay at 0.50, the conclusion is empirically supported; if not, the abstract and conclusions need revision. The paper is transparent about the grid limits, and I see no reason to question the authors' integrity. The recommendation remains CONDITIONAL, matching the reader's verdict, so the verdict is unchanged.","tokens_in":19109,"tokens_out":5695,"duration_ms":53414,"concrete_test":"Extend the Section 3 grid with accreted masses of 0.60, 0.70, 0.80, and 1.00 Msun for secondary initial masses of 0.40-2.00 Msun at the metallicities of the strong-Ba sample (Z about 0.006-0.010), using the same STARS plus FRUITY prescriptions, and re-run the maximum-likelihood fit for the strong Ba stars. Compare the resulting best-fit accreted-mass histogram and final-mass distribution with the Escorza et al. (2017) double-Gaussian (2.5 +/- 0.18 and 3.0 +/- 1 Msun). If the best-fit accreted masses remain pinned at 0.50 Msun despite the wider grid, the 'must accrete more than 0.50 Msun' claim survives; if they move above 0.50 and the final masses shift upward toward 2-3 Msun, the central claim is a grid-boundary artifact.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3 sets the accretion-mass grid to 0.05, 0.10, 0.20, 0.30, 0.40, and 0.50 Msun, so 0.50 is a hard upper boundary. Section 4 reports that the strong-Ba best fits are 'strongly peaked at high accretion masses, with most fits showing 0.50 Msun.' A maximum-likelihood estimate that lands on the edge of a discrete grid is not evidence that the true value is at or above that edge; it is evidence that the grid cannot represent larger values. The abstract and conclusions convert this boundary effect into the central claim that strong Ba stars 'must accrete more than 0.50 Msun' and that this limit 'is unable to reproduce the observed mass distribution of strong Ba stars.' The paper explicitly concedes the point in Section 6: 'Further modeling considering higher accretion masses may reveal a scenario that better describes the final mass distribution of the strong Ba stars.' Since the final mass is the initial mass plus the accreted mass, increasing the maximum accreted mass could permit larger final masses, e.g. near the 2.5 Msun peak of Escorza et al. (2017), while still leaving enough accreted s-process material on the surface to match the observed high [hs/Fe]. Without those models, the claim that the standard binary scenario fails for strong Ba stars is an extrapolation beyond the computed parameter space rather than a demonstrated result.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a grid of about 2700 binary accretion models computed with the STARS code, using FRUITY AGB yields, and compares them via a maximum-likelihood chi-squared analysis to observed Ba, CH, and CEMP-s stars. The authors report consistent AGB donor masses around 2-3 Msun, moderate accreted masses (≤0.5 Msun) for weak Ba stars, large accreted masses (≥0.5 Msun) for strong Ba stars, and low accreted masses (~0.1 Msun) for CH and CEMP-s stars. They also compare derived stellar masses with the observed mass distribution and discuss orbital properties and accretion efficiencies. The central conclusion is that strong Ba stars must accrete more than 0.5 Msun to explain their abundances, and that at this grid boundary the model cannot reproduce the observed mass distribution of strong Ba stars.","tokens_in":19330,"tokens_out":4682,"duration_ms":48343,"significance":"If established, the claim that strong Ba stars require accretion above 0.5 Msun would identify a concrete tension between standard binary mass-transfer models and observations, motivating new modeling efforts and possibly new mass-transfer mechanisms. The paper has notable strengths: a large and reusable model grid, use of independent FRUITY yields, explicit treatment of dilution and first dredge-up, and successful reproduction of the mass and abundance distributions for weak Ba, CH, and CEMP-s stars. The individual fits (e.g., PV UMa, HD 123949, CS 29512-073) are informative. However, the headline conclusion about strong Ba stars rests on an extrapolation beyond the computed grid, since the maximum accreted mass in the grid is 0.5 Msun. The paper also concedes in Section 6 that higher-accretion models may reveal a better scenario, which undercuts the strength of the abstract's claim. The neglect of thermohaline mixing and the acknowledged correlation problem in the chi-squared comparison add further caveats, especially for the metal-poor populations.","major_comments":[{"comment":"The central claim that strong Ba stars 'must accrete more than 0.50 Msun' is not established by the presented grid. The accreted mass is capped at 0.50 Msun by construction in Section 3, and the best-fit models piling at this boundary only demonstrate that the grid cannot represent larger values; a maximum-likelihood estimate at a discrete grid edge does not imply the true value lies at or above that edge. Moreover, since final mass is the sum of initial mass plus accreted mass, allowing larger Delta M could move the final masses upward, potentially resolving the reported mismatch with the observed mass distribution of Escorza et al. (2017). The paper's own Section 6 concedes that 'further modeling considering higher accretion masses may reveal a scenario that better describes the final mass distribution of the strong Ba stars,' which is in direct tension with the definitive wording in the abstract. Please either extend the grid to larger accreted masses for at least the strong Ba stars, or reframe the conclusion as a provisional extrapolation rather than a demonstrated failure of the standard scenario.","section":"Section 3, Section 4.1, Fig. 5, Abstract/Conclusions"},{"comment":"The neglect of thermohaline mixing is acknowledged but its impact on the inferred parameters is not quantified. The paper states in Section 5.1 that without thermohaline mixing, surface abundances remain nearly constant until first dredge-up, making differentiation between best-fit models difficult. For the metal-poor CEMP-s and CH stars, thermohaline mixing on the main sequence can significantly dilute light elements (Stancliffe et al. 2007), yet the low accreted masses (~0.05-0.1 Msun) are derived under the assumption of negligible thermohaline mixing. This could systematically bias the inferred accretion masses for these populations. Please provide a quantitative assessment of how the best-fit parameters would change if thermohaline mixing (or a range of mixing efficiencies) were included, or soften the population-level conclusions for CEMP-s and CH stars.","section":"Section 3, Section 5.1"},{"comment":"The chi-squared comparison assumes that the fitted quantities are uncorrelated, but the paper explicitly notes that the heavy-element abundances are highly correlated with one another and with the surface parameters. This is not a minor caveat: it affects the interpretation of the differences in chi-squared between competing models and could change the ranking of best-fit models, especially when the fits differ only subtly in abundance pattern. Please discuss the likely direction and magnitude of this effect, or implement a covariance-aware likelihood to demonstrate that the main results—particularly the strong-Ba preference for the maximum accreted mass—are robust to the treatment of correlated abundances.","section":"Section 3.1"}],"minor_comments":[{"comment":"The sentence 'The observed temperature The two best fit models...' is incomplete and should be rephrased.","section":"Section 4.1, PV UMa paragraph"},{"comment":"The variable 'Xorignial' appears to be a typo for 'Xoriginal'.","section":"Section 3, Eq. (1)"},{"comment":"The sentence beginning 'Stancliffe et al. (2007) found that ther-mohaline mixing is on the main sequence is effective...' contains a grammatical error; the phrase 'is on the main sequence is effective' should be corrected.","section":"Section 5.1"},{"comment":"Please add the sample sizes for each population to the caption, since the CH star sample is stated to be small and this affects the interpretation of the histograms.","section":"Figure 5"},{"comment":"The abstract states 'macc≥0.5 Msun' while Section 4.1 reports that 'most fits showing 0.5 Msun'; please ensure consistent phrasing across the paper","section":"Abstract and Section 4.1"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of A&A and represents a substantial modeling effort. The main issue is the extrapolative nature of the strong-Ba conclusion; the authors should either extend the grid or clearly temper the claim. The thermohaline mixing and chi-squared caveats are acknowledged in the manuscript but deserve more quantitative treatment, especially because they affect the reliability of the metal-poor results. I would not reject the paper, but the central claim needs to be placed on firmer ground before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nTwo things before you read. The paper builds a genuinely new 2700-model grid covering Ba, CH, and CEMP-s stars, and the headline claim—that strong Ba stars must accrete more than 0.5 solar masses and that the standard binary scenario then fails to match their mass distribution—is an extrapolation from the edge of that grid. The grid is valuable; the strong-Ba claim is not yet established.\n\nWhat it does well: it extends Stancliffe (2021) from a fixed metallicity and final mass to a wider grid across four populations, applies a clean maximum-likelihood comparison to published abundances, and recovers a consistent AGB donor mass range of 2–3 solar masses. The qualitative results for weak Ba, CH, and CEMP-s stars (modest accretion, roughly 0.05–0.4 solar masses) are plausible and consistent with prior work. The abundance fitting is not circular—observed abundances are compared against independent FRUITY yields. That is real progress.\n\nThe soft spots are real but concentrated. The strong-Ba fits pile at the hard upper boundary of the accretion grid (0.5 solar masses). A maximum-likelihood estimate landing on the grid edge shows the grid cannot represent larger values, not that the true values lie at or beyond the edge. The paper concedes as much in Section 6: 'Further modeling considering higher accretion masses may reveal a scenario that better describes the final mass distribution of the strong Ba stars.' That sentence undercuts the 'must accrete more' language in the abstract. Higher-accretion models could plausibly raise final masses toward the 2.5-solar-mass peak of Escorza et al. (2017) while still matching the high surface [hs/Fe], because the extra accreted material offsets the stronger dilution in a more massive convective envelope. Without those models, the claimed discrepancy is suggestive, not demonstrated.\n\nOther issues are minor to moderate. The fit excludes Mg, alpha elements, and Nb; the rationale (observed values systematically exceed AGB yield predictions) is defensible, but it narrows the comparison to C and s-process elements. Chi-squared is used despite the paper's own acknowledgement that heavy-element abundances are strongly correlated; that weakens the formal ranking, though it is unlikely to change the qualitative picture. Thermohaline mixing is omitted, which is a genuine concern for the low-metallicity CEMP-s and CH stars where accretion happens on the main sequence and is followed by first dredge-up; the paper's discussion shows the effect is metallicity-dependent and may not be negligible.\n\nWho is this for: stellar modelers and observers working on s-process enrichment and post-mass-transfer binaries. The grid and the population-level constraints deserve a serious referee, who should insist on extended accretion masses, a thermohaline treatment for the metal-poor classes, and a public release of the grid. I would not desk-reject.","headline":"Valuable new grid for post-accretion binaries whose headline strong-Ba claim overreaches the computed parameter space.","tokens_in":19975,"tokens_out":4718,"would_cite":true,"duration_ms":46494,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper argues that reproducing strong barium stars' abundances requires each star to accrete more than 0.5 solar masses, and at that limit the standard binary accretion models fail to reproduce their observed mass distribution.","keywords":["barium stars","CH stars","CEMP-s stars","s-process nucleosynthesis","AGB mass transfer","binary accretion","stellar evolution models","dilution factor"],"falsifier":"Compute the same grid with thermohaline mixing included and check whether the inferred accreted masses for CEMP-s and CH stars drop below about 0.1 solar masses and whether strong Ba stars can be fit with final masses above 2 solar masses. Alternatively, measure C, N, and O on the main sequence or subgiant branch of a sample of CEMP-s stars: if their surface C/N ratio declines before first dredge-up, thermohaline dilution is acting and the no-mixing models are false.","tokens_in":18802,"feed_emoji":"⭐","tokens_out":7827,"duration_ms":81831,"temperature":0.7,"pith_summary":"The paper tries to pin down the progenitors of barium, CH, and carbon-enhanced metal-poor s-process (CEMP-s) stars, binaries whose surfaces carry heavy elements made by the slow neutron capture process in an evolved AGB companion. Using a new grid of 2,700 accretion models spanning metallicities from [Fe/H] -2.15 to -0.15 and comparing them star-by-star to observed abundances by maximum likelihood, it concludes that AGB donors of about 2-3 solar masses can explain all four populations. The reproduction works well for weak Ba, CH, and CEMP-s stars, but fails for strong Ba stars: matching their high s-process abundances requires each to accrete more than 0.50 solar masses, and in that regime the models cannot reproduce the observed mass distribution of strong Ba stars. If this is right, strong Ba stars are the place to look for missing physics in how polluted binaries accrete.","feed_headline":"Strong barium stars demand giant accretions","feed_subtitle":"New 2,700-model grid reproduces weak Ba, CH, and CEMP-s stars but fails the strongest class.","key_machinery":"The argument is carried by a grid of about 2,700 binary evolution models in which each system is specified by metallicity, AGB donor mass, initial accretor mass, and accreted mass (0.05-0.50 solar masses). Accreted AGB ejecta, taken from published s-process yield tables, are deposited on the surface and tracked with a tracer composition so that the surface abundance at any time is the mass-weighted mix of accreted material and original stellar material. The crucial mixing quantity is the dilution factor $d = M_{\\mathrm{acc}}/M_{\\mathrm{mix}}$, the ratio of accreted mass to the mass of the convective envelope after the first dredge-up; it controls how much of the observed s-process enhancement survives. Comparing each model at each timestep to observed effective temperature, surface gravity, metallicity, carbon, and heavy-element abundances via a chi-square maximum-likelihood test selects the best-fit progenitor for each star.","core_discovery":"The central claim is that the same basic scenario, a 2-3 solar mass AGB donor transferring mass to a lower-mass companion, underlies weak Ba, strong Ba, CH, and CEMP-s stars, but the strength of the required transfer separates them. Weak Ba stars are best fit by moderate accretion (up to 0.5 solar masses) onto a roughly 2.0-2.5 solar mass star; CH and CEMP-s stars by small accretions (about 0.1 solar masses) onto roughly 1.0 solar mass stars; strong Ba stars by large accretions (at least 0.5 solar masses) onto roughly 1.0-2.0 solar mass stars. The paper's key negative result is that strong Ba stars, which require the highest accretion masses, cannot be made consistent with the observed mass distribution: the models push their final masses down near 1 solar mass, about a full solar mass below the peak found in earlier observational studies, because more massive accretors dilute the accreted s-process material more thoroughly. The paper therefore states that in the high-accretion limit it is unable to reproduce the observed mass distribution of strong Ba stars.","pith_inferences":["This reading suggests a testable extension: recomputing the grid with thermohaline mixing switched on for metal-poor accretors would shift the inferred small accreted masses for CEMP-s and CH stars, since the paper itself notes that thermohaline mixing can dilute light elements on the main sequence.","I infer that the strong-Ba mass discrepancy could also be read as evidence that the observed strong Ba sample is biased toward low-mass giants, or that a missing mixing process such as rotationally induced mixing changes how much dilution high-mass accretors experience.","A full binary population synthesis that feeds these best-fit progenitor parameters into initial binary distributions could predict how many strong versus weak Ba stars should exist, directly testing whether the required high-accretion channel is actually populated.","The authors' planned follow-up on long-period systems is a natural place to test whether wind Roche-lobe overflow or circumbinary discs can supply the angular momentum budget needed for large accretions."],"forward_implications":["If the grid is right, AGB donors of 2-3 solar masses are common progenitors of all four classes, which connects their chemical enrichment to the same nucleosynthetic production site.","The mass distributions of weak Ba, CH, and CEMP-s stars are reproduced, so current convective-mixing physics appears sufficient for those systems.","Strong Ba stars require accretion of at least 0.5 solar masses, and wind mass-transfer models cannot deliver that much material, so their formation likely needs a different mass-transfer channel.","The recovered final masses of strong Ba stars sit near 1.0 solar mass, about one solar mass below observational estimates, so either the observed strong Ba sample or the treatment of dilution is incomplete.","Strong Ba systems require accretion efficiencies around 25 percent or higher, above the wind-accretion efficiencies predicted by three-dimensional hydrodynamical models."],"supporting_citations":[{"why":"Supplies the earlier fixed-metallicity barium-star model grid and the dilution-factor definition that this paper generalizes.","marker":"Stancliffe (2021)"},{"why":"Provides the observed strong and weak Ba star sample and abundances used for the fits.","marker":"de Castro et al. (2016)"},{"why":"Supplies the homogeneous Ba, CH, and CEMP-s sample with stellar parameters and abundances.","marker":"Dimoff et al. (2024)"},{"why":"Gives the observed Ba star mass distribution that the strong-Ba models fail to match.","marker":"Escorza et al. (2017)"},{"why":"Provides the AGB s-process yields that set the composition of the accreted material.","marker":"Cristallo et al. (2015)"},{"why":"Establishes how dilution and thermohaline mixing affect CEMP-s surface abundances, the mixing framework this paper builds on.","marker":"Matrozis & Stancliffe (2016)"},{"why":"Gives three-dimensional wind-accretion efficiencies used as the physical comparison for inferred accretion efficiencies.","marker":"Liu et al. (2017)"},{"why":"Supplies the critical-rotation accretion limit used to argue that wind accretion cannot deliver the Ba-star accretion masses.","marker":"Matrozis et al. (2017)"},{"why":"Is the paper's own citation for main-sequence thermohaline dilution in metal-poor stars, which the models neglect.","marker":"Stancliffe et al. (2007)"}],"fun_headline_variants":["Strong Ba stars need big accretions but models miss their masses","Accretion paradox: heavy s-process delivery fails for strong Ba stars","Massive accretion can't explain strong barium stars' masses","Model grid reproduces weak Ba, CH, CEMP-s but not strong Ba"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result stands on the assumption that thermohaline mixing is negligible in the accretors, so surface abundances stay unchanged from the end of accretion until first dredge-up; if thermohaline mixing is significant, especially for metal-poor CEMP-s and CH stars, the inferred accreted masses would be systematically wrong.","fun_headline_variants_meta":{"raw":{"variants":["Strong Ba stars need big accretions but models miss their masses","Accretion paradox: heavy s-process delivery fails for strong Ba stars","Massive accretion can't explain strong barium stars' masses","Model grid reproduces weak Ba, CH, CEMP-s but not strong Ba"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000534,"raw_usage":{"total_tokens":2696,"prompt_tokens":1202,"completion_tokens":1494,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":818,"completion_tokens_details":{"reasoning_tokens":1418}},"tokens_in":818,"tokens_out":1494,"duration_ms":11485,"temperature":1.0,"reasoning_tokens":1418,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T13:12:34.725419+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the same grid with thermohaline mixing included and check whether the inferred accreted masses for CEMP-s and CH stars drop below about 0.1 solar masses and whether strong Ba stars can be fit with final masses above 2 solar masses. Alternatively, measure C, N, and O on the main sequence or subgiant branch of a sample of CEMP-s stars: if their surface C/N ratio declines before first dredge-up, thermohaline dilution is acting and the no-mixing models are false.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the earlier fixed-metallicity barium-star model grid and the dilution-factor definition that this paper generalizes."},{"cited_title":"J., Hansen, C","cited_arxiv_id":null,"evidence_quote":"Supplies the homogeneous Ba, CH, and CEMP-s sample with stellar parameters and abundances."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the observed Ba star mass distribution that the strong-Ba models fail to match."},{"cited_title":"2015, ApJS, 219, 40","cited_arxiv_id":null,"evidence_quote":"Provides the AGB s-process yields that set the composition of the accreted material."},{"cited_title":"& Stancliffe, R","cited_arxiv_id":null,"evidence_quote":"Establishes how dilution and thermohaline mixing affect CEMP-s surface abundances, the mixing framework this paper builds on."},{"cited_title":"J., Abate, C., & Matrozis, E","cited_arxiv_id":null,"evidence_quote":"Gives three-dimensional wind-accretion efficiencies used as the physical comparison for inferred accretion efficiencies."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the critical-rotation accretion limit used to argue that wind accretion cannot deliver the Ba-star accretion masses."}],"review_version":1}