{"id":"935b98fd-f3b3-4293-94db-815263ed482c","arxiv_id":"2411.13489","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Actively accreting stellar binaries with well-measured orbits lie on the equilibrium eccentricity curve predicted by circumbinary disk theory.","lead":"This paper checks whether young double stars still surrounded by a shared gas disk have the stretched, oval orbits that accretion theory predicts. It reports that the aligned systems in a published catalog sit on the predicted equilibrium curve, while tilted systems sort into the predicted aligned or polar states.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Claimed fit to the attractor curve is never quantified; without eccentricity error bars, the equal-q pair TWA 3A and HD 200775 alone may refute the fit, so the central result is unestablished.","rationale":"The paper is a short observational check of a real theoretical prediction, and the classification of polar versus aligned systems using the Martin and Lubow critical-inclination curve is a useful independent consistency test. I read the central claim as narrow: for the subset of actively accreting binaries with well-measured parameters, planar alignment, and negligible tidal circularization, the observed e-q pairs should lie near the simulated attractor curve. That claim would be genuinely interesting because field binaries show no such relation. The weakest point is not the existence of the theory curve but the evidence that the seven points actually agree with it. The text asserts the agreement and even calls it precise, but provides no error bars, no fit statistic, and no discussion of scatter. The TWA 3A and HD 200775 pair exposes the problem: two systems at nearly the same q differ by 0.33 in e, making it impossible for a single-valued increasing e_attr(q) to be close to both unless the quoted or unquoted errors are very large. Since the paper does not supply those errors, the reader cannot distinguish 'consistent with theory' from 'scattered within large errors.' The missing uncertainties are therefore the load-bearing issue; the accretion-history inference flagged by the Reader is downstream of it. If the test I propose confirms small residuals, the verdict should be ACCEPT with the caveat that a statistical presentation is still needed; if it shows large outliers, the claim should be REJECTED as an overinterpretation. I do not see fraud or circular reasoning; the issue is standard quantitative rigor. The Reader's CONDITIONAL verdict is appropriate, and I would not change it.","tokens_in":7185,"tokens_out":7587,"duration_ms":84523,"concrete_test":"Reconstruct Figure 2 with the full Czekala et al. (2019) uncertainties for e and q for the seven planar-aligned systems. Compute residuals (e_obs - e_attr(q))/sigma_e and a chi-square or RMS statistic, including the interpolation uncertainty in e_attr(q) from Siwek et al.'s table. Specifically check TWA 3A and HD 200775: if both lie within 1 sigma of the same e_attr(q) curve, the visual fit is confirmed; if one or both fall more than 2 sigma away, the claim fails. Also repeat the test against a constant-e null model to show that the promised q-dependent curve is actually preferred.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is the claimed quantitative agreement in Figure 2 between the planar-aligned systems and the Siwek et al. (2023) attractor curve e_attr(q). This agreement is asserted visually: the paper reports no observational uncertainties for e or q, no residual statistic, and no uncertainty for the interpolated theoretical curve. The omission is not cosmetic. TWA 3A (e=0.628, q=0.84) and HD 200775 (e=0.30, q=0.82) have nearly identical mass ratios yet differ in eccentricity by 0.33; since e_attr(q) is a single-valued curve that increases with q in the Siwek et al. results, no such curve can pass within a small tolerance of both systems. Either their eccentricity error bars are large enough (about 0.15 or more) that the match carries little information, or at least one of these 'aligned' systems is actually inconsistent with the predicted curve. The paper also refers to 'eight aligned systems' while Table 1 lists only seven planar-aligned systems, excluding the two tidally circularized e=0 systems, so the internal counting reinforces that the fit was not independently checked. Consequently, the central claim that observations 'follow' the equilibrium curve is not yet supported by the data as presented; it requires a quantitative goodness-of-fit test.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper assembles the well-characterized binary–disk systems from Czekala et al. (2019) and classifies them by mutual binary–disk inclination into tidally circularized, planar aligned, hierarchical trinary, and polar aligned categories, using the Martin–Lubow critical-inclination criterion. For the planar aligned systems, the authors compare the observed eccentricity and mass ratio against the equilibrium (attractor) eccentricity curve of Siwek et al. (2023) and claim that the systems lie on that curve, suggesting that they have accreted enough gas for circumbinary-disk torques to set their eccentricity, and that eccentricity evolution resumes only after disk dispersal.","tokens_in":7407,"tokens_out":2705,"duration_ms":32399,"significance":"If the central comparison is established, this would be a valuable observational test of circumbinary-disk eccentricity theory that was developed largely for supermassive black hole binaries, and it would add an interesting post-disk evolutionary constraint for young stellar binaries. The paper has genuine strengths: the theoretical curve is taken from published simulations rather than fitted to the data, the aligned/polar classification uses an independent published criterion, and the sample is small but carefully selected from a homogeneous catalog. The main significance, however, hinges entirely on the quantitative agreement in Figure 2, and that agreement is currently asserted visually rather than demonstrated.","major_comments":[{"comment":"The central claim that the aligned systems lie on the Siwek et al. (2023) attractor curve is not quantified. The paper reports no observational uncertainties for eccentricity or mass ratio, no residual statistic, and no uncertainty on the interpolated theoretical curve. This is not a cosmetic omission: TWA 3A (e=0.628, q=0.84) and HD 200775 (e=0.30, q=0.82) have nearly identical mass ratios but differ in eccentricity by 0.33. Since e_attr(q) is a single-valued, monotonically increasing curve over this q range in Siwek et al., no curve of that shape can pass within a small tolerance of both systems. Either the error bars on e are large enough (about 0.15 or more) that the match carries little information, or at least one of these two systems is inconsistent with the predicted curve as drawn. The authors should provide eccentricity and mass-ratio errors, a measure of scatter or a goodness-of-fit statistic, and the uncertainty on the interpolated attractor curve, so that the claimed agreement can actually be evaluated.","section":"§4, Figure 2"},{"comment":"There is an internal counting inconsistency that suggests the fit was not independently checked: the text in §4 states that 'the eight aligned systems lie approximately upon the attractor eccentricity curve,' but Table 1 lists only seven planar aligned systems (HD 131511, αCrB, TWA 3A, HD 200775, Ak Sco, DQ Tau, UZ Tau E). The two tidally circularized systems (V4046 Sgr, CoRoT 2239) have e=0 and are explicitly placed in a separate category, so they cannot be the eighth aligned system. The paper should correct the count and clarify exactly which systems are used in the Figure 2 fit.","section":"Table 1 and §4"},{"comment":"The interpretation that the aligned systems 'have accreted enough gas to significantly affect their orbit' is inferred from the match itself, because the sample selection does not include an accretion history or age criterion that would independently identify systems that have reached the equilibrium. As written, the agreement could also arise if these young binaries were simply born near the curve or if the catalog selection favors such systems. To make the evolutionary claim load-bearing, the authors should either add an independent check (for example, comparing accretion indicators or estimated accreted mass against the predicted equilibration timescale) or explicitly reframe the conclusion as a statement about the current e–q distribution only.","section":"§2 and §4"}],"minor_comments":[{"comment":"The wording oscillates between 'consistent with' (abstract), 'lie approximately upon' (§4), and 'precisely obeys' (§5); the final wording is too strong given the lack of quantified residuals and should be harmonized with the actual statistical support.","section":"§4 and §5"},{"comment":"Equation (1) is typeset in a garbled way in the manuscript text, and the displayed formula is hard to read; it should be reset so that the argument of the arccosine and the eccentricity dependence are unambiguous.","section":"§3, Equation (1)"},{"comment":"The sentence describing the interpolation of points where ė_b=0 and the interpolation for e_b=0.7 should specify whether the interpolation is linear in the (q, e) plane and whether the resulting curve has any systematic uncertainty from the coarse tabulation in Siwek et al. (2023).","section":"§3"},{"comment":"The procedure of setting i=0 and using the upper limit as the uncertainty for systems with i<10° is reasonable, but the same treatment is not applied to the eccentricity and mass-ratio columns, which are listed without any uncertainties; providing the original error bars from Czekala et al. (2019) in a supplementary table or in the caption would greatly improve the transparency of the comparison.","section":"§2, Table 1"}],"recommendation":"major_revision","confidential_remarks":"The central scientific idea is attractive and the paper is short and readable, but the main quantitative claim is currently unsupported. The most serious issue is the lack of uncertainties and a goodness-of-fit test for the Figure 2 comparison; the TWA 3A / HD 200775 pair alone shows that the claimed agreement is not self-evident. I would encourage the editor to request a revision that adds error bars, a residual analysis, and a corrected count of aligned systems. On the citation front, the author overlap with the theory papers (Duffell et al. 2020; D'Orazio & Duffell 2021) is transparent, and the theory curve is not fitted to the data here, so I do not see a circularity problem, but the referee should be aware of the self-citation when judging novelty."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick read of Murray & Duffell. The new bit: nobody has put the Czekala et al. stellar CBD systems on the Siwek et al. attractor curve before, and the paper gives a clean statement of what an observational confirmation would look like. It also does something useful before the main event: the polar/aligned split using the Martin-Lubow critical inclination works for the non-triple systems, which is a nice independent check. Credit where due: the paper is short, readable, and takes the theory seriously without trying to build a new framework.\n\nThe soft spot is the load-bearing comparison. The central claim—that the aligned systems lie on e_attr(q)—is asserted visually from seven points. Table 1 lists no uncertainties for e or q, and there is no residual statistic or goodness-of-fit test. That is not cosmetic. TWA 3A (e=0.628, q=0.84) and HD 200775 (e=0.30, q=0.82) have nearly identical mass ratios and eccentricities that differ by 0.33. If the attractor curve is single-valued and increasing in q, no such curve can pass within a small tolerance of both. Either the eccentricity uncertainties are around 0.15 or larger, in which case the fit carries little information, or one of these \"aligned\" systems is actually off the curve. The authors do not tell us which. The text also says \"eight aligned systems\" while Table 1 lists seven planar-aligned (excluding the two tidally circularized e=0 systems), which suggests the fit was never carefully checked. The conclusion that these systems \"precisely obey\" the predictions goes beyond what the data show.\n\nThere is also the selection-function issue: no model for how the Czekala catalog was assembled, and the inference that these stars \"have accreted enough gas to significantly affect their orbit\" is drawn from the very match being claimed. That is a history assumption, not a result.\n\nWho is this for? People working on circumbinary disk evolution and stellar binary formation, and anyone using Siwek-style curves for SMBHB population predictions. The paper is a legitimate first attempt and deserves referee time, but the fixable flaws are exactly the ones that make the difference between \"appears to fit\" and \"follows.\" I would send it to review and ask for a quantitative fit, error bars, and a selection discussion. As published, I'd treat the central claim as unconfirmed.","headline":"First stellar test of the Siwek attractor curve is a good idea, but the claimed fit rests on seven points with no error bars and a pair of systems that already look inconsistent.","tokens_in":7952,"tokens_out":2987,"would_cite":true,"duration_ms":32914,"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":"Actively accreting stellar binaries that are aligned with their circumbinary disks lie on the predicted equilibrium eccentricity curve.","keywords":["circumbinary disks","binary eccentricity","accreting binaries","equilibrium eccentricity","stellar accretion disks","binary-disk interaction","mass ratio","protostellar binaries"],"falsifier":"A larger sample of young circumbinary-disk binaries, selected without eccentricity or inclination priors, would refute the claim if the planar-aligned, non-tidally-circularized systems scatter far from the $\\dot{e}_b=0$ curve; specifically, if a majority of a dozen new systems fall more than about 0.1 in eccentricity from the Siwek et al. attractor, the equilibrium interpretation would be unsupported.","tokens_in":6945,"feed_emoji":"🌟","tokens_out":18082,"duration_ms":160050,"temperature":0.7,"pith_summary":"Accreting binaries that still have circumbinary disks occupy eccentricity–mass-ratio positions that match the equilibrium curve predicted by disk–binary hydrodynamics. The authors select from a catalog of binary–disk systems those with well-measured orbital parameters, no third body, mutual inclination consistent with zero, and separations too wide for tidal circularization; this aligned subset lies on the predicted curve. These are the first observational points tied to the predicted attractor eccentricity, rather than to the flat eccentricity distribution of the general Milky Way binary population. The paper concludes that eccentricity in such systems is set by disk accretion, and that some later mechanism must act after the disk dissipates to erase that preference.","feed_headline":"Accreting binaries match the predicted disk-torque eccentricity curve","feed_subtitle":"Aligned young binary stars with circumbinary disks sit on the equilibrium eccentricity curve predicted by theory.","key_machinery":"The load-bearing object is the equilibrium eccentricity curve: the locus of points in the $(q,e_b)$ plane at which the simulated eccentricity derivative vanishes, $\\dot{e}_b(q,e_b)=0$, interpolated from the fixed-binary, two-dimensional hydrodynamic grid of Siwek et al. (2023). Combined with the mass-ratio evolution formula $\\dot{q}(q)$ from Duffell et al. (2020), this curve yields a vector map showing binaries migrating onto the curve and then along it toward higher mass ratio. A second criterion, the Martin & Lubow (2019) critical inclination $i_{\\min}=\\arccos\\sqrt{5e_b^2/[4(1+e_b^2)]}$, separates systems that should evolve to planar versus polar alignment, and the paper uses it to define the aligned subsample to which the eccentricity curve applies.","core_discovery":"The central discovery is that the observed eccentricities of planar-aligned accreting stellar binaries trace the zero-derivative curve $\\dot{e}_b(q,e_b)=0$ computed by Siwek et al. (2023), meaning the binaries have evolved to the eccentricity at which disk torques stop changing it. The aligned systems also sit below the Martin & Lubow (2019) critical inclination curve, as they must if they are to remain aligned, while the two polar-aligned systems sit above the curve and the hierarchical triples scatter outside the simple binary–disk predictions. The authors read the match as evidence that the disks have been present long enough to drive the binaries onto the attractor; the same plot explains why the general field population, whose disks have long since vanished, shows no preferred eccentricity.","pith_inferences":["Because the aligned subsample contains only a handful of systems drawn from a catalog not designed for this test, part of the agreement could reflect selection effects; a blind, volume-limited survey of accreting binaries would sharpen the comparison.","The theoretical curve comes from two-dimensional, fixed-binary simulations, while real protoplanetary disks are three-dimensional and the binary's orbit responds to the disk; the level of agreement may be partly forgiving of those simplifications.","If the equilibrium is genuine, the same eccentricity–mass-ratio attractor should appear in other circumbinary-disk populations, such as post-asymptotic-giant-branch binaries, where a bimodal eccentricity distribution has already been predicted."],"forward_implications":["A real match would show that circumbinary disk torques reshape stellar binary orbits on observable timescales, providing a local proxy for the same physics predicted for supermassive black hole binaries.","Because the merged Milky Way binary population has a flat eccentricity distribution, some mechanism must act after the disk dissipates to erase the preferred eccentricity–mass-ratio pattern.","The position of a system like Ak Sco, with $q\\approx 1$, $e_b\\approx 0.47$, and an age of 18 Myr, is naturally read as a 'twin' binary that has accreted substantially, supporting the predicted mass-ratio evolution.","A proposed test using astrometric surveys of young binaries with no third body should recover the same attractor eccentricity if the equilibrium interpretation is correct."],"supporting_citations":[{"why":"It supplies the catalog of binary–disk systems with measured eccentricity, mass ratio, and mutual inclination, from which the aligned subsample is drawn.","marker":"I. Czekala et al. 2019"},{"why":"It predicts the equilibrium eccentricity curve $\\dot{e}_b=0$ as a function of mass ratio that the observed aligned systems are compared against.","marker":"M. Siwek et al. 2023"},{"why":"It provides earlier simulation results finding an attractor eccentricity near 0.4–0.5 for equal-mass binaries, establishing the equilibrium picture.","marker":"D. J. D'Orazio & P. C. Duffell 2021"},{"why":"It finds attractor eccentricities and a bimodal eccentricity preference for post-AGB binaries, supporting disk-driven eccentricity evolution.","marker":"J. Zrake et al. 2021"},{"why":"It provides the mass-ratio evolution formula used with the eccentricity curve to build the $q$–$e_b$ vector map.","marker":"P. C. Duffell et al. 2020"},{"why":"It gives the critical inclination criterion used to classify systems as aligned versus polar.","marker":"R. G. Martin & S. H. Lubow 2019"},{"why":"It shows the flat eccentricity distribution in the general binary population, the demographic contrast that motivates the post-disk evolution explanation.","marker":"A. M. Price-Whelan et al. 2020"},{"why":"It checks consistency between independent binary–disk simulation studies, supporting the reliability of the eccentricity predictions used here.","marker":"R. Valli et al. 2024"}],"fun_headline_variants":["Binary eccentricities follow disk-torque equilibrium curve","Accreting binaries trace predicted disk-torque curve","Disk-torque equilibrium curve matches observed binary orbits","Aligned binaries sit on predicted eccentricity attractor","Observation confirms disk-torque equilibrium for binaries"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The aligned systems have already accreted enough gas for disk torques to have driven their eccentricities to the equilibrium curve, rather than being born near it by chance or selected into the catalog because they lie near it.","fun_headline_variants_meta":{"raw":{"variants":["Binary eccentricities follow disk-torque equilibrium curve","Accreting binaries trace predicted disk-torque curve","Disk-torque equilibrium curve matches observed binary orbits","Aligned binaries sit on predicted eccentricity attractor","Observation confirms disk-torque equilibrium for binaries"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000453,"raw_usage":{"total_tokens":2217,"prompt_tokens":825,"completion_tokens":1392,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":441,"completion_tokens_details":{"reasoning_tokens":1316}},"tokens_in":441,"tokens_out":1392,"duration_ms":10909,"temperature":1.0,"reasoning_tokens":1316,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T16:21:12.690221+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A larger sample of young circumbinary-disk binaries, selected without eccentricity or inclination priors, would refute the claim if the planar-aligned, non-tidally-circularized systems scatter far from the $\\dot{e}_b=0$ curve; specifically, if a majority of a dozen new systems fall more than about 0.1 in eccentricity from the Siwek et al. attractor, the equilibrium interpretation would be unsupported.","supporting_citations":[{"cited_title":"M., et al","cited_arxiv_id":null,"evidence_quote":"It supplies the catalog of binary–disk systems with measured eccentricity, mass ratio, and mutual inclination, from which the aligned subsample is drawn."},{"cited_title":"2021, ApJL, 909, L13 4 The Astrophysical Journal, 982:113 (4pp), 2025 April 1 Murray & Duffell","cited_arxiv_id":null,"evidence_quote":"It finds attractor eccentricities and a bimodal eccentricity preference for post-AGB binaries, supporting disk-driven eccentricity evolution."},{"cited_title":"C., D ’Orazio, D., Derdzinski, A., et al","cited_arxiv_id":null,"evidence_quote":"It provides the mass-ratio evolution formula used with the eccentricity curve to build the $q$–$e_b$ vector map."},{"cited_title":"G., & Lubow, S","cited_arxiv_id":null,"evidence_quote":"It gives the critical inclination criterion used to classify systems as aligned versus polar."},{"cited_title":"M., Hogg, D","cited_arxiv_id":null,"evidence_quote":"It shows the flat eccentricity distribution in the general binary population, the demographic contrast that motivates the post-disk evolution explanation."},{"cited_title":"2024, A&A, 688, A128","cited_arxiv_id":null,"evidence_quote":"It checks consistency between independent binary–disk simulation studies, supporting the reliability of the eccentricity predictions used here."}],"review_version":1}