Nonperturbative condensate input of the QCD operator product expansion: four-quark condensates and the accuracy of vacuum saturation

In his 1976 study of effective quark masses in the chiral limit, Politzer used the operator product expansion (OPE) to show that spontaneous chiral-symmetry breaking enters short-distance quark two-point functions only through nonperturbative vacuum matrix elements, and acknowledged that these condensates are external input: the two-quark condensate was not computed in the framework, and four-quark condensates were estimated only through an uncontrolled factorization (vacuum saturation) hypothesis. Later literature resolved part of this agenda: the two-quark condensate is now precisely determined by lattice QCD, and renormalon theory clarified the asymptotic character of the OPE at finite momentum. However, the surviving core remains open: the dimension-six four-quark condensates of real QCD (with color count $N_c=3$) have no controlled first-principles determination, and the degree of violation of vacuum saturation, conventionally encoded in factors $\kappa$, is unknown, with channel-by-channel estimates scattered from $\kappa\approx1$ to $\kappa\approx7$. The problem asks for a controlled determination of the four-quark condensate input of the OPE at phenomenologically relevant scales, and a quantification of factorization violation, in real QCD.

Background

In quantum chromodynamics (QCD), the theory of the strong interaction, the light quarks ($u$, $d$, $s$) carry small 'current' masses that appear in the Lagrangian, yet inside hadrons they behave as if they had 'constituent' masses of a few hundred MeV. The difference signals spontaneous chiral symmetry breaking: in the chiral limit (vanishing current quark masses) the Lagrangian is nearly invariant under a symmetry that the quantum vacuum does not respect. A measurable signature is the quark condensate $v(M)=\langle 0|\bar\psi\psi|0\rangle$, the vacuum expectation value of the quark bilinear $\bar\psi\psi$ renormalized at scale $M$, which is nonzero and of order $\Lambda_{\mathrm{QCD}}^3$, where $\Lambda_{\mathrm{QCD}}$ is the intrinsic mass scale of the strong interaction.

Politzer's 1976 paper asked how the two sources of quark mass — explicit mass terms versus dynamical, spontaneous-symmetry-breaking effects — are separated at short distances. The tool is the operator product expansion (OPE): the statement that a product of local operators at short separation expands into a sum of local operators whose momentum-dependent coefficient functions are perturbatively computable at large spacelike momentum $Q$, multiplied by vacuum expectation values ('condensates') of those operators, which carry the nonperturbative physics. Applied to the quark two-point function, this defines an effective, scale-dependent quark mass $m(Q)$ whose leading spontaneous contribution behaves as $m_{\mathrm{spont}}(Q)\sim v(M)(\log Q^2)^{-15/27}/Q^2$ — power-suppressed relative to the running explicit mass $m_0(Q)\sim m_0(M)(\log Q^2)^{-12/27}$. Explicit masses therefore dominate at large $Q$, and ratios of effective quark masses approach ratios of bare masses even in the presence of spontaneous chiral symmetry breaking. Politzer also proposed to represent current-current spectral functions by perturbative graphs with $m(Q)$ substituted for the quark masses, as a practical mnemonic.

The paper itself acknowledged two gaps in this construction, which became the agenda of a whole subfield. First, the condensates are not computable within the OPE; they must be supplied from outside (experiment, models, or nonperturbative calculation). Second, beyond leading order, four-quark operators $\bar\psi\Gamma\psi\,\bar\psi\Gamma\psi$ (mass dimension six) appear, and their vacuum matrix elements were estimated by the factorization or vacuum saturation hypothesis,

$$\langle 0|\bar\psi\Gamma_1\psi\,\bar\psi\Gamma_2\psi|0\rangle \;\simeq\; \langle 0|\bar\psi\psi|0\rangle^{2}\times(\text{known spin-color factor}),$$

obtained by inserting a complete set of states between the two bilinears and keeping only the vacuum contribution. Shifman, Vainshtein and Zakharov (1979) systematized the OPE-plus-condensates idea into QCD sum rules, the standard framework for extracting hadron properties from short-distance expansions, in which dimension-six terms enter prominently (for example in nucleon sum rules). Factorization is exact in the limit of a large number of colors, $N_c\to\infty$, where vacuum expectation values of products of color-singlet operators factorize; but real QCD has $N_c=3$, and the accuracy of the approximation for four-quark condensates — conventionally encoded in correction factors $\kappa$ multiplying the factorized estimate — has never been established. A second refinement came from renormalon theory (Beneke, 1999): the perturbative coefficients of the OPE grow factorially, so the OPE is an asymptotic rather than convergent expansion, and the values of high-dimension condensates carry ambiguities that cancel against those of the perturbative series. This sharpened the source paper's observation that the expansion is controlled only at large $Q^2$, and that higher-dimension operators and nonperturbative effects limit its local predictive power at low momenta.

Meanwhile the two-quark condensate itself has been determined precisely: through the Gell-Mann–Oakes–Renner (GMOR) relation, which ties the product of the light-quark mass and condensate to the pion mass and decay constant (already available to the source paper), and from first-principles lattice QCD (QCD discretized and simulated numerically) — via the Banks–Casher relation, which connects the condensate to the density of near-zero eigenvalues of the Dirac operator (Fukaya et al., 2010), and via heavy-light current-current correlators, i.e. correlation functions of currents coupling one heavy and one light quark (Davies et al., 2019). By contrast, no controlled determination of the vacuum four-quark condensates exists: the 2019 review by Gubler and Satow records that none of the four-quark condensates is 'well constrained in any meaningful way', that vacuum saturation is 'the only method presently known to obtain a concrete numerical value for them', and that 'it is not known to what degree it is violated in real QCD with $N_c=3$', with channel-by-channel $\kappa$ estimates ranging from close to 1 to $2$–$3$, up to $6$, and even $\sim7$ for strange quarks.

Problem Statement

Determine, with controlled and quantified uncertainties, the nonperturbative vacuum matrix elements of the dimension-six four-quark operators that enter the QCD operator product expansion of current correlation functions at the momentum scales used in phenomenology, at a stated renormalization scheme and scale, and thereby establish the accuracy of the vacuum saturation (factorization) hypothesis for four-quark condensates in real QCD with color count $N_c=3$.

This is the surviving core of the acknowledged limitations of Politzer (1976): the two-quark condensate part of that input problem has since been resolved by lattice QCD and GMOR-type determinations, while the four-quark condensates, and the factorization assumption used to estimate them, remain undetermined. The objective has two parts:

(a) a determination of the four-quark condensate combinations that enter standard OPE analyses, at a stated renormalization scheme and scale, with a complete error budget covering statistical, quark-mass, and operator-renormalization/mixing uncertainties, plus discretization and finite-volume uncertainties for lattice determinations;

(b) the corresponding deviation factors $\kappa$ from vacuum saturation for those same condensate combinations, with uncertainties, replacing the current practice of treating $\kappa$ as a free or assumed parameter.

The common objective is to convert the nonperturbative condensate input of the OPE from an uncontrolled assumption into controlled knowledge; (a) and (b) are logically linked and mutually supporting, not independent problems.

References

1. H. D. Politzer, 'Effective quark masses in the chiral limit', Nuclear Physics B 117 (1976) 397-406. DOI: 10.1016/0550-3213(76)90405-3, https://doi.org/10.1016/0550-3213(76)90405-3
2. M. Gell-Mann, R. J. Oakes, B. Renner, 'Behavior of current divergences under SU(3) x SU(3)', Physical Review 175 (1968) 2195-2199. DOI: 10.1103/PhysRev.175.2195, https://doi.org/10.1103/PhysRev.175.2195
3. M. A. Shifman, A. I. Vainshtein, V. I. Zakharov, 'QCD and resonance physics. Theoretical foundations', Nuclear Physics B 147 (1979) 385-447. DOI: 10.1016/0550-3213(79)90022-1, https://doi.org/10.1016/0550-3213(79)90022-1
4. V. A. Novikov, M. A. Shifman, A. I. Vainshtein, M. B. Voloshin, V. I. Zakharov, 'Use and misuse of QCD sum rules, factorization and related topics', Nuclear Physics B 237 (1984) 525-552. DOI: 10.1016/0550-3213(84)90006-3, https://doi.org/10.1016/0550-3213(84)90006-3
5. K. G. Chetyrkin, A. A. Pivovarov, 'Vacuum saturation hypothesis and QCD sum rules', Nuovo Cimento A 100 (1988) 899-906. DOI: 10.1007/BF02789008, https://doi.org/10.1007/BF02789008; arXiv: hep-ph/0105093, https://arxiv.org/abs/hep-ph/0105093
6. S. Leupold, 'Factorization and non-factorization of in-medium four-quark condensates', Physics Letters B 616 (2005) 203-207. DOI: 10.1016/j.physletb.2005.04.046, https://doi.org/10.1016/j.physletb.2005.04.046; arXiv: hep-ph/0502061, https://arxiv.org/abs/hep-ph/0502061
7. R. Thomas, S. Zschocke, B. Kämpfer, 'Evidence for in-medium changes of four-quark condensates', Physical Review Letters 95 (2005) 232301. DOI: 10.1103/PhysRevLett.95.232301, https://doi.org/10.1103/PhysRevLett.95.232301; arXiv: hep-ph/0510156, https://arxiv.org/abs/hep-ph/0510156
8. S. Schael et al. (ALEPH Collaboration), 'Branching ratios and spectral functions of tau decays: Final ALEPH measurements and physics implications', Physics Reports 421 (2005) 191-284. DOI: 10.1016/j.physrep.2005.06.007, https://doi.org/10.1016/j.physrep.2005.06.007; arXiv: hep-ex/0506072, https://arxiv.org/abs/hep-ex/0506072
9. R. Thomas, T. Hilger, B. Kämpfer, 'Four-quark condensates in nucleon QCD sum rules', Nuclear Physics A 795 (2007) 19-46. DOI: 10.1016/j.nuclphysa.2007.08.012, https://doi.org/10.1016/j.nuclphysa.2007.08.012; arXiv: 0704.3004, https://arxiv.org/abs/0704.3004
10. E. Shintani, S. Aoki, S. Hashimoto, H. Matsufuru, J. Noaki, T. Onogi, N. Yamada, 'Strong coupling constant and four-quark condensates from vacuum polarization functions with dynamical overlap fermions', Proceedings of Science, PoS(LATTICE 2008) (2009) 134. DOI: 10.22323/1.066.0134, https://doi.org/10.22323/1.066.0134
11. A. Gómez Nicola, J. R. Peláez, J. Ruiz de Elvira, 'Nonfactorization of four-quark condensates at low energies within chiral perturbation theory', Physical Review D 82 (2010) 074012. DOI: 10.1103/PhysRevD.82.074012, https://doi.org/10.1103/PhysRevD.82.074012; arXiv: 1005.4370, https://arxiv.org/abs/1005.4370
12. H. Fukaya, S. Aoki, S. Hashimoto, T. Kaneko, J. Noaki, T. Onogi, N. Yamada (JLQCD Collaboration), 'Determination of the chiral condensate from 2+1-flavor lattice QCD', Physical Review Letters 104 (2010) 122002. DOI: 10.1103/PhysRevLett.104.122002, https://doi.org/10.1103/PhysRevLett.104.122002; arXiv: 0911.5555, https://arxiv.org/abs/0911.5555
13. M. Beneke, 'Renormalons', Physics Reports 317 (1999) 257-286. DOI: 10.1016/S0370-1573(98)00130-6, https://doi.org/10.1016/S0370-1573(98)00130-6; arXiv: hep-ph/9807443, https://arxiv.org/abs/hep-ph/9807443
14. F. L. Braghin, F. S. Navarra, 'Factorization breaking of four-quark condensates in the Nambu-Jona-Lasinio model', Physical Review D 91 (2015) 074008. DOI: 10.1103/PhysRevD.91.074008, https://doi.org/10.1103/PhysRevD.91.074008; arXiv: 1404.4094, https://arxiv.org/abs/1404.4094
15. P. Gubler, D. Satow, 'Recent progress in QCD condensate evaluations and sum rules', Progress in Particle and Nuclear Physics 106 (2019) 1-58. DOI: 10.1016/j.ppnp.2019.02.005, https://doi.org/10.1016/j.ppnp.2019.02.005; arXiv: 1812.00385, https://arxiv.org/abs/1812.00385
16. C. T. H. Davies, K. Hornbostel, J. Komijani, J. Koponen, G. P. Lepage, A. T. Lytle, C. McNeile, 'Determination of the quark condensate from heavy-light current-current correlators in full lattice QCD', Physical Review D 100 (2019) 034506. DOI: 10.1103/PhysRevD.100.034506, https://doi.org/10.1103/PhysRevD.100.034506; arXiv: 1811.04305, https://arxiv.org/abs/1811.04305
17. Zhi-Gang Wang, 'Does vacuum saturation work for the higher dimensional vacuum condensates in the QCD sum rules?', International Journal of Modern Physics A 36 (2021) 2150246. DOI: 10.1142/S0217751X21502468, https://doi.org/10.1142/S0217751X21502468; arXiv: 2107.09328, https://arxiv.org/abs/2107.09328
18. M. Ding, C. D. Roberts, S. M. Schmidt, 'Emergence of hadron mass and structure', Particles 6 (2023) 57-120. DOI: 10.3390/particles6010004, https://doi.org/10.3390/particles6010004; arXiv: 2211.07763, https://arxiv.org/abs/2211.07763
19. J. Kim, P. Gubler, S. H. Lee, 'Phi meson properties in nuclear matter from QCD sum rules with chirally separated four-quark condensates', Physical Review D 105 (2022) 114053. DOI: 10.1103/PhysRevD.105.114053, https://doi.org/10.1103/PhysRevD.105.114053; arXiv: 2204.11440, https://arxiv.org/abs/2204.11440
20. M. Mariño, R. Miravitllas, 'Trans-series from condensates', SciPost Physics 18 (2025) 101. DOI: 10.21468/SciPostPhys.18.3.101, https://doi.org/10.21468/SciPostPhys.18.3.101; arXiv: 2402.19356, https://arxiv.org/abs/2402.19356
21. Y. Liu, M. Mariño, 'Trans-series from condensates in the non-linear sigma model', SciPost Physics 20 (2026) 156. DOI: 10.21468/SciPostPhys.20.6.156, https://doi.org/10.21468/SciPostPhys.20.6.156; arXiv: 2507.02605, https://arxiv.org/abs/2507.02605