Continuum scaling and universality checks
Separate sampling, stochastic-time, volume, and spatial-cutoff limits while recording enough evidence to reproduce the extrapolation.
On this page
The question for today
What you’ll understand
- Match physical length scales while refining a lattice.
- Distinguish a finite-volume continuum limit from a thermodynamic limit.
- Assess cutoff fits and universality claims with explicit assumptions and reproducible analysis.
A statistical field simulation has several limits, and increasing the configuration count addresses only one of them. Finite spacing changes the action, finite volume changes long-distance fluctuations, numerical stochastic time can bias an uncorrected sampler, and finite trajectory length limits statistical information. A continuum claim must specify which physical quantities are held fixed as bare parameters and lattice size change. We use a real scalar action with dimensionless lattice variables and restore length dimensions explicitly, while keeping simulation time separate from any physical time coordinate.
Restore dimensions and follow a line of matched physics
For a continuum Euclidean scalar action with a canonical gradient term, use units in which its statistical or quantum action scale has been absorbed so the action is dimensionless. In d dimensions the physical field has length dimension one minus half d. The lattice field rescaling then makes the nearest-neighbor difference term dimensionless. The quadratic and quartic bare coefficients contain powers of lattice spacing shown below. Their renormalized counterparts need not equal these tree-level bare expressions in an interacting theory.
A continuum limit requires the correlation length in lattice spacings to grow while the physical reference length stays fixed. Near an interacting critical point this generally means tuning the bare mass parameter toward a coupling-dependent critical value, not merely setting its classical value to zero. Set the scale using a specified observable and match appropriate renormalized couplings along the refinement sequence. If ratios of independent physical observables drift without control, the simulations may be following different physical theories instead of one theory with a smaller cutoff. Existence and nontriviality of the limiting theory are additional questions, not consequences of refining a grid.
Volume and cutoff limits ask different questions
With L sites in one direction, the physical box length is a times L. Refining at fixed physical box therefore increases L. This can define a continuum theory in a finite box, but it does not make the box infinite relative to the correlation length. A thermodynamic or infinite-volume result requires a separately controlled ratio of box size to physical correlation length. At a critical point where that length diverges, finite-size scaling replaces a simple assumption that every box is already much larger than the correlations.
For a nearest-neighbor free scalar, the lattice momentum kernel differs from the continuum squared momentum by terms quadratic in spacing at fixed low momentum. This gives a useful analytic cutoff benchmark, with the small parameter being momentum times spacing. An interacting observable can have logarithms, other correction exponents, or action-dependent leading powers, so an arbitrary straight line in spacing squared is not universally justified. At fixed physical volume the number of lattice variables grows as inverse spacing to the d, and critical slowing can add further cost.
Universality needs matched observables and a reproducible record
Compare distinct lattice discretizations only after matching their symmetries, dimensionality, relevant couplings, and physical scales. In an appropriate common scaling regime, critical exponents and suitably normalized dimensionless relations can agree while bare critical couplings and nonuniversal amplitudes differ. Binder crossing values also depend on shape and boundary conventions. Algorithmic dynamic exponents can differ even when both algorithms sample the same static universality class. Agreement of one exponent or a good-looking collapse is therefore supporting evidence rather than a complete identification.
Preserve the action and lattice geometry, update schedule, thermalization decision, measurement definitions, random generator and stream construction, code version, precision, and raw or adequately blocked data. A seed alone is insufficient when generators, parallel ordering, or analysis change. Record covariance handling, fit ranges, correction assumptions, and scale-setting correlations so the extrapolation can be repeated. Independent implementations and exact tiny-system benchmarks support this record. The planar matrix lab below illustrates one-axis cell-area bookkeeping only; it has no spacing-dependent field action or continuum solver.
Work it through
Refine while keeping both physical lengths fixed
A simulation has physical correlation length 2 fm, spacing 0.1 fm, and L = 200 sites along each direction. Refine to spacing 0.05 fm while retaining the physical box and correlation length. Find both lattice correlation lengths, the new L, and the variable-count factor in two dimensions.
Initially and the box length is . The box-to-correlation-length ratio is ten.
After halving a, matched physics requires and . The ratio L divided by lattice correlation length remains ten.
In two dimensions the number of sites grows by . This counts variables only; force costs, autocorrelation, and warmup can change as well. The sequence refines a fixed physical box and does not by itself take an infinite-volume limit.
The lattice correlation length grows from 20 to 40, L grows from 200 to 400, and two-dimensional site count increases fourfold while the physical box-to-correlation ratio stays ten.
The living notebook
Explore the model
Interactive lab · matrix
A cell-area change is separate from field dynamics
Reset and set diagonal scale , shear ; predict selected final-matrix determinant one. Repeat from Reset with both controls , ; predict determinant 0.5, a one-axis planar cell-area reduction. At Reset the drawn geometry is identity; Play or Step interpolates it, while determinant and singular-value readouts always describe the selected final matrix. This fixed triangular family is not a lattice discretization, field eigensolver, or continuum simulation.
- Area scale · det A
- 1.4
- Singular values · σmax, σmin
- 1.597, 0.876
- Condition number · κ₂(A)
- 1.823
The selected matrix is A = [[a,b],[0,1]]. Playback interpolates geometry between the identity and A; readouts always describe A. Its determinant is a and its singular values come from AᵀA. Positive a makes every supported matrix invertible. This two-dimensional family is not a general eigensolver or a learned network. 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 forRefining geometric cells changes area bookkeeping; matching field parameters, physical volume, and continuum observables requires additional physics.
Pause & reason
Check your understanding
Your turn
Try it yourself
A dimensionless matched-volume observable obeys exactly the toy ansatz . Measurements are 1.08 at and 1.02 at 0.1. Find Q₀ and c₂, and state what these two points cannot test.
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
Subtract the two equations before solving for the intercept; both measurements use the same physical volume and matched theory.
Compare with the solution
The difference is , so . Substituting gives . The two points determine this assumed two-parameter form exactly, but cannot test an omitted quartic term, a logarithm, mismatched physical scales, finite-volume error, or underestimated sampling uncertainty. Additional spacings and volume checks are needed before interpreting the intercept as a controlled physical prediction.
Keep the idea
Take it with you
- Spatial refinement must match physical scales and renormalized physics, not just increase lattice dimensions.
- Finite-volume continuum, infinite-volume, sampling, and stochastic-time limits control different effects.
- Universality comparisons and reproducible extrapolations require shared conventions and preserved statistical evidence.
Read further
Original sources and deeper explanations to follow at your own pace.
- Kari Rummukainen, Simulation Methods in Physics, Quantum field theory on the lattice, pp. 1–20: Euclidean measure, scalar discretization, and volume versus continuum limits
- Kari Rummukainen, Simulation Methods in Physics, Finite size scaling: corrections, universal quantities, and finite-volume critical analysis
Optional progress, saved in this browser when available. Revisit any lesson at any time.
Lesson not marked complete.