Gauge Connections

A connection on a principal bundle is the geometric object underlying all gauge theories in physics. The connection 1-form generalizes the electromagnetic potential, the curvature 2-form generalizes the field strength, and gauge transformations are the symmetries that leave physics invariant.

Historical Context

Charles Ehresmann formalized the notion of a connection on a principal bundle in 1950, unifying the Levi-Civita connection of Riemannian geometry with the gauge connections of physics. Meanwhile, Chen-Ning Yang and Robert Mills (1954) independently discovered non-abelian gauge theories from the physics side, seeking to generalize electromagnetism to the strong nuclear force.

The deep connection between these mathematical and physical developments was recognized in the 1970s, when physicists learned that the Yang-Mills field is a connection on a principal bundle and the field strength is its curvature. This unification was championed by Wu and Yang (1975) in their "dictionary" between mathematics and physics terminology.

The non-abelian Stokes theorem, instanton solutions (Belavin, Polyakov, Schwartz, Tyupkin, 1975), and the Gribov ambiguity showed that the global geometry of the bundle space has profound physical consequences for quantum field theory.

Derivation 1: The Connection 1-Form

A connection on a principal $G$-bundle $P \to M$ is a$\mathfrak{g}$-valued 1-form $\omega \in \Omega^1(P, \mathfrak{g})$ satisfying:

$R_g^*\omega = \text{Ad}_{g^{-1}}\omega \quad \text{(equivariance under right G-action)}$

$\omega(\tilde{X}) = X \quad \text{for all } X \in \mathfrak{g} \text{ (reproduces generators on vertical vectors)}$

Local Gauge Potential

Given a local section $\sigma: U \to P$ (a gauge choice), the local gauge potential is the pullback:

$A = \sigma^*\omega \in \Omega^1(U, \mathfrak{g})$

In coordinates: $A = A_\mu^a T_a\,dx^\mu$ where $T_a$ are generators of $\mathfrak{g}$. For electromagnetism ($G = U(1)$), this is simply $A = A_\mu dx^\mu$—the 4-potential.

Horizontal distribution: The connection defines a splitting$T_pP = H_pP \oplus V_pP$ into horizontal and vertical subspaces, with$H_pP = \ker\omega_p$. Horizontal lifts of curves in $M$ to $P$ define parallel transport in the bundle.

Derivation 2: The Curvature 2-Form

The curvature of a connection $\omega$ is the$\mathfrak{g}$-valued 2-form on $P$:

$\Omega = d\omega + \frac{1}{2}[\omega, \omega]$

This is the structure equation (Cartan). The local field strength is:

$\boxed{F = dA + A \wedge A = \sigma^*\Omega}$

In components:

$F_{\mu\nu} = \partial_\mu A_\nu - \partial_\nu A_\mu + [A_\mu, A_\nu]$

For abelian $G = U(1)$: the commutator vanishes and $F = dA$ is the electromagnetic field tensor. For non-abelian $G$: the $[A,A]$ term produces self-interaction of gauge bosons.

Bianchi Identity

$D\Omega = d\Omega + [\omega, \Omega] = 0$

Locally: $D_\mu F_{\nu\rho} + D_\nu F_{\rho\mu} + D_\rho F_{\mu\nu} = 0$. For electromagnetism, this gives two of Maxwell's equations: $\nabla \cdot \mathbf{B} = 0$ and Faraday's law.

Derivation 3: Gauge Transformations

A gauge transformation is a change of local section$\sigma \to \sigma' = \sigma \cdot g$ where $g: U \to G$. Under this change:

$A \to A^g = g^{-1}Ag + g^{-1}dg$

Connection transforms inhomogeneously

$F \to F^g = g^{-1}Fg$

Curvature transforms homogeneously (adjoint representation)

Derivation

Starting from $A' = (\sigma')^*\omega = (\sigma \cdot g)^*\omega$, use the equivariance of $\omega$:

$(\sigma \cdot g)^*\omega = \text{Ad}_{g^{-1}}(\sigma^*\omega) + g^*\theta$

where $\theta = g^{-1}dg$ is the Maurer-Cartan form on $G$. This gives the transformation law $A' = g^{-1}Ag + g^{-1}dg$.

Physical gauge invariance: Observable quantities must be gauge-invariant. The Yang-Mills action $S = \frac{1}{4g^2}\int \text{Tr}(F \wedge *F)$ is gauge-invariant because $\text{Tr}(g^{-1}Fg \wedge *g^{-1}Fg) = \text{Tr}(F \wedge *F)$ by the cyclicity of trace.

Derivation 4: The Yang-Mills Equations

The Yang-Mills action on a 4-manifold $(M, g)$ is:

$S_{YM}[A] = -\frac{1}{2g^2}\int_M \text{Tr}(F \wedge *F)$

Varying with respect to $A$ yields the Yang-Mills equations:

$\boxed{D*F = D_\mu F^{\mu\nu} = 0}$

Combined with the Bianchi identity $DF = 0$, these are the non-abelian generalization of Maxwell's equations. For $G = U(1)$, they reduce exactly to source-free Maxwell equations.

Self-Dual and Anti-Self-Dual Connections

In 4D, the Hodge star satisfies $**F = F$ on 2-forms, so $F$ splits into self-dual and anti-self-dual parts. A connection satisfying:

$*F = \pm F \quad \text{(instantons / anti-instantons)}$

automatically satisfies the Yang-Mills equations (since $D*F = \pm DF = 0$ by Bianchi) and minimizes the action in its topological class.

Derivation 5: Wilson Loops and Holonomy

The holonomy of a connection around a closed loop$\gamma$ is defined by parallel transport:

$\text{Hol}(\gamma) = \mathcal{P}\exp\left(-\oint_\gamma A\right) \in G$

where $\mathcal{P}$ denotes path ordering. The Wilson loop is the trace:

$W_R(\gamma) = \text{Tr}_R\,\mathcal{P}\exp\left(-\oint_\gamma A\right)$

in representation $R$. This is gauge-invariant and captures the physical content of the gauge field.

Non-Abelian Stokes Theorem

For a small loop bounding a surface element $\Sigma$:

$\text{Hol}(\gamma) \approx \mathbf{1} - F_{\mu\nu}\,\delta\Sigma^{\mu\nu} + O(\delta\Sigma^2)$

The curvature is the infinitesimal holonomy. For abelian $U(1)$, this simplifies to Stokes's theorem: $\oint A = \int F$ (Aharonov-Bohm effect).

Confinement criterion (Wilson, 1974): In a confining gauge theory (like QCD), the Wilson loop obeys an area law: $\langle W(\gamma)\rangle \sim e^{-\sigma A}$where $\sigma$ is the string tension and $A$ the minimal area bounded by $\gamma$. This signals that separating quarks costs energy proportional to distance (linear confinement).

Interactive Simulation

This simulation computes the Dirac monopole connection and field strength, the BPST instanton profile for various size parameters, and the Wilson loop phase for a U(1) gauge field. It verifies flux quantization and the topological nature of the monopole charge.

Gauge Connections: Connection 1-Forms, Curvature & Wilson Loops

Python

script.py182 lines

#!/usr/bin/env python3
"""
gauge_connections.py - Connection 1-Forms, Curvature 2-Forms, and Gauge Transformations
Computes gauge field strengths, Wilson loops, and Yang-Mills configurations
"""
import numpy as np
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt
import base64
from io import BytesIO

# ============================================================
# 1. U(1) connection (electromagnetic potential)
# ============================================================
print("=" * 72)
print("  GAUGE CONNECTIONS: THE GEOMETRY OF FORCES")
print("=" * 72)
print()

print("  U(1) CONNECTION: ELECTROMAGNETISM")
print("  " + "-" * 50)
print()
print("  Connection 1-form: A = A_mu dx^mu")
print("  Curvature 2-form:  F = dA = (1/2) F_{mu nu} dx^mu ^ dx^nu")
print("  F_{mu nu} = partial_mu A_nu - partial_nu A_mu")
print()

# Magnetic monopole: A = g(1 - cos(theta)) dphi (Dirac string at south pole)
theta_vals = np.linspace(0.01, np.pi - 0.01, 200)
A_phi_north = (1 - np.cos(theta_vals))  # g=1
A_phi_south = -(1 + np.cos(theta_vals))  # south chart

print("  Dirac monopole (north chart):")
print("    A = (1 - cos(theta)) dphi")
print("  Dirac monopole (south chart):")
print("    A = -(1 + cos(theta)) dphi")
print("  Gauge transformation on overlap:")
print("    A_N - A_S = 2 dphi = d(2*phi)")
print("    g_{NS} = exp(2i*phi)")
print("    Winding number = 1 (Dirac quantization)")
print()

# Field strength
B_r = 1.0 / (np.sin(theta_vals) + 1e-10)  # B = sin(theta) dtheta ^ dphi / r^2
# Actually on S^2: F = sin(theta) dtheta ^ dphi
F_component = np.sin(theta_vals)
flux = np.trapz(F_component, theta_vals) * 2 * np.pi
print("  Magnetic flux through S^2:")
print("    Phi = integral F = {:.6f}".format(flux))
print("    = 4*pi = {:.6f}  (first Chern number c_1 = 1)".format(4*np.pi))
print()

# ============================================================
# 2. Non-abelian connection (Yang-Mills)
# ============================================================
print("  NON-ABELIAN CONNECTION: YANG-MILLS THEORY")
print("  " + "-" * 50)
print()
print("  Connection: A = A_mu^a T_a dx^mu (Lie algebra valued)")
print("  Curvature:  F = dA + A ^ A")
print("    F_{mu nu} = partial_mu A_nu - partial_nu A_mu + [A_mu, A_nu]")
print()
print("  For SU(2): A_mu = A_mu^a (sigma_a / 2i)")
print("  Generators satisfy [T_a, T_b] = epsilon_{abc} T_c")
print()

# SU(2) BPST instanton profile function
r_vals = np.linspace(0.01, 10.0, 300)
rho = 1.0  # instanton size parameter
f_inst = r_vals**2 / (r_vals**2 + rho**2)

print("  BPST Instanton (SU(2) on R^4):")
print("    A_mu^a = f(r) * eta^a_{mu nu} x^nu / r^2")
print("    f(r) = r^2 / (r^2 + rho^2)")
print("    Action: S = 8*pi^2 / g^2 (one instanton)")
print()

for rv in [0.5, 1.0, 2.0, 5.0]:
    fv = rv**2 / (rv**2 + rho**2)
    print("    r = {:.1f}: f(r) = {:.6f}".format(rv, fv))
print()

# ============================================================
# 3. Gauge transformations
# ============================================================
print("  GAUGE TRANSFORMATIONS")
print("  " + "-" * 50)
print()
print("  Under g: M -> G:")
print("    A -> g^{-1} A g + g^{-1} dg  (connection transforms)")
print("    F -> g^{-1} F g              (curvature transforms homogeneously)")
print()
print("  Gauge-invariant quantities:")
print("    Tr(F ^ *F)  (Yang-Mills action density)")
print("    Tr(F ^ F)   (topological charge density)")
print("    Wilson loop: W(C) = Tr P exp(oint_C A)")
print()

# ============================================================
# 4. Wilson loop computation
# ============================================================
print("  WILSON LOOP FOR U(1)")
print("  " + "-" * 50)
print()

# For constant magnetic field B through a loop of radius R
B_field = 1.0
R_loop_vals = np.linspace(0.1, 5.0, 200)
flux_vals = np.pi * R_loop_vals**2 * B_field
wilson = np.cos(flux_vals)  # Re(exp(i*flux)) for U(1)

print("  Uniform B field, circular Wilson loop of radius R:")
print("    W(C) = exp(i * B * pi * R^2)")
print("    |W| = 1  (unitary)")
print()
for rv in [0.5, 1.0, 2.0, 3.0]:
    phi_val = np.pi * rv**2 * B_field
    print("    R = {:.1f}: flux = {:.4f}, W = exp(i * {:.4f})".format(rv, phi_val, phi_val))
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: Monopole connection
ax1 = axes[0, 0]
ax1.plot(theta_vals, A_phi_north, color='#2dd4bf', linewidth=2.5, label='A_N (north chart)')
ax1.plot(theta_vals, A_phi_south, color='#f472b6', linewidth=2.5, label='A_S (south chart)')
ax1.axvline(x=np.pi/2, color='#f59e0b', linestyle='--', alpha=0.6, label='equator')
ax1.set_xlabel('theta', color='#5eead4')
ax1.set_ylabel('A_phi', color='#5eead4')
ax1.set_title('Dirac Monopole Connection', color='#99f6e4', fontsize=12, fontweight='bold')
ax1.legend(facecolor='#111111', edgecolor='#134e4a', labelcolor='#99f6e4', fontsize=8)

# Plot 2: Curvature (field strength)
ax2 = axes[0, 1]
ax2.plot(theta_vals, F_component, color='#34d399', linewidth=2.5)
ax2.fill_between(theta_vals, F_component, alpha=0.15, color='#34d399')
ax2.set_xlabel('theta', color='#5eead4')
ax2.set_ylabel('F_{theta phi}', color='#5eead4')
ax2.set_title('Monopole Field Strength on S^2', color='#99f6e4', fontsize=12, fontweight='bold')

# Plot 3: BPST instanton profile
ax3 = axes[1, 0]
for rho_val in [0.5, 1.0, 2.0, 4.0]:
    f_val = r_vals**2 / (r_vals**2 + rho_val**2)
    ax3.plot(r_vals, f_val, linewidth=2,
             label='rho = {:.1f}'.format(rho_val))
ax3.set_xlabel('r', color='#5eead4')
ax3.set_ylabel('f(r) = r^2/(r^2 + rho^2)', color='#5eead4')
ax3.set_title('BPST Instanton Profile', color='#99f6e4', fontsize=12, fontweight='bold')
ax3.legend(facecolor='#111111', edgecolor='#134e4a', labelcolor='#99f6e4', fontsize=8)

# Plot 4: Wilson loop phase
ax4 = axes[1, 1]
ax4.plot(R_loop_vals, np.cos(flux_vals), color='#2dd4bf', linewidth=2, label='Re(W)')
ax4.plot(R_loop_vals, np.sin(flux_vals), color='#f472b6', linewidth=2, label='Im(W)')
ax4.set_xlabel('Loop radius R', color='#5eead4')
ax4.set_ylabel('Wilson loop', color='#5eead4')
ax4.set_title('U(1) Wilson Loop in Uniform B', 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="Gauge Connections 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

Connection 1-Form

A Lie-algebra-valued 1-form on the principal bundle that defines parallel transport and horizontal subspaces. Locally: the gauge potential $A = A_\mu^a T_a dx^\mu$.

Curvature 2-Form

$F = dA + A \wedge A$: the field strength. Measures non-commutativity of parallel transport. The Bianchi identity $DF = 0$ gives half of the generalized Maxwell equations.

Gauge Transformations

Changes of local trivialization. The connection transforms inhomogeneously ($A \to g^{-1}Ag + g^{-1}dg$) but the curvature transforms covariantly.

Wilson Loops

Gauge-invariant observables defined by the trace of parallel transport around closed loops. They encode the physical content of gauge theories and diagnose confinement.

← Fiber Bundles Next: Characteristic Classes →