{"id":"0f711d59-b938-4c55-8fbb-4db6341acad9","arxiv_id":"2412.01025","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Binary evolution reduces C and s-process yields from low- and intermediate-mass stars by roughly 20-25% at a binary fraction of 0.7, while N and O yields remain nearly unchanged.","lead":"This paper models how having a stellar companion changes the chemical elements that Sun-like and slightly heavier stars release as they die. It finds that a realistic population with 70% binaries ejects about 20-25% less carbon and heavy s-process elements than a population of only single stars.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed 20-25% reduction in s-process yields relies on extrapolating K16 intershell abundances to stars below 1.5 Msun, an extrapolation the authors tie to their overproduction of high-[Ce/Y] Ba stars; this systematic is not covered by the paper's sensitivity tests.","rationale":"The reader's weakest_assumption was the calibration of TDU parameters against a CSLF assumed to be dominated by intrinsic carbon stars. That is a real concern, but the paper's own sensitivity analysis shows that varying the TDU parameters around the best fit changes the single-to-binary yield ratios very little (C, N, and Ba ratios of 0.76, 0.93, and 0.67 at fb=1.0 remain almost identical across TDU, wind, and mass-transfer variations). Thus the relative reduction, which is the headline, appears fairly robust to the calibration uncertainty. The more serious soft spot is the s-process intershell extrapolation for masses below 1.5 Msun. This is not a mere parameter variation; it is an acknowledged physical mismatch with the K16 nucleosynthesis models and with the observed [Ce/Y] distribution of Ba stars, and it affects the very mass range that dominates the Sr, Ba, and Pb yields. The paper reports the mismatch but does not test how the population-level s-process reduction changes when a realistic low-mass intershell composition is used. Because the central claim bundles C with s-process elements in a single 20-25% reduction statement, the s-process portion is less secure than the C portion. The paper remains publishable as a conditional contribution: the methodology is sound, the code is available, and the authors are transparent about limitations, but the s-process yields need a targeted test before the abstract's quantitative claim can be taken at face value. This does not change the reader's CONDITIONAL verdict; it sharpens the condition.","tokens_in":34601,"tokens_out":6338,"duration_ms":64290,"concrete_test":"Recompute the Sr, Ba, and Pb population yields for binary fractions 0.0 and 0.7 using a modified heavy-element intershell table in which stars with initial mass below 1.5 Msun produce no s-process elements (or are assigned the systematically lower yields of a model with convective 13C burning, e.g., from FRUITY at 1.2-1.3 Msun). If the binary-to-single yield ratios for Ba and Pb at fb=0.7 change by more than about 5 percentage points from the 0.75 and 0.77 values in Table 4, the s-process part of the headline claim is not robust to this extrapolation; if the ratios remain within 5%, the concern is mitigated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim for Sr, Ba, and Pb in Table 4 depends on s-process yields from AGB stars with initial masses as low as about 1.2 Msun (Section 3.1), yet Section 2.2 states that for stars below 1.5 Msun the heavy-element intershell abundances are taken from the 1.5 Msun K16 model. K16 does not produce a 13C pocket at those masses, and Section 5.2 explicitly notes that a 1.2 Msun star likely burns 13C convectively, yielding fewer neutrons and less s-process production than the 1.5 Msun model. The authors also report that 63% of the systems overproducing [Ce/Y] in their Ba-star comparison have AGB companions below 1.5 Msun, directly implicating this extrapolation. Because stars with masses below about 2 Msun eject roughly 80% of the Ba in the population, the absolute s-process yields and, potentially, the single-to-binary yield ratios are biased by this abundance-table assumption. The sensitivity analysis in Section 5.1 varies TDU efficiency parameters, the common-envelope parameter, mass-loss, and mass-transfer prescriptions, but it never varies the intershell abundances for low-mass stars, so the robustness of the s-process reduction to this specific systematic is untested. In contrast, the carbon yield reduction is largely independent of this s-process extrapolation and is better supported by the model variations.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper upgrades the binary population synthesis code binary_c by introducing s-process intershell abundances fitted to the Karakas and Lugaro (2016) models and calibrating the third dredge-up parameters (Delta_Mc,min = -0.13 Msun, lambda_min = 0.45) against the Galactic carbon-star luminosity function of Abia et al. (2022). It then computes population-weighted yields of C, N, O, Sr, Ba, and Pb for a grid of 1000 single-star and 640,000 binary models at solar metallicity, as a function of binary fraction. The central result is that a population with binary fraction 0.7 ejects about 18% less C and about 23-25% less Sr, Ba, and Pb than a single-star-only population, while N and O yields change by only a few percent. The paper also compares predicted Ba-star surface abundances with observations, finding an overproduction of high-[Ce/Y] stars, and reports rare predicted channels such as massive Ba stars and O-rich naked He stars.","tokens_in":34926,"tokens_out":6442,"duration_ms":63309,"significance":"If the central result is robust, it is an important contribution to Galactic chemical evolution: it quantifies a 20-25% reduction in AGB yields for a realistic binary fraction and identifies which elements are affected. The study has notable strengths: a large and systematic model grid, an explicit population-weighting scheme, a public code version, a calibration against observed carbon stars, and a sensitivity analysis that varies several binary and AGB parameters. The relative single-to-binary yield ratios for C, N, and Ba are shown to be stable under the tested variations of alpha_CE, third dredge-up parameters, mass-loss prescriptions, and mass-transfer treatments. The main weakness is that the s-process part of the central claim rests on an untested extrapolation of intershell abundances to stars below 1.5 Msun, which the paper itself identifies as problematic in Section 5.2.","major_comments":[{"comment":"The Sr, Ba, and Pb yield ratios in Table 4 are not robust to the assumed intershell abundances for AGB stars below 1.5 Msun. Section 2.2 states that for stars of mass <1.5 Msun the heavy-element intershell abundances are taken from the 1.5 Msun K16 model, while Section 5.2 notes that a 1.2 Msun star likely burns 13C convectively, yielding fewer neutrons than the 1.5 Msun model, and that 63% of the overproduced high-[Ce/Y] Ba stars have AGB companions below 1.5 Msun. The sensitivity tests in Section 5.1 vary alpha_CE, third dredge-up parameters, mass-loss prescriptions, RLOF, and WRLOF, but they never vary the low-mass intershell abundances. Since Section 5.1 also reports that stars below about 2 Msun eject roughly 80% of the Ba, this systematic can affect both the absolute s-process yields and the single-to-binary ratios. A concrete test would be to recompute the population ratios after setting the below-1.5-Msun heavy intershell abundances to zero or to a lower-neutron-density model; without such a test, the s-process part of the central claim is not yet fully supported.","section":"§5.2 and Table 4"},{"comment":"The abstract and the conclusions state that at a binary fraction of 0.7 there is a 20-25% reduction in C, Sr, Ba, and Pb yields, but Table 3 gives a C yield ratio of 0.82 at fb=0.7, i.e., an 18% reduction, and the 24% C reduction quoted in Table 3 and Section 5.1 corresponds to fb=1.0 rather than fb=0.7. The text should be corrected to distinguish the C reduction (about 18% at fb=0.7) from the heavier s-process reductions (23-25% at fb=0.7), or explicitly state the fb=1.0 values when quoting the larger reduction.","section":"Abstract, §6, and Table 3"},{"comment":"The third dredge-up calibration assumes that the observed carbon stars used for the CSLF fit are mostly intrinsic AGB stars rather than extrinsic binary-polluted stars. The paper explicitly acknowledges this assumption and filters out low-luminosity GB carbon stars, but it does not quantify the sensitivity of the fitted parameters, and therefore of the absolute yields, to a possible extrinsic AGB contamination. A useful check would be to construct a CSLF from the binary grid with a range of assumed extrinsic fractions and determine how much the best-fit Delta_Mc,min and lambda_min shift; this would test whether the carbon yield reduction is robust to the central calibration assumption.","section":"§2.1 and §3.1"}],"minor_comments":[{"comment":"The phrase 'approximately 50% of the abundances of the elements heavier than iron' is standard but should be attributed to a specific reference in the introduction, where it appears later; the abstract currently reads as a bald claim.","section":"Abstract"},{"comment":"The notation [X/Y] = log10(X/Y)_star - log10(X/Y)_sun is clear, but the solar ratio is written as 'log10(X/Y)_sun' only in the text; the equation would benefit from an explicit definition of the solar standard as Lodders (2003) at the point of use.","section":"Equation (6)"},{"comment":"The right-hand panel of Figure 14 does not show the low-[Ce/Y] tail of the observed sample; the authors state 3% of Ba stars have [Ce/Y] < -0.4, but this is not visible in the plotted range, making the claimed discrepancy harder to assess.","section":"Figure 14"},{"comment":"The sentence 'The most extreme O producers ... similar to an R Coronae Borealis star but O-rich instead of C-rich' introduces a class of objects not previously predicted; it would help to state explicitly in the text that no such object has yet been observed, as is done later.","section":"Section 3.3.1"}],"recommendation":"major_revision","confidential_remarks":"This is a solid population-synthesis study with a transparent code release and a careful parameter-sensitivity analysis. The main obstacle to acceptance is the untested extrapolation of K16 s-process abundances to stars below 1.5 Msun, which directly affects the Sr/Ba/Pb part of the headline result. If the authors can show that the single-to-binary yield ratios are insensitive to this extrapolation, or if they revise the claims to exclude or qualify the s-process reduction, the paper would be suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this paper if you care about binary effects on AGB yields. The headline result — at binary fraction 0.7, C and s-process yields drop ~20–25% relative to single stars — is backed by a large grid and a genuine attempt to anchor the model to observations. What's new is the update of binary_c's intershell abundances to Karakas & Lugaro (2016), the calibration of the third dredge-up to the carbon star luminosity function, and the population-level yields for C,N,O,Sr,Ba,Pb as a function of binary fraction. The Ba star comparison is a good sanity check, even if it reveals problems.\n\nThe paper is honest about its soft spots. The most serious one is the s-process extrapolation to stars below 1.5 Msun. For these, the heavy-element intershell abundances are taken from the 1.5 Msun K16 model, which doesn't produce a 13C pocket at those masses. The authors note in Sec 5.2 that a 1.2 Msun star likely burns 13C convectively, yielding fewer neutrons. They also report that 63% of the systems overproducing [Ce/Y] in the Ba star comparison have AGB companions below 1.5 Msun. The sensitivity analysis covers TDU efficiency, alpha_CE, mass-loss, and mass transfer, but not this intershell extrapolation. So the s-process yield reduction is less robust than the carbon reduction, which is largely independent of this assumption. This doesn't sink the paper, but it means the absolute s-process yields are more uncertain than the sensitivity tests suggest.\n\nThe CSLF calibration also assumes the observed carbon stars are mostly intrinsic single stars. The authors filter out some extrinsic candidates, but it's still an assumption. Given that the central C reduction is a relative comparison within the same code, this is less worrying, but it's worth noting.\n\nThe paper also overproduces high [Ce/Y] Ba stars and predicts rare massive Ba stars >10 Msun. The authors call these rare, which is fair, but it suggests the model is a bit too efficient at making s-process material.\n\nAll in all, the paper is a serious contribution to the field. It deserves peer review, and the referee should focus on the low-mass s-process treatment and ask for the data tables to be made available. I'd bring it to a reading group.","headline":"A serious binary population synthesis paper with a robust carbon yield result, but the s-process yields rest on an untested low-mass extrapolation.","tokens_in":35497,"tokens_out":2889,"would_cite":true,"duration_ms":26682,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A stellar population with a binary fraction of 0.7 ejects about 20–25% less carbon, strontium, barium, and lead than a single-star-only population, while N and O yields barely change.","keywords":["asymptotic giant branch stars","binary population synthesis","stellar yields","third dredge-up","s-process nucleosynthesis","carbon stars","barium stars","solar metallicity"],"falsifier":"If high-precision astrometry and radial velocities show that a substantial fraction (say more than 20%) of the carbon stars in the luminosity-function sample used for calibration are extrinsic binary-polluted stars rather than intrinsic AGB stars, the calibrated third dredge-up parameters and the resulting 20–25% yield reduction are biased.","tokens_in":34367,"feed_emoji":"⭐","tokens_out":10345,"duration_ms":81230,"temperature":0.7,"pith_summary":"AGB stars are major suppliers of carbon, nitrogen, and roughly half of the elements heavier than iron, but standard yield calculations treat them as single stars even though most such stars have companions. This paper uses a modified binary population synthesis code, with third dredge-up calibrated to the luminosity function of Galactic carbon stars and intershell abundances fitted to 328 isotopes, to compute population yields for binary fractions from 0 to 1 at solar metallicity. The central result is that a population with a binary fraction of 0.7 ejects roughly 20–25% less carbon and $s$-process elements such as Sr, Ba, and Pb than a single-star-only population, because binary interactions truncate or prevent the thermally pulsing AGB phase that manufactures and dredges up these elements. Nitrogen and oxygen yields change by only a few percent, with the N decrease and O increase nearly cancelling across the population.","feed_headline":"Binaries cut carbon and s-process yields by a quarter","feed_subtitle":"Stellar populations with a 0.7 binary fraction eject about a fifth less carbon and s-process material.","key_machinery":"The machinery is a synthetic binary population synthesis code with three modifications that carry the argument. First, the third dredge-up parameters are recalibrated to the observed Galactic carbon-star luminosity function, yielding a minimum core-mass reduction of $\\Delta M_{c,\\min} = -0.13\\,M_\\odot$ and a minimum dredge-up efficiency $\\lambda_{\\min}=0.45$. Second, the helium intershell abundance table is refitted to detailed AGB models covering 328 isotopes, with the heavy-element table keyed to the number of third dredge-up events rather than the number of thermal pulses, so $s$-process production (the slow neutron-capture chain that builds elements up to lead) is tied to actual dredge-up. Third, a grid of 640,000 binary systems and 1,000 single stars is weighted by a standard initial mass function and period/mass-ratio distributions to give population yields per solar mass of star-forming material. The third dredge-up is the convective mixing event that carries carbon and $s$-process products from the helium intershell to the stellar surface; without it the yields of C, Sr, Ba, and Pb are drastically reduced.","core_discovery":"The paper's central discovery is that binary companions suppress, rather than enhance, the net chemical return of low- and intermediate-mass stars. In the weighted population with a binary fraction of 0.7, the C yield falls to 0.82 of the single-star value, N to 0.95, and O rises to 1.02; Sr, Ba, and Pb fall to 0.76, 0.75, and 0.77 respectively (Tables 3 and 4). The mechanism is that only 60% of binary systems produce a TP-AGB star with at least five thermal pulses, versus 78% of single stars; common-envelope ejection and Roche-lobe overflow cut the AGB phase short before the third dredge-up can deliver carbon and $s$-process elements to the wind. There are rare counterexamples: low-mass binaries can overproduce carbon through extra thermal pulses after accretion or through helium/CO white-dwarf mergers of the R Coronae Borealis type, and some merged stars sustain hot-bottom burning longer and overproduce nitrogen. The predicted Ba star population reproduces most observed [Ce/Y] values but overproduces the fraction with [Ce/Y] > +0.2, leading the authors to conclude their models are over-efficient at making $s$-process elements in low-mass stars.","pith_inferences":["If the 20–25% reduction is real, chemical evolution models that take AGB yields from single stars alone are systematically over-injecting carbon and $s$-process elements; the binary fraction becomes a needed parameter in galactic chemical evolution calculations.","The overproduction of high-[Ba/Fe] and high-[Ce/Y] stars suggests the assumed $^{13}$C pocket size or dredge-up efficiency in low-mass models is too generous; a larger observed sample of Ba stars with measured white-dwarf masses could quantify the mismatch.","Extending the same grid to lower metallicity (e.g. Z = 0.001) would test whether the binary suppression factor changes with metallicity, since the $s$-process efficiency and mass-transfer rates both depend on Z.","The predicted massive Ba stars from intermediate-mass binaries without a common-envelope phase are a falsifiable population: a radial-velocity survey of s-process-enhanced early-type stars near the predicted orbital periods of 3000–8000 days could find or exclude them."],"forward_implications":["At a binary fraction of 0.7, the population ejects about 18% less C and about 24–25% less Sr, Ba, and Pb than a single-star-only population; at binary fraction 1.0 the reductions reach 24% for C and 33% for Ba.","Binary interactions reduce the number of TP-AGB stars by 23% and the number of hot-bottom burning systems by 24%, which propagates into lower carbon-star formation: 40% of binary primaries become C stars versus 51% of single stars.","The N yield drops only 5% and O rises 2% at a binary fraction of 0.7, so binary evolution does not strongly alter the CNO return of AGB populations even though it changes the distribution of [C/O], [N/O], and [C/N] among individual systems.","The model predicts 40% of Ba stars with [Ce/Y] > +0.2 versus 24% observed, and a rare population of about 30 massive (>10 $M_\\odot$) Ba stars in the Milky Way that has not yet been observed."],"supporting_citations":[{"why":"Supplies the synthetic AGB and binary evolution code on which all model grids are built.","marker":"Izzard et al. (2004)"},{"why":"Provides the detailed AGB models whose 328-isotope intershell abundances are fitted in the updated code.","marker":"Karakas and Lugaro (2016)"},{"why":"Provides the observed Galactic carbon-star luminosity function used to calibrate the third dredge-up parameters.","marker":"Abia et al. (2022)"},{"why":"Motivates the assumption that the low-luminosity carbon stars used in the calibration are mostly intrinsic.","marker":"Izzard and Tout (2004)"},{"why":"Supplies the initial mass function used to weight the synthetic grid into a physical stellar population.","marker":"Kroupa (2001)"},{"why":"Supplies the TP-AGB wind mass-loss prescription used to compute the ejected yields.","marker":"Vassiliadis and Wood (1993)"},{"why":"Provides the observed Ba star [Ce/Y] and [Fe/H] sample used to test the predicted Ba star abundances.","marker":"Cseh et al. (2018)"},{"why":"Provides observed Ba star masses and orbital periods used to test the predicted Ba star systems.","marker":"Jorissen et al. (2019)"}],"fun_headline_variants":["Binaries cut carbon and s-process yields by up to a quarter","Stellar companions reduce AGB carbon and s-process yields","With 70% binaries, carbon and s-process yields fall ~25%","Binary evolution lowers carbon and s-element yields in AGB stars"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The third dredge-up calibration assumes that the Galactic carbon stars used for the fit are mostly intrinsic single AGB stars rather than binary-polluted extrinsic carbon stars; if a substantial fraction are extrinsic, the calibrated dredge-up parameters and the resulting yield reductions are biased.","fun_headline_variants_meta":{"raw":{"variants":["Binaries cut carbon and s-process yields by up to a quarter","Stellar companions reduce AGB carbon and s-process yields","With 70% binaries, carbon and s-process yields fall ~25%","Binary evolution lowers carbon and s-element yields in AGB stars"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000849,"raw_usage":{"total_tokens":3774,"prompt_tokens":1109,"completion_tokens":2665,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":725,"completion_tokens_details":{"reasoning_tokens":2591}},"tokens_in":725,"tokens_out":2665,"duration_ms":18680,"temperature":1.0,"reasoning_tokens":2591,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T04:45:17.853375+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"If high-precision astrometry and radial velocities show that a substantial fraction (say more than 20%) of the carbon stars in the luminosity-function sample used for calibration are extrinsic binary-polluted stars rather than intrinsic AGB stars, the calibrated third dredge-up parameters and the resulting 20–25% yield reduction are biased.","supporting_citations":[],"review_version":1}