Chiral-Flavor Vacuum in QCD and Field Mixing

• Chiral-flavor vacuum is a dynamically structured ground state in QCD models that integrates chiral condensates with flavor and anomaly information.
• It arises in varied contexts—from fermionic zero-point contributions to Bogoliubov condensates—impacting symmetry breaking and phase transitions.
• Understanding this vacuum guides improvements in modeling criticality, field mixing, and anomaly resolution in quantum field theory.

Chiral-flavor vacuum denotes a ground-state structure in which chiral order parameters and flavor degrees of freedom are dynamically inseparable. In QCD-motivated effective theories it can mean the combined classical mesonic background and fermionic Dirac sea that determines the chiral condensate and the order of chiral restoration; in finite-$\theta$, finite-isospin, or background-field problems it can mean a vacuum reoriented in flavor space; and in quantum-field-theoretic treatments of mixing it can mean a Bogoliubov condensate vacuum associated with operators of definite chirality and flavor rather than definite mass or energy (Skokov et al., 2010, Sasaki et al., 2011, Adhikari et al., 2019, Blasone et al., 13 Jul 2025). The unifying point is that the vacuum is not merely a flavor-blind scalar background: it carries explicit information about anomaly structure, condensate alignment, zero-point fluctuations, or inequivalent Fock representations.

1. Definitions and conceptual range

In two-flavor QM and PQM models, “vacuum” includes the fermionic zero-point energy

$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$

so that the chiral condensate is determined jointly by the mesonic potential and the Dirac sea. The corresponding “no-sea approximation” omits this term even though it is the one-loop fermion determinant at $T=0$, not a mesonic vacuum expectation value (Skokov et al., 2010).

In three-flavor PNJL and EPNJL at finite $\theta$, the vacuum is specified by flavor-diagonal scalar and pseudoscalar condensates,

$$\sigma_f=\bar q_f q_f,\qquad \eta_f=\bar q_f i\gamma_5 q_f,$$

together with $\Phi,\Phi^*$. The thermodynamic vacuum is obtained from

$$\frac{\partial \Omega}{\partial X}=0,\qquad X=\Phi,\Phi^*,\sigma_f,\eta_f,$$

with the selected solution required to be the global minimum (Sasaki et al., 2011).

In finite-isospin $\chi$PT, the vacuum is an $SU(2)$ matrix field

$$\Sigma_\alpha=\cos\alpha+i\hat\phi_i\tau_i\sin\alpha,$$

whose orientation in flavor space is itself the order parameter. The normal vacuum corresponds to $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$0, while the pion-condensed vacuum is a rotated state in the charged-pion plane (Adhikari et al., 2019).

In QFT treatments of flavor-chiral oscillations, the vacuum is the state annihilated by flavor-chiral ladder operators. It is explicitly a condensate of particle-antiparticle pairs with definite masses and helicities, rather than the standard Dirac mass vacuum (Blasone et al., 13 Jul 2025).

Taken together, these usages indicate that the expression is context-dependent. A plausible synthesis is that “chiral-flavor vacuum” names the ground state once one has fixed which chiral and flavor charges are diagonal, which condensates are nonzero, and which fluctuation sectors are retained.

2. Fermion sea, effective potentials, and chiral criticality

The two-flavor QM/PQM literature makes the fermion vacuum term decisive. In the PQM model,

$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$1

with $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$2. In the chiral limit, dropping $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$3 yields a Landau form

$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$4

so the $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$5 term drives the quartic coefficient effectively negative at small $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$6, generating two competing minima and a barrier. The explicit conclusion is that in this no-sea approximation the chiral transition in the chiral limit is always first order. Including the fermion vacuum term cancels the thermal logarithm and restores an ordinary quartic Landau form,

$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$7

with second-order behavior for sufficiently large cutoff; numerically, the quoted threshold is $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$8 MeV. In dimensional regularization the renormalized vacuum piece becomes

$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$9

and the full mean-field potential is scale independent after parameter running (Skokov et al., 2010).

This vacuum restructuring propagates into flavor-sensitive thermodynamics. The generalized susceptibilities

$T=0$0

are strongly cutoff dependent if the fermion vacuum is included only in a cutoff scheme, and the sharp peak in $T=0$1 near $T=0$2 at small $T=0$3 is interpreted as an unphysical relic of the underlying first-order tendency. In the renormalized treatment, $T=0$4 develops a kink and $T=0$5 a rapid drop near the crossover; in the chiral limit these become nonanalytic with

$T=0$6

while the Polyakov loop modifies only the thermal sector and does not eliminate the need to retain the fermion sea (Skokov et al., 2010).

One-loop OPT studies reach a related conclusion from a different reorganization of perturbation theory. The free energy

$T=0$7

changes qualitatively when different zero-point pieces are kept. In one setup, keeping only $T=0$8 gives a first-order transition with $T=0$9 MeV; adding $\theta$0 changes the transition to second order with little change in $\theta$1; adding $\theta$2 instead gives a first-order transition with $\theta$3 MeV; and including all one-loop terms gives first order with $\theta$4 MeV because bosonic and fermionic vacuum terms partially cancel (Khan et al., 2010). In the related two-flavor QM analysis with explicit symmetry breaking, the full one-loop calculation yields a crossover in the whole $\theta$5-$\theta$6 plane at the physical point, while in the chiral limit the full one-loop result remains first order for all $\theta$7 considered (Andersen et al., 2011). A recurrent controversy in this class of models is therefore not whether the vacuum matters, but which vacuum terms are essential and how they must be renormalized.

3. Anomaly, $\theta$8 dependence, and topological vacua

At finite $\theta$9, the three-flavor PNJL/EPNJL vacuum is controlled by the $$\sigma_f=\bar q_f q_f,\qquad \eta_f=\bar q_f i\gamma_5 q_f,$$0-breaking KMT determinant. After the axial rotation

$$\sigma_f=\bar q_f q_f,\qquad \eta_f=\bar q_f i\gamma_5 q_f,$$1

the determinant becomes $$\sigma_f=\bar q_f q_f,\qquad \eta_f=\bar q_f i\gamma_5 q_f,$$2-independent and the light-quark mass sector becomes

$$\sigma_f=\bar q_f q_f,\qquad \eta_f=\bar q_f i\gamma_5 q_f,$$3

This makes the main $$\sigma_f=\bar q_f q_f,\qquad \eta_f=\bar q_f i\gamma_5 q_f,$$4-vacuum effect transparent: increasing $$\sigma_f=\bar q_f q_f,\qquad \eta_f=\bar q_f i\gamma_5 q_f,$$5 primarily reduces the effective even light mass $$\sigma_f=\bar q_f q_f,\qquad \eta_f=\bar q_f i\gamma_5 q_f,$$6, while the odd part remains parametrically small. At $$\sigma_f=\bar q_f q_f,\qquad \eta_f=\bar q_f i\gamma_5 q_f,$$7, $$\sigma_f=\bar q_f q_f,\qquad \eta_f=\bar q_f i\gamma_5 q_f,$$8 is even in $$\sigma_f=\bar q_f q_f,\qquad \eta_f=\bar q_f i\gamma_5 q_f,$$9 and develops a cusp at $\Phi,\Phi^*$0, with $\Phi,\Phi^*$1 and $\Phi,\Phi^*$2 jumping there; this is identified as the Dashen phenomenon, i.e. spontaneous parity breaking. At $\Phi,\Phi^*$3, parity restoration with increasing $\Phi,\Phi^*$4 is second order in PNJL, with $\Phi,\Phi^*$5 MeV, but first order in EPNJL, with $\Phi,\Phi^*$6 MeV. The paper further states that for sufficiently large $\Phi,\Phi^*$7 the chiral transition becomes first order even at $\Phi,\Phi^*$8 in EPNJL (Sasaki et al., 2011).

The Stern phase pushes the anomaly discussion in a different direction. It assumes

$\Phi,\Phi^*$9

while still allowing

$$\frac{\partial \Omega}{\partial X}=0,\qquad X=\Phi,\Phi^*,\sigma_f,\eta_f,$$0

through higher-order condensates. In the $$\frac{\partial \Omega}{\partial X}=0,\qquad X=\Phi,\Phi^*,\sigma_f,\eta_f,$$1, $$\frac{\partial \Omega}{\partial X}=0,\qquad X=\Phi,\Phi^*,\sigma_f,\eta_f,$$2 effective theory the $$\frac{\partial \Omega}{\partial X}=0,\qquad X=\Phi,\Phi^*,\sigma_f,\eta_f,$$3-vacuum energy behaves, in the macroscopic limit, as

$$\frac{\partial \Omega}{\partial X}=0,\qquad X=\Phi,\Phi^*,\sigma_f,\eta_f,$$4

If $$\frac{\partial \Omega}{\partial X}=0,\qquad X=\Phi,\Phi^*,\sigma_f,\eta_f,$$5, this produces two first-order transitions in $$\frac{\partial \Omega}{\partial X}=0,\qquad X=\Phi,\Phi^*,\sigma_f,\eta_f,$$6, generally away from $$\frac{\partial \Omega}{\partial X}=0,\qquad X=\Phi,\Phi^*,\sigma_f,\eta_f,$$7. For

$$\frac{\partial \Omega}{\partial X}=0,\qquad X=\Phi,\Phi^*,\sigma_f,\eta_f,$$8

the vacuum minimizes at $$\frac{\partial \Omega}{\partial X}=0,\qquad X=\Phi,\Phi^*,\sigma_f,\eta_f,$$9 and enters an exotic Aoki-like phase that breaks vector flavor symmetry,

$\chi$0

with two exactly massless modes (Kanazawa, 2015). This directly contradicts the common assumption that spontaneous chiral symmetry breaking necessarily implies a nonzero bilinear condensate.

A broader topological formulation starts from nonzero topological susceptibility in the fermion-free theory,

$\chi$1

or its gravitational analogue. In that framework, adding $\chi$2 massless fermions forces dynamical breaking of chiral flavor symmetry down to a maximal anomaly-free subgroup, makes the axial $\chi$3 pseudo-Goldstone massive through three-form screening, and removes massless fermions from the low-energy spectrum without assuming confinement (Dvali, 2017). This suggests that topology can be an organizing principle of the chiral-flavor vacuum, rather than a perturbation around an otherwise fixed vacuum manifold.

4. Vacuum reorientation by $\chi$4 and by magnetic fields

At nonzero isospin chemical potential, the two-flavor chiral vacuum rotates continuously in flavor space. The LO static potential is

$\chi$5

after minimizing over the direction in isospin space, which yields $\chi$6. The stationary condition gives either $\chi$7 or

$\chi$8

so the rotated solution exists only for $\chi$9. After expressing the theory in terms of physical quantities, the critical point is

$SU(2)$0

The Landau-Ginzburg expansion

$SU(2)$1

together with $SU(2)$2, shows that the vacuum-to-pion-condensed transition is second order. Above threshold, the vacuum lies in the charged-pion plane and spontaneously breaks $SU(2)$3, producing a Goldstone mode (Adhikari et al., 2019).

In a magnetized three-flavor NJL model with KMT interaction, the vacuum is written explicitly as a Bogoliubov state of quark-antiquark pairs,

$SU(2)$4

with flavor-dependent condensate functions and constituent masses

$SU(2)$5

This makes the vacuum simultaneously flavor resolved and flavor mixed. With the quoted parameter set, the vacuum constituent masses are $SU(2)$6 GeV and $SU(2)$7 GeV. In a magnetic field, the enhancement factors at $SU(2)$8 are approximately $SU(2)$9, $$\Sigma_\alpha=\cos\alpha+i\hat\phi_i\tau_i\sin\alpha,$$0, and $$\Sigma_\alpha=\cos\alpha+i\hat\phi_i\tau_i\sin\alpha,$$1 for $$\Sigma_\alpha=\cos\alpha+i\hat\phi_i\tau_i\sin\alpha,$$2, $$\Sigma_\alpha=\cos\alpha+i\hat\phi_i\tau_i\sin\alpha,$$3, and $$\Sigma_\alpha=\cos\alpha+i\hat\phi_i\tau_i\sin\alpha,$$4, respectively. Because $$\Sigma_\alpha=\cos\alpha+i\hat\phi_i\tau_i\sin\alpha,$$5 while $$\Sigma_\alpha=\cos\alpha+i\hat\phi_i\tau_i\sin\alpha,$$6, the magnetic field splits the flavor response, whereas the KMT determinant tends to reduce that splitting (Chatterjee et al., 2011).

The corresponding lattice picture in $$\Sigma_\alpha=\cos\alpha+i\hat\phi_i\tau_i\sin\alpha,$$7-flavor QCD at $$\Sigma_\alpha=\cos\alpha+i\hat\phi_i\tau_i\sin\alpha,$$8 is charge resolved rather than flavor symmetric. The renormalized condensates

$$\Sigma_\alpha=\cos\alpha+i\hat\phi_i\tau_i\sin\alpha,$$9

show magnetic catalysis, with the up condensate growing faster than the down condensate. For $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$00 GeV$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$01, the fitted slopes are

$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$02

while at smaller fields the growth is described by $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$03 with $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$04 and $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$05. Neutral pseudoscalar masses decrease monotonically and then saturate at nonzero values, neutral decay constants increase, and an approximate $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$06 scaling is observed for the $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$07- and $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$08-flavor components of the neutral pion mass, decay constant, and condensates. The correction to the GMOR relation is less than $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$09 for the neutral pion and less than $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$10 for the neutral kaon at $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$11 GeV$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$12 (Ding et al., 2020).

Three-flavor $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$13PT to $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$14 translates this response into a vacuum-free-energy problem. The field-dependent free energy $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$15 generates

$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$16

A technical result is that all two-loop charged-charged vacuum graphs cancel identically, leaving a simplified $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$17 structure. The light-quark condensate shift at $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$18 is reported to agree better with recent $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$19-flavor lattice data than the $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$20 result, while the renormalized magnetization remains uncertain because of the LEC dependence (Adhikari et al., 2023).

5. Flavor content, critical surfaces, and vacuum diagnostics

A direct way to vary the vacuum’s flavor content is to vary $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$21. In a coupled DSE analysis of massless QCD, increasing $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$22 weakens the infrared gluon interaction through quark-loop screening until the scalar quark self-energy disappears at

$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$23

The condensate scales as

$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$24

or equivalently $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$25 with $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$26, which is interpreted as a second-order chiral transition. Above $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$27, the vacuum is chirally symmetric but still retains a gluonic infrared scale, so the postcritical regime is described as walking rather than fully conformal (Chen et al., 16 Feb 2026).

The low Dirac spectrum offers a complementary diagnostic of vacuum flavor content. The spectral density

$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$28

is nearly constant below $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$29 MeV and then rises approximately linearly up to $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$30 MeV. This is interpreted as an infrared regime that is effectively $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$31, with the strange quark decoupled at the deepest scales, followed by a region where the strange sea becomes visible. Fitting NLO PQ$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$32PT yields

$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$33

and

$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$34

so the $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$35 and $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$36 chiral limits are quantitatively different. The same analysis shows that sea-quark masses can be constrained from the spectral density with $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$37 statistical uncertainties (Liang et al., 2021).

At finite temperature, the flavor dependence of the chiral vacuum is encoded in the critical surface. In a $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$38-flavor QM model with FRG, the vacuum order parameters are the light and strange condensates $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$39 and $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$40. A fixed-UV parameter choice can unphysically eliminate spontaneous chiral symmetry breaking in the chiral limit, whereas a fixed-$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$41 scheme is proposed as a heuristic continuation away from the physical point. With anomaly present, the FRG gives a very small first-order region near the three-flavor chiral limit, including $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$42 MeV and a tricritical point on the $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$43-axis near $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$44 MeV, much smaller than the eMFA values $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$45 MeV and $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$46 MeV (Resch et al., 2017). An independent many-flavor lattice method finds a finite critical heavy-flavor parameter,

$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$47

in the two-flavor chiral limit, which is taken to suggest that massless two-flavor QCD is second order rather than first order (Yamada et al., 2016).

Hadronic spectroscopy and current correlators also reflect the vacuum’s chiral-flavor structure. In a three-flavor heavy-baryon chiral model, the vacuum expectation value

$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$48

splits a $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$49 multiplet into a positive-parity $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$50 and a negative-parity $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$51, identified as chiral partners (Kawakami et al., 2019). At the level of precision correlators, three-loop two-flavor $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$52PT for hadronic vacuum polarization computes the low-virtuality flavor-current fluctuation

$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$53

to NNNLO; the amplitude requires elliptic master integrals, and its renormalizability depends on master-integral relations not implied by standard IBP reduction (Lellouch et al., 14 Oct 2025). These are not alternative definitions of the chiral-flavor vacuum, but they are precision probes of it.

6. Bogoliubov vacua, field mixing, and flavor-chiral charges

In QFT treatments of chirality, the left and right charges

$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$54

are not diagonal in the standard massive Dirac basis because the mass term couples $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$55 and $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$56. Diagonalizing them requires a Bogoliubov transformation mixing particle annihilation and antiparticle creation operators. The resulting vacuum is a time-dependent condensate of particle-antiparticle pairs and is orthogonal to the standard Dirac vacuum in the infinite-volume limit. The expectation value of the chiral charge on an initially left-chiral one-particle state reproduces the standard chiral-oscillation formula,

$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$57

but now as a vacuum-sensitive QFT result rather than a first-quantized wave-mechanics formula (Bittencourt et al., 2024).

The explicitly combined flavor-chiral construction for two-flavor mixing proceeds in two stages. First one diagonalizes chirality with a Bogoliubov map; then one rotates in flavor space,

$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$58

The massive chiral vacuum is

$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$59

This state is a condensate of particle-antiparticle pairs with definite mass $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$60, momentum, and helicity. A distinctive result is that the flavor rotation $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$61 leaves the standard mass vacuum invariant, so the vacuum annihilated by flavor-chiral operators coincides with the massive chiral vacuum: $$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$62 Expectation values of chiral-flavor charges recover the first-quantized flavor-chiral oscillation formulas, but with the vacuum structure made explicit (Blasone et al., 13 Jul 2025).

A more abstract algebraic treatment reaches a closely related conclusion from spontaneous symmetry breaking. Starting from

$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$63

dynamical field mixing requires a two-step breaking pattern,

$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$64

with an off-diagonal order parameter in the mass basis,

$$\Omega_{q\bar q}^{\rm vac}=-2N_fN_c\int\frac{d^3p}{(2\pi)^3}E_q\,\theta(\Lambda^2-\vec p^{\,2}),$$65

The resulting flavor vacuum condensate is coherent, entangled, and unitarily inequivalent to the mass vacuum, and Ward-Takahashi identities give the correct number of Nambu-Goldstone modes (Blasone et al., 2018).

This operator-theoretic line of work provides the most literal use of the expression “chiral-flavor vacuum.” It does not refer primarily to a thermodynamic effective potential, but to a ground state defined by inequivalent representations of the CAR algebra, by condensates that are off-diagonal in flavor or chirality, and by the nontrivial fact that the physically relevant charges are diagonal only after Bogoliubov transformation. In that sense, the chiral-flavor vacuum is not simply the vacuum of a theory with chiral and flavor symmetries; it is the vacuum selected once those symmetries are dynamically broken, mixed, or both.