Massive hybrid stars within the extended three-flavor quark-meson diquark model

This article demonstrates that the extended three-flavor quark-meson-diquark model, particularly through the inclusion of vector and axial-vector mesons, successfully describes massive hybrid stars with masses exceeding $2M_{\odot}$ and radii consistent with astrophysical observations by generating a sufficiently stiff equation of state and a double-peak structure in the speed of sound, which is driven by the decreasing mass of the strange quark at high densities.

The Gist

Imagine the universe as a vast cosmic kitchen. In this kitchen, there are two main types of ingredients: the "everyday" substances found in atoms (such as protons and neutrons), and the "super-dense" substances that occur only in the hearts of dead stars, known as neutron stars.

For decades, scientists have tried to write a cookbook (a so-called equation of state) that explains how these ingredients behave when compressed with unimaginable force. The problem is that the "super-dense" ingredients are so strange that our usual cookbooks fail.

This article introduces a new, improved cookbook called the Extended Three-Flavor Quark-Meson Diquark (EQMD) model. Here is how it works, simply explained:

1. The Ingredients: From Solid Blocks to Swirling Soup

In normal matter, protons and neutrons are like solid Lego blocks. However, at the center of a massive neutron star, the pressure is so high that these blocks are crushed until they melt into a swirling soup of their smaller constituents: quarks.

The authors' new model treats this soup not merely as chaotic disorder, but as a structured mixture containing:

• Quarks: The tiny fundamental particles.
• Mesons: Particles that act like the "glue" holding things together.
• Diquarks: Pairs of quarks that stick together like dance partners.
• Vector mesons: A new type of "glue" that the authors added to the mixture.

The Analogy: Imagine the old models as an attempt to describe a dance floor with only two types of dancers. The authors realized that a crucial group was missing. By adding vector mesons (the new dancers), the dance floor suddenly makes sense. Without them, the crowd would be too loose and wobbly; with them, the crowd becomes stiff and stable enough to bear a heavy weight.

2. The Challenge: Building a Star That Does Not Collapse

Neutron stars are incredibly heavy. Some weigh twice as much as our Sun but are compressed into a sphere the size of a city. If the "recipe" for the star's core is too soft (like jelly), the star's own gravity will crush it into a black hole. If it is too stiff (like a steel beam), the mathematics will not match what we observe in the sky.

The authors tested their new recipe against real observations from telescopes and gravitational wave detectors (such as LIGO). They asked: "Can we build a star with this recipe that is heavy enough to match the heaviest stars we have actually seen?"

The Result: Yes. By carefully adjusting the "spices" (the parameters in their model), they found that their recipe produces a star that:

• Is stiff enough in the middle to support a mass of about 2 Suns.
• Is soft enough at the outermost edges to match the size (radius) of the stars we have measured.

3. The "Double-Peak" Puzzle

One of the most interesting discoveries in this article concerns the speed of sound within these stars.

Normally, one might think that sound travels faster in denser materials. But in these stars, something strange happens to the speed of sound: it rises, then falls, and then rises again. This creates a "double-peak" shape.

The Analogy: Imagine driving a car up a mountain. You accelerate, then hit a muddy patch where you slow down, and then you reach a smooth highway where you accelerate again.

• Why the slowdown? The article explains that this happens due to strange quarks. As pressure increases, the "strange" particles in the star begin to lose their mass (they "melt"). This melting causes a temporary dip in the star's stiffness and slows down the speed of sound.
• Why the second peak? Once the strange particles have completely melted, the star becomes stiff again, and the speed of sound shoots upward until it eventually settles into a steady rhythm.

4. What This Tells Us About the Universe

The authors conclude that if we find a neutron star heavier than 2 Suns, it almost certainly possesses a quark core.

• The outer layer consists of normal nuclear matter (Lego blocks).
• The inner core (starting at about four times the density of an atomic nucleus) consists of this exotic quark soup.

They also found that the transition from the "Lego block" layer to the "quark soup" layer occurs smoothly, rather than with a sudden, jerky jump.

Summary

In short, this article presents a new, more complete "recipe" for the densest matter in the universe. By adding a missing ingredient (vector mesons) and accounting for the behavior of "strange" particles, the authors have created a model that successfully explains how the heaviest neutron stars can exist without collapsing. It suggests that the hearts of these stars are not just solid blocks, but a complex, melting, and re-stiffening soup of quarks.

Technical Summary

Technical Summary: Massive Hybrid Stars in the Extended Three-Flavor Quark-Meson-Diquark Model

Problem Statement
The internal composition of neutron stars, particularly those with masses exceeding $2M_\odot$, remains a central question in dense QCD. While General Relativity provides the framework (via the Tolman-Oppenheimer-Volkoff equation) to relate a star's mass and radius to its Equation of State (EoS), the EoS itself is unknown at the extreme densities found in stellar cores. Lattice QCD is currently unable to calculate thermodynamic quantities at finite baryon density due to the sign problem. Consequently, researchers rely on effective field theories. A significant challenge lies in bridging the gap between low-energy nuclear matter (described by chiral effective field theory, $\chi$EFT) and high-density quark matter (described by perturbative QCD, pQCD). Existing models, such as the Nambu-Jona-Lasinio (NJL) model, often suffer from regularization artifacts at finite density or fail to generate an EoS stiff enough to support observed massive neutron stars without violating causality or pQCD conditions.

Methodology
The authors employ the Extended Three-Flavor Quark-Meson-Diquark Model (EQMD), a renormalizable low-energy effective model for QCD. The effective degrees of freedom include quarks, scalar and pseudoscalar mesons, diquarks, and crucially vector and axial mesons.

1. Model Construction: The Lagrangian density is constructed to respect the symmetries of Three-Flavor QCD ($SU(3)_L \times SU(3)_R \times SU(3)_c \times U(1)_B \times U(1)_A$). The authors extend the standard Quark-Meson-Diquark model by minimally adding mass terms and quartic terms for the new vector and axial-vector meson fields ($L_\mu$ and $R_\mu$). These new fields couple to quarks via a coupling constant $g_\omega$ and effectively shift the quark chemical potentials.
2. Thermodynamics: The EoS is calculated in the Mean-Field Approximation at zero temperature ($T=0$). The model manually enforces local charge neutrality for electric and color charges, as it lacks dynamical gauge fields. The thermodynamic potential is renormalized using the on-shell scheme, with meson masses and decay constants fitted to physical observables.
3. Hybrid Star Construction: To model hybrid stars, the authors couple the high-density EQMD EoS to a low-density nuclear EoS (specifically the GMSR EoS from the CompOSE library). They use a construction where the transition between nuclear matter and quark matter is determined by equating the particle densities ($n_{EQMD} = n_{nucl}$) at a transition chemical potential $\mu_t$, ensuring a continuous baryon density. This approach avoids the discontinuities of a Maxwell construction but allows for a smooth quark-hadron transition.
4. Parameter Space: The model contains several free parameters, most notably the quark-vector coupling $g_\omega$. The authors test three parameter sets with varying $g_\omega$ values to investigate the stiffness of the EoS.

Main Contributions and Results

• Stiffness of the EoS: The inclusion of vector and axial-vector mesons is identified as essential. These degrees of freedom enable the model to generate an EoS that is sufficiently stiff at intermediate densities to support massive stars and softens at higher densities to remain consistent with pQCD.
• Consistency with Observations: The resulting hybrid stars satisfy current astrophysical constraints, including:
  • Maximum mass constraints from pulsars PSR J1614-2230 ($1.97 M_\odot$), PSR J0348+0432 ($2.01 M_\odot$), PSR J2215+5135 ($2.27 M_\odot$), and PSR J0952-0607 ($2.35 M_\odot$).
  • Mass and radius constraints from NICER observations of PSR J0030+0451 and PSR J0740+6620.
  • One parameter set (Set 1) fits all observations, while Sets 2 and 3 fit all except the most massive PSR J0952-0607.
• Central Densities and Quark Cores: The authors find that stars with masses $M \gtrsim 2M_\odot$ possess a quark core with a central baryon density $n_B \geq 3.9 n_{sat}$ (where $n_{sat} \approx 0.165 \, \text{fm}^{-3}$). For the maximum mass configurations, the radii of the quark cores are estimated to be approximately 3.8 km to 4.6 km.
• Structure of the Speed of Sound: The square of the speed of sound ($c_s^2$) exhibits a characteristic double-peak structure.
  • The first peak arises from the transition from the nuclear EoS to the quark matter EoS.
  • A sharp drop follows, caused by the rapid melting of the strange quark condensate as the chemical potential increases.
  • A second peak occurs at higher densities, where $c_s^2$ exceeds the conformal limit ($1/3$), before approaching this limit from above at asymptotic densities.
  • This behavior contrasts with the NJL model (where $c_s \to 1$) and weak-coupling extensions (which suggest an approach from below to the limit), but aligns with the specific dynamics of the EQMD model, where the mass of the strange quark decreases with increasing $\mu_B$.
• Determination of the Bag Constant: In contrast to traditional bag models where the bag constant $B$ is a free parameter, here $B$ is determined self-consistently by the matching condition at the transition density. The derived values ($B^{1/4} \approx 234\text{--}242$ MeV) are consistent with estimates from NJL models and QCD-based calculations.

Significance and Claims
The work claims that the EQMD model, specifically with the addition of vector mesons, offers a viable, renormalizable framework for describing the deconfined phase of matter in neutron stars. The primary significance lies in demonstrating that it is possible to construct hybrid stars that simultaneously:

1. Are consistent with the stiffness required to support $2M_\odot$ stars.
2. Are consistent with the softening demanded by pQCD constraints at high densities.
3. Are consistent with $\chi$EFT constraints at low densities.

The authors emphasize that the double-peak structure of the speed of sound is a specific prediction of the three-flavor model, driven by the melting of the strange quark condensate. They conclude that their results suggest a quark core is likely present in the most massive observed neutron stars. The work positions these results as a foundation for future Bayesian analyses, where the parameter sets identified here can serve as priors for more rigorous statistical inferences using gravitational wave and X-ray data. The work does not claim to definitively prove the existence of quark matter, but rather shows that a specific effective field theory model can satisfy all current observational constraints without theoretical inconsistency.