Codex @codex  0

The correlation-function susceptibility sum rule and $G(r)=r^{-(D-2+\eta )}{f}_G(r/\xi )$ give ${\chi }_s\propto {\xi }^{2-\eta }$ when $f_G(0)$ is finite, its large-distance tail is integrable and $\eta <2$. With correlation-length critical exponent $\nu$, this gives the displayed relation for the magnetic-susceptibility critical exponent. Microscopic distances contribute a regular background rather than the divergent critical power.

First interpret the displayed critical correlation function as the connected correlation functionAbove the transition at zero field the subtraction is zero; in an ordered phase it is essential. For a dimensionless source coupled to the total magnetization $M=\int d^{d}xm(x)$, differentiating $Z$ twice gives $\chi =V^{-1}(⟨M^2⟩-⟨M{⟩}^2)$. Translation invariance then proves the correlation-function susceptibility sum ruleInsert the scaling form, use spherical coordinates and set $s=r/\xi$:For $\eta <2$, finite nonzero $g(0)$ and an integrable large-$s$ tail, the last integral tends to a finite nonzero constant. Microscopic separations add an analytic response background. Thus the singular magnetic susceptibility obeys ${\chi }_s\asymp {\xi }^{2-\eta }$, and at zero field $\xi \asymp |t{|}^{-\nu }$ gives the Fisher scaling relationA dimensional magnetic field would add a nonsingular inverse-temperature prefactor. Integrating the full unsubtracted ordered-phase correlation function instead would yield a volume-divergent term $V⟨m{⟩}^2$, which is not the intrinsic susceptibility appearing in this exponent identity.

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)$.

Write ${\beta }_T=(k_{B}T{)}^{-1}$ to distinguish inverse temperature from an order-parameter critical exponent. For the Ising model, define the full two-point correlation function $G(r)=⟨{\sigma }_0{\sigma }_r⟩$ and its connected correlation function $G_c(r)=G(r)-M^2$ in a translation-invariant phase. Away from a thermodynamic critical point, the large-distance envelope decays exponentially, possibly multiplied by a power of distance:This defines the correlation length in physical distance units when the limit exists. At a continuous phase transition, $\xi \to \infty$ and the critical connected correlation function typically has algebraic decay $G_c(r)\sim |r{|}^{-(D-2+\eta )}$, where $\eta$ is the anomalous dimension. Below the transition the full $G(r)$ tends to $M^2$, so the decay definition must use $G_c$.

Let $S={\sum }_r{\sigma }_r$ and $F=-{\beta }_T^{-1}\log Z$ be the extensive Helmholtz free energy. Holding the field-independent constant $C$ and other physical couplings fixed,Because the Hamiltonian contains $-hS$, ${\partial }_h\log Z={\beta }_T⟨S⟩$ and ${\partial }_h^2\log Z={\beta }_T^2(⟨S^2⟩-⟨S{⟩}^2)$. Translation invariance then proves the correlation-function susceptibility sum rule:In the ordered zero-field phase these statements refer to a selected pure thermodynamic phase; averaging equally over opposite magnetizations introduces the separate spin-mixture contribution to zero-field susceptibility.

A normalized blocking kernel $K_b({\sigma }^{\prime },\sigma )\ge 0$ satisfies ${\sum }_{{\sigma }^{\prime }}{K}_b({\sigma }^{\prime },\sigma )=1$. For example, assign one block spin to each block of $b^D$ microscopic spins using majority sign, with equal probabilities for a tie. Define the blocked Hamiltonian, retaining the entire generated operator space, byThe field-independent contribution $N{g}_E$ belongs to the identity operator. With $N^{\prime }=N{b}^{-D}$ andsumming over ${\sigma }^{\prime }$ proves $Z(u,C,N)=Z(u^{\prime },C^{\prime },N^{\prime })$ exactly. Iterating gives $N_p=N{b}^{-pD}$ and the same partition function with $(u_p,C_p,N_p)$. Observables at large scales are reproduced by appropriately blocked observables and sources. Truncating the generated interactions makes this an approximation; dropping the additive constant already destroys the equality of free energies even when normalized spin expectations remain correct.

The extensive Helmholtz free energy is unchanged: $F(u_0,C_0,N)=F(u_p,C_p,N_p)$. In contrast, its value per site changes because the number of sites changes. Define the reduced zero-constant free-energy density $f_0(u)={\beta }_{T}F(u,0,N)/N$ in the thermodynamic limit. The preceding equality yieldsand iteration proves the additive free-energy recursion under blocking:Equivalently $C_p=b^{pD}[C_0+{\sum }_{j=0}^{p-1}{b}^{-jD}{g}_E(u_j)]$. Thus the inhomogeneous source records the identity-operator contribution to renormalization-group free energy from eliminated short-distance degrees of freedom.

To express this recursion for the singular part, write $f_0=a+f_s$, where $a$ is the regular analytic background. Then the same displayed recursion holds with $f_s$ andThis is the analytic subtraction of an inhomogeneous renormalization recursion. A purely analytic source can often be absorbed into $a$, after which it does not control the leading singular powers. It is not legitimate simply to discard an analytic source in the full free energy, where it can dominate the regular background. Moreover an analytic monomial $c{t}^m{h}^n$ contributesIf $m{l}_t+n{l}_h=D$, the sum is $p$ and hence generates a logarithm when $p\sim -\log |t|/(l_t\log b)$: a renormalization-group free-energy resonance. Such a term cannot be discarded when describing that logarithmic singularity. Neglect of the inhomogeneous term in the leading power-law scaling therefore assumes that its regular part has been subtracted, that any remaining source only changes finite scaling amplitudes, and that there are no relevant resonances or marginal logarithms. Omitting an irrelevant operator further requires absence of a dangerously irrelevant coupling, whose vanishing limit can make a scaling function singular.

A renormalization-group fixed point obeys $R_b{u}_{\ast }=u_{\ast }$. Linearization gives $\delta u^{\prime }=A\delta u$; in eigenvector coordinates $w_i^{\prime }=b^{l_i}{w}_i$. Positive $l_i$ label relevant operators, negative $l_i$ label irrelevant operators, and zero $l_i$ labels a marginal operator requiring nonlinear analysis. The critical surface is the stable manifold obtained by tuning all relevant scaling fields to zero. A repulsive renormalization-group trajectory leaves the fixed point along a relevant direction as the observation length increases. Locally a typical zero-field section has $dw/d\ell =-\omega w$, $dt/d\ell =l_{t}t$, $\ell =\log b$, so the critical surface $t=0$ attracts along $w$, while trajectories with either sign of $t$ depart from it. A magnetic scaling field adds another relevant direction.

Under the stated homogeneous-scaling assumptions and after fixing irrelevant fields at their limiting values, two relevant scaling fields obey $t^{\prime }=b^{l_t}t$, $h^{\prime }=b^{l_h}h$. These fields are analytic coordinates near the renormalization-group fixed point, proportional to the physical reduced temperature and conjugate field to leading order. The singular free-energy density consequently obeysHere $F_s$ denotes the singular free energy per original site or unit volume, with regular dimensional factors absorbed; the extensive singular part is $N{F}_s$. Iterate and choose $b^p=|t{|}^{-1/l_t}$ to obtain the scaling hypothesis for critical phenomena:The sign labels the two sides of the transition: $f_{+}$ describes $t>0$, the disordered side, while $f_{-}$ describes $t<0$, the ordered side. The two functions need not coincide; $f_{-}$ has the zero-field one-sided derivative appropriate to spontaneous magnetization. The numbers $l_t,l_h$ are scaling eigenvalues in logarithmic form: the corresponding discrete RG multipliers are $b^{l_t},b^{l_h}$.

The correlation length obeys $\xi (t,h)=b\xi (b^{l_t}t,b^{l_h}h)$. The same scale choice gives $\nu =1/l_t$. Differentiating $F_s$ once and twice with respect to the conjugate field, and twice with respect to temperature, givesAt $t=0$, instead choose the scale $b^p=|h{|}^{-1/l_h}$. Then $F_s(0,h)\sim |h{|}^{D/l_h}$ and $M\sim \mathrm{sgn}(h)|h{|}^{D/l_h-1}$, so $\delta =l_h/(D-l_h)$. Thus differentiation and scale matching establishThe first is the Widom scaling relation and the second the hyperscaling relation. The hyperscaling relation requires precisely the regular homogeneous scaling used here. Above the ordinary upper critical dimension, the quartic coupling is dangerously irrelevant in the ordered phase and this derivation cannot discard it: the mean-field critical exponents $\nu =1/2$, $\alpha =0$ do not obey $\alpha =2-D\nu$ for $D>4$. At marginal dimensions logarithmic factors likewise require care beyond the displayed pure powers.

Use ${\beta }_{\mathrm{th}}=1/(k_{B}T)$ for inverse temperature to distinguish it from the order-parameter exponent. The decaying two-point correlation function needed for the response identity is connected:Away from the critical point in a short-range Ising model pure phase, the correlation length controls the exponential long-distance falloff, $G(r)\sim A(r)e^{-|r|/\xi }$ with a possible algebraic prefactor. The full two-spin expectation below the transition instead tends to $M^2$. The spin-mixture contribution to zero-field susceptibility explains why a symmetric mixture is unsuitable for defining a finite connected correlation length in the ordered phase.

Let $F=-{\beta }_{\mathrm{th}}^{-1}\log Z$ be the total Helmholtz free energy. Take magnetization and susceptibility per site, soDifferentiating the Boltzmann weights at fixed temperature givesTranslation invariance then yields the correlation-function susceptibility sum ruleTotal magnetization and total susceptibility multiply these expressions by $N$; per physical volume they have an additional factor $a^{-D}$ relative to the per-site convention. These normalizations do not change the critical powers.

For the subsequent scaling requests, let $F_s$ denote singular free energy per site; the fixed extensive normalization is suppressed. The labels $+$ and $-$ denote $t>0$ and $t<0$, with distinct scaling functions and amplitudes on the two sides. Below the transition a one-sided derivative at $h=0$ selects a spontaneous-magnetization branch. Writing $x=h/|t{|}^{\Delta }$, differentiation of the scaling hypothesis for critical phenomena givesThe field-independent $I_{\pm }$ makes no contribution to these field derivatives, although it can change heat-capacity amplitudes. Assuming nonzero leading amplitudes and the indicated pure-phase limits, comparison with the exponent definitions givesThis is the Rushbrooke scaling relation. A temperature-independent critical isotherm requires the large-$|x|$ magnetization scaling function to grow as $\mathrm{sgn}(x)|x{|}^{\beta /\Delta }$: otherwise $M_s=|t{|}^{\beta }{\mathcal{M}}_{\pm }(x)$ would retain a temperature power as $t\to 0$ at fixed small $h$. Therefore $1/\delta =\beta /\Delta$. Since the preceding two equalities also give $\Delta =\beta +\gamma$, the Widom scaling relation follows:For the hyperscaling relation, additionally assume that one correlation volume supplies an order-one singular free-energy contribution and that there is a single isotropic correlation length, with no dangerously irrelevant coupling spoiling this argument. Then $F_s\sim {\xi }^{-D}\sim |t{|}^{D\nu }$. Two temperature derivatives give the singular heat capacity, soHyperscaling is an additional physical assumption, not a consequence of the two-variable equation alone, and generally fails above the upper critical dimension.

Finally the correlation scaling function $f_G(x)$ falls exponentially, up to algebraic factors, for $x\to \infty$ in a massive short-range pure phase. At $x\to 0$ it approaches a finite nonzero value in the scaling regime. Approximating the lattice sum at distances much larger than $a$, and denoting the area of the unit sphere by $S_{D-1}$, givesIf $\eta <2$ and the integral has a finite nonzero limiting value, microscopic distances contribute only a regular background and the leading divergence is ${\xi }^{2-\eta }$. Thus the Fisher scaling relation isThese statements concern leading singular powers; logarithmic corrections and cancellations of amplitudes require separate treatment.

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.

Use ${\beta }_{\mathrm{th}}=1/(k_{B}T)$ for inverse temperature, to distinguish it from the order-parameter exponent. In a translation-invariant pure thermodynamic phase, the connected correlation function isThe subtraction is essential below the transition. Its exponential correlation length is defined by ${\xi }^{-1}=-{\lim }_{r\to \infty }{r}^{-1}\log |G(r)|$, with direction specified if the lattice is anisotropic. Away from criticality, $G(r)$ has a tail proportional to $e^{-r/\xi }$ times an algebraic prefactor. In the massive continuum scalar approximation that prefactor is $r^{-(D-1)/2}$; at criticality the exponential cutoff disappears and $G(r)\asymp r^{-(D-2+\eta )}$ instead. The large-distance limit here is $r\gg \xi$, as the PDF shows.

Let $F_{\mathrm{tot}}=-{\beta }_{\mathrm{th}}^{-1}\log Z$. With the energy convention $-h{\sum }_r{\sigma }_r$, the magnetization per site and magnetic susceptibility areWriting $S={\sum }_r{\sigma }_r$, direct differentiation of the partition function gives ${\partial }_h⟨S⟩={\beta }_{\mathrm{th}}(⟨S^2⟩-⟨S{⟩}^2)$. Thus the correlation-function susceptibility sum rule isIn physical continuum coordinates the lattice sum becomes $a^{-D}\int d^{D}rG(r)$, so omitting the site-density factor would change the normalization. A derivative with respect to the dimensionless source ${\beta }_{\mathrm{th}}h$ instead absorbs the explicit inverse-temperature factor.

An exact renormalization-group transformation with length factor $b>1$ can be defined by a nonnegative normalized blocking kernel $K_b({\sigma }^{\prime },\sigma )$ satisfying ${\sum }_{{\sigma }^{\prime }}{K}_b({\sigma }^{\prime },\sigma )=1$. A deterministic majority rule is one example; randomized tie-breaking preserves its normalization. Define the blocked energy, including its constant, bySumming over ${\sigma }^{\prime }$ proves $Z(u^{\prime },C^{\prime },N^{\prime })=Z(u,C,N)$. All generated operators must be retained for this identity to remain exact; a finite-coupling truncation is an approximation. After $p$ iterations, $N_p=N/b^{pD}$ and the spacing is $a_p=b^{p}a$. Long-distance physics is preserved because this is a marginalization of microscopic probabilities, with observables and sources transformed along with the energy. It integrates short-distance variables rather than changing their effects on retained observables. Physical lengths remain fixed even though the correlation length in units of the new spacing is reduced by $b^p$.

For clarity, now use dimensionless free energy per site $F(u,C)=-N^{-1}\log Z=f(u)+{\beta }_{\mathrm{th}}C$. Exact equality of the partition functions givesA constant microscopic energy shifts every blocked weight by the same factor, so $C^{\prime }=b^{D}C+c(u)$. Substitution at one step yields the additive free-energy recursion under blockingInduction now gives the requested relationThe function $g$ is the identity-operator contribution to renormalization-group free energy: it records eliminated-mode entropy, integration constants and normalization factors. A free-energy function satisfying this inhomogeneous equation need not be purely singular until its analytic background has been subtracted.

At a renormalization-group fixed point, $R_b{u}_{\ast }=u_{\ast }$. Linearize in scaling coordinates $v_i$ so $v_i^{\prime }=b^{{\lambda }_i}{v}_i$. Here ${\lambda }_i$ is a scaling exponent, not the discrete multiplier itself. A positive ${\lambda }_i$ gives a relevant operator, a negative one an irrelevant operator, and zero a marginal operator. The critical surface is the stable manifold: tune every relevant scaling field to zero and the irrelevant coordinates flow toward $u_{\ast }$. A repulsive renormalization-group trajectory leaves along a relevant direction; it can still have irrelevant coordinates contracting toward that manifold. The plot shows the $h=0$ section when temperature and field are relevant: $t=0$ is the critical line in this section, while an independent nonzero field is another direction away from criticality.

The inhomogeneous term cannot simply be discarded from the full free energy. If an analytic background $A(u)$ solves $A(u)=b^{-D}A(R_{b}u)+g(u)$, its subtraction leaves $F_s=f-A$ obeying homogeneous scaling. For example, in linear thermal and field coordinates, a source monomial $g_{mn}{t}^m{h}^n$ is removed byWhen the denominator vanishes, the source is resonant and can instead produce a logarithm. Thus the homogeneous singular equation is justified after analytic subtraction when there is no relevant resonance affecting the leading singularity, no marginal logarithm that has been omitted, and no singular dependence on a dangerously irrelevant coupling. Merely saying that $g$ is regular is insufficient to discard the whole weighted sum without this analysis.

Under these assumptions, with two relevant scaling fields $t,h$, choose $p$ so that $b^{p{\lambda }_t}|t|$ is of order one. The scaling hypothesis for critical phenomena then givesWriting ${\Phi }_{\pm }=f_{\pm }+I_{\pm }$ recovers the PDF's notation:The sign labels the two matching branches $t>0$ and $t<0$, which enter different phases when the flow leaves the critical neighborhood. The constants $I_{\pm }$ retain a branch-dependent identity contribution or a chosen split of the scaling function. They are determined by the same matching problem, not independently adjustable nonuniversal free-energy offsets. An arbitrary analytic offset from $C_0$ belongs to the subtracted background. The formula describes the singular part; adding that background recovers the full physical free energy.