Beyond Protons and Neutrons: Unveiling the Secrets of Quarkyonic Matter

Author: Denis Avetisyan

New research explores the exotic state of matter where quarks become dominant, bridging the gap between traditional nuclear physics and the realm of quark-gluon plasmas.

In symmetric nuclear matter, quark momentum distributions-calculated using the quarkyonic Quantum Monte Carlo model for a density parameter of <span class="katex-eq" data-katex-display="false"> r_p = 0.7 </span> fm-exhibit a shift at densities exceeding saturation, indicating full occupation of low-momentum quark states and reflecting the underlying many-body physics of nucleon interactions. — In symmetric nuclear matter, quark momentum distributions-calculated using the quarkyonic Quantum Monte Carlo model for a density parameter of $r_p = 0.7$ fm-exhibit a shift at densities exceeding saturation, indicating full occupation of low-momentum quark states and reflecting the underlying many-body physics of nucleon interactions.

A relativistic quark model reveals how nuclear interactions influence the equation of state and saturation density of quarkyonic matter.

The conventional understanding of nuclear matter at extreme densities remains incomplete, particularly concerning the transition between hadronic and quark-dominated phases. This is addressed in ‘A dual description of quarks and baryons: Quarkyonic matter within a relativistic quark model’, which investigates quarkyonic matter-a phase where both quarks and baryons coexist-using a relativistic quark model combined with the quark-meson coupling model. The study reveals that nuclear interactions significantly influence the equation of state, leading to an earlier onset of quark saturation and a stiffer high-density behavior compared to non-interacting models. How do these findings refine our understanding of neutron star interiors and the properties of matter under the most extreme conditions in the universe?

Unveiling the Density Puzzle: Beyond Conventional Models

The intense gravitational forces within neutron stars compress matter to densities far exceeding that of atomic nuclei, presenting a significant challenge to conventional hadronic models. These models, which describe matter as composed of protons and neutrons interacting through established forces, begin to falter under such extreme conditions. Predictions derived from these models regarding the star’s equation of state – the relationship between pressure and density – diverge significantly from observational constraints gleaned from gravitational wave detections and electromagnetic observations of neutron star mergers. This discrepancy suggests that the underlying assumptions of hadronic models, particularly regarding the behavior of strongly interacting particles at ultra-high densities, require refinement or even a fundamental reassessment to accurately capture the physics at play within these cosmic laboratories.

Current theoretical models of matter at extreme densities, while successful in many regimes, falter when predicting the equation of state necessary to accurately describe neutron stars and other compact objects. This equation of state, which defines the relationship between pressure and density, remains poorly constrained under the immense pressures found within these celestial bodies. Consequently, interpreting observational data – such as the masses and radii of neutron stars derived from gravitational wave detections or electromagnetic radiation – becomes significantly challenging. Discrepancies between predicted and observed properties suggest that existing models may be oversimplified, failing to capture crucial aspects of the behavior of matter at these densities and necessitating exploration of more complex and nuanced descriptions of the strong nuclear force.

The behavior of matter at densities exceeding that of atomic nuclei-conditions found within neutron stars-presents a significant challenge to physicists attempting to model the transition from familiar hadronic matter to exotic states. Current theoretical frameworks struggle to precisely delineate this phase transition, potentially involving quarkyonic matter – a phase where quarks and gluons are deconfined but still exhibit collective behavior – or even fully deconfined quark matter. This uncertainty stems from the complexities of quantum chromodynamics (QCD) at high densities, where perturbative calculations fail and non-perturbative approaches are required. Accurately mapping this transition is crucial, as the equation of state-the relationship between pressure and density-is directly influenced by the matter’s composition and profoundly impacts the observed properties of neutron stars, such as their mass-radius relationship and cooling rates. Determining the precise nature of this transition remains a central goal in nuclear astrophysics, demanding innovative theoretical models and stringent tests against observational data.

The quarkyonic quantum Monte Carlo (QQMC) model reveals that both bulk <span class="katex-eq" data-katex-display="false">k_b^<i></span> and shell <span class="katex-eq" data-katex-display="false">k_s^</i></span> nucleon momenta exhibit density dependence. — The quarkyonic quantum Monte Carlo (QQMC) model reveals that both bulk $k_b^<i>$ and shell $k_s^</i>$ nucleon momenta exhibit density dependence.

A Hybrid Phase Emerges: The Nature of Quarkyonic Matter

Quarkyonic matter is theorized to be an intermediate phase of quantum chromodynamics (QCD) matter, existing between hadronic matter and quark matter. Unlike typical hadronic matter where quarks are confined within baryons and mesons, and unlike quark matter where quarks are deconfined, quarkyonic matter features a coexistence of both. Specifically, it’s characterized by the presence of large baryon masses and a momentum distribution where both nucleon-like (hadronic) and quark-like momentum components contribute significantly. This means that while baryons still exist, their internal quark structure begins to exhibit properties of a deconfined quark-gluon plasma, leading to a hybrid state with unique thermodynamic and transport characteristics. The existence of this phase is predicated on density ranges where quark degrees of freedom begin to emerge without fully dissolving the hadronic structure.

The momentum distribution of particles in quarkyonic matter deviates from traditional single-component models due to the simultaneous presence of nucleons and liberated quarks. Specifically, the distribution exhibits a broadening and complex structure resulting from the superposition of hadronic contributions – typically characterized by a relatively narrow momentum spread – and the wider distribution arising from deconfined quarks with higher average momenta. This dual nature manifests as a distribution that is no longer solely describable by a Fermi-Dirac function, but rather requires a combination of functions or a more complex analytical form to accurately represent the contributions from both baryon and quark degrees of freedom. Quantitative analysis of this momentum distribution, and its evolution with density, is essential for extracting the equation of state and transport properties of matter in this intermediate regime.

Precise modeling of the dual momentum distribution in quarkyonic matter is essential for accurately predicting the equation of state at intermediate densities, specifically between approximately 200 MeV and 600 MeV. This density range represents a transitional region where hadronic structures begin to dissolve and quark degrees of freedom emerge, but a complete transition to fully deconfined quark matter has not yet occurred. The equation of state in this region significantly impacts macroscopic properties like neutron star mass-radius relationships and the dynamics of heavy-ion collisions; inaccuracies in the momentum distribution will propagate to these observable quantities. Current theoretical efforts focus on developing effective models that capture the interplay between nucleon and quark momenta, employing techniques such as holographic duality and effective field theories to constrain the parameters governing this complex behavior.

Calculations using both Quantum Monte Carlo (QMC) and Quadratic Quantum Monte Carlo (QQMC) models demonstrate that binding energy per nucleon <span class="katex-eq" data-katex-display="false">\varepsilon_b</span> decreases with increasing baryon density ratio <span class="katex-eq" data-katex-display="false">\rho_N/\rho_0</span>. — Calculations using both Quantum Monte Carlo (QMC) and Quadratic Quantum Monte Carlo (QQMC) models demonstrate that binding energy per nucleon $\varepsilon_b$ decreases with increasing baryon density ratio $\rho_N/\rho_0$ .

Unveiling Substructure: The Quark-Meson Coupling Model

The Quark-Meson Coupling (QMC) model deviates from traditional nuclear matter calculations by treating nucleons not as fundamental particles, but as emergent phenomena arising from underlying quark constituents. This approach explicitly incorporates the relativistic quark substructure of nucleons – specifically, utilizing a constituent quark model where each nucleon is composed of three valence quarks – and calculates the properties of dense matter based on interactions between these constituent quarks, as well as between the nucleons themselves. This contrasts with non-relativistic approaches and allows for the investigation of how the internal structure of nucleons influences the macroscopic properties of dense matter, particularly at high densities where traditional nuclear models may break down. The model’s framework is built upon relativistic quantum field theory, enabling the investigation of phenomena governed by $QCD$ at a scale accessible through effective nuclear interactions.

The Quark-Meson Coupling (QMC) model represents nucleons as composed of quarks bound by a potential derived from a relativistic framework. Specifically, the internal wavefunctions of the quarks are described by a Gaussian function, mathematically formulated as a solution to the Dirac equation-the relativistic analogue of the Schrödinger equation-subject to a harmonic oscillator potential. This potential incorporates both scalar and vector components, representing the confining force between quarks and ensuring Lorentz covariance. The spatial distribution of quarks within the nucleon is thus defined by $\psi(\mathbf{r}) \propto e^{-(\mathbf{r}^2)/2\sigma^2}$ , where σ is a parameter defining the size of the nucleon, determined by solving the Dirac equation with the specified potential. This relativistic Gaussian wavefunction forms the foundation for calculating nucleon properties and interactions within the QMC model.

The Quark-Meson Coupling (QMC) model enables the calculation of macroscopic properties by explicitly accounting for both nucleon and quark degrees of freedom. Effective nucleon mass is determined by considering the interactions between quarks within nucleons and the surrounding meson fields, leading to a mass modification relative to the free nucleon mass. The equation of state (EOS) is then derived from a statistical ensemble incorporating both nucleons and quarks, utilizing Fermi-Dirac statistics and accounting for interactions mediated by meson exchange. This approach allows for the determination of the pressure as a function of density, and consequently, the stiffness of the nuclear matter EOS, which is crucial for understanding the properties of neutron stars and the dynamics of heavy-ion collisions. Calculations incorporate contributions to the energy density from both the nucleons, treated as composite objects, and the explicit quark degrees of freedom, ensuring a consistent description of dense baryonic matter.

Calculations within the Quark-Meson Coupling (QMC) model indicate a substantial quantitative contribution from nuclear interactions to the stiffening of the equation of state (EOS) in the quarkyonic regime. Specifically, the repulsive character of the ω-meson exchange, and to a lesser extent the ρ-meson exchange, between nucleons significantly increases the pressure at high densities. This effect is not solely attributable to quark degrees of freedom; rather, the residual nucleon-nucleon interactions, modeled through meson exchange potentials, provide a dominant contribution to the overall stiffness of the EOS at comparable densities, and impacts predictions for neutron star radii and the collective flow in heavy-ion experiments. $P = K(\rho - \rho_0)^{\gamma}$

The Quark-Meson Coupling (QMC) model provides a unified framework for describing the properties of dense nuclear matter across a range of compositions. Calculations are performed for both symmetric nuclear matter, consisting of equal numbers of protons and neutrons, and for pure neutron matter, relevant to the core of neutron stars. This adaptability is achieved through consistent treatment of nucleon and quark degrees of freedom within the model’s Lagrangian, allowing for the prediction of equation of state parameters and nucleon properties as a function of density for both systems. The model’s applicability to both compositions facilitates a broader understanding of the behavior of matter under extreme conditions, independent of proton-to-neutron ratio.

Sound velocity <span class="katex-eq" data-katex-display="false">\sqrt{s^2}</span> exhibits a dependence on the baryon density ratio <span class="katex-eq" data-katex-display="false">\rho_N/\rho_0</span>, differing between the QQMC and GQ models across proton fractions of 0.6, 0.7, and 0.8 fm. — Sound velocity $\sqrt{s^2}$ exhibits a dependence on the baryon density ratio $\rho_N/\rho_0$ , differing between the QQMC and GQ models across proton fractions of 0.6, 0.7, and 0.8 fm.

From Theory to Observation: Bridging the Gap to Neutron Stars

The Quantum Monte Carlo (QMC) model posits that quarkyonic matter – a hypothesized state existing within neutron stars – exhibits a distinct relationship between its sound velocity and binding energy, fundamentally dictating its stability and observable properties. This model predicts a specific sound velocity, essentially a measure of how quickly disturbances propagate through the matter, and a corresponding binding energy, which determines how tightly the quarks are bound together. Crucially, these parameters aren’t static; they vary with density and influence the equation of state, impacting the overall pressure and resistance to compression within the neutron star. A precise understanding of this interplay, as revealed by QMC simulations, is essential for connecting theoretical predictions to astronomical observations and ultimately deciphering the exotic physics occurring at extreme densities, potentially revealing clues about the behavior of matter beyond the limits of terrestrial experiments.

Recent investigations utilizing quantum Monte Carlo methods suggest a compelling relationship between nuclear and quark matter densities within neutron stars. These calculations demonstrate that the quark saturation density $\rho_{sat}$ , representing the point at which quark degrees of freedom become significant, exceeds the nuclear saturation density $\rho_0$ . This finding isn’t merely a numerical observation; it’s intrinsically linked to the relativistic treatment of particles within the model and the complex interplay of nuclear interactions. The increased $\rho_{sat}$ implies that quarkyonic matter, a phase bridging hadronic and quark matter, emerges only at densities exceeding those typically found in atomic nuclei. Consequently, this provides constraints on the equation of state of dense matter, impacting models of neutron star structure and potentially influencing observations of these exotic stellar remnants.

The predicted behavior of sound velocity ( $v_s^2$ ) within quarkyonic matter reveals a critical point at the quark saturation density ( $\rho_{sat}$ ). Calculations demonstrate that as matter approaches this density, $v_s^2$ exhibits a singularity, indicating a dramatic shift in the equation of state. This stiffening isn’t merely a numerical quirk; it suggests a fundamental change in the material’s resistance to compression. Prior to reaching $\rho_{sat}$ , the matter behaves as a relatively fluid medium, but beyond this point, it becomes increasingly resistant to further compression, potentially influencing the maximum mass and radius of neutron stars and providing insights into the exotic phases of matter present at extreme densities.

Investigations into quarkyonic matter reveal a significant increase in binding energy $ε_b$ as the nucleon density $ρ_N$ surpasses the quark saturation density $\rho_{sat}$ . This phenomenon indicates a transition to a more stable, tightly bound state of matter. Notably, this rise in binding energy is considerably more pronounced in pure neutron matter compared to symmetric nuclear matter, suggesting that the neutron-rich environment within neutron stars plays a crucial role in driving this transition. The disparity implies that the strong interaction’s influence differs between neutron-dominated and balanced nuclear compositions at these extreme densities, offering a potential pathway to constrain the equation of state of matter found within the cores of these celestial objects.

Current quarkyonic matter calculations, while providing valuable insights, stand to gain significantly from methodological advancements, particularly through integration with models like IdyllIq. This approach introduces a dual momentum space distribution, allowing for a more nuanced representation of particle interactions at extreme densities – conditions found within neutron stars. By refining the model’s capacity to accurately capture these interactions, researchers anticipate enhanced predictive power regarding the equation of state of matter at these densities. Such improvements are crucial for correlating theoretical calculations with observational data from neutron star observations, ultimately deepening the understanding of the behavior and stability of these enigmatic celestial objects and the exotic states of matter they contain.

The exploration of quarkyonic matter, as detailed in this study, necessitates a careful consideration of underlying structures and their interactions. It’s a process demanding patience, as quick conclusions regarding the equation of state or the onset of quark saturation can mask structural errors. This aligns with Galileo Galilei’s observation: “You cannot teach a man anything; you can only help him discover it himself.” The research doesn’t present definitive answers, but rather a framework for understanding how nuclear interactions contribute to the stiffness of matter at extreme densities, encouraging further investigation and refinement of the relativistic quark model.

Looking Ahead

The exploration of quarkyonic matter, as presented, reveals a landscape where the familiar boundaries of hadronic physics begin to dissolve. However, the model’s reliance on a specific relativistic quark model introduces a necessary caveat. The equation of state, while demonstrating a plausible path towards quark saturation, remains sensitive to the underlying assumptions concerning quark interactions and constituent quark masses. Future refinement demands a systematic investigation of these parameters, alongside comparisons with independent theoretical frameworks – perhaps those less anchored in constituent quark pictures.

A compelling avenue for further study lies in bridging the gap between this model-dependent approach and emergent phenomena observed in heavy-ion collisions. The calculated sound velocity, a critical parameter linking microscopic interactions to macroscopic behavior, warrants detailed comparison with experimental data. Discrepancies, rather than being dismissed as noise, should be treated as opportunities to refine the model’s description of strong interactions at finite density. The observed stiffening of the equation of state, attributable to nuclear interactions, deserves scrutiny in the context of neutron star mergers – a setting where extreme densities offer a unique testing ground.

Ultimately, the quest to understand quarkyonic matter is not merely about identifying a specific phase of matter. It is about probing the fundamental nature of confinement and deconfinement, and the subtle interplay between quarks and hadrons. Each refinement of the model, each comparison with experiment, is a step towards deciphering the patterns hidden within the seemingly chaotic realm of strong interactions.

Original article: https://arxiv.org/pdf/2603.19839.pdf

Contact the author: https://www.linkedin.com/in/avetisyan/

Unveiling the Density Puzzle: Beyond Conventional Models

A Hybrid Phase Emerges: The Nature of Quarkyonic Matter

Unveiling Substructure: The Quark-Meson Coupling Model

From Theory to Observation: Bridging the Gap to Neutron Stars

Looking Ahead

See also: