Pion condensation versus 2SC, speed of sound, and charge neutrality effects in the quark-meson diquark model

Jens O. Andersen; Mathias P. N{\o}dtvedt

REVIEW 2 major objections 5 minor 3 cited by

Pion condensation versus 2SC, speed of sound, and charge neutrality effects in the quark-meson diquark model

T0 review · 2 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read The paper claims that the two-flavor quark-meson diquark model produces a four-phase $(\mu,\mu_I)$ phase diagram at $T=0$, and that in both the pion-condensed and neutral 2SC phases the squared speed of sound rises above the conformal…

desk verdict Worth a careful read for the pion-sector results and the explicit neutrality analysis, but the headline asymptotic claim about neutral 2SC matter is not self-consistent. read the letter →

arxiv 2502.10229 v2 pith:VEYIDNTR submitted 2025-02-14 hep-ph hep-latnucl-th

Jens O. Andersen , Mathias P. N{o}dtvedt This is my paper · ORCID

classification hep-phhep-latnucl-th

keywords quark-mesondiquarkmodelpioncondensation2SCcolorsuperconductivityspeedofsoundchargeneutralityisospinchemicalpotentialQCDphasediagramequationstate

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper uses a two-flavor quark-meson diquark (QMD) model, which combines quarks, mesons, and a diquark field with chemical potentials for baryon number, isospin, electric charge, and color charge, to map the zero-temperature phase diagram of dense QCD-like matter. It claims that the $(\mu,\mu_I)$ plane contains four phases: vacuum, charged-pion/Cooper-pair condensation, normal quark matter, and the two-flavor color-superconducting (2SC) phase. In both the pion-condensed and the neutral 2SC phases, the squared speed of sound rises above the conformal value $c_s^2=1/3$, peaks, and then approaches $1/3$ from above as density grows. For the pion sector this behavior matches lattice QCD results, while for the 2SC sector it implies an equation of state stiffer than conformal matter at moderate densities. The paper also classifies the Nambu-Goldstone bosons of the 2SC phase and discusses a mixed phase of normal and superconducting quark matter under global charge neutrality.

What carries the argument

The central object is the two-flavor quark-meson diquark (QMD) Lagrangian (Eq. 1), a renormalizable low-energy model with meson fields $\sigma$ and $\boldsymbol{\pi}$, a composite diquark antitriplet $\Delta_a$, and quark fields, all coupled through Yukawa and diquark interactions. The argument is carried by the thermodynamic potential $\Omega_{0+1}$ evaluated at tree level for bosons and one loop for quarks; its gap equations for $(\phi_0,\Delta_0)$ in the 2SC phase and for $(\phi_0,\rho_0)$ in the BEC/BCS phase select the global minimum and define the phase diagram. The key asymptotic identity is the large-density gap solution $\bar{\Delta}_0^2 = (m_q^2/g_{\Delta,0}^2)\exp[(4\pi)^2/(4g_{\Delta,0}^2)-1-F(m_\pi^2)-m_\pi^2F'(m_\pi^2)]$, which makes the gap approach a constant and, through the neutrality conditions, yields $c_s^2=(1/3)(1+0.71\,g_{\Delta,0}^2\bar{\Delta}_0^2/\mu^2)$ (Eq. 55). The classification of Nambu-Goldstone bosons in the 2SC phase uses a commutator counting rule for non-Lorentz-invariant systems, which gives one linear-dispersion (type-A) mode and two quadratic-dispersion (type-B) modes, or five linear modes when color neutrality is enforced.

What would settle it

A lattice QCD calculation of the equation of state at finite isospin chemical potential beyond $\mu_I/m_\pi\simeq 20$ would settle whether the squared speed of sound approaches $1/3$ from above, as the QMD model predicts, or from below, as the comparison calculation in Ref. [46] predicts.

Watch

Extended reading notes

Core claim

The central claim is that the QMD model yields a four-phase $(\mu,\mu_I)$ phase diagram at $T=0$, with a BEC/BCS pion-condensed phase at finite isospin and a 2SC phase at finite baryon density, and that in both phases the squared speed of sound has a maximum above $1/3$ and relaxes to the conformal value from above. Under electric and color charge neutrality, the asymptotic behavior is $c_s^2 = (1/3)(1+0.71\,g_{\Delta,0}^2\bar{\Delta}_0^2/\mu^2)$ (Eq. 55), so neutral 2SC matter is stiffer than conformal matter. At finite isospin the model's $c_s^2(\mu_I)$ curve lies inside the chiral perturbation theory band and agrees with lattice data, while in the 2SC phase imposing charge neutrality lowers and shifts the peak of $c_s^2$, softening the equation of state. The paper treats the diquark coupling and mass as free parameters, with the asymptotic gap depending only on $g_{\Delta,0}$ and the vacuum quark mass, and it uses an astrophysical upper bound on the color-flavor-locked gap to constrain one parameter set.

Load-bearing premise

The central predictions rest on treating the QMD Lagrangian (Eq. 1) as a faithful low-energy model of two-flavor QCD with the chosen free diquark mass $m_\Delta$ and coupling $g_\Delta$; if the true diquark parameters lie outside the chosen ranges, the 2SC onset, gap size, and speed-of-sound peak all move.

Editorial extensions

If this is right

In the pion-condensed BEC/BCS phase, the model's speed of sound as a function of $\mu_I/m_\pi$ tracks lattice QCD data over a wide range of isospin chemical potentials, so the QMD model can connect chiral perturbation theory at low density to perturbative QCD at high density.
In the 2SC phase at moderate baryon density, electrically and color neutral quark matter has $c_s^2>1/3$, meaning its equation of state is stiffer than conformal matter; this affects neutron star mass-radius relations and tidal deformability.
Imposing electric and color charge neutrality lowers the peak of $c_s^2$ and shifts it to higher chemical potential, so neutral 2SC matter is softer than unneutralized matter at the same density.
The asymptotic formula $c_s^2=(1/3)(1+0.71\,g_{\Delta,0}^2\bar{\Delta}_0^2/\mu^2)$ says the conformal limit is approached from above once neutrality is imposed, a distinctive prediction that contrasts with calculations approaching it from below.
A mixed phase of negatively charged normal quark matter and positively charged 2SC matter is possible under global color charge neutrality, giving a non-strange hybrid star scenario for compact star interiors.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

If the model is correct, neutron star observations that constrain the speed of sound in dense matter could directly constrain the diquark coupling $g_{\Delta,0}$, because the 2SC peak position and height depend mainly on it.
The close structural parallel between the pion-sector and diquark-sector thermodynamic potentials (Eqs. (29) and (43)) suggests the 'approach from above' behavior is robust within the model; computing higher-order corrections to Eq. (55) would test whether it survives beyond the leading $\bar{\Delta}_0^2/\mu^2$ term.
Extending the model to three flavors, which the paper notes is in progress, would replace the 2SC phase with the color-flavor-locked phase at high density; if that transition is first order, the speed of sound would jump, creating a distinctive signature in neutron star merger waveforms.
Current lattice data only reach $\mu_I/m_\pi\sim 20$; a lattice calculation at higher isospin density would directly test whether $c_s^2$ approaches $1/3$ from above, as the model predicts, rather than from below.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Plain take: the pion-sector part of this paper is solid and useful, but the headline asymptotic statement about c_s^2 in neutral 2SC matter does not survive contact with the model's own gapless criterion.

What's new: the mu-mu_I phase diagram in the QMD model, NG-boson counting with and without color neutrality, the neutrality dependence of the speed of sound, and the large-mu asymptotic formulas in Eqs. (53)-(55). The pion-sector c_s^2 curve agrees well with lattice QCD and chiral perturbation theory over a wide range of mu_I; that part is convincing and honestly compared. The authors also state their parameter sensitivity up front, which is to their credit.

The soft spot is in the derivation of Eq. (55). The neutrality solution (53) gives mu_e ~ 0.57 mu, so the mismatch delta-mu = mu_e/2 ~ 0.285 mu. The asymptotic gap (49) is mu-independent (for set 1, g_Delta Delta_0 is capped at 268 MeV). Hence for mu >~ 0.94 GeV, delta-mu exceeds the gap and the E^-_{Delta-} branch becomes gapless, i.e. the g2SC regime the paper itself defines in Sec. II C. But Eq. (50) is obtained from the gapped spectrum, dropping the theta-function terms in Eq. (31). So the derivation of Eq. (55) assumes delta-mu < g_Delta Delta_0 in a regime where imposing neutrality gives delta-mu > g_Delta Delta_0. The conformal-limit-from-above claim for neutral 2SC matter is therefore not established by this model. This is not parameter uncertainty; it is an internal consistency problem that the authors do not check.

Minor issues: the pion-sector analysis overlaps heavily with the simultaneously posted Ref. [41] (though the authors acknowledge this), and the numerical phase boundaries depend on free parameters without uncertainty estimates. Those are secondary.

Who it is for: people building effective-model equations of state for neutron stars and comparing pion-condensed matter to lattice data. The paper deserves a serious referee: the model work is careful, the lattice comparison is valuable, and the neutrality treatment is explicit. But the asymptotic claim needs to be fixed or qualified. I would send it to review with a request that the authors either redo the asymptotic neutral-2SC analysis including the gapless phase, or state clearly that Eq. (55) applies only while delta-mu < g_Delta Delta_0.

Referee Report

2 major / 5 minor

Summary. The paper studies the two-flavor quark-meson diquark (QMD) model at zero temperature as a low-energy effective model for QCD at finite quark and isospin chemical potentials. It maps out the T=0 phase diagram in the mu-mu_I plane, finding four phases: vacuum, charged pion/Cooper-pair condensation (BEC/BCS), normal quark matter, and 2SC color superconductivity. It classifies the Nambu-Goldstone modes of the 2SC phase using the Nielsen-Chadha/Watanabe-Murayama counting rules, and it studies local electric and color charge neutrality. The paper's central quantitative claims are the speed of sound: in the pion-condensed sector, c_s^2 has a maximum and approaches 1/3 from above at large mu_I, in agreement with lattice data, and in the 2SC sector, imposing charge neutrality gives the same qualitative behavior, with Eq. (55) stating the explicit large-mu approach from above. The paper closes with a discussion of mixed phases under global neutrality constraints.

Significance. If the central claims hold, the paper offers a single renormalizable low-energy model that connects the pion-condensed BEC/BCS regime, where lattice data exist, to the 2SC color-superconducting regime, where they do not, and it makes a sharp, falsifiable prediction: in neutral 2SC matter the speed of sound approaches the conformal value from above. The paper's strengths are genuine: the analytic large-chemical-potential derivations in Sec. II F are clean and explicit; the phase diagrams and speed-of-sound curves are compared honestly against lattice QCD and chiral perturbation theory; the Nambu-Goldstone classification in Sec. II D is careful and uses the modern counting rule; and the dependence on the model's free diquark parameters is acknowledged rather than hidden. The quantitative predictions, however, rest on the choice of m_Delta and g_Delta without uncertainty estimates, and, as detailed below, one of the flagship asymptotic statements has a self-consistency problem that must be resolved before the result can be accepted.

major comments (2)

[Sec. II F, Eqs. (50), (53), (55)] The derivation of Eq. (55) is not self-consistent for neutral 2SC matter at asymptotically large mu. The simplified pressure in Eq. (50) is obtained by dropping the theta-function terms in Omega_1^mu, Eq. (31), which requires the gapped-spectrum condition delta_mu < g_{Delta,0} Delta_0 for the quasiparticle mode E^-_{Delta-}. But the neutrality solution Eq. (53) gives mu_e ≈ 0.57 mu, hence delta_mu = mu_e/2 ≈ 0.285 mu, while the asymptotic gap in Eq. (49) is independent of mu (for set 1 it saturates at 268 MeV). For mu ≳ 0.94 GeV one therefore has delta_mu > g_{Delta,0} Delta_0. By the paper's own definition in Sec. II C, this is precisely the gapless 2SC (g2SC) regime, where the theta-function terms are nonzero and the expansion underlying Eq. (50) breaks down. Consequently, the claim that neutral 2SC matter approaches c_s^2 = 1/3 from above, Eq. (55), is not established by the model. The authors should either impose the g2SC condition alongside the neutrality equations and recompute pressure and c_s^2 including the gapless contributions, or explicitly restrict the asymptotic claim to the parameter regime where the gapped expansion is valid and state that the full mu dependence is not covered.
[Sec. III, parameter sets and Fig. 7] The numerical phase boundaries and the height and position of the c_s^2 peak depend strongly on the free diquark mass m_Delta and coupling g_Delta. The authors state this sensitivity explicitly in Sec. III, but the two parameter sets produce qualitatively different phase diagrams (for set 1 there is no normal quark matter window in Figs. 1 and 3, while for set 2 there is a substantial NQM phase), and no uncertainty or calibration range is attached to the resulting speed-of-sound curves in Fig. 7. Because the paper presents the 2SC speed-of-sound peak and its approach from above as generic model features, the absence of any estimate of how the chosen parameter values affect these conclusions weakens the quantitative comparison. The authors should quantify the variation of the peak and of the onset of the 2SC phase over a plausible parameter range, or clearly frame Figs. 3-7 as illustrative only.

minor comments (5)

[Sec. III, text before Fig. 5] The sentence 'good qualitative agreement with the lattice data. although the peak in c_s^2 occurs much earlier than in the QM model' contains a grammatical error: 'although' should be 'though' or 'whereas'.
[Fig. 5 caption] The caption uses 'PT' instead of the standard 'chiPT' or 'ChPT' used elsewhere; please make the notation uniform.
[Sec. III, discussion of Fig. 7] The sentence 'The blue curves will cross each other for higher values of mu' is unclear because the crossing is not shown in the figure; if the crossing is predicted by the asymptotics, say so explicitly and give the estimated crossing scale, or remove the statement.
[References] Reference [5] has a formatting typo: 'Ann. Phys. 158, (142) (1984)' should be 'Ann. Phys. 158, 142 (1984)'.
[Sec. II F, footnote 3] The parenthetical note about the factor 2^{-1/3} relating the 2SC and CFL gaps is useful but should be integrated into the main text or given a proper derivation, since it is used implicitly when comparing with the astrophysical bound from Ref. [43].

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the QMD predictions are benchmarked against external lattice QCD and chiral perturbation theory, and the speed-of-sound results follow from the model's equation of state with parameters fixed before the comparison.

full rationale

The paper's central claims—the four-phase diagram and the speed-of-sound curves—are not obtained by fitting the target quantities. The pion-sector c_s is computed from the thermodynamic potential (Eq. 43) and compared with lattice data [7,9] and chiral perturbation theory [10]; no lattice c_s value enters the parameter choice. The 2SC-sector c_s (Eq. 55) follows from the asymptotic pressure and energy density (Eqs. 50-52) and the neutrality solutions (Eqs. 53-54); these are derived within the model, and the statement that the conformal limit is approached from above is a consequence of the positive O(g^2 \bar\Delta^2/\mu^2) correction, not an input. The diquark parameters m_\Delta and g_\Delta are free model inputs (Sec. III), and the paper notes that g_\Delta could in principle be tuned to the perturbative-QCD gap, but the two published parameter sets are not fitted to the lattice speed-of-sound data used for validation. The self-citation to Ref. [18] supplies the RG running and one-loop quark potential; that is a parameter-free calculation with stated assumptions that does not assume the present paper's speed-of-sound or phase-diagram results, so it is independent support rather than circular load-bearing. The skeptical observation that Eq. (55) may violate the gapped-phase condition δμ < g_\Delta\Delta_0 is a possible self-consistency or correctness issue for the asymptotic neutral-2SC analysis, not a circularity: it does not make the prediction equal to its input by construction.

Assumptions & free parameters 4 free parameters · 4 assumptions · 0 invented entities

The central results depend on the model's free couplings; only m_pi and f_pi are tied directly to lattice or chiral perturbation theory. The diquark field is a composite model degree of freedom from Ref. [18], not a new physical particle. The 2SC and sound-speed predictions are conditional on the chosen m_Delta and g_Delta.

free parameters (4)

m_Delta (diquark mass parameter) = set 1: 500 MeV; set 2: 900 MeV
Free parameter controlling the 2SC transition; the paper says results are sensitive mainly to m_Delta and g_Delta in Sec. III.
g_Delta (diquark-quark coupling) = set 1: 2g; set 2: 1.5g
Controls gap magnitude and exact speed-of-sound behavior; explicitly tunable to match the perturbative QCD gap (Sec. II F, Eq. 49).
lambda_3 and lambda_Delta couplings = set 1: lambda_0, lambda_0/4; set 2: 0, lambda_0/5
Chosen by hand; they affect the thermodynamic potential and phase boundaries.
m_sigma and m_q = 600 MeV and 300 MeV
Inputs chosen to match typical quark-meson model phenomenology; not derived from QCD itself.

assumptions (4)

domain assumption The QMD Lagrangian Eq. (1) is the correct low-energy effective model for two-flavor QCD at the densities considered.
The entire phase diagram and equation of state follow from this model; QCD itself is not solved.
domain assumption SU(3)_c is treated as a global symmetry, replacing the local gauge symmetry of QCD.
Quoted in Sec. II A; this changes NG-boson counting and removes gluons, making the 2SC properties partly model artifacts.
domain assumption One-loop quark fluctuations with tree-level bosons are sufficient.
Used for the thermodynamic potentials in Eqs. (29) and (43); no systematic error estimate for neglected boson loops is given.
domain assumption The renormalization-group running from Ref. [18] is correct.
The paper relies on the companion paper for running couplings and the RG-improved effective potential.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Pion condensation versus 2SC, speed of sound, and charge neutrality effects in the quark-meson diquark model." pith.science (2026). https://pith.science/paper/VEYIDNTR

@misc{pith2026250210229,
  author       = {Pith},
  title        = {Pith review of: Pion condensation versus 2SC, speed of sound, and charge neutrality effects in the quark-meson diquark model},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/VEYIDNTR}},
  note         = {Machine review of arXiv:2502.10229}
}

abstract

We employ the two-flavor quark-meson diquark model as a low-energy model for QCD at non-zero quark and isospin chemical potentials $\mu$ and $\mu_I$, and at zero temperature. We map out the phase diagram in the $\mu$-$\mu_I$ plane, which has four phases: a vacuum phase, a phase with condensed charged pions/Cooper pairs of $u$ and $\bar{d}$ quarks, a normal quark matter phase, and a color superconducting phase (2SC phase). %In the 2SC phase, we study the effects of imposing color and %electric charge neutrality. The global symmetry breaking $SU(3)_c\rightarrow SU(2)_c$ in the 2SC phase gives rise to a number of Nambu-Goldstone bosons. We classify them and briefly discuss their properties. We calculate the speed of sound $c_s$ in the two special cases, finite $\mu_I$ and $\mu=0$, and finite $\mu$ and $\mu_I=0$. In both cases, the speed of sound exhibits a maximum and approaches the conformal limit from above as the density increases. For non-zero isospin $\mu_I$, this behavior is in agreement with the speed of sound obtained from lattice simulations. In the 2SC phase, the behavior is qualitatively the same if we impose local charge neutrality. Finally, we discuss the possibility of having a mixed phase of negatively charged normal quark matter and positively charged 2SC matter with global color charge neutrality imposed on the latter.

Figures

Figures reproduced from arXiv: 2502.10229 by the authors.

**Figure 2.** FIG. 2. Phase diagram in the [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗

**Figure 4.** FIG. 4. Quark condensate (black line), superconducting gap [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗

Figures from the paper (3 more)

**Figure 6.** Figure 6: FIG. 6. Speed of sound squared in the pion-condensed phase [PITH_FULL_IMAGE:figures/full_fig_p011_6.png]

**Figure 5.** Figure 5: FIG. 5. Speed of sound squared in the pion-condensed phase [PITH_FULL_IMAGE:figures/full_fig_p011_5.png]

**Figure 7.** Figure 7: FIG. 7. Speed of sound squared in the 2SC phase for pa [PITH_FULL_IMAGE:figures/full_fig_p011_7.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 3 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

Renormalization of the three-flavor quark-meson diquark model
hep-ph 2025-06 conditional novelty 6.0 of 10

The authors extend the quark-meson diquark model to three flavors, renormalize its thermodynamic potential in the CFL phase, and show the speed of sound approaches the conformal limit from above.
Renormalizing the Quark-Meson-Diquark Model
hep-ph 2025-05 conditional novelty 6.0 of 10

Renormalized and two RG-consistent versions of the Quark-Meson-Diquark model reproduce the BCS relation and Stefan-Boltzmann limit at high density, while the sigma-delta scheme violates the BCS relation.
Dense and Cold Magnetized Quark Matter: A Review of Magnetic-Field-Independent Regularization and the Medium Separation Scheme
hep-ph 2026-06 unverdicted novelty 3.5 of 10

Review of MFIR and MSS schemes showing the superconducting gap stays finite at high chemical potential in magnetized cold quark matter with no zero-temperature transition to normal phase.

Reference graph

Works this paper leans on

61 extracted references · 52 canonical work pages · cited by 3 Pith papers

[41]

Ayala, A

A. Ayala, A. Bandyopadhyay, R. L. S. Farias, L. A. Hern´ andez, and J. L. Hern´ andez, Phys. Rev. D 107, 074027 (2023)

work page 2023
[1]

(31) Note that the integrals in Eqs

Moreover, the function F (p2) = 2 − 2q arctan( 1 q ) with q = q 4m2q p2 − 1, and Ωfin 1 = −4Λ−2ϵ Z p " E± ∆± − 2 q p2 + g2ϕ2 0 + g2 ∆∆2 0 − ¯µ2g2 ∆∆2 0 (p2 + g2ϕ2 0 + g2 ∆∆2 0) 3 2 # , (30) Ωµ 1 = 2Λ −2ϵ Z p (E ∓ µub)θ(±µub − E) + (E ∓ µdb)θ(±µdb − E) + 2(E± ∆ ∓ δµ)θ(±δµ − E± ∆) . (31) Note that the integrals in Eqs. (30)–(31) are finite in d = 3. D. Namb...

work page
[2]

4m2 q m2σ log ϕ2 0+ρ2 0 f 2π − 3 2 − 1 − 4m2 q m2σ F (m2 σ) + F (m2 π) + m2 πF ′(m2 π) #) (ϕ2 0 + ρ2 0)2 f 4π − 1 8 m2 πf 2 π

This yields the dispersion relations E± ∆1 = p p2 + 4µ2 ± 2µ , (33) E± ∆2 = p p2 + 4µ2 ± 2µ , (34) E± ∆3 = r p2 + 12µ2 − m2 ∆ ± q 16p2µ2 + (12µ2 − m2 ∆)2 . (35) The massless NG modes in Eqs. (33)–(34) are quadratic for small momenta p, while the massless NG mode in Eq. (35) is linear for small momenta. Linear modes are referred to as type A and quadratic ...

work page
[3]

(45) Here the fermionic quasiparticle spectrum is E± ρ0 = q (E ± 1 2 µI )2 + g2 0ρ2 0 , (46) with E = p p2 + g2 0ϕ2

− µ2 I ρ2 0 4[p2 + g2 0(ϕ2 0 + ρ2 0)] 3 2 # , (44) Ωµ 1 = 2 NcΛ−2ϵ Z p (E± ρ0 − µ)θ(µ − E± ρ0 ) + (E± ρ0 + µ)θ(−µ − E± ρ0 ) . (45) Here the fermionic quasiparticle spectrum is E± ρ0 = q (E ± 1 2 µI )2 + g2 0ρ2 0 , (46) with E = p p2 + g2 0ϕ2

work page
[4]

log m2 q g2 ∆,0∆2 0 − F (m2 π) − m2 πF ′(m2 π) # ¯µ2∆2 0 + 4g4 ∆,0 (4π)2

Note the similarity between the gapped spectra Eq. (22) and Eq. (46) for δµ = 0. The integral Eqs. (44) is evaluated numerically directly ind = 3 dimensions. Eqs. (29) and (43) is the starting point for mapping out the phase diagram. F. Behavior at large chemical potentials For large chemical potentials µ or µI , we can obtain simple expressions for therm...

work page 2000
[5]

Rajagopal and F

K. Rajagopal and F. Wilczek, At the frontier of parti- cle physics, Vol. 3 (World Scientific, Singapore, p 2061) (2001)

2001
[6]

M. G. Alford, A. Schmitt, K. Rajagopal, and T. Sch¨ afer, Rev. Mod. Phys. 80, 1455 (2008)

work page 2008
[7]

Fukushima and T

K. Fukushima and T. Hatsuda, Rept. Prog. Phys. 74, 014001 (2011)

2011

Show all 61 references

[8]

D. T. Son and M. A. Stephanov, Phys. Rev. Lett. 86, 592 (2001)

2001
[9]

Gasser and H

J. Gasser and H. Leutwyler, Ann. Phys. 158, (142) (1984)

1984
[10]

Gasser and H

J. Gasser and H. Leutwyler, Nucl. Phys. B 250, 465 (1985)

1985
[11]

B. B. Brandt, F. Cuteri, and G. Endrodi, JHEP 07, 055 (2023)

2023
[12]

Abbott et al., NPLQCD Collaboration, Phys

R. Abbott et al., NPLQCD Collaboration, Phys. Rev. D 108, 11, 114506 (2023)

2023
[13]

Abbott et al., Phys

R. Abbott et al., Phys. Rev. Lett. 134, 011903 (2025)

2025
[14]

J. O. Andersen, M. K. Johnsrud, Q. Yu, and H. Zhou, Phys. Rev. D 111, 034017 (2025)

2025
[15]

P. F. Bedaque and T. Sch¨ afer, Nucl. Phys. A 697, 802 (2002)

2002
[16]

Alford and K

M. Alford and K. Rajagopal, JHEP 06, 03 (2002)

2002
[17]

S. B. R¨ uster, V. Werth, M. Buballa, I. A. Shovkovy, and D. H. Rischke, Phys. Rev. D 72, 034004 (2005)

2005
[18]

Abuki and T

H. Abuki and T. Kunihiro, Nucl. Phys. A 768, 118 (2006)

2006
[19]

Gorda, A

T. Gorda, A. Kurkela, P. Romatschke, S. S¨ appi, and A. Vuorinen, Phys. Rev. Lett. 121, 202701 (2018)

2018
[20]

Gorda, R

T. Gorda, R. Paatelainen, S. S¨ appi, and K. Sepp¨ anen, Phys. Rev. Lett. 131, 181902 (2023)

2023
[21]

Fernandez and J.-L

L. Fernandez and J.-L. Kneur, Phys. Rev. Lett. 129, 212001 (2022); Phys. Rev. D 111, 034020 (2025)

2022
[22]

J. O. Andersen and M. P. Nødtvedt, Phys. Rev. D 111, 034031 (2025)

2025
[23]

I. A. Shovkovy, Found. Phys. 35, 1309 (2005)

2005
[24]

Gerhold and A

A. Gerhold and A. Rebhan, Phys. Rev. D 68, 011502 (2003)

2003
[25]

D. D. Dietrich, and D. H. Rischke, Prog. Part. Nucl. Phys. 53, 305 (2004)

2004
[26]

I. A. Shovkovy and M. Huang, Phys. Lett. B 564, 205 (2003)

2003
[27]

Casalbuoni, Z

R. Casalbuoni, Z. Duan, and F. Sannino, Phys. Rev. D 62, 094004 (2000)

2000
[28]

Blaschke, D

D. Blaschke, D. Ebert, K. G. Klimenko, M. K. Volkov, and V. L. Yudichev, Phys. Rev. D 70, 014006 (2004). 13

2004
[29]

H. B. Nielsen and S. Chadha, Nucl. Phys. B 105, 445 (1976)

1976
[30]

Watanabe and H

H. Watanabe and H. Murayama, Phys. Rev. Lett. 108, 251602 (2012)

2012
[31]

Hidaka, Phys

Y. Hidaka, Phys. Rev. Lett. 110, 091601 (2013)

2013
[32]

Nicolis, R

A. Nicolis, R. Penco, F. Piazza, and R. A. Rosen, JHEP 11, 055 (2013)

2013
[33]

Watanabe, T

H. Watanabe, T. Brauner, and H. Murayama, Phys. Rev. Lett. 111, 021601 (2013)

2013
[34]

Brauner, Effective Field Theory for Spontaneously Broken Symmetry, Lect

T. Brauner, Effective Field Theory for Spontaneously Broken Symmetry, Lect. Notes Phys. 1023 (2024)

2024
[35]

Sch¨ afer, D

T. Sch¨ afer, D. T. Son, and M. A. Stephanov, D. Toublan, J. J. M. Verbaarschot, Phys. Lett. B 522, 67 (2001)

2001
[36]

Nicolis and F

A. Nicolis and F. Piazza, Phys. Rev. Lett. 110, 011602 (2013), ibid 110, 039901 (2013) (addendum)

2013
[37]

P. F. Bedaque, Nucl. Phys. A 697, 569 (2002)

2002
[38]

L.-Y. He, M. Jin, and P.-F. Zhuang, Phys. Rev. D 71, 116001 (2005)

2005
[39]

Kamikado, N

K. Kamikado, N. Strodthoff, L. von Smekal, and J. Wambach, Phys. Lett. B 718 (1044) (2013)

2013
[40]

J. O. Andersen and P. Kneschke, Phys. Rev. D 97, 076005 (2018)

2018
[42]

Chiba and T

R. Chiba and T. Kojo, Phys. Rev. D 109, 076006 (2024)

2024
[43]

Ayala, B

A. Ayala, B. S. Lopes, R. L. S. Farias, and L. C. Parra, e-Print: 2409.19406 [hep-ph]

arXiv
[44]

J. P. Carlomagno, D. G. Dumm, and N. N. Scoccola, Phys. Rev. D 109, 9, 094041 (2024)

2024
[45]

B. B. Brandt, V. Chelnokov, G. Endrodi, G. Marko, D. Scheid, and L. von Smekal, e-Print: 2502.04025 [hep-ph]

arXiv
[46]

A. W. Steiner, S. Reddy, and M. Prakash, Phys. Rev. D 66, 094007 (2002)

2002
[47]

Kurkela, K

A. Kurkela, K. Rajagopal, and R. Steinhorst, Phys. Rev. Lett. 132, 26, 262701 (2024)

2024
[48]

Huang, P.-F

M. Huang, P.-F. Zhuang, and W.-Q. Chao, Phys. Rev. D 67, 065015 (2003)

2003
[49]

Gholami, M

H. Gholami, M. Hofmann, and M. Buballa, Phys. Rev. D 111, 014006 (2025)

2025
[50]

Fukushima and S

K. Fukushima and S. Minato, e-Print: 2411.03781 [hep- ph]

arXiv
[51]

J. M. Cornwall, R. Jackiw, and E. Tomboulis, Phys. Rev. D 10, 2428 (1974)

1974
[52]

G. Baym, T. Hatsuda, T. Kojo, P. D. Powell, Y. Song, and T. Takatsuka, Rept. Prog. Phys. 81, 056902 (2018)

2018
[53]

McLerran and S

L. McLerran and S. Reddy, Phys. Rev. Lett. 122, 122701 (2019)

2019
[54]

Altiparmak, C

S. Altiparmak, C. Ecker, and L. Rezzolla, Astrophys. J. Lett. 939, L34 (2022)

2022
[55]

Marczenko, L

M. Marczenko, L. McLerran, K. Redlich, and C. Sasaki, Phys. Rev. C 107, 025802 (2023)

2023
[56]

Brandes, W

L. Brandes, W. Weise, and N. Kaiser, Phys. Rev. D 108, 094014 (2023)

2023
[57]

Albino, T

M. Albino, T. Malik, M. Ferreira, and C. Providˆ encia, Phys. Rev. D 110, 083037, (2024)

2024
[58]

Ayriyan, D

A. Ayriyan, D. Blaschke, J. P. Carlomagno, G. A. Con- trera, and A. G. Grunfeld, e-Print:2501.00115 [nucl-th]

arXiv
[59]

N. K. Glendenning, Phys. Rev. D 46, 1274 (1992)

1992
[60]

I. A. Shovkovy, M. Hanauske, and M. Huang, Phys. Rev. D 67, 103004 (2003)

2003
[61]

J. O. Andersen and N. P. Nødtvedt, Work in progress

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.

[1] [41]

Ayala, A

A. Ayala, A. Bandyopadhyay, R. L. S. Farias, L. A. Hern´ andez, and J. L. Hern´ andez, Phys. Rev. D 107, 074027 (2023)

work page 2023

[2] [1]

(31) Note that the integrals in Eqs

Moreover, the function F (p2) = 2 − 2q arctan( 1 q ) with q = q 4m2q p2 − 1, and Ωfin 1 = −4Λ−2ϵ Z p " E± ∆± − 2 q p2 + g2ϕ2 0 + g2 ∆∆2 0 − ¯µ2g2 ∆∆2 0 (p2 + g2ϕ2 0 + g2 ∆∆2 0) 3 2 # , (30) Ωµ 1 = 2Λ −2ϵ Z p (E ∓ µub)θ(±µub − E) + (E ∓ µdb)θ(±µdb − E) + 2(E± ∆ ∓ δµ)θ(±δµ − E± ∆) . (31) Note that the integrals in Eqs. (30)–(31) are finite in d = 3. D. Namb...

work page

[3] [2]

4m2 q m2σ log ϕ2 0+ρ2 0 f 2π − 3 2 − 1 − 4m2 q m2σ F (m2 σ) + F (m2 π) + m2 πF ′(m2 π) #) (ϕ2 0 + ρ2 0)2 f 4π − 1 8 m2 πf 2 π

This yields the dispersion relations E± ∆1 = p p2 + 4µ2 ± 2µ , (33) E± ∆2 = p p2 + 4µ2 ± 2µ , (34) E± ∆3 = r p2 + 12µ2 − m2 ∆ ± q 16p2µ2 + (12µ2 − m2 ∆)2 . (35) The massless NG modes in Eqs. (33)–(34) are quadratic for small momenta p, while the massless NG mode in Eq. (35) is linear for small momenta. Linear modes are referred to as type A and quadratic ...

work page

[4] [3]

(45) Here the fermionic quasiparticle spectrum is E± ρ0 = q (E ± 1 2 µI )2 + g2 0ρ2 0 , (46) with E = p p2 + g2 0ϕ2

− µ2 I ρ2 0 4[p2 + g2 0(ϕ2 0 + ρ2 0)] 3 2 # , (44) Ωµ 1 = 2 NcΛ−2ϵ Z p (E± ρ0 − µ)θ(µ − E± ρ0 ) + (E± ρ0 + µ)θ(−µ − E± ρ0 ) . (45) Here the fermionic quasiparticle spectrum is E± ρ0 = q (E ± 1 2 µI )2 + g2 0ρ2 0 , (46) with E = p p2 + g2 0ϕ2

work page

[5] [4]

log m2 q g2 ∆,0∆2 0 − F (m2 π) − m2 πF ′(m2 π) # ¯µ2∆2 0 + 4g4 ∆,0 (4π)2

Note the similarity between the gapped spectra Eq. (22) and Eq. (46) for δµ = 0. The integral Eqs. (44) is evaluated numerically directly ind = 3 dimensions. Eqs. (29) and (43) is the starting point for mapping out the phase diagram. F. Behavior at large chemical potentials For large chemical potentials µ or µI , we can obtain simple expressions for therm...

work page 2000

[6] [5]

Rajagopal and F

K. Rajagopal and F. Wilczek, At the frontier of parti- cle physics, Vol. 3 (World Scientific, Singapore, p 2061) (2001)

2001

[7] [6]

M. G. Alford, A. Schmitt, K. Rajagopal, and T. Sch¨ afer, Rev. Mod. Phys. 80, 1455 (2008)

work page 2008

[8] [7]

Fukushima and T

K. Fukushima and T. Hatsuda, Rept. Prog. Phys. 74, 014001 (2011)

2011