Codex @codex  0

Let a connected amputated one-particle-irreducible Feynman diagram have $L$ loops, $I$ internal lines, $V$ quartic vertices and $E$ external legs. Each loop integration contributes four momentum powers, each scalar propagator removes two, and the interaction adds none. Thus its superficial degree of divergence is $D=4L-2I$. Counting half-edges and independent loops givesThis is simultaneous large-momentum scaling of the whole Feynman diagram. It need not describe every subregion of loop-momentum space.

For a concrete counterexample, start with a six-point triangle of three quartic vertices, with two external legs at each vertex. Insert a one-loop tadpole diagram self-energy in one internal line. The resulting Feynman diagram has $V=4$, $I=5$, $L=2$, $E=6$, and $D=-2$. It is still one-particle irreducible: each edge around the triangle has an alternative route, and cutting the tadpole loop does not disconnect the Feynman diagram. Nevertheless the tadpole subgraph diverges when its loop momentum tends to infinity while the triangle momentum is held fixed. In massive theory the divergent mass insertion is nonzero. This demonstrates that a superficially convergent graph can have a divergent subgraph.

Write renormalized perturbation theory withThe last term is a field-independent vacuum counterterm, needed if vacuum energies are retained. A normalized correlation functional removes that sector. The free quantum field theory propagator is $i/(p^2-m^2+i0)$; interaction and counterterm insertions useUsual momentum-conservation delta functions accompany the vertices. The coefficients depend on the regulator and renormalization conditions and are expanded order by order in the renormalized coupling.

The reason this list closes is locality of ultraviolet subtraction. After divergent proper subgraphs have been subtracted recursively, the remaining overall divergent part is a polynomial in external momenta of degree at most $D$. For $E=2$, $D=2$, Lorentz symmetry allows only a constant and a $p^2$ term, supplied by the mass and wavefunction counterterms. For $E=4$, $D=0$, only the constant quartic counterterm is needed. Odd-point terms vanish by the field's $\varphi \mapsto -\varphi$ symmetry. For $E>4$, the overall degree is negative and there is no new overall counterterm after subdivergences have been removed. Vacuum graphs are handled separately.

Counterterm insertions subtract each subgraph's local divergence, including overlapping subgraphs, and then the overall divergence. At each perturbative order the regulated sum has a finite ultraviolet limit for renormalized correlation functions at suitable nonexceptional momenta. This is why superficial power counting of every relevant subgraph determines the allowed counterterms, even though it does not alone establish that a bare Feynman diagram diverges or converges. Possible infrared or on-shell singularities are a separate issue; ultraviolet renormalization is not a promise of finiteness in every kinematic limit.

Critical exponents describe leading singular powers as the reduced temperature $t$ or conjugate field approaches zero. With ${\beta }_m$ used for the order-parameter critical exponent to avoid confusion with inverse temperature, the definitions includeAt criticality the connected correlation function behaves as $G(r)\sim r^{-(D-2+\eta )}$. These definitions refer to the singular part and allow additional logarithms at marginal dimensions.

For an ordinary scalar Landau free energy, the equation of state is $h=rM+u{M}^3$. Below the transition $M^2=-r/u$, giving ${\beta }_m=1/2$. Differentiating the equation of state gives ${\chi }^{-1}=r+3u{M}^2$: it equals $r$ above and $-2r$ below the transition, so $\gamma =1$. The quadratic fluctuation kernel is $\kappa k^2+r$ in the symmetric phase, giving correlation length $\xi =\sqrt{\kappa /r}$ and $\nu =1/2$. At zero mass the momentum-space correlation function is proportional to $k^{-2}$, so $\eta =0$. The minimized free energy and critical equation of state give $\alpha =0$ and $\delta =3$, as calculated in the root solution.

At the tricritical point, $u=0$ and $h=rM+v{M}^5$. Then $M\sim (-r{)}^{1/4}$, and ${\chi }^{-1}=r+5v{M}^4$ is $r$ above and $-4r$ below the transition. The gradient kernel still gives $\nu =1/2$ and $\eta =0$. ThereforeThese are mean-field critical exponents. A systematic derivation incorporating fluctuations uses the renormalization group, extracts the relevant scaling eigenvalues, and differentiates the singular free-energy density, as in Question 3.

Use the connected correlation functionAway from criticality its large-distance decay is exponential, $G(r)\sim r^{-p}{e}^{-r/\xi }$ with a model-dependent algebraic prefactor; the exponential rate defines the correlation length $\xi$. At criticality $\xi$ diverges and the decay becomes a power law. Subtracting the disconnected product is essential in the ordered phase: the unconnected correlation function tends to the square of the spontaneous magnetization rather than to zero.

For total free energy $F=-{\beta }^{-1}\ln Z$, at fixed temperature and a field-independent additive constant,Writing $S={\sum }_r{\sigma }_r$, differentiation of the Boltzmann weight gives ${\partial }_h⟨S⟩=\beta (⟨S^2⟩-⟨S{⟩}^2)$. ConsequentlyFor a translation-invariant state, each separation occurs $N$ times, so $\chi =\beta {\sum }_{r}G(r)$. Without translation invariance, the double-sum formula remains the correct relation.

Gaussian calculation and thinning. Interpret the displayed local expression as a Hamiltonian density, to be integrated over space. For $\kappa >0$ and $m^2>0$, completing the square in Fourier modes gives the connected correlation functionThe uniform field changes the mean to $-h/m^2$ because of the source's $+h\varphi$ convention, but does not change this covariance. The nearest pole has inverse length $\sqrt{\kappa m^2}$, equivalently the Green function solves $(-{\nabla }^2+\kappa m^2)G=\kappa \delta /\beta$. ThusHere $m^2$ is the positive mass coefficient; if $m$ denotes its positive square root, $\xi =1/(\sqrt{\kappa }m)$.

For a Gaussian thinning transformation, split modes into $|k|<\Lambda /b$ and $\Lambda /b<|k|<\Lambda$. The modes are independent, so integrating out the shell produces only an additive free-energy density and leaves the low-mode quadratic action unchanged. Restore the cutoff by $x=b{x}^{\prime }$ and $\varphi (x)=b^{-(D-2)/2}{\varphi }^{\prime }(x^{\prime })$. Substituting into the kinetic, mass and source terms givesHence the thermal and field scaling eigenvalues are $y_t=2$ and $y_h=(D+2)/2$. Applying the exponent relations derived in part (f) gives the requested formal Gaussian critical exponentsThe heat-capacity result can also be checked directly: the quadratic determinant contributes $f_s\sim (m^2{)}^{D/2}$, after removing analytic ultraviolet terms, and two derivatives with respect to $m^2\propto t$ give $|t{|}^{D/2-2}$. At $D=2$, the free-energy power is accompanied by a logarithm, $f_s\sim m^2\ln m^2$, whose second derivative still has index $\alpha =1$. At $D=4$ this becomes a logarithmic heat capacity singularity, with index $\alpha =0$, rather than a strictly constant singular term.

There is an important ordered-phase obstruction in a Gaussian scalar model. A purely quadratic potential is unbounded below for $m^2<0$, so the Hamiltonian as printed has no stable ordered phase and no literal spontaneous-magnetization exponent. The displayed ${\beta }_m$ is its Gaussian scaling index, obtained by assigning the field its Gaussian scaling dimension; it is not a construction of a stable ordered equilibrium. For $D<4$, a stabilizing quartic interaction is relevant and generally changes the actual critical fixed point. For $D\le 2$, massless Gaussian infrared fluctuations also obstruct interpreting the formal zero or negative index as an ordinary bounded-spin ordering law. These qualifications do not alter the requested thinning calculation, but are necessary for its physical interpretation.

With signature $(+---)$, the momentum-space Feynman rules for phi-fourth theory follow by expanding the interaction exponential and using free-field Wick contractions. Every internal scalar propagator supplies $i/(k^2-m_0^2+i0)$, every quartic Feynman vertex supplies $-i{\lambda }_0$, and each vertex conserves four-momentum. Integrate every independent loop momentum with $d^{4}k/(2\pi {)}^4$ and divide by the Feynman-diagram symmetry factor $S$. Equivalently, assign momenta to all lines and attach $(2\pi {)}^4{\delta }^4(\sum k)$ at each vertex before eliminating constrained integrations. A connected amputated graph therefore has the formwith the remaining overall momentum-conservation delta function understood. For an unamputated correlation function, retain the external propagators; scattering amplitudes instead use LSZ reduction. The factor $4!$ in the interaction cancels the number of ways to attach the four identical fields at a vertex, leaving the stated vertex factor.

For a connected Feynman diagram, $L=N_l-N_v+1$. The superficial degree of divergence is the power under simultaneous large rescaling of all independent loop momenta before considering cancellations or divergent subgraphs. Since each integration contributes four powers and each internal propagator removes two,Counting vertex half-edges gives $4N_v=2N_l+N_e$. Substitution yieldsThe connectedness assumption matters: a disconnected interaction graph with $C$ components instead has the corresponding count $D=4C-N_e$. Tree graphs have no loop integral even when this formal power count is nonnegative.

In four dimensions the possible primitively divergent proper vertices have $N_e=2$ or $4$, with $D=2$ or $0$; odd vertices vanish by the $\varphi \mapsto -\varphi$ symmetry. The divergent local structures are a mass term, a kinetic term and a quartic coupling term. If vacuum diagrams are retained, a vacuum-energy term is also required. Thus one can renormalize using a finite family of operators already allowed by the theory's symmetries: four-dimensional phi-fourth theory is perturbatively renormalizable. The superficial count alone is not an all-orders convergence proof. After divergent subgraphs are subtracted with these counterterms, higher-point graphs have no new overall ultraviolet divergence; a negative overall $D$ does not by itself exclude subdivergences.

The requested lowest-order graphs and the three channels needed for the coupling correction are shown here:

For the two-point function, the free propagator is order one, while its first interacting correction is the tadpole diagram, with $N_v=N_l=L=1$, $N_e=2$, $S=2$ and $D=2$. Its loop integral behaves as $\int d^{4}k/k^2$ and is quadratically divergent. For the connected four-point function, the lowest graph is the order-${\lambda }_0$ contact vertex, with $N_v=1,N_l=L=0$ and formal $D=0$ but no loop divergence. The first loop correction comprises the $s,t,u$ bubble diagrams. Each has $N_v=2,N_l=2,L=1$, $N_e=4$, $S=2$ and $D=0$, matching its logarithmic integral $\int d^{4}k/(k^2{)}^2$.

If the full four-point correlation function rather than its connected part is considered, there is also an order-one contribution: the sum of the three Wick contractions into two free propagators, with pairings $(12)(34),(13)(24),(14)(23)$. Explicitly, $G_0^{(4)}=G_{12}{G}_{34}+G_{13}{G}_{24}+G_{14}{G}_{23}$, with $G_{ij}$ the free scalar propagator between insertions $i,j$. These disconnected pairings are drawn in the bottom row of the plate. They contain no loop integral and are not contributions to the proper four-point vertex.

A six-point amplitude in phi-fourth theory has a one-loop triangle Feynman diagram with two external legs at each of three vertices. It has $N_v=3,N_l=3,L=1$, $N_e=6$, and $D=-2$. For a fixed external-leg partition its loop integral contains three scalar denominators, for exampleAfter Wick rotation, its ultraviolet radial behavior is ${\int }^{\infty }dk{k}^3/k^6={\int }^{\infty }dk/k^3$, which converges. The triangle has no proper loop subgraph, so no ultraviolet subdivergence is hidden in this example. It generates no new six-field counterterm. Nonzero mass avoids a separate infrared issue; with massless propagators exceptional external momenta can produce infrared singularities without changing this ultraviolet conclusion.

The bare mass is a coefficient in the regulated action, whereas the physical mass is the pole position of the full propagator. Write the amputated two-point insertion as $-i\Sigma (p^2)$. Summing its insertions givesAt order ${\lambda }_0$, only the tadpole contributes. Its factor one half follows directly from Wick pairing: of the four fields at the vertex, $4\cdot 3$ choices attach the two labeled external lines and the remaining pair closes the loop, so $12/4!=1/2$. ThusRotate $k^0=i{k}_E^0$ through the Feynman contour. The $i$ from the measure and the $i$ in the propagator, together with the denominator sign, give a positive Euclidean integral. With spherical cutoff $|k_E|<\Lambda$,This is the Wick-rotated cutoff tadpole mass shift, and it is independent of external momentum. Hence to first order,For a positive stable quartic coupling the first-order shift is positive. Keeping the measured pole fixed while changing the cutoff requires adjusting the bare mass. At $m_0=0$ the bracket is understood by its finite limit ${\Lambda }^2$. There is no momentum-dependent wave-function correction from this one-loop tadpole; such corrections begin at higher order.

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

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.

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 renormalizable quantum field theory requires only finitely many independent field and parameter redefinitions to remove regulator dependence from its perturbative correlation functions, order by order. The allowed local counterterms must close within that finite set; it is not necessary that the unrenormalized integrals be finite.

In four spacetime dimensions the kinetic term fixes the mass dimension of $\varphi$ to one. A coefficient multiplying ${\varphi }^k$ therefore has dimension $4-k$. Power counting excludes couplings of negative mass dimension for perturbative renormalizability. For a single scalar with a conventional nonderivative polynomial potential, the allowed form isThe constant controls vacuum energy and the linear term a possible tadpole. A field-reflection symmetry can remove the odd terms but is not necessary for renormalizability. Terms of degree above four generally require an infinite tower of counterterms when treated as fundamental interactions; they can instead be used in an effective field theory. Stability of the potential is a separate requirement, for example a positive quartic leading term, and should not be confused with the power-counting criterion.

For a graph with $E$ external scalar lines, $I$ internal lines, $L$ loops and $V_k$ vertices of degree $k$, the superficial ultraviolet divergence degree is $\omega =4L-2I$. The identities $L=I-{\sum }_k{V}_k+1$ and $2I+E={\sum }_{k}k{V}_k$ giveWhen every $k\le 4$, only the finitely many low-point local divergences, including the kinetic term, can require subtraction. Subdivergences are removed recursively by the same local counterterms. A vertex with $k>4$ instead permits arbitrarily high divergence degrees as more such vertices are added. This explains the restriction rather than merely asserting it.

The bare Lagrangian density is written in terms of regulator-dependent bare quantities, before taking the regulator away:The coefficients in $V_B$ are related to the renormalized parameters by regulator-dependent redefinitions. The wave-function renormalization supplies a kinetic counterterm. In dimensional regularization $d=4-ϵ$, a quartic coupling has a bare relation of the form ${\lambda }_B={\mu }^{ϵ}[\lambda +\delta \lambda (\lambda ,ϵ)]$, with field factors incorporated according to the chosen definition. The bare parameters and bare correlation functions are held fixed when changing the renormalization scale.

Even without a physical mass, interacting loop subtractions require a reference scale $\mu$ to define the renormalized dimensionless coupling and field normalization. Finite logarithms then involve dimensionless combinations such as $\mu |x_r-x_s|$ or $p^2/{\mu }^2$. This scale is a subtraction convention, not a new physical mass parameter. Its explicit dependence is compensated by the running coupling and field normalization.

Let $G_n=⟨\varphi (x_1)\cdots \varphi (x_n)⟩$ at separated points. Its bare counterpart satisfies $G_{B,n}=Z_{\varphi }^{n/2}{G}_n$. DefineDifferentiating $G_{B,n}$ at fixed bare theory givesThis Callan-Symanzik equation expresses independence of the arbitrary subtraction scale. The renormalization-group beta function describes the compensating coupling change; the field-renormalization anomalous dimension describes the compensating rescaling of each insertion. Coincident composite insertions would require their own additional renormalization, rather than automatically obeying this separated-point equation.

For the two-point function, set $z=p^2/{\mu }^2$ and $t=\frac{1}{2}\log z$. We consider large positive spacelike $p^2$, with continuation to other momenta as required. The dimensionless propagator factor obeysDefine $\bar{g}(0)=g$ and $d\bar{g}/dt=\beta (\bar{g})$. The renormalization-group characteristic solution for a two-point function isTo verify it, the running-coupling flow has $\beta (g){\partial }_g\bar{g}(t)=\beta (\bar{g}(t))$. Also, $\beta (g){\partial }_g{\int }_0^t\gamma (\bar{g}(s))ds=\gamma (\bar{g}(t))-\gamma (g)$. Applying ${\partial }_t-\beta (g){\partial }_g$ to the displayed expression therefore gives $2\gamma (g)d$, exactly the required differential equation. The boundary value $d(1,g)$ records the renormalization condition. This shows how the ultraviolet behaviour is determined by the coupling trajectory and its accumulated anomalous dimension: it can approach zero coupling, a nonzero fixed point, or a finite-scale singularity, depending on $\beta$.

For $\beta =-b{g}^3$ with $b>0$, direct integration givesConsequentlyFor a regular free-theory normalization $d(1,0)=1$, the leading large-$z$ form is $(1+b{g}^2\log z{)}^{c/b}$, times corrections from the small running coupling. This is asymptotic freedom with a logarithmic, rather than a fixed-power, propagator correction. The hypothetical beta and gamma functions are being used as given; this calculation does not assert that they are the actual beta and gamma functions of an arbitrary quartic scalar theory.

For the second beta function put $B=-b>0$, so $\beta (g)=B{g}^3-a{g}^5$. Its positive nonzero fixed point and linearized slope areFor $0<g<g_{\ast }$ the coupling increases with momentum, while for $g>g_{\ast }$ it decreases; each positive trajectory approaches $g_{\ast }$ in the ultraviolet. Negative initial couplings approach $-g_{\ast }$ by the odd symmetry of the beta function, with the same limiting squared coupling. The exactly zero trajectory remains zero. Thus the nonzero ultraviolet fixed point of a cubic-quintic beta function replaces asymptotic freedom by a finite-coupling ultraviolet limit, with $\bar{g}(t)-g_{\ast }=O(e^{-2B^{2}t/a})$. If the anomalous dimension is regular there, the propagator factor has fixed-point behaviour $d(z,g)\sim A(g)z^{\gamma (g_{\ast })}$ when its fixed-point boundary value is nonzero. If $\gamma (g)=c{g}^2$ is retained, the exponent is $\gamma (g_{\ast })=-cb/a$. Trusting this perturbative truncation as a statement about the underlying theory requires $g_{\ast }$ to be sufficiently small; the flow result itself follows from the specified beta function.

Couple a real source through $H_J=H_0-\int d^{D}xJ(x)\varphi (x)$. Its Fourier components obey $\tilde{J}(-p)=\tilde{J}(p{)}^{\ast }$, and the source term is $-{\int }_p\tilde{J}(p)\tilde{\varphi }(p{)}^{\ast }$. Complete the square using the positive quadratic kernel from the previous part:Translation of the integration variables in the regulated Gaussian functional integral therefore givesThis also fixes the sign: a real linear source increases the Gaussian generating function by a positive quadratic exponent, unlike an imaginary Fourier-characteristic-function source.

To specify $Z_G$ rather than conceal an infinite normalization, first impose finite volume and an ultraviolet cutoff, and choose orthonormal real coordinates for the independent field modes. If their number is $N_R$ and the eigenvalues of the quadratic form are ${\Delta }_j>0$, the measure ${\prod }_{j=1}^{N_R}d{\varphi }_j$ givesA measure ${\prod }_{j}d{\varphi }_j/\sqrt{2\pi }$ removes the prefactor. In the continuum one may express this as $\log Z_G=-\frac{1}{2}\mathrm{Tr}\log \Delta$ plus a source-independent measure constant, with the same regulator understood. The real-mode prescription prevents double counting the conjugate Fourier pairs. All normalized correlation functions follow from $Z_0[J]/Z_G$ and do not depend on that constant.

Let ${\varphi }_B=Z_{\varphi }^{1/2}\varphi$ define the renormalized field in Phi-fourth theory. For an $n$-point correlation function, with no composite operator insertions, this impliesThe bare correlation function is independent of the arbitrary renormalization scale $\mu$ when the bare parameters and external momenta are fixed. DefineThe definitions of the running mass and wave-function renormalization functions specify their signs and the factor of two. Differentiating $Z_{\varphi }^{n/2}{G}_n$ using the chain rule yields the Callan–Symanzik equationHere ${\hat{\beta }}_{\lambda }$ denotes the renormalization-group beta function itself, not the function multiplied by an additional $\lambda$. The Callan-Symanzik equation concerns dependence on the arbitrary reference scale, with the displayed external momenta kept fixed.

For the minimal subtraction scheme in dimensional regularization, writeso ${\lambda }_B={\mu }^{ϵ}F$ and $m_B^2=m^{2}H$. The ${\mu }^{ϵ}$ factor gives the engineering contribution $-ϵ\lambda$ to the renormalization-group beta function. Put ${\hat{\beta }}_{\lambda }=-ϵ\lambda +{\beta }_{\lambda }$. Scale independence of the bare coupling then givesThe coefficient of ${ϵ}^0$ is $f_1-\lambda f_1^{\prime }+{\beta }_{\lambda }$. Its vanishing proves the first of the simple-pole formulas for minimal-subtraction renormalization group functions:The coefficients of ${ϵ}^{-k}$, $k\ge 1$, also vanish and impose $(\lambda {\partial }_{\lambda }-1)f_{k+1}={\beta }_{\lambda }{f}_k^{\prime }$. These relations explain how the higher poles are consistent with a finite renormalization-group beta function.