research30 min studyLesson 946 in Physics

Finite-size critical-exponent fitting

Extract logarithmic size ratios while distinguishing designed power-law examples from controlled asymptotic exponent estimates.

On this page

What you’ll understand

  • Define susceptibility and magnetization observables with consistent units.
  • Extract exponent ratios and nu from size quotients.
  • Diagnose correction-to-scaling bias and correlated-fit uncertainty.

A straight line on a log-log plot is a beginning of a scaling analysis, not its conclusion. Critical exponents describe asymptotic behavior near a particular fixed point, while measured finite-lattice observables contain regular backgrounds, correction terms, and sampling errors. The fitted coupling may itself be uncertain. We use L as the number of sites along one direction of a fixed-shape d-dimensional box and keep all lattice observables dimensionless. Synthetic data below are constructed to make the algebra visible; their fitted powers are not asserted to be measured exponents of a physical universality class.

Choose the response and order parameter precisely

For the Ising source H appearing in the dimensionless action, differentiating the mean magnetization density gives the susceptibility N times the variance of m. If the source instead is the energy-valued physical field h, an extra inverse-temperature factor enters. At exact zero field in a finite symmetric ensemble, the signed mean magnetization vanishes, so absolute magnetization or its root-mean-square can serve as a nonzero finite-size amplitude. A susceptibility with absolute magnetization subtracted is another observable and must be labeled rather than silently substituted.

At a continuous transition in the ordinary scaling regime, susceptibility grows with exponent gamma over nu, while absolute magnetization shrinks with the magnetization exponent divided by nu. A dimensionless Binder slope grows with inverse nu. We write the magnetization exponent as beta mag to avoid confusing it with inverse temperature. These laws assume the correct critical coupling, fixed shape and boundaries, and appropriate control of irrelevant fields. First-order coexistence, upper-critical-dimensional logarithms, or dangerous irrelevant variables can require modified forms.

χH=N(⟨m2⟩−⟨m⟩2),χH(Kc,L)∼Lγ/ν,⟨∣m∣⟩∼L−βmag/ν,∂U∂t∣0∼L1/ν\chi_H=N(\langle m^2\rangle-\langle m\rangle^2),\qquad \chi_H(K_c,L)\sim L^{\gamma/\nu},\quad\langle|m|\rangle\sim L^{-\beta_{\mathrm{mag}}/\nu},\quad \left.\frac{\partial U}{\partial t}\right|_0\sim L^{1/\nu}
All displayed lattice observables and reduced coupling t are dimensionless. Response to h = k_BT H has additional inverse-energy units; beta_mag is a critical exponent, not the inverse temperature.↔ Scroll wide equations horizontally.

Ratios remove an amplitude, not every correction

A ratio between sizes L and sL cancels a leading amplitude. Taking its logarithm and dividing by logarithm of s gives an effective exponent. For an exact power law this equals the exponent at every size; with a correction proportional to L raised to minus omega, the quotient drifts. A regular additive background can produce a different drift, especially when the singular contribution is not yet dominant. A visually stable pair of adjacent ratios can result from cancellation between corrections.

Evaluate several sizes and vary the minimum size retained in the fit. Fit corrections only when the data can constrain their amplitudes and powers; otherwise an apparently flexible model can conceal the absence of an asymptotic regime. Reweighting to a pseudocritical crossing can reduce some tuning uncertainty but introduces a size-dependent coupling trajectory. Its scaling interpretation must match the chosen estimator. A slope exponent is sensitive to how the reduced coupling is defined, although a constant coupling rescaling changes its amplitude rather than its asymptotic power.

peff(L,s)=ln⁡[O(sL)/O(L)]ln⁡s=p+ln⁡[(1+bs−ωL−ω)/(1+bL−ω)]ln⁡sp_{\mathrm{eff}}(L,s)=\frac{\ln[O(sL)/O(L)]}{\ln s}=p+\frac{\ln[(1+bs^{-\omega}L^{-\omega})/(1+bL^{-\omega})]}{\ln s}
This expression assumes O(L) = A L^p(1 + b L^(−omega)) with positive measured O and no additional terms. It isolates one known correction form rather than guaranteeing that a real data set follows it.↔ Scroll wide equations horizontally.

Fit covariance and test physical consistency

Measurements at multiple couplings derived from the same trajectory are correlated. Bootstrap or jackknife resamples should repeat the full moment construction, reweighting, coupling selection, and exponent fit, using blocks that resolve the slowest relevant measurement mode. A poorly estimated covariance matrix can make a correlated fit unstable; blindly inverting it does not make the analysis more reliable. Independent replicas and an explicit covariance-treatment choice help separate statistical limitations from an inadequate scaling model.

Scaling relations provide consistency checks only under their hypotheses. For short-range systems below the upper critical dimension with ordinary hyperscaling, twice beta mag plus gamma equals d times nu. Boundary conditions, anisotropy, long-range forces, and finite imaginary-time extent in quantum models require further analysis. Agreement of one fitted ratio with a famous number cannot establish a universality class. Compare several dimensionless quantities and exponents, while preserving a record of discarded sizes and fit choices. The quadratic instrument tests one supplied power and performs no exponent inference.

2βmag+γ=dνunder ordinary hyperscaling2\beta_{\mathrm{mag}}+\gamma=d\nu\qquad\text{under ordinary hyperscaling}
This relation is a conditional scaling hypothesis with a specified regime, not a numerical identity for every finite system or for arbitrary mean-field and upper-critical-dimensional limits.↔ Scroll wide equations horizontally.

Recover designed exponent ratios

At the stipulated critical point of a synthetic three-dimensional scaling example, let susceptibilities at L = 8, 16, 32 be 16, 64, 256, while Binder slopes at L = 8 and 16 are 8 and 16. Extract gamma/nu and nu, then infer beta_mag/nu only if ordinary hyperscaling is assumed.

  1. Each susceptibility doubles its size argument while increasing by a factor four. The logarithmic quotient is ln⁡4/ln⁡2=2\ln4/\ln2=2, so the designed data have γ/ν=2\gamma/\nu=2.

  2. The slope ratio is two, giving 1/ν=ln⁡2/ln⁡2=11/\nu=\ln2/\ln2=1 and therefore nu equal to one. This uses the stipulated exact critical coupling and exact power laws.

  3. Under ordinary hyperscaling in three dimensions, 2βmag/ν+γ/ν=32\beta_{\mathrm{mag}}/\nu+\gamma/\nu=3. Hence βmag/ν=1/2\beta_{\mathrm{mag}}/\nu=1/2. These are algebraic properties of supplied synthetic data, not a numerical determination of the three-dimensional Ising exponents.

The result

The designed ratios give gamma/nu = 2 and nu = 1; conditional hyperscaling gives beta_mag/nu = 1/2, without identifying a physical universality class.

Explore the model

Interactive lab

Interactive lab · function

A supplied quadratic finite-size law

Reset and set coefficient a=1a=1, position x=1x=1; predict height one. Repeat from Reset with both controls a=1a=1, x=2x=2; predict height four. Let x represent a normalized size only for the supplied law O∝L2O\propto L^2. This parabolic function illustrates the quotient four at size ratio two; it has no lattice data, critical-coupling tuning, correction fit, or statistical uncertainty.

A supplied quadratic finite-size lawThis lab uses only y = ax². Playback traces the curve within the fixed integration interval from 0 to x; it does not estimate numerical error or change the chosen point. Current readings: y at chosen x: 0.608; Tangent slope · 2ax: 1.35; Signed area · ax³/3: 0.182. All parameters are adjustable with the labeled controls below.−2.4−1.201.22.4−9−4.504.59xy(x, ax²)
y = ax²TangentSigned integral
y at chosen x
0.608
Tangent slope · 2ax
1.35
Signed area · ax³/3
0.182
t = 0 s

This lab uses only y = ax². Playback traces the curve within the fixed integration interval from 0 to x; it does not estimate numerical error or change the chosen point. 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 forAn instrument can verify a declared power-law quotient while leaving the physical exponent-estimation problem entirely open.

Check your understanding

If a positive observable increases by a factor eight when L doubles and follows an exact power law, what is its exponent?

Try it yourself

Suppose the exact finite-size law is χ(L)=AL7/4(1+4/L)\chi(L)=A L^{7/4}(1+4/L) with positive A. Find the effective exponent from L = 8 and 16, and compare it with the asymptotic value.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

The correction factors are 3/2 and 5/4; their ratio multiplies the leading size quotient.

Compare with the solution

The susceptibility quotient is 27/4(5/6)2^{7/4}(5/6). Therefore peff=7/4+ln⁡(5/6)/ln⁡2≃1.48697p_{\mathrm{eff}}=7/4+\ln(5/6)/\ln2\simeq1.48697, below the asymptotic value 1.75. The correction is known exactly in this toy example, so no statistical uncertainty is involved. A two-size fit that omits it would return a biased effective power despite perfectly precise data; larger sizes or an adequately constrained correction model are necessary.

Take it with you

  • Observable definitions and coupling units precede any exponent fit.
  • Size ratios cancel amplitudes but retain correction and tuning effects.
  • Covariance, fit-range stability, and conditional scaling identities support a universality assessment together.

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.