{"id":"5bef0972-dba7-4881-83d4-71581a29a64f","arxiv_id":"2608.11940","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Past mass gainers in massive binary systems are identifiable by their CNO surface abundance pattern, and an analytic model can reconstruct their accretion history from observed abundances.","lead":"This paper shows that stars that once swallowed material from a binary companion carry a distinctive chemical fingerprint on their surfaces: high nitrogen with ordinary oxygen. The authors built an analytic tool that reads this fingerprint to reconstruct the original binary masses and how much material was transferred, and they use it to reinterpret several well-known hot stars.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Single-star baseline omits internal gravity wave mixing, so the claimed exclusivity of the CN-eq.+dilution branch for mass gainers is not established; the paper cites the relevant process but does not model it.","rationale":"The reader's CONDITIONAL verdict is driven by the same weakest assumption: the exclusivity of the CN-eq.+dilution branch depends on the completeness of the single-star comparison set. I agree with that assessment. The paper is otherwise strong: the analytic CNO lines are derived from conservation of CNO nuclei, the mock-star test recovers the injected accretion parameters, and the MESA grid and input files are publicly available. The algebra is internally consistent, and I found no flaw in Eqs. 4-15. The issue is external: the diagnostic's power depends on the absence of single-star processes that produce the same signature. The cited ref 18 directly suggests such a process exists. The paper's own Extended Data Fig. 1 shows that different single-star rotational mixing treatments already shift single stars toward the gainer region in the HeN diagram, and only the CNO diagram is claimed to break the degeneracy; but that break was tested only against two rotational-mixing recipes. The caveat paragraph's statement that gainers remain distinguishable 'regardless of rotational mixing treatment' is not evidence against IGW mixing, which is a wave-driven mixing process in radiative zones and is not calibrated by the same physics. Therefore the central fingerprint is conditional on a missing baseline. The proposed computation is feasible and would settle the concern: if IGW-mixed single stars remain below or beside the branch, the fingerprint survives; if they populate it, the identified stars can no longer be uniquely classified as gainers, and the derived accretion histories (f_CNO, Y_CNO, M_acc_CNO, initial masses) would be ambiguous. Verdict unchanged: CONDITIONAL is appropriate.","tokens_in":27710,"tokens_out":4448,"duration_ms":43759,"concrete_test":"Add an internal gravity wave mixing prescription to the single-star baseline: recompute the Jin+2024 single-star grid (same MESA setup, nuclear network, overshooting, rotation) with an IGW diffusivity calibrated to reproduce the B-type binary/single CNO abundances in Brinkman+2025, and plot surface log(N/C) vs log(N/O) over the full main sequence and core helium burning. The exclusivity claim fails if any single-star model with initial mass 5-50 Msun, vsini<400 km/s enters the branch region (e.g., log(N/O) between -0.5 and +0.5 with log(N/C)>1) during core hydrogen or helium burning. A weaker but informative check is to fit the four observed stars with IGW-mixed single-star models and compare chi-squared to the mass-gainer interpretation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the branch near the 'CN-eq.+dilution'-line in Fig. 1 is 'exclusively populated by mass gainers and is avoided by mass donors and single stars.' This exclusivity is what turns the analytic CNO framework (Eqs. 7, 11, 15) into an identification tool for apparent single stars such as gamma Col, HD 48279, HD 93840, and zeta Oph. The evidence against single-star occupancy is the blue contours from Jin+2024 and, in Extended Data Fig. 1, Ekstrom+2012. Neither grid includes mixing by internal gravity waves, although ref 18 (Brinkman+2025) is cited as showing that IGW mixing can explain CNO surface abundances of B-type detached eclipsing binaries and single stars. IGW mixing transports CN-processed material from the radiative region just above the convective core without necessarily bringing CNO-equilibrium core material to the surface; that is precisely the chemical configuration the authors identify with the gainer branch (high N/C at moderate N/O). The paper's caveat paragraph argues that mass gainers still show higher N/C than single stars 'regardless of rotational mixing treatment,' but that statement addresses rotational mixing, not internal gravity waves. Since the analytic framework and all applications inherit the exclusivity assumption, omission of an IGW single-star baseline is the weakest link.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that the surface CNO abundances of core-hydrogen-burning massive stars can identify past mass gainers in binary systems. The authors use a large MESA binary grid to show that gainers occupy a branch near the analytic 'CN-eq. + dilution' line in the log(N/C)-log(N/O) plane, a region they claim is avoided by mass donors and single stars. They derive an analytic mixing/dilution framework (Eqs. 7, 10, 11, 15) to invert observed CNO and He/N abundances for the accreted CNO-equilibrium mass and its helium content, and apply it to gamma Columbae, HD 48279, HD 93840, zeta Ophiuchi, and SN 1987A. A mock-star test based on a 22.4+7.8 Msun binary is used as validation.","tokens_in":27872,"tokens_out":6784,"duration_ms":54777,"significance":"If the exclusivity of the gainer branch is established, the paper offers a genuinely useful diagnostic: surface CNO ratios are routinely measured, and a model-independent inversion for accreted mass and composition would constrain mass-transfer physics and identify binary products that appear single. The algebraic derivations in the Methods are internally consistent and respect CNO conservation, the analytic lines are clearly useful coordinate tools, and the public availability of the MESA grid input files is a concrete strength. However, the central fingerprint claim depends on the completeness of the single-star comparison set, and the mock-star validation is a self-consistency test rather than an independent test. These two points determine whether the applications to gamma Col, HD 48279, HD 93840, and zeta Oph are as secure as the text suggests.","major_comments":[{"comment":"The central claim that the branch near the 'CN-eq. + dilution' line is 'exclusively populated by mass gainers and is avoided by mass donors and single stars' is not established by the comparison set shown. The single-star loci in Fig. 1 and Extended Data Fig. 1 come from the Jin+2024 and Ekstrom+2012 grids, neither of which includes mixing by internal gravity waves. The manuscript itself cites Brinkman+2025 (ref. 18) for the result that IGW mixing can explain the CNO surface abundances of B-type detached eclipsing binaries and single stars, and the configuration invoked there (CN-processed material mixed upward from just above the convective core without exposing CNO-equilibrium core material) is precisely the high-N/C, moderate-N/O signature assigned to the gainer branch. The caveat paragraph stating that mass gainers show higher N/C than single stars 'regardless of rotational mixing treatment' addresses rotational mixing only, not IGW mixing. To support the fingerprint claim, the authors should add an IGW-mixing single-star baseline, or analytically bound its locus in Fig. 1 and show that it does not enter the gainer branch. Without this, the identifications of gamma Col, HD 48279, HD 93840, and zeta Oph as mass gainers inherit an unproven exclusivity assumption.","section":"Results from detailed binary evolution models (Fig. 1; Extended Data Fig. 1)"},{"comment":"The mock-star validation is a self-consistency test rather than an independent validation. The test star is drawn from the same MESA grid that motivated the framework, so recovering its input quantities demonstrates internal consistency of the inversion, not agreement with independent physics. In addition, the recovered constraints in Extended Data Table 2 are very broad: for the Mock star, M1,i = 4.8-45.2 Msun versus the true 22.4 Msun, M2,i = 3.2-7.9 Msun versus 7.8 Msun, and beta = 0.03-0.92 versus 0.05. The statement that the method 'successfully reproduces the key properties of the accreted material ... and the initial binary configuration' therefore overstates the precision demonstrated. The authors should present this as a consistency check with the reported widths, or validate the inversion against models computed with different mixing prescriptions.","section":"Methods: 'The analytic framework and its application'; Extended Data Table 2"},{"comment":"The phrase 'an analytic framework which is independent of specific evolutionary models' is too strong. The inversion uses single-star model inputs for the donor's H/He gradient mass M1,CNO(M1,i), for the envelope mass Menv, and for the evolutionary mass Mevol, all taken from specific single-star grids. The framework is independent of the binary evolution models, which is valuable, but it is not independent of evolutionary models in general. This distinction should be stated explicitly so that the model dependence of the inferred initial masses and accretion efficiencies is not underestimated.","section":"Abstract and Methods: 'The analytic model', 'Constraints on the initial primary mass', 'Evolutionary mass'"}],"minor_comments":[{"comment":"The displayed formula reads 'log Mspec/Msun = log L/Lsun - log L/Lsun', which is a tautology as printed and cannot be the intended mass-luminosity relation; the missing numerical relation or spectroscopic luminosity term should be supplied.","section":"Methods: 'Spectroscopic mass'"},{"comment":"The SN 1987A point in Fig. 1 is shown after a correction for non-solar LMC CNO ratios, but this correction is discussed only in the caption and the Supplement. A sentence in the main text stating that the plotted SN 1987A position is the LMC-corrected value would prevent readers from misinterpreting the raw abundance measurement.","section":"Fig. 1 caption and Supplementary Information Section A"},{"comment":"The caption notes that the Ekstrom+2012 models have 'slightly different initial abundances compared to ours' but does not quantify the difference; given that the limiting lines are anchored to the initial CNO ratios, the comparison would be cleaner if the figure stated the initial C, N, and O values for both grids.","section":"Extended Data Figure 1 caption"}],"recommendation":"major_revision","confidential_remarks":"The algebraic framework and the population-level trends are a solid contribution, but the paper's headline claim of an exclusive mass-gainer fingerprint needs a stronger empirical baseline. I would ask the authors to confront the IGW-mixing channel in a quantitative way before publication, and to reframe the mock-star test as a consistency check with its actual widths. The paper is otherwise well within the scope of the journal and should be revisable without changing its central methodology."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should read this one. Jin & Langer present a comprehensive MESA binary grid, a new analytic CNO-abundance framework, and applications to gamma Col, HD 48279, HD 93840, zeta Oph, and SN 1987A. The central idea: mass gainers land on a distinct branch in the log(N/C)-log(N/O) plane, near the 'CN-eq.+dilution' line, because they accrete CN-processed material without dragging CNO-equilibrium core matter to the surface. If this holds, surface abundances become a practical population tool for finding the large unresolved population of mass gainers.\n\nWhat's new and good: the corrected CNO-cycling limiting line (Eq. 7), the CN-processing dashed-line family (Eq. 11), and the CN-equilibrium max-line (Eq. 15) are correctly derived and internally consistent. The mock-star test recovers the input. The single-star and binary grids are public on Zenodo. The paper is clearly written, and the applications are plausible, including a more natural mass-gainer explanation for gamma Col than the earlier stripped-star proposal.\n\nThe soft spot is the load-bearing exclusivity claim. The paper says the branch near 'CN-eq.+dilution' is 'exclusively populated by mass gainers and is avoided by mass donors and single stars.' That conclusion rests on two rotating single-star model grids, neither of which includes mixing by internal gravity waves. The paper cites Brinkman+2025 as a known mechanism that can produce high N/C at moderate N/O in single stars, but it does not model that mechanism as a baseline. IGW mixing can transport CN-processed material from the radiative region above the convective core without bringing CNO-equilibrium core material to the surface; that is exactly the configuration the authors identify with the gainer branch. The caveat paragraph addresses rotational mixing, not IGW. So the exclusivity claim is not yet established. This weakens but does not break the paper: the analytic framework remains useful, and the numerical models do show gainers occupying that region, but 'exclusively' is too strong given the current comparison set.\n\nOther soft spots are minor. The mock star is drawn from the same MESA grid that motivated the framework, so its recovery is a self-consistency check rather than an independent validation. The 'model-independent' framework still leans on single-star tracks for evolutionary and envelope masses. Error propagation into derived parameters is incomplete.\n\nWho should read it: anyone working on massive binary evolution, OB star surface abundances, or SN 1987A's progenitor. It deserves a serious referee. I'd send it to review with a strong request to either include an IGW-inclusive single-star baseline or soften the exclusivity language. If that's done, this will be a solid contribution.","headline":"Solid, genuinely new analytic framework for identifying binary mass gainers via CNO abundances, but the 'exclusive' branch claim is not yet established because the single-star baseline omits internal gravity wave mixing.","tokens_in":28536,"tokens_out":6288,"would_cite":true,"duration_ms":53104,"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":"Chemical fingerprint exposes stars that once gained mass.","keywords":["massive stars","binary mass transfer","mass gainers","CNO cycle","surface abundances","stellar mergers","SN 1987A","thermohaline mixing"],"falsifier":"A decisive check would be to find a main-sequence massive star with a well-established single-star history (no companion, no past accretion, constant radial velocity) whose measured N/C and N/O place it on or above the CN-eq.+dilution branch; alternatively, a single-star model grid that includes internal gravity-wave mixing and populates the branch would falsify the claimed exclusivity.","tokens_in":27387,"feed_emoji":"⭐","tokens_out":8261,"duration_ms":79339,"temperature":0.7,"pith_summary":"The paper claims that the carbon, nitrogen, and oxygen (CNO) surface abundances of core-hydrogen-burning massive stars carry a fingerprint of past binary mass transfer. In the diagnostic N/C versus N/O diagram, the branch near the CN-equilibrium-plus-dilution line is populated exclusively by mass gainers (stars that accreted material from a binary companion), while mass donors and single stars stay away from it. If this is right, many apparently single OB stars with elevated nitrogen are not merely rapidly rotating single stars but former accretors whose chemical surface patterns record an earlier binary interaction. The authors provide an analytic framework, independent of the details of their evolutionary models, that converts an observed position in this diagram into the amount and composition of accreted material, and they apply it to several well-studied stars and to SN 1987A.","feed_headline":"Chemical fingerprint exposes stars that once gained mass","feed_subtitle":"Nitrogen and carbon ratios separate mass gainers from single stars and reveal how much material was swallowed.","key_machinery":"The diagnostic CNO abundance plane, with $\\log(N/O)$ on one axis and $\\log(N/C)$ on the other, is framed by two analytic dilution lines: mixing pristine envelope matter with CNO-equilibrium matter (the 'CNO-eq.+dilution' line, Eq. 7) and mixing with CN-equilibrium matter (the 'CN-eq.+dilution' line, Eq. 10). Superposed are dashed tracks, Eq. 11, describing CN-cycling of an already diluted mixture, which end at the full-CN-processing line of Eq. 15. These lines do the work of the argument: an observed star's position fixes the dilution factor $f_{\\rm CNO}$ by projection onto the CNO-eq.+dilution line, and the helium-nitrogen diagram fixes $Y_{\\rm CNO}$, giving a model-independent reconstruction of the accretion history.","core_discovery":"On the paper's own terms, the central discovery is that stars which have accreted matter in a binary can be recognized long after the interaction by their surface CNO ratios alone. In a large grid of binary evolution models, mass gainers occupy a characteristic branch of the log(N/C) versus log(N/O) diagram that mass donors and single stars do not enter: matter accreted from the donor's hydrogen/helium-gradient zone is in CNO equilibrium, and later slow mixing plus CN-processing in the gainer's envelope raises N/C at roughly constant N/O. The paper argues that this branch is a unique chemical fingerprint, and that the amount of accreted CNO-equilibrium material (the dilution factor $f_{\\rm CNO}$) and the helium content of that material ($Y_{\\rm CNO}$) can be read off from an observed star's position using analytic mixing lines.","pith_inferences":["If the fingerprint survives comparison with a wider set of single-star mixing processes, the method becomes a population tool: the distribution of $f_{\\rm CNO}$ values among field stars would map the mass-transfer efficiency distribution across initial binary parameter space.","The same dilution-plus-CN-cycling formalism should extend to isotope ratios such as $^{13}{\\rm C}/^{12}{\\rm C}$ or $^{15}{\\rm N}/^{14}{\\rm N}$, which could break remaining degeneracies between accretion and rotational mixing in stars where element ratios alone are ambiguous.","A testable prediction of the accretion picture is that stars on the gainer branch should show signatures of accretion-induced spin-up or mixing, such as unusually rapid rotation or surface helium enrichment, more often than single stars of similar mass and age; a targeted survey of the branch could check this.","The exclusivity of the branch could be probed directly by computing single-star models with internal gravity-wave mixing as a baseline; if such models enter the branch, the fingerprint would need to be redefined."],"forward_implications":["Apparent single OB stars with high N/C at moderate N/O can be classified as former mass gainers, turning existing abundance surveys into a census of past binary accretion.","For each identified gainer, the analytic framework yields the mass of accreted CNO-equilibrium material, the average helium content of that material, and allowed ranges for the initial masses of both binary components and the mass-transfer efficiency.","The method separates stable mass transfer from mergers: merger products sit near the CNO-eq.+dilution line, whereas gainers from stable transfer can rise along the CN-cycling tracks, so SN 1987A's ring abundances point to a diluted CNO-equilibrium mixture consistent with a post-main-sequence merger.","For gamma Columbae, the paper concludes it is a mass gainer whose companion likely exploded as a stripped-envelope supernova, replacing the earlier interpretation of the star as a recently stripped object.","The constraints on initial mass ratio and accretion efficiency derived this way provide empirical benchmarks for future binary evolution models, independent of the models' assumed accretion physics."],"supporting_citations":[{"why":"Supplies the comprehensive binary evolution grid of mass gainer and donor surface abundances on which the two-branch pattern is identified.","marker":"(34)"},{"why":"Provides the comparison grid of rotating single-star models that the gainer branch is claimed to avoid.","marker":"(19)"},{"why":"Establishes the analytic dilution-line formalism for CNO surface abundances that the paper extends and corrects.","marker":"(17)"},{"why":"Prior analytic limiting lines for CN/CNO cycling; the paper shows its CNO-cycling line deviates from theirs at high N/O.","marker":"(27)"},{"why":"Source of gamma Columbae's surface abundances and the stripped-star interpretation that the paper argues against.","marker":"(21)"},{"why":"Circumstellar ring abundances of SN 1987A used to test the merger application.","marker":"(24)"},{"why":"Provides the critical mass-ratio criterion for stable mass transfer used to constrain the initial binary configurations.","marker":"(40)"},{"why":"Alternative single-star mixing process, internal gravity waves, that could in principle populate the gainer branch and threaten the fingerprint's uniqueness.","marker":"(18)"}],"fun_headline_variants":["Surface CNO ratios unmask stellar mass gainers","Chemical fingerprint reveals stars that gained mass in binaries","Gamma Col is a mass gainer, not a stripped star","Stellar chemistry exposes hidden binary mass transfer"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The claim that the branch near the CN-equilibrium-plus-dilution line belongs only to mass gainers depends on the comparison single-star models being complete: if some ordinary single-star process, such as internal gravity-wave mixing, can push a non-accreting star into the same region, the fingerprint would no longer be unique.","fun_headline_variants_meta":{"raw":{"variants":["Surface CNO ratios unmask stellar mass gainers","Chemical fingerprint reveals stars that gained mass in binaries","Gamma Col is a mass gainer, not a stripped star","Stellar chemistry exposes hidden binary mass transfer"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000339,"raw_usage":{"total_tokens":1862,"prompt_tokens":924,"completion_tokens":938,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":540,"completion_tokens_details":{"reasoning_tokens":876}},"tokens_in":540,"tokens_out":938,"duration_ms":9913,"temperature":1.0,"reasoning_tokens":876,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:22:45.750856+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive check would be to find a main-sequence massive star with a well-established single-star history (no companion, no past accretion, constant radial velocity) whose measured N/C and N/O place it on or above the CN-eq.+dilution branch; alternatively, a single-star model grid that includes internal gravity-wave mixing and populates the branch would falsify the claimed exclusivity.","supporting_citations":[],"review_version":1}