Index Theorems

The Atiyah-Singer index theorem is one of the greatest mathematical achievements of the 20th century. It computes the analytical index of an elliptic operator (the difference between the dimensions of its kernel and cokernel) in terms of topological data (characteristic classes). This result unifies the Gauss-Bonnet theorem, the Hirzebruch signature theorem, and the Riemann-Roch theorem, and explains chiral anomalies in quantum field theory.

Historical Context

Michael Atiyah and Isadore Singer announced their index theorem in 1963, providing a formula for the index of any elliptic differential operator on a compact manifold in terms of characteristic classes. The theorem unified several classical results: the Gauss-Bonnet theorem, the Hirzebruch signature theorem (1954), and the Riemann-Roch theorem (1850s/1954).

The connection to physics came in the late 1970s and 1980s. The chiral anomaly in quantum field theory (Adler-Bell-Jackiw, 1969) was recognized as a consequence of the index theorem for the Dirac operator. Anomaly cancellation in the Standard Model and in string theory (Green-Schwarz, 1984) relies fundamentally on index-theoretic calculations.

Atiyah and Singer received the Abel Prize in 2004 for their index theorem, which Atiyah called "the single most important theorem in mathematics of the 20th century."

Derivation 1: The Analytical Index

An elliptic operator $D: \Gamma(E) \to \Gamma(F)$ between sections of vector bundles over a compact manifold has finite-dimensional kernel and cokernel. The analytical index is:

$\text{ind}(D) = \dim\ker(D) - \dim\ker(D^*) = \dim\ker(D) - \dim\text{coker}(D)$

Key Properties

Topological invariance: The index is unchanged under continuous deformations of $D$

Stability: Adding a compact (lower-order) perturbation does not change the index

Additivity: $\text{ind}(D_1 \oplus D_2) = \text{ind}(D_1) + \text{ind}(D_2)$

Physical meaning: For the Dirac operator coupled to a gauge field,$\ker(D)$ contains the zero-energy fermion modes of one chirality, and $\ker(D^*)$ the other. The index counts the net number of chiral zero modes—this is the chiral anomaly.

Derivation 2: The Atiyah-Singer Index Theorem

For an elliptic operator $D$ on a compact manifold $M$, the index theorem equates the analytical index to the topological index:

$\boxed{\text{ind}(D) = \int_M \text{ch}(\sigma(D)) \wedge \text{Td}(TM \otimes \mathbb{C})}$

where $\text{ch}(\sigma(D))$ is the Chern character of the symbol of $D$ and$\text{Td}$ is the Todd class. For specific operators, this simplifies beautifully:

Special Cases

de Rham ($d + d^*$): $\text{ind} = \chi(M) = \int_M e(TM)$

Gauss-Bonnet-Chern theorem

Signature ($d + d^*$ on $\Omega^\pm$): $\text{ind} = \sigma(M) = \int_M L(TM)$

Hirzebruch signature theorem

Dirac ($\not{D}$): $\text{ind} = \int_M \hat{A}(TM) \wedge \text{ch}(E)$

The Dirac index theorem (chiral anomaly)

Dolbeault ($\bar{\partial}$): $\text{ind} = \chi(M, \mathcal{O}(E)) = \int_M \text{ch}(E) \wedge \text{Td}(TM)$

Hirzebruch-Riemann-Roch theorem

Derivation 3: The Dirac Index and $\hat{A}$-Genus

For the Dirac operator $\not{D}$ on a spin manifold $M^{2n}$ coupled to a vector bundle $E$:

$\boxed{\text{ind}(\not{D}_E) = \int_M \hat{A}(TM) \wedge \text{ch}(E)}$

The $\hat{A}$-genus (A-hat genus) is:

$\hat{A}(TM) = \prod_i \frac{x_i/2}{\sinh(x_i/2)} = 1 - \frac{p_1}{24} + \frac{7p_1^2 - 4p_2}{5760} + \cdots$

4-Dimensional Formula

In 4 dimensions, for $E$ a rank-$r$ bundle with instanton number $k$:

$\text{ind}(\not{D}_E) = \int_M \left(-\frac{p_1(TM)}{24} + c_2(E)\right) = -\frac{\sigma(M)}{8} + k$

On $S^4$ ($\sigma = 0$): $\text{ind}(\not{D}_E) = k$. An SU(2) instanton with charge $k$produces exactly $k$ zero modes of the Dirac operator.

Integrality constraint: Since $\hat{A}(M)$ must be an integer for spin manifolds, this constrains the topology. For example, $p_1(M^4)$ must be divisible by 24 for a spin 4-manifold with $E$ trivial. This rules out $\mathbb{CP}^2$ as a spin manifold (since $p_1 = 3$, not divisible by 24).

Derivation 4: Anomalies in Quantum Field Theory

The chiral anomaly in 4D gauge theory with massless fermions:

$\partial_\mu j^\mu_5 = \frac{1}{16\pi^2}\text{Tr}(F_{\mu\nu}\tilde{F}^{\mu\nu})$

This is a direct consequence of the index theorem. In an instanton background:

$\Delta Q_5 = n_+ - n_- = \text{ind}(\not{D}) = \frac{1}{8\pi^2}\int \text{Tr}(F \wedge F) = 2k$

Anomaly Cancellation

A gauge theory is consistent only if the gauge anomaly cancels. For a chiral fermion in representation $R$, the anomaly is proportional to:

$\mathcal{A} \propto \text{Tr}_R(T^a\{T^b, T^c\})$

The Standard Model is anomaly-free precisely because the hypercharges of quarks and leptons satisfy $\sum Y = 0$, $\sum Y^3 = 0$, and $\text{Tr}(T_a^2 Y) = 0$ within each generation. This is a highly nontrivial constraint on the particle content.

Green-Schwarz mechanism: In 10D type I and heterotic string theories, gravitational and gauge anomalies do not cancel individually but combine via the factorization$I_{12} = X_4 \wedge X_8$. A Chern-Simons coupling $B \wedge X_8$ cancels the residual anomaly. This works only for gauge groups $SO(32)$ and $E_8 \times E_8$.

Derivation 5: Families Index and Determinant Line Bundles

When the Dirac operator depends on parameters (e.g., a family of metrics or gauge fields), the index theorem extends to the families index theorem. The index becomes a virtual vector bundle over the parameter space:

$\text{ind}(D_t) = [\ker D_t] - [\text{coker}\,D_t] \in K(B)$

and the families index theorem computes its Chern character:

$\text{ch}(\text{ind}(D)) = \int_{M/B} \hat{A}(TM/B) \wedge \text{ch}(E)$

Determinant Line Bundle

The determinant line bundle $\mathcal{L} = \det(\text{ind}(D))$over the space of gauge connections is the object whose sections define the fermion path integral. Its curvature encodes the anomaly:

$c_1(\mathcal{L}) = \text{anomaly polynomial}$

A gauge anomaly means $\mathcal{L}$ is non-trivial over the gauge orbit space—the fermion determinant is not a well-defined function but a section of a line bundle.

Witten's global anomaly (1982): Even when the local anomaly cancels, there can be a global (non-perturbative) anomaly if $\pi_4(G)$ is nontrivial. For $SU(2)$ with an odd number of Weyl fermion doublets, $\pi_4(SU(2)) = \mathbb{Z}_2$implies the partition function changes sign under a topologically nontrivial gauge transformation, rendering the theory inconsistent.

Interactive Simulation

This simulation computes the Dirac index on $S^2$ with monopole backgrounds, verifies the index theorem on 4-manifolds with instantons, plots zero mode density distributions, and demonstrates the $\hat{A}$-genus as a function of the first Pontryagin class.

Index Theorems: Atiyah-Singer, Dirac Index & Anomalies

Python

script.py220 lines

#!/usr/bin/env python3
"""
index_theorems.py - Atiyah-Singer Index Theorem, Dirac Index, and Anomalies
Computes indices for Dirac operators, verifies the index theorem on S^2 and S^4
"""
import numpy as np
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt
import base64
from io import BytesIO

# ============================================================
# 1. Euler characteristic as an index
# ============================================================
print("=" * 72)
print("  INDEX THEOREMS: THE DEEPEST BRIDGE BETWEEN GEOMETRY AND ANALYSIS")
print("=" * 72)
print()

print("  GAUSS-BONNET AS AN INDEX THEOREM")
print("  " + "-" * 50)
print()
print("  The Euler characteristic is the index of (d + d*):")
print("    ind(d + d*) = dim ker(d + d*)|_even - dim ker(d + d*)|_odd")
print("               = sum(-1)^k b_k = chi(M)")
print()

manifolds = [
    ("S^2", 2, [1, 0, 1]),
    ("T^2", 0, [1, 2, 1]),
    ("S^4", 2, [1, 0, 0, 0, 1]),
    ("CP^2", 3, [1, 0, 1, 0, 1]),
    ("K3", 24, [1, 0, 22, 0, 1]),
]

for name, chi, betti in manifolds:
    betti_str = ", ".join(str(b) for b in betti)
    check = sum((-1)**k * b for k, b in enumerate(betti))
    print("    {}: chi = {}, b = ({}), check: {}".format(name, chi, betti_str, check))
print()

# ============================================================
# 2. Dirac index on S^2 with monopole background
# ============================================================
print("  DIRAC INDEX ON S^2 WITH MONOPOLE")
print("  " + "-" * 50)
print()
print("  On S^2 with U(1) monopole of charge n:")
print("    ind(D) = n+ - n- = (n+1) for n >= 0")
print("                     = 0 for n = 0")
print("    (from Atiyah-Singer: ind = integral ch(L) A-hat(TM))")
print()

for n in range(0, 8):
    # On S^2: ind(D_n) = n (number of zero modes of positive chirality minus negative)
    # For monopole charge n >= 0: n+ = n+1, n- = 1 if n=0, else 0
    if n == 0:
        n_plus, n_minus = 1, 1
    else:
        n_plus, n_minus = n, 0
    index = n_plus - n_minus
    print("    n = {}: n+ = {}, n- = {}, index = {}".format(n, n_plus, n_minus, index))
print()

# ============================================================
# 3. Atiyah-Singer for Dirac on 4-manifolds
# ============================================================
print("  DIRAC INDEX ON 4-MANIFOLDS")
print("  " + "-" * 50)
print()
print("  ind(D) = integral A-hat(TM) ch(E)")
print("  For 4-manifold with SU(2) bundle of instanton number k:")
print("    ind(D) = integral (-p_1/24 + c_2(E))")
print("           = -sigma/8 + k")
print()

manifolds_4d = [
    ("S^4", 0, [1, 2, 3, 4]),
    ("CP^2", 1, [1, 2, 3, 4]),
    ("K3", -2, [1, 2, 3, 4]),
]

for name, sigma_over_8, k_vals in manifolds_4d:
    print("    {} (sigma/8 = {}):".format(name, sigma_over_8))
    for k in k_vals:
        index = -sigma_over_8 + k
        print("      k = {}: ind(D) = {}".format(k, index))
    print()

# ============================================================
# 4. Anomaly computation
# ============================================================
print("  CHIRAL ANOMALY FROM INDEX THEOREM")
print("  " + "-" * 50)
print()
print("  In 4D QCD with SU(N) gauge group:")
print("    n+ - n- = integral (1/8pi^2) Tr(F ^ F) = k (instanton number)")
print()
print("  This means: in an instanton background with charge k,")
print("  the Dirac operator has k zero modes of one chirality.")
print("  This causes 't Hooft vertices violating B+L conservation.")
print()

# Anomaly coefficients for Standard Model
print("  Standard Model anomaly cancellation:")
print("    SU(3)^3: sum of cubic Casimirs = 0")
print("    SU(2)^2 U(1): sum = 0 (requires 3 generations)")
print("    U(1)^3: sum of Y^3 = 0")
print("    gravity-gauge: sum of Y = 0")
print()

# Per generation: quarks + leptons
# Q_L: (3, 2, 1/6), u_R: (3, 1, 2/3), d_R: (3, 1, -1/3)
# L_L: (1, 2, -1/2), e_R: (1, 1, -1)
# Sum Y: 6*(1/6) + 3*(2/3) + 3*(-1/3) + 2*(-1/2) + 1*(-1) = 1 + 2 - 1 - 1 - 1 = 0
print("  Per generation, sum of Y:")
print("    Q_L: 6 * (1/6) = 1")
print("    u_R: 3 * (2/3) = 2")
print("    d_R: 3 * (-1/3) = -1")
print("    L_L: 2 * (-1/2) = -1")
print("    e_R: 1 * (-1) = -1")
print("    Total: 1 + 2 - 1 - 1 - 1 = 0 [anomaly free]")
print()

# ============================================================
# 5. A-hat genus computation
# ============================================================
print("  A-HAT GENUS")
print("  " + "-" * 50)
print()
print("  A-hat(TM) = 1 - p_1/24 + (7p_1^2 - 4p_2)/5760 + ...")
print()

for name, p1, p2, dim in [("S^4", 0, 0, 4), ("CP^2", 3, 0, 4),
                           ("K3", -48, 0, 4), ("S^2 x S^2", 0, 0, 4)]:
    if dim == 4:
        a_hat = -p1 / 24.0
        print("    {}: p_1 = {}, A-hat = -p_1/24 = {:.4f}".format(name, p1, a_hat))
print()

# ============================================================
# Visualization
# ============================================================
fig, axes = plt.subplots(2, 2, figsize=(14, 11))
fig.patch.set_facecolor('#0a0a0a')

for ax in axes.flat:
    ax.set_facecolor('#111111')
    for spine in ax.spines.values():
        spine.set_color('#2dd4bf')
    ax.tick_params(colors='#99f6e4')
    ax.grid(True, color='#134e4a', alpha=0.4)

# Plot 1: Dirac index vs monopole charge on S^2
ax1 = axes[0, 0]
n_vals = np.arange(0, 10)
indices = n_vals.copy()
ax1.plot(n_vals, indices, 'o-', color='#2dd4bf', linewidth=2, markersize=8)
ax1.set_xlabel('Monopole charge n', color='#5eead4')
ax1.set_ylabel('Dirac index', color='#5eead4')
ax1.set_title('Dirac Index on S^2 with Monopole', color='#99f6e4', fontsize=12, fontweight='bold')

# Plot 2: Index vs instanton number on S^4
ax2 = axes[0, 1]
k_vals = np.arange(0, 8)
for sigma_label, sigma_val, color in [('S^4 (sig=0)', 0, '#2dd4bf'),
                                       ('CP^2 (sig=1)', 1, '#f472b6'),
                                       ('K3 (sig=-16)', -2, '#f59e0b')]:
    idx = -sigma_val + k_vals
    ax2.plot(k_vals, idx, 'o-', linewidth=2, markersize=6, label=sigma_label, color=color)
ax2.set_xlabel('Instanton number k', color='#5eead4')
ax2.set_ylabel('Dirac index', color='#5eead4')
ax2.set_title('Dirac Index on 4-Manifolds', color='#99f6e4', fontsize=12, fontweight='bold')
ax2.legend(facecolor='#111111', edgecolor='#134e4a', labelcolor='#99f6e4', fontsize=8)

# Plot 3: Zero mode density for monopole on S^2
ax3 = axes[1, 0]
theta_zero = np.linspace(0.01, np.pi - 0.01, 300)
for n in [1, 2, 3, 4]:
    # Zero mode density ~ sin^(n-1)(theta/2) cos^(n-1)(theta/2) sin(theta)
    density = np.sin(theta_zero)**(n) * (np.sin(theta_zero/2))**(2*n - 2)
    density = density / (np.trapz(density, theta_zero) + 1e-10)
    ax3.plot(theta_zero, density, linewidth=2.5, label='n = {}'.format(n))
ax3.set_xlabel('theta', color='#5eead4')
ax3.set_ylabel('Zero mode density', color='#5eead4')
ax3.set_title('Dirac Zero Modes on S^2', color='#99f6e4', fontsize=12, fontweight='bold')
ax3.legend(facecolor='#111111', edgecolor='#134e4a', labelcolor='#99f6e4', fontsize=8)

# Plot 4: Anomaly polynomial structure
ax4 = axes[1, 1]
# Plot the A-hat characteristic class contribution vs p_1
p1_range = np.linspace(-60, 60, 200)
a_hat_vals = -p1_range / 24.0
ax4.plot(p1_range, a_hat_vals, color='#34d399', linewidth=2.5, label='A-hat_4 = -p_1/24')
ax4.axhline(y=0, color='#6b7280', linewidth=1, linestyle='--')
ax4.axvline(x=0, color='#6b7280', linewidth=1, linestyle='--')

# Mark known manifolds
known = [("S^4", 0), ("CP^2", 3), ("K3", -48)]
for name, p1v in known:
    av = -p1v / 24.0
    ax4.plot(p1v, av, 'o', markersize=10, color='#f59e0b')
    ax4.annotate(name, xy=(p1v, av), xytext=(p1v + 3, av + 0.2),
                 color='#f59e0b', fontsize=10)
ax4.set_xlabel('p_1(TM)', color='#5eead4')
ax4.set_ylabel('A-hat genus', color='#5eead4')
ax4.set_title('A-hat Genus vs First Pontryagin Class', color='#99f6e4', fontsize=12, fontweight='bold')
ax4.legend(facecolor='#111111', edgecolor='#134e4a', labelcolor='#99f6e4')

plt.tight_layout(pad=2.0)
buf = BytesIO()
plt.savefig(buf, format='png', dpi=150, bbox_inches='tight', facecolor='#0a0a0a')
buf.seek(0)
img_b64 = base64.b64encode(buf.read()).decode()
print('<img src="data:image/png;base64,' + img_b64 + '" alt="Index Theorems Visualization" />')
buf.close()
print()
print("Visualization complete.")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Summary

Atiyah-Singer Index Theorem

Equates the analytical index (dim ker - dim coker) of an elliptic operator to a topological integral of characteristic classes. Unifies Gauss-Bonnet, signature, and Riemann-Roch theorems.

Dirac Index

$\text{ind}(\not{D}_E) = \int \hat{A}(TM) \wedge \text{ch}(E)$. Counts the net number of chiral zero modes. In 4D: $\text{ind} = -\sigma/8 + k$ (instanton number minus signature contribution).

Chiral Anomaly

The index theorem explains why the axial current is not conserved in an instanton background:$\Delta Q_5 = 2k$. Anomaly cancellation constrains the particle content of gauge theories.

Families Index

Extends the index to parametrized families of operators. The determinant line bundle encodes anomalies, and global anomalies arise from nontrivial homotopy groups of the gauge group.

← Homology & Cohomology Course Overview →