Pith. sign in

REVIEW 4 major objections 6 minor 2 references

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.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review

2026-08-02 20:15 UTC pith:3WIQ7NV4

load-bearing objection 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. the 4 major comments →

arxiv 2602.23646 v1 pith:3WIQ7NV4 submitted 2026-02-27 astro-ph.HE gr-qcnucl-th

Realistic Equations of State Informing Neutron Star Post-Merger Gravitational-Wave Frequencies

Spencer J. Magnall , Nilaksha Barman , Debarati Chatterjee , Paul D. Lasky , Simon Goode This is my paper
classification astro-ph.HE gr-qcnucl-th
keywords neutron star mergerspost-merger gravitational wavesequation of statefinite temperaturef-mode oscillationsKAGRA high frequencyrelativistic mean fielddetector optimisation
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

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.

Core claim

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

What carries the argument

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

Load-bearing premise

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.

What would settle it

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.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

Share X Bluesky LinkedIn Reddit HN

If this is right

  • 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.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • 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.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

4 major / 6 minor

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.

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 (4)
  1. [Sec. 2.2, Eqs. (1)–(4)] 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
  2. [Secs. 2.2 and 3] 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.
  3. [Sec. 2.2, Fig. 3] 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.
  4. [Sec. 2.3, Eq. (6)] 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.
minor comments (6)
  1. [Eq. (5)] 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.
  2. [Sec. 4] 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.
  3. [Abstract / Sec. 3] 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.
  4. [References] 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.
  5. [Sec. 2.2] Typo: 'we make perform semi-analytic calculations' should read 'we perform semi-analytic calculations'.
  6. [Fig. 5 caption] 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.

Circularity Check

0 steps flagged

No significant circularity: f_peak and SNR claims are model predictions from externally calibrated quasi-universal relations applied to independent TOV solutions, not fits relabeled as predictions.

full rationale

The derivation chain is not circular. The paper computes f_peak by solving the TOV equations with its RMF-based hot EoSs, then inserting the resulting (M0, R0) into the Doneva et al. (2013) quasi-universal relations (Eqs. 1–4). Those relations are empirical fits calibrated on zero-temperature stellar models, not on the paper's own output frequencies; using an external fitted relation is a standard modeling step, not a self-definitional reduction. The thermal shift (200–600 Hz) arises because hot EoSs change M0 and R0, which then change the fitted frequency prediction; the paper does not fit any f_peak value to the data and then call it a prediction. The KAGRA SNR comparison is likewise a forward calculation: signals are injected at the model frequencies and evaluated against detector noise curves, so the 3 kHz recommendation follows from the model distribution rather than presupposing it. The self-citations to Barman et al. (2025) and Barman & Chatterjee (2025) provide the thermal EoS construction formalism, but they do not import the f_peak result or any uniqueness theorem; they are independent support for the EoS model, which is constrained by χEFT and astrophysical data rather than by post-merger frequencies. The paper openly flags the main validity caveats—"The Doneva et al. relations also assume uniform rotation and are calculated using zero-temperature equations of state" and "errors of approximately 10−30%" from the Cowling approximation—but these are accuracy and robustness limitations, not evidence that any output is equivalent by construction to its input. The central claims therefore have genuine independent content, even though they carry acknowledged systematic uncertainty.

Axiom & Free-Parameter Ledger

7 free parameters · 8 axioms · 0 invented entities

The paper contributes a large EoS ensemble and a detector-design comparison; it does not contribute a new derivation of the frequency-structure relation. Every absolute frequency in the paper is produced by empirical fits (Eqs. 1–4) applied to TOV-derived M0, R0. The thermal effect enters only through the EoS-dependent M0, R0 variation, not through a temperature-dependent oscillation calculation. The remnant mass prior and fixed signal parameters are hand-chosen. No new physical entities are postulated.

free parameters (7)
  • σ0 universal-relation coefficients (1.562, 1.151) = 1.562, 1.151
    Eq. (4): the non-rotating f-mode frequency scale σ0 [kHz] = 1.562 + 1.151 (M0/1.4M☉)^{1/2}(R0/10km)^{-3/2}. Coefficients are empirical fits from Doneva et al. 2013; every f_peak value in the paper inherits this fitted scale.
  • Kepler-frequency relation coefficients (1.716, −0.189) = 1.716, −0.189
    Eq. (2): ΩK from Doneva et al. 2013 fit; sets the rotation rate at which all remnants are assumed to spin.
  • Corotating-mode rotation coefficients (0.235, 0.358) = 0.235, 0.358
    Eq. (1): correction to σ_corot/σ0 at Kepler rotation; fitted to rotating-star oscillation data.
  • Mass-increase fit coefficients (0.991, 9.36e-3, 3.28) = 0.991, 9.36×10⁻³, 3.28
    Eq. (3): maps non-rotating mass M0 to Kepler-rotating remnant mass; determines the ~24% mass increase and the 1.24 M_TOV upper endpoint of the remnant mass range.
  • Signal model amplitude and damping time = A = 10⁻²², t_damp = 0.025 s
    Sec. 2.3: chosen by hand for the SNR injections and stated to be EoS-independent, so the SNR ratios reflect only noise-curve differences at the chosen f0, not amplitude physics.
  • Remnant mass prior (uniform in non-rotating M0) = uniform 0 ≤ M0 ≤ M_TOV; remnant masses ≈2.2 M☉ – 1.24 M_TOV
    Sec. 2.2: an assumed population distribution; this prior directly controls the width of the f_peak distributions in Fig. 3 and therefore the 'need broadband' conclusion. Not validated against binary mass-ratio distributions or merger simulations.
  • 'Fixed' nuclear parameter set = nsat=0.15 fm⁻³, Esat=−16 MeV, Ksat=240, Jsym=32, Lsym=60, m*/m=0.65
    Table 1: hand-picked central values used for the thermal-profile illustration (Fig. 1) and to isolate thermal effects.
axioms (8)
  • domain assumption Quasi-universal relations (Eqs. 1–4), calibrated on cold rotating/non-rotating neutron star models, hold for all ~3300 EoSs in the ensemble including finite-temperature ones.
    Sec. 2.2 uses Eqs. 1–4 for cold, warm, and hot EoSs; Sec. 4 concedes the relations assume uniform rotation and zero temperature.
  • domain assumption The thermal dependence of the f-mode frequency is fully captured by the thermal dependence of M0, R0; the σ0(M,R) and ΩK relations themselves are temperature-independent.
    This transfer is the whole mechanism that produces the 200–600 Hz thermal shift; no hot-star oscillation calculation is performed.
  • domain assumption Post-merger thermal structure is approximated by fixed entropy per baryon S/A = 1 (warm) and 2 (hot), with no entropy gradients and no thermal evolution during emission.
    Sec. 2.1 and Sec. 4; authors note the entropy profile evolves on the GW emission timescale.
  • domain assumption The remnant is a uniformly rotating star at exactly the Kepler (break-up) frequency.
    Sec. 2.2: Ω = ΩK assumed; differential rotation (the dominant support mechanism for hypermassive remnants) is neglected; authors cite Uryū+2017 that uniformity holds only in the core.
  • domain assumption EoSs are purely nucleonic (no hyperons, no quarks).
    Sec. 4: adding hyperons softens the EoS and raises f_peak; hot hyperonic EoSs reach only ~2.2 M☉, making the 1.24 M_TOV remnant masses unattainable under uniform rotation.
  • domain assumption A uniform prior over non-rotating masses (0 to M_TOV) fairly represents the population of post-merger remnant masses.
    Sec. 2.2: no justification from binary population synthesis or merger simulations; the 2.2 M☉ lower bound is chosen, not derived.
  • domain assumption Cowling-approximation f-modes (as encoded in the Doneva et al. relations) are accurate to 10–30% for the quantities used.
    Sec. 2.2 states 10–30%; Sec. 4 states 'at least 30%'; the quoted error is internally inconsistent and is not propagated.
  • standard math The external constraints used to filter EoSs (χEFT at low density; neutron star mass/radius and tidal observations at high density) are accepted as prior literature.
    Sec. 2.1: parameter sets inconsistent with Drischler+2016, Abbott+2017c, Annala+2018, Most+2018, Tong+2020 are discarded.

reviewed 2026-08-02 · how reviews work

0 comments
Cite this review

Pith. "Pith review of Realistic Equations of State Informing Neutron Star Post-Merger Gravitational-Wave Frequencies." pith.science (2026). https://pith.science/paper/3WIQ7NV4

@misc{pith2026260223646,
  author       = {Pith},
  title        = {Pith review of: Realistic Equations of State Informing Neutron Star Post-Merger Gravitational-Wave Frequencies},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/3WIQ7NV4}},
  note         = {Machine review of arXiv:2602.23646}
}
Share X Bluesky LinkedIn Reddit HN
read the original abstract

Binary neutron star mergers are thought to produce hot, rapidly rotating neutron stars with masses that can far exceed their Tolman-Oppenheimer-Volkoff mass. The gravitational-wave emission from such remnants provides a unique opportunity to measure the nuclear equation of state at densities and temperatures not available to terrestrial experiments. Current detector design is informed by gravitational-wave signals from general relativistic hydrodynamics simulations of neutron star mergers, typically with hybrid thermal treatments for the equation of state, where a cold equation of state is modified by adding a thermal component. We use realistic equations of state based on the relativistic mean field model with consistent treatment of thermal effects to compute the distribution of expected peak gravitational-wave frequencies. Marginalising over equation of state and progenitor neutron star masses, we show the peak frequency of emission ranges from $\sim2.5$ to 4 kHz. The width of this distribution suggests the need for broadband observatories with kHz sensitivity, and calls into question some of the so-called post-merger optimised configurations. We show the proposed KAGRA high-frequency design is well-suited to measuring post-merger remnants when compared to the KAGRA broadband design.

Figures

Figures reproduced from arXiv: 2602.23646 by Debarati Chatterjee, Nilaksha Barman, Paul D. Lasky, Simon Goode, Spencer J. Magnall.

Figure 1
Figure 1. Figure 1: Temperature as a function of baryon density for nucle￾onic matter for warm (S/A = 1) and hot (S/A = 2) neutron star configurations. urations. The mass-radius relations for non-rotating neutron stars using our EoS sets are shown in Appendix A. 2.2 fpeak calculations We estimate of the dominant GW emission frequency fpeak by calculating the emission of the co-rotating l = |m| = 2 f-mode. This ideally require… view at source ↗
Figure 2
Figure 2. Figure 2: Peak GW frequency as a function of mass of post-merger remnant at Kepler rotation for cold (T = 0; grey), warm (S/A = 1; green) and hot (S/A = 2; red) EoSs. 2 3 4 5 fpeak (kHz) 0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 p(fpeak) N(T = 0) N(S/A = 1) N(S/A = 2) [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Peak GW frequency distributions with zero and finite temperature EoSs for all post-merger remnant masses. Shown are results for cold (T = 0; grey), warm (S/A = 1; green) and hot (S/A = 2; red) remnants. frequency for a given post-merger mass configuration as f￾mode frequency observed by an inertial observer. It is given by Equation 5. 2.3 Post-Merger Optimised Detector Configurations To quantify the benefi… view at source ↗
Figure 4
Figure 4. Figure 4: Distributions for the peak frequency of a post-merger remnant of a binary neutron star merger for cold (grey), warm (green) and hot (red) EoSs compared to the sensitivity curves of a 20 km Cosmic Explorer (blue), LIGO at A♯ sensitivity (orange), and various configurations of NEMO (blue, green) and a high frequency KAGRA detector (red, brown, purple). through the Ngarrgu Tindebeek / OzSTAR Australian na￾tio… view at source ↗
Figure 5
Figure 5. Figure 5: Violin plots of expected signal-to-noise ratios of post-merger remnants for our EoSs. The signal-to-noise ratios for the post￾merger optimised detectors are normalised against the standard or ‘broadband’ configuration of the detector. The top panel corresponds to a 2 kHz configuration of the KAGRA HF detector whereas the bottom panel corresponds to a 3 kHz configuration of the KAGRA HF detector. Each panel… view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

2 extracted references · 1 linked inside Pith

  1. [1]

    P., et al., 2017a, Classical and Quantum Gravity, 34, 044001 Abbott B

    Abac A., et al., 2025, The Science of the Einstein Telescope (arXiv:2503.12263),https://arxiv.org/abs/2503.12263 Abbott B. P., et al., 2017a, Classical and Quantum Gravity, 34, 044001 Abbott B. P., Abbott R., Abbott T. D., Acernese F., Ackley K., Adams C., Adams T., et al., 2017b, Phys. Rev. Lett., 119, 161101 Abbott B. P., et al., 2017c, Phys. Rev. Lett....

    Pith/arXiv arXiv 2025

  2. [2]

    The black dashed line represents a signal-to-noise ratio of 1, i.e, exactly the same detector response from the broadband and optimised detectors

    equations of state. The black dashed line represents a signal-to-noise ratio of 1, i.e, exactly the same detector response from the broadband and optimised detectors. MNRAS000, 000–000 (0000) L8Magnall et al. the Royal Astronomical Society, 517, 507 Lattimer J. M., 2012, Ann. Rev. Nucl. Part. Sci., 62, 485 Lattimer J. M., 2021, Ann. Rev. Nucl. Part. Sci.,...

    arXiv 2012

This paper was first reviewed by deepseek-v4-flash on August 2, 2026.