research30 min studyLesson 950 in Physics

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

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.

φx=a(d−2)/2ϕphys(x),r=a2m02,u=a4−dλ0,ξphys=a ξlat\varphi_x=a^{(d-2)/2}\phi_{\mathrm{phys}}(x),\qquad r=a^2m_0^2,\qquad u=a^{4-d}\lambda_0,\qquad \xi_{\mathrm{phys}}=a\,\xi_{\mathrm{lat}}
a is a length, phi_phys has length dimension (2−d)/2, m₀ has inverse-length units, and lambda₀ has length dimension d−4. varphi is the dimensionless field called phi in the sampling lesson. The parameter identifications refer to the stated canonical bare action with quartic normalization 1/4!.↔ Scroll wide equations horizontally.

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.

ℓphys=aL,p^2=4a2∑μ=1dsin⁡2(apμ2)=p2−a212∑μpμ4+O(a4p6)\ell_{\mathrm{phys}}=aL,\qquad \widehat p^2=\frac4{a^2}\sum_{\mu=1}^d\sin^2\left(\frac{ap_\mu}{2}\right)=p^2-\frac{a^2}{12}\sum_\mu p_\mu^4+O(a^4p^6)
p has inverse-length units, and the expansion assumes every |a p_mu| is small. The displayed kernel is a tree-level nearest-neighbor scalar result, not a universal interacting cutoff-error formula.↔ Scroll wide equations horizontally.

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.

Q(a,ℓphys)=Q0(ℓphys)+c2(aℓ∗)2+⋯Q(a,\ell_{\mathrm{phys}})=Q_0(\ell_{\mathrm{phys}})+c_2\left(\frac a{\ell_*}\right)^2+\cdots
Q is a dimensionless observable and ell_* a fixed physical reference length. This is a conditional cutoff ansatz for a matched finite-volume sequence with a justified leading quadratic error; Q₀ can still depend on volume.↔ Scroll wide equations horizontally.

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.

  1. Initially ξlat=2/0.1=20\xi_{\mathrm{lat}}=2/0.1=20 and the box length is 0.1×200=20 fm0.1\times200=20\,\mathrm{fm}. The box-to-correlation-length ratio is ten.

  2. After halving a, matched physics requires ξlat=2/0.05=40\xi_{\mathrm{lat}}=2/0.05=40 and L=20/0.05=400L=20/0.05=400. The ratio L divided by lattice correlation length remains ten.

  3. In two dimensions the number of sites grows by (400/200)2=4(400/200)^2=4. 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 result

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.

Explore the model

Interactive lab

Interactive lab · matrix

A cell-area change is separate from field dynamics

Reset and set diagonal scale a=1a=1, shear b=0b=0; predict selected final-matrix determinant one. Repeat from Reset with both controls a=0.5a=0.5, b=0b=0; 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.

A cell-area change is separate from field dynamicsThe 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. Current readings: Area scale · det A: 1.4; Singular values · σmax, σmin: 1.597, 0.876; Condition number · κ₂(A): 1.823. All parameters are adjustable with the labeled controls below.Matrix A with first row 1.4, 0.6 and second row zero, one; display interpolates between identity and A while final determinant remains 1.4.−4−2024−2−1012x · equal geometric scaleyA = [ 1.4 0.6 ; 0 1 ]Display blend: 0% · final values below
Interpolated geometryOriginal unit squareSelected final transform
Area scale · det A
1.4
Singular values · σmax, σmin
1.597, 0.876
Condition number · κ₂(A)
1.823
t = 0 s

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.

Check your understanding

Halving lattice spacing at fixed physical box length requires what change in the number L of sites along each direction?

Try it yourself

A dimensionless matched-volume observable obeys exactly the toy ansatz Q(a)=Q0+c2(a/ℓ∗)2Q(a)=Q_0+c_2(a/\ell_*)^2. Measurements are 1.08 at a/ℓ∗=0.2a/\ell_*=0.2 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 0.06=c2(0.04−0.01)0.06=c_2(0.04-0.01), so c2=2c_2=2. Substituting gives Q0=1.02−2(0.01)=1.00Q_0=1.02-2(0.01)=1.00. 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.

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.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.