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

We investigate the gravitational wave spectrum from first-order QCD phase transitions using the parity doublet model at finite baryon chemical potential. The model incorporates the chiral invariant mass $m_0$, representing the portion of nucleon mass that persists even when chiral symmetry is restored. Within the model, we identify two first-order phase transition regions: the nuclear liquid–gas transition and the chiral phase transition. By solving the bounce equation and computing the Euclidean action $S_3/T$, we obtain the gravitational wave spectra from both transitions. The liquid–gas transition yields $α\sim \mathcal{O}(1)$ and $β/H \sim \mathcal{O}(10)$–$\mathcal{O}(100)$ near the endpoint of the first-order line, producing signals with peak frequencies from the millihertz to the nanohertz band that can fit the existing data. In contrast, the chiral transition produces signals suppressed by approximately five orders of magnitude, well below the sensitivity of all current and planned detectors. These results connect the chiral invariant mass to the gravitational wave spectrum, offering a novel probe of the origin of nucleon mass through gravitational wave astronomy.

💡 Research Summary

The authors investigate the stochastic gravitational‑wave (GW) background that would be produced by first‑order QCD phase transitions in the early Universe, employing the parity‑doublet model of nuclear matter at finite baryon chemical potential. In this effective theory nucleons and their opposite‑parity partner N(1535) form a chiral doublet; the nucleon mass receives two contributions: a chiral‑invariant piece m₀ (which remains non‑zero when chiral symmetry is restored) and a component proportional to the chiral condensate ⟨ q̄q ⟩. By fitting the model parameters to nuclear saturation properties and fixing m₀≈800 MeV (consistent with neutron‑star constraints), the authors map out the low‑temperature phase diagram. Two distinct first‑order transition lines appear: (i) a low‑density nuclear liquid–gas transition and (ii) a high‑density chiral transition associated with the excitation of the parity partner.

The paper proceeds to compute the bubble nucleation dynamics using homogeneous thermal nucleation theory. The three‑dimensional Euclidean action S₃/T is obtained by solving the O(3) bounce equation for the σ field, which serves as the order parameter. Numerical solutions for the σ profile are presented for several values of the baryon chemical potential μ_B on both sides of each transition line. The action exhibits a characteristic “U‑shaped’’ temperature dependence: it is large at low T (high barrier), reaches a minimum at an intermediate temperature where thermal fluctuations most efficiently overcome the barrier, and diverges again as T approaches the critical temperature because the critical bubble radius diverges.

From the temperature dependence of S₃/T the authors extract the nucleation rate Γ∝exp(−S₃/T) and the dimensionless transition rate β/H. For the liquid–gas transition they find α≈𝒪(1) (the ratio of vacuum energy released to the radiation energy density) and β/H≈10–100 near the endpoint of the first‑order line. These values imply relatively slow transitions that generate strong GW signals. Using the standard formulas for bubble‑collision, sound‑wave, and magnetohydrodynamic turbulence contributions, the resulting GW spectrum peaks in the millihertz–nanohertz band, with amplitudes capable of fitting the stochastic background reported by recent pulsar‑timing‑array (PTA) collaborations (NANOGrav, EPTA, PPTA, etc.). Thus, a QCD liquid–gas transition occurring after a period of large baryon asymmetry (e.g., Affleck‑Dine baryogenesis followed by a brief low‑scale inflation) could be a viable source of the observed PTA signal.

In contrast, the chiral transition yields a much larger action, leading to β/H≈10⁴–10⁵ and α≈10⁻³. Consequently the GW power is suppressed by roughly five orders of magnitude relative to the liquid–gas case. The peak frequencies lie in the millihertz range, but the signal lies far below the sensitivity of planned space‑based interferometers such as LISA, Taiji, or DECIGO, as well as ground‑based detectors.

A key conceptual outcome is the explicit link between the chiral‑invariant mass m₀ and GW observables. Larger m₀ pushes the chiral transition to higher μ_B and makes the barrier steeper, which dramatically reduces the nucleation rate and weakens the GW signal. Therefore, a detection (or non‑detection) of a low‑frequency stochastic GW background could, in principle, constrain m₀ and shed light on the portion of the nucleon mass that does not arise from spontaneous chiral symmetry breaking—a long‑standing question in hadron physics.

The paper also discusses the broader relevance of the parity‑doublet framework: it has been successfully applied to neutron‑star phenomenology, and the same parameter set reproduces both astrophysical constraints and the QCD phase structure needed for GW calculations. By integrating nuclear‑matter modeling, finite‑temperature field theory, and cosmological GW production, the work provides a coherent bridge between microscopic QCD dynamics and macroscopic gravitational‑wave astronomy, opening a novel avenue to probe the origin of nucleon mass through upcoming GW observations.