⚛️ nuclear theory

Observable Signatures of a Quarkyonic Phase in Neutron Stars

This paper demonstrates through Bayesian inference that quarkyonic matter equations of state are consistent with current neutron star observations and identifies specific mass-radius signatures, such as curve curvature and central sound speed, that can distinguish neutron stars with quarkyonic cores from those without.

Original authors: Probit J Kalita, Tuhin Malik, Tianqi Zhao, Bharat Kumar, James M. Lattimer

Published 2026-08-04

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Probit J Kalita, Tuhin Malik, Tianqi Zhao, Bharat Kumar, James M. Lattimer

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

Imagine the universe as a giant cosmic kitchen, and inside it, there are some very special, incredibly dense ovens called neutron stars. These are the leftover cores of massive stars that have exploded, crushed down so tightly that a single teaspoon of their material would weigh as much as a mountain. For decades, scientists have been trying to figure out exactly what kind of "dough" makes up the inside of these stars. Is it just a soup of protons and neutrons (the building blocks of normal atoms), or does the pressure get so high that it breaks them apart into something even stranger, like a sea of free-floating quarks?

To solve this mystery, physicists look at two main clues: how heavy the star is and how big its radius (size) is. They also look at how "stiff" the material inside is. Think of stiffness like a mattress: a soft mattress sinks easily when you sit on it, while a stiff one pushes back hard. In a neutron star, if the material is too soft, the star collapses; if it's just right, it holds up a massive weight. Recently, we've learned that some neutron stars are surprisingly heavy (over twice the mass of our Sun) but also surprisingly small (less than 13.5 kilometers across). This combination is a puzzle because it suggests the material inside must get incredibly stiff very quickly as you go deeper, but then maybe relax a bit near the surface. The big question is: what causes this weird behavior?

This paper dives into a specific, exotic recipe for that cosmic dough called "quarkyonic matter." It's a bit like a hybrid dessert where you have a shell of normal particles (nucleons) hugging the outside, but the deep inside is filled with a crowd of non-interacting quarks. The authors wanted to know: if this quarkyonic stuff actually exists in real neutron stars, would we be able to spot it? They didn't just guess; they built a massive computer simulation using a method called "Bayesian inference," which is like a super-smart detective that tests millions of different possibilities against real-world data from telescopes and gravitational wave detectors. They found that yes, it is possible to create a model of quarkyonic matter that fits all the current observations perfectly. But more importantly, they discovered a unique "fingerprint" that could tell us if a star has this quarkyonic core.

The team found that if a neutron star has a quarkyonic core, it changes the shape of the relationship between the star's mass and its radius in a very specific way. Imagine drawing a line on a graph where the x-axis is the star's weight and the y-axis is its size. For normal stars, this line curves one way. But for stars with a quarkyonic core, the line gets steeper and curves differently, almost like a rollercoaster that suddenly shoots up. The authors identified three main ways to spot this difference in stars that weigh about 1.4 times the mass of our Sun (which is a standard, common weight for these stars).

First, they looked at the "sound speed" inside the star. In physics, sound speed tells you how fast a vibration travels through a material. In a quarkyonic star, this speed spikes up dramatically in the core, like a car suddenly hitting the gas pedal. Second, they looked at the "curvature" of that mass-radius line. Quarkyonic stars make this line bend in a distinct way that is different from stars made only of normal matter or stars with other types of exotic matter. Finally, they compared the sizes of two different stars: one weighing 2.0 times the Sun's mass and another weighing 1.4 times. In a quarkyonic star, the difference in their sizes is a specific value that helps separate them from other types.

The paper suggests that if future telescopes and gravitational wave detectors can measure the masses and sizes of neutron stars with high enough precision, we might be able to use these clues to confirm if quarkyonic matter exists. The authors point out that while we can't directly "hear" the sound speed inside a star yet, we can calculate it if we know the mass and radius well enough. They propose that by plotting the "curvature" of the mass-radius line against the size difference between heavy and medium stars, we could clearly separate the quarkyonic stars from the normal ones. It's like having a secret code: if a star falls into a specific zone on this graph, it's strong evidence that its core is made of this special, hybrid quarkyonic stuff.

However, the authors are careful to note that this isn't a guaranteed proof yet. Other things, like a sudden jump to a different type of quark matter, could potentially mimic some of these signs. But the combination of a steep slope, specific curvature, and that size difference is a very strong hint. The paper concludes that with the next generation of super-accurate detectors, we might only need to measure three unique neutron stars to get these numbers right and finally peek inside the heart of a neutron star to see if it's hiding a quarkyonic secret.

Technical Summary: Observable Signatures of a Quarkyonic Phase in Neutron Stars

Problem Statement
The internal composition of neutron stars (NSs) remains an unresolved mystery, particularly regarding the behavior of matter at extreme densities. Current astrophysical observations present a dual constraint: the existence of massive neutron stars ( $\gtrsim 2 M_\odot$ ) requires a stiff equation of state (EoS) at high densities to support the mass, while the compact radii inferred from the binary neutron star merger GW170817 suggest a relatively soft EoS at lower densities. This implies a rapid increase in the speed of sound ( $c_s$ ) within the core. While the "hyperon puzzle" (the softening of the EoS due to hyperon appearance) has been addressed by various models, the potential existence of a "quarkyonic" phase—a state where confined quarks dominate bulk thermodynamics while nucleons exist as color-singlet quasiparticles near the Fermi surface—remains to be observationally constrained. The primary challenge is determining whether quarkyonic matter can satisfy current observational bounds and, if so, identifying distinct signatures that differentiate quarkyonic cores from purely hadronic ones.

Methodology
The authors employ a Bayesian inference framework to constrain a modified $\beta$ -equilibrated quarkyonic matter model. The model, based on Ref. [13], incorporates protons and leptons, assuming that strongly interacting quarks near the Fermi surface form interacting nucleon triplets, while quarks deep in the Fermi sea remain non-interacting. The model is governed by seven free parameters: the transition density ( $\rho_t$ ), the QCD scale parameter ( $\Lambda_{QCD}$ ), symmetry energy ( $J$ ), slope of symmetry energy ( $L$ ), incompressibility ( $K$ ), binding energy ( $B$ ), and a potential interaction term ( $\gamma_1$ ).

The authors utilize nested sampling (via the UltraNest package) to explore the parameter space, generating EoSs that satisfy three progressively constrained sets of observational and theoretical data:

Set $\alpha$ : Includes pulsar mass/radius measurements from NICER (PSR J0030+0451, PSR J0740+6620), tidal deformability from GW170817, and low-density nuclear theory bounds ( $\chi$ EFT) up to nuclear saturation density.
Set $\beta$ : Adds NICER measurements for PSR J0437-4715 and PSR J0614-3329, with a maximum mass ceiling of $M_{max} \le 2.33 M_\odot$ .
Set $\gamma$ : Similar to Set $\beta$ but with a stricter maximum mass ceiling of $M_{max} \le 2.15 M_\odot$ .

All generated EoSs are required to be causal ( $c_s^2 \le 1$ ) and support a mass of at least $1.97 M_\odot$ . The authors solve the Tolman–Oppenheimer–Volkoff equations to derive mass-radius ( $M-R$ ) and mass-tidal deformability ( $M-\Lambda$ ) relations.

Key Contributions and Results

Viability of Quarkyonic EoS: The Bayesian inference demonstrates that quarkyonic matter equations of state can satisfy all current astrophysical observations (NICER, LIGO/Virgo) and theoretical constraints. The posterior distributions for parameters such as $\rho_t$ and $\Lambda_{QCD}$ are constrained, though $\rho_t$ remains partially unconstrained at high densities, indicating that existing observations are insufficient to fully determine the transition density.
Distinctive Signatures via Mass-Radius Geometry: The study identifies that the presence of a quarkyonic phase induces specific geometric changes in the $M-R$ $M - R$ relation, distinct from the softening caused by hyperons. Specifically:
- Curvature ( $\kappa_R$ ) and Slope ($dR/dM$): While hyperons typically reduce the slope and curvature (softening), quarkyonic matter leads to an enhanced slope and curvature due to rapid EoS stiffening.
- Central Sound Speed ( $c_{sc}^2$ ): The transition to a quarkyonic core causes a sharp rise in the central sound speed.
Bimodality and Separability: The posterior distributions for the central sound speed ( $c_{sc}^2$ $c_{sc}^{2}$ ) and the curvature ( $\kappa_R$ $κ_{R}$ ) of canonical ( $1.4 M_\odot$ $1.4 M_{⊙}$ ) neutron stars exhibit bimodality. The authors segregate the valid EoSs into "hadronic" (no quarkyonic phase) and "quarkyonic" (phase present) populations based on whether the central density exceeds $\rho_t$ $ρ_{t}$ .
- Separation Metrics: The study finds that the populations are separable based on the combination of:
  - The central sound speed ( $c_{sc}^2$ ) and the curvature ( $\kappa_R$ ).
  - The radius difference between a $2.0 M_\odot$ and a $1.4 M_\odot$ star ( $\Delta R = R_{2.0} - R_{1.4}$ ) and the curvature ( $\kappa_R$ ).
- Observational Feasibility: While the central sound speed is a strong discriminator, it is not directly observable. The authors propose that the combination of $\kappa_R$ and $\Delta R$ serves as a practical, observable proxy. The $1\sigma$ regions for hadronic and quarkyonic populations show minimal overlap in the $\kappa_R$ vs. $\Delta R$ plane.

Significance and Claims
The paper claims that observing a neutron star with signatures matching the quarkyonic population (specifically a combination of high curvature and specific radius differences) would provide strong evidence for the existence of a quarkyonic phase or a similar crossover transition in the core. The authors emphasize that while other mechanisms (such as first-order phase transitions to quark matter) could produce similar stiffening, the specific signatures identified are indicative of the quarkyonic mechanism.

Crucially, the authors note that distinguishing between a quarkyonic core and a core with massive quark matter (e.g., NJL models) may not be possible with these specific signatures alone, as NJL models can occupy overlapping regions in the parameter space. However, identifying a star in the quarkyonic region would still confirm the presence of quarks in the core. The study concludes that with the anticipated improvements in mass and radius measurement precision from next-generation gravitational wave detectors and radio telescopes, measuring the curvature and radius differences of just three unique neutron stars could be sufficient to infer the nature of the NS core composition.

More like this