Solution (source code)

= Solution

The <renormalization-group beta function> is
$$
\boxed{\beta(\lambda)=\left.\mu\frac{d\lambda}{d\mu}\right|_{\text{bare parameters fixed}}.}
$$
At one-loop accuracy, integrating $d\lambda/d\log\mu=b\lambda^2$ gives
$$
\boxed{\frac1{\lambda(\mu)}=\frac1{\lambda(\mu_0)}-b\log\frac\mu{\mu_0}.}
$$
For a positive weak coupling and $b>0$, the coupling increases toward the ultraviolet and decreases toward the infrared. Extrapolation predicts a <Landau pole> at $\mu=\mu_0\exp[1/(b\lambda(\mu_0))]$, where the weak-coupling approximation fails. For $b<0$, the positive coupling decreases toward zero as $\mu\to\infty$: this is <asymptotic freedom>. Toward the infrared it instead grows, so perturbation theory eventually fails. The zero coupling is a fixed point. If negative couplings are allowed mathematically, their direction of motion still satisfies $d\lambda/d\log\mu=b\lambda^2$, but the side on which the pole occurs reverses: for example $b>0,\lambda<0$ approaches zero from below in the ultraviolet. Negative quartic coupling does not give a stable scalar potential.

For the loop calculation use $d=4-\epsilon$, $\epsilon>0$, and the PDF's vertex normalization $-\lambda\phi^4/4!$. A <Feynman parameter> combines the denominators. Shift the Euclidean <loop momentum> to $\ell=k+xp$, and put $M_x^2=m^2+x(1-x)p^2$. The regulated integral without the vertex prefactor is
$$
I(p)=\mu^\epsilon\int_0^1dx\int\frac{d^d\ell}{(2\pi)^d}\frac1{(\ell^2+M_x^2)^2}.
$$
Using <Schwinger parameterization>, $q^{-2}=\int_0^\infty t e^{-tq}\,dt$, and the momentum <Gaussian integral> gives
$$
I(p)=\frac{\mu^\epsilon}{(4\pi)^{d/2}}\Gamma\left(2-\frac d2\right)
\int_0^1dx\,(M_x^2)^{d/2-2}.
$$
For $d=4-\epsilon$, $\Gamma(\epsilon/2)=2/\epsilon+O(1)$, while the other factors have finite limits. Thus the <massive scalar bubble pole in four dimensions> is
$$
\boxed{I(p)_{\rm div}=\frac1{8\pi^2\epsilon},\qquad
\left(\frac{i\lambda^2}{2}I(p)\right)_{\rm div}
=\frac{i\lambda^2}{16\pi^2\epsilon}.}
$$
This pole is independent of external <momentum> and mass and hence is canceled by a local quartic <counterterm>. There are three channels, corresponding to the three pairings of four external momenta. Their poles add, and a counterterm vertex $-i\delta\lambda$ cancels the sum when
$$
\boxed{\delta\lambda=\frac{3\lambda^2}{16\pi^2\epsilon},\qquad
\mathcal L_{\rm ct}^{(4)}=-\mu^\epsilon\frac{\delta\lambda}{4!}\phi^4.}
$$
This is the total counterterm for the three four-point diagrams specified here. The <minimal subtraction scheme> removes their poles only. The one-loop tadpole is momentum independent and requires a mass counterterm, but does not produce one-loop <wave-function renormalization>; that is why it does not alter this one-loop coupling calculation.

The dimensionful <bare coupling> is therefore
$$
\lambda_B=\mu^\epsilon\left[\lambda+\frac{a\lambda^2}{\epsilon}+O(\lambda^3)\right],
\qquad a=\frac3{16\pi^2}.
$$
Its independence of $\mu$ gives
$$
0=\epsilon\left(\lambda+\frac{a\lambda^2}{\epsilon}\right)
+\beta_d\left(1+\frac{2a\lambda}{\epsilon}\right)+O(\lambda^3).
$$
Putting $\beta_d=-\epsilon\lambda+b\lambda^2+O(\lambda^3)$, the finite coefficient of $\lambda^2$ is $b-a$. Thus
$$
\boxed{\beta_d=-\epsilon\lambda+\frac{3\lambda^2}{16\pi^2}+O(\lambda^3),\qquad
\beta_{d=4}=\frac{3\lambda^2}{16\pi^2}+O(\lambda^3).}
$$
This is the <one-loop quartic scalar beta function>, positive for nonzero weak real coupling. The factors would look different with a regulator $d=4-2\epsilon$, but the four-dimensional beta function is unchanged. For massless theory take nonexceptional external momentum, or an infrared regulator, so an infrared singularity is not confused with the ultraviolet pole.