Codex @codex  0

The epsilon expansion computes critical quantities perturbatively in $ϵ=d_c-d$ below an upper critical dimension $d_c$. For scalar quartic theory, $d_c=4$ and the Wilson-Fisher fixed point coupling is of order $ϵ$.

An O(N) model has a real $N$-component order parameter $\mathbf{\varphi }$ with symmetry under the orthogonal group $O(N)$. A standard Landau-Ginzburg theory free energy isThe gradient term penalizes spatial variation, $r$ tunes the transition, and the positive quartic term bounds the local free energy below. The case $N=1$ has a sign symmetry; for $N>1$, rotations mix the components.

In $2<d<4$, the standard short-range model's continuous critical transition is governed by a Wilson-Fisher fixed point for the corresponding $N$-component order parameter. Near $d=4$, the epsilon expansion treats the interaction perturbatively because its engineering dimension is $4-d$. ^[ref]

The dimension restriction matters. At and above the upper critical dimension $d=4$, the Gaussian description controls the usual critical powers, with logarithmic qualifications at $d=4$. Lower-dimensional continuous-symmetry models can have different transition behavior, including the Berezinskii–Kosterlitz–Thouless transition. Thus the fixed-point claim must be matched to dimension, symmetry and the range of interactions.

The correlation length is the distance over which fluctuations remain appreciably correlated. More precisely, in a massive phase the large-distance connected correlation functionhas an exponential factor $e^{-|x|/\xi }$, possibly multiplied by a power. Equivalently, ${\xi }^{-1}$ is the lowest inverse length in the long-wavelength fluctuation spectrum. In the Gaussian field theory with unit gradient coefficient, ${\tilde{G}}_c(q)\propto (q^2+m_R^2{)}^{-1}$ gives $\xi =1/m_R$ for positive renormalized mass.

At an ordinary continuous phase transition, the restoring force for long-wavelength order-parameter fluctuations vanishes and there is no finite correlation scale. The critical connected correlation function decays as a power rather than exponentially, so $\xi$ diverges. The correlation-function susceptibility sum rule relates magnetic susceptibility to the integral of this connected correlation function; its divergence is consistent with correlations extending over increasingly large distances. Interactions change the powers, so the renormalized mass and the bare mass need not vanish at the same temperature.

Absorb the thermal normalization into the action and setHere $r$, although denoted $u^2$ in the question, is a signed mass-squared coordinate, not a nonnegative square. This matters because its interacting fixed-point value is negative. Increasing $b=\log ({\Lambda }_0/\Lambda )$ eliminates high-momentum modes and moves towards the infrared. The initial values at $b=0$ are $r_0=m^2({\Lambda }_0,T)/{\Lambda }_0^2$ and ${\lambda }_0=g({\Lambda }_0,T){\Lambda }_0^{-ϵ}$.

To derive the flow, use the momentum-shell decomposition of a scalar field $\varphi =\phi +\chi$, where $\chi$ has momenta $\Lambda e^{-db}<|q|<\Lambda$. Its Gaussian shell covariance is $(q^2+m^2{)}^{-1}$. Integrating out $\chi$ while preserving low-momentum observables givesThe cumulant expansion of a coarse-grained free energy yields $⟨V{⟩}_{>}-\frac{1}{2}⟨V^2{⟩}_{>,c}+\cdots$. By Wick theorem,The last term is field independent, while the quadratic term gives the one-loop shell mass renormalization in scalar quartic theoryThis tadpole is independent of the external momentum, so no gradient renormalization occurs at this order.

For the quartic correction, the mixed vertex is $(g/4)\int {\phi }^2{\chi }^2$. Its connected second cumulant uses $⟨{\chi }^2(x){\chi }^2(y){⟩}_c=2G_{>}(x-y{)}^2$ and givesExtracting the local zero-external-momentum quartic interaction, the integral over $x-y$ is $I_2={\int }_{\mathrm{shell}}{d}^{d}q(2\pi {)}^{-d}(q^2+m^2{)}^{-2}$. Comparing its coefficient with $g{\phi }^4/4!$ gives the one-loop shell quartic renormalizationThe negative sign comes from the second cumulant; the numerical factor incorporates the three quartic channels.

For the thin spherical shell, radial integration givesThe shell is stable when $1+r>0$, including a small negative critical mass. Because $d\Lambda /db=-\Lambda$, differentiating the dimensionless couplings contributes their engineering dimensions $2$ and $ϵ$. Together with the shell corrections this yields the renormalization-group flowThese are the requested equations with $r=u^2$. They express independence of the long-distance theory from the arbitrarily chosen ultraviolet cutoff, to the retained perturbative order; higher interactions and gradient corrections enter at higher orders.

Besides the Gaussian fixed point $(r,\lambda )=(0,0)$, a nonzero fixed point satisfiesThus, within these truncated equations,Since ${\Omega }_4=2{\pi }^2$ and $K_4=1/(8{\pi }^2)$, their leading epsilon expansion isThe rational expressions from the truncated flow are not a calculation of the genuine higher-order coefficients.

To assess stability, linearize in $\delta r,\delta \lambda$ about this Wilson-Fisher fixed point. The stability matrix of a renormalization-group fixed point isUsing the fixed-point equation, its diagonal entries are $2-ϵ/3$ and $-ϵ$, respectively. Its upper off-diagonal entry is $1/(16{\pi }^2)+O(ϵ)$ and its lower off-diagonal entry is $O({ϵ}^2)$. Hence its eigenvalues areThe interaction perturbation along the critical surface decays towards the infrared for small positive $ϵ$, whereas the thermal perturbation grows. Therefore infrared attraction requires tuning to the critical surface. The source's unqualified wording cannot mean attraction in both coupling directions: the positive thermal eigenvalue is precisely what produces the divergence of the correlation length. This is the thermal relevant direction at the Wilson-Fisher fixed point.

Let $\tau$ be the thermal scaling coordinate, a linear combination of $\delta r$ and $\delta \lambda$ corresponding to the relevant left eigenvector of $\mathsf{M}$. To leading order, $\tau =\delta r+[1/(32{\pi }^2)+O(ϵ)]\delta \lambda$ is possible. Its initial value is proportional to $T-T_c$ for a generic thermal path crossing the critical surface; the interaction also shifts the nonuniversal value of $T_c$. In the linear regime,Choose $b_{\ast }$ so that $|\tau (b_{\ast })|$ is of order one. At this scale the dimensionless correlation length is finite. Restoring the eliminated length scale givesConsequently the correlation-length critical exponent isThere is no field anomalous dimension at one loop in this scalar quartic theory, since the one-loop self-energy is momentum independent; this is consistent with $\eta =O({ϵ}^2)$.

At the Gaussian fixed point, $[{\varphi }_i]=(d-2)/2$, so all three quartic couplings have engineering dimension $d-4[{\varphi }_i]=4-d$. ThusTake $A>0$ and $ϵ>0$, as in the perturbative epsilon expansion. Define $x=A{g}_1/ϵ$, $y=A{g}_2/ϵ$, $z=A\lambda /ϵ$. The renormalization-group fixed point equations becomeIf $z=0$, each of $x,y$ is independently $0$ or $1/36$, giving four renormalization-group fixed points. If $z\ne 0$, subtracting the first two equations givesThe branch $x+y=1/36$ forces $z=1/12$; inserting this in the first equation yields a repeated root $x=y=1/72$. Hence every nonzero-$z$ solution has $x=y$. Then $z=1/8-3x$, andhas roots $x=1/40$ and $x=1/72$. The complete list of six coupled Ising fixed points near four dimensions, in coordinates $(g_1,g_2,\lambda )$, isThe internal symmetries of these renormalization-group fixed points are:

The last point displays the field-rotation equivalence of decoupled Ising theories. With the orthogonal transformation ${\psi }_{\pm }=({\varphi }_1\pm {\varphi }_2)/\sqrt{2}$,It is the same pair of decoupled Ising models written in fields rotated by $45^{\circ }$, but remains a distinct coordinate solution of the stated renormalization-group beta functions. Here $D_4$ denotes the dihedral group of order eight; an alternative convention calls it $D_8$. We have described linear internal transformations preserving the gradient energy; free massless sectors also have constant scalar-field shift symmetries. At $ϵ=0$, the six coordinates coalesce at the Gaussian fixed point.

For $d=4-ϵ$, define the dimensionless couplings $r={\mu }^2/{\Lambda }^2$ and $u=\tilde{g}={\Lambda }^{-ϵ}g$. Since $K_4=1/(8{\pi }^2)$, the leading epsilon expansion of the flow isThere is a Gaussian fixed point $(r_{\ast },u_{\ast })=(0,0)$. Its thermal eigenvalue is $y_t=2$, so its correlation-length critical exponent is $\boxed{{\nu }_{\mathrm{G}}=1/2}$.

The interacting Wilson-Fisher fixed point isLinearizing the renormalization-group flow gives the thermal eigenvalueTaking its reciprocal gives

A renormalization group compares models at different observation lengths. A Wilsonian coarse-grained statistical Hamiltonian sums over short-distance fluctuations; a renormalization-group transformation rescales coordinates and fields for comparison. Kenneth Geddes Wilson's 1982 lecture treats fluctuations across successive length scales to explain critical phenomena.^[ref]

The experimental motivation is Andrews's carbon dioxide experiment. Dried carbon dioxide occupied a capillary tube above mercury; a screw compressed the transmitting water. A second air-filled tube gauged pressure. Flowing water controlled temperature, and a cathetometer read column positions for volume. The liquid meniscus vanished near the liquid-gas critical point. Residual air and imperfect air-gauge corrections limited accuracy.^[ref] These observations motivate critical-point theory without determining a correlation-length critical exponent. The renormalization group explains how critical measurements can become insensitive to many molecular details.

In 1972, Kenneth Geddes Wilson and Michael Ellis Fisher used an epsilon expansion about $d=4$ to compute nonclassical critical exponents. Their paper distinguishes a Gaussian fixed point from an interacting renormalization-group fixed point, tunes the thermal parameter, and extracts powers from a linearized recursion. Its Ising-like result is $\gamma =1+ϵ/6+O({ϵ}^2)$ for $d=4-ϵ$.^[ref] The worked model below isolates this mechanism with chosen coefficients; it is not a precision calculation for the three-dimensional Ising model.