{"id":"ec57c36a-e36d-4c3d-910d-59b027adaf0a","arxiv_id":"2507.01156","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"RIFT can recover the masses, spins, and orbital energy and angular momentum of simulated hyperbolic black hole encounters when paired with the TEOBResumSDALI waveform model.","lead":"Researchers adapted an existing gravitational-wave analysis tool, RIFT, to estimate the properties of black holes that whip past each other in unbound hyperbolic orbits instead of spiraling in. The tool recovered simulated masses and orbit parameters for two events, a scatter and a plunge, seen by a future Cosmic Explorer detector.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"In both zero-noise injections, the reported 90% credible intervals for mass ratio q exclude the injected value q=1, directly contradicting the claimed accurate recovery of mass parameters.","rationale":"The reader's weakest assumption focused on the same-model injection and the TEOBResumSDALI model's fidelity, which is a legitimate concern for astrophysical applicability. The specific q exclusion is more immediately damaging because it is visible from the paper's own reported numbers: with zero noise and injection equal to template, the true parameter must lie inside a well-calibrated posterior. The reader did note the mass-ratio prior inconsistency in the rationale, so there is partial agreement, but the weakest_assumption field identified a different issue. My proposed check is deliberately narrow: correct the prior convention and re-run the same two injections. This isolates whether the q bias is a recoverable setup error or a genuine failure of the parameter-estimation scheme. The verdict remains CONDITIONAL because the paper is presented as a proof of principle and the flaw is potentially fixable, but the authors must address the q bias before the 'first comprehensive infrastructure' claim can be accepted. If q=1 remains outside the 90% interval after the prior correction, the central accuracy claim should be rejected or substantially weakened.","tokens_in":13193,"tokens_out":4588,"duration_ms":55741,"concrete_test":"Re-run both zero-noise injections with a corrected, internally consistent mass-ratio prior (for example q=m_low/m_high ∈ [0.1,1]) using the same RIFT configuration, then form the 90% credible interval of the marginalized q posterior. If q=1 is again excluded, the pipeline has a genuine bias and the abstract's accuracy claim is unsupported; if q=1 is included after the prior correction, the current result is a prior/convention error that must be stated and corrected.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that RIFT with TEOBResumSDALI \"accurately recovers the mass, spins, and hyperbolic orbit parameters,\" with two zero-noise injections as the sole quantitative support. In both cases the injected mass ratio q=1 is outside the reported 90% credible interval: the scatter gives q=0.8188^{+0.1574}_{-0.2586}, or [0.5602, 0.9762], and the plunge gives q=0.8113^{+0.1679}_{-0.2528}, or [0.5585, 0.9792]. Because these are zero-noise injections and the recovery uses the same waveform model used to generate the signal, a correct pipeline and prior should put q=1 inside the posterior credible region. Its exclusion indicates a prior inconsistency or pipeline error, not merely waveform-model systematics. This is supported by an internal inconsistency in the Methods: the prior is stated as q∈{1.0, 10.0} with m1≥m2, yet the reported posteriors are below q=1, and the definition q=m2/m1 with m1≥m2 would have q≤1. The abstract's accuracy claim is therefore not established for at least one mass parameter in either demonstrated case.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper adapts the RIFT parameter-estimation code to the TEOBResumSDALI effective-one-body waveform model for hyperbolic black-hole encounters. It describes the intrinsic and extrinsic parameter space, discusses data-conditioning challenges for scatter, capture, and plunge waveforms, and presents zero-noise injection-recovery tests for one scatter and one plunge event at Cosmic Explorer design sensitivity, each with SNR ~ 42 and only the dominant (2,±2) mode. The authors report good recovery of total mass, energy, angular momentum, and effective spin, but the recovered mass ratio is centered near 0.82 for equal-mass injections in both cases. They conclude that parameter estimation of generic hyperbolic waveforms is now possible with RIFT.","tokens_in":13482,"tokens_out":6271,"duration_ms":60961,"significance":"The adaptation of a mature grid-based PE code to hyperbolic-encounter waveforms is a useful step, and if validated it would provide a practical tool for third-generation detectors. The paper's strengths include the use of a physically broad EOB model, explicit attention to waveform-class boundaries such as the plunge/scatter separatrix, and a clean zero-noise testbed that isolates pipeline behavior. However, the validation is self-referential (the same model generates the injections and serves as the recovery template), it uses only two examples in a single optimal detector configuration, and the mass-ratio results conflict with the abstract's accuracy claim. The infrastructure may well be sound, but the quantitative evidence presented does not yet establish unbiased recovery of all mass parameters.","major_comments":[{"comment":"For both zero-noise injections with q=1, the reported 90% credible intervals for q exclude the true value: the scatter gives q=0.8188^{+0.1574}_{-0.2586}, i.e. [0.560,0.976], and the plunge gives q=0.8113^{+0.1679}_{-0.2528}, i.e. [0.559,0.979]. Since the injections are noise-free and the recovery model is identical to the injection model, a correct prior and pipeline should place the true value well inside the credible interval. The abstract's claim that RIFT 'accurately recovers the mass' is therefore not supported for the mass ratio in either demonstrated case; the authors should re-examine the prior and sampling, or substantially soften the claim.","section":"Results, mass-ratio paragraph and Fig. 2"},{"comment":"The text defines q ≡ m2/m1 with m1 ≥ m2, which implies q ≤ 1, but then states the prior 'q∈{1.0,10.0}'. This is internally inconsistent. The reported posteriors (q < 1, with upper 90% bounds below 1) are consistent with a prior truncated at q=1, not with the stated q∈[1,10]. Please correct the prior definition or the reported range, and rerun the injections; the mass-ratio bias may be an artifact of this inconsistency.","section":"Methods, prior specification"},{"comment":"Because the injected signals are generated with the same TEOBResumSDALI model used for recovery and contain only the dominant (2,±2) mode, the study tests the self-consistency of the pipeline, not whether the model is faithful to the true waveforms from hyperbolic encounters. The closing-remarks claim that 'parameter estimation of generic hyperbolic waveforms is now possible' should be qualified accordingly, particularly since the authors themselves note that the (2,0) mode can be important for scattering and that higher-order modes would help break the mass-ratio degeneracy.","section":"Methods and Closing remarks"}],"minor_comments":[{"comment":"The text reads 'Finallyplunge events' and should be 'Finally, plunge events'.","section":"Methods, first paragraph"},{"comment":"The caption reads 'denote the the 5% and 95% percentiles'; the duplicate 'the' should be removed.","section":"Fig. 3 caption"},{"comment":"Reference [74] is cited for the statement that the (2,0) mode can be important in scattering systems, but the listed reference is about GW190521 as a merger of Proca stars and does not appear to support this statement; please verify and, if needed, replace the citation.","section":"References"},{"comment":"The phrase 'hyperbolic orbit parameters: the system energy and angular momentum defined at a fiducial initial separation' has an awkward colon construction; consider rephrasing as 'namely, the system energy and angular momentum defined at a fiducial initial separation'.","section":"Abstract"}],"recommendation":"major_revision","confidential_remarks":"The mass-ratio prior inconsistency and the exclusion of q=1 from both posterior credible intervals are the key blockers. The authors should resolve the parameter convention, rerun the injections, and then reassess whether the abstract's accuracy claim is warranted. A single-detector, optimally oriented, zero-noise pair of injections is also a narrow evidentiary base for the 'first comprehensive infrastructure' claim; expanding the study or substantially qualifying the claim would strengthen the paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Hi [Colleague],\n\nQuick read on arXiv:2507.01156. The headline: this is a legit first step toward PE for hyperbolic encounters, but the paper's central accuracy claim is undercut by a mass-ratio problem the authors don't acknowledge.\n\nWhat's new: they get RIFT working with TEOBResumSDALI on scatter and plunge waveforms, including the practical fixes—tapering for asymmetric starts, a peak-finder to locate merger time for captures. That's real, useful infrastructure work, and the paper is honest about its scope: two zero-noise injections, one detector, dominant (2,±2) mode, same model used for injection and recovery. The discussion of the scatter/plunge separatrix and the degeneracies in E0/M and p_phi is sensible.\n\nThe soft spot is load-bearing. In both injections (q=1), the 90% credible interval for q excludes 1. For scatter, q=0.8188^{+0.1574}_{-0.2586}; for plunge, q=0.8113^{+0.1679}_{-0.2528}. In a zero-noise injection with the same model, a correct pipeline should put the true value inside the credible region. And the Methods state a prior q∈{1.0,10.0} with q≡m2/m1 and m1≥m2, which is self-contradictory—q cannot exceed 1 under that definition. So either the prior is mis-specified or there's a bug in the sampling or in the reported values. The abstract's claim that RIFT 'accurately recovers the mass, spins, and hyperbolic orbit parameters' is not supported for the mass ratio in either demonstrated case.\n\nEverything else—total mass, E0/M, p_phi, spin—recovers reasonably well, and the authors flag the known limitations (higher modes, the (2,0) mode). But the q inconsistency needs to be fixed and the recovery re-run before the accuracy claim holds. I'd also want at least a few noise realizations and a non-equal-mass injection.\n\nVerdict: worth sending to a referee, because the infrastructure is novel and the authors are clearly not hiding the scaffolding. But it needs major revision, not a quick accept.\n\nWho's this for? GW data-analysis people building PE for eccentric/hyperbolic populations, and people planning CE-era searches. I'd cite it once the q issue is resolved.\n\nRecommendation: send to peer review with a request to fix the prior, redo the injections, and soften the abstract to match what's actually demonstrated.","headline":"Promising proof-of-principle for hyperbolic-encounter PE with RIFT, but the mass-ratio prior inconsistency and biased recovery undermine the central accuracy claim.","tokens_in":13960,"tokens_out":2304,"would_cite":false,"duration_ms":24120,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C35","83C57"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper shows that RIFT, a grid-based Bayesian algorithm, can recover the masses, spins, and hyperbolic orbital parameters of black-hole scatter and plunge signals at Cosmic Explorer sensitivity when paired with the TEOBResumSDALI…","keywords":["gravitational waves","parameter estimation","hyperbolic black hole encounters","effective-one-body waveforms","RIFT","dynamical capture","Cosmic Explorer","black hole scattering"],"falsifier":"Inject waveforms from an independent numerical-relativity surrogate or Cauchy-characteristic extraction (not TEOBResumSDALI) with known parameters into the same Cosmic Explorer noise, run RIFT, and check whether the true values fall inside the 90% credible regions; a systematic offset would show that the reported accuracy is set by the waveform model rather than by the inference algorithm.","tokens_in":13038,"feed_emoji":"🌊","tokens_out":9377,"duration_ms":100479,"temperature":0.7,"pith_summary":"The paper shows that gravitational-wave parameter estimation can be extended from the familiar quasicircular merging binaries to hyperbolic black-hole encounters—short, one-shot bursts produced when two black holes pass close enough to scatter, briefly capture, or plunge directly into each other. The claim is that RIFT, a grid-based Bayesian inference algorithm, can accurately recover the masses, spins, and the two hyperbolic orbit parameters (system energy $E_0/M$ and angular momentum $p_0^\\phi$ at a fiducial separation) when paired with the effective-one-body waveform model TEOBResumSDALI. The demonstration uses zero-noise scatter and plunge injections at signal-to-noise ratio $\\sim42$ in a single detector with Cosmic Explorer design sensitivity, recovering the injected parameters within the stated credible regions. This matters because unbound black-hole encounters are expected in dense stellar environments and would open a new gravitational-wave channel, but current pipelines are optimized for quasi-circular inspirals.","feed_headline":"RIFT recovers orbit parameters of black-hole flyby signals","feed_subtitle":"Bayesian inference now extracts mass, spin, and orbit parameters from unbound black-hole bursts.","key_machinery":"RIFT (Rapid Iterative FiTting) is a grid-based Bayesian algorithm that first marginalizes over the extrinsic parameters—distance, sky location, inclination, polarization, time, and phase—then iteratively interpolates the marginalized likelihood across the intrinsic parameter grid, using Bayes' theorem to produce posteriors. Its key property here is model-agnostic parallel evaluation, which lets it drive the slower effective-one-body model TEOBResumSDALI. That model generates waveforms for the three unbound-orbit classes—scatter, dynamical capture (zoom-whirl), and direct plunge—from the intrinsic parameters $(m_1,m_2,E_0/M,p_0^\\phi,\\chi_{1,z},\\chi_{2,z})$, with $E_0/M$ the orbital energy and $p_0^\\phi$ the angular momentum at the fiducial initial separation $r_0$. The machinery also includes class-dependent data conditioning: tapering of scatter waveforms (whose pre-event and post-event strain differ) and a peak-finding algorithm to identify the merger time for captures, where the peak strain occurs at the first flyby rather than at merger.","core_discovery":"The central claim is that RIFT—a grid-based, iterative Bayesian parameter-estimation code that marginalizes over extrinsic parameters—can recover the intrinsic parameters of hyperbolic black-hole encounters when paired with the effective-one-body waveform model TEOBResumSDALI. The paper demonstrates this with two zero-noise injections at signal-to-noise ratio $\\sim42$ in a single detector with Cosmic Explorer design sensitivity: a scatter event ($E_0/M=1.01$, $p_0^\\phi=4.40$) and a plunge event ($E_0/M=1.05$, $p_0^\\phi=4.00$), both with equal masses of $20\\,M_\\odot$ and zero spin. For the scatter, the recovered total mass, energy, and angular momentum are tightly localized (the latter two spanning roughly 9% and 1% of their prior ranges); for the plunge, the energy and angular-momentum posteriors are wider and show a sharp likelihood cutoff interpreted as the physical separatrix between plunge and scatter. The paper concludes that parameter estimation of generic hyperbolic waveforms—scatters, dynamical captures, and plunges—is now possible with this infrastructure.","pith_inferences":["Beyond the paper, a direct test of model systematics would be to inject waveforms from an independent numerical-relativity surrogate rather than TEOBResumSDALI itself; the paper's injections and recoveries use the same model, so the reported widths do not include waveform-model error.","The sharp separatrix in the plunge posterior suggests that for real loud events near $p_0^\\phi\\approx4.6$, a hierarchical or mixture treatment across waveform classes could be necessary; the paper does not implement this classification step.","The energy degeneracy in plunges implies that $E_0/M$ information is carried mainly by the pre-merger peak; measuring high-energy plunges may therefore be limited by the low-frequency sensitivity of the detector, a testable prediction for Cosmic Explorer's band.","One could extend the demonstration to a population study: generate a realistic distribution of hyperbolic encounters, run the pipeline, and compare recovered versus injected population hyperparameters to see whether selection effects from the separatrix bias inferred scattering-rate distributions."],"forward_implications":["The same RIFT plus TEOBResumSDALI pipeline can be run on dynamical-capture (zoom-whirl) signals, not only on scatter and plunge events, because the model spans all three waveform classes in one continuous parameter space.","At design sensitivity, a single Cosmic Explorer detector can constrain the total mass to a few percent of the prior range for both scatter and plunge events at SNR $\\sim42$; comparable Advanced LIGO recovery would require SNR $85+$.","The plunge recovery shows a sharp likelihood boundary at $p_0^\\phi\\approx4.6$ that is not a sampling edge but the physical transition from plunge to scatter, so the posterior itself can be used to classify the waveform family.","Because the model also includes tidal deformability, the same infrastructure can be extended to hyperbolic neutron-star or neutron-star–black-hole encounters, which could have electromagnetic counterparts.","Higher-order multipoles, especially the $(2,0)$ mode known to matter in scattering, are expected to break degeneracies in energy and mass ratio that remain in the plunge case."],"supporting_citations":[{"why":"Introduces the TEOBResumSDALI effective-one-body model for dynamical captures and defines the waveform classes and the plunge-scatter separatrix used throughout the paper.","marker":"[31]"},{"why":"Previous parameter-estimation application of the same model to GW190521 as a dynamical capture; it supplies the baseline the paper extends to generic hyperbolic encounters.","marker":"[33]"},{"why":"Supplies high-order analytical information underlying the TEOBResumSDALI model for noncircularized, spin-aligned binaries.","marker":"[52]"},{"why":"Defines the multipolar effective-one-body waveform family for spin-aligned, eccentric, and hyperbolic binaries that TEOBResumSDALI implements.","marker":"[53]"},{"why":"Prior parameter-estimation study with a numerical-relativity surrogate for hyperbolic encounters; it provides the comparison showing why a more flexible infrastructure is needed.","marker":"[60]"},{"why":"Introduces RIFT's grid-based marginalized-likelihood scheme, the inference engine used here.","marker":"[61]"},{"why":"Documents RIFT's rapid inference for precessing binaries, including the aligned-spin prior convention adopted in this work.","marker":"[62]"},{"why":"Describes the highly parallelized RIFT implementation that makes evaluating the expensive hyperbolic waveform model feasible.","marker":"[63]"},{"why":"Provides the Cosmic Explorer design sensitivity curve used for the injections.","marker":"[38]"}],"fun_headline_variants":["RIFT recovers orbit parameters of hyperbolic black-hole flybys","Parameter estimation works for unbound black-hole encounters with RIFT","Black-hole scattering: new method estimates mass, spin, and orbit","Hyperbolic encounters: RIFT tackles generic orbit waveforms for black holes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the same waveform model used to invent the test signals also describes real gravitational waves from hyperbolic encounters, so any systematic error in the model, particularly near the plunge-scatter boundary or from omitting higher modes, would bias all recovered parameters.","fun_headline_variants_meta":{"raw":{"variants":["RIFT recovers orbit parameters of hyperbolic black-hole flybys","Parameter estimation works for unbound black-hole encounters with RIFT","Black-hole scattering: new method estimates mass, spin, and orbit","Hyperbolic encounters: RIFT tackles generic orbit waveforms for black holes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000234,"raw_usage":{"total_tokens":1484,"prompt_tokens":921,"completion_tokens":563,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":537,"completion_tokens_details":{"reasoning_tokens":490}},"tokens_in":537,"tokens_out":563,"duration_ms":6989,"temperature":1.0,"reasoning_tokens":490,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T20:58:30.155034+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Inject waveforms from an independent numerical-relativity surrogate or Cauchy-characteristic extraction (not TEOBResumSDALI) with known parameters into the same Cosmic Explorer noise, run RIFT, and check whether the true values fall inside the 90% credible regions; a systematic offset would show that the reported accuracy is set by the waveform model rather than by the inference algorithm.","supporting_citations":[{"cited_title":"GW190521 as a dynamical capture of two nonspinning black holes","cited_arxiv_id":null,"evidence_quote":"Previous parameter-estimation application of the same model to GW190521 as a dynamical capture; it supplies the baseline the paper extends to generic hyperbolic encounters."},{"cited_title":"Novel scheme for rapid parallel pa- rameter estimation of gravitational waves from compact binary coalescences","cited_arxiv_id":null,"evidence_quote":"Introduces RIFT's grid-based marginalized-likelihood scheme, the inference engine used here."},{"cited_title":"Cosmic Explorer: The U.S. Con- tribution to Gravitational-Wave Astronomy beyond LIGO","cited_arxiv_id":null,"evidence_quote":"Provides the Cosmic Explorer design sensitivity curve used for the injections."}],"review_version":1}