{"id":"9db4fff2-b1ff-4bfc-b7bd-b02a8f335c31","arxiv_id":"1908.11418","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A phase-coherent time-domain model for kilohertz postmerger gravitational waves from neutron star remnants is introduced, with injection studies indicating postmerger detection at SNR around 8.5 and radius constraints near one kilometer.","lead":"Researchers built a new waveform model for the gravitational-wave signal emitted a few milliseconds after two neutron stars collide, covering the hard-to-model kilohertz band. The model is calibrated on numerical relativity simulations and aims to let future detectors measure the properties of ultradense nuclear matter.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The SNR-8.5 detectability threshold and the kilometer-precision RTOV_max claim rest on non-converged NR postmerger waveforms, and the injection tests use the same NR family used to fit the model, so shared NR systematics may inflate the apparent accuracy.","rationale":"The reader's weakest_assumption correctly identifies the non-convergence of the underlying NR postmerger waveforms as the primary risk, and the paper itself is unusually candid about this limitation in Appendix C. My reading agrees that the model is structurally sound, clearly specified, and honestly reported, and that conditional acceptance is appropriate. The load-bearing concern is not an internal inconsistency but an external-validation gap: the model is tuned to a NR dataset whose own resolution errors are comparable to the model's mismatch, and the injection study is performed within that same dataset. This does not invalidate the methodological contribution, but it does mean the specific numbers in the abstract (SNR 8.5, kilometer precision) are self-consistency measures rather than demonstrated accuracies against true signals. A concrete high-resolution or independent-code test would settle whether the concern actually shifts the numbers; without such a test, the conditional verdict should remain. I find no basis for rejection: the fits are tabulated, the validation set is held out, the quasiuniversal relations are cross-checked against SACRA for f2, and the limitations are explicitly flagged. The proposed verification step is feasible with existing data and would directly test the paper's headline claims.","tokens_in":26409,"tokens_out":2467,"duration_ms":28520,"concrete_test":"Re-run the validation-set injection recoveries using the highest-resolution (h = 0.135 km) SLy4 1.30+1.30 waveform from Appendix C as the injected signal, and also a Richardson-extrapolated continuum estimate obtained from the VLR/LR/SR/HR sequence, then recover with NRPM at the same SNRs as in Table II. If the recovered f2 or RTOV_max shifts by more than the reported 90% credible intervals, or if the minimal detectable SNR rises above ~8.5, the headline claims are not robust to NR resolution systematics. In parallel, compute full time-domain mismatches of NRPM against independent SACRA waveforms across the same κT2 range; if those mismatches exceed the CoRe-based values of Fig. 4, the training-set validation is insufficient.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claims — minimal postmerger detection SNR ≈ 8.5 and RTOV_max to kilometer precision — are established by injecting CoRe NR waveforms and recovering them with NRPM, a model fit to the same CoRe database. This makes the test circular with respect to NR systematic error. Appendix C explicitly reports that the postmerger phase is monotonic with grid resolution only for a few milliseconds after merger, that long-term data are not in a convergence regime, and that the amplitude is non-monotonic with resolution. Fig. 4 shows that the mismatch between the model and NR is often comparable to the mismatch between NR waveforms at different resolutions, meaning the model's validation error budget is dominated by the very uncertainty that the injection tests assume away. If the true gravitational-wave signal differs from the non-converged NR waveforms by a systematic shift in f2 or in the decay time, then both the fitted quasiuniversal relations (Table I and Eq. 23) and the injection-recovery posteriors inherit that shift; the reported SNR threshold and the quoted ~1 km RTOV_max precision are therefore not validated against real signals. The prompt-collapse criterion (Eq. 16) is acknowledged as ad hoc for hyperonic or phase-transition EOS, but the more load-bearing issue is that the model's self-consistency with its own training data cannot certify accuracy against nature.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs NRPM, a time-domain analytical model of binary neutron star postmerger gravitational waveforms, calibrated to the CoRe numerical-relativity database and characterized by the parameters (mass ratio, total mass, and tidal polarizability kappa_T^2). NRPM is attached to TEOBResumS inspiral-merger waveforms to form a phase-coherent inspiral-merger-postmerger (IMPM) model. The authors validate NRPM against in-family and a small validation set of NR waveforms, report average mismatches of order 0.3, and perform zero-noise Bayesian injection studies to claim that postmerger signals can be detected at network SNR ~8.5, that prompt collapse versus remnant formation can be distinguished by model selection, that the radius of the maximum-mass nonrotating neutron star can be constrained to about 1 km via a recalibrated quasiuniversal relation Rmax(f2), and that inconsistencies between inspiral and postmerger inferences can indicate EOS softening at densities of ~3-5 rho_0.","tokens_in":26721,"tokens_out":4338,"duration_ms":41389,"significance":"If the central claims hold, this is a valuable contribution: it is the first phase-coherent time-domain inspiral-merger-postmerger model for BNS signals, it extends quasiuniversal relations to amplitudes and times, and it demonstrates a concrete Bayesian framework for postmerger parameter estimation and model selection. The authors are also commendably transparent about the limitations of the underlying NR waveforms, explicitly stating in Appendix C that postmerger phase convergence is not achieved on long timescales and that the amplitude is non-monotonic with resolution. The main quantitative results, however, are calibrated and tested on the same simulation family and therefore measure the model's self-consistency against its training data rather than its accuracy against true gravitational-wave signals; for this reason the stated SNR threshold and kilometer-level radius precision are not yet established as statements about real signals.","major_comments":[{"comment":"The abstract and Section V A claim that template-based detection is possible at SNR ~8.5 and that RTOV_max can be determined to about 1 km. These quantitative claims are supported only by zero-noise injections of CoRe NR waveforms, which are the same simulation family used to calibrate NRPM in Section II B, into a model fitted to those simulations. Appendix C states that the postmerger phase is monotonic with grid resolution only for a few milliseconds after merger and that the long-term data are not in a convergence regime, while Fig. 4 shows that NR resolution mismatches are often comparable to the model-NR mismatches. Under these conditions, the injection-recovery tests validate the model's self-consistency with its training data rather than its accuracy against true signals, and the quoted threshold and precision do not include the dominant systematic error. Please propagate NR resolution uncertainties into the injection studies (for example, by injecting waveforms at different resolutions and adding a resolution-based systematic error term) or explicitly restate the claims as self-consistency statements.","section":"Sec. V A; Appendix C; Fig. 4"},{"comment":"The validation set of Table II is not independent of the calibration data: the fits in Section II B use 148 CoRe simulations plus 24 additional simulations, and the ten validation simulations are drawn from the same CoRe collaboration, often with the same microphysical EOS and evolution code, as the training set. The only external check in Appendix A compares the f2 quasiuniversal relation with SACRA data; it does not validate the full time-domain waveform model. As a result, the average mismatch ~0.3 and the injection results in Table II cannot be interpreted as a test of the model's ability to represent waveforms produced by other codes or EOS models. Please perform a full-waveform validation on independent NR data (for example, the SACRA catalog) or explicitly label the current validation as in-family.","section":"Sec. II B; Table II; Appendix A"},{"comment":"Table II shows that several injections recover biased parameter values at or above the claimed detection threshold: for the SLy4 (1.364+1.364) injection the f2 posterior is bimodal with the dominant mode far from the injected value; for the DD2 (1.50+1.50) injection the recovered kappa_T^2 of 196^{+79}_{-68} is inconsistent with the injected value of 91.1 at more than one sigma; and for the H4 injection Rmax is overestimated. The text attributes these biases to NR fit inaccuracies and model systematics. Because these biases occur in exactly the injections used to support the detectability claim, the statement in Section V A that 'the posterior distributions of the physical parameters include the injected values within the 95% confidence regions' is misleading. Please quantify the fraction of injections that are unbiased and discuss how the SNR threshold and the Rmax precision claim change when these systematic biases are incorporated.","section":"Sec. V A; Table II"},{"comment":"The prompt-collapse criterion kappa_T^2 < 80 +/- 40 in Eq. (16) is an empirical fit to hadronic EOS simulations and is explicitly acknowledged in the text to fail for hyperonic or phase-transition EOS. The paper nonetheless claims in the abstract and in Section V B that the model can infer whether the merger outcome is prompt collapse or a remnant star. The demonstration in Table III uses only the 2B and BHB Lambda-phi injections, both of which are consistent with the calibrated criterion. Since the model-selection argument depends on the template family rather than on Eq. (16) alone, please either demonstrate prompt-collapse inference for a case in which Eq. (16) is violated or restrict the claim to the class of EOS for which the criterion was calibrated.","section":"Sec. II B 2; Sec. V B; Eq. (16)"}],"minor_comments":[{"comment":"The abstract states a 'minimal signal-to-noise ratios (SNR) of 8' while the conclusion states SNR ~8.5; please reconcile these numbers and correct the singular/plural agreement in 'signal-to-noise ratios'.","section":"Abstract and Conclusion"},{"comment":"There are typographical errors in the tables and text: 'mininum' and 'maxinum' in Table I should be 'minimum' and 'maximum', 'mantaining' in Section III should be 'maintaining', and 'quasinuniversal' should be 'quasiuniversal'.","section":"Table I and Table II"},{"comment":"The caption says the primary posterior peak at ~5.2 kHz is 'beyond the Nyquist limit, not in the plot', but the analysis band is stated as [1024, 4096] Hz, so the posterior maximum exceeds the analysis band rather than the Nyquist frequency; please clarify how a marginalized posterior can be produced outside the analyzed frequency band.","section":"Fig. 7 caption"},{"comment":"Equation (23) is called an approximate relation but the reported chi^2 = 7.4e-5 is given without the number of degrees of freedom or the covariance of the fitted coefficients; please report the fit residuals and the coefficient covariance matrix so that the uncertainty on Rmax can be propagated correctly.","section":"Sec. V C; Eq. (23)"},{"comment":"The injections are performed in zero noise, so the quoted detectability threshold does not include false-alarm statistics, noise realization effects, or calibration uncertainties; at minimum this limitation should be stated in Section V A rather than only implicitly in the methodology.","section":"Sec. V A"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know: this is the first time-domain postmerger model that attaches phase-coherently to an EOB inspiral-merger waveform, and it is fitted to a much larger NR catalog than previous models. That is a genuine step forward. The quasiuniversal relations are extended to amplitudes and times, the Rmax(f2) relation is recalibrated, and the Bayesian model-selection demonstration for prompt collapse versus remnant is well constructed. The fits are tabulated, the validation set is held out, and Appendix A compares against independent SACRA data. Credit is due there.\n\nThe soft spots are real but mostly acknowledged in the paper itself. The headline numbers—minimal postmerger detection SNR around 8.5 and RTOV_max to about a kilometer—come from zero-noise injections of CoRe waveforms recovered with a model fitted to the same CoRe database. That is a self-consistency test, not a measurement of accuracy against nature. The stress-test note gets this right, and the paper's own Appendix C supports it: the postmerger phase is monotonic with resolution only for a few milliseconds, long-term data are not in a convergence regime, and Fig. 4 shows model-NR mismatches comparable to NR-NR resolution mismatches. If the true signals differ systematically from the non-converged NR waveforms, both the fitted relations and the injection posteriors inherit that shift. So the SNR threshold and the kilometer precision should be read as conditional on the NR waveforms being representative.\n\nThere are also some known biased recoveries—DD2 f2 overestimated, SLy4 bimodal posteriors—and the prompt-collapse criterion kappa_T^2 < 80 +/- 40 is admittedly ad hoc for hyperonic or phase-transition EOS. The abstract says SNR 8 while the body and conclusion say 8.5; a minor internal inconsistency.\n\nNone of this makes the paper unserious. The authors are unusually transparent about the NR limitations, and the model is clearly specified and reproducible enough to be built upon. The central argument holds as a feasibility study, not as a validated measurement pipeline. This paper deserves a serious referee and likely conditional acceptance, with the quantitative claims softened or explicitly framed as calibration against current NR data. I would bring it to reading group and would cite it as the standard reference for this modeling approach.","headline":"The first phase-coherent time-domain postmerger model is a real advance and the paper is honestly reported, but the SNR-8.5 and kilometer-precision claims are self-consistency checks against non-converged NR waveforms, not validated predictions.","tokens_in":27279,"tokens_out":1611,"would_cite":true,"duration_ms":18073,"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":"NRPM is the first phase-coherent time-domain model spanning inspiral, merger, and kiloHertz postmerger neutron-star waveforms, and it makes postmerger detection feasible at SNR about 8.5.","keywords":["gravitational waves","neutron star mergers","postmerger waveforms","equation of state","numerical relativity","kiloHertz gravitational waves","template-based detection","quasiuniversal relations"],"falsifier":"A concrete test: run the same NRPM fits on a new set of error-controlled, higher-resolution postmerger simulations with microphysical equations of state and check whether the predicted $\\hat{f}_2$ and the $\\hat{R}_{\\max}(\\hat{f}_2)$ relation shift by more than the current fit uncertainties. A real detection with independently measured inspiral radius would also settle it: if the postmerger-inferred $R_{\\rm TOV}^{\\max}$ disagrees with the inspiral-based radius at high SNR, the quasiuniversal calibration is wrong or the postmerger remnant probes different physics.","tokens_in":26191,"feed_emoji":"🌊","tokens_out":7439,"duration_ms":66542,"temperature":0.7,"pith_summary":"This paper constructs the first phase-coherent time-domain model that joins the inspiral, merger, and postmerger stages of a binary neutron star gravitational-wave signal, extending effective-one-body inspiral waveforms into the kiloHertz postmerger regime. The model, called NRPM, is fitted to a large set of numerical-relativity simulations and parameterized by binary mass ratio, total mass, and the tidal polarizability $\\kappa^T_2$. The paper argues that with this model, a template-based search can identify postmerger signals from GW170817-like events at network signal-to-noise ratios around 8.5 using third-generation detectors, and that Bayesian model comparison can tell whether the remnant collapsed promptly to a black hole or formed a neutron star. For non-prompt-collapse remnants, the paper derives a quasiuniversal relation between the dominant postmerger frequency and the radius of the maximum-mass nonrotating neutron star, and claims single-event constraints at the kilometer level. Such measurements would probe matter at densities several times nuclear saturation, where the equation of state is least known.","feed_headline":"SNR 8.5 suffices to detect neutron-star postmerger waves","feed_subtitle":"A single third-generation detection could also measure the maximum-mass neutron-star radius to about a kilometer.","key_machinery":"The central object is NRPM, a time-domain postmerger waveform model built from a C1 piecewise frequency and amplitude construction. The frequency is assembled from three characteristic spectral frequencies, the dominant $\\hat{f}_2$ and sidebands $\\hat{f}_{2\\pm 0}$, with a cubic interpolation from the merger frequency and a sine-modulated oscillation; the amplitude interpolates between four extrema and then decays exponentially. All parameters are fitted to numerical-relativity data through quasiuniversal rational or linear relations in $\\xi=\\kappa^T_2+c(1-4\\nu)$. The quasiuniversal relations are the load-bearing machinery: they convert the model's few physical parameters into concrete waveform morphology, and they are also the instrument by which the paper translates a measured peak frequency into the radius $R_{\\rm TOV}^{\\max}$ and into indications of equation-of-state softening.","core_discovery":"The central claim is that the kiloHertz gravitational-wave signal from a binary neutron star remnant can be captured by an analytical time-domain waveform whose parameters are determined, through quasiuniversal fits, by the binary's total mass, symmetric mass ratio, and tidal polarizability $\\kappa^T_2$. When attached to an effective-one-body inspiral-merger waveform at the amplitude peak, the model keeps phase coherence across the full observed band and enables matched-filtering and Bayesian inference on the postmerger portion. The paper's validation shows mismatches averaging around 0.3 against numerical-relativity waveforms, with better fidelity for long-lived remnants and lower postmerger frequencies. From the recovered postmerger peak frequency, the paper derives a quasiuniversal mapping to $R_{\\rm TOV}^{\\max}$, the radius of the most compact nonrotating neutron star, and states that a single detection at minimal SNR constrains that radius to about one kilometer. It also reports that inconsistencies between inspiral and postmerger tidal-parameter posteriors can signal equation-of-state softening at extreme densities.","pith_inferences":["The quoted SNR threshold is a statement about how well NRPM matches the specific simulated waveforms used, not yet a guarantee about real sky signals; error-converged postmerger templates would be needed to confirm the threshold observationally.","The same peak-frequency to radius machinery could be turned around: if a future high-SNR detection yields a postmerger-inferred radius that disagrees with the inspiral-derived radius, the discrepancy would point to new physics at postmerger densities, though the paper's hyperon example shows that breaking the quasiuniversal relation alone cannot identify the microphysical cause.","Promoting the damping parameter $\\alpha$ to a free inference parameter should remove the distance-estimation bias the paper notes and would make the model more agnostic; this is an immediate, testable extension.","Using independent numerical-relativity catalogs to refit the quasiuniversal relations, as the paper begins to do for $\\hat{f}_2$, would quantify how much of the validation mismatch is model bias versus numerical-relativity uncertainty."],"forward_implications":["A template-based matched-filter search using NRPM can claim detection of postmerger signals at network SNR $\\sim 8.5$, roughly the level expected for GW170817-like events with third-generation detectors.","Bayesian model selection between inspiral-merger and inspiral-merger-postmerger hypotheses distinguishes prompt collapse to a black hole from a surviving neutron-star remnant, with decisive log Bayes factors in the two tested injections.","For remnants that do not promptly collapse, the quasiuniversal relation $\\hat{R}_{\\max}(\\hat{f}_2)$ converts the measured postmerger peak frequency into an estimate of the maximum-mass neutron star radius with uncertainty near one kilometer.","The inspiral-postmerger consistency check can flag equation-of-state softening at densities of about 3 to 5 times nuclear saturation, as in the hyperon case studied, already at SNRs of about 11.","The model's quasiuniversal fits are directly usable by other waveform construction strategies, extending earlier results on spectral peak relations."],"supporting_citations":[{"why":"Supplies the quasiuniversal relation between the postmerger peak frequency and tidal polarizability that NRPM's frequency fits build on.","marker":"[16]"},{"why":"Introduces the prompt-collapse criterion $\\kappa^T_2<80\\pm40$ adopted in the model.","marker":"[12]"},{"why":"Provides the numerical-relativity waveform database used for the fits and model construction.","marker":"[61]"},{"why":"Provides the effective-one-body tidal waveform used for the inspiral part and the hybrid attachment at merger.","marker":"[82]"},{"why":"Introduced the $R_{\\rm TOV}^{\\max}(f_2)$ relation that the paper recalibrates in Eq. (23).","marker":"[48]"},{"why":"Supplies hyperon and DD2 simulations used to demonstrate inference of equation-of-state softening at extreme densities.","marker":"[58]"},{"why":"Characterizes long-lived remnant morphology and the m=1 spiral mode relevant to postmerger modeling.","marker":"[29]"},{"why":"Provides an independent numerical-relativity catalog used to cross-check the quasiuniversal $\\hat{f}_2$ fit.","marker":"[51]"}],"fun_headline_variants":["Model catches postmerger waves at SNR 8, pins neutron-star radius","Postmerger waves at SNR 8 give km-level radius","A single detection at SNR 8 measures neutron-star radius to km","Time-domain model unlocks postmerger waves and stiff-matter tests","Postmerger waveform model offers km-radius constraint at SNR 8"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the numerical-relativity postmerger waveforms used for fitting are accurate enough to anchor the quasiuniversal fits; the paper itself reports those waveforms are not yet in a convergence regime after the first few milliseconds, and it also relies on an ad hoc prompt-collapse criterion, $\\kappa^T_2<80\\pm40$, that may fail for hyperonic or phase-transition equations of state.","fun_headline_variants_meta":{"raw":{"variants":["Model catches postmerger waves at SNR 8, pins neutron-star radius","Postmerger waves at SNR 8 give km-level radius","A single detection at SNR 8 measures neutron-star radius to km","Time-domain model unlocks postmerger waves and stiff-matter tests","Postmerger waveform model offers km-radius constraint at SNR 8"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000833,"raw_usage":{"total_tokens":3652,"prompt_tokens":975,"completion_tokens":2677,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":591,"completion_tokens_details":{"reasoning_tokens":2585}},"tokens_in":591,"tokens_out":2677,"duration_ms":18074,"temperature":1.0,"reasoning_tokens":2585,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:15:54.762677+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete test: run the same NRPM fits on a new set of error-controlled, higher-resolution postmerger simulations with microphysical equations of state and check whether the predicted $\\hat{f}_2$ and the $\\hat{R}_{\\max}(\\hat{f}_2)$ relation shift by more than the current fit uncertainties. A real detection with independently measured inspiral radius would also settle it: if the postmerger-inferred $R_{\\rm TOV}^{\\max}$ disagrees with the inspiral-based radius at high SNR, the quasiuniversal calibration is wrong or the postmerger remnant probes different physics.","supporting_citations":[{"cited_title":"The gamma-rays that accompanied GW170817 and the observational signature of a magnetic jet breaking out of NS merger ejecta","cited_arxiv_id":"1710.05897","evidence_quote":"Introduced the $R_{\\rm TOV}^{\\max}(f_2)$ relation that the paper recalibrates in Eq. (23)."}],"review_version":1}