Past exam of the mathematics course of the University of Cambridge / 2006 / iii / Paper 51 (Statistical Field Theory)

An order parameter distinguishes phases by a macroscopic quantity that changes at the transition. For an Ising model ferromagnet it is the magnetization per Ising spin, $M=N^{-1}{\sum }_n⟨{\sigma }_n⟩$. At zero magnetic field the statistical Hamiltonian is invariant under spin inversion symmetry. The disordered phase preserves this symmetry, whereas a selected ordered phase has $M\ne 0$ and a spin-reversed partner with $-M$. At finite volume with symmetric boundary conditions the expected value is zero even below the transition; spontaneous magnetization means taking the thermodynamic limit before sending a symmetry-breaking field to zero.

Landau-Ginzburg theory introduces a coarse-grained local order parameter $\varphi (x)$ and a symmetry-constrained free-energy functional. For a short-range system with scalar spin inversion symmetry, a useful expansion isThe coefficients vary smoothly with microscopic control parameters at the level of the Landau approximation; normally $r=A(T-T_c)$ with $A>0$. The gradient term penalizes spatial variation. The magnetic field is conjugate to the order parameter, and odd powers at $h=0$ are excluded by spin inversion symmetry. With $u>0$, the sextic term may be omitted locally; with $u<0$, a stabilizing positive $v$ is essential. The equilibrium uniform order parameter minimizes the local potential $V(M)$, givingThis equation of state follows from minimization, rather than being imposed as an independent assumption. Vector or tensor order parameters and their allowed invariants describe other symmetry classes.

Use the connected correlation function $G(x)=⟨\varphi (x)\varphi (0)⟩-⟨\varphi {⟩}^2$, so that the constant ordered contribution is removed. Away from a thermodynamic critical point its long-distance decay defines the correlation length: $-{\lim }_{R\to \infty }{R}^{-1}\log |G(R)|={\xi }^{-1}$ when this limit exists. In a massive rotation-invariant phase the exponential decay is typically multiplied by an algebraic factor.

Expand the Landau-Ginzburg theory free energy about a stable uniform minimum $M$, writing $\varphi =M+\delta \varphi$. To quadratic order,integration of the Gaussian field theory gives the Fourier-space connected correlation function $\tilde{G}(q)=k_{B}T/(c+\kappa q^2)$. Thus the Ornstein--Zernike correlation function hasFor an ordinary quartic transition, $M=0$ above the transition and $M^2=-r/u$ below, so $c=r$ and $c=-2r$, respectively. The correlation length diverges as $|T-T_c{|}^{-1/2}$ on both sides, with different amplitudes. At criticality the mass vanishes and the Gaussian correlation function is proportional to $R^{-(D-2)}$ for $D>2$. Interactions can change this power to $R^{-(D-2+\eta )}$.

At a continuous phase transition the equilibrium order parameter tends continuously to its critical value, there is no latent heat, and the correlation length diverges. For $u>0$ at zero field, minimizing $V(M)=r{M}^2/2+u{M}^4/4$ gives $M=0$ when $r\ge 0$ and $M=\pm \sqrt{-r/u}$ when $r<0$. These minima join continuously at $r=0$. The minimized singular free-energy density is zero above and $-r^2/(4u)$ below, so its first thermal derivative is continuous, while its second derivative has a finite jump in mean-field approximation.

A first-order phase transition instead has competing minima whose free energies cross while their order parameters remain distinct. The first derivatives of the equilibrium free energy generally jump; along a temperature-driven crossing this gives latent heat. For the symmetric sextic potential with $u<0,v>0$, let $y=M^2>0$. At phase coexistence, the stationary points and equality of the two free energies requireEliminating $r$ gives $-u{y}^2/4-v{y}^3/3=0$, henceAt coexistence $V(M)=v{M}^2(M^2+3u/(4v){)}^2/6$, proving that all three degenerate minima are global. The equilibrium order parameter jumps from zero to one of the two nonzero values even though $r_{\mathrm{coex}}>0$. The disordered minimum loses local stability at $r=0$, and the nonzero stationary points first appear at $r=u^2/(4v)$. These spinodal points delimit metastability and are not the equilibrium transition line. Cubic terms, when symmetry permits them, provide another common route to a first-order transition.

Universality means that systems with different microscopic interactions can have the same asymptotic critical behaviour. Under renormalization-group flow, irrelevant microscopic differences shrink, and systems whose critical flows approach the same renormalization-group fixed point belong to one universality class. Important determinants are spatial dimension, order-parameter symmetry and number of components, the range of interactions and any additional relevant constraints. Microscopic Bravais lattice shape alone does not normally change the short-range scalar universality class.

Critical exponents, properly normalized scaling functions and certain dimensionless amplitude ratios are universal within a universality class. For example, the universal specific-heat amplitude ratio compares the two singular amplitudes above and below the transition, canceling their common normalization factors. Scaling functions require choices of thermal, field and length metric factors before comparisons between microscopic systems make sense. The critical temperature, critical field, Bravais lattice spacing, individual correlation-length and magnetic susceptibility amplitudes, and analytic free-energy backgrounds are not generally universal. Thus two models can share critical exponents without sharing their transition temperature or the absolute size of their response.

Let $t=(T-T_c)/T_c$, and suppose $r=r_{t}t$ with $r_t>0$ while $u>0$ and $\kappa >0$ remain finite at the transition. Critical exponents describe the leading singular powers: $M\sim (-t{)}^{\beta }$ below the transition, $\chi \sim |t{|}^{-\gamma }$, $\xi \sim |t{|}^{-\nu }$, $C_{V,s}\sim |t{|}^{-\alpha }$, $M(t=0)\sim \mathrm{sgn}(h)|h{|}^{1/\delta }$, and $G(R,t=0)\sim R^{-(D-2+\eta )}$.

The quartic equation of state $h=rM+u{M}^3$ immediately gives $M=\sqrt{-r/u}$ at zero field below the transition, so $\beta =1/2$. Differentiating with respect to the conjugate field gives ${\chi }^{-1}=r+3u{M}^2$, which is $r$ above and $2|r|$ below; therefore $\gamma =1$. At $r=0$, $M=(h/u{)}^{1/3}$ with its real sign, giving $\delta =3$. The quadratic correlation function has $\xi \propto |r{|}^{-1/2}$ and critical Ornstein--Zernike correlation function proportional to $q^{-2}$, giving $\nu =1/2$ and $\eta =0$.

Finally the minimized potential below the transition is $V_s=-r^2/(4u)$, while it is zero above. Its second temperature derivative is a finite step, so the mean-field heat-capacity critical exponent is $\alpha =0$. In particular,These are derived mean-field critical exponents. In the fluctuation-dominated critical region of a phase transition, one instead derives exponents from the scaling eigenvalues of the renormalization-group fixed point, as in the next question. A zero heat-capacity power exponent can also accompany logarithms in other theories, so specifying the power alone does not completely describe that singularity.

A tricritical point occurs when the continuous and first-order transition lines meet. In the symmetric sextic Landau free energy this requires tuning two coefficients, $r=u=0$, while $v>0$ stabilizes the potential. For $u>0$, the zero-field transition lies on $r=0$ and is continuous. For $u<0$, the coexistence curve is $r=3u^2/(16v)$. The order-parameter jump satisfies $M^2=-3u/(4v)$ and shrinks to zero as the tricritical point is approached. The labelled phase diagram below uses $(u,r)$ as local coordinates; two independent physical controls, such as temperature and composition, can provide these coordinates.

Along the tricritical tuning $u=0$, the equation of state is $h=rM+v{M}^5$. At zero field the ordered solution is $M=(-r/v{)}^{1/4}$, and the critical isotherm has $M=\mathrm{sgn}(h)(|h|/v{)}^{1/5}$. The inverse magnetic susceptibility is $r$ above and $4|r|$ below. The minimized singular potential is $-|r{|}^{3/2}/(3\sqrt{v})$ below, so differentiating twice thermally gives a $|t{|}^{-1/2}$ singularity. These calculations giveThe correlation-length and correlation-function exponents follow from the same quadratic gradient term. A generic path with a nonzero positive quartic coefficient instead has ordinary critical behaviour; observing tricritical powers requires the extra tuning. The sextic coupling is marginal in three dimensions, giving tricritical upper critical dimension three and logarithmic corrections there.

The mean-field approximation replaces each site's fluctuating surroundings by a self-consistent average. On a cubic lattice each spin has coordination number of a lattice $z=2D$, even though summing only over positive basis directions counts each bond once. Write ${\sigma }_n=M+\delta {\sigma }_n$ and drop products of deviations. Each Ising bond becomes ${\sigma }_n{\sigma }_m\simeq M{\sigma }_n+M{\sigma }_m-M^2$, so the approximate statistical Hamiltonian isThe constant corrects counting each bond twice. The resulting independent-spin partition function is $Z_{\mathrm{MF}}=e^{-{\beta }_{T}NzJ{M}^2/2}[2\cosh ({\beta }_{T}zJM){]}^N$, where ${\beta }_T=1/(k_{B}T)$. The one-spin probabilities giveEquating this average with the assumed $M$ givesHere $J>0$ is assumed; a ferromagnetic order parameter is not the appropriate choice for an antiferromagnetic coupling. Put $A=T_c/T$. For $A\le 1$, $\tanh (AM)<AM\le M$ at any $M>0$, so only $M=0$ is possible. For $A>1$, the function $\tanh (AM)-M$ has positive slope at zero, is strictly concave on $M>0$ and is negative at $M=1$. It therefore has exactly one positive zero and a spin-reversed negative zero. Expanding at the bifurcation givesThe transition is continuous, with order-parameter critical exponent $\beta =1/2$. Stability can also be checked from the Bragg-Williams free energy of the Ising model per Ising spin:Its derivative is $-zJM+k_{B}T\mathrm{arctanh}M$, giving the same self-consistency equation. Near zero,The quartic coefficient is positive at $T_c$, the stable minima approach zero continuously, and the minimized free energy has no first-derivative jump. Thus mean-field approximation predicts a second-order transition at $T_c=2DJ/k_B$. This is the mean-field Ising critical bifurcation.

For short-range interactions, large dimension gives many neighbours and suppresses the relative fluctuations of their average. More quantitatively, Gaussian field theory fluctuations averaged over a correlation volume have variance of order ${\xi }^{2-D}$, while the squared ordinary mean-field order parameter is of order $|t|\sim {\xi }^{-2}$. Their ratio is proportional to ${\xi }^{4-D}$. The Ginzburg criterion therefore makes neglect of critical fluctuations asymptotically consistent above four dimensions, inconsistent below four, and marginal at four, where logarithms occur. Large dimension is a route to validity; sufficiently long-range or infinite-range interactions can also produce mean-field behaviour. In low-dimensional short-range systems even the existence and location of a transition can differ from this approximation.

Use ${\beta }_T=1/(k_{B}T)$ for inverse temperature, to distinguish it from the order-parameter critical exponent $\beta$. In a translation-invariant pure thermodynamic phase, the connected two-point correlation function isFor finite correlation length, at $|r|\gg \xi$ it decays exponentially, possibly with an algebraic prefactor; ${\xi }^{-1}$ is the asymptotic exponential decay rate. Below an ordered transition the pure phase is selected by taking the thermodynamic limit before $h\to 0^{+}$ or $h\to 0^{-}$. A zero-field symmetric mixture of the two ordered phases has a nondecaying connected term and does not describe the usual pure-phase magnetic susceptibility.

Let $\mathcal{F}=-k_{B}T\log Z$ denote extensive physical free energy and $S={\sum }_r{\sigma }_r$. Since the field contributes $-hS$ to the statistical Hamiltonian, differentiation givesExpanding the variance and using translation invariance proves the correlation-function susceptibility sum rule:The sum includes $r=0$. If the source is instead defined as the dimensionless quantity ${\beta }_{T}h$, the corresponding derivative omits the displayed inverse-temperature factor.

A real-space renormalization group step of length factor $b>1$ uses a local blocking kernel $B_b({\sigma }^{\prime },\sigma )\ge 0$ satisfying ${\sum }_{{\sigma }^{\prime }}{B}_b({\sigma }^{\prime },\sigma )=1$. For example, a deterministic majority-spin rule assigns one coarse Ising spin to each block; tied blocks can be resolved with equal probabilities. Define the blocked statistical Hamiltonian byThe spacing becomes $a^{\prime }=ba$. All generated interactions, including multispin terms and the field-independent constant, must be retained for an exact transformation. Summing over ${\sigma }^{\prime }$ uses normalization to prove $Z(u^{\prime },C^{\prime },N^{\prime })=Z(u,C,N)$. Iteration gives $N_p=b^{-pD}N$ and $Z(u_p,C_p,N_p)=Z(u_0,C_0,N)$.

The coarse probability law is the exact marginal law of these block variables. Indeed, for any coarse observable $A$, its expected value is the microscopic expected value of ${\sum }_{{\sigma }^{\prime }}{B}_b({\sigma }^{\prime },\sigma )A({\sigma }^{\prime })$. A blocking rule that retains the long-wavelength order-parameter mode, with source and observable transformations tracked, therefore retains long-distance response and correlation functions. Their field normalization may change, while their physical length scales and critical singularities do not. Partition-function equality alone would not justify this assertion for an arbitrary kernel that discards all slow variables. Truncating the generated operator family is likewise an approximation.

For the free-energy recursion use dimensionless free energy per site,in the thermodynamic limit. Exact partition-function invariance gives $F(u_0,C_0)=b^{-pD}F(u_p,C_p)$. If one step sends $C$ to $C^{\prime }=b^{D}C+c(u)$, the coefficient $c(u)$ records the identity-operator contribution from eliminated fluctuations and normalization. ConsequentlyRepeated substitution derives the additive free-energy recursion under blocking:In physical energy units, divide this convention by ${\beta }_T$. The additive term is essential for absolute free energies even when it does not affect the leading critical powers.

A renormalization-group fixed point satisfies $R_b(u_{\ast })=u_{\ast }$. Linearize in scaling coordinates $v_i$ so that $v_i^{\prime }=b^{y_i}{v}_i+O(v^2)$. A relevant operator has $y_i>0$, an irrelevant operator has $y_i<0$, and a marginal operator has $y_i=0$ with behaviour decided by nonlinear terms. Thus the eigenvalues of the discrete linear map are $b^{y_i}$; the sign criterion applies to the generator exponents $y_i$, not to the discrete multipliers. The critical surface is the stable manifold on which every relevant scaling field has been tuned away. Irrelevant coordinates flow to the renormalization-group fixed point along it. A repulsive renormalization-group trajectory leaves the renormalization-group fixed point along an unstable relevant direction, ultimately entering a noncritical phase or another scaling regime.

To separate the singular free-energy density from its regular background, suppose the eliminated-mode contribution $g$ is analytic near the renormalization-group fixed point. Look for an analytic background $a(u)$ satisfying $a(u)=b^{-D}a(R_{b}u)+g(u)$. Then $F_s=f-a$ satisfies $F_s(u)=b^{-D}{F}_s(R_{b}u)$. In linear thermal and field coordinates, an analytic monomial $g_{mn}{t}^m{h}^n$ has background coefficientprovided the denominator is nonzero. This explicitly explains when analytic subtraction of an inhomogeneous renormalization recursion works. The printed inhomogeneous equation may describe the free-energy function before this subtraction, or a singular part defined only up to analytic terms. The inhomogeneous sum need not be small compared with the full free energy: an analytic constant can dominate a vanishing singular contribution. It can be omitted in deriving leading singular powers after regular subtraction, provided no renormalization-group free-energy resonance, marginal logarithm or dangerously irrelevant coupling obstructs the homogeneous critical asymptotics. At a renormalization-group free-energy resonance $m{\lambda }_t+n{\lambda }_h=D$, additive logarithms can remain, so simply dropping the term would miss them.

Suppose there are exactly two relevant scaling fields, thermal $t$ and magnetic $h$, with positive exponents ${\lambda }_t,{\lambda }_h$, and the remaining fields approach a regular limit. Physical reduced temperature and magnetic field agree with these coordinates to leading order, up to nonuniversal metric factors. The homogeneous scaling hypothesis for critical phenomena isas an asymptotic scaling equation. Choose $\ell =|t{|}^{-1/{\lambda }_t}$ to make the thermal argument $+1$ or $-1$. Defining $f_{\pm }(z)=F_s(\pm 1,z)$ givesThe two functions describe respectively $t>0$ and $t<0$: approaching the transition from disordered and ordered sides need not give the same function. In particular the ordered zero-field function has one-sided derivatives associated with spontaneous magnetization. Metric factors multiply this formula in different microscopic normalizations.

The correlation length in Bravais lattice units obeys $\xi (t,h)=\ell \xi ({\ell }^{{\lambda }_t}t,{\ell }^{{\lambda }_h}h)$. At $h=0$, the same choice of $\ell$ gives $\nu =1/{\lambda }_t$. Two thermal derivatives of $F_s(t,0)$ give $2-\alpha =D/{\lambda }_t$. One and two field derivatives give, respectively, $\beta =(D-{\lambda }_h)/{\lambda }_t$ and $\gamma =(2{\lambda }_h-D)/{\lambda }_t$, assuming the relevant branch amplitudes are nonzero. At $t=0$, choose $\ell =|h{|}^{-1/{\lambda }_h}$; then $M\propto \mathrm{sgn}(h)|h{|}^{(D-{\lambda }_h)/{\lambda }_h}$, so $\delta ={\lambda }_h/(D-{\lambda }_h)$. HenceSimple algebra now gives the Rushbrooke scaling relation $\alpha +2\beta +\gamma =2$ and the hyperscaling relation $\alpha =2-D\nu$. Also $\beta \delta ={\lambda }_h/{\lambda }_t=\beta +\gamma$, giving the Widom scaling relation. These statements refer to leading singular contributions with the homogeneous two-field scaling assumptions. Above the upper critical dimension, the quartic coupling is dangerously irrelevant because it is still needed to stabilize the ordered phase; ordinary mean-field critical exponents then satisfy Rushbrooke and Widom but generally not hyperscaling. Marginal logarithms also require qualifications at the upper critical dimension.

For the connected-correlation scaling form $G(r)=|r{|}^{-(D-2+\eta )}{f}_G(|r|/\xi )$, the function $f_G(x)$ tends to a finite nonzero value as $x\to 0$ and decays exponentially at large $x$, up to powers of $x$, in a massive short-range phase. The scaling ansatz alone does not fix that prefactor. If the usual massive Gaussian field correlation tail $G(r)\propto r^{-(D-1)/2}{e}^{-r/\xi }$ applies, matching the powers gives $f_G(x)\sim B{x}^{(D-3)/2+\eta }{e}^{-x}$, with a correlation-length-dependent amplitude in $G$ fixed by matching the critical scaling form.

In physical coordinates the lattice sum has site-density factor $a^{-D}$. Its long-distance singular part therefore giveswhere ${\Omega }_{D-1}$ is the area of the unit sphere in $D$ dimensions. For $\eta <2$, finite $f_G(0)$, an integrable exponential tail and a nonzero integral, the last factor tends to a finite nonzero constant as $\xi \to \infty$. Microscopic distances contribute a regular background. Thus ${\chi }_s\propto {\xi }^{2-\eta }\propto |t{|}^{-(2-\eta )\nu }$, proving the Fisher scaling relationEquivalently the field has scaling dimension $(D-2+\eta )/2$, so its conjugate field has ${\lambda }_h=(D+2-\eta )/2$, which gives the same identity from the exponent formulas above.

A partition function sums the Boltzmann weights of all states: $Z={\sum }_s{e}^{-{\beta }_{T}H(s)}$ for discrete states, or the corresponding integral for continuous fields. It normalizes equilibrium probabilities. Its physical Helmholtz free energy is $\mathcal{F}=-k_{B}T\log Z$, and its free-energy density is $\mathcal{F}/V$. With a momentum ultraviolet cutoff ${\Lambda }_0$, a scalar-field example isThe normalization of the functional measure and field-independent energy terms must be tracked when discussing absolute free energy.

Choose a smaller cutoff $\Lambda$ and split the field into retained modes ${\varphi }_{<}$ with $|q|<\Lambda$ and eliminated modes ${\varphi }_{>}$ with $\Lambda <|q|<{\Lambda }_0$. Define the Wilsonian coarse-grained statistical Hamiltonian bySubstitution into the partition function proves $Z=\int \mathcal{D}{\varphi }_{<}{e}^{-{\beta }_T{H}_{\Lambda }[{\varphi }_{<}]}$ for every intermediate cutoff. A source retained in this definition similarly preserves long-wavelength connected correlation functions and response. Thus the coefficients in a Hamiltonian density must run with $\Lambda$: the dependence induced by eliminated modes compensates the dependence of the remaining field integral. Holding the original bare coefficients fixed while changing the cutoff would generally change the system. Exact integration generates all interactions allowed by symmetry, including higher powers, derivative terms and an additive constant; a local derivative expansion or a truncation to a few couplings is an extra approximation.

Cutoff independence makes arbitrary intermediate coarse-graining scales unphysical. It allows calculations at the scale of the observable, avoids repeatedly resolving irrelevant short-distance details and exposes the renormalization-group fixed points responsible for universality. A microscopic cutoff, such as a Bravais lattice spacing, can still label a physical microscopic model; taking it to infinity as a continuum limit requires tuning appropriate bare couplings. This is distinct from lowering an intermediate Wilsonian cutoff in a fixed model.

To interpret the zero-cutoff limit, define $M=V^{-1}\int \varphi (x)d^{D}x$, the spatially uniform mode of the order parameter. Keep this mode fixed and integrate the others, rather than integrating over $M$ too. Write $\varphi =M+\eta$, with $\int \eta =0$, and define the constrained order-parameter free energy byIn a finite periodic box, once $\Lambda$ is below the smallest nonzero momentum, all nonconstant modes have been eliminated. The remaining effective energy is $H_{\Lambda }[M]=V{U}_{\Lambda }(M)$, so its defining integral is exactly $Z_M$. Therefore, with consistent normalization,When $\mathcal{H}(\Lambda ,M)$ denotes this uniform effective Hamiltonian density, this is the stated identification of $F(M)$ with its zero-cutoff limit. It requires a free-energy density constrained at fixed $M$, not the unconstrained extensive free energy; all integrated fluctuations and additive constants are included, the zero mode in field theory is retained, and a density for the uniform configuration is well defined. The finite-volume argument avoids an interchange of the thermodynamic and cutoff limits without justification. On a homogeneous pure-phase branch one may use a local potential description, while a thermodynamic limit allowing phase separation can convexify the constrained potential. A bare nonconvex Landau polynomial is accordingly a homogeneous-branch approximation, not a claim that the exact equilibrium Legendre potential stays nonconvex.

In a field $h$, the remaining integral is $Z(h)=\int dM{e}^{-{\beta }_{T}V[F_c(M)-hM]}$, up to the consistently chosen zero-mode measure. The thermodynamic saddle then gives $f(h)={\inf }_M\{F_c(M)-hM\}$ and selects the equilibrium magnetization. Landau theory assumes that the relevant uniform potential can be approximated by a regular analytic expansion in $M$, with symmetry-constrained coefficients, and minimizes it. At tree level the integrated fluctuations are neglected, recovering the familiar bare polynomial. With fluctuations retained, the coefficients run and can become singular near criticality; the zero-cutoff identification of constrained free energy alone does not justify mean-field critical exponents.

For an ordinary short-range scalar thermodynamic critical point, the leading interaction is quartic. Absorb inverse temperature and the gradient coefficient into field normalization and writeScale $x=b{x}^{\prime }$ and $\varphi (x)=b^{-(D-2)/2}{\varphi }^{\prime }(x^{\prime })$ to leave the gradient term invariant. The quartic term changes to $g{b}^{4-D}\int {\varphi }^{\prime 4}/4!$, so its engineering scaling exponent is $4-D$. It is irrelevant at the Gaussian fixed point for $D>4$, relevant for $D<4$, and marginal for $D=4$. HenceThe loop expansion makes the breakdown concrete. A quartic-vertex correction is proportional to $g^2{I}_2(r)$, whereFor $0<D<4$, rescaling $q=\sqrt{r}k$ gives $I_2(r)\sim C_D{r}^{(D-4)/2}$, diverging as $r\to 0$. At $D=4$, radial integration giveswhich is logarithmically divergent. Above four dimensions its infrared limit is finite; the leading cutoff-dependent constant can be absorbed into a coupling renormalization. The effective expansion parameter is thus of order $g{\xi }^{4-D}$ below four, or $g\log (\Lambda \xi )$ at four, since $\xi \sim r^{-1/2}$. It eventually becomes large for any fixed bare interaction in the lower-dimensional critical limit, defeating the uncorrected mean-field expansion. The Ginzburg criterion gives the same threshold by comparing correlation-volume fluctuations with the squared mean order parameter.

Below four dimensions an interacting renormalization-group fixed point can replace the Gaussian fixed point and change critical exponents, provided a continuous transition exists. At four dimensions the interaction is marginal and produces logarithmic corrections: the leading power exponents can remain their mean-field values, so “breakdown” there means that pure mean-field power laws omit those corrections, rather than necessarily that every power exponent changes. Above four the stabilizing quartic interaction is dangerously irrelevant for the ordered phase, explaining why mean-field thermodynamic exponents need not obey naive hyperscaling with the spatial dimension.