The published abstract of Abrikosov's 1981 letter on the instanton in the heated quark plasma establishes that, "independently of the temperature, in leading order of perturbation theory the instantons do not interact with the dense plasma of light quarks"; beyond leading order, whether higher-loop finite-temperature corrections produce such an interaction, suppress large instantons, or generate Debye-type screening is left open by the published result. Subsequent literature has settled the one-loop story: the Gross-Pisarski-Yaffe thermal determinant does suppress instantons at finite TT with a Debye-screening interpretation, and at quark chemical potential μ\mu a one-loop Pauli-blocking factor e−Nfμ2ρ2e^{-N_f\mu^2\rho^2} arises from the fermion zero modes (zero-energy solutions of the Dirac equation in the instanton background), as established by Abrikosov's own 1983 paper and adopted in the modern dense-QCD literature. What has never been computed in four-dimensional QCD is the two-loop finite-temperature correction to the caloron effective action and instanton size distribution - the leading beyond-one-loop instanton-medium interaction, with higher orders entering the required answer only where they are parametrically enhanced relative to it - even though modern work repeatedly identifies it as needed: lattice QCD finds the high-temperature instanton density roughly ten times larger than the one-loop dilute-instanton-gas prediction, and axion-cosmology applications depend on precisely this regime. The open problem is to determine that correction.

Background

An instanton is a self-dual solution of the classical Euclidean Yang-Mills equations (the Belavin-Polyakov-Schwartz-Tyupkin solution) carrying topological charge Q=1Q=1; it mediates quantum tunneling between gauge vacua of different winding number, and the fermion zero modes (zero-energy solutions of the Dirac equation) localized in its core generate the anomalous multi-fermion interaction known as the 't Hooft vertex, the microscopic origin of the U(1)AU(1)_A anomaly and of anomalous chiral processes in electroweak and QCD plasmas. At finite temperature TT, Euclidean time is compactified to a circle of circumference β=1/T\beta=1/T with periodic boundary conditions for the gauge fields and antiperiodic boundary conditions for the quarks, so that fields are expanded in discrete thermal frequencies (Matsubara frequencies), ωn=2πnT\omega_n=2\pi nT for bosons and (2n+1)πT(2n+1)\pi T for fermions; the instanton generalizes to the Harrington-Shepard caloron, a periodic self-dual solution. In the dilute instanton gas approximation (DIGA), instantons of size ρ\rho are treated as a non-interacting gas with a size distribution n(ρ,T)∝d(ρ)exp⁡[−Seff(ρ,T)]n(\rho,T)\propto d(\rho)\exp[-S_{\mathrm{eff}}(\rho,T)], where Seff=8π2/g2(ρ)+ΔS(ρ,T)S_{\mathrm{eff}}=8\pi^2/g^2(\rho)+\Delta S(\rho,T) is the classical action (with gg the QCD gauge coupling, running with the scale ρ\rho) plus the medium correction ΔS\Delta S from quantum fluctuations. Integrating n(ρ,T)n(\rho,T) over ρ\rho gives the high-temperature topological susceptibility χtop(T)\chi_{\mathrm{top}}(T), the variance of the total topological charge, which is the key input to the cosmological history of the QCD axion (a hypothetical light pseudoscalar particle).

The one-loop medium correction at finite temperature was computed by Pisarski and Yaffe and summarized by Gross, Pisarski and Yaffe:

S(λ)=λ2(2Nc+Nf)/3+12A(λ)(1+(Nc−Nf)/6),λ=πρT,S(\lambda)=\lambda^2(2N_c+N_f)/3+12A(\lambda)(1+(N_c-N_f)/6),\quad \lambda=\pi\rho T,

where NcN_c is the number of colors, NfN_f the number of quark flavors, and A(λ)A(\lambda) a known function with A(λ)≈−λ2/36A(\lambda)\approx-\lambda^2/36 at small λ\lambda. This factor exponentially suppresses large instantons at high TT. Because instantons carry long-range color-electric fields, the dominant ρ2mD2\rho^2m_D^2-type piece is interpretable as Debye screening - the screening of static electric color fields in the plasma, with electric screening mass mD2=(Nc/3+Nf/6)g2T2m_D^2=(N_c/3+N_f/6)g^2T^2 at one loop. Korthals Altes and Sastre later rewrote this one-loop thermal determinant in compact form and clarified its infrared structure: the gluonic (periodic) determinant is logarithmically infrared divergent at large ρT\rho T, while the fermionic (antiperiodic) determinant is infrared finite.

Abrikosov's 1981 letter approached the same medium from a different decomposition of the thermal fermion propagator: in a thermal ensemble at temperature TT and quark chemical potential μ\mu, the free propagator splits into a vacuum part and an on-shell medium part whose phase-space weights are Fermi-Dirac occupation factors. Computing this medium part in the caloron background at one loop in SU(2) QCD with one massless quark flavor (plus the first nontrivial small-quark-mass term), he found that the quark number density equals that of the plasma without the caloron and concluded that instantons do not interact with the on-shell quark plasma at one loop, independently of TT. The analysis is restricted to one loop and to the dilute-gas setting; since the published result holds only "in leading order of perturbation theory", whether two-loop and higher finite-temperature corrections generate an instanton-plasma interaction remains open.

The subsequent arc frames the open question. Abrikosov's own two-loop computation in a cold dense plasma at T=0T=0 found logarithmically enhanced interaction terms that decrease the effective action as the instanton size grows. His 1983 long paper established that at μ≠0\mu\neq 0 the fermion zero modes produce a one-loop Pauli-blocking suppression (the filled Fermi sea blocks occupation of the zero-mode fermion states) ΔS=Nfμ2ρ2\Delta S=N_f\mu^2\rho^2, with total small-ρT\rho T action 8π2/g2(ρ)+Nfμ2ρ2+(π2/18)ρ2T2(11Nc+7Nf−6)8\pi^2/g^2(\rho)+N_f\mu^2\rho^2+(\pi^2/18)\rho^2T^2(11N_c+7N_f-6), the last term being the small-λ\lambda limit of the Gross-Pisarski-Yaffe factor. Morris, Ross and Sachrajda computed the two-loop instanton density at T=0T=0 and found it renormalization-scheme dependent. Modern work - Boccaletti and Nogradi's precision re-evaluation of the one-loop finite-TT distribution, Pisarski and Rennecke's observation that the lattice instanton density exceeds the one-loop dilute-gas result by a prefactor of order ten that a complete two-loop thermal computation might ameliorate, and Högl and Moore's one-loop caloron-gas treatment of the finite-TT topological susceptibility with heavy quarks - all operates at one loop. Whether the instanton-quark-plasma interaction survives beyond one loop at finite temperature, and with what sign, remains unknown.

Problem Statement

Abrikosov (1981) concluded, per the published abstract, that "independently of the temperature, in leading order of perturbation theory the instantons do not interact with the dense plasma of light quarks"; the full text derives this at one loop in SU(2) QCD with one massless quark flavor, from the on-shell medium part of the thermal fermion propagator in a caloron background leaving the quark number density unchanged. Beyond leading order, whether two-loop and higher finite-temperature corrections generate a net interaction between instantons and the heated plasma is left open. The one-loop part of the question has since been settled - the Gross-Pisarski-Yaffe determinant suppresses instantons at finite TT (with a Debye-screening interpretation of the dominant ρ2mD2\rho^2m_D^2-type term, where ρ\rho is the instanton size and mDm_D the electric color-screening mass), and a finite quark chemical potential μ\mu adds the Pauli-blocking factor e−Nfμ2ρ2e^{-N_f\mu^2\rho^2} (NfN_f being the number of quark flavors) - but the beyond-one-loop part has never been computed in four-dimensional QCD.

Determine the two-loop finite-temperature correction to the caloron effective action and to the instanton size distribution n(ρ,T)n(\rho,T) of a thermal quark-gluon medium; higher-order terms enter the required answer only where they are parametrically enhanced relative to the two-loop result. This determination includes the interaction between the caloron and the thermal quark medium, understood as the (ρ,T)(\rho,T)-dependent part of the two-loop correction to the medium term ΔS(ρ,T)\Delta S(\rho,T) of the effective action, and must establish whether that interaction survives beyond one loop and whether it suppresses large instantons, acting as Debye-type screening beyond the one-loop factor, or instead decreases the effective action with growing instanton size ρ\rho, as Abrikosov's two-loop analysis of the cold dense plasma suggests.

The answer should be stated in a definite renormalization scheme together with its scheme dependence (the T=0T=0 two-loop density of Morris-Ross-Sachrajda is already scheme dependent), should specify its domain of validity in ρT\rho T (the semiclassical dilute-gas regime versus the infrared-sensitive large-instanton regime), and should be consistent with the known one-loop infrared structure (logarithmic divergence of the periodic gluon determinant, finiteness of the antiperiodic fermion determinant). The two-loop correction to ΔS(ρ,T)\Delta S(\rho,T) and to n(ρ,T)n(\rho,T), together with the sign of the interaction, is subject to these scheme, domain-of-validity, and infrared-consistency conditions.

References

  1. A. A. Abrikosov, Instanton in the heated quark plasma, Phys. Lett. B 99 (1981) 421-424. DOI: 10.1016/0370-2693(81)90561-X; https://doi.org/10.1016/0370-2693(81)90561-X
  2. A. A. Abrikosov, Instantons in the quark plasma the two-loop approximation, Nucl. Phys. B 182 (1981) 441-461. DOI: 10.1016/0550-3213(81)90129-2; https://doi.org/10.1016/0550-3213(81)90129-2
  3. A. A. Abrikosov, Instantons in the quark plasma, Nucl. Phys. B 218 (1983) 459-483. DOI: 10.1016/0550-3213(83)90375-9; https://doi.org/10.1016/0550-3213(83)90375-9
  4. A. A. Belavin, A. M. Polyakov, A. S. Schwartz and Yu. S. Tyupkin, Pseudoparticle solutions of the Yang-Mills equations, Phys. Lett. B 59 (1975) 85-87. DOI: 10.1016/0370-2693(75)90163-x; https://doi.org/10.1016/0370-2693(75)90163-x
  5. B. J. Harrington and H. K. Shepard, Periodic Euclidean solutions and the finite-temperature Yang-Mills gas, Phys. Rev. D 17 (1978) 2122-2125. DOI: 10.1103/PhysRevD.17.2122; https://doi.org/10.1103/PhysRevD.17.2122
  6. R. D. Pisarski and L. G. Yaffe, The density of instantons at finite temperature, Phys. Lett. B 97 (1980) 110-112. DOI: 10.1016/0370-2693(80)90559-6; https://doi.org/10.1016/0370-2693(80)90559-6
  7. D. J. Gross, R. D. Pisarski and L. G. Yaffe, QCD and Instantons at Finite Temperature, Rev. Mod. Phys. 53 (1981) 43-80. DOI: 10.1103/RevModPhys.53.43; https://doi.org/10.1103/RevModPhys.53.43
  8. V. Baluni, Screening of instantons in the quark gas, Phys. Lett. B 106 (1981) 491-496. DOI: 10.1016/0370-2693(81)90264-1; https://doi.org/10.1016/0370-2693(81)90264-1
  9. T. R. Morris, D. A. Ross and C. T. Sachrajda, Higher-order quantum corrections in the presence of an instanton background field, Nucl. Phys. B 255 (1985) 115-148. DOI: 10.1016/0550-3213(85)90131-2; https://doi.org/10.1016/0550-3213(85)90131-2
  10. I. I. Balitsky and V. M. Braun, Instanton-induced contributions of fractional twist to the cross section of hard gluon-gluon scattering in QCD, Phys. Rev. D 47 (1993) 1879-1888. DOI: 10.1103/PhysRevD.47.1879; https://doi.org/10.1103/PhysRevD.47.1879
  11. T. Schäfer, Instanton effects in QCD at high baryon density, Phys. Rev. D 65 (2002) 094033. DOI: 10.1103/PhysRevD.65.094033; arXiv:hep-ph/0201189; https://arxiv.org/abs/hep-ph/0201189
  12. C. P. Korthals Altes and A. Sastre, Thermal instanton determinant in compact form, Phys. Rev. D 90 (2014) 125002. DOI: 10.1103/PhysRevD.90.125002; arXiv:1406.6911; https://arxiv.org/abs/1406.6911
  13. R. D. Pisarski and F. Rennecke, Multi-instanton contributions to anomalous quark interactions, Phys. Rev. D 101 (2020) 114019 (preprint version: arXiv:1910.14052, The Myriad Uses of Instantons). DOI: 10.1103/PhysRevD.101.114019; https://arxiv.org/abs/1910.14052
  14. A. Boccaletti and D. Nogradi, The semi-classical approximation at high temperature revisited, JHEP 03 (2020) 045. DOI: 10.1007/JHEP03(2020)045; arXiv:2001.03383; https://arxiv.org/abs/2001.03383
  15. Y. Fujimoto, Interplay between the weak-coupling results and the lattice data in dense QCD, arXiv:2408.12514; https://arxiv.org/abs/2408.12514
  16. B. Högl and G. D. Moore, Finite-T topological susceptibility with heavy quarks, JHEP 12 (2025) 089. DOI: 10.1007/JHEP12(2025)089; arXiv:2501.11990; https://arxiv.org/abs/2501.11990