First-order pion condensation at the QCD epoch: tricritical behavior at high isospin chemical potential and the primordial gravitational-wave imprint

Whether the early universe, under a large but observationally allowed lepton flavor asymmetry, entered a pion-condensation phase through a first-order phase transition during the QCD epoch — and what gravitational-wave and cosmological imprint such a transition would leave — remains unestablished. In the two-flavor quark-meson model at mean field, the pion-condensation boundary contains a tricritical point, so that for $|l_e + l_\mu| \gtrsim 0.1$ the cosmic trajectory enters the pion-condensed phase through a first-order transition; but this feature is model-dependent (it is absent in a renormalization-group-invariant mean-field formulation) and is unconfirmed by lattice QCD at high isospin chemical potential. The detailed gravitational-wave prediction of the first-order transition — bubble-nucleation history, the inverse-duration parameter $\beta/H$, latent-heat partition, spectral templates, and distinguishability from a quark-hadron transition in pulsar timing array data — has been explicitly deferred to future work, most recently by a March 2026 follow-up that computed only the complementary speed-of-sound imprint on the gravitational-wave background.

Background

The total lepton asymmetry of the universe is tightly constrained by the cosmic microwave background and Big Bang Nucleosynthesis (BBN), but the asymmetries of individual lepton flavors, $l_f$ (the difference between the number densities of leptons and antileptons of flavor $f$, including the corresponding neutrinos, normalized to the entropy density), can in principle be large, $|l_f| \sim 0.01$–$1$, if they nearly cancel in the sum. Because the primordial plasma is electrically neutral and baryon and lepton numbers are conserved back to the QCD epoch, a large imbalance $l_e + l_\mu$ induces a large electric-charge chemical potential $\mu_Q$ at temperatures $T \sim 150$ MeV. For two massless light-quark flavors the charge chemical potential maps to an isospin chemical potential $\mu_I \simeq \mu_Q/2$; when $\mu_I$ exceeds half the pion mass, $m_\pi/2$, the lowest charged-pion mode softens and a Bose-Einstein condensate of $\pi^\pm$ forms ("pion condensation" or pion superfluidity). This onset is a second-order transition and is firmly established by chiral perturbation theory (the effective field theory of low-energy pions) and by lattice quantum chromodynamics (QCD, the first-principles numerical simulation of the strong interaction on a discretized spacetime grid). What is not established is the order of the transition deeper in the condensed phase at high $\mu_I$, where effective models disagree. The quark-meson (QM) model — a chiral effective theory of the $\sigma$ field and pions coupled to quarks, here treated at mean field with the no-sea approximation (Dirac-sea/vacuum contributions to the thermodynamic potential neglected) — predicts that the second-order pion-condensation line turns first order at a tricritical point, located by Ferreira et al. at $T_{tc} \approx 124$ MeV, $\mu_{Qc} \approx 186$ MeV for $\mu_B = 100$ MeV. The Polyakov loop (an order-parameter-like variable tracing deconfinement, used to extend the model to Polyakov-quark-meson type), three-flavor extensions, mesonic fluctuations beyond mean field, and renormalization-group-invariant (RGI) treatments of the vacuum terms all shift or may remove this tricritical point: an RGI mean-field formulation finds no tricritical point at the $\mu_Q$ values investigated. On the lattice, the pion-condensation boundary is confirmed second order at least up to $\mu_I \approx 104.5$ MeV, and large-isospin-density ground states have been simulated at unphysical pion mass, but the existence of a first-order segment at higher $\mu_I$ has neither been confirmed nor excluded. Tracing the expanding universe through the $(T, \mu_B, \mu_Q)$ phase diagram (the "cosmic trajectory", fixed by conservation of baryon number, lepton flavor numbers, and charge neutrality), Vovchenko et al. showed with a lattice-anchored hadron-resonance-gas equation of state that for $|l_e + l_\mu| \gtrsim 0.1$ the trajectory crosses the pion-condensed phase, and Cao, He, and Zhang found reentrant pion superfluidity along trajectories computed in the Polyakov–Nambu–Jona-Lasinio (PNJL) model, a chiral quark model with a local four-quark (Nambu–Jona-Lasinio, NJL) interaction supplemented by the Polyakov loop. Ferreira et al. then showed that in the QM model the entry can be first order (via the tricritical point) with exit second order, making pion condensation a candidate source of primordial gravitational waves (GWs): a first-order transition proceeds by bubble nucleation, and the resulting GW spectrum is controlled by the latent heat $L$ released (they estimate $L \approx 2000$ MeV/fm$^3$ at $|l_e + l_\mu| = 0.2$), the inverse transition duration $\beta/H$ (bubble-nucleation rate versus Hubble expansion), and the partition of latent heat into bulk fluid motion and reheating. At $T \sim 150$ MeV the signal peaks in the nanohertz range, where pulsar timing arrays (PTAs) — collaborations timing millisecond pulsars to detect a stochastic GW background, e.g. the NANOGrav (North American Nanohertz Observatory for Gravitational Waves) 15-year signal — are most sensitive. The same physics leaves other fingerprints: the speed of sound exceeds the conformal bound $c_s^2 \leq 1/3$ inside the condensed phase (seen also on the lattice at high isospin density), and the altered equation of state can shift weakly interacting massive particle (WIMP) dark-matter freeze-out, the primordial black hole (PBH) mass spectrum, and possibly nucleate compact self-bound pion stars.

Problem Statement

As posed by Ferreira et al. (Phys. Rev. D 112, 094009, 2025), the problem is to establish the quantitative robustness and the observational consequences of pion condensation driven by large lepton flavor asymmetry during the QCD epoch. Its components jointly serve one objective — to determine whether the universe actually entered the pion-condensed phase through a first-order phase transition and what its observable imprint would be — and are linked as follows. (1) Transition order and the tricritical point (the existence condition for everything downstream): determine whether the pion-condensation phase boundary contains a first-order segment emanating from a tricritical point at high charge/isospin chemical potential once the systematics of the two-flavor mean-field no-sea quark-meson treatment are removed — under the three-flavor extension, Polyakov-loop coupling, inclusion of vacuum terms (renormalization-group-invariant formulation), and mesonic fluctuations beyond mean field — and obtain a lattice-QCD confirmation or exclusion of such a first-order segment at large $\mu_I$ (at $\mu_B = 0$, where the fermion determinant is real and positive, so that the sign problem — a complex determinant that blocks Monte Carlo importance sampling — does not obstruct the simulation), including the effect of small finite $\mu_B$ along realistic cosmic trajectories. (2) Gravitational-wave phenomenology (the observable prediction, conditional on the first-order entry found in component 1): compute the complete GW signal of the first-order pion-condensation transition along the cosmic trajectory — the bubble-nucleation history, the inverse-duration parameter $\beta/H$, the fraction of latent heat partitioned into bulk fluid motion versus reheating, the resulting spectral templates $\Omega_{GW}(f)$, their distinguishability from a first-order quark-hadron transition and from other nanohertz sources, and their quantitative compatibility with pulsar timing array data such as the NANOGrav 15-year signal. (3) Remaining observational consequences: quantify the effects of the pion-condensation equation of state and of the first-order transition on the relic abundance of WIMP dark matter (freeze-out), on the primordial black hole mass spectrum, and on the possible formation of compact pion stars during nucleation. Component 1 sets the existence and strength of the transition; component 2 converts it into a GW prediction; component 3 enumerates the remaining cosmological fingerprints. The existence question of component 1 is answered by the same standard of evidence whether the outcome is positive or negative, and components 2 and 3 follow from that outcome: a rigorous exclusion of a first-order segment answers component 1 negatively, in which case the gravitational-wave prediction and the remaining cosmological consequences take the form appropriate to the absence of a first-order transition.

References

1. O. Ferreira, E. S. Fraga, M. Hippert, J. Schaffner-Bielich, 'Chiral symmetry breaking and pion condensation in the early Universe', Phys. Rev. D 112, 094009 (2025). DOI: 10.1103/bcz6-xxn8. URL: https://arxiv.org/abs/2507.06518
2. V. Vovchenko, B. B. Brandt, F. Cuteri, G. Endrődi, F. Hajkarim, J. Schaffner-Bielich, 'Pion Condensation in the Early Universe at Nonvanishing Lepton Flavor Asymmetry and Its Gravitational Wave Signatures', Phys. Rev. Lett. 126, 012701 (2021). DOI: 10.1103/PhysRevLett.126.012701. URL: https://arxiv.org/abs/2009.02309
3. O. Ferreira, E. S. Fraga, J. Schaffner-Bielich, 'Signatures from pion condensation and lepton flavor asymmetries in the cosmological gravitational wave background', arXiv:2603.26946 (2026). DOI: 10.48550/arXiv.2603.26946. URL: https://arxiv.org/abs/2603.26946
4. B. B. Brandt, V. Chelnokov, G. Endrodi, G. Marko, D. Scheid, L. von Smekal, 'Renormalization group invariant mean-field model for QCD at finite isospin density', arXiv:2502.04025 (2025). DOI: 10.48550/arXiv.2502.04025. URL: https://arxiv.org/abs/2502.04025
5. R. Abbott, W. Detmold, F. Romero-López, Z. Davoudi, M. Illa, A. Parreño, R. J. Perry, P. E. Shanahan, M. L. Wagman, 'Lattice quantum chromodynamics at large isospin density', Phys. Rev. D 108, 114506 (2023). DOI: 10.1103/PhysRevD.108.114506. URL: https://arxiv.org/abs/2307.15014
6. B. B. Brandt, G. Endrodi, S. Schmalzbauer, 'QCD phase diagram for nonzero isospin-asymmetry', Phys. Rev. D 97, 054514 (2018). DOI: 10.1103/PhysRevD.97.054514. URL: https://arxiv.org/abs/1712.08190
7. G. Cao, 'Kaon superfluidity in the early Universe', Phys. Rev. D 110, 034004 (2024). DOI: 10.1103/PhysRevD.110.034004. URL: https://arxiv.org/abs/2405.07208
8. G. Cao, L. He, P. Zhang, 'Reentrant pion superfluidity and cosmic trajectories within a PNJL model', Phys. Rev. D 104, 054007 (2021). DOI: 10.1103/PhysRevD.104.054007. URL: https://arxiv.org/abs/2105.08932
9. J. B. Kogut, D. Toublan, 'QCD at small nonzero quark chemical potentials', Phys. Rev. D 64, 034007 (2001). DOI: 10.1103/PhysRevD.64.034007. URL: https://arxiv.org/abs/hep-ph/0103271
10. B. B. Brandt, G. Endrodi, E. S. Fraga, M. Hippert, J. Schaffner-Bielich, S. Schmalzbauer, 'New class of compact stars: Pion stars', Phys. Rev. D 98, 094510 (2018). DOI: 10.1103/PhysRevD.98.094510. URL: https://arxiv.org/abs/1802.06685
11. B. B. Brandt, G. Endrodi, G. Markó, 'Equation of state of isospin asymmetric QCD with small baryon chemical potentials', arXiv:2411.12918 (2024). DOI: 10.48550/arXiv.2411.12918. URL: https://arxiv.org/abs/2411.12918
12. R. F. Basta, B. B. Brandt, F. Cuteri, G. Endrődi, O. Philipsen, 'Pion condensation at non-zero isospin chemical potential with Wilson fermions', arXiv:2502.05051 (2025). DOI: 10.48550/arXiv.2502.05051. URL: https://arxiv.org/abs/2502.05051
13. F. Di Clemente, A. Drago, L. Formaggio, C. Ratti, V. Vovchenko, G. Yadav, 'Upper Bound on the Cosmic Baryon Chemical Potential from Lepton-Flavor Asymmetry', arXiv:2511.11995 (2025). DOI: 10.48550/arXiv.2511.11995. URL: https://arxiv.org/abs/2511.11995
14. H. Shah, T. Gyure, A. Leon, F. Di Clemente, M. Hippert, C. Ratti, V. Vovchenko, 'QCD critical surface from constant entropy contours', arXiv:2606.28282 (2026). DOI: 10.48550/arXiv.2606.28282. URL: https://arxiv.org/abs/2606.28282
15. G. Agazie et al. (NANOGrav Collaboration), 'The NANOGrav 15 yr Data Set: Evidence for a Gravitational-wave Background', Astrophys. J. Lett. 951, L8 (2023). DOI: 10.3847/2041-8213/acdac6. URL: https://arxiv.org/abs/2306.16213