{"id":"fec21155-70b0-4996-b911-7f5cda96fd1c","arxiv_id":"2610.27164","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The paper identifies the dominant sources of discrepancy between black-hole perturbation theory and numerical relativity across the waveform, showing that inspiral errors are largely dissipative, precessing systems have a new nonlinear dependence, and merger errors come from the change in…","lead":"Using a modular time-domain Teukolsky framework, the paper isolates which finite-mass-ratio corrections matter most in binary black hole inspirals and mergers. For nonspinning systems, a second-order dissipative flux removes 87% of the inspiral rate error; for precessing systems, a new ν² sin²ι dependence reveals missing precessional dynamics, while the merger offset arises from the change in spin of the final black hole.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek/deepseek-v4-flash","headline":"The identification of ν² sin²ι dependence as a finite-mass-ratio correction to precessional dynamics rests solely on a fit to 0PA tier waveforms without testing whether including available conservative corrections (C) modifies this dependence.","rationale":"The paper is thorough and the central claims for nonspinning and aligned-spin binaries are well-supported. The precessing result is the least tested component of the central narrative; it is based on a single-tier fit without cross-checking at higher perturbative order. The proposed test would clarify robustness. However, the paper's language is appropriately tentative ('suggesting'), and the core contributions (framework, catalog, validation, merger explanation) are unaffected. The verdict remains ACCEPT.","tokens_in":30237,"tokens_out":3672,"duration_ms":39677,"concrete_test":"For 5–10 precessing binaries spanning the ι and q range used in the fit, repeat the β̂ measurement using the 0PA-H-C tier (including horizon absorption and conservative correction). If the fitted A₃ coefficient changes by more than its stated uncertainty (3.74 ± 0.33), the interpretation as a robust finite-mass-ratio precession correction would require revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper claims that for precessing binaries, a ν² sin²ι dependence in the inspiral time rescaling β̂ is not explained by leading orbit-averaged PN terms, 'suggesting finite-mass-ratio corrections to precessional dynamics' (b4, b98). This interpretation is based entirely on the bare adiabatic (0PA) tier (Eq. 5), which omits the conservative correction C, horizon absorption H, and secondary spin. The paper has these ingredients available for generic inclined Kerr orbits (b31: metric reconstruction for redshift). If the ν² sin²ι dependence arises from truncation error in the dynamics rather than genuine precessional corrections, the claim would be misattributed. The paper does not test whether including C (or H) shifts the fitted A₃ coefficient. Since the conservative correction is known to alter β̂ by ~0.333/ν for aligned spins (b103), a similar effect on the precessing coefficient is plausible but unexamined. This does not invalidate the paper's core results for nonspinning and aligned-spin binaries, but it weakens the precessing component of the central narrative.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek/deepseek-v4-flash","summary":"This paper develops a modular time-domain Teukolsky framework (gwKerrDynamics + gwTeukolsky) that isolates finite-mass-ratio corrections to binary-black-hole dynamics while keeping the radiation at first perturbative order. By comparing 1571 time-domain waveforms (the Vishveshwara Catalog) with numerical relativity across 348 binaries, the authors identify which perturbative ingredients matter at each stage of the coalescence. For nonspinning binaries, the second-order dissipative flux removes 87% of the inspiral rate discrepancy; for precessing binaries, a nu^2 sin^2 iota dependence in the time rescaling not explained by leading orbit-averaged PN terms suggests finite-mass-ratio corrections to precessional dynamics; near merger, the offset is quantitatively accounted for by the difference between the initial and remnant Kerr backgrounds. The paper includes extensive external validation against independent codes (pybhpt, WaSABI, published self-force data) and resolution studies.","tokens_in":31803,"tokens_out":1204,"duration_ms":35540,"significance":"If the central claims hold, this work provides a clear physical picture of where and why black-hole perturbation theory begins to differ from numerical relativity across the mass-ratio and spin parameter space. The finding that the dominant missing physics changes character from second-order dissipative forcing during inspiral to the background Kerr-spacetime change near merger is conceptually important for waveform modeling and for understanding the limits of perturbative descriptions. The modular framework and the public catalog of tiered waveforms are valuable community resources. The honest discussion of limitations (Schwarzschild F1 applied to Kerr, non-independent 1PA comparison, etc.) strengthens the paper's credibility.","major_comments":[{"comment":"The claim that the nu^2 sin^2 iota dependence in the inspiral time rescaling beta_hat suggests finite-mass-ratio corrections to precessional dynamics is based solely on the 0PA tier (bare adiabatic dynamics). The authors acknowledge elsewhere that the conservative correction C can alter beta_hat by about 0.333/nu for aligned spins. Because the framework has C, H, and other ingredients available for generic inclined orbits (via metric reconstruction, b31), the paper should either (a) test whether including C (or H) shifts the fitted A3 coefficient, or (b) explicitly discuss this as an unresolved limitation of the precessing part of the analysis. As written, the evidence for the precessional interpretation is weaker than for the nonspinning and aligned-spin results, and the manuscript should address this asymmetry.","section":"Eq. (5) and surrounding text (b4, b93-b98)","authors_acknowledge":null,"authors_acknowledge_p":null}],"minor_comments":[{"comment":"The caption says 'Symbols mark the configurations used in the analysis, with filled circles in the right panel showing the measured SXS merger values.' It would be helpful to also state what the open symbols represent (e.g., BHPT values, or configurations without NR data).","section":"Fig. 2 caption","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The closed-form fit coefficients a_k, c_k, p_k in Eq. (C1) are defined but their numerical values are not given. While the paper states these are closed-form approximations and production waveforms use direct interpolation, providing the coefficient values (even in a table) would be useful for reproducibility.","section":"Appendix C","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The 1PA comparison with WaSABI is described as not an independent validation because the inputs are derived from WaSABI data (b250). This is correctly noted, but the sentence could be moved to the figure caption for immediate transparency.","section":"Appendix D, Fig. 6","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The abstract states that for precessing binaries a nu^2 sin^2 iota dependence suggests finite-mass-ratio corrections to precessional dynamics. The qualifier 'suggesting' is appropriate, but given the major comment above, the paper should be more explicit that this inference rests on the 0PA tier and has not been tested with conservative corrections.","section":"Abstract and b4","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The manuscript is generally well organized and clearly written. The use of AI-assisted development is transparently documented in Appendix G. A few very extended sentences could be broken for readability (e.g., the 10-line sentence in the third paragraph of the Introduction).","section":"General","authors_acknowledge":null,"authors_acknowledge_p":null}],"recommendation":"minor_revision","confidential_remarks":"The skeptic's concern about the precessing claim is legitimate but proportionate. The paper's central results for nonspinning and aligned-spin binaries are robust and well validated. The precessing claim is already tentative, but the manuscript should either test it or explicitly flag the missing check. Overall this is a strong paper that merits publication after addressing the single major comment. The level of validation and the public release of the catalog set a high standard. The AI-development discussion is unusually transparent and helpful."},"author_rebuttal":null,"desk_editor":{"model":"deepseek/deepseek-v4-flash","letter":"This is a solid, well-executed paper that delivers what it promises: a clean modular decomposition of where finite-mass-ratio corrections enter binary-black-hole dynamics. The central result—that 87% of the inspiral-rate discrepancy for nonspinning binaries comes from the second-order dissipative flux—is convincingly demonstrated, and the explanation of merger offsets via the change in Kerr background is physically clear. The framework itself (gwKerrDynamics + gwTeukolsky) is extensively validated against independent codes (pybhpt, WaSABI) and analytic limits, and the 1571-waveform catalog is a genuine community resource for LISA-era modeling. The paper is also refreshingly honest about its limitations: the Schwarzschild F1 flux applied to Kerr, the non-independent 1PA comparison, and the convergence floor near merger are all stated plainly.\n\nThe strongest claim—that precessing binaries show a ν² sin²ι dependence suggesting finite-mass-ratio corrections to precessional dynamics—rests on a fit to the bare adiabatic (0PA) tier only. The stress-test concern is real: the paper has the conservative correction C available for generic inclined orbits (metric reconstruction for redshift), but never tests whether including C shifts the fitted A₃ coefficient. For aligned spins, C changes β̂ by ~0.333/ν (b103), so a similar effect on the precessing coefficient is plausible. This does not invalidate the core results for nonspinning and aligned-spin binaries, but it means the precessing narrative is weaker than the rest. The paper should either perform that test or explicitly temper the claim.\n\nOverall, this is a significant contribution that deserves serious peer review. For nonprecessing systems the evidence is strong; the precessing part needs one additional cross-check before the conclusions can be trusted at the same level. I would accept the paper conditional on the author adding a test of whether including conservative corrections modifies the ν² sin²ι fit for precessing binaries.\n\nFor the reading group: yes, this will generate good discussion about modularity and the limits of BHPT.\n\nI would cite this in my own work on waveform modeling.","headline":"A thorough and honest modular decomposition of finite-mass-ratio effects; the precessing claims are less supported than the rest.","tokens_in":31110,"tokens_out":1351,"would_cite":true,"duration_ms":15794,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"pith_extraction":{"msc":[],"pacs":["04.30.-w","04.25.dg","04.25.Nx","04.70.Bw"],"model":"deepseek/deepseek-v4-flash","headline":"A modular Teukolsky framework isolates the distinct physical mechanisms that drive the departure of black-hole perturbation theory from numerical relativity across inspiral, precession, and merger.","keywords":["black-hole perturbation theory","numerical relativity","gravitational waves","Teukolsky equation","self-force","binary black holes","time-domain waveforms"],"falsifier":"A full Kerr second-order self-force computation for an equal-mass, aligned-spin binary with χ1z = χ2z = 0.6 (the configuration of SXS:BBH:2509) would either confirm that the residual dephasing of order 3 rad after including the Schwarzschild F1 flux is indeed removed by the Kerr second-order flux, or else show that other corrections (such as higher-order conservative effects) are needed.","tokens_in":31735,"feed_emoji":"🌀","tokens_out":2678,"duration_ms":28233,"temperature":0.7,"pith_summary":"This paper aims to identify which finite-mass-ratio corrections cause black-hole perturbation theory (BHPT) to deviate from numerical relativity (NR) at different stages of a binary black-hole merger. The author develops a modular time-domain Teukolsky framework that adds individual corrections—horizon absorption, conservative self-force, second-order dissipation, secondary-spin effects—to the orbital dynamics while keeping the radiative field at first perturbative order. By comparing the resulting time-domain waveforms with NR for 348 binaries, the paper shows that the dominant missing physics changes character across the waveform: during inspiral it is the dissipative forcing, for precessing systems it is a new nonlinear dependence on tilt that is not captured by standard post-Newtonian theory, and near merger it is the difference between the initial Kerr background and the remnant black hole. The work provides a unified picture of where and why BHPT departs from full relativity, and it releases 1571 waveforms as the Vishveshwara Catalog.","feed_headline":"Three distinct nonlinearities limit black-hole perturbation theory","feed_subtitle":"A modular Teukolsky framework reveals separate missing physics for inspiral, precession, and merger, each dominant in its own regime.","key_machinery":"The paper's central tool is a modular time-domain Teukolsky pipeline (gwKerrDynamics + gwTeukolsky) that allows finite-mass-ratio corrections to be added independently to the orbital dynamics while holding the first-order radiative calculation fixed. The key diagnostic is the time rescaling parameter β, defined by t_NR = β t_Teuk, which measures the rate error relative to NR. The analysis also uses the second-order dissipative flux computed in Schwarzschild (from the Black Hole Perturbation Toolkit), the conservative redshift invariant, horizon-absorption fluxes, and a Newtonian mass-asymmetry factor for each mode. The decomposition into separate, switchable corrections lets the author identify the dominant missing physics at each stage by its effect on β and the waveform phase.","core_discovery":"The central claim is that the discrepancy between first-order time-domain Teukolsky waveforms and numerical relativity has a physically distinct structure at each stage of a binary black-hole coalescence. During the inspiral of nonspinning binaries, adding the known second-order Schwarzschild dissipative flux removes 87% of the rate discrepancy, bringing it to within about 1% of NR. For precessing binaries, a $ν^{2}$ $sin^{2}$ ι dependence appears in the inspiral time rescaling that is not explained by the leading orbit-averaged post-Newtonian terms, pointing to finite-mass-ratio corrections to precessional dynamics. Near merger, the same second-order dissipation removes only 1–2% of the offset; instead, 96% of the predicted frequency shift arises from the change in spin between the initial and remnant Kerr black holes, as the Teukolsky field continues to ring at the frequency of the initial background. These three regimes—forcing, precession, and background mismatch—together account for the dominant nonlinearity that limits the perturbative description.","pith_inferences":["The success of the Schwarzschild second-order flux for nonspinning binaries, together with its partial failure for spinning primaries, implies that a full Kerr second-order self-force calculation is likely the single most impactful near-term target for improving perturbative waveform accuracy at comparable masses.","The result that the merger discrepancy is dominated by the spin difference between initial and remnant Kerr backgrounds suggests that a simple analytic correction—shifting the Teukolsky ringdown frequency to match the remnant's quasinormal mode—could rescue the late-time waveform without a full second-order evolution.","If the precessional ν^2 sin^2 ι term is confirmed by future NR simulations at higher mass ratios, it would indicate that orbit-averaged post-Newtonian theory systematically underestimates the strength of finite-mass-ratio precession effects, with implications for parameter estimation of precessing intermediate-mass-ratio binaries.","The modular approach itself could be applied to other problems in gravitational-wave astronomy—for example, isolating the effects of eccentricity or gas drag in extreme-mass-ratio inspirals—by keeping the field evolution fixed and varying only the orbital dynamics."],"forward_implications":["If the inspiral discrepancy is largely captured by a single O(ν) dissipative correction, then extending the second-order self-force calculation to Kerr backgrounds—especially for spinning primaries—should bring BHPT waveforms for aligned-spin binaries into percent-level agreement with NR across the whole inspiral.","The ν^2 sin^2 ι dependence in precessing binaries suggests that a finite-mass-ratio correction to the precession of the orbital plane about the total angular momentum is needed; including this effect could significantly improve perturbative models for generic misaligned binaries.","Near merger, the dominant correction is not a self-force effect but rather a change of background: the Teukolsky field must be evolved on the evolving Kerr spacetime that eventually becomes the remnant. This motivates using a time-dependent Kerr background or a transition–plunge prescription that tracks the black hole's mass and spin evolution.","The released Vishveshwara Catalog of 1571 time-domain waveforms provides a systematic testbed for waveform models and for training surrogate or machine-learned corrections that interpolate between BHPT and NR across the mass-ratio–spin plane."],"supporting_citations":[{"why":"Provides the second-order dissipative flux in Schwarzschild that is the key correction for nonspinning inspiral.","marker":"[23]"},{"why":"Supplies the second-order energy flux for nonspinning binaries, used as the F1 correction.","marker":"[25]"},{"why":"Previous surrogate model for nonspinning binaries calibrated to NR; provides baseline for comparison and validation.","marker":"[34]"},{"why":"The SXS catalog supplies the numerical-relativity waveforms against which all comparisons are made.","marker":"[92]"},{"why":"The third SXS catalog extends the NR sample, used for the merger and ringdown analysis.","marker":"[93]"},{"why":"Provides the gravitational self-force data for conservative corrections on generic Kerr orbits.","marker":"[71]"},{"why":"The Compère–Küchler transition-to-plunge prescription is used to connect inspiral to plunge.","marker":"[97]"},{"why":"Establishes the time-domain Teukolsky source-modeling approach that this work extends.","marker":"[107]"},{"why":"Provides the ISCO estimate for final mass and spin used to predict the merger frequency shift.","marker":"[135]"}],"fun_headline_variants":["Three distinct nonlinearities limit black-hole theory","Inspiral, precession, merger: three regimes of missing physics","Black-hole perturbation theory hits three separate walls","The three shapes of nonlinearity in binary-black-hole dynamics","Where perturbative gravity breaks down: three regimes"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The second-order dissipative flux is computed in Schwarzschild but applied to Kerr backgrounds by evaluating it at the matched orbital frequency, without a full Kerr second-order self-force calculation; this assumption directly affects the claimed 87% reduction for nonspinning binaries versus the weaker improvement for spinning primaries, and the interpretation of residual spin dependence as evidence for absent Kerr second-order physics.","fun_headline_variants_meta":{"raw":{"variants":["Three distinct nonlinearities limit black-hole theory","Inspiral, precession, merger: three regimes of missing physics","Black-hole perturbation theory hits three separate walls","The three shapes of nonlinearity in binary-black-hole dynamics","Where perturbative gravity breaks down: three regimes"]},"model":"deepseek/deepseek-v4-flash","effort":"low","cost_usd":6.6e-05,"raw_usage":{"cost":6.5508352e-05,"is_byok":false,"cost_details":{"upstream_inference_cost":6.5508352e-05,"upstream_inference_prompt_cost":3.3503232e-05,"upstream_inference_completions_cost":3.200512e-05},"total_tokens":1368,"prompt_tokens":992,"completion_tokens":376,"prompt_tokens_details":{"audio_tokens":0,"video_tokens":0,"cached_tokens":256,"cache_write_tokens":0},"completion_tokens_details":{"audio_tokens":0,"image_tokens":0,"reasoning_tokens":300}},"tokens_in":992,"tokens_out":376,"cost_source":"openrouter_usage","duration_ms":4130,"image_pages":[],"temperature":1.0,"billed_attempts":1,"reasoning_tokens":300,"transport_provider":"openrouter","cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-10-21T15:08:32.182960+00:00","model_set":{"reader":"deepseek/deepseek-v4-flash"},"falsifier":"A full Kerr second-order self-force computation for an equal-mass, aligned-spin binary with χ1z = χ2z = 0.6 (the configuration of SXS:BBH:2509) would either confirm that the residual dephasing of order 3 rad after including the Schwarzschild F1 flux is indeed removed by the Kerr second-order flux, or else show that other corrections (such as higher-order conservative effects) are needed.","supporting_citations":[{"cited_title":"Gravitational-Wave Extraction from an Inspiraling Configuration of Merging Black Holes,","cited_arxiv_id":null,"evidence_quote":"Provides the second-order dissipative flux in Schwarzschild that is the key correction for nonspinning inspiral."},{"cited_title":"Gravitational waveforms for compact binaries from second-order self-force theory","cited_arxiv_id":"2112.12265","evidence_quote":"Supplies the second-order energy flux for nonspinning binaries, used as the F1 correction."},{"cited_title":"A Surrogate Model for Gravitational Wave Signals from Comparable- to Large- Mass-Ratio Black Hole Binaries","cited_arxiv_id":"1910.10473","evidence_quote":"Previous surrogate model for nonspinning binaries calibrated to NR; provides baseline for comparison and validation."},{"cited_title":"Im- proved filters for gravitational waves from inspiralling com- pact binaries,","cited_arxiv_id":null,"evidence_quote":"The SXS catalog supplies the numerical-relativity waveforms against which all comparisons are made."},{"cited_title":"The SXS Collaboration catalog of binary black hole simulations","cited_arxiv_id":"1904.04831","evidence_quote":"The third SXS catalog extends the NR sample, used for the merger and ringdown analysis."},{"cited_title":"The mass and angular momentum of reconstructed metric perturbations","cited_arxiv_id":"1702.00969","evidence_quote":"Provides the gravitational self-force data for conservative corrections on generic Kerr orbits."},{"cited_title":"Transition from adiabatic inspiral to plunge into a spinning black hole","cited_arxiv_id":"1101.3749","evidence_quote":"The Compère–Küchler transition-to-plunge prescription is used to connect inspiral to plunge."},{"cited_title":"An Approach to gravita- tional radiation by a method of spin coefficients,","cited_arxiv_id":null,"evidence_quote":"Establishes the time-domain Teukolsky source-modeling approach that this work extends."},{"cited_title":"LIGO Algorithm Library - LALSuite,","cited_arxiv_id":null,"evidence_quote":"Provides the ISCO estimate for final mass and spin used to predict the merger frequency shift."}],"review_version":1}