← Back to Part IX: From Perelman to the Standard Model

Instantons, Topology Change, and Perelman Surgery

Self-dual gauge configurations that tunnel between vacua — and their striking parallel with Perelman's surgical resolution of singularities

The BPST Instanton

The Belavin-Polyakov-Schwarz-Tyupkin (BPST) instanton is a self-dual solution to the Yang-Mills equations on $\mathbb{R}^4$ (or equivalently $S^4$by conformal compactification). Self-duality means:

$$F_{\mu\nu} = \tilde{F}_{\mu\nu} \equiv \frac{1}{2}\,\epsilon_{\mu\nu\rho\sigma}\, F^{\rho\sigma}$$

The self-duality condition implies the Yang-Mills equations $D_\mu F^{\mu\nu} = 0$automatically (via the Bianchi identity $D_\mu \tilde{F}^{\mu\nu} = 0$).

Explicit Solution

In the singular gauge, the BPST instanton with size $\rho$ centred at the origin is:

$$A^a_\mu = \frac{2\rho^2}{x^2(x^2 + \rho^2)}\, \bar\eta^a_{\mu\nu}\, x^\nu$$

where $\bar\eta^a_{\mu\nu}$ are the 't Hooft symbols that intertwine the Lorentz and colour indices. The field strength is:

$$F^a_{\mu\nu} = -\frac{4\rho^2}{(x^2 + \rho^2)^2}\, \bar\eta^a_{\mu\nu}$$

Instanton Action

The action of the instanton is topologically quantised:

$$S_{\text{inst}} = -\frac{1}{2g^2}\int \text{tr}(F_{\mu\nu} F^{\mu\nu})\, d^4x = \frac{8\pi^2}{g^2}$$

This is independent of the instanton size $\rho$ and position — reflecting the conformal symmetry of classical Yang-Mills theory in four dimensions. The Bogomolny bound guarantees that self-dual solutions minimise the action in a given topological sector:

$$S \geq \frac{8\pi^2}{g^2}\,|\nu|, \qquad \nu = \frac{g^2}{32\pi^2}\int F^a_{\mu\nu}\tilde{F}^{a\mu\nu}\, d^4x$$

The Bogomolny bound is saturated exactly when $F = \pm\star F$ (self-dual or anti-self-dual). This follows from the algebraic identity:

$$\int |F \mp \star F|^2\, d^4x = 2\int |F|^2\, d^4x \mp 2\int F \wedge F \geq 0$$

Topology Change via Instantons

The instanton interpolates between topologically distinct vacua. At Euclidean time$\tau \to -\infty$, the gauge field is in sector $|n\rangle$; at $\tau \to +\infty$, it is in sector $|n+1\rangle$:

$$\langle n+1\, |\, e^{-H T}\, |\, n\rangle \sim e^{-S_{\text{inst}}} = e^{-8\pi^2/g^2}$$

The tunnelling amplitude is exponentially suppressed at weak coupling, making instantons non-perturbative objects invisible to any finite order of perturbation theory.

In the dilute instanton gas approximation, the vacuum energy depends on$\theta$ as:

$$E(\theta) = -2K\, e^{-8\pi^2/g^2}\cos\theta, \qquad K = \text{det}^\prime(-D^2)^{-1/2}$$

where $K$ is the one-loop fluctuation determinant around the instanton background.

Instanton Moduli Space

A single instanton in $SU(2)$ has 8 collective coordinates (moduli): 4 for position $x_0^\mu$, 1 for size $\rho$, and 3 for global gauge orientation. The instanton measure on moduli space is:

$$d\mu_{\text{inst}} = C_N\, \frac{d^4 x_0\, d\rho}{\rho^5}\, \left(\frac{8\pi^2}{g^2}\right)^{2N} e^{-8\pi^2/g^2}\, \rho^{4N-5}\, d\Omega_{SU(N)}$$

The factor $\rho^{4N-5}$ comes from the one-loop determinant. For $N=2$this gives $\rho^3$, meaning small instantons are suppressed. For$N \geq 3$ with enough light fermions, the integral over $\rho$is infrared-divergent — a signal of strong coupling dynamics.

The 't Hooft Vertex

In the presence of $N_f$ massless quark flavours, each instanton produces a $2N_f$-fermion vertex ('t Hooft vertex):

$$\mathcal{L}_{\text{'t Hooft}} \sim e^{-8\pi^2/g^2}\, \prod_{f=1}^{N_f} \bar\psi_{L,f}\, \psi_{R,f} + \text{h.c.}$$

This vertex violates the axial $U(1)_A$ symmetry, resolving the$U(1)_A$ problem: the $\eta'$ meson is heavy because instantons explicitly break the would-be $U(1)_A$ symmetry.

Parallel with Perelman Surgery

Perelman's proof of the Poincaré conjecture uses surgery to handle singularities that develop under Ricci flow. The parallel with instantons is structural:

Perelman Surgery

Neck-pinch singularity develops

Cut along the neck

Cap with standard spherical caps

Topology changes: one manifold becomes two

Instanton Tunneling

Energy barrier between vacua

Instanton path through the barrier

Self-dual configuration is the “cap”

Topology changes: winding number shifts by 1

Both processes are controlled by an action or entropy functional. For instantons, the suppression is:

$$\text{Instanton:} \quad e^{-S_{\text{inst}}} = e^{-8\pi^2/g^2}$$

For Perelman surgery, the analogous control parameter is the reduced volume:

$$\text{Surgery:} \quad \widetilde{V}(\tau) = \frac{1}{(4\pi\tau)^{n/2}} \int_M e^{-\ell(q,\tau)}\, d\text{vol}_g$$

In both cases, the topology change is not catastrophic but controlled: the action/entropy bounds the rate and ensures the process improves a global functional. Ricci flow decreases the Perelman $\mathcal{W}$-entropy, while instantons minimise the Yang-Mills action in their topological sector.

The Perelman $\mathcal{W}$-entropy that controls the surgery is:

$$\mathcal{W}(g, f, \tau) = \int_M \left[\tau\,(R + |\nabla f|^2) + f - n\right] \frac{e^{-f}}{(4\pi\tau)^{n/2}}\, d\text{vol}_g$$

Under Ricci flow $\partial_t g_{ij} = -2R_{ij}$, this functional is monotonically non-decreasing:

$$\frac{d\mathcal{W}}{dt} = 2\tau \int_M \left|R_{ij} + \nabla_i\nabla_j f - \frac{g_{ij}}{2\tau}\right|^2 \frac{e^{-f}}{(4\pi\tau)^{n/2}}\, d\text{vol}_g \geq 0$$

Compare with the Yang-Mills action, which is minimised by the instanton:

$$S_{\text{YM}} = \frac{1}{2g^2}\int |F|^2\, d^4x \geq \frac{8\pi^2}{g^2}\,|\nu|$$

The Bogomolny bound for Yang-Mills parallels the logarithmic Sobolev inequality that underlies Perelman's monotonicity formula. Both express the idea that topology constrains the minimum of a functional.

Deep Analogy:

Both the instanton and the Perelman surgery are “geometric transitions” — the manifold or field configuration passes through a critical point where the topology changes, controlled by the variational principle. The instanton is the gauge-theory analogue of the neck-pinch, and the self-dual equation $F = \star F$ is the analogue of the soliton equation that governs the cap geometry.

Simulation: BPST Instanton Profiles and Tunneling Rates

Python

script.py77 lines

import numpy as np
import matplotlib.pyplot as plt

# BPST Instanton: action density and topological charge density
# as a function of distance from center, for several instanton sizes rho

r = np.linspace(0.01, 6, 500)

# Instanton sizes
rho_values = [0.5, 1.0, 1.5, 2.0, 3.0]
colors_action = ['#fb7185', '#f97316', '#fbbf24', '#4ade80', '#67e8f4']
colors_charge = ['#fda4af', '#fdba74', '#fde68a', '#86efac', '#a5f3fc']

fig, axes = plt.subplots(1, 3, figsize=(14, 4.5))
fig.patch.set_facecolor('#0a0a1a')

for ax in axes:
    ax.set_facecolor('#0a0a1a')
    ax.tick_params(colors='#94a3b8')
    ax.spines['bottom'].set_color('#334155')
    ax.spines['left'].set_color('#334155')
    ax.spines['top'].set_visible(False)
    ax.spines['right'].set_visible(False)

# Panel 1: Action density |F|^2(r)
# For BPST instanton: |F|^2 = 192 * rho^4 / (r^2 + rho^2)^4
# In 4D, the radial measure is r^3 dr, so the density per unit r is:
# f(r) = |F|^2 * r^3 = 192 * rho^4 * r^3 / (r^2 + rho^2)^4

for rho, c in zip(rho_values, colors_action):
    F_sq = 192 * rho**4 / (r**2 + rho**2)**4
    axes[0].plot(r, F_sq, color=c, linewidth=2, label=f'rho = {rho}')

axes[0].set_xlabel('r', color='#94a3b8')
axes[0].set_ylabel('|F|^2(r)', color='#94a3b8')
axes[0].set_title('Action Density', color='#e2e8f0')
axes[0].legend(fontsize=8, facecolor='#1e293b', edgecolor='#334155', labelcolor='#e2e8f0')
axes[0].set_yscale('log')
axes[0].set_ylim(1e-4, 1e4)

# Panel 2: Topological charge density * r^3 (radial profile in 4D)
for rho, c in zip(rho_values, colors_charge):
    q_r = 192 * rho**4 * r**3 / (r**2 + rho**2)**4
    axes[1].plot(r, q_r, color=c, linewidth=2, label=f'rho = {rho}')

axes[1].set_xlabel('r', color='#94a3b8')
axes[1].set_ylabel('q(r) * r^3', color='#94a3b8')
axes[1].set_title('Radial Charge Profile (4D)', color='#e2e8f0')
axes[1].legend(fontsize=8, facecolor='#1e293b', edgecolor='#334155', labelcolor='#e2e8f0')

# Panel 3: Tunneling amplitude as function of coupling
g_values = np.linspace(0.3, 3.0, 200)
# Instanton suppression: exp(-8*pi^2 / g^2)
tunneling = np.exp(-8 * np.pi**2 / g_values**2)

axes[2].semilogy(g_values, tunneling, color='#c084fc', linewidth=2.5)
axes[2].axvline(x=1.0, color='#fbbf24', linestyle='--', alpha=0.5, label='g = 1')
axes[2].axvline(x=0.5, color='#fb7185', linestyle='--', alpha=0.5, label='g = 0.5 (perturbative)')

# Mark QCD-like coupling
axes[2].fill_betweenx([1e-60, 1], 0.3, 0.6, alpha=0.1, color='#67e8f4', label='Perturbative regime')
axes[2].fill_betweenx([1e-60, 1], 1.5, 3.0, alpha=0.1, color='#fb7185', label='Non-perturbative')

axes[2].set_xlabel('Coupling g', color='#94a3b8')
axes[2].set_ylabel('exp(-8 pi^2 / g^2)', color='#94a3b8')
axes[2].set_title('Instanton Tunneling Rate', color='#e2e8f0')
axes[2].legend(fontsize=7, facecolor='#1e293b', edgecolor='#334155', labelcolor='#e2e8f0', loc='center right')
axes[2].set_ylim(1e-60, 10)

plt.tight_layout()
plt.savefig('output.png', dpi=150, bbox_inches='tight', facecolor='#0a0a1a')
plt.show()
print("Plotted BPST instanton profiles and tunneling rates.")
print("Left: action density peaks at r=0, broader for larger rho.")
print("Center: radial charge profile peaks at r ~ rho (most of the topological charge).")
print("Right: tunneling is exponentially suppressed at weak coupling.")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

← CS & Anomalies Next: Higgs as Bundle →