{"id":"36960a53-1ed2-476c-aea2-ea366fad1ada","arxiv_id":"2507.13515","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"From 77 accreting binaries, a fitted power-law family for angular momentum loss is proposed to replace the failed magnetic braking model.","lead":"This paper measures long-term orbital period changes for 25 X-ray binaries and combines them with 52 cataclysmic variables to fit new empirical laws for angular momentum loss. If these laws hold, the standard magnetic braking model of binary evolution is wrong by orders of magnitude and modelers get a new empirical replacement.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Universal-law Mdot exponent may be an artifact of regressing the subtraction residual on the same noisy Mdot used in Eq. 10; synthetic-error test can settle it.","rationale":"Fair-reading summary: the O−C compilation and the demonstration that MBM predictions fail by orders of magnitude are valuable and seem largely independent of the universal-law fit. The concern targets only the step from Pdot_AML residuals to the empirical power-law family. The reader's weakest assumption is the same one: the accretion rates and wind efficiencies used to construct Pdot_AML must be accurate enough and independent of the residual, because they are also the input regressors. I agree with that identification. The mechanism is specific: Eq. 10 makes Pdot_mt ≈ 3P(1−q) Mdot/M_comp. Using a noisy estimate of this term to build the dependent variable and then regressing the result on the same estimate injects a correlated error. Standard regression treats the independent variables as fixed; here they are not, and the dependent variable was built from them. This is not resolved by the paper's adopted 0.25-dex systematic error, which only inflates the scatter and preserves any central-value bias. The synthetic test I propose is cheap and decisive. If it shows the bias is negligible, the universal law's exponents are credible and the conditional acceptance can stand. If it reproduces δ≈0.43 with δ_true=0, the central claim is an artifact and the paper should be revised to re-fit after decoupling the subtraction, for example by using independent Mdot estimates or Bayesian joint modeling, before any 'universal' law is claimed. Since the test has not been run, I keep the reader's CONDITIONAL verdict rather than moving to REJECT; the condition should explicitly include this error-correlation test.","tokens_in":61161,"tokens_out":7165,"duration_ms":81212,"concrete_test":"Run the manuscript's pipeline on synthetic data drawn from Table 6. Take the tabulated P, M_prim, M_comp, and true Mdot as inputs; set the true AML law to have δ_true = 0 (e.g., Pdot_AML = −C P^1.29 M_prim^2.75 M_comp^−1.0) and generate Pdot_obs = Pdot_GR + Pdot_mt(true Mdot, true masses) + Pdot_AML + Gaussian noise at the quoted Pdot errors. Then multiply each Mdot and mass by independent log-normal factors with the factor-of-ten widths stated in Section 5, recompute Pdot_mt from the noisy values, form Pdot_AML_est via Eq. 2, and fit Eq. 14 with the same weighting and the same treatment of non-parabolic O−C systems as Pdot=0. If the recovered δ exceeds 0.2 (or approaches the reported 0.43), the fitted Mdot exponent is an artifact of the subtraction/regression correlation. Repeat with true Pdot_AML=0 to check whether the fitted law is entirely spurious.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—Eqs. 15–16 / Table 8 as a 'universal' AML law—requires that the fitted exponents describe a real extra angular-momentum-loss process. The weakest link is that the dependent variable is manufactured from the same quantities that later serve as regressors. Eq. 2 defines Pdot_AML = Pdot − Pdot_GR − Pdot_mt, with Pdot_mt linear in Mdot for RLOF (Eq. 10) and depending on q = M_comp/M_prim; the same Mdot and mass values are then inserted into Eq. 14, log(−Pdot_AML) = log C + α log P + β log M_prim + γ log M_comp + δ log Mdot. Section 5 concedes Mdot uncertainties of order a factor of ten and poorly known wind efficiencies ϵ. If Mdot is overestimated, the subtracted Pdot_mt is too large and Pdot_AML moves negative; if underestimated, it moves positive. Regressing this error-contaminated residual on log Mdot creates a built-in positive correlation, so a nonzero δ (and shifted β, γ through q) can appear even when the true AML has no dependence on Mdot. The fit's reduced chi-square near unity after adding an adopted 0.25-dex systematic error does not protect against this, because the bias is in the central value, not the scatter. Thus the main new result is not established as a measurement of AML unless this correlation is shown to be negligible.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper assembles eclipse and minimum-light timings for 25 X-ray binaries, measures or collects their secular orbital period derivatives, and combines them with 52 cataclysmic variables from the author's prior work. After subtracting gravitational-radiation and mass-transfer contributions (Eqs. 2, 9-11), it forms 84 Pdot_AML values for 77 systems, uses them to test magnetic braking and related AML prescriptions, and finds large discrepancies with the MBM. The paper then fits power laws in P, Mprim, Mcomp, and Mdot to Pdot_AML (Eq. 14), reports three 'universal' empirical AML laws (Eqs. 15-16 and Table 8), and claims these describe the actual evolution of all 77 XRBs and CVs.","tokens_in":61479,"tokens_out":8309,"duration_ms":88602,"significance":"The O-C measurements are a substantial and welcome contribution, especially the long-baseline timing for Sco X-1, Her X-1, and V4641 Sgr, the TESS-based timings, and the careful discussion of systematic jitter in several systems. If the AML-law result survives scrutiny, it would provide an empirically based replacement for the MBM. However, the headline universal law is a least-squares fit to the same systems it claims to describe, and the construction of Pdot_AML from the same Mdot that later serves as a regressor creates a serious correlated-error concern. The paper is therefore significant as a measurement paper and as a falsification test of MBM, but the universal-law claim is not yet established.","major_comments":[{"comment":"The headline Mdot exponent (δ=0.43) may be an artifact of regressing a constructed residual on the same noisy Mdot used to construct it. Pdot_AML is defined by subtracting a mass-transfer term that is linear in Mdot (Eq. 10), and log Mdot is then one of the regressors in Eq. (14). With Mdot uncertainties of 'typically like a factor of ten' as stated in Section 5, an overestimate of Mdot makes the subtracted Pdot_mt too large and drives Pdot_AML more negative, while an underestimate drives it less negative; this built-in anti-correlation can produce a positive δ and can also bias β and γ through q=Mcomp/Mprim even when the true AML has no Mdot dependence. Adding a 0.25-dex systematic error changes the scatter but not the central-value bias, so a reduced chi-square near unity does not validate the law. Please add a synthetic-error test (simulate data with a true Pdot_AML independent of Mdot, add factor-ten Mdot errors, and show that the fitting procedure recovers δ=0) or an errors-in-variables / orthogonal regression treatment.","section":"Sec. 8.2 / Eq. (16), with Sec. 5, Eqs. (2), (10), (14)"},{"comment":"The claim of a 'universal' law is not supported by a fit to the same 77 systems from which the law was derived. No out-of-sample or cross-validated test is presented, so the statement that Eqs. (15)-(16) are 'the best representations of the actual evolution for all 77' is a goodness-of-fit restatement, not a predictive test. Additionally, the abstract promises a third law for P>1.0 day, but Table 8 contains only the below-gap fit, the 0.13-1.0 day fit, the HMXB/IMXB fit, and an all-systems fit; long-period CVs such as U Sco, V394 CrA, and T CrB appear only in the 'All' row. The paper either needs to present the dedicated P>1.0-day fit or revise the abstract's claim of three laws.","section":"Table 8 / Sec. 8.2"},{"comment":"The statistical treatment of the fit is under-specified. The text states that a systematic error of 0.25 dex is adopted for all binary groups to bring reduced chi-square near unity, but for the below-gap fit the reduced chi-square is 0.4 after this addition, indicating that the error budget is overestimated rather than calibrated. Table 6 reports the last-column acceptable ranges for Pdot_AML as asymmetric, yet the fitting section does not say how these asymmetric errors are converted into sigma, nor how the factor-ten Mdot uncertainties and the poorly known wind-capture efficiencies epsilon are propagated into the parameter errors quoted in Table 8. Please replace the ad hoc systematic with a transparent likelihood or Monte Carlo propagation, and report the sensitivity of α, β, γ, δ to the assumed Mdot and epsilon errors.","section":"Sec. 8.1 / Table 6"}],"minor_comments":[{"comment":"The units of Pdot_AML, P, masses, and Mdot are only given in the text; the abstract's 'in appropriate units' is too terse for a headline equation.","section":"Abstract / Eq. (16)"},{"comment":"The 'k' and 'kk' shorthand in the Pdot columns is nonstandard and should be defined in the table caption or replaced by explicit powers of ten.","section":"Table 6"},{"comment":"The reduction of Eq. (11) to Eq. (10) for epsilon=1 is not shown; a short check would help readers verify that the wind and RLOF cases are consistent.","section":"Sec. 5, Eq. (11)"},{"comment":"The term 'jerks' is introduced for fast O-C kinks; since this is not a standard term, define it at first use.","section":"Sec. 8.1"},{"comment":"Chi-square values are quoted without the number of degrees of freedom in the table; the text gives some reduced values, but the table should state the dof explicitly.","section":"Table 8"},{"comment":"The comparison of 50 measures for 44 systems should state explicitly whether the multiple inter-eruption intervals are treated as independent in the statistical tests.","section":"Sec. 7.2"}],"recommendation":"major_revision","confidential_remarks":"The strength of this paper is the hard-won O-C dataset for 25 XRBs and the careful discussion of individual timing systematics. The weakness is the regression-on-residual construction of the universal AML law, which creates a real risk that the reported Mdot dependence is an artifact of correlated errors. I would not accept the universal-law claim without a synthetic-error test and an out-of-sample or cross-validated check. The paper also overstates 'universal' relative to what Table 8 actually fits."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read it for the data, not for the law. Schaefer has put together a genuinely useful catalog of orbital period changes for 25 XRBs, with carefully built O-C diagrams and new timings from TESS, the Harvard plates, and AAVSO. Combining those with 52 CVs and testing the magnetic braking model is legitimate, and the MBM's unique-track claim is convincingly dead for these systems. The GR and mass-transfer subtraction is standard, the worked examples are instructive, and the citation pattern is fine: he engages Knigge et al. 2011, Patterson, and the individual system papers.\n\nThe soft spot is exactly where your reader put it. Equation 2 defines Pdot_AML by subtracting a term proportional to Mdot, then Equation 14 regresses that residual on Mdot. With Mdot uncertainties 'typically like a factor of ten' (the paper's own Section 5), correlated errors between the subtracted term and the regressor can produce a positive delta even when the true AML has no Mdot dependence. The same contamination shifts beta and gamma through q. Adding a 0.25-dex systematic to force reduced chi-square near unity does not cure bias in the central value; it only inflates the error bars. So the headline delta ~ 0.43 is not established as a property of nature. The law is also fit to the same 77 systems it claims to describe, with no out-of-sample test. The below-gap fit is especially fragile: huge exponent uncertainties, a reduced chi-square well below unity, and the arbitrary treatment of non-parabolic O-C systems as Pdot=0 with ad hoc uncertainties.\n\nNone of this undercuts the Pdot measurements or the MBM failure. But the universal-law claim is too strong. The paper deserves a serious referee. The referee should ask for a synthetic-error test: inject noisy Mdot, show whether the regression recovers delta=0 when the true law has delta=0. An out-of-sample check on a few systems would also help. Without that, the universal law is a fitting result, not a measurement.\n\nFor binary evolution modelers and CV/XRB observers, the Pdot table is worth having. I'd bring the paper to reading group and would cite the table; I would not yet cite the law.","headline":"A valuable Pdot catalog that refutes magnetic braking, but the claimed universal AML law is a same-sample fit whose Mdot exponent is likely contaminated by correlated errors.","tokens_in":62042,"tokens_out":4018,"would_cite":true,"duration_ms":46149,"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":"This paper measures orbital period changes for 25 X-ray binaries and, combined with 52 cataclysmic variables, derives an empirical power-law law for angular momentum loss that it claims is universal across accreting binaries.","keywords":["X-ray binaries","cataclysmic variables","orbital period change","angular momentum loss","magnetic braking","O-C diagrams","binary evolution","empirical power law"],"falsifier":"Take a system with period near 0.5 day whose masses are measured by eclipses and radial velocities and whose accretion rate is measured independently from X-ray luminosity and a Gaia distance. If its $\\dot P$ differs from the Equation 16 prediction by more than the propagated uncertainties, or if a decade-long change in $\\dot M$ does not produce the predicted power-law change in $\\dot{P}_{\\rm AML}$, the universal law would be falsified.","tokens_in":60890,"feed_emoji":"💫","tokens_out":4863,"duration_ms":53213,"temperature":0.7,"pith_summary":"This paper tries to establish that the slow orbital decay of accreting binaries—both X-ray binaries and cataclysmic variables—is governed by a single empirical angular-momentum-loss law, and that the long-standing magnetic braking model is ruled out by direct measurements. It reports new period-change measures for 25 X-ray binaries from O−C diagrams spanning decades and combines them with 52 cataclysmic variables. After subtracting gravitational radiation and mass-transfer contributions, the residual period change $\\dot{P}_{\\rm AML}$ is fitted to a power law in orbital period, stellar masses, and accretion rate, yielding one law for periods 0.13–1.0 days and analogous laws for shorter and longer periods. If correct, the three laws together describe the evolution of all 77 systems and give population-synthesis modelers a replacement for the magnetic braking recipe.","feed_headline":"77 binaries obey one empirical orbital-decay law","feed_subtitle":"Period-change data from X-ray binaries and cataclysmic variables reveal a universal angular-momentum-loss law.","key_machinery":"The central tool is the O−C diagram: plots of observed minus calculated eclipse or minimum times fitted with parabolas, whose curvature gives the steady period change $\\dot P$. Equations 1–11 then subtract the well-known gravitational-radiation contribution $\\dot P_{\\rm GR}$ and the mass-transfer contribution $\\dot P_{\\rm mt}$ from the measured $\\dot P$ to isolate the residual $\\dot P_{\\rm AML}$. Equation 14 converts the resulting 84 measures into a chi-square fit of a power law in $P$, $M_{\\rm prim}$, $M_{\\rm comp}$, and $\\dot M$, producing the fitted exponents in Table 8. The load-bearing step is treating the residual $\\dot P_{\\rm AML}$ as a real extra angular-momentum-loss signal rather than an artifact of uncertain accretion rates.","core_discovery":"The central claim is that the dominant angular momentum loss in accreting binaries is empirically $\\dot{P}_{\\rm AML} = -1500\\times 10^{-12}\\, P^{1.29} M_{\\rm prim}^{2.75} M_{\\rm comp}^{-1.00} \\dot{M}_{-8}^{0.43}$ for periods from 0.13 to 1.0 days, with separate fitted power laws for binaries below the period gap and for binaries with $P>1$ day. The same 77 systems show that magnetic braking predictions are wrong by more than an order of magnitude for most systems. The paper concludes that the unknown AML mechanism is controlled by the accretion process, because $\\dot{P}_{\\rm AML}$ rises with accretion rate, and that the empirical family of laws can stand as a universal description of binary evolution until a physical mechanism is identified.","pith_inferences":["If Equation 16 is causal, systems whose accretion rate changes sharply, like the recurrent nova U Sco after its 2010 eruption, should show a corresponding power-law change in $\\dot{P}_{\\rm AML}$; the paper reports such a jump but leaves it unexplained, so watching the next eruption would directly test the $\\dot{M}^{0.43}$ term.","The steep positive exponent on $M_{\\rm prim}$ predicts that, at equal period and accretion rate, binaries with more massive white dwarfs or neutron stars lose orbital angular momentum faster; this is a testable ranking within existing eclipsing systems.","If the universal law holds, the minimum period and period gap of cataclysmic variables should be derivable from Equations 15–16 rather than from magnetic braking physics, making the observed period distribution an independent check on the fitted exponents."],"forward_implications":["The magnetic braking model's single evolutionary track is contradicted: 7 of 8 X-ray binaries with main-sequence companions and most cataclysmic variables deviate from its predictions by orders of magnitude.","Evolution and population-synthesis calculations can replace the magnetic braking recipe with Equations 15–16, which reproduce the observed $\\dot{P}_{\\rm AML}$ scatter to about 0.33 dex.","Because $\\dot{P}_{\\rm AML}$ depends on the accretion rate, the dominant loss mechanism must live in the accretion flow, stream, or boundary layer, not in the companion's magnetic wind.","Below the period gap, at least 6 of 18 systems show a non-zero $\\dot{P}_{\\rm AML}$, so gravitational radiation alone does not drive those binaries.","The period gap and minimum period, long cited as successes of magnetic braking, are no longer evidence for that mechanism."],"supporting_citations":[{"why":"Supplies the revised magnetic braking model predictions and the single-evolutionary-track structure that the paper tests and rejects.","marker":"Knigge et al. (2011)"},{"why":"Provides the 52 cataclysmic variable $\\dot P$ measures that make up half of the 77-system sample.","marker":"Schaefer (2024)"},{"why":"Original magnetic braking formulation with the adopted gamma=3 power law that the paper compares against.","marker":"Rappaport, Verbunt, & Joss (1983)"},{"why":"Equation 8 gives the magnetic braking $\\dot P_{\\rm mb}$ prediction used in the Section 7.2 comparison.","marker":"Paxton et al. (2015)"},{"why":"CAML and e-CAML conjectured angular-momentum-loss laws that the paper fits and strongly rejects.","marker":"Schreiber, Zorotovic, & Wijnen (2016)"},{"why":"Prior collection of X-ray binary $\\dot P$ values and wind accretion rates that the paper extends and uses for the HMXB sample.","marker":"Falanga et al. (2015)"},{"why":"Standard prescription for how mass transfer and gravitational radiation contribute to the observed period change.","marker":"Frank, King, & Raine (2002)"}],"fun_headline_variants":["Universal orbital decay law from 77 X-ray binaries","Empirical law links angular momentum loss in 77 binaries","77 binaries reveal universal angular momentum loss law","One empirical law governs orbital decay in 77 binaries","Magnetic braking fails; empirical orbital decay law wins"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The fitted law assumes that the accretion rates and wind-capture efficiencies in Table 6 are accurate enough, despite uncertainties of roughly a factor of ten in $\\dot M$ and poorly known wind efficiencies, that the residual $\\dot{P}_{\\rm AML}$ after subtracting GR and mass transfer is a real signal rather than an artifact of those errors.","fun_headline_variants_meta":{"raw":{"variants":["Universal orbital decay law from 77 X-ray binaries","Empirical law links angular momentum loss in 77 binaries","77 binaries reveal universal angular momentum loss law","One empirical law governs orbital decay in 77 binaries","Magnetic braking fails; empirical orbital decay law wins"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000578,"raw_usage":{"total_tokens":2809,"prompt_tokens":1111,"completion_tokens":1698,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":727,"completion_tokens_details":{"reasoning_tokens":1623}},"tokens_in":727,"tokens_out":1698,"duration_ms":12556,"temperature":1.0,"reasoning_tokens":1623,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T16:23:30.098642+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a system with period near 0.5 day whose masses are measured by eclipses and radial velocities and whose accretion rate is measured independently from X-ray luminosity and a Gaia distance. If its $\\dot P$ differs from the Equation 16 prediction by more than the propagated uncertainties, or if a decade-long change in $\\dot M$ does not produce the predicted power-law change in $\\dot{P}_{\\rm AML}$, the universal law would be falsified.","supporting_citations":[{"cited_title":"2011, ApJS, 194, 28 (K2011)","cited_arxiv_id":null,"evidence_quote":"Supplies the revised magnetic braking model predictions and the single-evolutionary-track structure that the paper tests and rejects."}],"review_version":1}