advanced31 min studyLesson 36 in Physics

Fluctuations and equilibrium response

Relate measurable equilibrium variances to heat capacity and magnetic susceptibility, while separating static identities from dynamical response.

On this page

What you’ll understand

  • Derive the static response of an observable to its conjugate field.
  • Compute magnetization and susceptibility for independent two-state moments.
  • Identify the equilibrium, linearity, and finite-size conditions behind fluctuation-response reasoning.

Noise is not always a nuisance to average away. In equilibrium, the size of spontaneous fluctuations can reveal how strongly a system responds to a small external change. Energy noise determines heat capacity; number noise determines number susceptibility. A magnetic two-state model lets us derive the corresponding identity exactly and see where a linear approximation stops working. We will keep static equilibrium response separate from the additional assumptions required to predict the time dependence after a perturbation.

Couple a field to the measured quantity

Write a classical Hamiltonian as Hh=H0−hMH_h=H_0-hM, where the generalized field hh has the energy units needed to make hMhM an energy. Assume that the observable MM itself has no explicit hh dependence. Differentiating its canonical expectation changes both the Boltzmann numerator and its normalizing partition function. The two terms combine into a variance multiplied by β\beta. Thus the response is nonnegative at positive temperature for a stable canonical equilibrium with this conjugate-field convention.

The same simple formula holds for a quantum observable commuting with the Hamiltonian and entering through the same linear coupling. Noncommuting quantum observables require an imaginary-time correlation expression instead of blindly substituting an equal-time variance. Likewise, an observable with explicit field dependence contributes an extra derivative term. Keeping these restrictions visible prevents a compact classical identity from being mistaken for a universal formula for every optical or quantum response.

Hh=H0−hM,χM=(∂⟨M⟩∂h)T=β(⟨M2⟩−⟨M⟩2)H_h=H_0-hM,\qquad \chi_M=\left(\frac{\partial\langle M\rangle}{\partial h}\right)_T=\beta\bigl(\langle M^2\rangle-\langle M\rangle^2\bigr)
This equal-time variance identity assumes classical equilibrium, or the stated commuting quantum case, a field-independent M, and a Hamiltonian linear in its conjugate field. The nonnegative-response statement assumes T > 0.↔ Scroll wide equations horizontally.

Independent moments expose the scale

For NN independent variables si=±1s_i=\pm1, define M=∑isiM=\sum_i s_i and energy −hM-hM. A single moment has partition factor 2cosh⁡(βh)2\cosh(\beta h), mean tanh⁡(βh)\tanh(\beta h), and variance sech⁡2(βh)\operatorname{sech}^2(\beta h). Independence makes means and variances additive. At zero field the mean magnetization vanishes, but the magnetization variance is NN. A zero mean therefore does not imply that every microscopic configuration is unmagnetized.

At small field, the total mean is approximately NβhN\beta h. This linear response requires ∣βh∣≪1|\beta h|\ll1; at larger field the exact hyperbolic tangent saturates instead of growing without bound. Distinguish the extensive total susceptibility from susceptibility per moment. A physical magnetic moment would introduce additional factors converting a laboratory magnetic field into the energy parameter hh. Our dimensionless spins deliberately keep that unit conversion separate from the statistical derivation.

⟨M⟩=Ntanh⁡(βh),Var⁡(M)=Nsech⁡2(βh),χM=Nβsech⁡2(βh)\langle M\rangle=N\tanh(\beta h),\qquad \operatorname{Var}(M)=N\operatorname{sech}^2(\beta h),\qquad \chi_M=N\beta\operatorname{sech}^2(\beta h)
Spins are dimensionless and h is an energy per spin. At zero field χM=N/(kBT), while the susceptibility per spin is χM/N.↔ Scroll wide equations horizontally.

Correlations, dynamics, and the observation window

For interacting moments, the total variance contains a sum of connected correlations between different sites, not only individual-site variances. Long-range correlations near a critical point can strongly amplify susceptibility and disrupt ordinary inverse-square-root relative-noise scaling. In a finite system the correlations and variance remain finite. A fixed-magnetization ensemble would instead constrain the total variance, illustrating why one must match the ensemble and field contact before interpreting an observed fluctuation.

The static identity compares equilibria after a small parameter change; it gives no relaxation time by itself. Dynamical fluctuation-dissipation relations connect frequency-dependent response to time correlations under specified equilibrium and microscopic assumptions. A system driven far from equilibrium, aging over the measurement, or observed before it relaxes need not obey the equilibrium formula. Even in equilibrium, a short experimental trace can miss slow fluctuations and underestimate a susceptibility. The measurement window and equilibration protocol belong beside the algebra, rather than being treated as irrelevant details.

Predict a weak-field response from zero-field noise

There are 100 independent dimensionless moments at thermal energy kBT=2k_BT=2 in a chosen energy unit. Find their zero-field variance and susceptibility, then compare the linear and exact magnetization at h=0.1h=0.1 in that unit.

  1. At zero field, each spin has variance one and the independent sum has Var⁡(M)=100\operatorname{Var}(M)=100. Also β=1/2\beta=1/2 in inverse energy units.

  2. The susceptibility is χM=βVar⁡(M)=50\chi_M=\beta\operatorname{Var}(M)=50 per energy unit. Linear response predicts ⟨M⟩≈χMh=5\langle M\rangle\approx\chi_Mh=5.

  3. The exact result is 100tanh⁡(0.05)≈4.99584100\tanh(0.05)\approx4.99584. The small difference is consistent with the first nonlinear correction being cubic in βh\beta h; the parameter controlling the approximation is 0.05, not the raw field number alone.

The result

A static noise measurement predicts a small-field equilibrium response. It does not predict how quickly that magnetization is reached.

Explore the model

Interactive lab

Interactive lab · probability

Convert excitation-count noise into spin noise

Interpret each success as an up spin among n independent moments, so M = 2K − n for binomial success count K. The displayed count variance np(1−p) becomes spin variance 4np(1−p). Set p to one half to reproduce zero-field fluctuations.

Convert excitation-count noise into spin noiseThese are exact binomial probabilities for independent trials with one fixed success probability. Playback tours the distribution; it does not resample outcomes. Current readings: Mean · np: 5; Variance · np(1−p): 2.5; Total probability: 1. All parameters are adjustable with the labeled controls below.03581000.070.150.220.3Successes kP(X = k)k = 0 · P = 9.77e-4
Exact probabilityHighlighted outcomeMean np
Mean · np
5
Variance · np(1−p)
2.5
Total probability
1
t = 0 s

These are exact binomial probabilities for independent trials with one fixed success probability. Playback tours the distribution; it does not resample outcomes. 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 forThe lab calculates an exact independent-spin distribution. It does not model interacting correlations, susceptibility dynamics, or a field-to-probability conversion unless that relation is supplied from the equilibrium formulas.

Check your understanding

Which extra statement is needed before replacing a quantum static susceptibility by β times an equal-time variance?

Try it yourself

At zero field, double the number of independent spins from NN to 2N2N while holding temperature fixed. How do total magnetization variance, total susceptibility, and the standard deviation of magnetization per spin change?

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 addition of variances and divide the total magnetization by the new system size before comparing relative noise.

Compare with the solution

The variance doubles from NN to 2N2N, and total susceptibility doubles from βN\beta N to 2βN2\beta N. The standard deviation of M/NM/N is initially 1/N1/\sqrt N and becomes 1/2N1/\sqrt{2N}, a factor 1/21/\sqrt2 smaller. This reduction assumes independence; long-range correlations can change the size dependence.

Take it with you

  • Static equilibrium response follows from differentiating a normalized probability distribution.
  • Total susceptibility can be extensive even while fluctuations per particle shrink.
  • Quantum noncommutation, correlated criticality, and nonequilibrium measurement require additional care.

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.