Crossover Equation of State Constrained by Astronomical Observations and pQCD

This paper investigates a hadron-quark crossover equation of state for neutron stars by constraining Nambu--Jona-Lasinio model parameters with perturbative QCD and astronomical observations, finding that such crossover models significantly enhance maximum neutron star masses and offer distinct radial oscillation signatures that could serve as probes for quark matter interiors.

The Gist

Imagine a neutron star as the ultimate cosmic pressure cooker. It's the collapsed core of a dead star, so dense that a single teaspoon of its material would weigh a billion tons on Earth. Inside this pressure cooker, nature is trying to figure out what happens when you squeeze matter so hard that the usual rules break down.

This paper is like a team of cosmic detectives trying to solve a mystery: What is the "recipe" for the stuff inside a neutron star?

Here is the breakdown of their investigation, explained in simple terms:

1. The Two Main Ingredients: The "Legos" and the "Soup"

Inside a neutron star, matter exists in two potential states:

• The Hadronic Phase (The Legos): At lower densities, the star is made of protons and neutrons. Think of these as individual Lego bricks stacked together. This is the "normal" state of matter we know.
• The Quark Phase (The Soup): If you squeeze hard enough, those Lego bricks might melt. The protons and neutrons break apart into their smaller components (quarks), turning into a hot, dense "soup" of free-floating particles.

The big question is: Does the transition happen like a sudden explosion (Legos shattering instantly), or is it a smooth melting process? The authors argue for the "smooth melting" (called a crossover).

2. The Detective Work: Using Clues from Everywhere

To figure out the recipe, the scientists couldn't just build a neutron star in a lab (it's too heavy!). Instead, they used clues from three different places:

• Clue A: The High-Density Theory (The "Far Future" Map):
  At incredibly high densities (densities we can't reach on Earth), a theory called pQCD (perturbative Quantum Chromodynamics) works like a reliable map. It tells us what the "soup" should look like if we squeeze it enough. The authors used this map to set a "speed limit" and "weight limit" for their recipe.
• Clue B: The Real-World Observations (The "Police Report"):
  Astronomers have been measuring real neutron stars. They know how heavy some are (some are twice as heavy as our Sun) and how big they are. If their recipe predicts a star that collapses under its own weight, the recipe is wrong.
• **Clue C: The "Speed Limit" of Sound (The "Traffic Rule"):
  In physics, nothing can travel faster than light. This means the "speed of sound" inside the star can't get too crazy. If their recipe makes sound travel too fast, it breaks the laws of physics.

3. The Big Discovery: The "Stiffness" Knob

The authors built a mathematical model that blends the "Legos" and the "Soup." They found that the transition isn't just about what the matter is, but how stiff it is.

Imagine the matter inside the star is like a mattress.

• A soft mattress (soft equation of state) squishes down easily. If the star is made of this, it might collapse into a black hole before it gets very heavy.
• A stiff mattress (stiff equation of state) resists squishing. This allows the star to hold up more weight.

The "Knob" they turned:
They found a specific "knob" in their math (called the vector coupling) that controls how stiff the quark soup is.

• If they turn the knob too low, the soup is too soft, and the star collapses.
• If they turn it too high, the soup gets so stiff it violates the "Far Future" map (pQCD).
• The Sweet Spot: They found a very narrow range where the knob is just right. This allows the star to be heavy enough to match the observations (like the 2-solar-mass stars) without breaking the laws of physics.

4. The "Melting" Effect

Here is the coolest part: When they added this "smooth melting" (crossover) to their models, it acted like a structural reinforcement.

• For stars that were already made of "stiff Legos," the melting didn't change much.
• But for stars made of "soft Legos," the transition to quark soup actually stiffened the whole structure. It was like adding steel beams to a weak building. This allowed the "soft" stars to suddenly support much more weight, solving a puzzle about how heavy neutron stars can get.

5. Listening to the Star's Heartbeat

Finally, the paper suggests a new way to test this theory. Neutron stars don't just sit there; they vibrate.

• Think of a bell. If you hit a bell made of soft rubber, it makes a low, dull sound. If you hit a bell made of hard steel, it rings high and clear.
• The authors calculated that if a neutron star has this "smooth melting" inside, its vibration frequency (its "ring") will change dramatically compared to a star made only of Legos.
• Specifically, as the star gets heavier, the "ring" might suddenly jump up or show a weird "double peak." This could be the "smoking gun" that tells astronomers, "Hey, there's quark soup in there!"

Summary

In short, this paper says:

1. Neutron stars likely contain a mix of normal matter and a "quark soup."
2. The transition between them is smooth, not sudden.
3. By using clues from high-energy physics and real telescope data, they found the exact "recipe" that keeps these stars from collapsing.
4. This recipe suggests that if we can listen to the "ringing" of neutron stars (via gravitational waves), we might finally hear the sound of quarks melting inside them.

It's like finally figuring out the secret recipe for a cosmic cake that is so dense it defies gravity, and realizing that the way it "rings" when you tap it could tell us exactly what's inside.

Technical Summary

1. Problem Statement

The internal composition of neutron stars (NSs) at densities exceeding nuclear saturation density ($n_0 \approx 0.15 \text{ fm}^{-3}$) remains a major open question in nuclear physics and astrophysics. Key challenges include:

• The Hadron-Quark Transition: It is unclear whether the transition from hadronic matter to deconfined quark matter is a sharp first-order phase transition (Maxwell/Gibbs construction) or a smooth crossover.
• The "Hyperon Puzzle": The inclusion of hyperons or other exotic degrees of freedom often softens the Equation of State (EOS), preventing the formation of massive NSs ($>2M_\odot$), which contradicts observational data.
• Model Constraints: While perturbative QCD (pQCD) is reliable at extremely high densities ($>40n_0$), it is difficult to connect these high-density constraints to the intermediate densities found in NS cores. Conversely, phenomenological models (like Relativistic Mean Field or NJL models) lack rigorous constraints from fundamental QCD at high densities.
• Parameter Ambiguity: In the Nambu–Jona-Lasinio (NJL) model used for quark matter, the vector coupling ($G_v$) and diquark coupling ($H$) are poorly constrained, leading to large uncertainties in the stiffness of the quark EOS.

2. Methodology

The authors construct a unified hadron-quark crossover EOS by combining three distinct theoretical frameworks and constraining them with multi-messenger data:

• Hadronic Phase (Low Density): Three Relativistic Mean-Field (RMF) models are employed to describe the hadronic phase ($n_B < 2n_0$):
  1. Nonlinear RMF (NLRMF) with the IUFSU parameter set.
  2. Density-dependent RMF (DDRMF) with the DDVTD parameter set.
  3. Density-dependent point-coupling (DDPC) with the DDPC1 parameter set.
• Quark Phase (High Density): The Nambu–Jona-Lasinio (NJL) model (three-flavor) is used for the quark phase ($n_B > 5n_0$). It includes:
  • Scalar four-quark interaction ($G_s$) for chiral symmetry breaking.
  • 't Hooft determinant interaction ($K$) for $U(1)_A$ anomaly.
  • Vector interaction ($G_v$) for repulsion.
  • Diquark interaction ($H$) for color superconductivity (Cooper pairing).
  • Density Dependence: To better satisfy pQCD constraints, the vector coupling $G_v$ is extended to be density-dependent ($G_v(n_q)$), introducing rearrangement terms to ensure thermodynamic consistency.
• Crossover Region ($2n_0 - 5n_0$): A smooth interpolation is performed in the pressure-chemical potential ($P-\mu_B$) space using a fifth-order polynomial. Continuity of pressure, baryon density, and compressibility is enforced at the boundaries.
• Constraints:
  1. pQCD Constraints: Using a scale-averaging likelihood approach, the EOS is matched to pQCD calculations (up to next-to-next-to-leading order) at high chemical potential ($\mu_H = 2.6$ GeV), assuming a Color-Flavor Locked (CFL) superconducting phase.
  2. Astrophysical Observations: Constraints from massive pulsars (PSR J0740+6620, PSR J0348+0432), radius measurements (PSR J0030+0451, PSR J0437–4715, PSR J0614–3329), and tidal deformability from GW170817.
  3. Causality: The speed of sound ($c_s^2$) must not exceed the speed of light ($c_s^2 \leq 1$).

3. Key Contributions

• Quantitative Coupling Constraints: The study provides tight constraints on the NJL model parameters. The diquark coupling is constrained to $H \simeq 1.5 G_s$, while the vector coupling is restricted to $G_v \lesssim 1.1 G_s$ when combining pQCD and astrophysical data.
• Density-Dependent Vector Coupling: The authors introduce a density-dependent vector coupling $G_v(n_q)$ in the NJL model. This modification allows the EOS to satisfy pQCD constraints more effectively and generates a distinct "stiffening" behavior in the crossover region compared to constant coupling models.
• Unified Framework: The paper successfully bridges the gap between low-density nuclear physics (RMF), intermediate-density crossover, and high-density pQCD, creating a self-consistent EOS that satisfies causality and thermodynamic stability.

4. Key Results

• Maximum Mass Enhancement: The hadron-quark crossover significantly enhances the maximum mass of neutron stars, particularly for softer hadronic EOSs (like DDVTD and IUFSU).
  • For the DDVTD set, the maximum mass increases by nearly 20% (from $1.85 M_\odot$ to $2.24 M_\odot$) when the crossover is included with optimized vector coupling.
  • This resolves the tension where pure hadronic models with these parameter sets fail to support $2 M_\odot$ stars.
• Speed of Sound ($c_s^2$): The crossover region exhibits a pronounced peak in the speed of sound.
  • For density-independent couplings, the peak shifts to lower densities as $G_v$ increases.
  • For density-dependent couplings ($G_v(n_q)$), $c_s^2$ reaches its maximum rapidly and then decreases, staying within the causal limit ($c_s^2 \leq 1$).
• Tidal Deformability ($\Lambda$): The inclusion of quark matter and the strength of the vector coupling significantly affect $\Lambda$.
  • Weak vector coupling reduces $\Lambda$, while strong coupling increases it.
  • All constructed EOSs satisfy the GW170817 constraint ($\Lambda_{1.4} \lesssim 800$).
• Radial Oscillation Frequencies ($\nu$):
  • Pure hadronic EOSs show a frequency peak around $0.5 M_\odot$.
  • Crossover EOSs shift the peak to higher masses (around $1.8 M_\odot$) and can exhibit a double-peak structure in the $\nu-M$ plane.
  • The fundamental radial oscillation frequency for intermediate-mass stars ($1.4 - 2.0 M_\odot$) is significantly higher in crossover models (increasing by 20–40%) compared to pure hadronic models.
• Trace Anomaly: The trace anomaly ($\Delta = 1/3 - P/\epsilon$) decreases rapidly in the crossover region for stiffer quark EOSs, approaching the conformal limit ($\Delta \to 0$) at high densities.

5. Significance

• Evidence for Quark Matter: The results suggest that the existence of massive neutron stars ($>2 M_\odot$) with relatively small radii (consistent with PSR J0614–3329) strongly favors a crossover transition to quark matter, which provides the necessary stiffening at intermediate densities.
• New Observational Probes: The study highlights that fundamental radial oscillation frequencies are a sensitive probe of the internal composition. The predicted double-peak structure and the rapid rise in frequency for intermediate-mass stars offer a potential signature to distinguish hybrid stars from pure neutron stars using future gravitational-wave detectors.
• Theoretical Robustness: By rigorously constraining the NJL parameters with pQCD likelihoods, the paper reduces the theoretical uncertainty in quark matter modeling, providing a more reliable framework for interpreting multi-messenger astronomical data.
• Future Directions: The framework sets the stage for studying finite-temperature effects, rotation, and non-radial modes, which are crucial for understanding binary mergers and post-merger remnants.

In conclusion, this work demonstrates that a hadron-quark crossover EOS, constrained by both high-energy QCD theory and precise astronomical observations, can simultaneously explain the existence of massive neutron stars and predict unique dynamical signatures (oscillation frequencies) that may soon be detectable.