{"id":"3a028081-62c2-4ee6-aaff-bcba51d15f76","arxiv_id":"2602.23646","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Across ~3300 consistent finite-temperature neutron-star equations of state, the post-merger gravitational-wave peak frequency is predicted at 2.5–4 kHz (hot median ≈3.0 kHz), favoring a 3 kHz-detuned KAGRA configuration over broadband.","lead":"This paper predicts the gravitational-wave 'ringing' frequency of neutron stars after they collide, using thousands of realistic hot-nuclear-matter models, and finds the peak lies between about 2.5 and 4 kHz. The results recommend where future detectors like a high-frequency KAGRA should aim their sensitivity, and warn against tuning too narrowly.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central thermal-shift and frequency-range claims rest on untested assumption that zero-temperature f-mode relations hold for hot stars; 200–600 Hz shift is smaller than stated 10–30% systematic error.","rationale":"The paper's unique contribution is the use of consistent finite-temperature RMF EoSs; everything quantitative flows through the Doneva et al. relations. The assumption that these cold-calibrated relations are temperature-independent is the sole mechanism by which thermal effects feed into f_peak. It is untested, and the stated 10–30% Cowling error is not propagated, making the thermal shift comparable to the method's noise. This is precisely the reader's weakest_assumption, and I see no reason to move away from a CONDITIONAL verdict. A focused numerical test on the same EoS sets would settle whether the assumption lands.","tokens_in":11015,"tokens_out":7816,"duration_ms":74176,"concrete_test":"Select ~20 EoSs spanning the parameter table (Table 1). For each, compute the S/A=2 TOV configuration and calculate its l=2 f-mode frequency using a relativistic Cowling-approximation code (or a full perturbative code). Compare these frequencies with the predictions of Eq. 4 using the same M0,R0. If the median deviation exceeds ~200 Hz and scales with temperature, the thermal shift claim is undermined. A second check: compare with numerical-relativity f_peak fits (e.g., Bauswein et al. 2013; Radice et al. 2017) for a few EoSs to constrain the systematic error.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claims — the 200–600 Hz thermal reduction of f_peak and the 2.5–4 kHz prediction range — are produced entirely by inserting hot-EoS TOV masses and radii into the Doneva et al. (2013) quasi-universal relations (Eqs. 1–4), which were calibrated on zero-temperature, uniformly rotating neutron stars. The authors assume these relations are temperature-independent, but this is untested. At fixed M0 and R0, a hot star has a different density stratification and effective adiabatic index than a cold star, so the f-mode frequency need not be a unique function of M0 and R0. The cited thermal extensions (Pradhan et al. 2022; Barman et al. 2025) are not used to validate this bridge. The paper itself states in Sec. 4 that the relations are 'calculated using zero-temperature equations of state' and expects 'at least 30% error' from the Cowling approximation. No systematic error is propagated into the stated medians, credible intervals, or SNR ratios. A 30% error at 3.5 kHz is ~1000 Hz, larger than the reported 200–600 Hz thermal shift, so the shift lies within the method's noise. Consequently, the 3 kHz detector recommendation rests on an unvalidated absolute frequency scale.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs ~3300 nucleonic equations of state in a relativistic mean-field model with cold (T=0), warm (S/A=1), and hot (S/A=2) thermodynamics, computes non-rotating TOV masses/radii, and then uses the Doneva et al. (2013) quasi-universal relations (their Eqs. 1–5) to estimate the post-merger gravitational-wave peak frequency f_peak for remnants rotating at the Kepler frequency. Marginalising over EoS and remnant mass, it reports medians (and 90% credible intervals) of 3627+581/−389 Hz (cold), 3476+591/−397 Hz (warm), and 3039+633/−424 Hz (hot), a thermal reduction of roughly 200–600 Hz, and a spread of about 2.5–4 kHz. It then compares KAGRA high-frequency broadband, 2 kHz, and 3 kHz configurations and claims the 3 kHz detuned configuration gives ~2.5 times the signal-to-noise ratio of the broadband configuration. The paper is written as a Letter aimed at informing post-merger detector design.","tokens_in":11272,"tokens_out":9720,"duration_ms":95167,"significance":"If the thermal shift and absolute frequency band were robust, this would be a useful, systematic EoS-ensemble-based estimate for post-merger detector optimisation. The paper has several strengths: it uses a large, internally consistent set of finite-temperature EoSs; the computational pipeline is transparent and the quoted medians are reproducible from the stated equations; the physical direction of the thermal shift (hot stars are less compact, lower f_peak) is sensible; and the authors honestly enumerate their approximations in Sec. 4. The central quantitative claims, however, are not currently supported to the level needed for a detector-tuning recommendation: the frequency band and the thermal shift are generated by inserting hot-EoS masses/radii into zero-temperature quasi-universal fits, and the paper's own stated Cowling-approximation error ('at least 30%') is larger than the reported thermal shift.","major_comments":[{"comment":"The headline 200–600 Hz thermal shift and the 2.5–4 kHz band are generated entirely by substituting hot-EoS TOV masses/radii into the Doneva et al. (2013) quasi-universal relations, which were calibrated on cold, uniformly rotating neutron stars. The paper assumes these relations are temperature-independent, but this is untested: for a hot star the density stratification and effective adiabatic index differ from a cold star of the same (M0,R0), so the f-mode frequency need not be a unique function of (M0,R0). The cited finite-temperature f-mode calculations (Pradhan et al. 2022; Barman et al. 2025) are not used to validate this bridge. Because this assumption carries the central quantitative claims, the authors should either validate the hot-star application against those finite-temperature calculations, use thermal quasi-universal relations derived for hot EoSs, or explicitly re-frame t","section":"Sec. 2.2, Eqs. (1)–(4)"},{"comment":"Systematic uncertainties are not propagated into the quoted medians, credible intervals, or SNR ratios. The paper states that the Cowling approximation incurs 'at least 30% error'; at f≈3–3.5 kHz this is ~1000 Hz, which is larger than the reported 200–600 Hz thermal shift and comparable to the 90% credible intervals. A common systematic offset would not affect relative comparisons between EoSs, but it directly affects the absolute frequency scale used to choose between 2 kHz and 3 kHz detuned configurations. Please provide a conservative systematic error budget (e.g., show how the f_peak distributions and the SNR ratio change under a ±30% frequency scaling) or explain why the detector recommendation is insensitive to such a systematic.","section":"Secs. 2.2 and 3"},{"comment":"The remnant population is constructed by drawing non-rotating M0 uniformly up to M_TOV and then assuming exactly Kepler rotation. Both choices are strong, unphysical priors. The observed binary neutron star mass distribution peaks near 1.33 M⊙, so a uniform prior overweighting high masses inflates the high-frequency tail of the f_peak distribution. Similarly, post-merger remnants are differentially rotating and evolve/spin down; Kepler rotation is an upper limit rather than a representative value. The 2.5–4 kHz band and the relative merit of 2 kHz vs 3 kHz tuning should be tested against a realistic remnant mass prior and a sub-Kepler rotation range (e.g., Ω/Ωk ∈ [0.7, 1.0]) or a differential-rotation parametrisation.","section":"Sec. 2.2, Fig. 3"},{"comment":"The SNR comparison uses a fixed amplitude A=10^-22 and damping time tdamp=0.025 s for all EoSs and remnant masses. The damping time sets the spectral bandwidth and therefore directly affects whether a narrow detuned configuration or a broadband configuration performs better. Real post-merger spectra contain additional peaks (as the paper acknowledges in Sec. 3), and a fixed sine-Gaussian model may bias the optimal tuning. The claim that the 3 kHz configuration gives ~2.5 times the SNR should be shown to be robust to a factor-of-two variation in tdamp and to a multi-peak spectral model (even at lower amplitudes), since the conclusion is otherwise an artifact of the signal model.","section":"Sec. 2.3, Eq. (6)"}],"minor_comments":[{"comment":"Please specify the sign of the azimuthal number m. As written, with l=|m|=2 and m=+2, Eq. (5) gives negative inertial frequencies for the mass/radius range considered; presumably the prograde f-mode requires m=−2. Clarify the convention.","section":"Eq. (5)"},{"comment":"The text says the thermal reduction is '~200−600 Hz', but the median differences are 3627−3476=151 Hz (warm) and 3627−3039=588 Hz (hot). Either quote the full distribution-based range or correct the statement.","section":"Sec. 4"},{"comment":"The abstract's '~2.5 to 4 kHz' range is a generous union of the 90% credible intervals; the actual quoted intervals span roughly 2.6–4.2 kHz. Reconcile the numbers.","section":"Abstract / Sec. 3"},{"comment":"References Abbott et al. 2017b and 2017c both have the same journal/volume/page (PRL 119, 161101); verify the bibliographic data for the 2017c entry.","section":"References"},{"comment":"Typo: 'we make perform semi-analytic calculations' should read 'we perform semi-analytic calculations'.","section":"Sec. 2.2"},{"comment":"The caption says 'cold (S/A=0)', but Sec. 2.1 defines cold as T=0 and the warm/hot configurations as S/A=1 and 2; S/A=0 is not used. Use T=0 for consistency.","section":"Fig. 5 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper is a plausible and clearly written Letter, but its central quantitative claims are suspended on an untested extrapolation of zero-temperature quasi-universal relations to hot, Kepler-rotating remnants. The authors' own admission in Sec. 4 that the relations are zero-temperature and carry at least 30% error should have been accompanied by a systematic propagation and, ideally, a cross-check against the finite-temperature f-mode literature they cite. The additional dependencies on the uniform mass prior and the fixed signal model are also testable. These are substantial but addressable revisions rather than reasons for rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a competent, honest design-study Letter that does something genuinely useful — building a ~3300-member finite-temperature EoS ensemble and pushing it through semi-analytic f-mode fits to produce a predicted post-merger peak-frequency distribution and a KAGRA HF SNR comparison. The paper is transparent about its approximations, which is rare and worth appreciating. But the central quantitative claims — the 200–600 Hz thermal shift and the 2.5–4 kHz band — sit inside the paper's own stated systematic errors. The Doneva et al. relations are calibrated on cold, uniformly rotating stars; the paper inserts hot TOV masses and radii into them and assumes the relations are temperature-independent. That assumption is untested, and 30% of 3.5 kHz is ~1 kHz, much larger than the claimed thermal shift. So the thermal-shift detection is not yet supported.\n\nWhat the paper does well: it uses a consistent thermal EoS framework, applies astrophysical constraints to the ensemble, and is explicit about the limitations (Sec. 4 lists the Cowling error, the zero-temperature relations, the nucleonic-only assumption, fixed S/A, uniform rotation, fixed signal amplitude/damping). The medians and credible intervals are internally consistent with the equations. The signal model is simple but adequate for a design study.\n\nThe width of the distribution is also partly set by a uniform prior on non-rotating remnant mass, which is an arbitrary choice. The SNR ratio comparison (3 kHz detuned ≈2.5x broadband, 2 kHz ≈1.8–2x) is more robust to common-mode errors in the frequency map, but the recommendation of 3 kHz specifically depends on the absolute frequency scale, which carries roughly kHz-level uncertainty. So the direction of the recommendation (use a detuned high-frequency configuration rather than broadband) is probably right; the specific detune frequency is not well pinned down by this method.\n\nIf I were refereeing, I'd ask for: propagation of the 10–30% systematic error into the frequency distributions and SNR ratios; a sensitivity test against temperature-dependent f-mode calculations (or at least a check against numerical-relativity fits); and a discussion of how the mass prior affects the width. None of these are fatal — the paper is a legitimate contribution to detector planning — but they are needed before the numbers are used to commit to a narrowband configuration.\n\nReader's take aligns with mine; the stress-test concern lands squarely. This paper deserves a serious referee; I'd send it to review with requests for revision rather than desk-reject, and I'd cite it for the ensemble construction and the SNR comparison while being careful about the absolute frequencies.","headline":"A transparent, competent design-study Letter whose headline numbers (thermal shift, 3 kHz recommendation) are smaller than the paper's own stated systematic errors; the SNR-ratio comparison is more robust, but the absolute frequency scale is unvalidated.","tokens_in":11979,"tokens_out":2076,"would_cite":true,"duration_ms":21533,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Using realistic thermal equations of state, this paper predicts post-merger gravitational-wave peak frequencies of roughly 2.5–4 kHz, with finite-temperature effects lowering the peak by 200–600 Hz, favouring detectors tuned near 3 kHz.","keywords":["neutron star mergers","post-merger gravitational waves","equation of state","finite temperature","f-mode oscillations","KAGRA high frequency","relativistic mean field","detector optimisation"],"falsifier":"Compare the same ~3300 equations of state in full general-relativistic hydrodynamic merger simulations with consistent finite-temperature microphysics: if the emitted peak frequency falls outside the predicted 2.5–4 kHz range or the thermal downshift differs from 200–600 Hz, the semi-analytic prediction is falsified. A detected post-merger gravitational-wave signal from a binary neutron star merger would settle it observationally.","tokens_in":10744,"feed_emoji":"📡","tokens_out":5670,"duration_ms":50009,"temperature":0.7,"pith_summary":"The paper constructs a large ensemble of neutron-star equations of state from a relativistic mean field model, with consistent finite-temperature treatments at fixed entropy per baryon, and asks where the dominant post-merger gravitational-wave peak should lie. Marginalising over equations of state and progenitor masses, the predicted peak frequency spans about 2.5 to 4 kHz. Including finite-temperature effects, which puff up the remnant and reduce its compactness, lowers the peak frequency by roughly 200–600 Hz. The authors then compare detector configurations and conclude that a high-frequency design detuned near 3 kHz is best matched to these predictions, providing about 2.5 times the signal-to-noise of the broadband configuration. A sympathetic reader would care because this is a concrete, theory-driven input to how next-generation gravitational-wave observatories should be tuned.","feed_headline":"Post-merger GW peaks span 2.5–4 kHz; heat shifts them down 600 Hz","feed_subtitle":"Thermal equations of state point to ~3 kHz tuning, where a detuned high-frequency detector gains ~2.5x signal-to-noise.","key_machinery":"The load-bearing machinery is the combination of (i) a non-linear relativistic mean field (NL-RMF) model that generates zero- and finite-temperature equations of state at fixed entropy per baryon (S/A = 1 and 2) satisfying nuclear and astrophysical constraints; (ii) the quasi-universal f-mode relations of Doneva et al. 2013, which map a non-rotating star's mass and radius to the co-rotating l=|m|=2 f-mode frequency and include Kepler-frequency rotation, converting to the inertial frame by subtracting the rotation frequency; and (iii) the normalised signal-to-noise ratio comparing post-merger-optimised (detuned) and broadband detector noise curves. The f-mode relations carry the argument: fin","core_discovery":"The central claim is that realistic finite-temperature equations of state, built within a phenomenological relativistic mean field model constrained by chiral effective field theory and astrophysical observations, shift the expected post-merger gravitational-wave peak frequency downward relative to cold-matter predictions. For roughly 3300 accepted nuclear parameter sets, the median peak frequency falls from 3627 Hz (cold) to 3476 Hz (warm, entropy per baryon 1) and 3039 Hz (hot, entropy per baryon 2), with 90% credible intervals spanning roughly 2.5–4 kHz. Because thermal pressure expands the remnant, compactness decreases and the fundamental quadrupolar f-mode frequency drops. The authors","pith_inferences":["Because the same thermal puffing that lowers the dominant peak should also shift the secondary post-merger spectral peaks seen in simulations, detector designs targeting those features may need comparable re-tuning — an extension the paper does not make.","The temperature-independence assumption on the f-mode relations is testable: running the same ~3300 equations of state through full general-relativistic merger simulations would show whether the 200–600 Hz shift survives in a fully self-consistent treatment.","If hyperonic or quark degrees of freedom appear, they would soften the equation of state and raise f-mode frequencies, potentially pushing some remnants above the paper's 4 kHz upper bound; the authors note this but leave it to future work.","The uniform-rotation assumption may underestimate the maximum remnant mass; adding differential rotation would stretch the high-mass tail of the distribution and broaden the predicted frequency band further."],"forward_implications":["If the predicted 2.5–4 kHz band is correct, observatories need sensitivity across the full kHz range, not a narrow notch at one frequency.","If the ~200–600 Hz thermal downshift is real, post-merger optimised detectors should be detuned near 3 kHz rather than 2 kHz.","A high-frequency configuration detuned to ~3 kHz would give roughly 2.5 times the signal-to-noise of the broadband design for a typical remnant, improving the odds of a first post-merger detection.","The overlap of cold and warm distributions suggests temperature matters less at S/A = 1 but becomes significant at S/A = 2, so constraining the remnant's thermal state is itself informative."],"fun_headline_variants":["Heat shifts neutron-star merger GW peak down ~600 Hz","Realistic thermal EOS lower post-merger GW peak to ~3 kHz","Post-merger GW peaks span 2.5–4 kHz with realistic finite-T EOS","Finite-temperature EOS shift post-merger GW peak to ~3 kHz"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The quantitative predictions assume the empirical mass-radius-to-frequency relations, calibrated on zero-temperature stars with roughly 10–30% error and uniform rotation, hold unchanged for hot stars at fixed entropy per baryon; if temperature alters the relation itself rather than just the star's size, the 200–600 Hz thermal shift and the 2.5–4 kHz band would need revision.","fun_headline_variants_meta":{"raw":{"variants":["Heat shifts neutron-star merger GW peak down ~600 Hz","Realistic thermal EOS lower post-merger GW peak to ~3 kHz","Post-merger GW peaks span 2.5–4 kHz with realistic finite-T EOS","Finite-temperature EOS shift post-merger GW peak to ~3 kHz"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000566,"raw_usage":{"total_tokens":2521,"prompt_tokens":751,"completion_tokens":1770,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":495,"completion_tokens_details":{"reasoning_tokens":1687}},"tokens_in":495,"tokens_out":1770,"duration_ms":12290,"temperature":1.0,"reasoning_tokens":1687,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T20:15:15.046845+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compare the same ~3300 equations of state in full general-relativistic hydrodynamic merger simulations with consistent finite-temperature microphysics: if the emitted peak frequency falls outside the predicted 2.5–4 kHz range or the thermal downshift differs from 200–600 Hz, the semi-analytic prediction is falsified. A detected post-merger gravitational-wave signal from a binary neutron star merger would settle it observationally.","supporting_citations":[],"review_version":1}