Exploring Radial Oscillations in Slow Stable and Hybrid Neutron Stars

The Gist

The Big Picture: Squeezing Cosmic Sponges

Imagine neutron stars as the universe's most extreme sponges. They are the dead, crushed cores of massive stars, packed so tightly that a single teaspoon of their material would weigh as much as a mountain. Because they are so dense, they are like a laboratory for testing the laws of physics under extreme pressure.

This paper is about how these cosmic sponges bounce back when you squeeze them. Specifically, the authors are studying what happens when a neutron star pulsates (expands and contracts) like a heart beating. They want to know: Does the star bounce back safely, or does it collapse into a black hole?

The Twist: The "Frozen" vs. "Relaxed" State

Most previous studies assumed that when you squeeze a neutron star, the particles inside have enough time to rearrange themselves instantly to find a comfortable, balanced state. The authors call this the "Equilibrium" state. Think of this like a crowd of people in a room who instantly shuffle around to find the most comfortable spots as soon as the room starts shrinking.

However, the authors argue that in reality, the particles might not have time to shuffle. The "weak force" (a fundamental particle interaction) that allows them to change their identity is slow. So, when the star is squeezed quickly, the particles get "frozen" in their current arrangement. The authors call this the "Adiabatic" (or out-of-equilibrium) state.

The Analogy:

• Equilibrium: Imagine a jar of marbles. If you shake the jar slowly, the marbles settle into the tightest, most efficient packing immediately.
• Adiabatic (Frozen): Imagine shaking that jar very fast. The marbles don't have time to settle; they stay jumbled in the positions they were in before you started shaking. This "jumbled" state is actually stiffer and harder to compress than the settled state.

What They Did

The team built computer models of three different types of neutron stars:

1. ZL70: Made entirely of normal nuclear matter (protons and neutrons).
2. Gibbs40: A "hybrid" star where normal matter turns into quark matter (a soup of free-floating quarks) in a sharp, sudden transition.
3. KW55: Another hybrid star where the transition to quark matter is smooth and gradual.

They then simulated squeezing these stars and calculated two things:

1. Sound Speed: How fast a "pulse" of pressure travels through the star.
2. Stability: At what point the star stops bouncing and collapses.

Key Findings

1. The "Frozen" State is Smoother
When the authors looked at the "frozen" (adiabatic) state, they found that the speed of sound and the stiffness of the star changed more smoothly. In the "relaxed" (equilibrium) state, the transition to quark matter caused jagged spikes and sudden jumps in the data. In the "frozen" state, these jumps were smoothed out.

• Analogy: It's like driving over a bumpy road. In the equilibrium model, you hit a sudden, sharp pothole. In the adiabatic model, it's more like a gentle, rolling hill.

2. The "Slow Stable" Zone
This is the most exciting discovery. Usually, if a neutron star gets too heavy, it becomes unstable and collapses.

• The Old View: Once the star hits its maximum weight, it's done for.
• The New View: Because the "frozen" state is stiffer, the star can actually support more weight before collapsing.

The authors found a "Slow Stable" branch. Imagine a bridge that looks like it should collapse under a heavy truck. In the old model, it falls. In this new model, because the materials inside are "frozen" and stiff, the bridge holds up for a little longer, carrying a heavier load than expected.

3. Connecting to Real Stars (PSR J0740+6620)
There is a real neutron star called PSR J0740+6620 that is incredibly heavy (about 2 times the mass of our Sun) but surprisingly small (a radius of less than 11 km).

• The authors suggest that this star might be sitting on this new "Slow Stable" branch.
• If a star is this heavy but this small, it might be because its internal particles are "frozen" in a stiff configuration, allowing it to exist in a state that was previously thought to be unstable.

The Bottom Line

This paper suggests that we might have been underestimating how heavy and compact neutron stars can be. By accounting for the fact that particles inside these stars can't rearrange themselves instantly (the "frozen" effect), the stars become stiffer. This allows them to survive at higher masses and smaller sizes than previously thought, potentially explaining the existence of heavy, compact stars like PSR J0740+6620.

In short: Neutron stars might be tougher and more flexible than we thought, provided their insides are "frozen" in place.

Technical Summary

Technical Summary: Exploring Radial Oscillations in Slow Stable and Hybrid Neutron Stars

Problem Statement
Neutron stars (NSs) serve as unique laboratories for nuclear physics, where their internal structure and stability are governed by the Equation of State (EoS) of dense matter. While radial oscillations (quasi-normal modes) are critical for probing the EoS and stability of NSs, traditional analyses often assume fluid elements remain in chemical equilibrium (beta-equilibrium) during oscillations. This assumption holds when the relaxation time for weak interactions is significantly shorter than the oscillation period. However, at lower temperatures (below 1 MeV), neutrino relaxation times can exceed oscillation times, causing the chemical composition to remain "frozen" during compression and expansion. This non-equilibrium (adiabatic) scenario alters the adiabatic index and sound speed, potentially extending the range of stable stellar configurations. Furthermore, the presence of phase transitions (hadronic to quark matter) introduces complexities in sound speed and stability that differ between equilibrium and non-equilibrium treatments. This study addresses the gap in understanding radial oscillations for slow, stable hybrid stars under fixed chemical composition, specifically examining models without density discontinuities.

Methodology
The authors construct families of static neutron stars using three distinct EoSs:

1. ZL70: A purely nucleonic model based on the Zhao-Lattimer (ZL) formalism.
2. Gibbs40: A hybrid model incorporating nucleons, quarks, and leptons with a soft first-order phase transition described via Gibbs construction.
3. KW55: A hybrid model similar to Gibbs40 but utilizing a crossover phase transition description.

The study employs the Tolman-Oppenheimer-Volkoff (TOV) equations to determine hydrostatic equilibrium. For radial oscillations, the authors solve the linearized perturbation equations derived from general relativity (following Chandrasekhar). Crucially, they distinguish between two sound speed calculations:

• Equilibrium Sound Speed ($c_{eq}$): Derived from the total derivative of pressure with respect to baryon density, assuming instantaneous beta-equilibrium.
• Adiabatic Sound Speed ($c_{ad}$): Derived from partial derivatives at fixed particle fractions, representing the "frozen" composition scenario where weak processes are too slow to maintain equilibrium during oscillation.

The radial oscillation equations are cast as a Sturm-Liouville eigenvalue problem to determine the fundamental mode frequency ($f$-mode). Stability is assessed by identifying where the $f$-mode frequency vanishes ($f=0$). The authors compare the stability limits defined by the equilibrium case (dot markers) against the adiabatic case (star markers).

Key Contributions and Results

• Smoothness of Physical Parameters: The study demonstrates that considering out-of-equilibrium (adiabatic) effects prevents the singularities and discontinuities often observed in the equilibrium case, particularly near phase transitions. The adiabatic sound speed ($c_{ad}$) and adiabatic index ($\gamma$) exhibit smoother trends compared to their equilibrium counterparts ($c_{eq}$ and $\gamma_{eq}$).
• Extension of Stable Branches: A primary finding is that the adiabatic treatment extends the stable branches of stellar models beyond the maximum mass configuration defined by the equilibrium case. In the equilibrium scenario, stars with masses exceeding the maximum mass are unstable. However, under the adiabatic assumption, these "slow stable" configurations can sustain higher-order mass doublets.
• Frequency Behavior: The $f$-mode frequency vanishes at lower central energy densities in the equilibrium case compared to the adiabatic case. For the adiabatic case, the frequency persists to higher central densities, indicating a prolonged timescale for instability.
• Model Comparisons:
  • The purely nucleonic model (ZL70) yields the highest maximum mass ($2.103 M_\odot$ in equilibrium) and exhibits the stiffest behavior.
  • Among hybrid models, Gibbs40 is stiffer than KW55, resulting in a slightly higher maximum mass.
  • At a fixed mass of $1.4 M_\odot$, the Gibbs40 model exhibits a higher $f$-mode frequency than KW55 and ZL70.
• Radius Implications: The slow stable branch allows for configurations with masses similar to the maximum mass but with radii approximately 0.6 km smaller. This suggests that the lower bound of the NS radius can be extended in the adiabatic scenario.

Significance and Claims
The paper claims that incorporating out-of-equilibrium effects on chemical composition is essential for a comprehensive understanding of neutron star dynamics, particularly for slow stable hybrid stars. By demonstrating that the adiabatic approach extends the stable branches of stellar models, the authors suggest this mechanism could explain the existence and observation of high-mass neutron stars like PSR J0740+6620. Specifically, a radius measurement of $R < 11$ km for a two-solar-mass neutron star might indicate a star residing on a slow stable branch, whereas standard observations (e.g., NICER) suggest PSR J0740+6620 likely resides on the regular stable branch. The study provides insights into the dynamic variations of fundamental radial modes and offers a theoretical framework that bridges the gap between equilibrium and non-equilibrium descriptions of dense matter, relevant for the era of gravitational wave astronomy. The authors note that while their linearized analysis identifies these slow stable branches, higher-order perturbation terms could potentially render these stars unstable at large oscillation displacements.