NJL-Chiral Soliton and the Nucleon Equation of State at supra-saturation density: Impact of Chiral Symmetry Restoration

This paper proposes a method to construct the equation of state for supra-saturation density nuclear matter by modeling nucleons as topological solitons within a bosonized Nambu-Jona-Lasinio framework, demonstrating that the dynamical restoration of chiral symmetry leads to a stiffened equation of state compatible with neutron star observations.

The Gist

The Big Picture: What is this paper about?

Imagine the universe is made of two main types of "stuff": the stuff that makes up stars and planets (nuclear matter), and the stuff that makes up the very core of an atom (quarks). Usually, scientists think these are two completely different worlds.

This paper asks a fascinating question: What if the inside of a single proton (a nucleon) is actually a tiny, super-dense star?

The authors propose that if you squeeze matter hard enough—like the pressure found inside a neutron star—you don't need to invent new physics. You just need to look at what's already happening inside a single proton. They built a mathematical model to see how protons behave when they are squeezed together, and they found that this behavior perfectly explains how neutron stars hold themselves together.

---

The Main Characters and Tools

To understand their story, let's meet the cast of characters:

1. The Proton (The Nucleon): Think of a proton not as a solid marble, but as a soft, fluffy cloud surrounding a hard, dense rock.
   • The cloud is made of pions (particles that carry the "glue" of the strong force).
   • The rock is the "hard core" where the three main quarks live.
2. The NJL Model (The Blueprint): This is the set of rules the authors use to describe how quarks and gluons interact. It's like the "instruction manual" for building a proton.
3. Chiral Symmetry (The Elastic Band): This is a fancy physics term for a rule that keeps the proton's structure balanced. Imagine a rubber band holding the proton together. When things are normal, the band is tight. But when you squeeze the proton (increase density), the band starts to loosen. This "loosening" is called Chiral Symmetry Restoration.
4. The Soliton (The Shape-Shifter): In this paper, the proton is treated as a "soliton." Think of a soliton as a wave in a pond that keeps its shape even as it moves. The proton is a stable wave of energy that doesn't just fall apart.

---

The Story: Squeezing the Proton

1. The "Hard Core" Idea

The authors start with a bold idea: The equation of state (EoS) of a neutron star is just the pressure inside a single proton.

Imagine you have a giant room full of people (nuclear matter). If you squeeze them together, they push back. The authors say, "Let's ignore the people for a moment and just look at the pressure inside one person's chest." They calculated the pressure and energy density right in the center of a proton and found that it matches the pressure needed to support a massive neutron star.

2. The "Fluffy Cloud" vs. The "Hard Rock"

In normal conditions, a proton has a hard core (the rock) and a soft cloud (the pions).

• The Cloud: This part is very sensitive to the environment. When you squeeze protons together, this cloud gets squished and changes shape.
• The Rock: The hard core is tough. It doesn't change much until the pressure is extreme.

The authors used their model to see what happens when you increase the density. They found that as you squeeze the protons, the "rubber band" (chiral symmetry) starts to loosen. This causes the "hard rock" inside the proton to swell up slightly.

3. The "Swelling" Effect

Here is the magic trick:

• As the density increases, the proton's hard core gets bigger (it swells).
• Because the cores are swelling, they start to bump into each other sooner than we thought.
• When they bump, the "soft" parts of the protons merge, and the "hard" cores start to overlap.

This overlap is the key. The authors found that when the cores overlap, the matter becomes incredibly stiff (resistant to compression). This stiffness is exactly what is needed to prevent a neutron star from collapsing into a black hole.

---

The "Aha!" Moment: Why This Matters

For a long time, scientists had a problem. They knew that at very high densities, protons should turn into a soup of free quarks (quark-gluon plasma). But they didn't know how that transition happened or what the pressure looked like during the change.

This paper provides a bridge:

• Old View: Protons are solid; then suddenly, they explode into quark soup.
• New View (This Paper): Protons are like balloons. As you squeeze them, they swell, their internal structure changes, and they gradually merge into a giant, unified blob of quark matter.

The authors showed that by accounting for the "loosening of the rubber band" (chiral symmetry restoration), their model predicts a stiff equation of state. This means the matter is very hard to compress, which matches real-world observations of massive neutron stars (like those weighing twice as much as our Sun).

The Analogy of the "Crowded Dance Floor"

Imagine a crowded dance floor where everyone is dancing in a circle (the proton).

• Normal Density: Everyone has their own space. The "cloud" of their arms (the pion cloud) is waving around, but their bodies (the core) aren't touching.
• High Density: The DJ turns up the music and the room gets smaller. Everyone is pushed closer.
• The Result: The "arms" (clouds) get squished and disappear. The "bodies" (cores) start to touch.
• The Twist: In this paper, the authors realized that as the bodies touch, they don't just bounce off; they actually expand slightly because the rules of the dance (chiral symmetry) change. This expansion makes the crowd push back harder against the walls of the room.

Conclusion

In simple terms, this paper says: Neutron stars are just giant collections of protons that have been squeezed so hard that their internal structures have changed.

By understanding how a single proton changes when it's under extreme pressure (specifically, how its internal "symmetry" restores itself), the authors could predict how neutron stars behave. They found that this internal change makes the star's core incredibly strong, explaining why these stars don't collapse under their own weight. It's a beautiful example of how understanding the smallest building blocks of the universe helps us understand the largest objects in the sky.

Technical Summary

1. Problem Statement

The paper addresses the challenge of describing the Equation of State (EoS) of dense nuclear matter at supra-saturation densities (densities exceeding nuclear saturation density, $\rho_{sat} \approx 0.16 \text{ fm}^{-3}$), such as those found in the cores of neutron stars.

• The Gap: While perturbative QCD works at very high energies and Lattice QCD works at high temperatures/low densities, the intermediate high-density/low-temperature regime is difficult to access due to the "sign problem" in Lattice QCD.
• The Hypothesis: The authors investigate the conjecture that at sufficiently high baryon densities, the thermodynamic EoS of bulk nuclear matter can be identified with the mechanical EoS of the nucleon core.
• The Specific Challenge: Previous models (e.g., Ref. [51]) treated the nucleon as a topological soliton (Skyrmion) but approximated chiral symmetry restoration effects via simple rescaling of parameters. This paper aims to provide a self-consistent, dynamical treatment of chiral symmetry restoration within the nucleon structure itself to determine its impact on the EoS.

2. Methodology

The authors employ a framework combining the Nambu–Jona-Lasinio (NJL) model with Chiral Soliton (Skyrmion) theory.

A. Theoretical Framework

1. Bosonized NJL Model: The underlying quark-level NJL Lagrangian is bosonized using path-integral techniques. This generates an effective mesonic theory containing scalar ($S$), vector ($\omega, \rho$), and pseudoscalar ($\pi$) fields.
2. Chiral Soliton Ansatz: Nucleons are modeled as topological solitons stabilized by vector mesons ($\omega$ and $\rho$), avoiding the need for the ad-hoc Skyrme term. The soliton is described using the hedgehog ansatz for the pion field and specific ansatzes for the vector mesons (Wu-Yang-'t Hooft-Polyakov for $\rho$, time-component only for $\omega$).
3. In-Medium Dynamics:
   • Chiral symmetry restoration is implemented via a density-dependent scalar field $S$.
   • The scalar field $S$ acts as an order parameter related to the constituent quark mass. As density increases, $S$ decreases (approaching zero), signaling the restoration of chiral symmetry.
   • All meson parameters (pion decay constant $F_\pi$, meson masses $M_\pi, M_\sigma, M_v$, and couplings $g_v$) become functionals of the scalar field $S$.
4. Numerical Implementation:
   • The authors solve coupled non-linear differential equations for the meson profiles ($F, G, \omega$) using a relaxation method on a compactified coordinate system ($t = \tanh(r)$). This overcomes the instability and sensitivity issues of the traditional shooting method.
   • A phenomenological scaling factor $\chi$ is applied to match the soliton mass to the physical nucleon mass ($\sim 939$ MeV), correcting for the lack of explicit spin-isospin quantization in the static soliton approximation.

B. Construction of the Bulk EoS

To derive the bulk EoS from the single-nucleon profile:

1. Core Identification: The "hard core" of the nucleon is identified with the valence quark distribution.
2. Mapping: For a given bulk density $\rho_N$, the corresponding scalar field $S(\rho_N)$ is determined.
3. Extraction: The authors map the bulk energy density to a specific radial coordinate $r_{QC}$ inside the soliton where the local energy density matches the bulk average. The pressure at this radius is then extracted to form the $P(\epsilon)$ relation.

3. Key Contributions

1. Self-Consistent Chiral Restoration: Unlike previous works that rescaled parameters post-hoc, this study dynamically evolves the soliton structure (mass, radii, profiles) as the scalar field $S$ changes, capturing the interplay between scalar and vector sectors.
2. Numerical Robustness: The implementation of the relaxation method with coordinate compactification provides a stable and accurate solution for the stiff differential equations governing the soliton profiles, validated against existing literature.
3. Microscopic Derivation of Stiffening: The paper demonstrates that the progressive restoration of chiral symmetry naturally leads to a stiffening of the nucleon-based EoS, making it compatible with astrophysical constraints (e.g., $2M_\odot$ neutron stars).
4. Hard Deconfinement Threshold: The study identifies specific density thresholds where soliton cores overlap, signaling the transition from hadronic matter to quark matter ("hard deconfinement").

4. Key Results

A. Vacuum Properties and Parameter Sensitivity

• The model reproduces vacuum soliton properties (mass $\sim 1450$ MeV, baryon charge radius $\sim 0.5$ fm) consistent with Ref. [51].
• Vector Meson Mass ($M_v$): Increasing $M_v$ localizes the $\omega$-field, reducing the isoscalar charge radius but leaving the baryon charge radius (hard core) relatively robust.
• Scalar Coupling: The nucleon-scalar coupling $g_S$ is found to be small, and the scalar susceptibility is negative, a feature attributed to the lack of explicit confinement in the current model.

B. Impact of Chiral Symmetry Restoration

• Soliton Swelling: As chiral symmetry restores (scalar field $s$ increases), the quark distribution delocalizes. The baryon charge radius (hard core) increases, causing nucleon cores to overlap at lower densities than predicted by vacuum radii.
• EoS Stiffening:
  • The in-medium EoS becomes significantly stiffer as chiral symmetry restores.
  • With the scaling factor $\chi$ (matching physical mass), the EoS reaches the stiffness required by neutron star observations at scalar fields $s \sim 0.2-0.3 f_\pi$.
  • Without scaling, stiffness is reached at higher fields ($s \sim 0.5 f_\pi$).
  • This stiffening is a result of the competition between the scalar sector (which tends to soften) and the vector sector (which stiffens), with the self-consistent treatment capturing the net effect correctly.

C. Density-Dependent EoS and Phase Transition

• Overlap Density: The soliton-based EoS remains valid until the hard cores overlap. This occurs at:
  • $\rho_N \approx 8-10 \rho_{sat}$ for different NJL parameter sets.
  • This overlap coincides with the "breakdown" of the soliton description, interpreted as the onset of hard deconfinement (transition to quark matter).
• Comparison with NS Models: The derived EoS ($p_{QC, \chi}$ vs $\epsilon_{QC, \chi}$) aligns well with established high-density neutron star models (SLy4, QHC18) and satisfies the $2M_\odot$ mass constraint.
• Speed of Sound: The squared speed of sound ($c_s^2$) in the derived EoS remains consistent with causality and observational constraints.

5. Significance

• Unified Description: The work provides a microscopic bridge between the internal structure of a single nucleon and the macroscopic EoS of neutron stars, validating the hypothesis that the nucleon core EoS dictates bulk matter behavior at high densities.
• Chiral Dynamics: It demonstrates that chiral symmetry restoration is not merely a background effect but actively modifies the mechanical properties (pressure and energy density) of the nucleon, leading to a stiffer EoS that supports massive neutron stars.
• Phase Transition Mechanism: It offers a physically transparent mechanism for the hadron-to-quark phase transition: as density increases, chiral symmetry restoration swells the nucleons until their hard cores overlap, triggering a percolation-like transition to a deconfined quark phase.
• Future Directions: The authors note that future work should include axial mesons, explicit confinement potentials, and spin-isospin quantization to refine the scalar susceptibility and the precise density of the phase transition.

In summary, this paper successfully constructs a self-consistent NJL-chiral soliton model that dynamically incorporates chiral symmetry restoration, demonstrating that this mechanism naturally produces a stiff Equation of State compatible with modern neutron star observations and provides a clear microscopic picture of the onset of quark matter.