Chapter 1: Classical Field Theory

Before quantizing fields, we must understand them classically. Classical field theory extends the Lagrangian mechanics of particles to continuous systems described by fields $\phi(x^\mu)$ defined at every spacetime point. The action principle, Noether's theorem, and the equations of motion for scalar, spinor, and vector fields form the foundation upon which all of quantum field theory is built.

Derivation 1: Lagrangian Field Theory

In classical mechanics, the dynamics of a system with $N$ degrees of freedom is governed by the Lagrangian $L(q_i, \dot{q}_i, t)$. For a field, we promote discrete coordinates to a continuous field $\phi(x^\mu)$ defined over all spacetime.

The Action Principle

The action functional is defined as the integral of the Lagrangian density $\mathcal{L}$ over all spacetime:

$S[\phi] = \int d^4x \, \mathcal{L}(\phi, \partial_\mu \phi)$

Here $d^4x = dt \, dx \, dy \, dz$ and $\mathcal{L}$ depends on the field and its first derivatives $\partial_\mu \phi = \frac{\partial \phi}{\partial x^\mu}$. The Lagrangian density replaces the discrete sum over degrees of freedom with a continuous integral:

From particles to fields: In the continuum limit, discrete coordinates $q_i(t)$ become $\phi(\mathbf{x}, t)$, where $\mathbf{x}$ plays the role of the label $i$. The Lagrangian becomes:

$L(t) = \int d^3x \, \mathcal{L}(\phi, \partial_\mu \phi)$

Euler-Lagrange Equations for Fields

We demand that the action is stationary under small variations $\phi \to \phi + \delta\phi$, where $\delta\phi$ vanishes on the boundary. Computing $\delta S = 0$:

Step 1: Vary the action

$\delta S = \int d^4x \left[ \frac{\partial \mathcal{L}}{\partial \phi} \delta\phi + \frac{\partial \mathcal{L}}{\partial(\partial_\mu \phi)} \delta(\partial_\mu \phi) \right]$

Step 2: Since $\delta(\partial_\mu \phi) = \partial_\mu(\delta\phi)$, integrate by parts:

$\delta S = \int d^4x \left[ \frac{\partial \mathcal{L}}{\partial \phi} - \partial_\mu \frac{\partial \mathcal{L}}{\partial(\partial_\mu \phi)} \right] \delta\phi + \text{boundary terms}$

Step 3: Boundary terms vanish; since $\delta\phi$ is arbitrary:

$\boxed{\frac{\partial \mathcal{L}}{\partial \phi} - \partial_\mu \left( \frac{\partial \mathcal{L}}{\partial(\partial_\mu \phi)} \right) = 0}$

This is the Euler-Lagrange equation for fields — the field-theoretic analog of $\frac{d}{dt}\frac{\partial L}{\partial \dot{q}} = \frac{\partial L}{\partial q}$. It determines the classical dynamics of any field given its Lagrangian density.

Conjugate Momentum and Hamiltonian Density

The canonical momentum conjugate to $\phi$ is:

$\pi(x) = \frac{\partial \mathcal{L}}{\partial \dot{\phi}(x)}$

The Hamiltonian density is obtained by Legendre transformation:

$\mathcal{H} = \pi \dot{\phi} - \mathcal{L}$

Derivation 2: Noether's Theorem for Fields

Emmy Noether's theorem is one of the most profound results in theoretical physics: every continuous symmetry of the action corresponds to a conserved current. This connects the symmetries of nature to its conservation laws.

Statement and Proof

Consider an infinitesimal transformation that leaves the action invariant:

$\phi(x) \to \phi(x) + \alpha \, \Delta\phi(x)$

where $\alpha$ is an infinitesimal parameter. If the Lagrangian changes by at most a total divergence $\mathcal{L} \to \mathcal{L} + \alpha \, \partial_\mu K^\mu$, then we can construct a conserved current. The variation of $\mathcal{L}$ is:

$\delta\mathcal{L} = \frac{\partial \mathcal{L}}{\partial \phi} \delta\phi + \frac{\partial \mathcal{L}}{\partial(\partial_\mu \phi)} \partial_\mu(\delta\phi)$

Using the Euler-Lagrange equation for the first term, this becomes:

$\delta\mathcal{L} = \partial_\mu \left( \frac{\partial \mathcal{L}}{\partial(\partial_\mu \phi)} \delta\phi \right)$

Setting $\delta\mathcal{L} = \alpha \, \partial_\mu K^\mu$ and $\delta\phi = \alpha \, \Delta\phi$, we define the Noether current:

$\boxed{j^\mu = \frac{\partial \mathcal{L}}{\partial(\partial_\mu \phi)} \Delta\phi - K^\mu}$

with $\partial_\mu j^\mu = 0$ (conservation law)

The conserved charge is:

$Q = \int d^3x \, j^0(x), \qquad \frac{dQ}{dt} = 0$

Energy-Momentum Tensor

Spacetime translations $x^\mu \to x^\mu + a^\mu$ give $\Delta\phi = -a^\nu \partial_\nu \phi$and $\delta\mathcal{L} = -a^\nu \partial_\nu \mathcal{L} = -a^\nu \partial_\mu(\delta^\mu_\nu \mathcal{L})$. With $K^\mu = -a^\nu \delta^\mu_\nu \mathcal{L}$, the four conserved currents are:

$\boxed{T^{\mu\nu} = \frac{\partial \mathcal{L}}{\partial(\partial_\mu \phi)} \partial^\nu \phi - \eta^{\mu\nu} \mathcal{L}}$

$\partial_\mu T^{\mu\nu} = 0$ — the canonical energy-momentum tensor

Conserved charges from $T^{\mu\nu}$:

• Energy: $E = \int d^3x \, T^{00}$ (from time translation invariance)
• Momentum: $P^i = \int d^3x \, T^{0i}$ (from spatial translation invariance)

Angular Momentum

Under Lorentz transformations $x^\mu \to \Lambda^\mu_{\ \nu} x^\nu$, with infinitesimal parameter $\omega_{\rho\sigma}$, the conserved current leads to the angular momentum tensor:

$M^{\mu\rho\sigma} = x^\rho T^{\mu\sigma} - x^\sigma T^{\mu\rho} + S^{\mu\rho\sigma}$

The first two terms give the orbital angular momentum, and $S^{\mu\rho\sigma}$ is the intrinsic spin angular momentum. For scalar fields $S^{\mu\rho\sigma} = 0$; for spinor fields, this term gives rise to the spin-1/2 property.

Derivation 3: The Klein-Gordon Equation

The Klein-Gordon equation describes a free relativistic scalar field. It is the simplest relativistic wave equation, governing spin-0 particles like the Higgs boson.

Lagrangian for a Real Scalar Field

We write the Lagrangian density for a free real scalar field of mass $m$:

$\mathcal{L} = \frac{1}{2}(\partial_\mu \phi)(\partial^\mu \phi) - \frac{1}{2}m^2 \phi^2$

In the mostly-minus metric convention $\eta^{\mu\nu} = \text{diag}(+1, -1, -1, -1)$:

$(\partial_\mu \phi)(\partial^\mu \phi) = \dot{\phi}^2 - |\nabla\phi|^2$

This ensures Lorentz invariance: the kinetic term is a scalar under Lorentz transformations.

Deriving the Equation of Motion

We apply the Euler-Lagrange equation. Computing each term:

Term 1:

$\frac{\partial \mathcal{L}}{\partial \phi} = -m^2 \phi$

Term 2: Using the chain rule on the kinetic term,

$\frac{\partial \mathcal{L}}{\partial(\partial_\mu \phi)} = \partial^\mu \phi$

Therefore:

$\partial_\mu \frac{\partial \mathcal{L}}{\partial(\partial_\mu \phi)} = \partial_\mu \partial^\mu \phi = \Box \phi$

Combining via the Euler-Lagrange equation:

$\boxed{(\Box + m^2)\phi = 0} \quad \text{where} \quad \Box = \partial_\mu \partial^\mu = \frac{\partial^2}{\partial t^2} - \nabla^2$

Plane Wave Solutions

We seek solutions of the form $\phi(x) = A \, e^{-ik \cdot x}$ where $k \cdot x = k_\mu x^\mu = \omega t - \mathbf{k} \cdot \mathbf{x}$. Substituting:

$(-k_\mu k^\mu + m^2) \phi = 0 \implies k_\mu k^\mu = m^2$

$\omega^2 = |\mathbf{k}|^2 + m^2$

This is the relativistic dispersion relation. Note that it admits both positive and negative frequency solutions:

$\omega_{\mathbf{k}} = \pm \sqrt{|\mathbf{k}|^2 + m^2}$

The general solution is a superposition over all momenta:

$\phi(x) = \int \frac{d^3k}{(2\pi)^3} \frac{1}{\sqrt{2\omega_{\mathbf{k}}}} \left[ a(\mathbf{k}) \, e^{-ik \cdot x} + a^*(\mathbf{k}) \, e^{ik \cdot x} \right]$

The requirement $\phi = \phi^*$ (real field) forces the coefficients of positive and negative frequency modes to be complex conjugates.

Energy-Momentum Tensor for Klein-Gordon

Applying our general formula for $T^{\mu\nu}$:

$T^{\mu\nu} = \partial^\mu \phi \, \partial^\nu \phi - \eta^{\mu\nu} \left[ \frac{1}{2}(\partial_\alpha \phi)(\partial^\alpha \phi) - \frac{1}{2}m^2\phi^2 \right]$

The energy density is:

$T^{00} = \frac{1}{2}\dot{\phi}^2 + \frac{1}{2}|\nabla\phi|^2 + \frac{1}{2}m^2\phi^2$

All three terms are positive-definite, ensuring the energy is bounded below. This is a crucial requirement for stability of the theory.

Derivation 4: The Dirac Equation

The Klein-Gordon equation is second-order in time, leading to difficulties with probability interpretation (negative probabilities). Dirac sought a first-order relativistic wave equation.

Dirac's Brilliant Ansatz

Start from the relativistic energy-momentum relation:

$E^2 = |\mathbf{p}|^2 c^2 + m^2 c^4$

Dirac wanted a linear equation $i\hbar \frac{\partial \psi}{\partial t} = H\psi$ with$H = c\boldsymbol{\alpha} \cdot \mathbf{p} + \beta m c^2$, where $\boldsymbol{\alpha}$ and $\beta$are to be determined. Squaring:

$H^2 = c^2 \sum_{i,j} \frac{1}{2}(\alpha_i \alpha_j + \alpha_j \alpha_i) p_i p_j + mc^3 \sum_i (\alpha_i \beta + \beta \alpha_i) p_i + \beta^2 m^2 c^4$

For $H^2 = E^2 = |\mathbf{p}|^2 c^2 + m^2 c^4$, we need:

$\{\alpha_i, \alpha_j\} = 2\delta_{ij} \mathbf{I}, \quad \{\alpha_i, \beta\} = 0, \quad \beta^2 = \mathbf{I}$

These anticommutation relations cannot be satisfied by ordinary numbers — they require matrices!

Gamma Matrices

Defining $\gamma^0 = \beta$ and $\gamma^i = \beta \alpha_i$, the anticommutation relations become the Clifford algebra:

$\boxed{\{\gamma^\mu, \gamma^\nu\} = 2\eta^{\mu\nu} \mathbf{I}_{4 \times 4}}$

The minimum dimension satisfying this algebra in 4D spacetime is $4 \times 4$ matrices. In the Dirac (standard) representation:

$\gamma^0 = \begin{pmatrix} \mathbf{I} & 0 \\ 0 & -\mathbf{I} \end{pmatrix}, \quad \gamma^i = \begin{pmatrix} 0 & \sigma^i \\ -\sigma^i & 0 \end{pmatrix}$

where $\sigma^i$ are the $2 \times 2$ Pauli matrices.

The Covariant Dirac Equation

In natural units ($c = \hbar = 1$), multiplying by $\beta/c$ and using the gamma matrices:

$\boxed{(i\gamma^\mu \partial_\mu - m)\psi = 0}$

The Dirac equation — using the Feynman slash notation: $(i\not{\partial} - m)\psi = 0$

Spinor Solutions and Antiparticles

The Dirac equation has plane-wave solutions $\psi = u(\mathbf{p}) e^{-ip \cdot x}$and $\psi = v(\mathbf{p}) e^{ip \cdot x}$. For a particle at rest ($\mathbf{p} = 0$):

Substituting into the Dirac equation gives $(i\gamma^0 \partial_0 - m)\psi = 0$:

$\gamma^0 E \, u = m \, u \implies E = \pm m$

The four-component spinor has two positive-energy solutions (spin up/down) and two negative-energy solutions. Dirac interpreted the negative-energy states as the Dirac sea, whose holes are antiparticles (positrons). The modern interpretation uses the Feynman-Stuckelberg approach: negative-energy solutions propagating backward in time are antiparticles propagating forward.

Positive-energy spinors ($E = +\sqrt{|\mathbf{p}|^2 + m^2}$):

$u(\mathbf{p}) = \begin{pmatrix} \sqrt{E + m} \, \chi \\ \frac{\boldsymbol{\sigma} \cdot \mathbf{p}}{\sqrt{E + m}} \, \chi \end{pmatrix}$

where $\chi$ is a two-component spinor (spin up or down).

At low momentum ($|\mathbf{p}| \ll m$), the lower components are suppressed by $\sim |\mathbf{p}|/m$ — recovering the non-relativistic Pauli equation. In the ultra-relativistic limit ($|\mathbf{p}| \gg m$), upper and lower components become equal, and chirality becomes a good quantum number.

Derivation 5: Maxwell Equations from a Lagrangian

Maxwell's equations, the crown jewel of classical electrodynamics, emerge naturally from a gauge-invariant Lagrangian for the vector potential $A^\mu$.

The Field Strength Tensor

Define the electromagnetic field strength tensor:

$F_{\mu\nu} = \partial_\mu A_\nu - \partial_\nu A_\mu$

This antisymmetric tensor encodes the electric and magnetic fields:

$F_{0i} = -E_i$ (electric field), $F_{ij} = -\epsilon_{ijk} B_k$ (magnetic field)

The Maxwell Lagrangian

$\mathcal{L} = -\frac{1}{4} F_{\mu\nu} F^{\mu\nu}$

Expanding in terms of $\mathbf{E}$ and $\mathbf{B}$:

$\mathcal{L} = \frac{1}{2}(\mathbf{E}^2 - \mathbf{B}^2)$

The electric field contributes kinetic energy; the magnetic field contributes potential energy.

Deriving Maxwell's Equations

Apply the Euler-Lagrange equation with $A_\nu$ as the field variable:

Since $\mathcal{L}$ does not depend on $A_\nu$ directly (only through $\partial_\mu A_\nu$):

$\frac{\partial \mathcal{L}}{\partial A_\nu} = 0$

Computing the derivative with respect to $\partial_\mu A_\nu$:

$\frac{\partial \mathcal{L}}{\partial(\partial_\mu A_\nu)} = -F^{\mu\nu}$

The Euler-Lagrange equation gives:

$\boxed{\partial_\mu F^{\mu\nu} = 0}$

These are the source-free Maxwell equations ($\nabla \cdot \mathbf{E} = 0$ and $\nabla \times \mathbf{B} = \frac{\partial \mathbf{E}}{\partial t}$)

The other two Maxwell equations ($\nabla \cdot \mathbf{B} = 0$ and Faraday's law) follow identically from the definition of $F_{\mu\nu}$ (the Bianchi identity):

$\partial_\mu \tilde{F}^{\mu\nu} = 0 \quad \text{where} \quad \tilde{F}^{\mu\nu} = \frac{1}{2}\epsilon^{\mu\nu\rho\sigma} F_{\rho\sigma}$

Gauge Invariance

The Lagrangian is invariant under the gauge transformation:

$A_\mu \to A_\mu + \partial_\mu \Lambda(x)$

Since $F_{\mu\nu} \to F_{\mu\nu} + \partial_\mu \partial_\nu \Lambda - \partial_\nu \partial_\mu \Lambda = F_{\mu\nu}$. This gauge redundancy means $A_\mu$ has only two physical (transverse) degrees of freedom — the two polarizations of light. The photon is massless because a mass term $\frac{1}{2}m^2 A_\mu A^\mu$would break gauge invariance.

Coupling to matter: Adding a source current$j^\nu$, the Lagrangian becomes $\mathcal{L} = -\frac{1}{4}F_{\mu\nu}F^{\mu\nu} - j^\nu A_\nu$, giving $\partial_\mu F^{\mu\nu} = j^\nu$. Current conservation $\partial_\nu j^\nu = 0$follows from the antisymmetry of $F^{\mu\nu}$ and is guaranteed by gauge invariance via Noether's theorem.

Simulation: Dispersion Relations & Dirac Spinors

This simulation visualizes the Klein-Gordon dispersion relation for various masses, the mass gap, the group velocity (always subluminal), and the momentum dependence of Dirac spinor components showing the transition from non-relativistic to ultra-relativistic behavior.

Classical Field Theory: Dispersion & Spinor Analysis

Python

script.py148 lines

#!/usr/bin/env python3
"""
classical_fields.py - Klein-Gordon Dispersion & Dirac Spinor Components
Visualizes relativistic dispersion relations and spinor structure
"""
import numpy as np
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt
import base64
from io import BytesIO

# ============================================================
# Part 1: Klein-Gordon Dispersion Relation
# E^2 = p^2 c^2 + m^2 c^4  (natural units c=hbar=1: omega^2 = k^2 + m^2)
# ============================================================

m_values = [0.0, 0.5, 1.0, 2.0]
k = np.linspace(-5, 5, 500)

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

# --- Plot 1: KG dispersion for various masses ---
ax1 = axes[0, 0]
ax1.set_facecolor('#1a1025')
for spine in ax1.spines.values():
    spine.set_color('#4c1d95')
ax1.tick_params(colors='#c4b5fd')
ax1.grid(True, color='#4c1d95', alpha=0.3)

for m in m_values:
    omega_plus = np.sqrt(k**2 + m**2)
    omega_minus = -np.sqrt(k**2 + m**2)
    label_p = f'm = {m:.1f} (+)' if m > 0 else 'massless (+)'
    label_m = f'm = {m:.1f} (-)' if m > 0 else 'massless (-)'
    line, = ax1.plot(k, omega_plus, linewidth=2, label=label_p)
    ax1.plot(k, omega_minus, linewidth=2, color=line.get_color(), linestyle='--', alpha=0.6)

ax1.set_xlabel('Momentum k', color='#a78bfa', fontsize=11)
ax1.set_ylabel('Energy omega', color='#a78bfa', fontsize=11)
ax1.set_title('Klein-Gordon Dispersion: omega^2 = k^2 + m^2', color='#e9d5ff', fontsize=13, fontweight='bold')
ax1.legend(facecolor='#1a1025', edgecolor='#4c1d95', labelcolor='#c4b5fd', fontsize=9)
ax1.axhline(y=0, color='#6d28d9', linewidth=0.8)

# --- Plot 2: Mass gap visualization ---
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)

m = 1.0
omega = np.sqrt(k**2 + m**2)
ax2.fill_between(k, -m, m, alpha=0.15, color='#f97316', label=f'Mass gap |omega| < m = {m}')
ax2.plot(k, omega, color='#a78bfa', linewidth=2.5, label='Positive branch')
ax2.plot(k, -omega, color='#a78bfa', linewidth=2.5, linestyle='--', label='Negative branch')
ax2.axhline(y=m, color='#f97316', linewidth=1, linestyle=':')
ax2.axhline(y=-m, color='#f97316', linewidth=1, linestyle=':')
ax2.set_xlabel('Momentum k', color='#a78bfa', fontsize=11)
ax2.set_ylabel('Energy omega', color='#a78bfa', fontsize=11)
ax2.set_title('Mass Gap in Klein-Gordon Spectrum', color='#e9d5ff', fontsize=13, fontweight='bold')
ax2.legend(facecolor='#1a1025', edgecolor='#4c1d95', labelcolor='#c4b5fd', fontsize=9)

# --- Plot 3: Group velocity v_g = k / sqrt(k^2 + m^2) ---
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)

k_pos = np.linspace(0.01, 10, 500)
for m in [0.1, 0.5, 1.0, 2.0]:
    vg = k_pos / np.sqrt(k_pos**2 + m**2)
    ax3.plot(k_pos, vg, linewidth=2, label=f'm = {m}')

ax3.axhline(y=1.0, color='#ef4444', linewidth=1.5, linestyle='--', label='c = 1 (speed of light)')
ax3.set_xlabel('Momentum k', color='#a78bfa', fontsize=11)
ax3.set_ylabel('Group velocity v_g', color='#a78bfa', fontsize=11)
ax3.set_title('Group Velocity: v_g = d omega / dk < c', color='#e9d5ff', fontsize=13, fontweight='bold')
ax3.legend(facecolor='#1a1025', edgecolor='#4c1d95', labelcolor='#c4b5fd', fontsize=9)
ax3.set_ylim(0, 1.15)

# --- Plot 4: Dirac Spinor Components ---
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)

# Positive energy u-spinor: u(p) ~ (sqrt(E+m) * chi, sqrt(E-m) * sigma.p_hat * chi)
# For spin-up along z, chi = (1,0), p along z
m_dirac = 1.0
p = np.linspace(0, 5, 500)
E = np.sqrt(p**2 + m_dirac**2)

# Upper component amplitude |u_1|^2 + |u_2|^2 proportional to E+m
upper = np.sqrt(E + m_dirac)
# Lower component amplitude proportional to sqrt(E-m) for p along z
lower = np.sqrt(E - m_dirac)

# Normalize
norm = np.sqrt(upper**2 + lower**2)
upper_frac = (upper / norm)**2
lower_frac = (lower / norm)**2

ax4.plot(p, upper_frac, color='#22d3ee', linewidth=2.5, label='Upper (particle)')
ax4.plot(p, lower_frac, color='#f472b6', linewidth=2.5, label='Lower (antiparticle)')
ax4.axvline(x=m_dirac, color='#fbbf24', linewidth=1, linestyle=':', label=f'p = m = {m_dirac}')
ax4.set_xlabel('Momentum p', color='#a78bfa', fontsize=11)
ax4.set_ylabel('Fraction of spinor norm', color='#a78bfa', fontsize=11)
ax4.set_title('Dirac Spinor: Upper vs Lower Components', color='#e9d5ff', fontsize=13, fontweight='bold')
ax4.legend(facecolor='#1a1025', edgecolor='#4c1d95', labelcolor='#c4b5fd', fontsize=9)

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="Classical Field Theory Plots" />')
buf.close()

print()
print("=" * 72)
print("  CLASSICAL FIELD THEORY - NUMERICAL ANALYSIS")
print("=" * 72)
print()
print("Klein-Gordon Dispersion Relation: omega^2 = k^2 + m^2")
print("-" * 50)
for m in [0.0, 0.5, 1.0, 2.0]:
    print(f"  m = {m:.1f}:  omega(k=0) = {np.sqrt(m**2):.4f},  omega(k=3) = {np.sqrt(9+m**2):.4f}")
print()
print("Dirac Spinor Analysis (m = 1.0, spin-up along z):")
print("-" * 50)
for pval in [0.0, 0.5, 1.0, 2.0, 5.0]:
    Eval = np.sqrt(pval**2 + 1.0)
    up = np.sqrt(Eval + 1.0)
    lo = np.sqrt(max(Eval - 1.0, 0))
    nm = np.sqrt(up**2 + lo**2)
    print(f"  p = {pval:.1f}: E = {Eval:.4f}, upper = {(up/nm)**2:.4f}, lower = {(lo/nm)**2:.4f}")
print()
print("At rest (p=0): purely upper component (particle at rest)")
print("Ultra-relativistic (p >> m): equal upper and lower (chiral limit)")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Complex Scalar Field and Internal Symmetries

A complex scalar field $\Phi = (\phi_1 + i\phi_2)/\sqrt{2}$ has the Lagrangian:

$\mathcal{L} = (\partial_\mu \Phi^*)(\partial^\mu \Phi) - m^2 \Phi^* \Phi$

This possesses a global $U(1)$ symmetry $\Phi \to e^{i\alpha}\Phi$, whose Noether current is:

$j^\mu = i[\Phi^* (\partial^\mu \Phi) - (\partial^\mu \Phi^*) \Phi]$

The conserved charge $Q = \int d^3x \, j^0$ is interpreted as electric charge (or particle number). The complex field describes charged particles — it has distinct particle and antiparticle excitations, unlike the real scalar field.

Pattern recognition: The structure of classical field theory is deeply systematic. Given any symmetry group $G$, one can construct fields transforming in representations of $G$, build invariant Lagrangians, and derive the resulting equations of motion and conservation laws. This is the template used to construct the entire Standard Model of particle physics, where $G = SU(3) \times SU(2) \times U(1)$.

Summary: The Classical Field Theory Toolkit

Lagrangian Framework

The action principle $\delta S = 0$ yields the Euler-Lagrange equations, providing equations of motion for any field once the Lagrangian density is specified.

Noether's Theorem

Symmetries give conservation laws: translations $\to$ energy-momentum, Lorentz invariance $\to$ angular momentum, internal symmetries $\to$ charges.

Relativistic Wave Equations

Klein-Gordon for spin-0, Dirac for spin-1/2, Maxwell for spin-1. Each follows from a Lorentz-invariant Lagrangian with the appropriate field content and symmetries.

Gauge Invariance

The electromagnetic Lagrangian possesses gauge symmetry $A_\mu \to A_\mu + \partial_\mu \Lambda$, which ensures the photon is massless and constrains the form of interactions. This principle generalizes to non-Abelian gauge theories (the Standard Model).

← Course Overview Next: Canonical Quantization →