Radial oscillations of pulsating neutron stars: The UCIa equation-of-state case

This paper investigates how the $\sigma$-cut stiffening scheme applied to the UCIa equation of state affects the equilibrium properties and radial oscillation frequencies of nonrotating neutron stars, demonstrating that the modified model supports stable configurations up to $\sim 2M_\odot$ while systematically increasing mode frequencies, thereby validating radial spectra as a complementary probe of high-density nuclear matter.

The Gist

Imagine a neutron star as the ultimate cosmic weightlifter. It's a dead star so dense that a single teaspoon of its material would weigh a billion tons on Earth. Inside these stars, matter is squeezed to a point where atoms are crushed together, creating a "soup" of subatomic particles.

The big mystery for physicists is: How stiff is this cosmic soup?

If you push on a sponge, it squishes easily. If you push on a steel beam, it barely moves. Neutron stars are somewhere in between, but we don't know exactly how "squishy" or "stiff" they are at their core. This "stiffness" is called the Equation of State (EoS).

This paper is like a stress test for these cosmic weightlifters. The authors are asking: If we change the rules of how matter behaves at extreme pressures, does the star stay standing, or does it collapse?

Here is a breakdown of their work using simple analogies:

1. The Problem: The "Too Soft" Star

Physicists have a baseline model for how neutron star matter behaves (called the UCIa model). Think of this as a standard recipe for a cake.

• The Issue: Recent observations have found neutron stars that are incredibly heavy (twice the mass of our Sun). The standard recipe might be too "soft" to support that much weight without the star collapsing into a black hole.
• The Fix: The authors tried to "stiffen" the recipe. They added a special ingredient called a $\sigma$-cutoff.
  • Analogy: Imagine the standard recipe uses a soft gelatin. The new ingredient is like adding a rigid steel frame inside the gelatin. It doesn't change the taste at the edges (normal density), but in the very center (high density), it stops the material from squishing down as easily.

2. The Experiment: The "Pulsing" Test

How do you know if your new "stiffened" recipe actually works? You don't just look at the star; you make it pulse.

• The Analogy: Think of a drum. If you hit a drum with a loose skin (soft), it makes a low, dull thud. If you hit a drum with a tight skin (stiff), it makes a high, sharp ping.
• The Science: The authors simulated "radial oscillations." This means they imagined the entire star expanding and contracting like a breathing lung.
  • They calculated the frequency (pitch) of this "breathing."
  • Soft Star: Low pitch, slow breathing.
  • Stiff Star: High pitch, fast breathing.

3. The Results: The "Safety Check"

The authors compared two versions of the star:

1. The Original (UCIa): The standard, softer recipe.
2. The Stiffened (UCIa + fs=0.58): The recipe with the steel frame added.

What they found:

• Holding the Weight: The stiffened version could easily support the heavy 2-solar-mass stars that astronomers have actually observed. The original version was struggling.
• The Pitch: The stiffened stars "sang" at higher frequencies. Just like the tight drum, the stiffer the star, the faster it vibrates.
• Stability: This is the most important part. If a star is too heavy for its own structure, it becomes unstable and collapses. The authors checked the "pitch" of the fundamental vibration.
  • If the pitch is real and positive, the star is stable.
  • If the pitch turns imaginary (mathematically speaking), the star is about to collapse.
  • Result: The stiffened models remained stable even at the heaviest observed masses. They passed the stress test.

4. Why This Matters: The "Cosmic Seismologist"

In the past, scientists mostly looked at how big a star is (Radius) and how heavy it is (Mass) to guess its internal structure. It's like trying to guess what's inside a wrapped gift just by weighing it.

This paper introduces Asteroseismology (star-quakes).

• Analogy: Instead of just weighing the gift, they are listening to the sound it makes when you tap it.
• By listening to the "song" of the star (its vibration frequencies), they can confirm if the internal "recipe" is physically possible.

The Bottom Line

The authors showed that by adding a specific "stiffening" rule to their physics model, they could create neutron stars that:

1. Are heavy enough to match real observations.
2. Are stiff enough to not collapse.
3. Vibrate at specific, predictable frequencies that future telescopes might one day hear.

It's a way of saying: "We have a theory for how the universe's densest matter works. We tweaked the math to make it stronger, and it turns out that stronger version is the one that actually survives in our universe."

Future Outlook:
Currently, our radio telescopes and gravitational wave detectors (like LIGO) aren't sensitive enough to hear these high-pitched "pulses" from neutron stars yet. But the authors are preparing the map for the future. When the next generation of super-sensitive detectors (like the Einstein Telescope) comes online, they will be able to listen to these stars and say, "Ah, that pitch matches the 'Stiffened' model!" confirming exactly what is happening inside these cosmic giants.

Technical Summary

1. Problem Statement

Neutron stars (NS) serve as unique laboratories for studying cold, ultra-dense matter. While the Equation of State (EoS) of matter at nuclear saturation density is relatively well-constrained, the behavior of matter at supranuclear densities (several times the saturation density, $n_0 \approx 0.16 \text{ fm}^{-3}$) remains uncertain.

• The Conflict: Observational data from multimessenger astronomy (NICER X-ray timing and LIGO/Virgo gravitational waves) impose strict constraints. Models must support massive neutron stars ($\sim 2M_\odot$), which requires a "stiff" EoS at high densities, while simultaneously satisfying constraints on radii ($R_{1.4}$) and tidal deformabilities ($\Lambda_{1.4}$) from canonical $1.4M_\odot$ stars, which favor softer EoSs.
• The Gap: Most studies focus on static observables (Mass-Radius relations, Tidal Deformability). There is a need to test whether EoS modifications designed to satisfy static constraints also yield dynamically stable configurations under radial oscillations. Radial oscillations provide a sharp diagnostic for dynamical stability in General Relativity (GR).

2. Methodology

The authors employ a multi-step theoretical framework combining Relativistic Mean-Field (RMF) theory with General Relativistic perturbation theory.

A. Microphysical Model (RMF with $\sigma$-cutoff)

• Baseline: The study uses the UCIa parameter set, a calibrated RMF model for pure nucleonic matter (neutrons and protons) in $\beta$-equilibrium.
• Modification: To stiffen the EoS at high densities without spoiling the calibration near saturation, the authors implement a $\sigma$-cutoff scheme.
  • A regulator potential $U_{cut}(\sigma)$ is added to the scalar sector of the Lagrangian.
  • This potential suppresses the growth of the attractive scalar $\sigma$ field at supranuclear densities.
  • Physical Effect: Suppressing the scalar attraction effectively increases the effective baryon masses and raises the pressure at high densities, stiffening the EoS.
• Parameters: Two cases are compared:
  1. $f_s = 0$: The original UCIa (no cutoff).
  2. $f_s = 0.58$: The stiffened UCIa with the regulator active.

B. Macroscopic Structure (TOV Equations)

• The authors solve the Tolman-Oppenheimer-Volkoff (TOV) equations to generate equilibrium sequences for non-rotating, spherically symmetric stars.
• They compute Mass-Radius ($M-R$) relations, maximum masses ($M_{max}$), and tidal deformabilities ($\Lambda$) for both EoS variants.

C. Dynamical Stability (Radial Oscillations)

• The study derives and solves the linearized general-relativistic radial pulsation equations (Sturm-Liouville eigenvalue problem).
• Variables: Radial displacement $\xi(r)$ and pressure perturbation $\eta(r)$.
• Stability Criterion: The stability of a configuration is determined by the sign of the squared eigenfrequency of the fundamental mode ($\omega_0^2$).
  • $\omega_0^2 > 0$: Stable.
  • $\omega_0^2 = 0$: Marginally stable (onset of instability).
  • $\omega_0^2 < 0$: Unstable (collapse).
• Numerical Method: A shooting method is used to find discrete eigenfrequencies ($\omega_n$) and eigenfunctions for the fundamental mode ($n=0$) and overtones ($n=1$ to $10$).

3. Key Contributions

1. Self-Contained Derivation: The paper provides a transparent derivation of the modified RMF field equations in the presence of the $U_{cut}(\sigma)$ potential, clarifying how the scalar response is altered at high density.
2. Dynamical Consistency Check: It moves beyond static constraints by demonstrating that radial oscillation spectra serve as an independent filter for EoS models. A model satisfying $2M_\odot$ constraints must also remain radially stable up to that mass.
3. Systematic Frequency Analysis: The authors compute the full spectrum of radial modes (fundamental + 10 overtones) and analyze the "large frequency separation" ($\Delta \nu$), a key observable in asteroseismology.

4. Key Results

A. Equation of State and Static Properties

• Pressure Stiffening: The stiffened model ($f_s=0.58$) matches the original UCIa up to $n \approx 2.65 n_0$. Beyond this density, the pressure increases significantly due to the suppression of the scalar field.
• Mass-Radius Relations:
  • The stiffened model supports a higher maximum mass ($M_{max}$) compared to the baseline.
  • Both models are consistent with NICER radius constraints for PSR J0030+0451, PSR J0740+6620, and PSR J0437-4715.
  • Both models satisfy the $2M_\odot$ requirement and the tidal deformability constraints from GW170817.
• Tidal Deformability: The stiffening shifts $\Lambda_{1.4}$ to lower values, remaining within the bounds inferred from gravitational wave observations.

B. Radial Oscillation Spectra

• Frequency Shifts: For a fixed stellar mass, the stiffened model ($f_s=0.58$) exhibits systematically higher eigenfrequencies than the baseline model. This reflects the increased pressure support and stiffness of the fluid.
• Stability Limit: Using the $\omega_0^2 = 0$ criterion, the authors identify that the stiffened models remain radially stable up to the observed $\sim 2M_\odot$ mass scale.
• Mode Systematics:
  • The eigenfunctions ($\xi, \eta$) display standard nodal structures (number of nodes equals the overtone number $n$).
  • The large frequency separation ($\Delta \nu_n = \nu_{n+1} - \nu_n$) tends toward a constant value at high overtones, consistent with asymptotic expectations in asteroseismology.

C. Numerical Data (Table I)

• Frequencies for the fundamental mode ($n=0$) range from $\sim 2.3$ kHz to $3.2$ kHz depending on mass and EoS.
• The stiffened model yields higher frequencies (e.g., at $1.4M_\odot$, $\nu_0 \approx 3.20$ kHz for $f_s=0.58$ vs. $3.13$ kHz for $f_s=0$).

5. Significance and Conclusion

• Complementarity: The study demonstrates that radial oscillation spectra provide a complementary consistency test to static and tidal observables. An EoS that is "stiff enough" to support $2M_\odot$ stars must also be dynamically stable; the radial spectrum verifies this link between microphysics (scalar sector regulation) and macroscopic stability.
• Future Prospects: While current ground-based detectors (LIGO) lack the sensitivity to detect these high-frequency ($\sim$ kHz) radial modes, future third-generation detectors (Einstein Telescope, Cosmic Explorer) are expected to reach the necessary sensitivity.
• Implications: If future observations detect radial oscillations, the measured frequency spectrum could directly constrain the high-density behavior of the EoS and the validity of scalar-field regulators like the $\sigma$-cut scheme.

In summary, this paper validates the UCIa + $\sigma$-cutoff model as a viable description of dense matter that satisfies both static multimessenger constraints and dynamical stability requirements, establishing radial oscillations as a critical theoretical tool for future neutron star asteroseismology.