Quarkyonic Quark-Meson Coupling Model for Nuclear and Neutron Matter

arXiv:2512.04505 · nucl-th, astro-ph.HE, hep-ph · Submitted 2025-12-04 · Read on arXiv

Listen

Radio episode about this paper

Transcript

Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.

Vera: I'm Vera, and with me are Jocelyn and Subrahmanyan, guest researcher.

Jocelyn: Today's paper: "Quarkyonic Quark-Meson Coupling Model for Nuclear and Neutron Matter".

Vera: A novel nuclear model, called the quarkyonic quark-meson coupling (QQMC) model,

Jocelyn: First, who's behind it and why it matters.

Title and authors: Vera: So we're starting with this paper titled "Quarkyonic Quark-Meson Coupling Model for Nuclear and Neutron Matter," and the authors are Koichi Saito, Tsuyoshi Miyatsu, and Myung-Ki Cheoun. It sounds like they're trying to build a model that bridges the gap between how we understand nuclear matter at low densities and what happens in extremely dense environments.

Jocelyn: I was looking at the title, and it immediately makes me think about unifying different ways of looking at this problem, which is exactly what the authors seem to be doing by combining dual quarkyonic and quark-meson coupling approaches. It suggests a comprehensive look across all density regimes.

Subrahmanyan: From a theoretical perspective, uniting these two distinct frameworks is ambitious; it implies that the physics governing matter transitions from being purely hadronic to being dominated by quarks can be described coherently within one structure, which is very important for connecting terrestrial experiments to astrophysical observations.

Vera: Exactly! The authors are taking something that covers low density up to the crossover region and trying to make it quantitative enough so we can actually test it against real-world data from both heavy-ion collisions and neutron star observations.

Jocelyn: And I wonder what the direct implications are for our understanding of matter under extreme conditions; they seem focused on how these quark degrees of freedom influence the bulk properties of nuclear matter.

Subrahmanyan: The implication is that we might be able to get a more robust picture of the Equation of State in that tricky crossover region, which is where many current models struggle to provide a consistent description.

The paper's summary: Vera: So, what the paper actually summarizes is this novel model, the QQMC model, which incorporates Pauli blocking at the quark level alongside how the nucleon structure changes in a medium. It’s designed to handle everything from low density right up to that crossover point where hadronic matter starts turning into something else.

Jocelyn: That sounds like they are addressing some real theoretical headaches that plague nuclear physics; it's not just about fitting data at one point, but creating a framework that works across a whole density spectrum.

Subrahmanyan: The key summary point here is the unification aspect; by combining the dual quarkyonic model with the QMC model, they aim to describe physical quantities like pressure and sound velocity consistently using these quark degrees of freedom.

Vera: And they specifically mention using a relativistic, gaussian quark wavefunction for nucleon structure instead of some other functions used in other models, which is a methodological choice they're making to describe those nucleons in the medium.

Jocelyn: So, when you look at the results described in the summary, it seems like they are tackling the inherent issues where simpler models often show discontinuities or divergences at quark saturation density.

Subrahmanyan: That difficulty with singularities is central; they introduce an infrared regulator to smooth out that behavior and then combine everything into this final QQMC model to handle those issues quantitatively.

The paper's improvements: Vera: The paper details some specific improvements they made, like introducing a regulator to smear the sharp Fermi surface using something like the theta function in nucleon momentum distributions. This is a direct fix for those singular behaviors we talked about earlier.

Jocelyn: That regulator seems pretty practical; it's an adjustment that allows them to keep the physical quantities continuous at saturation density, which is a significant step forward from the naive ideal Fermi gas picture they started with.

Subrahmanyan: The authors also found that when they include the nuclear interaction within this QQMC framework, specifically by considering nucleon size in matter, the quark saturation density actually shifts lower compared to earlier calculations.

Vera: That shift is interesting; they state that for a nucleon radius of rp = zero point six or zero point eight fm, the quark saturation density becomes about three point six or one point five times rho0 in symmetric nuclear matter, which is a specific quantitative result they derived from this refinement <ref:2512.04505#pg1>.

Jocelyn: So it shows that the input parameters, like the nucleon size, have a substantial effect on how high the quark saturation density ends up being in the model.

Subrahmanyan: And they've pointed out that this entire framework is important because it provides a unified way to incorporate both Pauli blocking at the quark level and medium modification of nucleon structure, which is what makes it suitable for describing matter in that crossover region between hadronic and quark degrees of freedom <ref:2512.04505#pg1>.

Conclusion: Vera: So to wrap up the QQMC model, the authors conclude that this unified framework successfully incorporates both quark-level Pauli blocking and medium modification of nucleon structure, which they say is essential for describing dense matter in that crossover region between hadronic and quark degrees of freedom.

Jocelyn: They've shown that with these additions, they can achieve a quantitative description of physical quantities like sound velocity consistent with neutron star observations, and pressure consistent with heavy-ion collision data.

Subrahmanyan: The implication for the cosmic picture is that this model offers a way to link the microscopic structure of quarks and mesons to macroscopic observables like neutron star stiffness, which is what we need when we look at those massive pulsars.

Vera: It's certainly a strong connection; they've even shown how the sound velocity in neutron star matter reaches a maximum and then decreases gradually, which aligns with Bayesian inference analysis results.

Jocelyn: It’s really encouraging to see this level of detail in connecting these disparate areas of physics; it makes the model feel more grounded in observable constraints.

Subrahmanyan: I just think the work highlights how incorporating nuclear interactions is important to get quantitative accuracy; for instance, the sound velocity in neutron star matter can be explained by the QQMC model where it reaches a maximum and then decreases gradually, consistent with Bayesian inference analysis results <ref:2512.04505#pg2>.

Koichi Saito, Tsuyoshi Miyatsu, Myung-Ki Cheoun

Department of Physics and Astronomy, Tokyo University of Science · Department of Physics and OMEG Institute, Soongsil University

nucl-th, astro-ph.HE, hep-ph

Submitted: 2025-12-04

Updated: 2026-07-30

Comments: 50 pages, 18 figures, 3 tables, accepted for publication in PRC

Journal ref: Phys. Rev. C 114, 045201 (2026)

DOI: 10.1103/s71s-85r5

License: http://creativecommons.org/licenses/by/4.0/

Importance score: 74/100

The gist: A novel nuclear model, called the quarkyonic quark-meson coupling (QQMC) model, is constructed by uniting dual quarkyonic and quark-meson coupling approaches to describe nuclear matter from low

Key concepts

Quarkyonic Quark-Meson Coupling (QQMC) Model: A novel nuclear model that combines two approaches: dual quarkyonic and quark-meson coupling. It describes nuclear matter from low density up to the crossover region by treating nucleons as confined quarks interacting with scalar and vector mesons, allowing for a quantitative description of physical properties.
Pauli Blocking at the Quark Level: This effect occurs when quarks are confined within nucleons. Pauli blocking means that quarks cannot occupy states that are already filled, which modifies the quark momentum distribution. This is crucial for accurately describing dense matter and ensuring physical consistency in the model.
Quark Saturation Density ($ ho_{sat}$): This is a critical density where quarks with momenta between 0 and $q_b$ become fully occupied, forming a 'bulk Fermi sea.' The value of this density is sensitive to the chosen quark wavefunction; using the relativistic Gaussian function results in a higher $ ho_{sat}$ than in simpler models.
Medium Modification of Nucleon Structure: This refers to how the structure and properties of individual nucleons change when embedded in dense nuclear matter. In QQMC, mean fields from scalar and vector mesons modify the effective quark mass and single-particle energy, reflecting this change in the medium.

Terminology

Bayesian Inference — 7× in this paper · explained in 5 episodes
Unified Framework — 3× in this paper · explained in 5 episodes
Core — 2× in this paper · explained in 2 episodes
Accuracy — 1× in this paper · explained in 1 episode
Consistency — 1× in this paper · explained in 3 episodes
Degeneracy — 1× in this paper · explained in 3 episodes

Summary

A novel nuclear model, called the quarkyonic quark-meson coupling (QQMC) model, is constructed by uniting dual quarkyonic and quark-meson coupling approaches to describe nuclear matter from low density up to the crossover region. This framework incorporates both Pauli blocking at the quark level and medium modification of nucleon structure, allowing for a quantitative description of physical quantities like sound velocity consistent with neutron star observations and pressure consistent with heavy-ion collision data.

The gist

A novel nuclear model based on the quark degrees of freedom, which can cover a wide range of nuclear densities, from low density to the crossover region.

Model Construction and Components

The model unifies the dual quarkyonic model with the quark-meson coupling (QMC) model to construct a novel nuclear model based on the quark degrees of freedom.
It uses a relativistic, gaussian quark wavefunction to describe nucleon structure, which is adopted in this paper as a gaussian function instead of the Yukawa-type function used in the IdylliQ model.
The QMC model is extended to include quark degrees of freedom by considering mean fields of scalar and vector mesons interacting with confined quarks, where the effective quark mass and single-particle energy are modified by these fields.
The resulting framework, the QQMC model, includes the effect of Pauli blocking at the quark level as well as the effect of scalar polarizability of nucleon in medium.

Quarkyonic Phase and Momentum Distributions

The quark saturation density, denoted as quark saturation density, emerges above a critical density where quarks with momenta between 0 and qb are fully occupied, leading to the bulk Fermi sea.
In the ideal Fermi gas picture (GQ model), the quark momentum distribution is described by relations like Eq. (12) and Eq. (13), ensuring duality through the sum rule, Eq. (2).
The model introduces characteristic momenta in the quarkyonic phase, "kb and ks, where kb defines the under-occupied bulk part at low momentum and ks gives the upper bound of the shell structure at high momentum."

Addressing Singularities and Interactions

In the naive GQ model, physical quantities like chemical potential and pressure are discontinuous, and sound velocity diverges at quark saturation density.
A minimal correction is introduced by a regulator that smears the sharp Fermi surface using a function like θ(kb − k) in nucleon momentum distributions to remove singular behavior.
The boundary conditions for determining kb and ks are taken from the IdylliQ model, which are then modified using the regulator, ensuring that the physical quantities (except v2s) are continuous at ρsat.

Physical Predictions and Consistency Checks

The QQMC model predicts a sound velocity consistent with data inferred from observed neutron star data by using a neural network model or Bayesian inference analysis.
The pressure calculated by the QQMC model lies within the range deduced from experimental flow data in heavy-ion collisions at high energy and kaon production data.
The choice of nucleon size, specifically rp = 0.6–0.8 fm, seems most suitable for describing dense nuclear matter in the QQMC model, as various physical quantities depend on it accordingly.
The inclusion of nuclear interaction is important to describe physical quantities quantitatively; for instance, the sound velocity in neutron star matter can be explained by the QQMC model where it reaches a maximum and then decreases gradually, consistent with Bayesian inference analysis results.

Summary of Findings

The quark saturation density is very sensitive to the quark wavefunction; for the relativistic gaussian wavefunction used here, the value of ρsat turns out to be higher than that in the NR case.
When nuclear interaction is included (QQMC model), the quark saturation density is reduced because the nucleon size is swollen in matter, leading to a value of ρ∗sat = 1.46ρ0 in SNM for rp = 0.8 fm, which is lower than the value of 1.63 by about 10%.
The model provides a unified framework that incorporates both quark-level Pauli blocking and medium modification of nucleon structure, which is essential for describing dense matter in the crossover region between hadronic and quark degrees of freedom.

Future Directions

Further investigation is needed to determine if the boundary condition used to find kb and ks is sufficient, and a general rule for constructing the quark momentum distribution, fQ(q), in the postsaturation region satisfying Fermi statistics and energy minimization must be found.
The singular behavior around quark saturation density should fundamentally be resolved by QCD dynamics rather than relying on ad hoc regulators like the one introduced here.
It may be possible to include temperature and apply the model to EoSs for neutron stars with hyperons in future work.

Improvements for AI systems

As a fastidious and diligent researcher, I have analyzed the provided scientific paper, Quarkyonic Quark-Meson Coupling Model for Nuclear and Neutron Matter. This work proposes a novel nuclear model (the QQMC model) that unifies the dual quarkyonic picture with the quark-meson coupling (QMC) model.

Here are specific improvements to AI systems based on the theoretical framework presented in this paper, detailing what the improved AI system can achieve:

Based on the QQMC Model and its underlying physics, here are specific improvements for AI systems:

The improved AI system can accurately simulate and predict the Equation of State (EoS) of dense nuclear matter across a wide density range (from low density to crossover). It can specifically model the transition from baryonic to quark matter using the quarkyonic phase concept.
The system can perform high-precision calculations of transport coefficients, such as pressure, energy density, chemical potential, and sound velocity, in both symmetric nuclear matter (SNM) and pure neutron matter (PNM).
The AI can predict neutron star properties by calculating the sound velocity profile as a function of density. Crucially, it can simulate the stiffening of the EoS slightly above saturation density via quark degrees of freedom, which is essential for matching constraints from heavy-ion collisions and gravitational wave observations (like GW170817).
The AI can incorporate relativistic effects and nucleon structure modifications (via the Gaussian quark wavefunction and in-medium mass/polarizability) into its calculations, allowing it to predict how nuclear properties change when nucleons are embedded in a dense medium.
The system can be trained to perform Bayesian inference on observational data from neutron stars (e.g., mass and radius constraints) to constrain the model parameters (like the nucleon radius parameter, rp), thereby testing the model's predictive power against astrophysical data.
It can analyze and predict discrepancies between different theoretical approaches (e.g., comparing results from the naive Gaussian Quarkyonic (GQ) model versus the fully coupled QQMC model) to identify which physical ingredients—like Pauli blocking or nuclear interactions—are most critical for quantitative accuracy at high densities.

In summary, an AI system based on this paper would transition from being a general physics simulator to a specialized tool capable of:

Accurately modeling the phase transition region in dense matter.
Generating realistic EoS curves consistent with multi-messenger astronomy (neutron stars).
Providing quantitative predictions for observables used in terrestrial experiments (heavy-ion collisions).

Abstract

We unite the dual quarkyonic model with the quark-meson coupling (QMC) model to construct a novel nuclear model based on the quark degrees of freedom, which can cover a wide range of nuclear densities, from low density to the crossover region. In the model, the relativistic, gaussian quark wavefunction is used to describe the nucleon structure. We first evaluate the energy density, chemical potential, pressure and sound velocity within the ideal Fermi gas picture. In this case, those physical quantities are discontinuous or divergent at the quark saturation density, where the quarkyonic phase emerges. To remove such singular behavior, we next introduce an infrared regulator, and combine the dual quarkyonic model and the QMC model to include the nuclear interaction -- we call it the quarkyonic quark-meson coupling (QQMC) model. In this model, the quark saturation density depends strongly on the nucleon size. For example, when r p = 0.6, (0.8) fm, where r p is the root-mean-square radius of the proton, the quark saturation density is about 3.6,(1.5) times ρ 0 in symmetric nuclear matter, where ρ 0 is the nuclear saturation density. Furthermore, the nuclear interaction plays an important role in considering physical quantities quantitatively. In fact, the QQMC model can produce the sound velocity which is consistent with that inferred from the observed data of several neutron stars. Furthermore, pressure in symmetric or pure neutron matter deduced from the experiments of heavy-ion collisions at high energy can be explained by the QQMC model as well. We discuss in detail the formulation for the QQMC model and the physical quantities calculated by the model.