advanced32 min studyLesson 340 in Physics

Finite-size scaling and data collapse

Use sample size as a controlled infrared cutoff and assess whether rescaled observations support a stated critical scaling hypothesis.

On this page

What you’ll understand

  • Relate finite-size rounding to the ratio of correlation length and sample size.
  • Construct scaled thermal and susceptibility coordinates consistently.
  • Evaluate Binder ratios and explain why visually attractive collapse is insufficient by itself.

Every laboratory sample and computer lattice is finite, so its thermodynamic functions are rounded rather than truly singular at positive temperature. That limitation becomes a tool when size is varied systematically. Finite-size scaling predicts how the rounding window, response peaks, and dimensionless ratios change as a continuous critical point is approached. A data collapse is a quantitative test of those shared predictions, not permission to tune arbitrary smooth curves until they overlap.

Size cuts off the critical length scale

Suppose a conventional short-range continuous transition has bulk correlation length ξ∼∣t∣−ν\xi\sim|t|^{-\nu}. A sample of linear size L cannot behave as an infinite collection of independent correlated regions once xi is comparable with L. The crossover condition gives a rounding window of order L−1/νL^{-1/\nu}. The relevant size is linear extent, not total site count; for a d-dimensional cubic sample, N grows as LdL^d.

The usual finite-size ansatz replaces a bulk divergence by a power of L times a scaling function of tL1/νtL^{1/\nu}. Boundary conditions, shape, and normalization conventions must be kept consistent across the size sequence. A susceptibility peak is a finite-size pseudocritical feature, and its location can shift with L. It should not be identified automatically with the infinite-system transition temperature or assigned a universal shift amplitude.

x=tL1/ν,χH(t,L)=Lγ/νX(x),⟨∣m∣⟩=L−β/νM(x),∣t∣round∼L−1/νx=tL^{1/\nu},\qquad \chi_H(t,L)=L^{\gamma/\nu}\mathcal X(x),\qquad \langle|m|\rangle=L^{-\beta/\nu}\mathcal M(x),\qquad |t|_{\mathrm{round}}\sim L^{-1/\nu}
L is measured in lattice-spacing units, t=(T−Tc)/Tc, and m=M/N is dimensionless. chiH=Var(M)/N is per-site response to dimensionless H=h/(kBT). These leading forms assume conventional scaling below the upper critical dimension, zero ordering field, a fixed shape and boundary condition, and negligible corrections.↔ Scroll wide equations horizontally.

Ratios can remove a leading unknown amplitude

For a scalar reversal-symmetric order parameter, define a Binder quantity from its second and fourth moments. The leading magnetization scale cancels in the ratio. With the convention used here, a narrow symmetric pair of peaks gives two thirds and a centred Gaussian gives zero. Some references multiply this definition by three halves to make the ordered limit one; numerical crossing values cannot be compared without translating that convention.

At an ordinary critical point, a dimensionless Binder ratio or xi divided by L can approach a size-independent value for fixed geometry. Crossings for successive size pairs therefore help estimate the critical temperature. Irrelevant operators generate corrections such as L−ωL^{-\omega}, causing crossing drift and imperfect overlap. The limiting ratio may be universal within a specified shape and boundary class, but it is not one number common to every geometry or every order-parameter dimension.

UL=1−⟨m4⟩3⟨m2⟩2,UL(t)=U(tL1/ν)+L−ωU1(tL1/ν)+⋯ ,ω>0U_L=1-\frac{\langle m^4\rangle}{3\langle m^2\rangle^2},\qquad U_L(t)=\mathcal U(tL^{1/\nu})+L^{-\omega}\mathcal U_1(tL^{1/\nu})+\cdots,\quad \omega>0
UL is dimensionless and uses a scalar order parameter with nonzero second moment. omega is a leading irrelevant correction exponent when a power correction applies. The Gaussian and two-peak limits refer to this convention; vector order parameters need different normalizations.↔ Scroll wide equations horizontally.

A credible collapse survives changes of analysis window

Choose the critical temperature and exponents from a physically justified model, then plot the scaled ordinate against the scaled thermal variable for all sizes. Test whether a single smooth scaling function describes statistically compatible data. Correlated Monte Carlo measurements need appropriate uncertainty estimates, and uncertainty in the critical temperature propagates into the horizontal coordinates. Agreement for one or two sizes is weak evidence because correction terms and flexible fitting functions can absorb systematic errors.

Repeat the analysis while removing the smallest sizes, narrowing the thermal window, and including justified correction terms. Compare other observables and dimensionless crossings rather than optimizing visual overlap alone. Above the upper critical dimension, dangerous irrelevant variables can modify the naive size powers; at a BKT transition logarithmic corrections and an essential correlation length require a different scaling treatment. First-order coexistence also has its own volume-dependent rounding, so it should not be forced into a conventional continuous-transition collapse.

Match scaled points in a two-dimensional Ising test

For an illustrative zero-field two-dimensional Ising scaling test use ν=1\nu=1 and γ=7/4\gamma=7/4. One point has L=16, t=0.020, and χH=100\chi_H=100. A second has L=32 and t=0.010. Find the second susceptibility required for leading scaling collapse and the common scaled coordinates.

  1. The horizontal variables are x1=0.020(16)=0.32x_1=0.020(16)=0.32 and x2=0.010(32)=0.32x_2=0.010(32)=0.32. Matching temperatures without this L rescaling would compare different critical regimes.

  2. The common vertical coordinate from the first point is Y1=χH/L7/4=100/167/4=100/128=0.78125Y_1=\chi_H/L^{7/4}=100/16^{7/4}=100/128=0.78125, because 167/4=27=12816^{7/4}=2^7=128.

  3. The second susceptibility must be larger by 27/4≃3.36362^{7/4}\simeq3.3636, giving about 336.36. Dividing it by 327/432^{7/4} gives the same 0.78125. This is a conditional prediction of the leading ansatz, not evidence from an independent data set.

The result

Leading collapse requires chiH≈336.36 at the second point; both scaled points are (0.32,0.78125), assuming corrections and backgrounds are negligible.

Explore the model

Interactive lab

Interactive lab · sampling

Separate sampling precision from physical size scaling

For each fresh run, Reset and set BOTH Success probability p=0.50 and Independent samples n. Compare n=25 and n=100, predicting sample-mean standard deviations 0.10 and 0.05. The instrument displays exact binomial distributions, not Monte Carlo resampling. Its n is a number of independent measurements, not an Ising lattice length L.

Separate sampling precision from physical size scalingExact binomial masses for the mean of n independent Bernoulli trials, compared with normal mass between adjacent half-count boundaries; normal tails outside the finite support remain outside the bars. Playback highlights outcomes, not random samples. The Berry–Esseen bound applies to standardized CDFs, not individual bars or continuity-corrected errors. At p = 0 or 1, normal comparison and that bound are disabled. Normal-CDF evaluation uses a numerical approximation with absolute error below 7.5×10⁻⁸ (at most 1.5×10⁻⁷ per interval, before rounding). Current readings: Mean sample proportion: 0.5; Sample-mean variance: 0.0125; Uniform CDF error bound: ≤ 0.1062. All parameters are adjustable with the labeled controls below.Exact distribution of the mean of 20 independent Bernoulli variables with success probability 0.5; playback highlights outcome 0.00.250.50.75100.0580.1150.1730.23Sample mean · k / nProbability massNormal bin masses are an approximationk = 0 · exact probability = 9.54e-7
Exact binomial massNormal interval massHighlighted outcome / mean
Mean sample proportion
0.5
Sample-mean variance
0.0125
Uniform CDF error bound
≤ 0.1062
t = 0 s

Exact binomial masses for the mean of n independent Bernoulli trials, compared with normal mass between adjacent half-count boundaries; normal tails outside the finite support remain outside the bars. Playback highlights outcomes, not random samples. The Berry–Esseen bound applies to standardized CDFs, not individual bars or continuity-corrected errors. At p = 0 or 1, normal comparison and that bound are disabled. Normal-CDF evaluation uses a numerical approximation with absolute error below 7.5×10⁻⁸ (at most 1.5×10⁻⁷ per interval, before rounding). 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 forMore independent measurements improve statistical precision. That operation differs from enlarging a correlated physical sample, and this lab does not perform a critical data collapse.

Check your understanding

Why can a visually good collapse of two small-system curves fail to establish a universality class?

Try it yourself

Two symmetric scalar magnetization distributions both have ⟨m2⟩=0.040\langle m^2\rangle=0.040. Distribution A is centred Gaussian, with fourth moment 0.0048. Distribution B consists of narrow peaks at m=±0.20, with fourth moment 0.0016. Calculate the Binder quantity for each and explain the distinction.

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

Give me a hint

Use the stated definition without the optional three-halves normalization. Identical second moments do not force identical fourth moments or distribution shapes.

Compare with the solution

The common denominator is 3(0.040)2=0.00483(0.040)^2=0.0048. For A, UA=1−0.0048/0.0048=0U_A=1-0.0048/0.0048=0. For B, UB=1−0.0016/0.0048=2/3U_B=1-0.0016/0.0048=2/3. The second moments alone therefore do not distinguish a broad Gaussian around zero from a symmetric pair with fixed nonzero magnitude; the fourth-moment ratio does. These limiting values diagnose distribution shape in the stated convention. A finite observed Binder value or one crossing does not alone prove an infinite-system phase transition: size dependence, boundary conditions, and corrections still need to be examined.

Take it with you

  • Finite-size scaling treats L as an infrared cutoff with specified geometry and boundaries.
  • Binder ratios remove leading magnetization scale but retain convention and correction dependence.
  • Data collapse must test a physical scaling hypothesis with uncertainty and stability checks.

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.