Lattice QCD Basics

Core Idea

Lattice QCD is the non-perturbative formulation of quantum chromodynamics on a discrete spacetime lattice, where the path integral is evaluated numerically using Monte Carlo methods. It provides first-principles calculations of hadron masses, decay constants, form factors, and other quantities that are inaccessible to perturbation theory. Lattice QCD is essential for extracting fundamental parameters (quark masses, CKM elements, alpha_s) from experimental data.

Explainer

Lattice QCD was proposed by Kenneth Wilson in 1974 as a way to define QCD non-perturbatively by discretizing spacetime into a four-dimensional hypercubic lattice. The continuum path integral Z = integral D[A] D[psi] D[psi-bar] exp(i*S_QCD) is replaced by a well-defined (finite-dimensional) integral over link variables U and quark fields, evaluated in Euclidean spacetime (after Wick rotation). The lattice provides both an ultraviolet cutoff (the lattice spacing a) and an infrared cutoff (the lattice volume L⁴), making the theory fully regulated.

The numerical evaluation uses importance sampling Monte Carlo: gauge field configurations are generated with probability proportional to exp(-S_lattice), and observables are estimated as averages over these configurations. The inclusion of dynamical quarks (quark loops in the vacuum) requires computing the fermion determinant, which is the most computationally expensive part. Modern algorithms (Hybrid Monte Carlo with mass preconditioning and deflation) have made full dynamical calculations routine on current supercomputers. A typical state-of-the-art calculation uses lattice spacings a = 0.06-0.12 fm, volumes (6 fm)³, physical quark masses (m_pi ~ 135 MeV), and 2+1+1 dynamical quark flavors (u, d, s, c).

The flagship results of lattice QCD include: (1) the light hadron spectrum, computed to ~1% precision and matching experiment perfectly; (2) the strong coupling constant alpha_s(M_Z) = 0.1179 +/- 0.0009, the most precise determination; (3) quark masses (m_c, m_b to sub-percent precision; m_s to ~1%; m_u, m_d to ~5%); (4) CKM matrix elements from B and K meson form factors (f_B, f_{B_s}, B -> pi l nu, B -> D(*) l nu form factors), which are essential inputs for the unitarity triangle; (5) hadronic vacuum polarization and light-by-light contributions to the muon g-2, currently a major focus due to the tension between the experimental measurement and the Standard Model prediction.

The limitations of lattice QCD are primarily in processes involving real-time dynamics (scattering amplitudes, transport coefficients), which require Minkowski spacetime and are not directly accessible in Euclidean lattice calculations. Progress has been made using methods like the Luscher finite-volume formalism (relating discrete energy levels in a box to scattering phase shifts) and the Backus-Gilbert method for spectral function reconstruction. Multi-hadron states, resonances, and transition form factors at large momentum transfer remain challenging. Despite these limitations, lattice QCD is the only systematic, improvable, first-principles method for computing non-perturbative QCD quantities, and its results underpin much of the precision flavor physics program.