Lattice and nonperturbative methods
Understand a lattice regulator, extract a mass from Euclidean decay, and identify statistical and systematic limits.
On this page
The question for today
What you’ll understand
- Derive the leading lattice correction to a free kinetic term.
- Extract a ground-state energy from a Euclidean correlator under explicit assumptions.
- Separate continuum, volume, sampling, and sign-problem limitations.
Strong coupling does not make a theory meaningless; it makes a particular expansion unhelpful. A Euclidean lattice supplies a finite regulator with which many nonperturbative observables can be estimated. The task is then to remove or control the regulator and sampling errors. A smooth animation is not evidence for a continuum limit, so we will derive a simple discretization error explicitly.
Spacing is a regulator, not a physical crystal by default
Replace continuum points by sites separated by spacing and derivatives by finite differences. For a free scalar with a nearest-neighbor kinetic term, each momentum component enters through . Expanding the sine at small recovers the continuum momentum with corrections of order . This is a statement about a specified discretization; improved actions can cancel leading errors, and interacting observables can have additional logarithmic factors.
Gauge theories place group-valued parallel transporters on links rather than simply replacing every gauge potential by an unconstrained site vector. Products around elementary plaquettes give gauge-invariant local terms in the Wilson action. Its conventional parameter for SU(N) is a bare lattice parameter, not the renormalization-group beta function. To approach a continuum theory one tunes bare parameters while holding chosen physical quantities fixed as decreases.
Euclidean time filters the spectrum
Insert energy eigenstates between an operator and its adjoint in a zero-temperature Euclidean correlator. After subtracting any vacuum expectation value and setting the vacuum energy to zero, the connected correlator is a sum of decaying exponentials. If the operator overlaps the lightest state in the selected channel, sufficiently long Euclidean time suppresses heavier states relative to it. The decay rate therefore exposes an energy gap without identifying a literal particle trajectory through sampled fields.
Taking the logarithm of neighboring correlator values defines an effective mass. A plateau can indicate ground-state dominance, but it must be assessed with correlated statistical errors and changes of the fitting window or operator basis. Finite temporal extent introduces backward-propagating contributions, often requiring a hyperbolic-cosine fit. At long times the signal can decay faster than its uncertainty, so simply moving to later times is not a cost-free cure for excited-state contamination.
What sampling can and cannot guarantee
When the Euclidean weight is nonnegative, importance sampling can estimate expectation values without enumerating every field configuration. Consecutive samples are correlated, so the effective sample count can be much smaller than the number stored. Thermalization, autocorrelation, topology freezing, finite spatial volume, and finite lattice spacing all require checks. Statistical precision improves with more suitable samples, but it does not remove a systematic error caused by the wrong volume or an unextrapolated regulator.
At nonzero baryon chemical potential and in other settings, the weight can become complex or change sign. Sampling its absolute value and reweighting leaves a small average phase that can fall exponentially with volume, producing the sign problem. There is no general efficient classical remedy for arbitrary theories. Tensor networks, quantum simulation, contour methods, and specialized expansions address particular regimes with their own limits. None should be advertised as a universal replacement for continuum extrapolation and a transparent error budget.
Work it through
Extract a gap from an ideal single exponential
Suppose a connected correlator is exactly C(τ) = A exp(−mτ), and C(τ+a)/C(τ) = exp(−0.4). What gap is determined?
Insert the single exponential into the ratio. The unknown overlap A cancels.
Take the logarithm of the inverse ratio to obtain am = 0.4.
The dimensionful mass is m = 0.4/a. Converting to physical units requires an independently calibrated lattice spacing; the ratio alone measures only the dimensionless lattice gap.
A perfect exponential gives an exact effective-mass plateau in this toy case. Real data need excited-state, finite-time, scale-setting, and uncertainty analyses before a physical mass is quoted.
The living notebook
Explore the model
Interactive lab · spectral
Observe what a finite mode cutoff omits
Increase N in the finite sine sum and compare the result at two decay powers p. Higher modes add shorter spatial features. This illustrates resolution and truncation only; it is not a Monte Carlo chain, a lattice gauge action, a continuum extrapolation, or a simulated sign problem.
- Sum at x = π/2
- 0.867
- Smallest wavelength
- 2π / 5
- Highest-mode coefficient
- 0.2
A finite sum of sin(nx − t)/nᵖ. Playback changes t; increasing the cutoff includes finer spatial scales. This is not a proof about a continuum limit or renormalization. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.
Look forRegulator refinement changes resolution, but a physical continuum claim requires matched parameters and convergence evidence.
Pause & reason
Check your understanding
Your turn
Try it yourself
For one small momentum component p, use the free kinetic expansion to estimate the fractional leading error. What happens if a is halved at fixed physical p?
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
Divide the correction by p² before changing a.
Compare with the solution
The leading fractional change is . Halving reduces its magnitude by a factor of four in this small- nearest-neighbor example. An interacting observable need not follow this exact coefficient, so the scaling must still be tested with actual matched calculations.
Keep the idea
Take it with you
- The continuum limit changes the regulator while preserving calibrated physics.
- Euclidean decay encodes energy gaps subject to overlap and finite-time assumptions.
- Sampling error, regulator bias, finite volume, and the sign problem are different obstacles.
Read further
Original sources and deeper explanations to follow at your own pace.
Optional progress, saved in this browser when available. Revisit any lesson at any time.
Lesson not marked complete.