Rate function for magnetization

The large-deviation rate function governing probabilities of atypical magnetization values under a Gibbs measure; its minimizers describe typical phases and its shape encodes coexistence and metastability.
Rate function for magnetization

Let σi{1,+1}\sigma_i \in \{-1,+1\} be spins in a volume Λ\Lambda with Λ=N|\Lambda|=N. The (intensive) magnetization is

mN  =  1NiΛσi. m_N \;=\; \frac{1}{N}\sum_{i\in\Lambda}\sigma_i.

Under an equilibrium Gibbs measure (see ), mNm_N often satisfies a large deviation principle (LDP) at speed NN:

P(mNm)    exp(NI(m)), \mathbb{P}(m_N \approx m)\;\asymp\;\exp\big(-N\, I(m)\big),

where II is the rate function. See and .

Connection to the pressure (log-partition density)

In the canonical ensemble with external field hh (see and ), exponential tilting of magnetization corresponds to shifting the field.

Indeed,

Eβ,h ⁣[etNmN]=ZN(β,h+t/β)ZN(β,h). \mathbb{E}_{\beta,h}\!\left[e^{t\,N m_N}\right] ={} \frac{Z_N(\beta,\,h+t/\beta)}{Z_N(\beta,\,h)}.

If the pressure

p(β,h)  =  limN1NlogZN(β,h) p(\beta,h) \;=\; \lim_{N\to\infty}\frac{1}{N}\log Z_N(\beta,h)

exists (see ), then the scaled cumulant generating function is

Λβ,h(t)=p(β,h+t/β)p(β,h). \Lambda_{\beta,h}(t)=p(\beta,h+t/\beta)-p(\beta,h).

When Λβ,h\Lambda_{\beta,h} is differentiable, the Gartner–Ellis mechanism yields a convex rate function given by the Legendre–Fenchel transform

Iβ,h(m)=suptR{tmΛβ,h(t)}. I_{\beta,h}(m)=\sup_{t\in\mathbb{R}}\Big\{t m-\Lambda_{\beta,h}(t)\Big\}.

Non-differentiability of p(β,h)p(\beta,h) in hh is a signature of phase coexistence and can produce flat pieces in Iβ,hI_{\beta,h}.

Interpreting the shape of I(m)I(m)

  • Single phase (no coexistence): I(m)I(m) has a unique minimizer at the typical magnetization.
  • Symmetry breaking / coexistence: at h=0h=0 below critical temperature in Ising-type systems, I(m)I(m) may have multiple minimizers (or a flat interval of minimizers at speed NN), reflecting competing phases; see and .
  • Metastability: local (non-global) minima of I(m)I(m) correspond to metastable magnetizations in a static large-deviation sense; see .

A key caution at coexistence: intermediate magnetizations may be realized by phase separation with probability costs scaling like surface area (not volume). In that case, the speed-NN rate function can develop a plateau (zero cost on an interval), while the sharper cost is governed by interfaces and surface tension; see .

Explicit example: Curie–Weiss (mean-field) magnetization rate function

For the Curie–Weiss model (see ), the magnetization obeys an LDP with an explicit rate function.

Let

I0(m)=1+m2log1+m2+1m2log1m2+log2for m[1,1], I_0(m)=\frac{1+m}{2}\log\frac{1+m}{2}+\frac{1-m}{2}\log\frac{1-m}{2}+\log 2 \quad\text{for } m\in[-1,1],

the Cramér rate function of i.i.d. ±1\pm1 spins (see ). Then for coupling JJ and field hh, define the mean-field “free-energy-like” function

ϕβ,h(m)=I0(m)βJ2m2βhm. \phi_{\beta,h}(m)= I_0(m) - \frac{\beta J}{2}m^2 - \beta h\, m.

The magnetization rate function can be written as

Iβ,h(m)=ϕβ,h(m)infm[1,1]ϕβ,h(m). I_{\beta,h}(m)=\phi_{\beta,h}(m)-\inf_{m'\in[-1,1]}\phi_{\beta,h}(m').

Its minimizers are exactly the stable equilibrium magnetizations, and its local minima encode mean-field metastability (compare and ).