On the treatment of thermal effects in the equation of… — Plain-Language Explanation

Imagine two neutron stars—cosmic behemoths so dense that a teaspoon of their material would weigh a billion tons on Earth—spiraling toward each other. They crash, merge, and create a chaotic, super-hot remnant. This event is one of the most violent in the universe, and it sings a song in the fabric of spacetime called a gravitational wave.

Scientists want to listen to this song to understand the "recipe" of neutron stars. But to do that, they need to run computer simulations. The problem is, how you cook the "thermal ingredients" (heat and pressure) in your simulation changes the outcome.

This paper is like a culinary showdown between two different ways of cooking the same dish: The "Full Recipe" (Tabulated EOS) vs. The "Shortcut Recipe" (Hybrid EOS).

The Two Cooking Methods

The Shortcut (Hybrid EOS): This is the method most scientists have used for a long time. Imagine you are making a complex stew. Instead of tracking every single vegetable and spice, you say, "Okay, the cold part is a solid block, and the hot part is just like steam." It's a simplification. It's fast, easy to calculate, and usually good enough.
The Full Recipe (Tabulated EOS): This is the new, more rigorous method. Here, the computer tracks the actual physics of the hot, dense matter. It knows exactly how the "spices" (particles) interact when they are super-heated and squeezed. It's computationally expensive (like cooking a 12-course meal), but it's much more accurate.

What Happened When They Compared Them?

The researchers ran simulations of neutron star crashes using both methods and waited to see what happened. Here are the key findings, translated into everyday terms:

1. The "Survival of the Fittest" (Will it collapse?)

The Shortcut: Sometimes, the "shortcut" recipe makes the remnant star collapse into a black hole too quickly. It's like the stew losing its structure and turning into a puddle because the heat wasn't calculated right.
The Full Recipe: The more accurate simulations showed that the remnant could actually hold itself together longer. The heat provided enough "push" (pressure) to keep the star from collapsing immediately.
The Lesson: If you use the shortcut, you might think a star dies when it actually survives for a while. This matters because a surviving star sings a different song than a collapsing one.

2. The "Cosmic Song" (Gravitational Waves)

When the merged star vibrates, it emits gravitational waves. Think of this like a bell ringing.

Early Ringing: Right after the crash, both methods heard the same high-pitched "ding" (the fundamental frequency). They agreed on the initial sound.
The Late Ringing: As time went on (after 50–100 milliseconds), the songs diverged. The "Full Recipe" stars started singing a lower, deeper note that the "Shortcut" stars missed or got wrong.
Why it matters: Future detectors (like the "Cosmic Explorer") will be sensitive enough to hear these deep, late notes. If we use the shortcut recipe to interpret the data, we might misread the size or composition of the star.

3. The "Boiling Pot" (Convection)

Inside the hot remnant, the material isn't still; it's churning like a pot of boiling water. This is called convection.

The Discovery: Both methods showed that the star churns, but the pattern of the churning was different.
The "Inertial Modes": This churning excites a special type of vibration called "inertial modes." Think of it like swirling a cup of coffee; the swirl creates its own rhythm. The paper confirms that these swirls happen even with the super-accurate "Full Recipe," but they happen differently than the "Shortcut" predicted.
The Twist: These swirls create a unique gravitational wave signal that could be the "smoking gun" for third-generation detectors.

The Big Takeaway

For years, scientists have been using the "Shortcut" because it was the only thing that fit in their computers. This paper says: "The shortcut is okay for a quick taste, but if you want to know the true flavor of the universe, you need the full recipe."

The differences might seem small at first, but in the world of neutron stars, a small difference in heat calculation can mean the difference between a star that lives for 100 milliseconds and one that dies instantly. As our detectors get better, we need to stop guessing with shortcuts and start using the full, complex physics to decode the universe's most violent songs.

In short: The universe is more complex than our simple models. To hear its true song, we need to listen with the most accurate instruments we have.

1. Problem Statement

The fate of the remnant formed after a binary neutron star (BNS) merger—whether it forms a stable neutron star, a hypermassive neutron star (HMNS), or promptly collapses into a black hole (BH)—is heavily dependent on the Equation of State (EOS) of dense nuclear matter. A critical factor in this evolution is the treatment of finite-temperature effects.

Current numerical relativity simulations predominantly use hybrid EOSs, which approximate thermal effects by adding an ideal-gas-like thermal pressure term ( $P_{th} = \kappa_{th}\rho^{\Gamma_{th}}$ ) to a cold, piecewise-polytropic EOS. While computationally efficient, this approach simplifies complex microphysics. Conversely, fully tabulated, finite-temperature EOSs derived from nuclear physics models offer higher accuracy by explicitly including dependencies on temperature, density, and composition.

The central problem addressed is whether the simplified thermal treatment in hybrid EOSs significantly alters the dynamics of the post-merger remnant, the excitation of gravitational wave (GW) modes (specifically inertial modes), and the validity of "quasi-universal relations" used to infer neutron star properties from GW data. Previous studies suggesting the late-time excitation of inertial modes were based solely on hybrid EOSs; this paper investigates if these findings hold when using more consistent tabulated EOSs.

2. Methodology

The authors performed long-term, high-resolution numerical-relativity simulations of equal-mass BNS mergers ( $M_0 = 1.4 M_\odot$ ) using the Einstein Toolkit and the ILLINOISGRMHD code.

EOS Models: Four microphysical EOSs were selected to span a range of stiffness and thermal properties: SLy4, DD2, HShen, and LS220.
Comparative Approach: For each EOS, two representations were simulated:
1. Tabulated (tab): Fully tabulated, finite-temperature EOSs in $\beta$ -equilibrium.
2. Hybrid (hyb): A cold piecewise-polytropic component (7-piece and 10-piece fits) combined with a thermal ideal-gas term ( $\Gamma_{th} = 1.8$ ).
Initial Data: Single-star equilibrium configurations were constructed to ensure consistency, though the authors noted inherent discrepancies between the tabulated and hybrid representations even for cold, spherical stars.
Simulation Duration: Simulations extended up to $\gtrsim 150$ ms post-merger, allowing for the study of late-time dynamics.
Analysis Tools:
- GW Spectral Analysis: Used Prony's method (ESPRIT) to extract dominant oscillation frequencies ( $f_{max}$ , $f_{2,i}$ , $f_2$ ) and inertial modes.
- Convective Stability: Employed the Ledoux criterion and the Solberg-Høiland criterion (combining Brunt-Väisälä frequency $N^2$ and epicyclic frequency $K^2$ ) to analyze convective instabilities.
- Thermal Analysis: Compared thermal energy distributions and temperature definitions, highlighting the ambiguity of defining temperature in hybrid models.

3. Key Contributions

Direct Comparison of EOS Treatments: This is one of the first studies to directly compare long-term BNS merger evolutions using fully tabulated finite-temperature EOSs against their hybrid counterparts for the same set of microphysical models.
Validation of Inertial Modes: The study confirms that the late-time excitation of inertial modes (previously reported only in hybrid simulations) persists in simulations using fully tabulated EOSs, validating their physical relevance for future GW detectors.
Re-evaluation of Quasi-Universal Relations: The authors derive new quasi-universal relations for post-merger GW frequencies and demonstrate that deviations from established relations become significant at late post-merger phases ( $>30-40$ ms), depending heavily on the EOS treatment.
Temperature Definition Critique: The paper provides a rigorous analysis of temperature definitions, showing that the standard hybrid approximation ( $T_{hyb} \propto \epsilon_{th}$ ) significantly underestimates the true temperature in degenerate matter regimes compared to tabulated EOSs.

4. Key Results

A. Remnant Fate and Dynamics

Collapse vs. Stability: Hybrid EOSs tend to produce more compact remnants that collapse to black holes earlier than their tabulated counterparts. For the soft LS220 EOS, the hybrid model collapsed at $\sim 66$ ms, while the tabulated model remained stable for the full simulation duration.
Thermal Pressure: Tabulated EOSs provide more accurate thermal pressure support, preventing premature collapse. The hybrid approximation often underestimates thermal pressure in high-density, degenerate regions.
Density Evolution: Tabulated models generally reach lower maximum densities than hybrid models, indicating different responses to gravitational compression due to thermal effects.

B. Gravitational Wave Signals

Early Post-Merger ( $<10$ ms): The fundamental quadrupolar $f_{2,i}$ mode frequencies are largely consistent between hybrid and tabulated models, though soft EOSs (SLy4, LS220) show larger discrepancies in merger times.
Late Post-Merger ( $>30$ ms): Significant deviations emerge.
- Frequency Shifts: The late-time $f_2$ mode and inertial modes show frequency shifts of up to 390 Hz (e.g., in the HShen EOS) between tabulated and hybrid models.
- Inertial Modes: Both approaches excite inertial modes (frequencies $< f_{2}$ ) driven by convection. However, the timing and amplitude of these modes differ. Tabulated models show distinct convective patterns that persist beyond 100 ms, triggering sustained inertial mode excitation.
Quasi-Universal Relations: While relations for the merger peak frequency ( $f_{max}$ ) remain robust, relations for intermediate and late-time frequencies ( $f_2$ , inertial modes) deviate significantly from established fits when using tabulated EOSs, suggesting that thermal modeling is crucial for late-time GW inference.

C. Convective Instabilities

Stability Criteria: Using the Solberg-Høiland criterion ( $\omega^2 = N^2 + K^2 > 0$ ), the authors found that while the core is generally stable, unstable rings ( $\omega^2 < 0$ ) form in the outer layers of the HMNS.
Convection Patterns: These unstable regions drive convection, which excites inertial modes. The convective patterns differ substantially between tabulated and hybrid models, with hybrid models sometimes showing more apparent convection due to inaccurate thermal stratification.
Detection Potential: The excited inertial modes have frequencies potentially within the sensitivity band of third-generation GW detectors (Einstein Telescope, Cosmic Explorer).

D. Thermal Energy and Temperature

Temperature Underestimation: The study demonstrates that the hybrid definition of temperature ( $T_{hyb} = (\Gamma_{th}-1)\epsilon_{th}$ ) underestimates the actual temperature in degenerate matter (densities $>10^{12}$ g/cm $^3$ ).
Thermal Profiles: Tabulated EOSs show higher average thermal energy in the bulk and different spatial distributions compared to hybrid models, directly influencing the remnant's size and rotational profile.

5. Significance

Observational Implications: The findings suggest that using hybrid EOSs to interpret late-time GW signals from BNS mergers may lead to systematic errors in inferring neutron star properties (radius, tidal deformability) and the EOS stiffness.
Future Detectors: The confirmation of late-time inertial modes in tabulated simulations strengthens the case for targeting these frequencies with next-generation detectors.
Modeling Standards: The paper argues that for accurate long-term simulations of BNS remnants, especially those investigating thermal evolution and stability, fully tabulated finite-temperature EOSs are necessary. Hybrid approximations, while useful for initial explorations, fail to capture critical microphysical interactions that dictate the remnant's fate and GW emission characteristics.
Theoretical Advancement: By linking convective instabilities to inertial mode excitation in a consistent thermal framework, the study bridges the gap between hydrodynamic stability theory and observable GW signatures.

In conclusion, the paper establishes that the treatment of thermal effects is not merely a minor numerical detail but a fundamental driver of post-merger dynamics, significantly affecting the remnant's lifetime, stability, and the gravitational wave spectrum detectable by future observatories.

On the treatment of thermal effects in the equation of state on neutron star merger remnants