Gravitational wave spectrum from first-order QCD phase transitions based on a parity doublet model

Using the parity doublet model, this paper demonstrates that gravitational waves from the nuclear liquid-gas phase transition in QCD could produce detectable signals in the millihertz to nanohertz bands, whereas those from the chiral phase transition are too weak to be observed, thereby linking the chiral invariant mass to potential gravitational wave signatures of nucleon mass origins.

The Gist

Imagine the early universe as a giant, super-hot pot of soup. For a tiny fraction of a second after the Big Bang, this soup was made of tiny, free-floating particles called quarks and gluons (the "quark-gluon plasma"). As the universe cooled down, this soup was supposed to "freeze" into the protons and neutrons that make up everything we see today (like stars, planets, and us).

Usually, physicists think this freezing happens smoothly, like water slowly turning into ice. But this paper asks: What if it happened like a sudden, violent explosion of bubbles instead?

Here is the story of that research, broken down into simple concepts:

1. The Two Types of "Freezing"

The researchers used a special mathematical recipe (called the Parity Doublet Model) to simulate this cooling process. They discovered that depending on how much "stuff" (density) was in the soup, the universe could freeze in two very different ways:

• The "Bubble Tea" Transition (Liquid-Gas): At lower densities, the universe might have undergone a transition similar to water boiling into steam or gas condensing into liquid. This creates big, energetic bubbles.
• The "Deep Freeze" Transition (Chiral): At extremely high densities, the particles change their internal structure (their "chirality"). This is like a deep freeze where the particles rearrange themselves completely.

2. The Cosmic "Pop" and the Sound

When a first-order phase transition happens, it's not quiet. Imagine a pot of water superheated past its boiling point. Suddenly, bubbles of steam erupt violently.

• The Bubbles: These bubbles of the new "frozen" state expand at nearly the speed of light.
• The Collision: When these bubbles crash into each other, they create a massive shockwave.
• The Ripples: Just as a stone thrown into a pond creates ripples, these cosmic bubble crashes create ripples in space and time called Gravitational Waves.

3. The Big Discovery: One Loud, One Whisper

The researchers calculated the "volume" (amplitude) and "pitch" (frequency) of these gravitational waves for both scenarios.

• The "Bubble Tea" Scenario (Liquid-Gas):

  • The Sound: This transition creates a loud, clear signal.
  • The Pitch: The waves are very low frequency (in the "nanohertz" range).
  • The Match: This signal fits perfectly with the mysterious "hum" that astronomers have recently detected using giant radio telescopes (Pulsar Timing Arrays). It suggests that the early universe might have had a violent "bubble tea" moment that we can actually hear today.

• The "Deep Freeze" Scenario (Chiral):

  • The Sound: This transition is incredibly quiet. The signal is about 100,000 times weaker than the first one.
  • The Pitch: Even if we could hear it, it's so faint that our current and future detectors are like trying to hear a whisper in a hurricane. It is effectively invisible to us right now.

4. The Secret Ingredient: The "Ghost Mass"

Why is there such a difference? The paper points to a mysterious property called Chiral Invariant Mass ($m_0$).

Think of a proton's mass like a cake.

• Part of the cake comes from the ingredients mixing together (spontaneous symmetry breaking).
• But this model says there is a "ghost ingredient" ($m_0$) that gives the proton its weight even before the mixing happens.

The researchers found that the size of this "ghost ingredient" dictates how loud the cosmic bubble pops are. By listening to the gravitational waves, we might be able to measure this "ghost mass," helping us solve the mystery of where the weight of matter actually comes from.

The Bottom Line

This paper suggests that the universe might have had a loud, explosive "bubble tea" moment in its infancy. If we listen closely with our gravitational wave detectors, we might finally hear that echo. This sound wouldn't just tell us about the early universe; it would act like a cosmic X-ray, revealing the fundamental building blocks of mass itself.

However, if the universe had a "deep freeze" moment at high densities, it was so quiet that we likely won't hear it with our current technology. We have to look for the loud "bubble tea" signals to learn the most.

Technical Summary

1. Problem Statement

The nature of the Quantum Chromodynamics (QCD) phase transition at finite baryon chemical potential ($\mu_B$) remains an open question due to the sign problem in lattice QCD calculations. While the transition at low density is a smooth crossover, it is hypothesized that a first-order phase transition may occur at high densities or specific conditions in the early Universe. Such transitions can generate a stochastic gravitational wave (GW) background.

However, previous studies using effective models (e.g., Friedberg-Lee, Quark-Meson) often predict transition rates ($\beta/H$) that are too large, pushing GW signals to frequencies detectable only by space-based interferometers (like LISA), or fail to account for the specific origin of nucleon mass. Furthermore, the connection between the chiral invariant mass ($m_0$)—the portion of nucleon mass persisting even when chiral symmetry is restored—and the resulting GW spectrum has not been systematically explored.

The authors aim to:

1. Investigate GW spectra from two distinct first-order QCD phase transitions: the nuclear liquid–gas transition and the chiral phase transition.
2. Utilize the Parity Doublet Model, which naturally incorporates $m_0$, to calculate the transition dynamics.
3. Determine if these signals are detectable by current or planned GW observatories (PTAs, LISA) and how they relate to the origin of nucleon mass.

2. Methodology

A. Theoretical Framework: Parity Doublet Model

The authors employ the Parity Doublet Model, where nucleons ($N$) and their parity partners ($N(1535)$) form a chiral doublet. The nucleon mass arises from two sources:

• Chiral Invariant Mass ($m_0$): A mass term that persists even when chiral symmetry is restored.
• Spontaneous Symmetry Breaking: A term proportional to the chiral condensate $\langle \bar{q}q \rangle$.

The thermodynamic potential $\Omega$ is constructed including mean fields for the scalar ($\sigma$) and vector ($\omega$) mesons. The model parameters are fixed by fitting to nuclear saturation properties ($n_0 = 0.16 \text{ fm}^{-3}$) and vacuum hadron masses. The study fixes $m_0 = 800 \text{ MeV}$, consistent with recent neutron star constraints (NICER, LIGO/Virgo).

B. Phase Diagram Analysis

The model predicts two distinct first-order phase transition regions at finite temperature ($T$) and chemical potential ($\mu_B$):

1. Liquid–Gas Transition: Occurs at low $\mu_B$ ($\sim 920 \text{ MeV}$).
2. Chiral Phase Transition: Occurs at high $\mu_B$ ($\sim 1635 \text{ MeV}$), associated with the excitation of the $N(1535)$ state.

C. Bubble Nucleation Dynamics

To calculate the GW spectrum, the authors solve the bounce equation (Euclidean action) to determine the nucleation rate of true vacuum bubbles within the metastable false vacuum.

• Action: The 3D Euclidean action $S_3/T$ is computed by solving the equation of motion for the $\sigma$ field profile $\sigma(r)$ under spherical symmetry.
• Nucleation Temperature ($T_n$): Defined where the nucleation rate $\Gamma \sim H^4$ (one bubble per Hubble volume per Hubble time).
• Transition Parameters:
  • Inverse Duration ($\beta/H$): Calculated as $T_n \frac{d(S_3/T)}{dT}$.
  • Transition Strength ($\alpha$): Defined carefully to account for the baryon number contribution at finite $\mu_B$:
    $$\alpha = \frac{-\Delta(\varepsilon - \mu_B n_B) + 3\Delta p}{4\varepsilon_r}$$
    where $\varepsilon_r$ is the radiation energy density. This subtraction is crucial as rest mass energy changes do not drive GWs.

D. Gravitational Wave Spectrum Calculation

The GW spectrum is computed considering three sources:

1. Sound Waves (SW): Dominant contribution from bulk fluid motion.
2. MHD Turbulence: Estimated as $\sim 5\%$ of the bulk kinetic energy.
3. Bubble Wall Collisions: Negligible for relativistic walls in this model.

The authors assume Jouguet detonation velocities and use numerical fits for efficiency factors ($\kappa_v$) to convert vacuum energy into kinetic energy.

3. Key Results

A. Phase Transition Characteristics

• Liquid–Gas Transition:

  • Near the endpoint of the first-order line ($\mu_B \approx 921.8 \text{ MeV}$), the transition strength $\alpha$ reaches $\mathcal{O}(1)$.
  • The inverse duration $\beta/H$ drops to $\mathcal{O}(10) - \mathcal{O}(100)$.
  • Result: This combination produces a strong GW signal with peak frequencies spanning the millihertz to nanohertz band.
  • Observation: The spectrum for $\mu_B \approx 921.8 \text{ MeV}$ provides a good fit to the NANOGrav 15-year data, suggesting a first-order nuclear liquid–gas transition could be the source of the observed stochastic GW background.

• Chiral Phase Transition:

  • Occurs at much higher $\mu_B$ ($\sim 1635 \text{ MeV}$).
  • The transition strength $\alpha$ is suppressed by approximately five orders of magnitude ($\alpha \sim 10^{-5}$) compared to the liquid–gas case.
  • Reason: At high $\mu_B$, the background energy density ($\varepsilon_r$) is dominated by the baryon rest mass energy. Since the latent heat is normalized by this large background, the relative transition strength $\alpha$ becomes vanishingly small.
  • Result: The resulting GW amplitude is well below the sensitivity of all current (PTAs) and planned (LISA, Taiji) detectors.

B. Dependence on Chiral Invariant Mass ($m_0$)

The study highlights that $m_0$ is a critical parameter. It determines the location and strength of the phase transitions. A different $m_0$ would shift the transition lines in the $T-\mu_B$ plane, altering the nucleation dynamics and the resulting GW spectrum.

4. Significance and Implications

1. Origin of Nucleon Mass: The paper establishes a novel link between the chiral invariant mass ($m_0$) and the gravitational wave spectrum. By observing (or constraining) GWs from the early Universe, one could potentially infer the value of $m_0$, offering a cosmological probe of the origin of hadron mass that complements neutron star observations and heavy-ion collisions.
2. Interpretation of NANOGrav Data: The results support the hypothesis that the stochastic GW background detected by PTAs could originate from a first-order nuclear liquid–gas transition in the early Universe, provided the initial baryon asymmetry was large and subsequently diluted (e.g., via "little inflation").
3. Detection Limits for High-Density QCD: The study reveals a fundamental limitation in GW astronomy for probing the high-density QCD phase diagram. Transitions occurring at very high baryon densities (like the chiral transition) produce signals too weak to be detected because the transition strength is diluted by the massive background energy density.
4. Degeneracy Breaking: While GW spectra alone cannot uniquely identify the microscopic order parameter (liquid-gas vs. chiral vs. deconfinement), the amplitude serves as a diagnostic. A strong signal in the nanohertz band favors low-density transitions, whereas the absence of a signal does not rule out high-density transitions.

5. Conclusion

The authors demonstrate that within the Parity Doublet Model, the nuclear liquid–gas transition is a viable source for the observed nanohertz gravitational wave background, while the chiral phase transition at high density is effectively invisible to current detectors due to extreme suppression of the transition strength. This work underscores the potential of multi-messenger astronomy (combining GWs, neutron star data, and collider physics) to constrain the chiral invariant mass and understand the phase structure of QCD.