Dilute Instanton Ensembles and Theta Dependence

A dilute gas of charge $+1$ instantons and charge $-1$ anti-instantons turns the sector phase $e^{i\theta Q}$ into a periodic vacuum energy. The familiar cosine is not a universal strong-coupling theorem: it follows when individual event weights are finite and the ensemble is dilute enough that connected interactions can be neglected.

Required background. Instanton measures, zero modes, and determinants supplies the one-event weight. Theta parameters, theta states, and sector sums supplies the phase $e^{i\theta Q}$ and its periodicity.

Helpful background. Topological susceptibility and vacuum response defines the response observable extracted below.

Independent topological events

Let $\mathcal K$ be the fully integrated one-instanton density per Euclidean spacetime volume after every bosonic and fermionic zero mode relevant to the vacuum amplitude has been saturated. Assume:

• the instanton and anti-instanton weights are equal at $\theta=0$;
• their charges are $+1$ and $-1$;
• the mean separation is large compared with the event core;
• interactions and overlap corrections are negligible at the desired order.

In a Euclidean volume $\mathcal V$, a configuration with $n_+$ instantons and $n_-$ anti-instantons has $Q=n_+-n_-$. Independent center integrations and indistinguishability give

$$\frac{Z(\theta)}{Z_{\rm pert}} = \sum_{n_+,n_-=0}^{\infty} \frac{(\mathcal K\mathcal V)^{n_++n_-}}{n_+!\,n_-!} e^{i\theta(n_+-n_-)}.$$

The sums factorize:

$$\frac{Z(\theta)}{Z_{\rm pert}} = \exp\!\left(\mathcal K\mathcal V e^{i\theta}\right) \exp\!\left(\mathcal K\mathcal V e^{-i\theta}\right) = \exp\!\left(2\mathcal K\mathcal V\cos\theta\right).$$

The product in the middle line is meant literally; combining its exponents gives the final expression. Dividing by the $\theta=0$ partition function removes the common normalization:

$$\frac{Z(\theta)}{Z(0)} = \exp\!\left[2\mathcal K\mathcal V(\cos\theta-1)\right].$$

Thus the vacuum-energy density obeys

$$\mathcal E(\theta)-\mathcal E(0) = 2\mathcal K(1-\cos\theta),$$

and the dilute-gas susceptibility is

$$\chi = \left.\frac{\partial^2\mathcal E}{\partial\theta^2}\right|_{\theta=0} =2\mathcal K.$$

The derivation is exact under the independent-event assumptions. Whether those assumptions hold is a separate dynamical question. The corresponding gauge-instanton sum is developed in Mariño 2015, § 4.5, pp. 139–146.

Probability and cumulants

At $\theta=0$, $n_+$ and $n_-$ are independent Poisson variables with mean $\mathcal K\mathcal V$. Their difference $Q$ therefore follows a Skellam distribution,

$$P(Q) = e^{-2\mathcal K\mathcal V} I_{|Q|}(2\mathcal K\mathcal V),$$

where $I_n$ is a modified Bessel function. The cumulant-generating function is

$$\log\left\langle e^{i\theta Q}\right\rangle = 2\mathcal K\mathcal V(\cos\theta-1).$$

Consequently,

$$\frac{\langle Q^2\rangle_c}{\mathcal V}=2\mathcal K, \qquad \frac{\langle Q^4\rangle_c}{\mathcal V}=2\mathcal K,$$

while derivatives with respect to $i\theta$ alternate signs in the Taylor coefficients of $\mathcal E(\theta)$. Equality of all even charge cumulants per volume is a diagnostic of the ideal unit-charge gas, not a general property of gauge theory.

Where the one-event density comes from

In a weakly coupled regime, $\mathcal K$ is obtained by integrating the one-instanton measure over its non-position moduli. The construction is displayed in the moduli-to-measure chain. The underlying one-instanton weight and its fermion factors were derived in ‘t Hooft 1976, §§ III–VII, pp. 3436–3449. For four-dimensional Yang–Mills theory on $\mathbb R^4$, the size integration typically reaches $\rho\Lambda\sim1$, so $\mathcal K$ is not calculable by simply extending the one-loop BPST measure through the infrared.

Fermion zero modes can eliminate the vacuum contribution altogether. If at least one Dirac fermion is exactly massless and no external fermion fields are inserted, the one-event Grassmann integrals are unsaturated, so $\mathcal K_{\rm vac}=0$. With small masses, the measure contains $\prod_f(m_f\rho)$. This is consistent with removing $\theta$ by an anomalous axial change of variables when a quark is massless; it is not captured by inserting a mass-independent $\mathcal K$ into the cosine.

The comparison of canonical saddle calculations places the loop parameter, diluteness parameter, and evidence status side by side. Here the two essential controls are $g^2(1/\rho)\ll1$ over the contributing sizes and a small packing fraction.

Interactions and higher harmonics

When events interact, the partition function is organized by connected clusters rather than independent Poisson factors. A neutral instanton–anti-instanton pair changes the $\theta$-independent normalization, while a connected cluster of net charge $k$ contributes a harmonic proportional to $\cos(k\theta)$ when CP is unbroken at $\theta=0$:

$$\mathcal E(\theta)-\mathcal E(0) = \sum_{k\ge1}c_k\,[1-\cos(k\theta)].$$

The ideal gas retains only $c_1=2\mathcal K$. Higher harmonics can arise from correlated clusters or elementary fractional events under different boundary conditions. Periodicity alone does not determine the coefficients, their signs, or the convergence of this expansion. The finite-temperature instanton gas and its breakdown regimes are reviewed from the semiclassical measure in Gross, Pisarski, and Yaffe 1981, pp. 43–80.

The order of limits matters. The thermodynamic limit $\mathcal V\to\infty$ is taken after constructing the intensive cluster expansion. At finite volume, $Z(\theta)$ is analytic. Nonanalytic branch changes can emerge only in an infinite-volume limit and cannot be inferred from the elementary cosine without additional dynamics.

Common pitfalls

Calling the cosine universal. It is the leading result of an ideal unit-charge gas. Strongly coupled Yang–Mills theory need not have a calculable BPST fugacity or a single cosine vacuum energy.

Ignoring fermion saturation. The bosonic action does not determine whether an instanton contributes to $Z$. Unsaturated Grassmann zero modes make the vacuum contribution vanish.

Using event number as the dilution test. The mean number grows with volume. Diluteness is a local packing condition comparing event size and mean separation.

Exercises

1. Derive the ideal-gas susceptibility directly from charge fluctuations at $\theta=0$.

Solution

For independent Poisson variables, $\operatorname{Var}(n_\pm)=\mathcal K\mathcal V$. Since $Q=n_+-n_-$,

$$\langle Q^2\rangle_c =\operatorname{Var}(n_+)+\operatorname{Var}(n_-) =2\mathcal K\mathcal V.$$

Dividing by $\mathcal V$ gives $\chi=2\mathcal K$, in agreement with two derivatives of the vacuum energy.

2. Suppose elementary events have charges $\pm1/N$, equal fugacity $\zeta$, and are independently summed within a compactified effective theory. What local harmonic do they generate, and why is this not by itself a $2\pi N$-periodic microscopic theory?

Solution

The independent sum gives $2\zeta[1-\cos(\theta/N)]$ on one branch. A microscopic $2\pi$-periodic theory must also include the branch relabeling or global constraint that restores $2\pi$ periodicity. The local fractional-event expression is valid only with its compactification and branch data specified.

References

• Gross, David J., Robert D. Pisarski, and Laurence G. Yaffe. “QCD and Instantons at Finite Temperature.” Reviews of Modern Physics 53 (1981): 43–80. DOI.
• Mariño, Marcos. Instantons and Large N: An Introduction to Non-Perturbative Methods in Quantum Field Theory. Cambridge: Cambridge University Press, 2015. DOI.
• ‘t Hooft, Gerard. “Computation of the Quantum Effects Due to a Four-Dimensional Pseudoparticle.” Physical Review D 14 (1976): 3432–3450; erratum 18 (1978): 2199. DOI.