⚛️ nuclear theory

Quarkyonic matter suppresses neutron-star $g$ modes and reverses their mass trend

This paper demonstrates that quarkyonic matter suppresses neutron-star $g$ -mode frequencies by nearly an order of magnitude and reverses their mass dependence compared to standard nucleonic models, due to the strong-interaction-maintained chemical equilibrium that eliminates quark thermodynamic freedom and drastically reduces core stratification.

Original authors: Probit J Kalita, Bharat Kumar

Published 2026-09-23

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Probit J Kalita, Bharat Kumar

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Deep in the heart of a neutron star, where matter is crushed to densities that defy everyday intuition, the laws of physics behave in ways that are both strange and precise. These stellar remnants, the collapsed cores of massive stars, are laboratories for the densest matter in the universe. For decades, astronomers have tried to understand what lies inside them. They know that as you go deeper, the pressure rises, and the material becomes stiffer, resisting further compression. But they do not know exactly what that material is made of. Is it a soup of protons and neutrons, or does it break down into its fundamental building blocks, quarks? The answer matters because the composition of the star determines how it vibrates. Just as a bell rings with a specific tone depending on its shape and material, a neutron star oscillates at frequencies that reveal its internal structure. Among these vibrations, a special type called a gravity mode is unique. Unlike other vibrations that depend on how stiff the star is, the gravity mode depends entirely on how the star's composition changes with depth. If the layers of the star are distinct, this mode rings clearly; if the layers are mixed or uniform, the sound fades away.

A team of researchers in India has now used this principle to explore a specific, theoretical form of matter called quarkyonic matter. This idea suggests that inside a neutron star, quarks exist in a sea, but the protons and neutrons are not dissolved into it. Instead, they are pushed into a thin shell, like a crust on a ball of liquid. This arrangement creates a star that is incredibly stiff and can support a very heavy mass, solving a long-standing puzzle about how some neutron stars can be so massive without collapsing. However, the researchers found that this specific arrangement has a dramatic effect on the star's vibrations. By simulating ten different models of these stars, they discovered that the presence of this quarkyonic core essentially silences the gravity mode. The frequency of the vibration drops significantly, and more surprisingly, the way the vibration changes as the star gets heavier is the exact opposite of what happens in ordinary stars.

To reach this conclusion, the scientists built detailed computer models of neutron stars. They started with the known physics of protons and neutrons and then added the rules for quarkyonic matter, where quarks and nucleons are in a state of constant, rapid chemical balance enforced by the strong nuclear force. This balance happens so quickly that during the brief moment a star vibrates, the quarks cannot move independently; they are locked in step with the protons and neutrons. The researchers calculated how the speed of sound travels through this material. In a normal star, the speed of sound changes depending on whether you squeeze the material slowly or quickly, and this difference creates a restoring force that allows the gravity mode to exist. In their quarkyonic models, this difference almost vanished. The quarks acted like a buffer, absorbing the chemical imbalances that would normally create the vibration. As a result, the core of the star became nearly uniform in its composition, losing the stratification needed to support the gravity mode.

The effect was profound. In ordinary neutron star models, the gravity mode vibrates at frequencies between 158 and 522 hertz. In the quarkyonic models, this frequency dropped to a range of 81 to 254 hertz. But the most striking finding was how the frequency behaved as the star's mass increased. In a normal star, as you add more mass, the star becomes more compact and the gravity mode vibrates faster. In the quarkyonic models, the opposite happened: as the star grew heavier, the vibration slowed down. This reversal occurs because the heavy mass pushes the transition to quarkyonic matter outward, leaving the vibration trapped in a thinner and thinner shell of ordinary matter on the outside. The core, now vast and unstratified, moves as a solid block, while the vibration is confined to the narrow layer surrounding it.

The researchers also looked at how these stars would radiate energy as gravitational waves. Because the core of the quarkyonic star does not participate in the vibration, it does not radiate energy efficiently. The damping time, or how long the vibration lasts, becomes much longer, and the signal reaching Earth would be much weaker. This makes detecting such a star difficult, but it also provides a clear signature. If astronomers can measure the vibration frequency of a neutron star and see how it changes as the star's mass changes, they can distinguish between a normal star and one with a quarkyonic core. A normal star will show a rising frequency with mass; a quarkyonic star will show a falling one.

This work does not prove that quarkyonic matter exists in the universe, but it provides a concrete way to test for it. The study rules out the idea that a simple phase transition, where matter abruptly changes from one state to another, would produce the same result. In those scenarios, the vibration frequency would jump up, not down. The specific behavior found here—a drop in frequency and a reversal of the mass trend—is unique to the continuous, shell-like structure of quarkyonic matter. While the calculations rely on computer simulations and make certain simplifying assumptions, such as ignoring the star's rotation and magnetic fields, the core finding is robust. The sign of the trend, the direction in which the frequency moves, remains clear even when accounting for potential errors in the models.

Ultimately, this paper offers a new lens through which to view the densest matter in the cosmos. It suggests that the silence of a gravity mode, or its strange, slowing rhythm, could be the tell-tale sign of a quarkyonic core. As gravitational wave detectors become more sensitive, they may soon be able to listen to the vibrations of these distant stars. If they hear a frequency that drops as the star gets heavier, it would be a powerful confirmation that the interior of a neutron star is not just a soup of particles, but a complex, layered structure where quarks and nucleons dance in a delicate, silent balance.

Technical Summary: Quarkyonic Matter Suppresses Neutron-Star g Modes and Reverses Their Mass Trend

Problem Statement
The bulk properties of neutron stars (mass, radius, tidal deformability) constrain the equation of state (EOS) of dense matter but suffer from degeneracy: interiors of vastly different composition can yield identical macroscopic observables. Oscillation spectra offer a potential resolution, as mode frequencies depend on both the stiffness of matter and its internal layering (composition gradients). Specifically, gravity ( $g$ ) modes are driven by the difference between the equilibrium sound speed ( $c_e$ ) and the adiabatic sound speed ( $c_s$ ), a difference arising from the inability of slow weak interactions to re-equilibrate lepton fractions during an oscillation. While hybrid stars with first-order phase transitions are known to spike $g$ -mode frequencies due to a discontinuous drop in $c_e$ , the behavior of $g$ modes in quarkyonic matter—a phase where quarks are deconfined but confined to a Fermi sea while nucleons are pushed into a momentum shell—remains uncomputed. The central question is how the unique thermodynamic constraints of quarkyonic matter (strong-interaction equilibrium between quarks and nucleons) affect the buoyancy and oscillation spectrum.

Methodology
The authors construct a family of ten equations of state based on the Zhao–Lattimer model, which imposes genuine chemical equilibrium between quarks and nucleons via the strong interaction.

Models: Seven models feature quarkyonic matter with varying symmetry-energy slopes ( $L$ ), shell-transparency scales ( $\Lambda$ ), and transition densities ( $n_t$ ). Three purely nucleonic models (Zhao–Lattimer functionals) with identical low-density parameters serve as controls.
Thermodynamics: A key theoretical contribution is the derivation of the adiabatic sound speed ( $c_s$ ) for quarkyonic matter. The authors argue that because quark–nucleon equilibration is governed by the strong interaction (timescale $\sim 10^{-23}$ s), it is instantaneous relative to $g$ -mode periods (ms). Conversely, lepton fractions change only via weak interactions (seconds to hours). Thus, $c_s$ is calculated at fixed lepton fractions with strong equilibrium re-imposed at every step. This leads to a closed-form expression for the buoyancy-driving term $c_s^2 - c_e^2$ (Eq. 20), showing that quarks do not act as independent thermodynamic degrees of freedom but rather "slave" to nucleons, absorbing chemical imbalances.
Oscillation Solver: The relativistic Cowling approximation is used to solve the $l=2$ oscillation equations. The code integrates the first-order system for radial and tangential displacements. Extensive verification is performed against independent code paths, finite-difference checks, and analytic limits, confirming the formalism's accuracy to $10^{-6}$ in frequency.
Analysis: Frequencies are decomposed into a dynamical frequency factor $(GM/R^3)^{1/2}$ and a dimensionless remainder to separate structural effects (mass-radius relation) from compositional effects.

Key Contributions

Derivation of Sound Speeds: The paper provides the first derivation of $c_s$ and $c_e$ for quarkyonic matter, demonstrating that the strong-interaction constraint forces the quark sea to buffer chemical imbalances, drastically reducing the buoyancy factor.
Buoyancy Collapse: It is shown that in quarkyonic matter, the product $c_e^2 c_s^2$ jumps by orders of magnitude at the transition, causing the buoyancy factor $N^2 \propto (c_e^{-2} - c_s^{-2})$ to collapse by a factor of 9.5 to 19 (and up to 284 in soft cases), rather than spiking as in first-order transitions.
Mode Confinement: The analysis reveals that the $g$ mode is effectively expelled from the quarkyonic core (which becomes nearly unstratified) and confined to the thin nucleonic shell surrounding it.
Decomposition of Trends: By separating structural and compositional terms in the mass derivative of frequencies, the authors isolate the specific signature of quarkyonic matter.

Results

Frequency Suppression: The fundamental $g_1$ mode frequency in quarkyonic stars is significantly lower than in nucleonic controls. Across the parameter space, frequencies drop from a range of 158–522 Hz (nucleonic) to 81–254 Hz (quarkyonic).
Reversal of Mass Trend: A critical finding is the reversal of the mass trend. In nucleonic stars, $f_{g1}$ increases with mass. In quarkyonic stars with early transitions ( $n_t = 0.3$ fm $^{-3}$ ), $f_{g1}$ decreases monotonically with mass. This occurs because the quarkyonic core expands with mass, thinning the buoyant nucleonic shell where the mode resides, an effect that outweighs the increase in local gravity.
Dimensionless Remainders: When frequencies are scaled by the dynamical frequency, the $f$ and $p_1$ modes of quarkyonic stars behave similarly to nucleonic stars (their trends are inherited from the anomalous mass-radius relation where quarkyonic stars expand with mass). However, the dimensionless $g_1$ remainder ( $f_{g1}/f_{dyn}$ ) is constant to within 2–9% for nucleonic sequences but falls by 25–33% for quarkyonic sequences.
Damping and Radiation: Due to the rigid-plug motion of the non-buoyant core, the energy density perturbation $\delta \epsilon$ is driven to zero in the core. Consequently, the gravitational-wave damping times are orders of magnitude longer, and the strain amplitude is significantly lower compared to hadronic controls.
Universal Relations: The $g$ -mode frequency does not collapse onto a universal relation with compactness or tidal deformability, unlike the $f$ and $p_1$ modes. The quarkyonic $g$ -mode deviates from nucleonic universal relations by factors of two, confirming its sensitivity to composition rather than just structure.

Significance and Claims
The paper claims that the $g$ mode is the only oscillation mode capable of distinguishing quarkyonic matter from standard hadronic matter or first-order phase transitions, primarily through the sign of its mass derivative and the suppression of its dimensionless frequency remainder.

Differentiation: A rising $f_{g1}$ with mass indicates a purely hadronic star. A rising $f_{g1}$ with a discontinuous jump at a specific mass indicates a Gibbs mixed-phase hybrid star. A falling $f_{g1}$ with mass (for $M \gtrsim 1.8 M_\odot$ ) is the unique signature of an early-onset quarkyonic transition.
Observational Strategy: The authors emphasize that a single frequency measurement is insufficient due to overlap in absolute values. The diagnostic power lies in measuring the frequency at two different masses (or one frequency with an independent mass measurement) to determine the slope.
Modesty: The authors explicitly state that their frequencies are lower bounds (due to the Cowling approximation, which underestimates $g$ -mode frequencies by $\sim 10\%$ ). They note that while the sign of the mass trend reversal survives a 10% correction, the smoothness of the trend requires full general-relativistic verification. They also caution that finite weak reaction rates and superfluidity could further suppress the mode or alter the damping, but the fundamental mechanism of buoyancy collapse remains robust. The paper concludes that the $g$ mode is the "sharpest statement" of what a $g$ mode measures: composition, not just structure.

More like this