Chapter 8: Anomalies

Anomalies are quantum violations of classical symmetries. The most important is the chiral (Adler-Bell-Jackiw) anomaly: the axial-vector current, conserved at the classical level, acquires a divergence through quantum effects. This is not a failure of the theory but a profound feature — anomalies constrain the particle content of consistent gauge theories, predict the $\pi^0 \to \gamma\gamma$ decay rate, and play essential roles in topology, index theorems, and the structure of the Standard Model.

The Chiral Anomaly

Consider a massless Dirac fermion coupled to an electromagnetic field. Classically, the theory has two conserved currents: the vector current $j^\mu = \bar{\psi}\gamma^\mu\psi$ (corresponding to charge conservation) and the axial-vector current $j^{5\mu} = \bar{\psi}\gamma^\mu\gamma^5\psi$ (corresponding to chiral symmetry). At the classical level:

$\partial_\mu j^\mu = 0 \quad \text{(vector current conservation)}$

$\partial_\mu j^{5\mu} = 0 \quad \text{(classical axial current conservation, for } m=0\text{)}$

However, quantum corrections — specifically the triangle diagram with one axial and two vector vertices — violate the axial current conservation. The Adler-Bell-Jackiw (ABJ) result is:

$\partial_\mu j^{5\mu} = \frac{e^2}{16\pi^2}F^{\mu\nu}\tilde{F}_{\mu\nu}$

where $\tilde{F}_{\mu\nu} = \frac{1}{2}\varepsilon_{\mu\nu\rho\sigma}F^{\rho\sigma}$ is the dual field strength tensor. This is a remarkable result: the right-hand side is a total derivative ($F\tilde{F} = \partial_\mu K^\mu$ for a certain current $K^\mu$), but cannot be removed by any regularization that preserves gauge (vector current) invariance.

The anomaly is exact — it receives no corrections beyond one loop. This is the Adler-Bardeen theorem. The coefficient $e^2/(16\pi^2)$ is fixed by the triangle diagram and is not modified by higher-order perturbative or non-perturbative effects.

The Triangle Diagram

The anomaly arises from the AVV (axial-vector-vector) triangle diagram: a fermion loop with one axial-vector vertex ($\gamma^\mu\gamma^5$) and two vector vertices ($\gamma^\nu, \gamma^\rho$). The amplitude is:

$T^{\mu\nu\rho}(k_1, k_2) = (-1)\int \frac{d^4\ell}{(2\pi)^4} \text{Tr}\left[\gamma^\mu\gamma^5 S(\ell)\gamma^\nu S(\ell - k_1)\gamma^\rho S(\ell - k_1 - k_2)\right] + (k_1 \leftrightarrow k_2, \nu \leftrightarrow \rho)$

The integral is linearly divergent and requires regularization. The key issue is that no regularization can simultaneously preserve both the vector Ward identity ($k_{1\nu}T^{\mu\nu\rho} = 0$) and the axial Ward identity ($(k_1 + k_2)_\mu T^{\mu\nu\rho} = 0$).

The calculation proceeds by introducing a Feynman parameter $x$, shifting the loop momentum, and evaluating the resulting integrals. The surface term that arises from the shift of a linearly divergent integral is the origin of the anomaly:

$\int \frac{d^4\ell}{(2\pi)^4}\frac{\partial}{\partial\ell^\mu}f(\ell) \neq 0 \quad \text{(for linearly divergent integrals)}$

This is the mathematical origin of the anomaly: the naive step of shifting integration variables in a divergent integral is not justified. The physical choice is to preserve vector current conservation (gauge invariance) at the expense of axial current conservation.

Anomaly Coefficient

For a fermion with charge $Q$ in a representation with color multiplicity $N_c$, the anomaly coefficient is:

$\mathcal{A} = \sum_f N_c^f Q_f^2$

The anomaly is a one-loop exact result. For multiple fermion species, the anomaly coefficients simply add.

Fujikawa Method: Path Integral Derivation

Fujikawa provided an elegant non-perturbative derivation of the chiral anomaly using the path integral. The key insight is that while the classical action is invariant under chiral rotations, the path integral measure is not.

The Argument

Consider an infinitesimal chiral transformation:

$\psi \to (1 + i\alpha\gamma^5)\psi, \quad \bar{\psi} \to \bar{\psi}(1 + i\alpha\gamma^5)$

The classical action $S[\psi, \bar{\psi}, A]$ is invariant (for $m = 0$), but the path integral measure transforms as:

$\mathcal{D}\psi\mathcal{D}\bar{\psi} \to \mathcal{D}\psi\mathcal{D}\bar{\psi} \exp\left(-2i\alpha\int d^4x\, \mathcal{A}(x)\right)$

The Jacobian factor $\mathcal{A}(x)$ is computed by expanding $\psi$ in eigenmodes of the Dirac operator $i\not{D}$:

$\mathcal{A}(x) = \sum_n \phi_n^\dagger(x)\gamma^5\phi_n(x) = \lim_{M\to\infty}\sum_n \phi_n^\dagger(x)\gamma^5 e^{-\lambda_n^2/M^2}\phi_n(x)$

where $\lambda_n$ are the eigenvalues of $i\not{D}$ and we have introduced a Gaussian regulator. Evaluating this regulated sum using heat kernel techniques gives:

$\mathcal{A}(x) = \frac{e^2}{32\pi^2}F^{\mu\nu}\tilde{F}_{\mu\nu}$

This reproduces the ABJ anomaly. The Fujikawa method makes clear that the anomaly is fundamentally topological in nature — it is connected to the index of the Dirac operator via the Atiyah-Singer index theorem:

$\text{index}(i\not{D}) = n_+ - n_- = \frac{e^2}{32\pi^2}\int d^4x\, F^{\mu\nu}\tilde{F}_{\mu\nu}$

where $n_+$ and $n_-$ are the numbers of zero modes with positive and negative chirality. This topological connection explains why the anomaly is exact and receives no higher-order corrections.

Anomaly Cancellation in the Standard Model

For a gauge theory to be consistent at the quantum level, all gauge anomalies must cancel. An anomalous gauge symmetry would destroy the Ward identities that guarantee unitarity and renormalizability. The anomaly cancellation conditions for the Standard Model gauge group $SU(3)_C \times SU(2)_L \times U(1)_Y$ are:

$[SU(3)]^2 \times U(1)_Y$

$\sum_\text{quarks} Y_q = 2 \times \frac{1}{6} + \frac{2}{3} + \left(-\frac{1}{3}\right) = 0$

The factor of 2 counts the $SU(2)$ doublet components. Cancellation requires that the sum of hypercharges over all colored fermions vanishes.

$[SU(2)]^2 \times U(1)_Y$

$\sum_\text{doublets} N_c Y = 3 \times \frac{1}{6} + 1 \times \left(-\frac{1}{2}\right) = 0$

The sum runs over all left-handed $SU(2)$ doublets, weighted by color multiplicity.

$[U(1)_Y]^3$

$\sum_\text{all} N_c \cdot n_{SU(2)} \cdot Y^3 = 3 \cdot 2 \cdot \left(\frac{1}{6}\right)^3 + 3 \cdot \left(\frac{2}{3}\right)^3 + 3 \cdot \left(-\frac{1}{3}\right)^3 + 2 \cdot \left(-\frac{1}{2}\right)^3 + (-1)^3 = 0$

This cubic anomaly condition places the tightest constraint on the hypercharge assignments.

$[\text{grav}]^2 \times U(1)_Y$

$\sum_\text{all} N_c \cdot n_{SU(2)} \cdot Y = 3 \cdot 2 \cdot \frac{1}{6} + 3 \cdot \frac{2}{3} + 3 \cdot \left(-\frac{1}{3}\right) + 2 \cdot \left(-\frac{1}{2}\right) + (-1) = 0$

The mixed gravitational-gauge anomaly must also cancel for consistency in curved spacetime.

All four conditions are satisfied in the Standard Model, but only because of the precise relationship between quark and lepton quantum numbers within each generation. This remarkable cancellation strongly suggests that quarks and leptons are related at a deeper level, pointing toward grand unification.

$\pi^0 \to \gamma\gamma$ Decay

The neutral pion decay $\pi^0 \to \gamma\gamma$ is the classic physical application of the chiral anomaly. The $\pi^0$ is the Goldstone boson associated with the spontaneously broken axial $SU(2)$ symmetry, and its coupling to two photons is determined entirely by the anomaly.

The anomaly-induced effective coupling is:

$\mathcal{L}_{\pi^0\gamma\gamma} = \frac{\alpha}{8\pi f_\pi}N_c(Q_u^2 - Q_d^2)\pi^0 F^{\mu\nu}\tilde{F}_{\mu\nu}$

The decay rate is:

$\Gamma(\pi^0 \to \gamma\gamma) = \frac{\alpha^2 m_{\pi}^3}{64\pi^3 f_\pi^2}\left[N_c\left(Q_u^2 - Q_d^2\right)\right]^2$

With $N_c = 3$, $Q_u = 2/3$, $Q_d = -1/3$, the anomaly factor is $N_c(Q_u^2 - Q_d^2) = 3(4/9 - 1/9) = 1$. This gives a predicted lifetime of $\tau \approx 8.4 \times 10^{-17}$ s, in excellent agreement with the experimental value of $8.5 \times 10^{-17}$ s.

This result has deep significance:

• Without the anomaly ($\partial_\mu j^{5\mu} = 0$), the decay $\pi^0 \to \gamma\gamma$ would be forbidden, contradicting experiment.

• The rate is proportional to $N_c^2$, providing direct experimental evidence that quarks come in three colors.

• The prediction requires no adjustable parameters — it follows entirely from symmetry and the anomaly coefficient.

Global vs Gauge Anomalies

It is crucial to distinguish between anomalies in global symmetries and anomalies in gauge symmetries, as they have very different physical consequences.

Anomalous Global Symmetries

When a global symmetry has an anomaly, the symmetry is simply broken at the quantum level. The classical conservation law $\partial_\mu j^\mu = 0$ is replaced by an anomalous Ward identity. This is a physical effect with observable consequences. The chiral anomaly $\partial_\mu j^{5\mu} \neq 0$ is of this type: the axial$U(1)_A$ symmetry of massless QCD is broken by the anomaly, explaining why the$\eta'$ meson is much heavier than the pion (the $U(1)_A$ problem).

Anomalous Gauge Symmetries

When a gauge symmetry has an anomaly, the theory is inconsistent. The Ward identities that guarantee unitarity and renormalizability are violated, and the theory produces nonsensical predictions (negative probabilities, non-unitary S-matrix). Therefore, gauge anomaly cancellation is a strict consistency requirement for any viable gauge theory. This is why the Standard Model must have exactly the right fermion content to cancel all gauge anomalies.

't Hooft Anomaly Matching

't Hooft's anomaly matching condition is a powerful non-perturbative constraint: the anomaly coefficients of a global symmetry, computed in terms of the fundamental (UV) degrees of freedom, must match those computed in terms of the composite (IR) degrees of freedom:

$\mathcal{A}_\text{UV}[G_\text{global}] = \mathcal{A}_\text{IR}[G_\text{global}]$

This condition constrains the possible infrared dynamics of strongly-coupled theories. For example, it played a key role in understanding the chiral symmetry breaking pattern in QCD and constraining the spectrum of composite (bound state) fermions.

Topological Aspects of Anomalies

The deep connection between anomalies and topology is one of the most beautiful aspects of modern quantum field theory. The integrated anomaly equation relates the change in fermion number to the topological charge of the gauge field:

$\Delta Q_5 = \int d^4x\, \partial_\mu j^{5\mu} = \frac{e^2}{16\pi^2}\int d^4x\, F\tilde{F} = 2\nu$

where $\nu \in \mathbb{Z}$ is the instanton number (Pontryagin index) of the gauge field configuration. This means that in the background of an instanton ($\nu = 1$), two units of chirality are created — one fermion of each chirality is transmuted into the other.

Instantons and Tunneling

Instantons are topologically non-trivial solutions of the Euclidean field equations. In QCD, they mediate tunneling between vacua with different Chern-Simons numbers, and the anomaly determines the fermion zero modes in the instanton background. The Atiyah-Singer index theorem states:

$n_+ - n_- = \nu$

where $n_\pm$ count the zero modes of positive/negative chirality. Each quark flavor contributes one zero mode per instanton, leading to the 't Hooft effective vertex that violates the anomalous $U(1)_A$ symmetry by $2N_f$ units. In the electroweak theory, the analogous process (the sphaleron) violates baryon plus lepton number $B + L$ while conserving $B - L$, providing a mechanism for baryogenesis in the early universe.

The $\theta$ Vacuum and Strong CP Problem

The existence of instantons implies that the QCD vacuum is characterized by a parameter $\theta$:

$\mathcal{L}_\theta = \frac{\theta g^2}{32\pi^2}G^a_{\mu\nu}\tilde{G}^{a\mu\nu}$

This term violates CP symmetry unless $\theta = 0$ or $\pi$. The experimental bound on the neutron electric dipole moment requires $|\theta| < 10^{-10}$. Why$\theta$ is so small is the strong CP problem. The most elegant solution is the Peccei-Quinn mechanism, which introduces a new symmetry whose spontaneous breaking produces a light pseudoscalar — the axion — that dynamically relaxes $\theta$ to zero. Axions are also a leading candidate for dark matter.

Non-Perturbative Anomalies

Beyond the perturbative ABJ anomaly, there exist non-perturbative (global) anomalies discovered by Witten. The Witten $SU(2)$ anomaly states that an $SU(2)$ gauge theory with an odd number of Weyl fermion doublets is inconsistent:

$Z \to -Z \quad \text{under a topologically non-trivial gauge transformation}$

This anomaly cannot be detected by perturbative Feynman diagram calculations — it arises from the fact that $\pi_4(SU(2)) = \mathbb{Z}_2$, meaning there exist gauge transformations that are not continuously connected to the identity.

In the Standard Model, the electroweak $SU(2)_L$ has an even number of doublets per generation (the quark doublet counts as three due to color), so the Witten anomaly cancels. This provides yet another non-trivial consistency check on the Standard Model fermion content.

The study of anomalies has become a central theme in modern theoretical physics, connecting quantum field theory to topology, differential geometry, and condensed matter physics. Anomaly inflow, anomaly polynomials, and their classification by cobordism theory continue to provide deep insights into the structure of quantum field theories and their consistent formulation.

Computational Analysis: Anomaly Triangle Diagram

We compute the anomaly coefficients for the Standard Model fermions, verify anomaly cancellation, calculate the $\pi^0 \to \gamma\gamma$ decay rate as a function of the number of colors, and explore the spectral flow of the Dirac operator in a background field.

Anomaly Triangle Diagram Contributions

Python

script.py357 lines

#!/usr/bin/env python3
"""
anomalies.py - Chiral Anomaly and Triangle Diagram Contributions
Computes anomaly coefficients, pi0 -> gamma gamma, and visualizes the triangle diagram
"""
import numpy as np
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt
import base64
from io import BytesIO

print("=" * 72)
print("  ANOMALIES: TRIANGLE DIAGRAMS & CHIRAL ANOMALY")
print("=" * 72)
print()

# ============================================================
# 1. Triangle Diagram: AVV Amplitude
# The anomalous amplitude for axial current coupling to two photons:
# T^{mu nu lambda}(k1, k2) = Tr[gamma^lambda gamma^5 S(q) gamma^mu S(q-k1) gamma^nu S(q-k1-k2)]
#                           + crossed diagram
# After integration, this gives the ABJ anomaly
# ============================================================
print("1. TRIANGLE DIAGRAM CALCULATION (AVV)")
print("-" * 50)

# The triangle diagram integral (simplified numerical model)
# We compute the tensor structure numerically for demonstration
# The key result is that the divergent part cannot be regularized
# to preserve both vector and axial-vector current conservation

# Pauli-Villars style computation:
# Anomaly coefficient = sum over fermions of Q_f^2 * T_f * N_c

# Standard Model fermion content
fermions = {
    'e':    {'Q': -1.0, 'T3': -0.5, 'Y': -1.0, 'Nc': 1},
    'nu_e': {'Q':  0.0, 'T3':  0.5, 'Y': -1.0, 'Nc': 1},
    'u':    {'Q':  2/3, 'T3':  0.5, 'Y':  1/3, 'Nc': 3},
    'd':    {'Q': -1/3, 'T3': -0.5, 'Y':  1/3, 'Nc': 3},
}

print("  Fermion anomaly coefficients (per generation):")
print(f"  {'Fermion':>8s}  {'Q':>6s}  {'T3':>6s}  {'Y':>6s}  {'Nc':>4s}  {'Nc*Q^2':>8s}  {'Nc*Y^3':>8s}")
total_Q2 = 0
total_Y3 = 0
for name, props in fermions.items():
    nc_q2 = props['Nc'] * props['Q']**2
    nc_y3 = props['Nc'] * props['Y']**3
    total_Q2 += nc_q2
    total_Y3 += nc_y3
    print(f"  {name:>8s}  {props['Q']:>6.2f}  {props['T3']:>6.2f}  {props['Y']:>6.2f}  {props['Nc']:>4d}  {nc_q2:>8.4f}  {nc_y3:>8.4f}")

print(f"  {'Total':>8s}  {' ':>6s}  {' ':>6s}  {' ':>6s}  {' ':>4s}  {total_Q2:>8.4f}  {total_Y3:>8.4f}")
print()

# ============================================================
# 2. Anomaly Cancellation in the Standard Model
# For a consistent quantum theory, gauge anomalies must cancel:
# Tr[T^a {T^b, T^c}] = 0 for gauge generators
# ============================================================
print("2. ANOMALY CANCELLATION IN THE STANDARD MODEL")
print("-" * 50)

# The anomaly cancellation conditions per generation:
# 1. [SU(3)]^2 U(1)_Y: sum Nc * Y_q for quarks = 0
# 2. [SU(2)]^2 U(1)_Y: sum Y_i for doublets = 0
# 3. [U(1)_Y]^3: sum Nc * Y^3 = 0
# 4. [grav]^2 U(1)_Y: sum Nc * Y = 0
# 5. [SU(2)]^2 [SU(2)]: automatically zero (SU(2) is safe)
# 6. [SU(3)]^3: automatically zero (real rep contribution)

# Full SM fermion content (one generation, left-handed)
sm_fermions = {
    'Q_L (u,d)':   {'SU3': 3, 'SU2': 2, 'Y':  1/6},   # quark doublet
    'u_R':         {'SU3': 3, 'SU2': 1, 'Y':  2/3},   # right-handed up
    'd_R':         {'SU3': 3, 'SU2': 1, 'Y': -1/3},   # right-handed down
    'L_L (nu,e)':  {'SU3': 1, 'SU2': 2, 'Y': -1/2},   # lepton doublet
    'e_R':         {'SU3': 1, 'SU2': 1, 'Y': -1.0},   # right-handed electron
}

print("  SM Weyl fermions (one generation, left-handed):")
print(f"  {'Field':>15s}  {'SU(3)':>6s}  {'SU(2)':>6s}  {'Y':>8s}")
for name, props in sm_fermions.items():
    print(f"  {name:>15s}  {props['SU3']:>6d}  {props['SU2']:>6d}  {props['Y']:>8.4f}")
print()

# Check anomaly cancellation conditions
# Condition 1: [SU(3)]^2 U(1)_Y
A_33Y = 0
for name, props in sm_fermions.items():
    if props['SU3'] == 3:
        A_33Y += props['SU2'] * props['Y']  # T(fund) = 1/2 for SU(3) fund
print(f"  [SU(3)]^2 x U(1)_Y anomaly: {A_33Y:.6f}  {'CANCELLED!' if abs(A_33Y) < 1e-10 else 'NOT CANCELLED'}")

# Condition 2: [SU(2)]^2 U(1)_Y
A_22Y = 0
for name, props in sm_fermions.items():
    if props['SU2'] == 2:
        A_22Y += props['SU3'] * props['Y']  # T(fund) = 1/2 for SU(2) fund
print(f"  [SU(2)]^2 x U(1)_Y anomaly: {A_22Y:.6f}  {'CANCELLED!' if abs(A_22Y) < 1e-10 else 'NOT CANCELLED'}")

# Condition 3: [U(1)_Y]^3
A_YYY = 0
for name, props in sm_fermions.items():
    A_YYY += props['SU3'] * props['SU2'] * props['Y']**3
print(f"  [U(1)_Y]^3 anomaly:         {A_YYY:.6f}  {'CANCELLED!' if abs(A_YYY) < 1e-10 else 'NOT CANCELLED'}")

# Condition 4: [grav]^2 U(1)_Y (gravitational anomaly)
A_gravY = 0
for name, props in sm_fermions.items():
    A_gravY += props['SU3'] * props['SU2'] * props['Y']
print(f"  [grav]^2 x U(1)_Y anomaly:  {A_gravY:.6f}  {'CANCELLED!' if abs(A_gravY) < 1e-10 else 'NOT CANCELLED'}")

print()

# ============================================================
# 3. Pi0 -> gamma gamma Decay Rate
# Gamma = (alpha^2 * m_pi^3) / (64 * pi^3 * f_pi^2) * |A|^2
# A = Nc * (Q_u^2 - Q_d^2) = 3 * (4/9 - 1/9) = 1
# ============================================================
print("3. PI0 -> GAMMA GAMMA DECAY RATE")
print("-" * 50)

alpha_em = 1.0 / 137.036
m_pi0 = 0.135  # GeV
f_pi = 0.092   # pion decay constant in GeV
Nc = 3

# Anomaly amplitude factor
A_anomaly = Nc * ((2.0/3)**2 - (1.0/3)**2)  # Nc * (Q_u^2 - Q_d^2)
print(f"  Anomaly factor A = Nc*(Q_u^2 - Q_d^2) = {A_anomaly:.4f}")

# Decay rate from anomaly
Gamma_pi0 = (alpha_em**2 * m_pi0**3) / (64 * np.pi**3 * f_pi**2) * A_anomaly**2
tau_pi0 = 1.0 / Gamma_pi0  # lifetime in natural units

# Convert to seconds (hbar = 6.582e-25 GeV*s)
hbar = 6.582e-25  # GeV*s
tau_pi0_sec = hbar / Gamma_pi0

print(f"  Decay rate: Gamma = {Gamma_pi0:.6e} GeV")
print(f"  Lifetime: tau = {tau_pi0_sec:.4e} seconds")
print(f"  Experimental value: tau = 8.5 x 10^-17 seconds")
print(f"  Ratio (theory/experiment): {tau_pi0_sec / 8.5e-17:.3f}")
print()
print("  The anomaly prediction agrees with experiment to ~3%!")
print("  This was a crucial early test that Nc = 3 (three colors).")
print()

# ============================================================
# 4. Fujikawa Method: Path Integral Measure Non-Invariance
# Under chiral rotation psi -> exp(i*alpha*gamma5) psi:
# D[psi]D[psi_bar] -> D[psi]D[psi_bar] * exp(-2i*alpha*A[A])
# where A[A] = (e^2/(16*pi^2)) * integral F_munu * F_tilde^munu
# ============================================================
print("4. FUJIKAWA METHOD: ANOMALY FROM PATH INTEGRAL MEASURE")
print("-" * 50)

# The Jacobian of the chiral transformation involves the index of the Dirac operator
# Index(D) = n_+ - n_- = (1/(16*pi^2)) * integral F * F_tilde
# (Atiyah-Singer index theorem)

# Demonstrate with a simple model: 2D fermion in background field
# The anomaly in 2D: d_mu j^5_mu = (e/(2*pi)) * epsilon^{mu nu} F_{mu nu}
# = (e/pi) * E (the electric field)

print("  Chiral transformation: psi -> exp(i*alpha*gamma_5) psi")
print("  The path integral measure transforms as:")
print("    D[psi]D[psi_bar] -> D[psi]D[psi_bar] exp(-2i*alpha * A[A])")
print()
print("  where the anomaly functional is:")
print("    A[A] = (e^2 / (16*pi^2)) * integral d^4x F_{mu nu} F_tilde^{mu nu}")
print()
print("  This gives the anomalous divergence:")
print("    partial_mu j^5_mu = (e^2 / (8*pi^2)) * F_{mu nu} F_tilde^{mu nu}")
print()

# Numerical demonstration: eigenvalue flow in a simple model
# Consider a 1D Dirac operator with a background gauge field
N_sites = 200
L = 10.0
dx = L / N_sites
x = np.linspace(-L/2, L/2, N_sites)

# Background gauge field (instanton-like profile)
A_bg_strengths = [0.0, 0.5, 1.0, 2.0, 3.0]

print("  Spectral flow of Dirac operator eigenvalues:")
print(f"  {'A_strength':>12s}  {'n_zero_modes':>14s}  {'spectral_asym':>15s}")
for A_str in A_bg_strengths:
    # Simple 1D Dirac-like operator: H = -i*d/dx + A(x)
    A_field = A_str * np.exp(-x**2 / 2)  # Gaussian profile

# Construct discrete Hamiltonian
    H = np.zeros((N_sites, N_sites))
    for i in range(N_sites):
        if i < N_sites - 1:
            H[i, i+1] = -1j / (2 * dx)
            H[i+1, i] = 1j / (2 * dx)
        H[i, i] = A_field[i]

eigenvalues = np.linalg.eigvalsh(H.real)
    n_near_zero = np.sum(np.abs(eigenvalues) < 0.1)
    spectral_asym = np.sum(np.sign(eigenvalues))
    print(f"  {A_str:>12.1f}  {n_near_zero:>14d}  {spectral_asym:>15.1f}")

print()

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

# --- Plot 1: Triangle diagram schematic (momentum flow) ---
ax1 = axes[0, 0]
ax1.set_facecolor('#1a1025')
for spine in ax1.spines.values():
    spine.set_color('#4c1d95')
ax1.tick_params(colors='#c4b5fd')

# Draw triangle diagram
triangle_x = [0, 2, 1, 0]
triangle_y = [0, 0, np.sqrt(3), 0]
ax1.plot(triangle_x, triangle_y, color='#a78bfa', linewidth=3)

# Vertices
ax1.plot([0, 2, 1], [0, 0, np.sqrt(3)], 'o', color='#22d3ee', markersize=12, zorder=5)

# Labels
ax1.annotate('gamma (k1)', xy=(-0.3, -0.3), color='#f472b6', fontsize=12, fontweight='bold')
ax1.annotate('gamma (k2)', xy=(1.8, -0.3), color='#f472b6', fontsize=12, fontweight='bold')
ax1.annotate('j^5_mu (q)', xy=(0.5, np.sqrt(3)+0.15), color='#fbbf24', fontsize=12, fontweight='bold')

# Momentum arrows on fermion lines
ax1.annotate('', xy=(1.0, 0), xytext=(0.0, 0),
    arrowprops=dict(arrowstyle='->', color='#34d399', lw=2))
ax1.annotate('', xy=(1.8, np.sqrt(3)*0.6), xytext=(2, 0),
    arrowprops=dict(arrowstyle='->', color='#34d399', lw=2))
ax1.annotate('', xy=(0.2, np.sqrt(3)*0.6), xytext=(1, np.sqrt(3)),
    arrowprops=dict(arrowstyle='->', color='#34d399', lw=2))

# Wavy lines for photons
for t in np.linspace(0, 1, 20):
    x_wave = -0.5 + t * 0.5
    y_wave = -0.5 + 0.1 * np.sin(t * 6 * np.pi)
    ax1.plot(x_wave, y_wave, '.', color='#f472b6', markersize=2)

for t in np.linspace(0, 1, 20):
    x_wave = 2 + t * 0.5
    y_wave = -0.5 + 0.1 * np.sin(t * 6 * np.pi)
    ax1.plot(x_wave, y_wave, '.', color='#f472b6', markersize=2)

ax1.set_xlim(-1, 3)
ax1.set_ylim(-1, 2.5)
ax1.set_title('AVV Triangle Diagram', color='#e9d5ff', fontsize=13, fontweight='bold')
ax1.set_aspect('equal')
ax1.axis('off')

# --- Plot 2: Anomaly coefficients per SM generation ---
ax2 = axes[0, 1]
ax2.set_facecolor('#1a1025')
for spine in ax2.spines.values():
    spine.set_color('#4c1d95')
ax2.tick_params(colors='#c4b5fd')
ax2.grid(True, color='#4c1d95', alpha=0.3)

anomaly_types = ['[SU(3)]^2 Y', '[SU(2)]^2 Y', '[Y]^3', '[grav]^2 Y']
lepton_contrib = [0, -0.5 + 0, -(-0.5)**3 - (-1)**3, -0.5*2 + (-1)]
quark_contrib = [2*(1/6), 3*2*(1/6), 3*2*(1/6)**3 + 3*(2/3)**3 + 3*(-1/3)**3, 3*2*(1/6) + 3*(2/3) + 3*(-1/3)]
# Recalculate properly
lepton_vals = [0, -1/2, 2*(-1/2)**3 + (-1)**3, 2*(-1/2) + (-1)]
quark_vals = [2*(1/6), 3*(1/6), 3*2*(1/6)**3 + 3*(2/3)**3 + 3*(-1/3)**3, 3*2*(1/6) + 3*(2/3) + 3*(-1/3)]
total_vals = [l+q for l, q in zip(lepton_vals, quark_vals)]

x_pos = np.arange(len(anomaly_types))
width = 0.25

bars1 = ax2.bar(x_pos - width, lepton_vals, width, color='#22d3ee', label='Leptons', alpha=0.8)
bars2 = ax2.bar(x_pos, quark_vals, width, color='#a78bfa', label='Quarks', alpha=0.8)
bars3 = ax2.bar(x_pos + width, total_vals, width, color='#34d399', label='Total', alpha=0.8)

ax2.set_xticks(x_pos)
ax2.set_xticklabels(anomaly_types, color='#c4b5fd', fontsize=9)
ax2.set_ylabel('Anomaly coefficient', color='#a78bfa', fontsize=11)
ax2.set_title('SM Anomaly Cancellation', color='#e9d5ff', fontsize=13, fontweight='bold')
ax2.legend(facecolor='#1a1025', edgecolor='#4c1d95', labelcolor='#c4b5fd', fontsize=10)
ax2.axhline(y=0, color='#ef4444', linewidth=1.5, linestyle='--')

# --- Plot 3: pi0 decay rate vs Nc ---
ax3 = axes[1, 0]
ax3.set_facecolor('#1a1025')
for spine in ax3.spines.values():
    spine.set_color('#4c1d95')
ax3.tick_params(colors='#c4b5fd')
ax3.grid(True, color='#4c1d95', alpha=0.3)

Nc_vals = np.arange(1, 8)
Gamma_vs_Nc = []
for nc in Nc_vals:
    A = nc * ((2.0/3)**2 - (1.0/3)**2)
    G = (alpha_em**2 * m_pi0**3) / (64 * np.pi**3 * f_pi**2) * A**2
    Gamma_vs_Nc.append(G)

Gamma_exp = hbar / 8.5e-17  # experimental decay rate in GeV

ax3.bar(Nc_vals, np.array(Gamma_vs_Nc) * 1e9, color='#a78bfa', alpha=0.8, label='Theory')
ax3.axhline(y=Gamma_exp * 1e9, color='#ef4444', linewidth=2, linestyle='--', label='Experiment')
ax3.set_xlabel('Number of colors Nc', color='#a78bfa', fontsize=11)
ax3.set_ylabel('Gamma(pi0 -> gamma gamma) [10^-9 GeV]', color='#a78bfa', fontsize=11)
ax3.set_title('pi0 Decay Rate vs Number of Colors', color='#e9d5ff', fontsize=13, fontweight='bold')
ax3.legend(facecolor='#1a1025', edgecolor='#4c1d95', labelcolor='#c4b5fd', fontsize=10)

# --- Plot 4: Spectral flow (eigenvalue shift with gauge field) ---
ax4 = axes[1, 1]
ax4.set_facecolor('#1a1025')
for spine in ax4.spines.values():
    spine.set_color('#4c1d95')
ax4.tick_params(colors='#c4b5fd')
ax4.grid(True, color='#4c1d95', alpha=0.3)

colors_spec = ['#22d3ee', '#a78bfa', '#f472b6', '#fbbf24', '#34d399']
for j, A_str in enumerate([0.0, 0.5, 1.0, 2.0, 3.0]):
    A_field = A_str * np.exp(-x**2 / 2)
    H = np.zeros((N_sites, N_sites))
    for i in range(N_sites):
        if i < N_sites - 1:
            H[i, i+1] = -1.0 / (2 * dx)
            H[i+1, i] = 1.0 / (2 * dx)
        H[i, i] = A_field[i]
    eigenvalues = np.sort(np.linalg.eigvalsh(H))
    # Plot low-lying eigenvalues
    low_idx = (np.abs(eigenvalues) < 3)
    eig_low = eigenvalues[low_idx]
    ax4.plot(np.full_like(eig_low, A_str), eig_low, '.', color=colors_spec[j],
             markersize=3, alpha=0.7)

ax4.set_xlabel('Background field strength', color='#a78bfa', fontsize=11)
ax4.set_ylabel('Eigenvalue', color='#a78bfa', fontsize=11)
ax4.set_title('Spectral Flow of Dirac Operator', color='#e9d5ff', fontsize=13, fontweight='bold')

plt.tight_layout(pad=2.0)

buf = BytesIO()
plt.savefig(buf, format='png', dpi=150, bbox_inches='tight', facecolor='#0f0a1a')
buf.seek(0)
img_b64 = base64.b64encode(buf.read()).decode()
print(f'<img src="data:image/png;base64,{img_b64}" alt="Anomalies and Triangle Diagrams" />')
buf.close()

print()
print("Key result: Anomaly cancellation in the SM requires the existence of")
print("complete generations (quarks + leptons). The pi0 -> gamma gamma rate")
print("provides direct experimental evidence for Nc = 3 colors.")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Summary: Anomalies

Chiral Anomaly

The axial current acquires an anomalous divergence $\partial_\mu j^{5\mu} = (e^2/16\pi^2)F\tilde{F}$. This is an exact one-loop result (Adler-Bardeen theorem) arising from the impossibility of simultaneously preserving vector and axial-vector gauge invariance.

Fujikawa Method

The anomaly has a beautiful path integral interpretation: under chiral rotations, the fermion measure is not invariant. The Jacobian is determined by the index of the Dirac operator, connecting the anomaly to the Atiyah-Singer index theorem and topology.

Anomaly Cancellation

All gauge anomalies cancel in the Standard Model, but only because of the precise quantum number assignments of quarks and leptons. This remarkable cancellation requires complete generations and hints at grand unification.

$\pi^0 \to \gamma\gamma$ Decay

The anomaly predicts $\Gamma(\pi^0 \to \gamma\gamma) \propto N_c^2 \alpha^2 m_\pi^3 / f_\pi^2$, in excellent agreement with experiment for $N_c = 3$. This provides direct evidence for three quark colors and is one of the most striking confirmations of quantum field theory.

← Previous: Running Couplings & RG Next: Part III — Non-Abelian Gauge Theories →