Hadronic description of nuclear matter and neutron star properties

This paper demonstrates that a hadronic quantum hadrodynamics model incorporating $\sigma, \omega, \rho$, and $a_0$ mesons, along with their interactions, can simultaneously explain nuclear matter properties and astrophysical observations of neutron stars, predicting a peak in the speed of sound at intermediate densities that results in small intermediate-mass stars and a maximum mass of approximately $2M_\odot$.

The Gist

Imagine the universe as a giant cosmic kitchen. Inside this kitchen, the most extreme chefs are Neutron Stars. These are the dead, collapsed cores of massive stars, so dense that a single teaspoon of their material would weigh as much as a mountain.

For decades, physicists have been trying to figure out exactly what "ingredients" make up these cosmic mountains. Is it just a giant ball of squeezed protons and neutrons (like a super-dense Lego brick)? Or does the pressure get so high that the bricks melt into a soup of free-floating quarks (the fundamental particles inside protons and neutrons)?

This paper, written by Yao Ma, Yong-Liang Ma, and Jia-Ying Xiong, tackles this mystery with a fresh approach. Here is the story of their discovery, explained simply.

1. The Great Cosmic Tension

Imagine you are trying to bake a cake that satisfies two very picky judges:

• Judge A (The Nuclear Physicist): Says, "The cake must taste exactly like our experiments with heavy atoms on Earth."
• Judge B (The Astronomer): Says, "The cake must be big enough to weigh 2 Suns, but small enough to fit inside a city."

For a long time, these two judges were fighting. The recipes (theories) that satisfied Judge A made cakes that were too big for Judge B. The recipes that fit Judge B were too light for Judge A. Scientists thought, "Maybe we need a secret ingredient, like a 'quark soup' or some exotic matter, to fix this."

2. The "Universal Recipe" (The GQHD Model)

The authors decided to test a bold idea: What if we don't need any exotic secret ingredients? What if the standard "Lego bricks" (protons and neutrons) plus a few standard "glue" particles (mesons) are enough?

They built the most comprehensive "recipe book" possible, called the General Quantum Hadrodynamics (GQHD) model. Think of this as a master chef adding every possible spice and interaction they know of (sigma, omega, rho, and a0 particles) to the mix.

They didn't just guess; they used a Bayesian Joint Analysis. Imagine a super-smart computer that tastes thousands of different cake variations at once, checking them against every piece of data we have from Earth labs and space telescopes simultaneously. It's like a blind taste test where the computer scores every recipe on a scale of 1 to 100 based on how well it fits reality.

3. The Surprise Discovery

The computer found something amazing: The "Standard Lego" recipe actually works!

They discovered that you don't need to melt the bricks into quark soup to explain the universe. A universe made entirely of hadrons (protons, neutrons, and their interactions) can satisfy both judges perfectly.

The Secret Sauce:
The magic happened because of a specific interaction between four particles (sigma, omega, rho, and a0). The authors call this the $\sigma\omega\rho a_0$ interaction.

To use an analogy: Imagine the neutrons in the star are people in a crowded elevator.

• At first, they are polite and keep their distance (low density).
• As the elevator gets packed (medium density), they start to push against each other, but then they suddenly find a way to "hug" or interact in a specific, complex way that changes the pressure.
• This creates a peak in the "Speed of Sound."

What is the Speed of Sound here?
In a star, the "speed of sound" isn't about noise; it's a measure of stiffness.

• Soft EoS (Equation of State): The star is squishy like a marshmallow. It's easy to compress.
• Stiff EoS: The star is hard like a rock. It resists compression.

The authors found that this special interaction makes the star behave like a hybrid material:

1. In the middle: It gets "squishy" (soft). This allows medium-sized neutron stars to shrink to a smaller radius, matching the observations of stars like PSR J0614-3329.
2. At the very core: It gets "rock hard" (stiff). This prevents the heaviest stars from collapsing into black holes, allowing them to reach the massive 2-sun weight limit.

4. Why This Matters

This discovery changes the game in two ways:

1. Simplicity Wins: It suggests we might not need to invent "exotic" physics (like quark matter) to explain neutron stars. The standard physics of protons and neutrons, if we just understand their complex interactions better, might be enough.
2. The "Goldilocks" Test: The paper argues that the key to solving the mystery is measuring medium-sized neutron stars.
   • If we find a medium-sized star that is too big, it means the "squishy" middle part isn't working, and maybe we do need exotic matter.
   • If it fits the "small radius" prediction, it confirms that our "Standard Lego" recipe is correct.

The Bottom Line

The authors are saying: "Don't look for aliens in the kitchen yet. The ingredients we already have might be enough to bake the perfect cosmic cake."

They found that the complex dance between standard particles creates a "sweet spot" in the star's density that naturally explains why some stars are small and others are massive. The next step? We need better telescopes to measure the size of those medium-sized stars to confirm if this "Standard Lego" theory is the final answer.

Technical Summary

1. Problem Statement

The composition of dense nuclear matter (NM) inside neutron stars (NSs) remains a fundamental open question in nuclear and astrophysics. While massive neutron stars ($\sim 2M_\odot$) suggest a stiff Equation of State (EoS) with a high speed of sound, recent measurements of intermediate-mass pulsars (e.g., PSR J0614-3329) indicate smaller radii, implying a softer EoS at intermediate densities. This creates a "tension" where theoretical models struggle to simultaneously satisfy:

• Nuclear constraints: Experimental data on nuclear matter properties around saturation density ($n_0 \approx 0.16 \text{ fm}^{-3}$).
• Astrophysical constraints: Mass-radius ($M-R$) relations of massive NSs and intermediate-mass NSs, as well as tidal deformability from gravitational wave events (GW170817).

Previous approaches often invoke exotic degrees of freedom (DoFs) such as quark matter, quarkyonic matter, or color-flavor-locked phases to resolve this tension. This paper investigates whether the simplest hadronic description (pure nucleons and mesons) can satisfy all constraints without resorting to exotic phases.

2. Methodology

The authors employ a General Quantum Hadrodynamics (GQHD) model combined with a Bayesian Joint Analysis (BJA).

• Theoretical Model (GQHD):

  • The Lagrangian includes nucleons ($N$) and all meson resonances below 1 GeV: scalar $\sigma$, vector $\omega$, isovector $\rho$, and the scalar isovector $a_0(980)$ (often denoted as $\delta$).
  • The model is constructed up to dimension-4, incorporating self-interactions of mesons and cross-interactions (e.g., $\sigma\omega$, $\sigma\rho$, $\omega\rho$, and crucially, the mixed $\sigma\omega\rho a_0$ interaction).
  • The pion is omitted as it vanishes in the Relativistic Mean Field (RMF) approximation.
  • The energy density and pressure are calculated using the standard RMF approach, and the NS structure is derived by solving the Tolman-Oppenheimer-Volkoff (TOV) equations.

• Statistical Framework (Bayesian Joint Analysis):

  • The authors utilize Bayes' theorem to constrain the 21-dimensional parameter space (coupling constants and meson masses) simultaneously against two datasets:
    1. Nuclear Matter Properties: Binding energy, pressure, incompressibility, symmetry energy, and its slope at and near saturation density.
    2. Astrophysical Observations: $M-R$ constraints from massive pulsars (PSR J0740+6620, etc.), intermediate pulsars (PSR J0614-3329), and tidal deformability ($\Lambda_{1.4}$) from GW170817.
  • A normalized likelihood function is used to weight observables with different magnitudes, allowing for a unified statistical assessment at the 1$\sigma$ level.

3. Key Contributions

• Unified Hadronic Description: The study demonstrates that a pure hadronic model, without invoking phase transitions to quark matter, can simultaneously saturate constraints from nuclear experiments and diverse astrophysical observations.
• Role of the $a_0$ Meson and Mixed Interactions: The inclusion of the $a_0(980)$ meson and, specifically, the $\sigma\omega\rho a_0$ mixed interaction term (parameter $b_8$) is identified as the critical mechanism enabling the model to reproduce the required EoS behavior.
• Sound Velocity Peak: The model naturally predicts a peak structure in the speed of sound ($v_s$) at densities of $\sim (2-3)n_0$. This non-monotonic behavior allows the EoS to be soft at intermediate densities (yielding small radii for intermediate-mass stars) and stiff at high densities (supporting massive stars).
• Systematic Comparison: The work provides a rigorous comparison between the proposed GQHD model and established Walecka-type models (TM1, NL1, NL3, FSUGold, etc.), highlighting the limitations of simpler models in satisfying the full set of constraints.

4. Key Results

• Model Performance:

  • The optimized GQHD parameter sets (GQHD1 and GQHD2) successfully reproduce nuclear matter properties (binding energy, symmetry energy, etc.) within experimental uncertainties.
  • Unlike standard Walecka models (e.g., NL1, NL3) which predict radii too large for intermediate-mass pulsars, the GQHD model predicts radii consistent with PSR J0614-3329 while maintaining a maximum mass of $\sim 2M_\odot$.
  • The Bayesian Joint Analysis (BJA) score for GQHD2 is significantly higher than for other models, indicating superior overall agreement with data.

• Equation of State (EoS) Characteristics:

  • The EoS derived from GQHD is soft at intermediate densities (due to attractive forces enhanced by the mixed interactions) but becomes stiff at high densities.
  • This stiffness transition is driven by the peak in the speed of sound ($v_s^2 = dP/d\varepsilon$), which reaches values consistent with causality and massive star observations.

• Tidal Deformability:

  • The model predicts a tidal deformability for a $1.4M_\odot$ star of $\Lambda_{1.4} \approx 320-440$.
  • While this falls within the broad constraints of GW170817, it is noted that the predicted values are on the lower end, suggesting that future, more precise measurements of tidal deformability could provide even stricter constraints on hadronic models.

• Mass-Radius Relations:

  • The model predicts that intermediate-mass neutron stars ($\sim 1.3-1.5 M_\odot$) have smaller radii compared to massive stars, a feature directly linked to the peak in the speed of sound.

5. Significance and Outlook

• Redefining the Necessity of Exotic Matter: The paper argues that there is currently no compelling motivation to introduce exotic degrees of freedom (like deconfined quarks) to explain current neutron star data. The "simplest" hadronic approach, when extended with the most general mesonic interactions, is sufficient.
• Diagnostic Power of Intermediate Mass Stars: The authors emphasize that the sequential measurement of mass and radius, particularly for intermediate-mass neutron stars, is the most crucial future test. These stars are sensitive to the "soft" part of the EoS and the speed of sound peak, serving as a discriminator between pure nucleonic stars and hybrid (exotic) stars.
• Future Directions:
  • The current model is phenomenological. Future work should incorporate fundamental QCD symmetries, specifically chiral symmetry, to derive the couplings from first principles rather than fitting them.
  • The study sets a benchmark for precision: distinguishing between pure hadronic and hybrid structures will require next-generation facilities to measure the mass-radius relation with higher precision, specifically targeting the intermediate mass range.

In conclusion, this work provides a robust theoretical framework showing that the complex interplay of standard hadronic resonances, particularly the $a_0$ meson and its mixed interactions, can naturally resolve the tension between nuclear and astrophysical data, offering a unified description of neutron star matter without requiring exotic phase transitions.