Use the wave-function renormalization coefficient and anomalous dimension definition stipulated in this part, keeping their signs and normalization. Put $A=1/[6(32{\pi }^2{)}^2]$. To the required order,At fixed bare ${\lambda }_4$, the leading change of the dimensionless coupling is $\mu d{g}_4/d\mu =-2ϵg_4+O(g_4^2)$, because ${\lambda }_4$ has mass dimension $2ϵ$. Therefore the chain rule, applied before setting $ϵ=0$, givesThusThe vanishing engineering contribution $-2ϵg_4$ multiplies a pole and leaves a finite answer. Dropping it at the start would incorrectly give zero. The numerical answer is conditional on the supplied two-loop coefficient. A convention defining the field anomalous dimension as ${\gamma }_{\varphi }=\frac{1}{2}\mu d\log Z/d\mu$ instead gives ${\gamma }_{\varphi }=-\gamma /2$ for the same $Z$; substituting that convention silently would answer a different question.

The change of field and parameters separates the Lagrangian density in Euclidean space intoThe two-leg counterterm interaction vertex is $-(\delta Z{p}^2+\delta m^2)$, with the overall momentum conservation delta function understood. With $D_0(p)=1/(p^2+m^2)$, the explicit first-order contribution to the propagator isThe minus sign follows either from expanding $e^{-S_{\mathrm{ct}}}$ or from expanding $(p^2+m^2+\delta Z{p}^2+\delta m^2{)}^{-1}$.

For $k=2$, the loop correction is $-{\tilde{\lambda }}_4{T}_D(m)D_0(p{)}^2/2$. Its ultraviolet pole from c isThere is no external $p^2$ dependence in the amputated tadpole diagram, so no one-loop wave-function renormalization is required. The minimal subtraction scheme therefore givesOnly the pole is subtracted: including $1-{\gamma }_E+\log 4\pi$ would impose a different finite subtraction prescription. The defined $\delta m^2$ is the coefficient in $Z{m}_0^2=m^2+\delta m^2$, not automatically the difference $m_0^2-m^2$ beyond this order.

The single-vertex scalar two-point graph at first order in $g_{2k}$ has one $2k$-valent interaction vertex. Attach the two external lines to it and pair its remaining $2k-2$ legs into $k-1$ tadpole diagrams. The graph is connected, has $L=k-1$ independent loop momenta, and is momentum independent after amputation apart from the overall Dirac delta function. In particular, first order in the coupling need not mean one loop.

The number of contractions with fixed external labels isDividing by the interaction factorial $(2k)!$ gives the inverse Feynman-diagram symmetry factor $1/[2^{k-1}(k-1)!]$. Use $m>0$ for the mass in the lowest-order propagator and define the tadpole integralThe linear-order correction to the amputated function is consequentlyThis is the correction to the paper's amputated connected correlator, not a definition of an inverse-propagator self-energy with an unspecified sign. At this order using $m_0$ instead of $m$ gives the same lowest-order integral.

For completeness, Schwinger parameterization and a Gaussian integral giveThis first converges for $0<D<2$ and then defines its dimensional regularization continuation. At $D=4-2ϵ$, the Gamma function expansion yieldsThus its highest ultraviolet pole isFor $k=2$ the remainder means finite terms. For $k>2$ the unrenormalized graph has a higher-order pole from its independent ultraviolet subdivergences, rather than only a simple pole. A massless tadpole integral is scaleless and vanishes in dimensional regularization; that does not justify discarding the massive ultraviolet divergence by setting $m=0$ before separating ultraviolet and infrared effects.

Expand the normalized Euclidean path integral with weight $e^{-S_E}$ and apply Wick's theorem. Use labelled external momenta and a fixed Fourier transform convention throughout. A line carrying momentum $q$ contributes $(q^2+m_0^2{)}^{-1}$. Each $2k$-valent interaction vertex contributes $-{\lambda }_{2k}=-{\mu }^{d_k}{g}_{2k}$: its $(2k)!$ possible contractions cancel the factorial in the Lagrangian density. A vertex also carries $(2\pi {)}^D{\delta }^{(D)}(\sum q)$ with all incident momenta incoming. Integrate each independent loop momentum with $d^{D}q/(2\pi {)}^D$ and divide by the Feynman-diagram symmetry factor for fixed external labels. There are no Minkowski factors of $i$.

For the full $N$-point Green's function, include every graph in which all vertices belong to components attached to external insertions. In particular, disconnected components carrying external insertions must be included, and the zeroth-order result is the sum of free pairings. Pure vacuum bubbles cancel against the denominator. The connected Green's function restricts to connected graphs, while a one-particle-irreducible vertex imposes the stronger restriction that removing one internal line cannot disconnect the graph. These are different objects. Odd $N$ vanishes in the symmetric vacuum because every interaction has even valence. At fixed perturbative order, multiply by the appropriate power of the coupling and sum all allowed graph topologies and external attachments.

If the overall momentum-conservation delta function is removed by convention, report the remaining function as its coefficient. The ordinary unamputated Green's function includes external scalar propagators; dividing them out defines an amputated function. The definition in c performs this amputation but does not, by itself, impose the usual one-particle-irreducibility condition.

Use mass dimensions with $[x]=-1$, $[{\partial }_{\mu }]=1$ and dimensionless action. The Lagrangian density in Euclidean space has dimension $D$, and the kinetic term fixes $2+2[{\varphi }_0]=D$. ThusA dimensionless coupling constant is thereforeHere $\mu$ is a positive regularization scale with mass dimension one. The interaction is marginal at its upper critical dimension $D_c=2k/(k-1)$, relevant below it and irrelevant above it by this engineering-dimension test. In particular, $d_2=4-D$, while $d_k=4-2k+2(k-1)ϵ$ near $D=4-2ϵ$. At this stage $g_{2k}$ is the dimensionless version of the displayed bare coupling; after introducing counterterms one may use the same scale factor for the renormalized coupling.

Take the free squared mass $m^2\ge 0$, as for the ordinary massive free theory. The quartic part of the scalar potential can be writtenIt is nonnegative for every real field value exactly whenNecessity follows from the axis ${\varphi }_2=0$ and the diagonal ${\varphi }_1={\varphi }_2$: a negative coefficient on either ray makes $V$ tend to minus infinity under rescaling. The displayed sum proves sufficiency, since $m^2({\varphi }_1^2+{\varphi }_2^2)/2$ is also nonnegative. This is a copositive matrix condition on the quadratic form in the nonnegative variables $({\varphi }_1^2,{\varphi }_2^2)$, rather than positive definiteness on all independent signed variables.

The strict inequalities $\lambda >0$, $\kappa >-\lambda$ make the quartic part grow in every nonzero direction. They are not necessary for a stable vacuum when $m^2>0$. For example, $\lambda =\kappa =0$ gives two stable free massive fields; $\lambda >0$, $\kappa =-\lambda$ has a flat quartic diagonal but a confining quadratic term. If $m^2=0$, boundary cases have flat vacuum directions and need not have a unique isolated minimum. If one instead allows negative squared mass, a flat quartic ray makes the full potential unbounded below; the stability conditions must then be strengthened. This distinguishes boundedness from strict quartic coercivity.

Under the generator of a, the variation isIt vanishes identically if and only if $\kappa =\lambda$. In that case $V_4=\lambda ({\varphi }_1^2+{\varphi }_2^2{)}^2$ preserves the entire global U(1) symmetry. More explicitly, the interacting Euler-Lagrange equations giveHenceA special solution lying on an axis does not establish an operator conservation law. Among the stable couplings, the symmetry-preserving line is $\kappa =\lambda \ge 0$.

The canonical momenta are ${\pi }_k={\dot{\varphi }}_k$, soDifferent species commute at equal times, so these products have no ordering ambiguity and $Q$ is Hermitian. Write ${\int }_{\mathbf{p}}=\int d^{3}p/(2\pi {)}^3$. Substituting the mode expansion of a free field and using the spatial Fourier transform identity givesSubtract the expression with species interchanged. The two-annihilator terms cancel after $\mathbf{p}\mapsto -\mathbf{p}$, as do the two-creator terms. Commuting different species in the remaining terms yieldsThere is no infinite vacuum contribution to the Noether charge to subtract: the contractions that would generate a species-diagonal constant are absent.

The bosonic canonical commutation relations giveDefine the creation operatorsTheir commutators with $Q$ are $[Q,c_{\pm }^{†}]=\pm c_{\pm }^{†}$. Since $Q|0⟩=0$, the one-particle states haveThe fixed-momentum states are delta-function normalized; a square-integrable wave packet of either species is a normalizable state with the same Noether charge eigenvalue. In this basis $Q={\int }_{\mathbf{p}}(c_{+}^{†}{c}_{+}-c_{-}^{†}{c}_{-})$, so it measures a difference of occupation numbers. The equal-time identity $i[Q,{\varphi }_1]={\varphi }_2$, $i[Q,{\varphi }_2]=-{\varphi }_1$ independently checks the sign against a.

Divide the infinitesimal variation by its constant parameter $\alpha$, so the transformation generator is $(\delta {\varphi }_1,\delta {\varphi }_2)=({\varphi }_2,-{\varphi }_1)$. The variation of the kinetic term isand that of ${\varphi }_1^2+{\varphi }_2^2$ is $2{\varphi }_1{\varphi }_2-2{\varphi }_2{\varphi }_1=0$. Thus the Lagrangian density is invariant pointwise, without a nonzero boundary current. The finite transformation is a rotation of the two real scalar fields, giving an SO(2) global symmetry in field theory.

Noether's theorem supplies the parameter-independent Noether currentIts divergence is ${\partial }_{\mu }{j}^{\mu }={\varphi }_2\square {\varphi }_1-{\varphi }_1\square {\varphi }_2$, since the first-derivative terms cancel. Both fields obey the same Klein-Gordon equation, so this vanishes on solutions. With spatial boundary conditions giving zero flux at infinity, the Noether charge isA finite box with periodic boundary conditions works as well. The spatial integral is conserved because the local continuity equation converts its time derivative into a boundary flux; neither vanishing local current nor constant individual fields is required. Equivalently, $\Phi =({\varphi }_1+i{\varphi }_2)/\sqrt{2}$ transforms as $\Phi \mapsto e^{-i\alpha }\Phi$, exhibiting the equivalent global U(1) symmetry and fixing the sign of the generator.

Retain the specified leading poles and write $A=a_1/ϵ$, $B=b_1/{ϵ}^2$. Solving $g=g_0+A{g}_0^2+B{g}_0^3$ recursively givesThe scale derivative must again be reexpressed in the renormalized coupling:The leading-pole consistency under coupling reparametrization shows that a finite renormalization-group beta function requires the two-loop double-pole consistency conditionThe leading two-loop pole is fixed by the one-loop simple pole; it is not an independent parameter. This is an instance of the simple-pole formulas for minimal-subtraction renormalization group functions and their higher-pole consistency relations.

The source's leading-pole truncation does not fix the finite two-loop coefficient of the beta function. For example, adding the omitted term $c_1{g}_0^3/ϵ$ leaves the double-pole condition unchanged but contributes $-2c_1{g}^3$ to the finite beta function. Thus cancelling the displayed $1/ϵ$ term is not a claim that the physical two-loop running vanishes.

At fixed regulator abbreviateThe coupling expansion is $g=g_0+A{g}_0^2+O(g_0^3)$, whose inverse is $g_0=g-A{g}^2+O(g^3)$. Using part c before removing the regulator,Therefore the regulator-free coefficients at the orders determined by the supplied expansions areThe canonical term $-ϵg$ vanishes, and the factors of $ϵ$ cancel the simple poles. The finite coefficients $a_2,z_2$ and higher positive powers of $ϵ$ do not enter these leading regulator-free terms. Inverting the coupling is essential: substituting $g_0=g$ too early would incorrectly give $-2a_1{g}^2$.

The calculation establishes finiteness only at the displayed perturbative orders. Higher-order pole cancellations require the higher-order coefficients, as the next part makes explicit.

Write ${\mathcal{D}}_0=\mu {\partial }_{\mu }+{\beta }_0(g_0){\partial }_{g_0}$. Allow a general bare coefficient ${\gamma }_0$ temporarily, so the starting equation is $({\mathcal{D}}_0-n{\gamma }_0/2){\hat{\Gamma }}^{(n)}=0$. Substitute ${\hat{\Gamma }}^{(n)}=Z_{\varphi }(g_0{)}^{-n/2}{\Gamma }^{(n)}(p,g(g_0),\mu )$ and apply the product and chain rules:Multiplying by $Z_{\varphi }^{n/2}$ identifiesThus the renormalization-group beta function transforms by the chain rule under a coupling reparametrization, while the field coefficient also records the scale dependence of the finite or divergent field normalization. For the actual bare parametrization here, substitute ${\beta }_0=-ϵg_0$ and ${\gamma }_0=0$.

The dimensionless bare parameter is $g_0={\lambda }_0{\mu }^{-ϵ}$, so direct differentiation at fixed ${\lambda }_0$ givesThe dimensionful bare vertex itself is scale independent. Reexpressing it as ${\hat{\Gamma }}^{(n)}(p,g_0,\mu )={\Gamma }_0^{(n)}(p,g_0{\mu }^{ϵ})$ introduces only the combination $g_0{\mu }^{ϵ}$. HenceThere is no multiplicative field rescaling in this bare parametrization, soThis does not assert that the interacting renormalized field has zero anomalous dimension. It says that the unrescaled bare field has no additional normalization factor under this change of the arbitrary scale.

At fixed external momenta, fixed regulator and fixed bare coupling, define the total scale derivativeThe bare vertex has no independent dependence on the arbitrary renormalization scale. Differentiating ${\Gamma }^{(n)}=Z_{\varphi }^{n/2}{\Gamma }_0^{(n)}$ therefore givesDefining the two scale-response functions as in the source yieldsThis is the Callan-Symanzik equation. Its minus sign and factor $n/2$ follow from the specified $Z_{\varphi }^{n/2}$ multiplying the bare vertex. In particular, the source's $\gamma =\mu d\log Z_{\varphi }/d\mu$ is twice the field anomalous dimension defined using $\log Z_{\varphi }^{1/2}$. A convention change must be made consistently in both definitions and equation.

To first order, invert $m_0^2=m^2[1+g/(8{\pi }^2ϵ)]$ and use $g={\lambda }_0{\mu }^{-ϵ}+O({\lambda }_0^2)$:The remainder is a perturbative statement at fixed regulator. When differentiating, hold the bare mass $m_0^2$ fixed as well as ${\lambda }_0$; otherwise the mass running would not be defined by a fixed bare theory. Since $\mu {\partial }_{\mu }{\mu }^{-ϵ}=-ϵ{\mu }^{-ϵ}$,Thus, with the source's extra factor of one half,The factor $-ϵ$ from the running dimensionless bare coupling cancels the counterterm pole. Differentiating as though $g$ were independent of $\mu$ would miss this finite result. Some conventions call $\mu d\log m^2/d\mu$ itself the mass anomalous dimension; that convention is twice the one defined here.

Source correction: equation 11 on the original page has $\Gamma (\alpha -1/2)$. The dimensionally regulated integral instead contains $\Gamma (\alpha -d/2)$. This is not an interchangeable convention: for $\alpha =2$ and $d=4-ϵ$, the correct factor $\Gamma (ϵ/2)$ has a logarithmic ultraviolet pole, whereas the printed $\Gamma (3/2)$ is finite.

To obtain the correct expression, shift the loop momentum and perform a Wick rotation. For positive Euclidean ${\mathcal{M}}^2$, introduce a Schwinger parameterization:The derivation first holds in the convergent range and defines the result elsewhere by analytic continuation. For integer $\alpha$, rotating back giveswith the branch inherited from $i0$. This is the Lorentzian continuation of the Euclidean massive loop integral.

Set $p=0$, $\alpha =1$, $d=4-ϵ$. Since $\Gamma (-1+ϵ/2)=-2/ϵ+O(1)$,Consequently the directed tadpole in part d givesThe counterterm is $-i{\delta }_{m^2}$, so cancelling only the pole in the minimal subtraction scheme requiresThis is the quartic complex-scalar mass renormalization result. The scale factors and logarithms contribute finite terms at this order and do not change the displayed pole coefficient. A real-field check gives a tadpole coefficient $\frac{1}{2}(3g/2+g/2)=g$, in agreement with the oriented complex-field computation.

Use renormalized propagators of mass $m$ and tree vertex $-ig{\mu }^{ϵ}$. The interaction counterterms have Feynman factors $-i{\delta }_{m^2}$ for a two-leg mass insertion and $-i{\mu }^{ϵ}g{\delta }_g$ for a four-leg insertion. At one loop the two-point one-particle-irreducible vertex isThe directed tadpole has $S=1$: interchanging its two ends would reverse the field labels, so the real-scalar tadpole's factor $1/2$ does not apply. It is independent of external momentum, so no wave-function counterterm is required at this order.

For the four-point vertex, label the two unstarred incoming fields by $1,2$ and the starred ones by $3,4$, with all momenta incoming. The tree contact, the three oriented bubbles and the contact counterterm giveThe first bubble has two equivalent arrows along the same direction between its vertices, giving $S=2$. In the other two channels the arrows oppose one another, and interchanging the lines does not preserve their oriented attachments, giving $S=1$. This is the oriented complex-scalar quartic graph factors rule.

The mass tadpole is order $g{m}^2$, while the four-point loop is order $g^2$. In the multiplicative coupling convention of the source,The distinction between ${\delta }_g$ and the whole additive four-leg coefficient $g{\delta }_g$ matters for the requested order counting.

The propagator carries a real-species label rather than a charge-flow arrow:Differentiating the interaction with respect to its four external real fields gives the compact tensor ruleThus four ${\varphi }_1$ legs or four ${\varphi }_2$ legs give $-3i{\lambda }_0/2$, while two of each give $-i{\lambda }_0/2$. Vertices with an odd number of either species vanish. For example, $4!({\lambda }_0/16)=3{\lambda }_0/2$, whereas $2!2!(2{\lambda }_0/16)={\lambda }_0/2$. These factors already account for assignments at one vertex; graph symmetry factors must still be included separately. Momentum conservation and loop measures are the same as in part a.

The two real fields have canonical normalization becauseThereforeThe interaction is invariant under $O(2)$ rotations and reflections of $({\varphi }_1,{\varphi }_2)$. The continuous rotations are the original complex-field $U(1)$ symmetry. The factor $1/16$ follows from both the original $1/4$ and the square of the $1/2$ in ${\varphi }^{\ast }\varphi$.

A free contraction must connect one $\varphi$ and one ${\varphi }^{\ast }$ field. In momentum space,Represent the nonzero propagator by a directed line so the two conjugate ends remain distinguishable; an undirected line alone would incorrectly allow same-charge contractions. Choose one consistent arrow convention, for example pointing from the starred end to the unstarred end.

The interaction contains two identical $\varphi$ fields and two identical ${\varphi }^{\ast }$ fields. Their $2!2!$ assignments cancel the coefficient $1/4$, givingThe vertex has two incoming and two outgoing arrows. Include $(2\pi {)}^d{\delta }^d(\sum p_i)$ at every vertex, integrate $d^{d}k/(2\pi {)}^d$ for each independent loop, and divide by the oriented graph's symmetry factor. There is no fermion-loop sign.

Let ${\mathcal{V}}_4=i{\Gamma }^{(4)}$ denote the amputated Feynman vertex, with tree term $v_4=-i\lambda$. A one-loop one-particle-irreducible Feynman diagram has $V_3+2V_4=4$, allowing three families: two quartic vertices, two cubic vertices and one quartic vertex, or four cubic vertices. DefineUsing $I_2,I_3$ from part c and writing $P_{ij}=p_i+p_j$, the complete one-loop result isIn each term of the middle sum, $k,l$ are the two complementary indices. Their order does not change the integral after a loop-momentum redefinition.

The first line contains three bubble channels, one for each partition of the four labelled external legs into two pairs. Each has $S=2$ from interchanging its two identical internal lines. In the second line, the quartic vertex carries a chosen pair of external legs, giving $\left(\genfrac{}{}{0ex}{}{4}{2}\right)=6$ mixed triangles. Each has $S=1$. In the last line, cyclic orderings of four labelled external legs are identified under rotation and reversal, leaving $(4-1)!/2=3$ boxes, again with $S=1$ each. The figure in part b shows one representative of each family and its multiplicity. External-leg corrections and tadpoles connected by bridges are excluded because ${\Gamma }^{(4)}$ is a one-particle-irreducible vertex. All integrals are left unevaluated, as appropriate here.