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

By employing realistic relativistic mean field equations of state with consistent thermal treatments, this study demonstrates that post-merger gravitational-wave peak frequencies span 2.5 to 4 kHz, thereby highlighting the necessity for broadband observatories with kilohertz sensitivity and validating the KAGRA high-frequency design over its broadband counterpart.

The Gist

Imagine two neutron stars—cosmic corpses so dense that a teaspoon of their material would weigh a billion tons on Earth—spiraling into each other and colliding. When they smash together, they don't just disappear; they often create a new, super-hot, super-spinning "baby" neutron star. This baby is a chaotic, vibrating mess, screaming out ripples in space-time called gravitational waves.

This paper is like a sound engineer's manual for the future. It asks: "If we want to hear this cosmic scream clearly, what kind of microphone (detector) should we build?"

Here is the breakdown of their findings using everyday analogies:

1. The Problem: The "Hot" vs. "Cold" Confusion

For a long time, scientists tried to predict the "pitch" (frequency) of this cosmic scream using simplified math. They treated the newborn neutron star like a cold, rigid rock.

• The Old Way: Imagine a cold, hard steel ball. If you hit it, it rings at a specific, high-pitched note.
• The New Reality: The paper argues that the newborn star isn't a cold rock; it's a giant, super-heated balloon. It's so hot (hotter than the center of the sun) that the heat makes it puff up and become less dense.

The Analogy: Think of a guitar string.

• A cold, tight string (the old model) vibrates at a high, sharp frequency.
• A hot, puffy string (the new model) is looser and thicker. It vibrates at a lower, deeper frequency.

The authors used complex physics to simulate thousands of these "hot balloons" and found that because they are puffed up by heat, they sing at a lower pitch than we previously thought.

2. The Discovery: The "Sweet Spot" is Lower

The team calculated that the "scream" from these collisions will likely happen between 2.5 kHz and 4 kHz (that's 2,500 to 4,000 vibrations per second).

• The Shift: Because of the heat, the "sweet spot" for listening has shifted down.
• The Implication: If we build a detector tuned to listen to the highest possible notes (like 4 kHz), we might miss the loudest part of the song because the star is actually singing a bit lower, around 3 kHz.

3. The Solution: Tuning the Microphone

Gravitational wave detectors (like LIGO or KAGRA) are like giant microphones. Currently, they are "broadband," meaning they can hear a wide range of sounds, but they aren't super-sensitive to any one specific note.

To hear the faint post-merger signal, scientists propose "detuning" the detector. This is like taking a radio and turning the dial to focus only on one specific station, making that station crystal clear while ignoring the static of others.

• The Test: The authors compared different "tunings" for the KAGRA detector (a Japanese gravitational wave observatory).
• The Winner: They found that tuning the detector to listen specifically at 3,000 Hz (3 kHz) is the best strategy.
  • Tuning it to 2 kHz was okay, but not as good.
  • Tuning it to 3 kHz made the signal 2.5 times louder than using the standard, broad setting.

4. Why This Matters

Right now, we haven't heard the "post-merger" scream yet. Our current microphones aren't sensitive enough. This paper tells us exactly how to build the next generation of microphones.

• The Warning: The authors caution that because the "pitch" can vary (from 2.5 to 4 kHz depending on how heavy and hot the star is), we can't tune our detector to a single perfect note. We need a "broadband" approach that is very sensitive across a wide range, or a detector that can be quickly retuned.
• The Verdict: The proposed KAGRA High-Frequency upgrade, specifically tuned around 3 kHz, is the best bet we have to finally hear these cosmic collisions and learn about the strange, hot matter inside them.

Summary in One Sentence

By realizing that newborn neutron stars are hot, puffy, and sing at a lower pitch than we thought, this paper proves that our future gravitational wave detectors need to be tuned to listen specifically around 3,000 Hz to catch the clearest sound of the universe's most violent collisions.

Technical Summary

1. Problem Statement

Binary neutron star (BNS) mergers produce hot, rapidly rotating remnants that emit gravitational waves (GWs). The frequency of the dominant emission mode (the fundamental $f$-mode) is a unique probe of the hot nuclear Equation of State (EoS) at densities and temperatures inaccessible to terrestrial experiments.

However, current detector designs and search strategies rely heavily on simulations using "hybrid" thermal treatments (a cold EoS plus an ideal gas thermal component with a fixed adiabatic index $\Gamma$). This approach is increasingly recognized as insufficient because:

• It fails to capture non-nucleonic degrees of freedom (e.g., hyperons) expected at merger conditions.
• It does not consistently treat thermal effects across the full parameter space of temperature ($T$), baryon density ($n_b$), and charge fraction ($Y_Q$).
• There is a lack of consensus on the expected peak frequency ($f_{peak}$) distribution when realistic finite-temperature EoSs are used, making it difficult to optimize future high-frequency GW detectors (which often use narrowband "detuned" configurations).

2. Methodology

The authors developed a comprehensive framework to predict post-merger GW frequencies using realistic EoSs and semi-analytic approximations.

A. Equation of State Construction

• Model: They utilized a Non-Linear Relativistic Mean Field (NL-RMF) model.
• Constraints: The model parameters were constrained using:
  • Nuclear saturation properties (density, energy, incompressibility, symmetry energy, etc.) within experimental uncertainty ranges.
  • Microscopic calculations from Chiral Effective Field Theory ($\chi$EFT) at low densities ($0.5 - 1.4 n_{sat}$).
  • Astrophysical observations (mass/radius constraints) at high densities ($>3 n_{sat}$).
• Thermodynamic Configurations: Three distinct states were analyzed to simulate the post-merger environment:
  1. Cold: $T = 0$ K.
  2. Warm: Entropy per baryon $S/A = 1$.
  3. Hot: Entropy per baryon $S/A = 2$ (representing temperatures $\sim 20-80$ MeV found in merger remnants).
• Dataset: Approximately 3,300 valid nuclear parameter sets were generated to construct EoSs for all three thermodynamic states.

B. Frequency Calculation ($f_{peak}$)

Instead of computationally expensive full numerical relativity simulations for every EoS, the authors employed semi-analytic relations (following Chakravarti & Andersson 2020; Doneva et al. 2013):

1. Static Properties: Solved the Tolman-Oppenheimer-Volkoff (TOV) equations to determine the mass ($M_0$) and radius ($R_0$) of non-rotating stars for each EoS.
2. Rotation Effects: Assumed the remnant is a hypermassive neutron star rotating uniformly at the Kepler (break-up) frequency ($\Omega_k$).
   • Calculated the co-rotating $f$-mode frequency ($\sigma_{corot}$) using relations dependent on $\Omega/\Omega_k$.
   • Converted to the inertial frame frequency ($\sigma_{inert}$) using $\sigma_{inert} = \sigma_{corot} - m\Omega$.
3. Mass Distribution: Assumed a uniform distribution of remnant masses ranging from $2.2 M_\odot$ up to $1.24 \times M_{TOV}$ (accounting for mass increase due to rotation).

C. Detector Performance Analysis

• Signal Model: Modeled the GW signal as an exponentially damped sinusoid (dominant $l=|m|=2$ $f$-mode) with fixed amplitude ($A=10^{-22}$) and damping time ($t_{damp}=0.025$ s).
• Detectors: Compared Broadband configurations against Post-Merger Optimized (Narrowband) configurations for:
  • KAGRA HF: Proposed high-frequency upgrade with detuned sensitivities at 2 kHz and 3 kHz.
  • NEMO: Neutron Star Extreme Matter Observatory.
  • Cosmic Explorer (CE): 20 km layout.
  • LIGO A$\sharp$: Future upgrade.
• Metric: Calculated the Signal-to-Noise Ratio (SNR) for $\sim 10^5$ injected signals across all sky positions and EoS sets. The improvement was quantified as the ratio of SNR in the optimized detector vs. the broadband detector ($\rho_{norm}$).

3. Key Results

A. Impact of Temperature on Frequency

• Thermal "Puffing": Finite temperature effects introduce thermal pressure, causing the neutron star to expand ("puff out"), thereby decreasing its compactness.
• Frequency Shift: Since $f$-mode frequency scales with the square root of average density, the reduced compactness leads to a systematic decrease in peak frequency as temperature/entropy increases.
  • Cold ($T=0$): Median $f_{peak} \approx 3627$ Hz.
  • Warm ($S/A=1$): Median $f_{peak} \approx 3476$ Hz.
  • Hot ($S/A=2$): Median $f_{peak} \approx 3039$ Hz.
• Overall Distribution: The total expected range of $f_{peak}$ spans approximately 2.5 kHz to 4 kHz.

B. Detector Optimization

• 3 kHz vs. 2 kHz: The 3 kHz detuned configuration of KAGRA HF provides the best match to the predicted frequency distribution, particularly for hot EoSs.
  • It offers an average SNR increase of $\sim 2.5\times$ compared to the broadband configuration.
  • The 2 kHz configuration offers a smaller improvement ($\sim 1.8-2\times$), which diminishes as the EoS becomes colder (since colder stars emit at higher frequencies, closer to 3 kHz).
• Comparison with Other Detectors:
  • NEMO and Cosmic Explorer (20 km) show similar sensitivity in the expected $f_{peak}$ range but are generally outperformed by the 3 kHz KAGRA HF configuration for this specific narrowband target.
  • The broad distribution of frequencies suggests that while narrowband optimization is beneficial, the bandwidth must be wide enough to capture the $\sim 1.5$ kHz spread.

4. Limitations and Approximations

The authors acknowledge several approximations that may introduce errors (estimated at 10–30%):

• Cowling Approximation: Used in frequency relations, neglecting metric perturbations.
• Uniform Rotation: Assumed the remnant rotates uniformly, ignoring differential rotation which is significant in the outer layers of merger remnants.
• Static Thermal Profile: Used a fixed entropy per baryon profile rather than a dynamic, evolving thermal gradient.
• Nucleonic Matter Only: The EoSs are purely nucleonic. The inclusion of hyperons or quark matter (which soften the EoS) would likely increase frequencies further and potentially lower the maximum mass supportable by uniform rotation.

5. Significance and Conclusion

• Detector Design: The study provides critical data for the design of next-generation GW detectors. It argues that high-frequency optimized detectors (specifically at $\sim 3$ kHz) are essential for detecting post-merger remnants.
• Broadband Necessity: The wide spread of predicted frequencies (2.5–4 kHz) cautions against overly narrow bandwidths. Detectors must balance high sensitivity at a specific frequency with sufficient bandwidth to capture the full distribution of possible remnants.
• Physics Constraints: By moving away from hybrid thermal treatments to realistic finite-temperature EoSs, the paper establishes a more robust baseline for interpreting future GW data. It highlights that ignoring thermal effects leads to an overestimation of the peak frequency, potentially misaligning detector tuning.
• Future Work: The authors call for future studies to incorporate differential rotation and non-nucleonic degrees of freedom (hyperons/quarks) to refine these frequency predictions further.

In summary, this paper demonstrates that realistic thermal treatments shift the expected post-merger GW signal to lower frequencies ($\sim 3$ kHz), strongly supporting the case for KAGRA's proposed 3 kHz high-frequency upgrade as a primary instrument for the first detection of post-merger gravitational waves.