Quarkyonic Matter in Dense QCD

• Quarkyonic matter is a proposed dense QCD phase characterized by a quark-filled Fermi sea and confined baryonic excitations at the Fermi surface.
• Effective models, including mean-field and holographic approaches, reveal a distinct shell structure with pressure scaling as O(Nc)μQ⁴ and a stiff EOS.
• This paradigm impacts neutron star physics by explaining high mass limits and delayed hyperon onset, while suggesting unique heavy-ion collision signatures.

Quarkyonic matter is a theoretically proposed phase of quantum chromodynamics (QCD) expected to arise at high baryon density and low temperature. It is characterized by the coexistence of a dense Fermi sea of quarks—governing the bulk thermodynamics—while excitations near the Fermi surface remain confined into color-singlet baryons. This regime appears in the large number–of–colors (large-$N_c$) limit of QCD and is conjectured to interpolate between ordinary hadronic matter at low density and fully deconfined quark matter at asymptotically high densities. The quarkyonic paradigm underlies a diverse set of recent theoretical frameworks addressing the equation of state (EOS) of dense nuclear matter, the phenomenology of neutron stars, and the structure of the QCD phase diagram (McLerran et al., 2018, Duarte et al., 2023, Torrieri et al., 2012, Fujimoto et al., 2024, Sun et al., 24 Jan 2026).

1. Defining Features and Theoretical Foundations

The concept of quarkyonic matter originated with McLerran and Pisarski, who showed that for large $N_c$, the physics of cold, dense QCD at quark chemical potential $\mu_Q$ satisfying $$\Lambda_\mathrm{QCD} \lesssim \mu_Q \lesssim \sqrt{N_c}\,\Lambda_\mathrm{QCD}$$ is distinct from both ordinary hadronic and deconfined quark-gluon plasma (QGP) phases (Torrieri et al., 2012, McLerran et al., 2018, Kojo, 2010). In this window:

• Bulk degrees of freedom: Deep inside the Fermi sea, thermodynamics is quark-like, with pressure $P\sim O(N_c)\,\mu_Q^4$ and baryon number density $n_B\sim O(N_c)\,\mu_Q^3$.
• Confinement: Near the Fermi surface, confinement persists. The Polyakov loop remains near zero and only color-singlet (baryonic) excitations are accessible.
• Distinct shell structure: The momentum distribution separates into:
  • An interior (low $k$) filled by quarks.
  • A "shell" at the Fermi surface (width $\Delta \sim \Lambda_\mathrm{QCD}$) containing baryons.
  • This structure persists up to a critical chemical potential/deconfinement threshold, at which Debye screening by quarks overcomes gluonic antiscreening and color charges become unscreened (Kojo, 2010, Torrieri et al., 2012, McLerran et al., 2018).
• Phase diagram: The resulting $T$–$\mu_B$ (baryon chemical potential) diagram features a hadronic region at low $\mu_B$, quarkyonic matter at intermediate $\mu_B$, and the deconfined QGP at high $\mu_B$ (Park et al., 2021).

This behavior is robust not only in schematic models but also in explicit treatments, such as mean-field theories, percolation models, and holographic QCD approaches (Duarte et al., 2023, Moss et al., 2024, Kovensky et al., 2020).

2. Microscopic Realizations and Effective Theories

Quarkyonic matter can be constructed in several complementary frameworks:

• Constituent quark models: A Hamiltonian incorporating color confinement and hyperfine (color–spin) interactions reveals that as density increases toward $n_B \sim 4\,n_0$ ($n_0\approx0.16\,\text{fm}^{-3}$), quark modes begin to fill the phase space, particularly freeing the $d$ (or $u$) quark from nucleons with the most attractive $(ud)$ diquark correlations remaining intact (Park et al., 2021).
• Mean-field field theories: Effective Lagrangians with quark, nucleon, and "ghost" (compensating) fields enforce the exclusion of nucleonic states already occupied by quarks in the Fermi sea, ensuring quark–nucleon duality and fully capturing the shell structure of quarkyonic matter (Duarte et al., 2023, Duarte et al., 2021).
• Quarkyonic duality models: Explicit mapping between baryon occupation numbers and quark distributions is imposed, so that quarks fill up to a Fermi momentum $k_F^Q$ and baryons reside in a shell $k_F^Q < k < k_F^B$, with the shell thickness $\Delta$ determined by density and duality conditions. This provides an analytically tractable model to demonstrate the Pauli-blocking constraints intrinsic to quarkyonic matter (Fujimoto et al., 2024, Koch et al., 2024).
• Quantum van der Waals and RMF extensions: RMF or quantum van der Waals models incorporating quark shell degrees of freedom, density-dependent masses, and excluded volume effects yield unified EOSs that interpolate between nuclear and quark–dominated matter at high density, often calibrated to neutron star and heavy-ion constraints (Sun et al., 24 Jan 2026, Moss et al., 2024, Xia et al., 2023).

In all constructions, the dynamical shell structure arises as a consequence of Pauli blocking by the quark Fermi sea, and the shell thickness (parametrically $\Delta \sim \Lambda_\mathrm{QCD}$) sets the scale for baryonic occupation (McLerran et al., 2018, Duarte et al., 2023).

3. Equations of State, Chiral Symmetry, and Phase Structure

A salient property of quarkyonic matter is its impact on the EOS at supranuclear densities:

• Stiffening of the EOS: The emergence of quark degrees of freedom just above saturation density ($n_0$) induces a rapid increase in the pressure $P(n_B)$ and sound speed $c_s^2 = dP/d\epsilon$. Typically, $c_s^2$ exhibits a pronounced peak $c_s^2 \sim 0.6$–1 near $2–4\,n_0$ (McLerran et al., 2018, Moss et al., 2024).
• Crossover behavior: The nucleon→quark transition is typically continuous, with the nucleon shell fraction decreasing smoothly as density increases past a critical threshold. For $n_B \gtrsim 3–4\,n_0$, the bulk is dominated by quarks, with only a thin baryonic shell (Duarte et al., 2023, Gao et al., 2024, Cao et al., 2020).
• Chiral symmetry restoration: As the shell shrinks and baryons fade from the spectrum, chiral condensates decrease and chiral symmetry is restored. In parity–doublet or generalized sigma models, the nucleon and its negative-parity partner become degenerate, and constituent quark masses approach their bare values (Gao et al., 2024).
• Deconfinement: Deconfinement is parametrically delayed to high densities, $$\mu_B^\mathrm{deconf} \sim N_c^{1/2}\Lambda_\mathrm{QCD}$$, since gluon antiscreening dominates over quark screening until very high densities (Torrieri et al., 2012).
• Holographic realizations: The Witten–Sakai–Sugimoto model realizes a chirally symmetric but confined quarkyonic phase at large $N_c$ and high chemical potential, with a first-order transition from baryonic to quarkyonic matter at $T=0$ (Kovensky et al., 2020).
• Chiral spirals: Quarkyonic matter supports spatially inhomogeneous chiral condensates—so-called chiral spirals—breaking translational, chiral, and parity symmetries at the microscopic level, with a large mass gap $\Delta \sim \Lambda_\mathrm{QCD}$ (Kojo, 2010, Kojo, 2011).

4. Phenomenology: Astrophysical and Experimental Consequences

Quarkyonic matter directly addresses several core phenomena in compact stars and heavy-ion physics:

• Neutron star structure: Quarkyonic EOSs accommodate large neutron star masses $M_\mathrm{max}\gtrsim2\,M_\odot$, radii $R_{1.4}\sim11$–13 km, and tidal deformabilities in line with multimessenger observations (GW170817, NICER) (McLerran et al., 2018, Duarte et al., 2023, Moss et al., 2024, Xia et al., 2023, Kumar et al., 2023).
• Hyperon puzzle: Pauli blocking by the quark Fermi sea delays the onset of hyperons to much higher densities ($n_B\gtrsim5–6\,n_0$), and strongly suppresses their low-momentum states. This yields only modest softening of the EOS above hyperon threshold and preserves large neutron star maximum masses (Fujimoto et al., 2024, Sun et al., 24 Jan 2026).
• Cooling and transport: The reduction of proton (and generally lepton) fractions in $\beta$-equilibrium, together with quenching of fast neutrino emission channels, alters predictions for neutron star cooling, generally suppressing direct URCA processes and favoring slower cooling rates (Margueron et al., 2021, McLerran et al., 2018).
• Superconductivity and magnetism: Extensions incorporating color superconductivity (CSQY matter) predict further stiffening of the EOS, superconformal sound speed excursions, and modifications to neutron star mass-radius sequences (Gärtlein et al., 3 Sep 2025). Quarkyonic matter also presents a unique environment for ferromagnetic instabilities, potentially enabling spontaneous large-scale magnetic fields in neutron star cores (Gao et al., 9 Jul 2025).
• Experimental observables: In heavy-ion collisions at FAIR energies, quarkyonic percolation should manifest as a dip in low-mass dilepton production rates—a “band-gap” feature—unlike the continuum of the quark–gluon plasma or the resonance structure of a hadron gas (Torrieri et al., 2012).

Context	Quarkyonic EOS Impact	References
Neutron stars	Large $M_\mathrm{max}$, stiff core	(McLerran et al., 2018, Duarte et al., 2023)
Hyperon puzzle	Delayed onset, mild softening	(Fujimoto et al., 2024, Sun et al., 24 Jan 2026)
Cooling/transport	Suppressed direct URCA, slow cooling	(Margueron et al., 2021, McLerran et al., 2018)
Heavy-ion collisions	Dilepton band-gap signature	(Torrieri et al., 2012)
Magnetism	Possible core ferromagnetism	(Gao et al., 9 Jul 2025)

5. Percolation, Duality, and Model-Independent Features

Beyond model-dependent details, several universal elements of quarkyonic matter emerge:

• Percolation transition: Dense baryonic matter can be recast as a percolation problem for overlapping quark wavefunctions. At high enough baryon density (one baryon per baryon volume), quark wavefunctions percolate throughout the system while color confinement persists locally, giving rise to "percolating but confined" quarkyonic matter (Torrieri et al., 2012).
• Quark-baryon duality: Pauli exclusion enforces that, as low-momentum quark states are filled, baryons are expelled from the core of the Fermi sphere and can only occupy a shell near the Fermi surface (Fujimoto et al., 2024, Duarte et al., 2021, Koch et al., 2024). This duality is enforced in effective theories by ghost fields or direct mapping of occupation numbers.
• Absence of sharp boundaries: The quark-hadron transition in the quarkyonic regime is generally a smooth crossover, with no strong first-order transition—contrasting with MIT bag model–style hybrid stars and enabling wide flexibility to match astrophysical constraints (Duarte et al., 2023, Kumar et al., 2023).
• Compatibility with chiral restoration: Quarkyonic matter naturally incorporates both broken and (partially) restored chiral symmetry, depending on the: chiral condensate, parity doublet splitting, and occupation structure (Gao et al., 2024, Kojo, 2010).

6. Outstanding Issues and Outlook

The quarkyonic framework provides a unified, QCD-grounded paradigm for high-density matter. Yet a number of open problems and developments remain:

• Quantitative calibration: Parameters controlling the shell width, density-dependent masses, and EOS transition need to be continuously tested against ab initio methods (e.g., lattice QCD at finite density, quantum Monte Carlo), nuclear experiment, and heavy-ion data (Duarte et al., 2023, Xia et al., 2023).
• Strangeness and finite-$N_c$ corrections: Realistic models must include hyperon and strange quark degrees of freedom as well as subleading $1/N_c$ corrections (Sun et al., 24 Jan 2026, Fujimoto et al., 2024).
• Nonzero temperature and dynamical transitions: Extending quarkyonic models to $T>0$, isospin asymmetry, and time-dependent processes relevant to mergers: these are being progressively handled in QvdW and RMF-based realizations (Moss et al., 2024, Xia et al., 2023).
• Distinctive observables: Identification of unambiguous astrophysical or experimental signals (e.g., sound speed peaks, magnetic instabilities, modified merger spectra, unique cooling curves, or band-gap dilepton dips) remains a priority for distinguishing quarkyonic matter from other QCD–inspired scenarios (Gärtlein et al., 3 Sep 2025, Torrieri et al., 2012, Kumar et al., 2023, Gao et al., 9 Jul 2025).
• Holographic and lattice QCD links: The emergence of quarkyonic matter in holographic QCD provides a top-down, strongly coupled perspective matching the semiclassical large-$N_c$ picture; continued interplay with nonperturbative lattice results is anticipated (Kovensky et al., 2020).
• Universal phase structure: The persistence of a shell structure and duality between quarks and baryons is robust across a wide range of models and is a plausible candidate for universal behavior in dense QCD.

The quarkyonic paradigm thus provides a coherent and flexible structure for exploring cold, dense strongly-interacting matter, with profound consequences for neutron star astrophysics, high-energy experiments, and QCD phenomenology (Torrieri et al., 2012, McLerran et al., 2018, Fujimoto et al., 2024, Duarte et al., 2023).