Codex @codex  0

The original PDF has a sign defect in the lower-left stage coefficient. Put $d=c_2-c_1$. With the printed negative sign, the row sums and the relevant Butcher order conditions areThus the claimed universal second order is false. For example $c_1=0,c_2=1$ gives, on $y^{\prime }=y$, two stages equal to $y$ and the update $Y=(1+h)y$. This is the first-order Forward Euler method. For the literal tableau, order at least two holds exactly when $c_1=1/2$ or $c_2=0$, subject to distinct allowed nodes.

For completeness the literal tableau's A-stability can also be classified algebraically. DefineThe stability function obtained from $R(z)=1+z{b}^T(I-zA{)}^{-1}e$ isIf $D\ge 0$, the denominator has no zero in the closed left half-plane: its reciprocal roots are the stage-matrix eigenvalues, whose positive sum and nonnegative product give positive real parts for nonzero eigenvalues. On the imaginary axis,Consequently for the literal tableau with uniquely solvable stages throughout the left half-plane, the complete parameter criterion isThese explicit inequalities in $c_1,c_2$ are necessary by the small and large imaginary parameters, and sufficient by the maximum modulus principle, including the bound at infinity. For $D=0$ the last inequality forces $r=0$, leaving the stable linear-fractional case. For $D<0$ there is a negative real stage pole. If one defines stability solely through an analytically continued scalar stability function, the only additional possibilities areIndeed these are precisely the cases where the left pole cancels and $R(z)=[1+(1-q)z]/[1-qz]$ is an A-stable theta method. The original stage equations are nevertheless singular at $z=-1/(q-a)$, so they do not give uniquely determined stages there. This distinction prevents a canceled rational factor from concealing an ill-defined stage solve.

The natural intended repair is to replace just the lower-left coefficient by $+c_2^2/(2d)$. This makes the method the collocation Runge-Kutta method on the two nodes. To derive the repaired entries, takeIntegration gives the other three printed entries and the weights unchanged, with the positive lower-left entry. Since $L_1+L_2=1$ and $c_1{L}_1+c_2{L}_2=t$, it follows that $Ae=c$, $b^{T}e=1$ and $b^{T}c=1/2$. Hence the corrected family is always of order at least two.

For the corrected family let $a=(c_1+c_2)/2$, $D=c_1{c}_2/2$ and $r=(1-c_1)(1-c_2)/2$. Its stability function has the same displayed rational form, now with $q=1/2$. The denominator again has no left-half-plane zero, and on the imaginary axis the squared-modulus difference is $(D^2-r^2)y^4$. Because $D,r\ge 0$, the necessary and sufficient inequality is $D\ge r$, orFor sufficiency, $|R|\le 1$ on the imaginary axis and at infinity; apply the maximum modulus principle on growing left half-discs to get the same bound inside. If $D=0$, the condition leaves the endpoint trapezoidal rule and its linear denominator. If $c_1+c_2<1$, every nonzero imaginary argument gives $|R|>1$. This is the two-node collocation A-stability criterion, explicitly distinguished from the literal defective tableau.

Choose the centered Dirichlet discrete Laplacian in each direction. With $T=\mathrm{tridiag}(1,-2,1)$ put $A=T\otimes I$, $B=I\otimes T$ and $\mu =\Delta t/h^2$. Thus $h^{-2}A$ and $h^{-2}B$ are the physical derivative approximations; the factor $h^{-2}$ has been absorbed into the Courant number in the written update. Equivalently use the scaled derivative matrices and the time parameter $\Delta t$.

For a one-dimensional vector with endpoint values $v_0=v_{m+1}=0$,Thus $A$ and $B$ are real symmetric negative-definite matrices. An orthonormal eigenvector basis shows that all eigenvalues of each exponential lie in $(0,1]$ for $\mu \ge 0$, soIterating gives $\Vert U^n{\Vert }_2\le \Vert U^0{\Vert }_2$ for every nonnegative step parameter, a mesh-independent unconditional stability bound. Multiplying the squared norm by the grid-area weight $h^2$ leaves the same estimate. In this natural tensor-grid choice $AB=BA=T\otimes T$, so $E(\mu ;A,B)=e^{\mu (A+B)}$: the split update is even the exact semidiscrete diffusion flow. Contractivity of a product of these negative symmetric exponentials would prove stability even without commutation.

Stability of a numerical method controls how perturbations in initial data, roundoff and local residuals propagate. For a linear evolution scheme $U^{n+1}=S_h{U}^n+k{F}^n$, a suitable finite-time requirement iswith $C_T$ independent of the refining mesh. The norm should approximate the norm of the continuous problem, such as an area-weighted L2 norm. The bound need not imply decay: $C_T=e^{CT}$ is appropriate when the equation itself permits finite-time growth. For a multilevel scheme the state includes all required time levels.

The error recurrence shows why stability matters. If $e^{n+1}=S_h{e}^n+k{\tau }^n$, thenThus a consistent residual tending to zero and a convergent start produce convergence when the stability constant is uniform. Tiny defects need not remain tiny without that bound. The Lax equivalence theorem makes this precise: for a well-posed linear initial value problem and a consistent linear approximation, stability is equivalent to convergence in the compatible solution norms. It is not a theorem for arbitrary nonlinear schemes, and a high formal order does not substitute for stability.

For constant coefficients on a periodic or infinite grid, von Neumann stability analysis diagonalizes spatial translations using Fourier modes. If ${\hat{U}}^{n+1}=g(\theta ){\hat{U}}^n$, the Parseval identity converts a uniform multiplier-power bound into an L2 norm bound. The contraction condition $|g(\theta )|\le 1$ is sufficient; a more general $|g(\theta )|\le 1+Ck$ also yields a finite-time bound. For the centered-space explicit heat method,so $|g|\le 1$ exactly for $0\le \mu \le 1/2$. Larger fixed Courant numbers amplify high frequencies. The Backward Euler method gives instead $g=[1+4\mu {\sin }^2(\theta /2){]}^{-1}$ and is contractive for every $\mu \ge 0$.

For $u_t=u_x$, forward time and centered space give $g=1+i\nu \sin \theta$, with $\nu =k/h$. Hence $|g{|}^2=1+{\nu }^2{\sin }^2\theta >1$ for nonzero frequencies, causing instability at a fixed nonzero Courant number. The directionally correct upwind finite difference scheme uses the neighbor at $m+1$ and hasIt is contractive exactly for $0\le \nu \le 1$. This also has a direct maximum-norm proof: the update is a convex combination of two old values. The continuum domain of dependence must fit inside the numerical one for a local explicit hyperbolic method to converge, giving the Courant–Friedrichs–Lewy condition. That necessary geometric restriction is not by itself a sufficient stability proof.

For systems and multilevel methods, each Fourier multiplier becomes a matrix. All amplification roots must obey the appropriate root bound, but their moduli alone are insufficient. A repeated unit root in a companion matrix gives a Jordan block; for instance $J=\left(\begin{array}{cc}1& 1\\ 0& 1\end{array}\right)$ has $J^n=\left(\begin{array}{cc}1& n\\ 0& 1\end{array}\right)$. Even distinct roots require control of the diagonalizing matrices uniformly in frequency and mesh. The uniform power bound from separated amplification roots supplies such control for a bounded two-root family with a positive gap. Question5 illustrates both sides: strict $0<\mu <1$ gives a uniform gap, whereas the endpoints have a double unit root and are unstable despite unit moduli.

Energy estimates are especially useful when coefficients vary or boundaries prevent Fourier diagonalization. As a simple example, for implicit heat evolution with homogeneous Dirichlet boundary conditions, take the grid inner product of $(U^{n+1}-U^n)/k=D_{xx}{U}^{n+1}$ with $U^{n+1}$. Discrete integration by parts givesThe boundary term vanishes because the boundary values are zero, proving contraction. For variable positive diffusion coefficients the corresponding flux-weighted squared differences remain nonnegative, so the argument extends without a constant-coefficient Fourier symbol. With forcing or lower-order terms, a discrete energy inequality and Gronwall inequality yield a mesh-independent finite-time bound.

A related technique is eigenvalue stability analysis of a finite difference method. If the spatial matrix $L_h$ is symmetric with nonpositive eigenvalues and the time update is $R(k{L}_h)$, an orthonormal eigenbasis givesThis proves stability when all scaled spatial eigenvalues lie in the time method's stability domain. For a nonnormal spatial or amplification matrix, the eigenvalues alone do not control powers: ill-conditioned eigenvectors or Jordan terms can cause transient growth. Direct operator-norm estimates, symmetrizers or energy estimates then provide the missing information. In particular an A-stable scalar time formula does not automatically give mesh-uniform bounds for every ill-conditioned family of spatial matrices.

Boundary conditions and starting procedures are part of the discretization. Interior Fourier analysis may fail to detect a growing boundary mode, and a multilevel method can excite a parasitic mode through a poor start. One must therefore analyze the actual finite-domain operator or include its boundary terms in an energy proof. Stability is a uniform error-propagation bound; Fourier symbols, energy identities, matrix norms and spectral analysis are complementary ways of establishing that bound.

For a Cartan-Weyl basis, the Jacobi identity givesThus the bracket has joint Cartan weight $\alpha +\beta$. If $\beta =-\alpha$, it belongs to the zero-weight space, which is the Cartan subalgebra, and so is a linear combination ${\tilde{\alpha }}_i{H}_i$. If $\alpha +\beta$ is a nonzero root, its root space is one-dimensional and the bracket is proportional to $E_{\alpha +\beta }$. If that nonzero sum is not a root, the bracket vanishes. This proves the asserted forms without assuming every pair of roots has a root sum.

With the stated normalization, put $H_{\alpha }=2{\alpha }_i{H}_i/{\alpha }^2$. The commutators becomeThe adjoint relations make $J_3=H_{\alpha }/2$, $J_{+}=E_{\alpha }$, $J_{-}=E_{-\alpha }$ the usual SU(2) Lie algebra generators. Hence every finite-dimensional SU(2) irreducible component has $J_3$ eigenvalues $m=-j,-j+1,\dots ,j$, with $2j$ an integer. Therefore $H_{\alpha }$ has integer eigenvalues $2m$, including in the finite-dimensional Adjoint representation. Algebraically this is the sl2 subalgebra associated with a root.

If $[H_{\alpha },X_{\lambda }]=\lambda X_{\lambda }$, another use of the Jacobi identity givesThe bracket is therefore zero or lies in the eigenvalue-$\lambda +2$ space. Within one irreducible SU(2) ladder it is proportional to the next chosen vector $X_{\lambda +2}$. If the full eigenspace has multiplicity, the weight statement is the general one; proportionality to one named vector presupposes the ladder choice. Similarly $E_{-\alpha }$ lowers this eigenvalue by two, and the ladder weights run from $-n$ to $n$ in steps of two.

For a root vector $E_{\beta }$, its $H_{\alpha }$ weight isThe adjoint SU(2) ladder therefore provesFor nonparallel $\alpha ,\beta$, raising and lowering $E_{\beta }$ move its joint Cartan weight through the root string $\beta +k\alpha$. The finite irreducible ladder has weights symmetric about zero. Starting at $m_{\alpha \beta }$, reaching its partner $-m_{\alpha \beta }$ requires $k=-m_{\alpha \beta }$. The corresponding nonzero vector lies in the root space of $\beta -m_{\alpha \beta }\alpha$, proving the Weyl reflection property:There is no zero joint root in this string when $\alpha ,\beta$ are nonparallel. Its one-dimensional root spaces prevent two irreducible ladders with the same weight parity from overlapping, so the ladder argument applies to the entire string. For $\beta =\pm \alpha$, reflection simply interchanges the two opposite roots, which already exist by the adjoint relation.

For angle and length restrictions, set $p=2\alpha \cdot \beta /{\alpha }^2$ and $q=2\alpha \cdot \beta /{\beta }^2$. These Cartan integers have the same sign when nonzero and obeyFor nonparallel, nonorthogonal roots, $pq$ is an integer strictly between zero and four. The possibilities are consequently:

Orthogonal roots have $p=q=0$ and $\theta =\pi /2$. Parallel roots have $pq=4$; the integer pairs would allow only length multiples $1,2,1/2$. The stated reducedness condition excludes the latter two, so parallel roots are equal or opposite, with angles zero or $\pi$ and equal lengths. ThusNot every root system realizes every listed angle.

The orthogonal pair alone gives no ratio through $q/p$. To finish the length restriction for a simple Lie algebra, use irreducibility of its root system. Let $U$ span the Weyl group orbit of $\alpha$. Both $U$ and $U^{\perp }$ are invariant under the orthogonal root reflections. If a root had nonzero projections in both, reflection of its projection in $U$ would have a nonzero $U^{\perp }$ component, violating invariance. Hence every root is entirely in one of these subspaces. A proper nonzero $U$ would split the root system into two orthogonal parts, contradicting irreducibility. This proves that the Weyl orbit spans an irreducible root space.

In particular, some orbit root $\gamma$ is nonorthogonal to any chosen $\beta$, and ${\gamma }^2={\alpha }^2$. Applying the previous integer-pair restriction to $\gamma ,\beta$ therefore also handles orthogonal $\alpha ,\beta$:There can be at most two different root lengths: three lengths would have ratios two and three relative to the shortest, but then the longest-to-middle squared-length ratio would be $3/2$, which is forbidden. A reducible collection of orthogonal rank-one factors could have independently chosen scales; the simple-algebra assumption excludes that exception.

Choose a generic linear functional to define positive roots. The simple roots are those positive roots which cannot be written as sums of two positive roots; they form a basis in which all root coefficients are integers of one sign. Distinct simple roots have nonpositive inner product: if their inner product were positive, their root string would contain their difference, and whichever difference is positive would decompose one of them. In rank two, irreducibility excludes the orthogonal case. The allowed simple-root angles are therefore $120^{\circ },135^{\circ },150^{\circ }$, giving the three irreducible diagrams in the classification of rank-two root systems: A2 root system, B2 root system and G2 root system. Reflections in the two simple roots generate their respective six, eight and twelve roots.

Here are explicit simple-root choices and positive roots; adding their negatives gives each complete diagram.

The nonnegative coefficients in each positive-root list verify the chosen simple roots explicitly. The orthogonal four-root system $A_1\times A_1$ is rank two but reducible, so it is not a diagram for the simple algebra assumed here.

Define the Cartan matrix in the column-coroot convention:Using the ordered roots above givesThe alternative row-coroot convention transposes these matrices. Stating the convention and root ordering avoids a spurious disagreement in the unequal-length cases.

Write ${\beta }_T=(k_{B}T{)}^{-1}$ for inverse temperature to distinguish it from the order-parameter critical exponent. A normalized blocking kernel $B_b({\sigma }^{\prime },\sigma )$ is a nonnegative conditional distribution of coarse spins given the fine spins, normalized byFor deterministic blocking it can impose a majority rule or block average, with a Dirac delta replacing the sum constraint for continuous coarse variables. Each block occupies $b^D$ original cells, so the new lattice spacing and number of sites areThe physical volume $N{a}^D$ is unchanged. After $p$ blockings, $a_p=b^{p}a$ and $N_p=b^{-pD}N$.

Define the coarse statistical Hamiltonian by the exact Boltzmann-weight identityThe operator family must include all generated interactions for this renormalization-group transformation to be exact; restricting it to a few couplings is an approximation. The scalar $g(u)$ has energy units and is the contribution per original site that is independent of ${\sigma }^{\prime }$. Restoring the original constant term $NC$ giveshenceThis convention measures $g$ per original site; defining a constant per new site would move a factor $b^D$ into its definition.

The identity-operator contribution to renormalization-group free energy $C_p$ records the configuration-independent normalization generated by the eliminated spins. It includes their contribution to free energy, including entropy and measure factors, and need not be independent of temperature. It has no effect on normalized probabilities of the retained spins, but is indispensable for the total partition function and its temperature derivatives. Iterating its recursion givesThus $N_p{C}_p$ is the original constant $N{C}_0$ plus the accumulated free-energy contributions of the eliminated degrees of freedom.

To obtain the additive free-energy recursion under blocking, use energy-valued free energy per current lattice site,The blocking kernel normalization impliesSubstituting the recursion for $C^{\prime }$ cancels $C$, givingA second step gives $f_0(u_0)=b^{-2D}{f}_0(u_2)+g(u_0)+b^{-D}g(u_1)$, and induction provesFor dimensionless free energy $-\log Z/N$, multiply $F,f_0,g,C$ throughout by ${\beta }_T$.

The singular free-energy density must be distinguished from the exact per-site representative $f_0$. If $a(u)$ is a regular background and $f_s=f_0-a$, thenThus the same displayed inhomogeneous recursion holds with $f_s,g_s$ after subtraction. When a suitable analytic background solves $a(u)=g(u)+b^{-D}a(Ru)$, it leaves a homogeneous singular recursion. Analytic subtraction of an inhomogeneous renormalization recursion explains this separation. A renormalization-group free-energy resonance can obstruct analytic subtraction and produce a logarithm. The original function $g$ originates in the short-distance spins that were summed out, rather than in an extra external interaction.

Near a renormalization-group fixed point $u_{\ast }$, linearize the transformation:An eigenvector of $\mathsf{A}$ with eigenvalue $b^{y_i}$ defines a scaling coordinate $t_i$ with $t_i^{\prime }=b^{y_i}{t}_i$ to leading order. A relevant operator has $y_i>0$: its coupling grows towards large distances, so it must be tuned to reach the fixed point. An irrelevant operator has $y_i<0$ and dies away; a marginal operator has $y_i=0$ and requires higher-order flow analysis. An operator of scaling dimension $x_i$ has $y_i=D-x_i$.

The coarse spins also require a multiplicative normalization. In continuum notation, writewhere ${\varphi }_{<}$ is the retained field and $\eta$ is its anomalous dimension. This restores the normalization of the critical long-distance connected correlation function. Since a uniform field conjugate to an order parameter enters through $\int d^{D}xH\varphi$, it has eigenvalue exponent $y_h=D-x_{\varphi }=(D+2-\eta )/2$.

Let $t$ be the thermal relevant coordinate with exponent $y_t>0$. After regular backgrounds have been removed, the leading singular free-energy density obeysfor a large total scale factor $B=b^p$. Choosing $B=|t{|}^{-1/y_t}$ gives the scaling hypothesis for critical phenomenaTwo temperature derivatives give the heat-capacity critical exponent; one conjugate-field derivative gives the order-parameter critical exponent; two such derivatives give the magnetic-susceptibility critical exponent. Scaling the correlation length gives $\xi (t)=B\xi (B^{y_t}t)$. Consequently,At $t=0$, instead choose $B=|H{|}^{-1/y_h}$; differentiating gives $\delta =y_h/(D-y_h)$. These scaling relations for critical exponents assume that irrelevant variables can be set to their fixed-point values without singular consequences. A dangerously irrelevant coupling or logarithmic correction requires keeping the corresponding dependence.

For the Gaussian field theory, integrate out the Fourier modes in $\Lambda /b<|q|<\Lambda$ and retain ${\varphi }_{<}$. The modes are independent at zero external field, so this Gaussian thinning transformation changes the normalization by a Gaussian determinant and leaves a quadratic statistical Hamiltonian for the retained modes. Restore the cutoff withIn the integrated statistical Hamiltonian, the volume element contributes $b^D$, each gradient contributes $b^{-1}$, and each field contributes $b^{-(D-2)/2}$. Thus the coefficient ${\kappa }^{-1}$ of the gradient energy is unchanged, whileIn particular $y_t=2$, $y_h=(D+2)/2$ and $\eta =0$. The Gaussian critical exponents and the Gaussian order-parameter scaling index therefore give the requested indicesThe sign convention $+h\varphi$ corresponds to the physical conjugate field $H=-h$ and does not affect these powers.

There are physical qualifications to these Gaussian indices. The strictly quadratic Gaussian field theory is stable for positive $m^2$ and has no stable ordered phase for negative $m^2$; its purported spontaneous $\beta$ is therefore the formal Gaussian order-parameter scaling index. A stabilizing quartic term gives an actual ordered phase. Below four dimensions that term is relevant and changes the exponents; above four dimensions it is a dangerously irrelevant coupling and gives the actual mean-field $\beta =1/2$, rather than the Gaussian hyperscaling value.

The heat capacity power can also be checked directly. Apart from constants, the zero-field Gaussian free-energy density isIf $m^2\propto t$, its singular second temperature derivative is controlled by $\int d^{D}q(q^2+\kappa m^2{)}^{-2}$. For $0<D<4$, rescaling $q=\sqrt{\kappa m^2}k$ gives a power $t^{D/2-2}$, as above. At $D=4$ there is a logarithm, and for $D>4$ an analytic background dominates the heat capacity. The negative formal Gaussian $\alpha$ then describes the subleading singular piece, with logarithmic resonances at appropriate even dimensions. This is the Gaussian specific-heat infrared threshold, and it prevents confusing singular indices with the entire observable.

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

Because $A$ is a real symmetric positive-definite matrix, its inverse exists and all its eigenvalues are positive. Completing the square givesThe translation $v=u-A^{-1}b$ has unit Jacobian determinant. Choose an orthogonal matrix $O$ with $O^{T}AO=\mathrm{diag}(a_1,\dots ,a_n)$, where $a_i>0$. The change $v=Oz$ also has absolute Jacobian one. The remaining integral is the product of $n$ one-dimensional Gaussian integrals, each equal to $\sqrt{2\pi /a_i}$. Consequently the linear-source multivariate Gaussian integral isThe determinant square root is the positive one, as required by convergence of the original real integral.

Keeping both gradient coefficients fixed requires $z^2=B{b}^2=B{c}^4$. HenceEquivalently the real-space field normalization is $z_x=z/B=b^{(3-d)/4}$. The transformed thermal and source couplings areThis is the Gaussian fixed point at $t=h=0$, with anisotropic effective volume exponent $D_{\mathrm{eff}}=(d+1)/2$. The stiffnesses remain positive parameters labeling the normalization of coordinates and fields. A strictly positive $t$ is a relevant perturbation away from the renormalization-group fixed point, not itself a critical renormalization-group fixed point.

The two correlation lengths also follow directly by balancing the kernel: $K{q}_{\parallel }^2\sim t$ gives ${\xi }_{\parallel }\sim (K/t{)}^{1/2}$, while $L{q}_{\perp }^4\sim t$ gives ${\xi }_{\perp }\sim (L/t{)}^{1/4}$. Their different exponents are consistent with the required anisotropic scaling. An omitted scalar quartic interaction has scaling eigenvalue $4-D_{\mathrm{eff}}=(7-d)/2$, so the quadratic renormalization-group fixed point alone does not establish the interacting smectic's critical behavior below physical dimension seven.

The scaling hypothesis for critical phenomena describes the singular part of the equilibrium free-energy density, rather than its analytic background. For an isotropic system in spatial dimension $D$, with thermal scaling field $t$ and conjugate field $h$, writeThis follows from coarse graining near a renormalization-group fixed point: a correlation volume is enlarged by $b^D$, while the relevant thermal and field perturbations grow with their eigenvalues. Choosing $b=|t{|}^{-1/y_t}$ givesThe two functions describe the two sides of the transition. Nonuniversal metric factors multiplying $t,h,f_s$ must be allowed when comparing microscopic models. Exponents and suitably normalized scaling functions or amplitude ratios characterize a universality class.