---
abstract: |
    The exact NSVZ relation between a $\beta$-function of ${\cal
    N}=1$ SQED and an anomalous dimension of the matter superfields is studied within the Slavnov higher derivative regularization approach. It is shown that if the renormalization group functions are defined in terms of the bare coupling constant, this relation is always valid. In the renormalized theory the NSVZ relation is obtained in the momentum subtraction scheme supplemented by a special finite renormalization. Unlike the dimensional reduction, the higher derivative regularization allows to fix this finite renormalization. This is made by imposing the conditions $Z_3(\alpha,\mu=\Lambda)=1$ and $Z(\alpha,\mu=\Lambda)=1$ on the renormalization constants of ${\cal N}=1$ SQED, where $\Lambda$ is a parameter in the higher derivative term. The results are verified by the explicit three-loop calculation. In this approximation we relate the $\overline{\mbox{DR}}$ scheme and the NSVZ scheme defined within the higher derivative approach by the finite renormalization.
author:
- |
    A.L.Kataev\
    \
    ,\
    \
    K.V.Stepanyantz\
    ,\
title: '**NSVZ scheme with the higher derivative regularization for ${\cal N}=1$ SQED**'
---

Keywords: higher covariant derivative regularization, supersymmetry, $\beta$-function, subtraction scheme.

Introduction
============

Supersymmetry is a beautiful concept of quantum field theory [@Ogievetsky:1975nu], which is studied both by theoreticians and by experimentalists. Investigations of SUSY models reveal interesting theoretical features (see e.g.[@West:1990tg; @Buchbinder:1998qv; @Diakonov:2010zzb]). Among them is the finiteness of the ${\cal N}=4$ supersymmetric Yang Mills (SYM) theory. As the consequence, this theory obeys the conformal symmetry even at the quantum level due to vanishing of the renormalization group (RG) $\beta$-function. This was proved in all orders of the perturbation theory in [@Grisaru:1982zh; @Mandelstam:1982cb; @Brink:1982pd; @Howe:1983sr] after the explicit analytical three-loop calculation of Ref. [@Avdeev:1980bh], made with the dimensional reduction (DRED) regularization [@Siegel:1979wq] (this calculation was also confirmed in [@Grisaru:1980nk; @Caswell:1980ru]). DRED is a modification of the dimensional regularization method [@'tHooft:1972fi; @Bollini:1972ui; @Ashmore:1972uj; @Cicuta:1972jf]. It is applied in SUSY theories since the dimensional regularization explicitly breaks SUSY [@Delbourgo:1974az]. Renormalization of SUSY theories is usually made using a variant of the modified minimal subtraction scheme $\overline{\mbox{MS}}$ [@Bardeen:1978yd][^1], namely in the $\overline{\mbox{DR}}$ scheme. However, it is well known that DRED is not mathematically consistent [@Siegel:1980qs]. Its inconsistency leads to the loss of explicit SUSY [@Avdeev:1981vf], which will be broken by quantum corrections in higher loops [@Avdeev:1982np; @Avdeev:1982xy]. Nevertheless, a recent analytical calculation [@Velizhanin:2010vw] demonstrates that the DRED regularization and the $\overline{\mbox{DR}}$ scheme are working in the ${\cal N}=4$ SYM theory up to the four-loop approximation and give the vanishing result for its $\beta$-function.

SUSY also leads to the absence of radiative corrections to the $\beta$-function of ${\cal N}=2$ SYM theories starting from the two-loop approximation [@Howe:1983wj]. In this case vanishing of the two- and three-loop contributions to the $\beta$-function was explicitly demonstrated in [@Avdeev:1981ew] using the $\overline{\mbox{DR}}$ scheme. However, the calculations made in Ref. [@Avdeev:1982np], even after correcting them in Refs. [@Jack:2007ni; @Velizhanin:2008rw], show that the $\beta$-function obtained from the fermion-fermion-vector, scalar-scalar-vector and ghost-ghost-vector vertexes vanishes, but obtaining the $\beta$-function from the fermion-fermion-scalar vertex gives a non-vanishing result in the three-loop approximation [@Velizhanin:2008rw]. This result reveals unsolved theoretical problems which appear if DRED is used for the regularization of SUSY theories.

${\cal N}=1$ SUSY leads to the absence of divergent quantum corrections to the superpotential [@Grisaru:1979wc]. Another interesting feature of ${\cal N}=1$ SYM models is existence of the exact relation between the $\beta$-function and the anomalous dimension of the matter superfields, derived in Refs. [@Novikov:1983uc; @Novikov:1985rd; @Shifman:1986zi]. This relation is usually called “the exact Novikov, Shifman, Vainshtein, and Zakharov (NSVZ) $\beta$-function”. It was first obtained in terms of the non-renormalized bare coupling constant $\alpha_0$ for SYM models in Refs. [@Novikov:1983uc; @Novikov:1985rd] and for the ${\cal N}=1$ supersymmetric electrodynamics (SQED) in Refs. [@Vainshtein:1986ja; @Shifman:1985fi]. For ${\cal N}=1$ SQED the NSVZ $\beta$-function is

$$\label{NSVZ}
\beta(\alpha_0) =
\frac{\alpha_0^2}{\pi}\Big(1-\gamma(\alpha_0)\Big).$$

However, a general problem arises, whether it is possible to specify a concrete subtraction scheme, which leads to the exact NSVZ $\beta$-function, if the RG functions are defined through the renormalized coupling constant. With DRED, supplemented by the $\overline{\mbox{DR}}$ subtraction scheme, this problem was studied in Refs. [@Jack:1996vg; @Jack:1996cn; @Jack:1998uj; @Jack:2007ni]. The results of the one- [@Ferrara:1974pu] and two-loop [@Jones:1974pg] calculations agree with the NSVZ $\beta$-function, because a two-loop $\beta$-function and a one-loop anomalous dimension are scheme independent in theories with a single coupling constant. In higher orders [@Jack:1996cn; @Jack:2007ni; @Harlander:2006xq] the exact NSVZ $\beta$-function for the RG functions defined in terms of the renormalized coupling constant can be obtained with the $\overline{\mbox{DR}}$ scheme after an additional finite renormalization. This finite renormalization should be fixed in each order of the perturbation theory, starting from the three-loop approximation [@Jack:1996vg; @Jack:1996cn; @Jack:1998uj; @Jack:2007ni]. However, there is no general prescription, how one should construct this finite renormalization using the $\overline{\mbox{DR}}$ scheme. Investigations of quantum corrections using other regularizations [@Shifman:1985tj; @Mas:2002xh] are usually made in one- and two-loop approximations.

In this paper we study in details, how the exact relation (\[NSVZ\]) can be obtained using the Slavnov higher derivative regularization [@Slavnov:1971aw; @Slavnov:1972sq], which is mathematically consistent and does not break the supersymmetry in all orders [@Krivoshchekov:1978xg; @West:1985jx]. Therefore, this regularization has theoretical advantages over the mathematically inconsistent DRED. Application of the higher derivative regularization to the evaluation of quantum corrections in ${\cal N}=1$ SUSY theories reveals one more interesting feature: integrals needed for obtaining a $\beta$-function defined in terms of the bare coupling constant are integrals of total derivatives [@Soloshenko:2003nc; @Pimenov:2009hv; @Stepanyantz:2011zz; @Stepanyantz:2011wq] and even double total derivatives [@Smilga:2004zr; @Stepanyantz:2011bz; @Stepanyantz:2012zz; @Stepanyantz:2012us]. This implies that one of the loop integrals can be calculated analytically, and a $\beta$-function in a $L$-loop approximation can be related with an anomalous dimension of the matter superfields in the $(L-1)$-loop approximation [@Stepanyantz:2011jy]. As a consequence, the NSVZ $\beta$-function can be naturally obtained for the RG functions defined in terms of the bare coupling constant $\alpha_0$ in the case of using the higher derivative regularization. However, if the RG functions are defined in terms of the renormalized coupling constant, the NSVZ $\beta$-function is obtained only in a special subtraction scheme, which is constructed in this paper for ${\cal
N}=1$ SQED regularized by higher derivatives. Unlike DRED, the higher derivative regularization allows to construct the NSVZ scheme by imposing the boundary conditions on the renormalization constants. The results are verified by an explicit three-loop calculation. In this approximation we also relate the $\overline{\mbox{DR}}$ scheme with the NSVZ scheme obtained with the higher derivatives regularization.

The paper is organized as follows: In Sect. \[Section\_Regularization\] we introduce the higher derivative regularization and remind how the NSVZ $\beta$-function can be obtained using this regularization, if the RG functions are defined in terms of the bare coupling constant and the Pauli–Villars masses are proportional to the parameter $\Lambda$ in the higher derivative term. However, after a rescaling of the Pauli–Villars masses which depends on the bare coupling constant the NSVZ relation for these RG functions is no longer valid. This is discussed in Sect. \[Section\_Rescaling\]. A standard definition of the RG functions is recalled in Sect. \[Section\_Boundary\_Conditions\]. In general, these RG functions (defined in terms of the renormalized coupling constant) do not satisfy the NSVZ relation. However, there is an NSVZ scheme in which the NSVZ relation is valid. In Sect. \[Section\_Boundary\_Conditions\] we demonstrate that this scheme is obtained by imposing a special boundary conditions on the renormalization constants, and the RG functions defined by two different ways coincide in this scheme. Changing the boundary conditions it is possible to change the relation between the $\beta$-function and the anomalous dimension. All these results are verified by an explicit three-loop calculation with the higher derivative regularization in Sect. \[Section\_Explicit\_Three\_Loop\]. In Sect. \[Section\_DRED\_NSVZ\] we compare the results obtained in the $\overline{\mbox{DR}}$-scheme with the results obtained with the higher derivative regularization in the NSVZ scheme. Some technical details are discussed in the Appendexes.

Derivation of the NSVZ $\beta$-function for ${\cal N}=1$ SQED
=============================================================

\[Section\_Regularization\]

In this section we recall how the NSVZ $\beta$-function can be obtained with the higher derivative regularization for ${\cal
N}=1$ SQED. It is convenient to describe this theory in terms of ${\cal N}=1$ superfields [@West:1990tg; @Buchbinder:1998qv]. Then in the massless limit the action of ${\cal N}=1$ SQED is written as

$$\label{Action}
S = \frac{1}{4e_0^2}\mbox{Re}\int d^4x\,d^2\theta\,W^a W_a +
\frac{1}{4} \int d^4x\,d^4\theta\,\Big(\phi^* e^{2V}\phi +
\widetilde\phi^* e^{-2V} \widetilde\phi\Big),$$

where $e_0$ is a bare coupling constant, $\phi$ and $\widetilde\phi$ are (bare) chiral superfields, and $V$ is a real gauge superfield.

In order to regularize this theory by higher derivatives, it is necessary to insert into the first term of Eq. (\[Action\]) a regularizing function $R$, such that $R(0)=1$ and $R(\infty)=\infty$, which contains higher derivatives [@Slavnov:1971aw; @Slavnov:1972sq]:

$$\frac{1}{4e_0^2}\mbox{Re}\int d^4x\,d^2\theta\,W^a W_a
\to\frac{1}{4e_0^2} \mbox{Re}\int d^4x\,d^2\theta W^a
R(\partial^2/\Lambda^2) W_a.$$

For example, it is convenient to choose this function as $R = 1+\partial^{2n}/\Lambda^{2n}$, where $\Lambda$ is a dimensionful parameter. Then the divergences remain only in the one-loop approximation [@Faddeev:1980be]. According to the standard prescription they should be regularized by inserting the Pauli–Villars determinants into the generating functional [@Slavnov:1977zf]:

$$\label{Generating_Functional}
Z[J,j,\widetilde j] = \int DV\, D\phi\,
D\widetilde\phi\,\prod\limits_{I=1}^n
(\det(V,M_I))^{c_I}\exp\Big(i S_{\mbox{\scriptsize reg}} + i
S_{\mbox{\scriptsize gf}} + i S_{\mbox{\scriptsize source}} \Big),$$

where $J$, $j$, and $\widetilde j$ are the sources, and $M_I$ are masses of the Pauli–Villars fields.

Let us describe all parts of this expression:

$$\begin{aligned}
\label{Regularized_Action}
&& S_{\mbox{\scriptsize reg}} = \frac{1}{4e_0^2}\mbox{Re}\int
d^4x\,d^2\theta\,W^a R(\partial^2/\Lambda^2) W_a + \frac{1}{4}
\int d^4x\,d^4\theta\,\Big(\phi^* e^{2V}\phi +
\widetilde\phi^* e^{-2V} \widetilde\phi\Big);\qquad\\
\label{PV_Action} && S_{\mbox{\scriptsize gf}} = - \frac{1}{64
e_0^2}\int d^4x\,d^4\theta\, \Big(V R(\partial^{2}/\Lambda^{2})
D^2 \bar D^2
V + V R(\partial^{2}/\Lambda^{2}) \bar D^2 D^2 V\Big);\\
\label{Sources_Action} && S_{\mbox{\scriptsize source}} = \int
d^4x\,d^4\theta\,V J + \Big(\int d^4x\,d^2\theta\, (\phi j +
\widetilde\phi\, \widetilde j ) +\mbox{c.c.}\Big).\end{aligned}$$

In the Abelian case it is not necessary to introduce ghost (super)fields. The Pauli–Villars determinants can be presented as functional integrals over the corresponding (chiral) Pauli–Villars superfields $\Phi$ and $\widetilde\Phi$:

$$\Big(\det(V,M)\Big)^{-1} = \int
D\Phi\,D\widetilde\Phi\,\exp(iS_{\mbox{\scriptsize PV}}),$$

where the action for the Pauli–Villars superfields is

$$\label{PV_Action1}
S_{\mbox{\scriptsize PV}} = \frac{1}{4}\int d^8x\,\Big(\Phi^*
e^{2V}\Phi + \widetilde\Phi^* e^{-2V}\widetilde\Phi \Big) +
\Big(\frac{1}{2} \int d^4x\,d^2\theta\,M \Phi \widetilde\Phi
+\mbox{c.c.}\Big).$$

For cancellation of remaining one-loop divergences the coefficients $c_I$ should satisfy the conditions [@Slavnov:1977zf]

$$\label{C_Conditions}
\sum\limits_I c_I =1;\qquad \sum\limits_I c_I M_I^2 =0.$$

We choose masses of the Pauli–Villars fields so that they will be proportional to the parameter $\Lambda$ in the higher derivative term [@Stepanyantz:2011wq; @Stepanyantz:2011jy]:

$$\label{Scheme}
M_I = a_I\Lambda,$$

where $a_I$ are arbitrary real constants which do not depend of the bare charge. Then the NSVZ relation is obtained exactly in all orders for the RG functions defined in terms of the bare coupling constant independently of a concrete renormalization prescription. Let us briefly describe, how this can be shown.

Due to the Ward identity the two-point Green function of the gauge superfield is transverse:

$$\label{D_Definition}
\Gamma^{(2)}_{V} - S_{\mbox{\scriptsize gf}} = - \frac{1}{16\pi}
\int
\frac{d^4p}{(2\pi)^4}\,d^4\theta\,V(\theta,-p)\,\partial^2\Pi_{1/2}
V(\theta,p)\, d^{-1}(\alpha_0,\Lambda/p),$$

where the supersymmetric transversal projector is given by

$$\partial^2\Pi_{1/2}  = -\frac{1}{8} D^a \bar D^2 D_a.$$

The function $d$ defined by Eq. (\[D\_Definition\]) coincides with the invariant charge $\alpha_{\mbox{\scriptsize
inv}}$ [@Bogolyubov:1980nc], related with the polarization operator $\Pi(\alpha_0,\Lambda/p)$ by the equation

$$d_0^{-1}(\alpha_0, \Lambda/p) = \frac{1}{\alpha_0} \Big(1
+\frac{\alpha_0}{\pi} \Pi_0(\alpha_0, \Lambda/p)\Big)=\frac{1}
{\alpha_{\mbox{\scriptsize inv}}(\alpha_0,\Lambda/p)}.$$

Similarly, a part of the effective action corresponding to the two-point Green function of the matter superfields can be written in the form

$$\label{Gamma_Matter}
\Gamma^{(2)}_\phi = \frac{1}{4}\int \frac{d^4p}{(2\pi)^4}
d^4\theta\, \Big(\phi^*(\theta,-p)\, \phi(\theta,p) +
\widetilde\phi^*(\theta,-p)\, \widetilde\phi(\theta,p) \Big)
G(\alpha_0,\Lambda/p).$$

Here we use the following renormalization procedure [@Bogolyubov:1980nc; @Collins:1984xc]: we define a renormalized coupling constant $\alpha(\alpha_0,\Lambda/\mu)$, where $\mu$ is a renormalization parameter, requiring that the inverse invariant charge $d^{-1}(\alpha_0(\alpha,\Lambda/\mu),\Lambda/p)$ is finite in the limit $\Lambda\to \infty$ (and, therefore, in this limit depends only on $\alpha$ and $\mu/p$). Then the renormalization constant $Z_3$ is defined according to the equation

$$\frac{1}{\alpha_0} \equiv \frac{Z_3(\alpha,\Lambda/\mu)}{\alpha}.$$

In order to renormalize the two-point Green function of the matter superfields, we construct the renormalization constant $Z$, requiring that the renormalized two-point Green function $Z
G$ is finite in the limit $\Lambda\to \infty$:

$$G_{\mbox{\scriptsize ren}}(\alpha,\mu/p) = \lim\limits_{\Lambda\to
\infty} Z(\alpha, \Lambda/\mu)G(\alpha_0,\Lambda/p).$$

Certainly, the renormalized coupling constant $\alpha$ and the renormalization constant $Z$ are not uniquely defined and depend on a choice of a renormalization scheme [@Vladimirov:1975mx].

The RG functions are usually defined in terms of the renormalized coupling constant. This definition is presented in Sect. \[Section\_Boundary\_Conditions\]. However, it is also possible to define the RG functions in terms of the bare coupling constant $\alpha_0$ according to the following prescription:

$$\begin{aligned}
\label{Beta_Definition1}
&& \beta\Big(\alpha_0(\alpha,\Lambda/\mu)\Big) \equiv \frac{d
\alpha_0(\alpha,\Lambda/\mu)}{d\ln\Lambda}
\Big|_{\alpha=\mbox{\scriptsize const}};\vphantom{\Bigg|}\\
\label{Gamma_Definition1} &&
\gamma\Big(\alpha_0(\alpha,\Lambda/\mu)\Big) \equiv - \frac{d \ln
Z(\alpha,\Lambda/\mu)}{d
\ln\Lambda}\Big|_{\alpha=\mbox{\scriptsize const}},\end{aligned}$$

where $\alpha$ and $\Lambda$ are considered as independent variables. Nevertheless, for finding these functions it is necessary to use the relation between the bare coupling constant $\alpha_0$ and the renormalized coupling constant $\alpha$. Really, the differentiation with respect to $\Lambda$ is made at a fixed value of $\alpha$. Therefore, before this differentiating we should express $\alpha_0$ in terms of $\alpha$ and $\Lambda/\mu$. This implies that expressions dependent on the renormalization scheme appear at intermediate steps of the calculation. However, final expressions for the RG functions (\[Beta\_Definition1\]) and (\[Gamma\_Definition1\]) are independent of a renormalization prescription. In order to see this, let us consider the function

$$d^{-1}(\alpha_0,\Lambda/p) = d^{-1}(\alpha,\mu/p) + \mbox{(terms
vanishing in the limit $p\to 0$)}$$

and differentiate it with respect to $\ln\Lambda$ (at a fixed value of $\alpha$). The derivative of the right hand side evidently vanishes in the limit $p\to 0$. Therefore,

$$\label{Vanishing_Derivative}
0 = \lim\limits_{p\to 0} \frac{d
d^{-1}(\alpha_0,\Lambda/p)}{d\ln\Lambda}
\Big|_{\alpha=\mbox{\scriptsize const}} = \lim\limits_{p\to
0}\Big(\frac{\partial d^{-1}(\alpha_0,\Lambda/p)}{\partial
\alpha_0} \beta(\alpha_0) - \frac{\partial
d^{-1}(\alpha_0,\Lambda/p)}{\partial\ln p} \Big),$$

where in the last equality $\alpha_0$ and $p$ are considered as independent variables. From this equation the $\beta$-function (\[Beta\_Definition1\]) is expressed through the expressions which evidently do not depend on a renormalization prescription. Therefore, this function does not depend on a choice of a renormalization scheme. Similarly, differentiating the equality

$$\ln G(\alpha_0,\Lambda/q) = \ln G_{\mbox{\scriptsize
ren}}(\alpha,\mu/q) - Z(\alpha,\Lambda/\mu) + \mbox{(terms
vanishing in the limit $q\to 0$)}$$

with respect to $\ln\Lambda$ at a fixed value of $\alpha$, in the limit $q\to 0$ we obtain

$$\gamma(\alpha_0) = \lim\limits_{q\to 0} \Big(\frac{\partial\ln
G(\alpha_0,\Lambda/q)}{\partial \alpha_0} \beta(\alpha_0) -
\frac{\partial\ln G(\alpha_0,\Lambda/q)}{\partial \ln q}\Big).$$

Therefore, the anomalous dimension (\[Gamma\_Definition1\]) does not depend on a choice of the function $\alpha_0(\alpha,\Lambda/\mu)$ as well.

In order to find an expression for the $\beta$-function (\[Beta\_Definition1\]) in the case of ${\cal N}=1$ SQED regularized by higher derivatives, it is possible to make the substitution

$$\label{V_Substitution}
V \to \bar \theta^a \bar\theta_a \theta^b \theta_b\equiv \theta^4$$

in the expression

$$\Delta\Gamma^{(2)}_V = \Gamma^{(2)}_V - S_{\mbox{\scriptsize gf}}
- S.$$

Then using Eq. (\[D\_Definition\]) we obtain

$$\begin{aligned}
\label{Calculate2}
&& \frac{d(\Delta\Gamma^{(2)}_{\mbox{\scriptsize
V})}}{d\ln\Lambda} \Big|_{V(x,\theta)=\theta^4} =
(2\pi)^3\delta^{4}(p) \frac{d}{d\ln \Lambda}\,
\Big(d^{-1}(\alpha_0,\Lambda/p)-\alpha_0^{-1}\Big)\nonumber\\
&&\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad = -
(2\pi)^3\delta^{4}(p) \frac{d\alpha_0^{-1}}{d \ln \Lambda} =
(2\pi)^3\delta^{4}(p) \frac{\beta(\alpha_0)}{\alpha_0^2},\qquad\end{aligned}$$

where (exactly as in Eqs. (\[Beta\_Definition1\]) and (\[Gamma\_Definition1\])) $\alpha$ and $\Lambda$ are considered as independent variables, and we take into account the first equality in Eq. (\[Vanishing\_Derivative\]).[^2]

The crucial observation is that if the higher derivative method is used for a regularization, the integrals obtained by calculating the left hand side of Eq. (\[Calculate2\]) are integrals of (double) total derivatives [@Soloshenko:2003nc; @Smilga:2004zr; @Stepanyantz:2011jy]. Therefore, one of the loop integrals can be calculated analytically. Using a rather complicated technique [@Stepanyantz:2011jy], it is possible to obtain

$$\label{Green_Function_Relation}
\frac{\beta(\alpha_0)}{\alpha_0^2} = \frac{d}{d\ln \Lambda}\,
\Big(d^{-1}(\alpha_0,\Lambda/p)-\alpha_0^{-1}\Big)\Big|_{p=0} =
\frac{1}{\pi}\Big(1-\frac{d}{d\ln\Lambda} \ln
G(\alpha_0,\Lambda/q)\Big|_{q=0}\Big)$$

exactly in all orders. The right hand side of this equation can be expressed through the anomalous dimension (\[Gamma\_Definition1\]):

$$\begin{aligned}
\label{NSVZ_Bare}
&& \frac{1}{\pi}\Big(1-\frac{d}{d\ln\Lambda} \ln
G(\alpha_0,\Lambda/q)\Big|_{q=0}\Big) = \frac{1}{\pi}-
\frac{1}{\pi}\frac{d}{d\ln\Lambda} \Big(\ln G_{\mbox{\scriptsize
ren}}(\alpha,\mu/q) - \ln
Z(\alpha,\Lambda/\mu) \Big)\Big|_{q=0}\nonumber\\
&& = \frac{1}{\pi} \Big(1 -
\gamma\Big(\alpha_0(\alpha,\Lambda/\mu)\Big)\Big),\end{aligned}$$

where we take into account that the function $Z G$ does not depend on $\Lambda$ in the limit $q\to 0$. Thus, we obtain

$$\frac{\beta(\alpha_0)}{\alpha_0^2}= \frac{1}{\pi} \Big(1 -
\gamma\Big(\alpha_0(\alpha,\Lambda/\mu)\Big)\Big).$$

This relation gives the NSVZ $\beta$-function (\[NSVZ\]), which was first obtained in [@Vainshtein:1986ja; @Shifman:1985fi] using a different method.

It is important that this relation between the RG functions defined in terms of the [*bare*]{} coupling constant is valid for an arbitrary choice of the functions $\alpha_0(\alpha,\Lambda/\mu)$ and $Z(\alpha,\Lambda/\mu)$, which specify a subtraction scheme. (Note that at intermediate steps these functions are needed for calculating the RG functions (\[Beta\_Definition1\]) and (\[Gamma\_Definition1\]).) Thus, with the higher derivative regularization the NSVZ $\beta$-function is naturally obtained in terms of the bare coupling constant.

Regularization dependence of the RG functions defined in terms of the bare charge
=================================================================================

\[Section\_Rescaling\]

The RG functions (\[Beta\_Definition1\]) and (\[Gamma\_Definition1\]) are related by Eq. (\[NSVZ\]) independently of a particular choice of counterterms, if the theory is regularized according to the prescription described in Sect. \[Section\_Regularization\]. However, it is possible to change the higher derivative regularization so that to modify the NSVZ relation. This can be done, for example, by a bare coupling constant dependent rescaling of the Pauli–Villars masses. Let us investigate, how the NSVZ relation for bare quantities is modified after this rescaling.

Actually, the regularization for which the NSVZ relation is obtained for the RG functions (\[Beta\_Definition1\]) and (\[Gamma\_Definition1\]) is singled out by the condition (\[Scheme\]). More exactly, the coefficients $a_I$ should not depend on $\alpha_0$. However, one-loop divergences can be also cancelled, if we choose

$$\label{MZ}
M_I = a_I \Lambda/z_0(\alpha_0),$$

where $z(\alpha_0)$ is an arbitrary finite function of the bare coupling constant $\alpha_0 = e_0^2/4\pi$. This is equivalent to the redefinition of the $S_{\mbox{\scriptsize reg}}$ and $S_{\mbox{\scriptsize PV}}$:

$$\begin{aligned}
\label{Z_Regularized_Action}
&& S_{\mbox{\scriptsize reg}} = \frac{1}{4e_0^2}\mbox{Re}\int
d^4x\,d^2\theta\,W^a W_a + \frac{1}{4}z_0(\alpha_0) \int
d^4x\,d^4\theta\,\Big(\phi^* e^{2V}\phi + \widetilde\phi^* e^{-2V}
\widetilde\phi\Big);\\
\label{Z_PV_Action} && S_{\mbox{\scriptsize PV}} =
\frac{1}{4}z_0(\alpha_0) \int d^8x\,\Big(\Phi^* e^{2V}\Phi +
\widetilde\Phi^* e^{-2V}\widetilde\Phi \Big) + \Big(\frac{1}{2}
\int d^4x\,d^2\theta\,M \Phi \widetilde\Phi
+\mbox{c.c.}\Big).\qquad\end{aligned}$$

These actions are obtained after rescaling of the intergration variables

$$\phi\to \sqrt{z_0(\alpha_0)}\,\phi;\qquad \widetilde\phi\to
\sqrt{z_0(\alpha_0)}\, \widetilde\phi;\qquad \Phi_I \to
\sqrt{z_0(\alpha_0)}\,\Phi_I;\qquad \widetilde\Phi_I \to
\sqrt{z_0(\alpha_0)}\, \widetilde\Phi_I.$$

in the generating functional (\[Generating\_Functional\]). In order to avoid anomalous contributions of the integration measure, discussed in [@ArkaniHamed:1997mj], the Pauli–Villars fields should be rescaled in the same way as the usual fields. Then $S_{\mbox{\scriptsize reg}}$ given by Eq. (\[Z\_Regularized\_Action\]) and $S_{\mbox{\scriptsize PV}}$ given by Eq. (\[Z\_PV\_Action\]) produce diagrams with insertions of

$$\frac{1}{4}(z_0(\alpha_0)-1) \int d^4x\,d^4\theta\,\Big(\phi^*
e^{2V}\phi + \widetilde\phi^* e^{-2V} \widetilde\phi\Big)$$

on the matter lines and similar diagrams with closed loops of the Pauli–Villars fields. In order to calculate a contribution of these diagrams, we note that Eq. (\[Green\_Function\_Relation\]) is also valid for all diagrams beyond the one-loop approximation and for all values of $z_0$. This can be proved repeating the calculations made in [@Stepanyantz:2011jy] taking into account that dependence on the Pauli–Villars masses is the same in both sides of this equation. However, a one-loop contribution to the right hand side of Eq. (\[Green\_Function\_Relation\]) is

$$\label{35} \frac{1}{\pi}\cdot \frac{d}{d\ln\Lambda}
\sum\limits_{I=1}^n c_I \ln \frac{M_I}{p}\Big|_{p=0}.$$

If the Pauli–Villars masses are given by Eq. (\[MZ\]), instead of Eq. (\[35\]) we obtain

$$\frac{1}{\pi} \sum\limits_{I=1}^n c_I \Big(1-\frac{d\ln
z_0}{d\ln\Lambda}\Big) = \frac{1}{\pi} \Big(1-\frac{d\ln
z_0}{d\ln\Lambda}\Big) = \frac{1}{\pi} \Big(1-\beta(\alpha_0)
\frac{d\ln z_0}{d\alpha_0}\Big).$$

For $z_0=1$ this gives $1/\pi$ in the right hand side of Eq. (\[Green\_Function\_Relation\]). Now it is evident that for an arbitrary function $z_0$ the generalization of Eq. (\[Green\_Function\_Relation\]) can be written as

$$\frac{d}{d\ln\Lambda} \Big(d^{-1}(\alpha_0,\Lambda/p) -
\alpha_0^{-1}\Big)\Big|_{p=0}  =
\frac{1}{\pi}\Big(1-\frac{d}{d\ln\Lambda} \ln
G(\alpha_0,\Lambda/q)\Big|_{q=0} - \frac{d\ln
z_0(\alpha_0)}{d\ln\Lambda} \Big).$$

From this equation we easily obtain that the $\beta$-function (\[Beta\_Definition1\]), which depends on the bare coupling constant $\alpha_0$, has the form

$$\beta(\alpha_0) = \frac{\alpha_0^2}{\pi} \cdot \frac{1 -
\gamma(\alpha_0)}{1 + \alpha_0^2\, (d\ln z_0/d\alpha_0)/\pi}.$$

Therefore, changing the Pauli–Villars masses according to Eq. (\[MZ\]) we modify the expression for the exact $\beta$-function expressed in terms of the bare coupling constant.

Also the NSVZ relation can be modified by a finite tuning of the bare charge $\alpha_0 \to \alpha_0'(\alpha_0)$. Combining this finite tuning with the rescaling of the Pauli–Villars masses, it is possible to obtain the RG functions dependent on the bare charge with coefficients equal to those obtained within the $\overline{\mbox{DR}}$ approach.

NSVZ scheme
===========

\[Section\_Boundary\_Conditions\]

Although the above prescriptions allow to obtain the NSVZ $\beta$-function in terms of the bare coupling constant, it is possible to study a problem whether the NSVZ relation is valid for the RG functions defined in terms of the renormalized coupling constant, if the theory is regularized by higher derivatives.

Let us consider the RG functions, defined in terms of the renormalized coupling constant:

$$\begin{aligned}
\label{Beta_Definition2} &&
\widetilde\beta\Big(\alpha(\alpha_0,\Lambda/\mu)\Big) \equiv
\frac{d\alpha(\alpha_0,\Lambda/\mu)}{d\ln\mu}\Big|_{\alpha_0=\mbox{\scriptsize const}};\\
\label{Gamma_Definition2} &&
\widetilde\gamma\Big(\alpha(\alpha_0,\Lambda/\mu)\Big) \equiv
\frac{d}{d\ln\mu}\ln
ZG(\alpha_0,\Lambda/\mu)\Big|_{\alpha_0=\mbox{\scriptsize const}}
= \frac{d\ln Z(\alpha(\alpha_0,\Lambda/\mu),
\Lambda/\mu)}{d\ln\mu}\Big|_{\alpha_0=\mbox{\scriptsize
const}},\qquad\end{aligned}$$

where $\alpha_0$ and $\mu$ are considered as independent variables. These functions depend on the arbitrariness of choosing $\alpha$ and $Z$. As a consequence, in general, these functions do not satisfy the NSVZ relation, which was originally derived for the bare quantities.

The arbitrariness of choosing $\alpha$ and $Z$ can be fixed by imposing certain boundary conditions. For example, the momentum subtraction (MOM) scheme is defined by the following requirements:

$$\label{MOM_Boundary_Conditions}
Z_{\mbox{\scriptsize MOM}}G(\alpha_{\mbox{\scriptsize MOM}},p=\mu)
= 1; \qquad d^{-1}(\alpha_{\mbox{\scriptsize MOM}},p=\mu)
=\alpha_{\mbox{\scriptsize MOM}}^{-1}.$$

In terms of the polarization operator the second condition can be rewritten as $\Pi(\alpha_{\mbox{\scriptsize MOM}},p=\mu) = 0$. In the MOM scheme the $\beta$-function (\[Beta\_Definition2\]) coincides with the Gell-Mann–Low function $\Psi(\alpha)$ (see [@Gorishnii:1990kd] for the detailed explanation).

The other schemes can be obtained by making finite renormalizations:

$$\label{Finite_Renormalization}
\frac{1}{\alpha} = \frac{z_3(\alpha_{\mbox{\scriptsize
MOM}})}{\alpha_{\mbox{\scriptsize MOM}}};\qquad
Z(\alpha,\Lambda/\mu) = z(\alpha_{\mbox{\scriptsize MOM}})\,
Z_{\mbox{\scriptsize MOM}}(\alpha_{\mbox{\scriptsize
MOM}},\Lambda/\mu),$$

where $z_3$ and $z$ are finite functions. In the three-loop approximation such a renormalization relating the NSVZ scheme and the MOM scheme is constructed in Appendix \[Appendix\_Renormalization\].

In this section we construct a scheme of the ultraviolet renormalization for which the NSVZ relation is also satisfied for the renormalized functions in case of using the higher derivative regularization. This scheme is related with the MOM scheme by the finite renormalization (\[Finite\_Renormalization\]). It is constructed by imposing an additional condition, which can be formulated as follows: There should be a point $x_0
=\ln\Lambda/\mu_0$ such that

$$\label{NSVZ_Scheme}
\alpha_0(\alpha_{\mbox{\scriptsize NSVZ}},x_0) =
\alpha_{\mbox{\scriptsize NSVZ}}; \qquad Z_{\mbox{\scriptsize
NSVZ}}(\alpha_{\mbox{\scriptsize NSVZ}},x_0)=1$$

for all values of $\alpha_{\mbox{\scriptsize NSVZ}}$. Equivalently, there is a point $x_0$ such that

$$\label{NSVZ_Scheme_Z3}
(Z_3)_{\mbox{\scriptsize NSVZ}}(\alpha_{\mbox{\scriptsize
NSVZ}},x_0) = 1;\qquad Z_{\mbox{\scriptsize
NSVZ}}(\alpha_{\mbox{\scriptsize NSVZ}},x_0)=1.$$

Although these conditions are similar to the requirement $Z_3=1$, which was used in Ref. [@Kataev:2013vua] to define the conformally invariant limit of perturbative QED [@Kataev:2010tm], Eqs. (\[NSVZ\_Scheme\]) and (\[NSVZ\_Scheme\_Z3\]) are imposed in a single point $x_0$ only. In Appendix \[Appendix\_Existence\] we prove that the scheme defined by Eq. (\[NSVZ\_Scheme\]) exists and without loss of generality it is possible to choose $x_0=0$. This result is also verified by an explicit calculation in the three-loop approximation.

Let us prove this statement. We assume that Eq. (\[NSVZ\_Scheme\]) is satisfied. It is convenient to define a variable $x = \ln\Lambda/\mu$. Then the RG functions (\[Beta\_Definition2\]) and (\[Gamma\_Definition2\]) can be written as

$$\begin{aligned}
&& \widetilde\beta\left(\alpha(\alpha_0,x)\right) =
-\frac{\partial
\alpha(\alpha_0,x)}{\partial x}; \nonumber\\
&& \widetilde\gamma\left(\alpha(\alpha_0,x)\right) = -\frac{d \ln
Z\left(\alpha(\alpha_0,x),x\right)}{d x} = - \frac{\partial \ln
Z(\alpha,x)}{\partial \alpha}\cdot
\frac{\partial\alpha(\alpha_0,x)}{\partial x} - \frac{\partial \ln
Z\left(\alpha(\alpha_0,x),x\right)}{\partial x},\qquad\end{aligned}$$

where the total derivative with respect to $x$ also acts on $x$ inside $\alpha$. Calculating these expressions at the point $x=x_0$ and taking into account Eq. (\[NSVZ\_Scheme\]) we obtain

$$\label{Equality}
\widetilde\beta(\alpha_0) = \beta(\alpha_0);\qquad
\widetilde\gamma(\alpha_0) = \gamma(\alpha_0).$$

Really, in order to prove the first equality we note that

$$\begin{aligned}
&& \alpha_0(\alpha,x) = \alpha + \beta(\alpha_0) (x-x_0) +
O\left((x-x_0)^2\right);\nonumber\\
&& \alpha(\alpha_0,x) = \alpha_0 - \widetilde\beta(\alpha) (x-x_0)
+ O\left((x-x_0)^2\right)\nonumber\\
&&\qquad\qquad\qquad\qquad\qquad = \alpha_0 -
\widetilde\beta(\alpha_0) (x-x_0) + O\left((x-x_0)^2\right).\qquad\end{aligned}$$

Similarly, the second equation in (\[Equality\]) can be obtained taking into account that

$$\frac{\partial \ln Z(\alpha,x_0)}{\partial \alpha} = 0$$

due to Eq. (\[NSVZ\_Scheme\]). According to [@Stepanyantz:2011jy] the functions (\[Beta\_Definition1\]) and (\[Gamma\_Definition1\]) satisfy the equation

$$\beta(\alpha_0) =
\frac{\alpha_0^2}{\pi}\Big(1-\gamma(\alpha_0)\Big),$$

if one uses the regularization described in Sect. \[Section\_Regularization\]. Therefore, from Eq. (\[Equality\]) we conclude that

$$\widetilde\beta(\alpha) =
\frac{\alpha^2}{\pi}\Big(1-\widetilde\gamma(\alpha)\Big),$$

if the boundary conditions (\[NSVZ\_Scheme\]) are imposed. In the three-loop approximation this is verified by an explicit calculation in Sect. \[Section\_Explicit\_Three\_Loop\].

The scheme dependence in the three-loop approximation
=====================================================

\[Section\_Explicit\_Three\_Loop\]

Let us reconsider the results of the three-loop calculation [@Soloshenko:2003nc] for ${\cal N}=1$ SQED and investigate the scheme dependence of the RG functions in this approximation. (The details of the three-loop calculation are described in Appendix \[Appendix\_Renormalization\].) We will start with the two-point Green function of the matter superfields. If the Pauli–Villars masses are chosen according to the prescription (\[Scheme\]), in the two-loop approximation this function is given by the following expression [@Soloshenko:2003sx]:[^3]

$$\begin{aligned}
\label{G(alpha0)}
&& G(\alpha_0,\Lambda/p) = 1 - \int \frac{d^4k}{(2\pi)^4}
\frac{2e_0^2}{k^2 R_k (k+p)^2} + \int \frac{d^4k}{(2\pi)^4}
\frac{d^4l}{(2\pi)^4} \frac{4e_0^4}{k^2 R_k l^2 R_l}
\Bigg(\frac{1}{(k+p)^2 (l+p)^2}\nonumber\\
&& + \frac{1}{(l+p)^2 (k+l+p)^2} - \frac{(k+l+2p)^2}{(k+p)^2
(l+p)^2 (k+l+p)^2} \Bigg) + \int \frac{d^4k}{(2\pi)^4}
\frac{d^4l}{(2\pi)^4} \frac{4 e_0^4}{k^2 R_k^2 (k+p)^2}
\nonumber\\
&& \times \Bigg(\frac{1}{l^2(k+l)^2} - \sum\limits_{I=1}^n c_I
\frac{1}{\left(l^2 + M_I^2\right)\left((k+l)^2+M_I^2\right)}\Bigg)
+ O(e_0^6) ,\end{aligned}$$

where $R_k\equiv R(k^2/\Lambda^2)$. (This expression is written in the Euclidean space after the Wick rotation.)

We split the bare coupling constant $\alpha_0 = e_0^2/4\pi$ into the renormalized coupling constant $\alpha$ and a counterterm. In the lowest approximation

$$\label{One-Loop_Coupling_Renormalization}
\frac{1}{\alpha_0} = \frac{1}{\alpha} - \frac{1}{\pi} \Big(
\ln\frac{\Lambda}{\mu} + b_1\Big) + O(\alpha),$$

where the logarithm compensates a one-loop divergence and a finite parameter $b_1$ is not fixed. Different choices of $b_1$ correspond to different subtraction schemes. Equivalently, Eq. (\[One-Loop\_Coupling\_Renormalization\]) can be rewritten in the form

$$\alpha_0 = \alpha\Big(1 + \frac{\alpha}{\pi}
\Big(\ln\frac{\Lambda}{\mu} + b_1\Big) + O(\alpha^2) \Big).$$

In order to cancel divergences in the two-point Green function of the matter superfields, the function $G(\alpha_0,\Lambda/p)$ should be multiplied by the renormalization constant $Z$. This corresponds to the following redefinition of the matter superfields:

$$\label{Rescaling_Of_Arguments}
\phi = Z^{1/2}(\alpha,\Lambda/\mu)\,\phi_{\mbox{\scriptsize
ren}};\qquad \widetilde \phi =
Z^{1/2}(\alpha,\Lambda/\mu)\,\widetilde\phi_{\mbox{\scriptsize
ren}},$$

where $\phi$ and $\widetilde\phi$ are the bare fields, and $\phi_{\mbox{\scriptsize ren}}$ and $\widetilde\phi_{\mbox{\scriptsize ren}}$ are the renormalized fields. It is easy to see that in the one-loop approximation the renormalization constant $Z$, which is chosen so that the renormalized Green function $Z G$ should be finite, is

$$Z = 1 + \frac{\alpha}{\pi}\Big(\ln\frac{\Lambda}{\mu}+g_1\Big) +
O(\alpha^2),$$

where $g_1$ is an arbitrary finite constant. In principle, without loss of generality the constant $g_1$ can be excluded by a redefinition of the parameter $\mu$. However, we will not do this, because in the next orders similar coefficients cannot be excluded by this way. Choosing a higher derivative term in the momentum space as

$$\label{Explicit_Regulator}
R_k = 1+\frac{k^{2n}}{\Lambda^{2n}},$$

it is possible to find the following explicit expression for the constant $Z$ in the next order [@Soloshenko:2003sx]:

$$\label{Two_Loop_Z}
Z = 1 + \frac{\alpha}{\pi}\Big(\ln\frac{\Lambda}{\mu}+g_1\Big)
+\frac{\alpha^2}{\pi^2}\ln^2\frac{\Lambda}{\mu}
-\frac{\alpha^2}{\pi^2}\ln\frac{\Lambda}{\mu}\Big(\sum\limits_{I=1}^n
c_I\ln a_I - b_1 + \frac{3}{2} - g_1\Big) + \frac{\alpha^2
g_2}{\pi^2} + O(\alpha^3),$$

where

$$a_I = M_I/\Lambda$$

and $g_2$ is another finite constant. Choosing $g_2$ we partially fix the subtraction scheme in the considered approximation. The subtraction scheme in the considered approximation is completely fixed by choosing $g_1$, $g_2$, and the finite constants in the renormalization of the coupling constant, $b_1$ and $b_2$ (defined below).

Now, let us calculate the anomalous dimension according to Eqs. (\[Gamma\_Definition1\]) and (\[Gamma\_Definition2\]).

In order to find the anomalous dimension according to the prescription (\[Gamma\_Definition1\]), we should differentiate

$$\begin{aligned}
&& \ln Z(\alpha,\Lambda/\mu) =
\frac{\alpha}{\pi}\Big(\ln\frac{\Lambda}{\mu}+g_1\Big)
+\frac{\alpha^2}{2\pi^2}\ln^2\frac{\Lambda}{\mu} \nonumber\\
&&\qquad\qquad\qquad\qquad\quad\ \
-\frac{\alpha^2}{\pi^2}\ln\frac{\Lambda}{\mu}\Big(\sum\limits_{I=1}^n
c_I\ln a_I - b_1 + \frac{3}{2}\Big) + \frac{\alpha^2}{\pi^2}
\Big(g_2-\frac{1}{2}g_1^2\Big) + O(\alpha^3)\qquad\end{aligned}$$

with respect to $\ln\Lambda$ and write the result in terms of $\alpha_0$. Then we obtain

$$\begin{aligned}
\label{Gamma_Answer1}
&& \gamma(\alpha_0) = - \frac{d\ln Z}{d\ln\Lambda} =
-\frac{\alpha}{\pi} -
\frac{\alpha^2}{\pi^2}\Big(\ln\frac{\Lambda}{\mu} -
\sum\limits_{I=1}^n c_I \ln a_I + b_1 -
\frac{3}{2}\Big) + O(\alpha^3)\nonumber\\
&&\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad =
-\frac{\alpha_0}{\pi} +
\frac{\alpha_0^2}{\pi^2}\Big(\frac{3}{2}+\sum\limits_{I=1}^n c_I
\ln a_I \Big) + O(\alpha_0^3).\qquad\end{aligned}$$

This expression is independent of the finite constants $g_i$ and $b_i$, which fix the subtraction scheme.

The anomalous dimension $\widetilde\gamma(\alpha)$ defined by Eq. (\[Gamma\_Definition2\]) can be constructed similarly. For this purpose we rewrite $\ln Z$ in terms of $\alpha_0$ using Eq. (\[One-Loop\_Coupling\_Renormalization\]):

$$\begin{aligned}
&& \ln Z(\alpha(\alpha_0,\Lambda/\mu),\Lambda/\mu) =
\frac{\alpha_0}{\pi}\Big(\ln\frac{\Lambda}{\mu}+g_1\Big)
-\frac{\alpha_0^2}{2\pi^2}\ln^2\frac{\Lambda}{\mu} \nonumber\\
&&\qquad\qquad\qquad
-\frac{\alpha_0^2}{\pi^2}\ln\frac{\Lambda}{\mu}\Big(\sum\limits_{I=1}^n
c_I\ln a_I +g_1 + \frac{3}{2}\Big) + \frac{\alpha_0^2}{\pi^2}
\Big(g_2-\frac{1}{2}g_1^2 - g_1 b_1\Big) + O(\alpha_0^3).\qquad\end{aligned}$$

Differentiating this expression with respect to $\ln\mu$ and writing the result in terms of $\alpha$ we obtain

$$\begin{aligned}
\label{Gamma_Answer2}
&& \widetilde\gamma(\alpha) = \frac{d\ln Z}{d\ln\mu} =
-\frac{\alpha_0}{\pi} +
\frac{\alpha_0^2}{\pi^2}\Big(\ln\frac{\Lambda}{\mu} +
\sum\limits_{I=1}^n c_I \ln a_I + g_1 +
\frac{3}{2}\Big) + O(\alpha_0^3)\nonumber\\
&&\qquad\qquad\qquad\qquad\qquad\qquad\quad  = -\frac{\alpha}{\pi}
+ \frac{\alpha^2}{\pi^2}\Big(\frac{3}{2}+\sum\limits_{I=1}^n c_I
\ln a_I  - b_1 + g_1\Big) + O(\alpha^3).\qquad\end{aligned}$$

Unlike Eq. (\[Gamma\_Answer1\]), this expression depends on the constants $g_1$ and $b_1$.

Let us study, under what conditions the NSVZ relation is valid. For this purpose we write the explicit expression for the three-loop renormalized coupling constant. The integrals defining the three-loop $\beta$-function are presented, for example, in Ref. [@Stepanyantz:2011jy]. They were obtained in [@Soloshenko:2003nc], but a more compact form for writing them was proposed later. If the higher derivative term is given by Eq. (\[Explicit\_Regulator\]), these integrals can be calculated explicitly using the results of [@Soloshenko:2003sx]. After this calculation we split the bare coupling constant into the renormalized coupling constant and a counterterm in such a way that the function $d^{-1}$ is finite. It is possible to find (see Appendix \[Appendix\_Renormalization\]) that the result can be written in the following form:

$$\begin{aligned}
\label{Alpha0}
&& \frac{1}{\alpha_0} = \frac{1}{\alpha} - \frac{1}{\pi} \Big(
\ln\frac{\Lambda}{\mu} + b_1\Big) - \frac{\alpha}{\pi^2}
\Big(\ln\frac{\Lambda}{\mu} + b_2\Big)\nonumber\\
&& \qquad\qquad\qquad -
\frac{\alpha^2}{\pi^3}\Big(\frac{1}{2}\ln^2 \frac{\Lambda}{\mu}
-\ln \frac{\Lambda}{\mu} \sum\limits_{I=1}^n c_I \ln a_I
-\frac{3}{2} \ln \frac{\Lambda}{\mu} + b_1 \ln\frac{\Lambda}{\mu}
+ b_3\Big) + O(\alpha^3),\qquad\end{aligned}$$

where $b_1$, $b_2$, and $b_3$ are finite constants, which define the subtraction scheme. (This three-loop expression coincides with Eq. (\[One-Loop\_Coupling\_Renormalization\]) at the one-loop level.)

Differentiating Eq. (\[Alpha0\]) with respect to $\ln\Lambda$ and writing the result in terms of $\alpha_0$ we obtain the $\beta$-function defined by Eq. (\[Beta\_Definition1\]):

$$\begin{aligned}
\label{Beta_Expression}
&& \frac{\beta(\alpha_0)}{\alpha_0^2} = \frac{1}{\pi} +
\frac{\alpha}{\pi^2}  + \frac{\alpha^2}{\pi^3}\Big(\ln
\frac{\Lambda}{\mu} - \sum\limits_{I=1}^n c_I \ln a_I -\frac{3}{2}
+ b_1\Big) +
O(\alpha^3)\nonumber\\
&& \qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad =
\frac{1}{\pi} + \frac{\alpha_0}{\pi^2} -
\frac{\alpha_0^2}{\pi^3}\Big( \sum\limits_{I=1}^n c_I \ln a_I
+\frac{3}{2} \Big) + O(\alpha_0^3).\qquad\end{aligned}$$

Thus, this $\beta$-function does not depend on the finite constants $g_i$ and $b_i$, as was proved in Eq. (\[Vanishing\_Derivative\]). Moreover, comparing this expression with Eq. (\[Gamma\_Answer1\]) we see that for all values of $g_i$ and $b_i$ the NSVZ relation is valid.

In order to calculate the $\beta$-function defined by Eq. (\[Beta\_Definition2\]), we re-express $\alpha$ in terms of $\alpha_0$:

$$\begin{aligned}
\label{Alpha}
&& \frac{1}{\alpha} = \frac{1}{\alpha_0} + \frac{1}{\pi} \Big(
\ln\frac{\Lambda}{\mu} + b_1\Big) + \frac{\alpha_0}{\pi^2}
\Big(\ln\frac{\Lambda}{\mu} + b_2\Big)\nonumber\\
&& \qquad + \frac{\alpha_0^2}{\pi^3}\Big(-\frac{1}{2}\ln^2
\frac{\Lambda}{\mu} -\ln \frac{\Lambda}{\mu} \sum\limits_{I=1}^n
c_I \ln a_I -\frac{3}{2} \ln \frac{\Lambda}{\mu}
-b_2\ln\frac{\Lambda}{\mu} - b_1 b_2 + b_3\Big) +
O(\alpha_0^3).\qquad\end{aligned}$$

Differentiating this expression with respect to $\ln
\mu$ and writing the result in terms of $\alpha$, we obtain

$$\begin{aligned}
\label{Beta_Answer2}
&& \frac{\widetilde\beta(\alpha)}{\alpha^2} = \frac{1}{\pi} +
\frac{\alpha_0}{\pi^2}  + \frac{\alpha_0^2}{\pi^3}\Big(-\ln
\frac{\Lambda}{\mu} - \sum\limits_{I=1}^n c_I \ln a_I -\frac{3}{2}
-b_2\Big) +
O(\alpha_0^3)\nonumber\\
&& \qquad\qquad\qquad\qquad\qquad\quad\ \ = \frac{1}{\pi} +
\frac{\alpha}{\pi^2} - \frac{\alpha^2}{\pi^3}\Big(
\sum\limits_{I=1}^n c_I \ln a_I +\frac{3}{2} - b_1 + b_2\Big) +
O(\alpha^3).\qquad\end{aligned}$$

In this case the $\beta$-function is related with the anomalous dimension of the matter superfields (\[Gamma\_Answer2\]) by Eq. (\[NSVZ\]) only if finite constants defining a subtraction scheme satisfy the condition

$$b_2 = g_1.$$

It is easy to see that this identity follows from the system of equations

$$\label{Three_Loop_Scheme}
\alpha_0 = \alpha(\alpha_0,\Lambda/\mu_0) + O(\alpha_0^4);\qquad 1
= Z(\alpha_0,\Lambda/\mu_0) + O(\alpha_0^2),$$

which in the considered approximation coincides with Eq. (\[NSVZ\_Scheme\]). Indeed, from the second equation we get

$$\ln\frac{\Lambda}{\mu_0} = -g_1.$$

Substituting this relation to the first equation of (\[Three\_Loop\_Scheme\]) one can see that $b_1 = g_1$ and $b_2 =
g_1$. Although the constant $b_1$ is also fixed by Eq. (\[Three\_Loop\_Scheme\]), it is evident that the relation between the $\beta$-function and the anomalous dimension is independent of this constant. For $b_1=g_1$, $b_2=g_1$ in the considered approximation $\gamma(\alpha) = \widetilde\gamma(\alpha)$ and $\beta(\alpha) = \widetilde\beta(\alpha)$. This confirms the results obtained in Sect. \[Section\_Boundary\_Conditions\]. In higher orders it is also necessary to consider terms with larger degrees of $\alpha_0$. For example, in the four-loop approximation we should also consider terms proportional to $\alpha_0^4$ in the first equation of (\[NSVZ\_Scheme\]) and terms proportional to $\alpha_0^2$ in the second equation of (\[NSVZ\_Scheme\]). This allows to obtain the next coefficients $b_3$ and $g_2$. If we set $g_1=0$, then the NSVZ scheme corresponds to $g_2=b_1=b_2=b_3=0$.

Relation between the NSVZ and $\overline{\mbox{DR}}$ schemes in the three-loop approximation {#Section_DRED_NSVZ}
============================================================================================

The three-loop $\beta$-function and the two-loop anomalous dimension for an arbitrary ${\cal N}=1$ SYM theory with matter in the $\overline{\mbox{DR}}$ scheme have been calculated in [@Jack:1996vg].[^4] The result does not satisfy the NSVZ relation. However, in [@Jack:1996vg] the authors obtained a special redefinition of the coupling constant after which the NSVZ relation is satisfied. For ${\cal N}=1$ SQED, considered here, this redefinition is written as

$$\label{Jones_Result}
\alpha\equiv \alpha_{\mbox{\scriptsize NSVZ}} \to \alpha'\equiv
\alpha_{\overline{\mbox{\scriptsize{DR}}}} = \alpha -
\frac{\alpha^3}{4\pi^2} + O(\alpha^4).$$

In general, under a finite renormalization

$$\alpha \to \alpha'(\alpha);\qquad Z'(\alpha',\Lambda/\mu) =
z(\alpha) Z(\alpha,\Lambda/\mu)$$

the $\beta$-function (\[Beta\_Definition2\]) and the anomalous dimension (\[Gamma\_Definition2\]) are changed according to the following rules:

$$\begin{aligned}
\label{Beta_Transformation}
&& \widetilde \beta'(\alpha') =
\frac{d\alpha'}{d\ln\mu}\Big|_{\alpha_0=\mbox{\scriptsize const}}
= \frac{d\alpha'}{d\alpha} \widetilde
\beta(\alpha);\\
\label{Gamma_Transformation} && \widetilde \gamma'(\alpha') =
\frac{d\ln Z'}{d\ln\mu}\Big|_{\alpha_0=\mbox{\scriptsize const}} =
\frac{d\ln z}{d\alpha}\cdot \widetilde \beta(\alpha) +
\widetilde\gamma(\alpha).\end{aligned}$$

Using these equations it is easy to see that if $\widetilde\beta(\alpha)$ and $\widetilde\gamma(\alpha)$ satisfy the NSVZ relation, then

$$\widetilde\beta'(\alpha') = \frac{d\alpha'}{d\alpha}\cdot
\frac{\alpha^2}{\pi}
\frac{1-\widetilde\gamma'(\alpha')}{1-\alpha^2 (d\ln z/d
\alpha)/\pi}\,\Big|_{\alpha=\alpha(\alpha')}.$$

&gt;From Eq. (\[Beta\_Transformation\]) we see that a shift $\alpha\to \alpha+\delta\alpha$ leads to the following change of the $\beta$-function:

$$\delta\widetilde\beta(\alpha) =
\widetilde\beta(\alpha)\frac{\partial
\delta\alpha}{\partial\alpha} - \delta\alpha
\frac{\partial\widetilde\beta}{\partial\alpha}.$$

Substituting the explicit expression for $\delta\alpha$ from Eq. (\[Jones\_Result\]) we reproduce the result of [@Jack:1996vg] for ${\cal N}=1$ SQED:

$$\delta\widetilde\beta(\alpha) = \widetilde
\beta_{\overline{\mbox{\scriptsize{DR}}}} - \widetilde
\beta_{\mbox{\scriptsize NSVZ}} = -\frac{\alpha^2}{\pi}\cdot
\frac{3\alpha^2}{4\pi^2} + \frac{\alpha^3}{4\pi^2} \cdot
\frac{2\alpha}{\pi} +O(\alpha^5) = -\frac{\alpha^4}{4\pi^3} +
O(\alpha^5).$$

In this paper we relate the results obtained with the higher derivative regularization with the results obtained in [@Jack:1996vg]. More exactly, we construct a finite renormalization relating the NSVZ scheme with the $\overline{\mbox{DR}}$ scheme and boundary conditions for which the functions $\widetilde\beta(\alpha)$ and $\widetilde\gamma(\alpha)$, depending on the renormalized coupling constant, coincide with the results obtained in the $\overline{\mbox{DR}}$ scheme.

The coefficient $b_1$ is found by comparing the two-loop expression for the anomalous dimension of the matter superfield

$$\label{Gamma_DRED}
\widetilde\gamma_{\overline{\mbox{\scriptsize{DR}}}}(\alpha) =
-\frac{\alpha}{\pi} + \frac{\alpha^2}{\pi^2} + O(\alpha^3)$$

with Eq. (\[Gamma\_Answer2\]). Two these expressions coincide, if

$$b_1 - g_1 = \sum\limits_{I=1}^n c_I \ln a_I + \frac{1}{2}.$$

The next coefficient $b_2$ can be found by comparing the results for the three-loop $\beta$-function. The result of [@Jack:1996vg] in our notation can be written as

$$\label{Beta_DRED}
\widetilde\beta_{\overline{\mbox{\scriptsize{DR}}}}(\alpha) =
\frac{\alpha^2}{\pi} + \frac{\alpha^3}{\pi^2} -
\frac{5\alpha^4}{4\pi^3} + O(\alpha^5).$$

Comparing this equation with (\[Beta\_Answer2\]) we obtain

$$b_2 - g_1 = \frac{1}{4}.$$

Therefore, for $x_0 = \ln\Lambda/\mu_0 =- g_1$

$$\begin{aligned}
&& Z_{\overline{\mbox{\scriptsize{DR}}}}(\alpha_0,x_0)=1 + O(\alpha_0^2); \nonumber\\
&&
\frac{1}{\alpha_{\overline{\mbox{\scriptsize{DR}}}}(\alpha_0,x_0)}
= \frac{1}{\alpha_0} + \frac{1}{\pi}\Big(\sum\limits_{I=1}^n c_I
\ln a_I + \frac{1}{2}\Big) + \frac{\alpha_0}{4\pi^2} +
O(\alpha_0^2).\end{aligned}$$

In the considered approximation these equations define a scheme in which the RG functions (\[Beta\_Definition2\]) and (\[Gamma\_Definition2\]) coincide with the corresponding RG functions obtained with the $\overline{\mbox{DR}}$ scheme. (It is expedient to compare these conditions with Eq. (\[NSVZ\_Scheme\]) defining the NSVZ scheme.)

It is easy to see that $Z_{\overline{\mbox{\scriptsize{DR}}}}$ and $\alpha_{\overline{\mbox{\scriptsize{DR}}}}$ can be obtained from $Z_{\mbox{\scriptsize NSVZ}}$ and $\alpha_{\mbox{\scriptsize
NSVZ}}$ (constructed with the higher derivative regularization) by the following finite renormalization:

$$\begin{aligned}
\label{NSVZ_DRED}
&&\hspace*{-7mm} Z_{\overline{\mbox{\scriptsize{DR}}}}
(\alpha_{\overline{\mbox{\scriptsize{DR}}}},\Lambda/\mu)=\Big(1 -
\frac{\alpha_{\mbox{\scriptsize NSVZ}}}{\pi}
\Big(\sum\limits_{I=1}^n c_I \ln a_I + \frac{1}{2}-b_1\Big) +
O(\alpha_{\mbox{\scriptsize NSVZ}}^2)\Big) Z_{\mbox{\scriptsize
NSVZ}}(\alpha_{\mbox{\scriptsize
NSVZ}},\Lambda/\mu); \nonumber\\
&&\hspace*{-7mm}
\frac{1}{\alpha_{\overline{\mbox{\scriptsize{DR}}}}} =
\frac{1}{\alpha_{\mbox{\scriptsize NSVZ}}} + \frac{b_1}{\pi} +
\frac{\alpha_{\mbox{\scriptsize NSVZ}}}{4\pi^2} -
\frac{\alpha_{\mbox{\scriptsize
NSVZ}}}{\pi^2}\Big(\sum\limits_{I=1}^n c_I \ln a_I + \frac{1}{2} -
b_1\Big) + O(\alpha_{\mbox{\scriptsize NSVZ}}^2),\end{aligned}$$

where $b_1$ is an arbitrary finite constant. For

$$b_1= \sum\limits_{I=1}^n c_I \ln a_I + \frac{1}{2}$$

Eq. (\[NSVZ\_DRED\]) corresponds to the result of Ref. [@Jack:1996vg], which is given by Eq. (\[Jones\_Result\]), because in [@Jack:1996vg] the authors construct the NSVZ scheme in which the anomalous dimension is fixed by the condition

$$Z_{\overline{\mbox{\scriptsize{DR}}}}
(\alpha_{\overline{\mbox{\scriptsize{DR}}}},\Lambda/\mu)=
Z_{\mbox{\scriptsize NSVZ}}(\alpha_{\mbox{\scriptsize
NSVZ}},\Lambda/\mu).$$

In the NSVZ scheme obtained in this paper at the one-loop level this condition can be satisfied by a special choice of the constant $b_1$. Fixing this additional constant we construct a finite renormalization after which not only the $\beta$-functions, but also the anomalous dimensions coincide. In the next orders the situation is the same: it is possible to redefine the coupling constant so that the NSVZ relation is valid for the RG functions defined in terms of the renormalized coupling constant. The finite renormalization of the matter superfields is not fixed by this requirement. This renormalization can be found requiring coincidence of the anomalous dimensions.

Conclusion
==========

For ${\cal N}=1$ SQED regularized by higher derivatives the NSVZ relation is automatically satisfied, if a $\beta$-function and an anomalous dimension are defined in terms of the bare coupling constant according to Eqs. (\[Beta\_Definition1\]) and (\[Gamma\_Definition1\]). If these definitions are used, then the regularization described in Sect. \[Section\_Regularization\] always gives the exact NSVZ $\beta$-function. In order to obtain a result different from the exact NSVZ $\beta$-function in this case, it is necessary to make a finite rescaling of the Pauli–Villars masses and a finite tuning of the bare charge $\alpha_0$. In particular, using this procedure it is possible to make coefficients of the RG functions defined in terms of the bare charge equal to the ones of the RG functions defined in terms of the renormalized charge in the $\overline{\mbox{DR}}$ scheme.

If a $\beta$-function and an anomalous dimension are defined in terms of the renormalized coupling constant, the NSVZ $\beta$-function is obtained only in a special (NSVZ) scheme, which is related with the MOM scheme by a finite renormalization. In case of using the DRED regularization the only way to construct this finite renormalization is to use the definition of the NSVZ scheme. However, if the theory is regularized by higher derivatives there is a concrete extra prescription: there should be a point $x_0$ in which

$$\label{NSVZ_Scheme_In_Conclusion}
Z_3(\alpha,x_0) = 1;\qquad Z(\alpha,x_0)=1.$$

In order to obtain a result different from the NSVZ $\beta$-function, it is possible to make a finite renormalization, which, in turn, changes the boundary conditions.

These results are verified by an explicit three-loop calculation. In particular, we relate the $\overline{\mbox{DR}}$ scheme and the NSVZ scheme obtained with the higher derivative regularization by a finite renormalization of both the coupling constant and the matter superfields similar to the result of [@Jack:1996vg]. However, an additional finite constant which gives the required value of the two-loop anomalous dimension should be fixed.

[**Acknowledgements.**]{}

We are grateful to V.A.Rubakov, I.V.Tyutin, and V.N.Velizhanin for valuable discussions. The work of AK is supported in part by the Grant NSh-5590.2012.2, Russian Foundation of Basic Research grants No. 11-01-00182 and No. 11-02-00112 and is done within the framework of the Ministry of Science and Education contract No. 8412. The work of KS was supported by Russian Foundation for Basic Research grant No. 11-01-00296.

[**Appendixes**]{}

Three-loop renormalization
==========================

\[Appendix\_Renormalization\]

In this appendix we describe the renormalization of ${\cal N}=1$ SQED, regularized by higher derivatives, at the three-loop level following Refs. [@Soloshenko:2003sx; @Soloshenko:2003nc]. We start with the expression (\[G(alpha0)\]) for the two-point Green function of the matter superfields. We split the bare coupling constant $\alpha_0 = e_0^2/4\pi$ into the renormalized coupling constant $\alpha$ and a counterterm. In the lowest approximation this can be done according to Eq. (\[One-Loop\_Coupling\_Renormalization\]), in which the logarithm compensate a one-loop divergence, and a finite parameter $b_1$ is not fixed. The expression (\[G(alpha0)\]) was calculated in Ref. [@Soloshenko:2003sx] for $R_k = 1+k^{2n}/\Lambda^{2n}$. The result can be written in the following form:

$$\begin{aligned}
&& G(\alpha_0,\Lambda/p) = 1 - \frac{\alpha_0}{\pi}
\ln\frac{\Lambda}{p} -\frac{\alpha_0}{2\pi} +
\frac{\alpha_0^2}{\pi^2} \ln{}^2 \frac{\Lambda}{p} +
\frac{\alpha_0^2}{\pi^2}\ln\frac{\Lambda}{p}\Big(\sum\limits_{I=1}^n
c_I \ln a_I +\frac{5}{2}\Big) + \frac{\alpha_0^2}{\pi^2}
c_2\nonumber\\
&& + \Big(\mbox{terms vanishing in the limit $\Lambda\to
\infty$}\Big) + O(\alpha_0^3),\end{aligned}$$

where the constant $c_2$ was not found in Ref. [@Soloshenko:2003sx]. (A finite part of this function in the one-loop approximation has been explicitly calculated.) The finite constants $a_I$ are related with the Pauli–Villars masses by Eq. (\[Scheme\]). In terms of the renormalized coupling constant the function $G$ is given by

$$\begin{aligned}
\label{G(alpha)}
&& G(\alpha,\Lambda/\mu,\Lambda/p) = 1 - \frac{\alpha}{\pi}
\Big(\ln \frac{\Lambda}{p} + \frac{1}{2}\Big)
+\frac{\alpha^2}{\pi^2} \Big[\ln{}^2 \frac{\Lambda}{p}
-\ln\frac{\Lambda}{\mu} \ln\frac{\Lambda}{p} - b_1 \ln
\frac{\Lambda}{p} - \frac{1}{2} \ln \frac{\Lambda}{\mu} -
\frac{1}{2} b_1\qquad\nonumber\\
&& + \ln\frac{\Lambda}{p} \Big(\sum\limits_{I=1}^n c_I \ln a_I
+\frac{5}{2}\Big) + c_2\Big] + (\mbox{terms vanishing in the limit
$\Lambda\to \infty$}) + O(\alpha^3).\end{aligned}$$

This expression should be multiplied by a factor $Z$ such that the expression $ZG$ is finite. This is true, if the renomalization constant $Z$ is given by Eq. (\[Two\_Loop\_Z\]):

$$Z = 1 + \frac{\alpha}{\pi}\Big(\ln\frac{\Lambda}{\mu}+g_1\Big)
+\frac{\alpha^2}{\pi^2}\ln^2\frac{\Lambda}{\mu}
-\frac{\alpha^2}{\pi^2}\ln\frac{\Lambda}{\mu}\Big(\sum\limits_{I=1}^n
c_I\ln a_I - b_1 + \frac{3}{2} - g_1\Big) + \frac{\alpha^2}{\pi^2}
g_2 + O(\alpha^3),$$

where $g_1$ and $g_2$ are undefined finite constants. Really, calculating the product $Z G$ in the limit $\Lambda\to \infty$, it is easy to see that all terms containing $\ln \Lambda$ cancel each other and

$$\begin{aligned}
&& G_{\mbox{\scriptsize ren}}(\alpha,\mu/p) =
\lim\limits_{\Lambda\to \infty} ZG = 1 - \frac{\alpha}{\pi}
\Big(\ln \frac{\mu}{p} - g_1 +\frac{1}{2}\Big) +
\frac{\alpha^2}{\pi^2} \Big[\ln{}^2 \frac{\mu}{p} +
\ln\frac{\mu}{p} \Big(\sum\limits_{I=1}^n
c_I \ln a_I +\frac{5}{2}\Big)\nonumber\\
&& - (b_1+g_1) \ln\frac{\mu}{p} + c_2 - \frac{1}{2} g_1  + g_2 -
\frac{1}{2} b_1\Big] + O(\alpha^3).\qquad\end{aligned}$$

This expression is finite for all finite values of $b_1$, $g_1$, and $g_2$. Therefore, they can be chosen in an arbitrary way (unlike the constant $c_2$, which is fixed, but unknown).

Values of the constants $b_i$ and $g_i$ can be found from boundary conditions, which specify the scheme. For example, the MOM scheme is defined by the boundary conditions (\[MOM\_Boundary\_Conditions\]):

$$Z_{\mbox{\scriptsize MOM}}G(\alpha_{\mbox{\scriptsize MOM}},p=\mu)
= 1; \qquad d^{-1}(\alpha_{\mbox{\scriptsize
MOM}},p=\mu)=\alpha_{\mbox{\scriptsize MOM}}^{-1}.$$

In the considered case the first condition gives

$$g_1 = \frac{1}{2};\qquad g_2 = - c_2 + \frac{1}{2} g_1
+\frac{1}{2} b_1 = - c_2 + \frac{1}{4} +\frac{1}{2} b_1.$$

(The constant $b_1$ is defined from the other boundary condition.) Therefore, in the subtraction scheme defined by the boundary conditions (\[MOM\_Boundary\_Conditions\]) the constants $b_i$ and $g_i$ are related with the finite parts of the Green function ($c_1 = -1/2$ and $c_2$).

However, it is possible to impose the different boundary conditions (\[NSVZ\_Scheme\_Z3\]):

$$Z_{\mbox{\scriptsize NSVZ}}(\alpha_{\mbox{\scriptsize
NSVZ}},\mu=\Lambda)=1; \qquad (Z_3)_{\mbox{\scriptsize
NSVZ}}(\alpha_{\mbox{\scriptsize NSVZ}},\mu=\Lambda)=1.$$

In this case

$$1 = Z_{\mbox{\scriptsize NSVZ}}(\alpha_{\mbox{\scriptsize
NSVZ}},\Lambda=\mu) = 1 + \frac{\alpha_{\mbox{\scriptsize
NSVZ}}}{\pi}g_1 + \frac{\alpha^2_{\mbox{\scriptsize NSVZ}}}{\pi^2}
g_2 + O(\alpha_{\mbox{\scriptsize NSVZ}}^3),$$

so that

$$g_1=0;\qquad g_2=0.$$

In this case the constants $g_i$ (and $b_i$, see below) are independent of the finite parts of the effective action (i.e. the constants $-1/2$ and $c_2$).

Let us proceed to the renormalization of the coupling constant. For this purpose we should write the function $d^{-1}$. The integrals defining this function in the three-loop approximation have been found in Ref. [@Soloshenko:2003nc]. Although it is very difficult to calculate their finite parts, we know [@Stepanyantz:2011jy] that this function satisfies Eq. (\[Green\_Function\_Relation\]):

$$\frac{d}{d\ln\Lambda} \Big(d^{-1}(\alpha_0,\Lambda/p) -
\alpha_0^{-1}\Big)\Big|_{p=0}  =
\frac{1}{\pi}\Big(1-\frac{d}{d\ln\Lambda} \ln
G(\alpha_0,\Lambda/q)\Big|_{q=0}\Big).$$

This equation gives the NSVZ relation in terms of the bare charge, i.e. for the RG functions (\[Beta\_Definition1\]) and (\[Gamma\_Definition1\]). The first term in the left hand side vanishes, because the function $d^{-1}$ depends only on $\alpha$ and $\mu/p$ and (expressed in terms of the renormalized coupling constant) is independent of $\Lambda$ in the limit $p\to 0$. The right hand side of this equation can be calculated differentiating the logarithm of Eq. (\[G(alpha)\]), taking into account that the limit $\Lambda\to \infty$ in the massless case is equivalent to the limit $q\to 0$ due to the dependence on $q/\Lambda$. Then it is easy to see that all terms containing the dependence on the momentum $q$ vanish. Expressing the result in terms of the bare coupling constant we obtain Eq. (\[Beta\_Expression\]):

$$- \frac{d}{d\ln\Lambda}(\alpha_0^{-1}) = \frac{1}{\pi} +
\frac{\alpha_0}{\pi^2}  - \frac{\alpha_0^2}{\pi^3}
\Big(\sum\limits_{I=1}^n c_I \ln a_I +\frac{3}{2}\Big) +
O(\alpha_0^3).$$

A renormalization of the coupling constant can be found by integrating this equation, which should be written in terms of the renormalized coupling constant:

$$- \frac{d}{d\ln\Lambda} (\alpha_0^{-1}) = \frac{1}{\pi} +
\frac{\alpha}{\pi^2} +\frac{\alpha^2}{\pi^3} \ln
\frac{\Lambda}{\mu} +\frac{\alpha^2}{\pi^3} b_1 -
\frac{\alpha^2}{\pi^3} \Big(\sum\limits_{I=1}^n c_I \ln a_I
+\frac{3}{2}\Big) + O(\alpha^3).$$

The result is given by Eq. (\[Alpha0\]):

$$\begin{aligned}
&& \frac{1}{\alpha_0} = \frac{1}{\alpha} - \frac{1}{\pi} \Big(
\ln\frac{\Lambda}{\mu} + b_1\Big) - \frac{\alpha}{\pi^2}
\Big(\ln\frac{\Lambda}{\mu} + b_2\Big)\nonumber\\
&& \qquad\qquad\qquad -
\frac{\alpha^2}{\pi^3}\Big(\frac{1}{2}\ln^2 \frac{\Lambda}{\mu}
-\ln \frac{\Lambda}{\mu} \sum\limits_{I=1}^n c_I \ln a_I
-\frac{3}{2} \ln \frac{\Lambda}{\mu} + b_1 \ln\frac{\Lambda}{\mu}
+ b_3\Big) + O(\alpha^3),\qquad\nonumber\end{aligned}$$

where $b_1$, $b_2$, and $b_3$ are integration constants in each order.

Although we have not calculated the inverse invariant charge $d^{-1}$, starting from this equation it is possible to restore its divergent part by requiring its finiteness after the substitution $\alpha_0 \to \alpha_0(\alpha,\Lambda/\mu)$:

$$\begin{aligned}
&& d^{-1}(\alpha_0,\Lambda/p) = \frac{1}{\alpha_0} +
\frac{1}{\pi}\Big(\ln\frac{\Lambda}{p} + d_1\Big) +
\frac{\alpha_0}{\pi^2}\Big(\ln\frac{\Lambda}{p} +
d_2\Big)\nonumber\\
&& + \frac{\alpha_0^2}{\pi^3}\Big(-\frac{1}{2}\ln^2
\frac{\Lambda}{p} -\ln \frac{\Lambda}{p} \sum\limits_{I=1}^n c_I
\ln a_I -\frac{3}{2} \ln \frac{\Lambda}{p} - d_2
\ln\frac{\Lambda}{p} + d_3\Big)\nonumber\\
&& + \mbox{(terms vanishing in the limit $\Lambda\to \infty$)} +
O(\alpha_0^3),\vphantom{\Big(}\end{aligned}$$

where $d_1$, $d_2$, and $d_3$ are finite constants, which should be found by calculating Feynman graphs. It is easy to see that for any values of $d_i$ and for an arbitrary choice of $b_i$ the above expression for $d^{-1}$ is a finite function of $\alpha$ and $\mu/p$:

$$\begin{aligned}
&& d^{-1}\Big(\alpha_0(\alpha,\Lambda/\mu),\Lambda/p\Big) =
\frac{1}{\alpha} + \frac{1}{\pi}\Big(\ln\frac{\mu}{p} + d_1 -
b_1\Big) + \frac{\alpha}{\pi^2}\Big(\ln\frac{\mu}{p} + d_2 -
b_2\Big)\nonumber\\
&& + \frac{\alpha^2}{\pi^3}\Big(-\frac{1}{2}\ln^2 \frac{\mu}{p}
-\ln \frac{\mu}{p} \sum\limits_{I=1}^n c_I \ln a_I -\frac{3}{2}
\ln \frac{\mu}{p} + (b_1 - d_2) \ln\frac{\mu}{p} + d_3 - b_3 + b_1
d_2\Big)\nonumber\\
&& + \mbox{(terms vanishing in the limit $\Lambda\to \infty$)} +
O(\alpha^3).\vphantom{\Big(}\end{aligned}$$

The coefficients $b_i$ can be found from boundary conditions. If the boundary conditions (\[MOM\_Boundary\_Conditions\]), defining the MOM scheme, are imposed, we obtain:

$$b_1 = d_1;\qquad b_2 = d_2;\qquad b_3 = d_3 + b_1 d_2 = d_3+ d_1
d_2.$$

Therefore, in this case the coefficients defining the scheme are expressed in terms of the finite parts of the Green function. In the limit $\Lambda\to \infty$ it is easy to see that in the MOM scheme

$$\begin{aligned}
&& Z_{\mbox{\scriptsize MOM}}G(\alpha_{\mbox{\scriptsize
MOM}},\mu/p)
= 1 - \frac{\alpha_{\mbox{\scriptsize MOM}}}{\pi} \ln \frac{\mu}{p}\\
&&\qquad\qquad\qquad\qquad\qquad + \frac{\alpha_{\mbox{\scriptsize
MOM}}^2}{\pi^2} \Big(\ln{}^2 \frac{\mu}{p} + \ln \frac{\mu}{p}
\Big( \sum\limits_{I=1}^n c_I \ln a_I + 2 -d_1\Big) \Big)
+ O(\alpha_{\mbox{\scriptsize MOM}}^3);\nonumber\\
&& d^{-1}(\alpha_{\mbox{\scriptsize MOM}},\mu/p) =
\frac{1}{\alpha_{\mbox{\scriptsize MOM}}} + \frac{1}{\pi}
\ln\frac{\mu}{p} + \frac{\alpha_{\mbox{\scriptsize MOM}}}{\pi^2} \ln\frac{\mu}{p}\\
&&\qquad\ + \frac{\alpha_{\mbox{\scriptsize
MOM}}^2}{\pi^3}\Big(-\frac{1}{2}\ln^2 \frac{\mu}{p} -\ln
\frac{\mu}{p} \sum\limits_{I=1}^n c_I \ln a_I -\frac{3}{2} \ln
\frac{\mu}{p} + (d_1 - d_2) \ln\frac{\mu}{p}\Big) +
O(\alpha_{\mbox{\scriptsize MOM}}^3).\qquad\nonumber\end{aligned}$$

If we impose the boundary conditions (\[NSVZ\_Scheme\_Z3\]), defining the NSVZ scheme, then the coefficients $b_i$ do not depend on the finite parts of the Green function:

$$\frac{1}{\alpha_{\mbox{\scriptsize NSVZ}}} =  \frac{1}{\alpha_0} =
\frac{1}{\alpha_{\mbox{\scriptsize NSVZ}}} - \frac{1}{\pi} b_1 -
\frac{\alpha_{\mbox{\scriptsize NSVZ}}}{\pi^2} b_2 -
\frac{\alpha_{\mbox{\scriptsize NSVZ}}^2}{\pi^3} b_3 +
O(\alpha_{\mbox{\scriptsize NSVZ}}^3),$$

so that

$$b_1=0; \qquad b_2=0; \qquad b_3=0.$$

In this case the renormalized Green functions in the limit $\Lambda\to \infty$ are written as

$$\begin{aligned}
&& Z_{\mbox{\scriptsize NSVZ}}G(\alpha_{\mbox{\scriptsize
NSVZ}},\mu/p) = 1 - \frac{\alpha_{\mbox{\scriptsize NSVZ}}}{\pi}
\Big(\ln \frac{\mu}{p}+\frac{1}{2}\Big)\\
&&\qquad\qquad\qquad\qquad\quad + \frac{\alpha_{\mbox{\scriptsize
NSVZ}}^2}{\pi^2}\Big[\ln{}^2 \frac{\mu}{p} + \ln\frac{\mu}{p}
\Big(\sum\limits_{I=1}^n c_I \ln a_I + \frac{5}{2}\Big) +
c_2\Big] + O(\alpha_{\mbox{\scriptsize NSVZ}}^3);\nonumber\\
&& d^{-1}(\alpha_{\mbox{\scriptsize NSVZ}},\mu/p) =
\frac{1}{\alpha_{\mbox{\scriptsize NSVZ}}} +
\frac{1}{\pi}\Big(\ln\frac{\mu}{p} + d_1\Big) +
\frac{\alpha_{\mbox{\scriptsize NSVZ}}}{\pi^2}\Big(\ln\frac{\mu}{p} + d_2\Big)\\
&&\qquad + \frac{\alpha_{\mbox{\scriptsize
NSVZ}}^2}{\pi^3}\Big(-\frac{1}{2}\ln^2 \frac{\mu}{p} -\ln
\frac{\mu}{p} \sum\limits_{I=1}^n c_I \ln a_I -\frac{3}{2} \ln
\frac{\mu}{p} - d_2 \ln\frac{\mu}{p} + d_3\Big) +
O(\alpha_{\mbox{\scriptsize NSVZ}}^3).\qquad\nonumber\end{aligned}$$

It is easy to verify that the MOM scheme defined by the boundary conditions (\[MOM\_Boundary\_Conditions\]) and the NSVZ scheme defined by the boundary conditions (\[NSVZ\_Scheme\_Z3\]) are related by the following finite renormalization:

$$\begin{aligned}
&& \frac{1}{\alpha_{\mbox{\scriptsize NSVZ}}} =
\frac{z_3(\alpha_{\mbox{\scriptsize
MOM}})}{\alpha_{\mbox{\scriptsize MOM}}} =
\frac{1}{\alpha_{\mbox{\scriptsize MOM}}}- \frac{1}{\pi} d_1 -
\frac{\alpha_{\mbox{\scriptsize MOM}}}{\pi^2} d_2 -
\frac{\alpha_{\mbox{\scriptsize MOM}}^2}{\pi^3} (d_3+ d_1
d_2) + O(\alpha_{\mbox{\scriptsize MOM}}^3);\nonumber\\
&& Z_{\mbox{\scriptsize NSVZ}}(\alpha_{\mbox{\scriptsize
NSVZ}},\Lambda/\mu) = z(\alpha_{\mbox{\scriptsize MOM}})
Z_{\mbox{\scriptsize MOM}}(\alpha_{\mbox{\scriptsize
MOM}},\Lambda/\mu)
\vphantom{\frac{1}{2}}\\
&&\qquad\qquad\qquad = \Big(1-\frac{\alpha_{\mbox{\scriptsize
MOM}}}{2\pi} - \frac{\alpha_{\mbox{\scriptsize
MOM}}^2}{\pi^2}\Big(\frac{1}{2} d_1-c_2\Big) +
O(\alpha_{\mbox{\scriptsize MOM}}^3)\Big) Z_{\mbox{\scriptsize
MOM}}(\alpha_{\mbox{\scriptsize MOM}},\Lambda/\mu),
\nonumber\qquad\end{aligned}$$

where the remaining coefficients can be found by calculating finite parts of the Feynman diagrams.

Existence of the NSVZ scheme
============================

\[Appendix\_Existence\]

In this section we prove that the scheme defined by Eq. (\[NSVZ\_Scheme\]) exists and without loss of generality it is possible to set $x_0=0$. For this purpose we note that the renormalized coupling constant $\alpha$ is defined so that the function $d^{-1}\Big(\alpha_0(\alpha,\Lambda/\mu),\Lambda/p\Big)$ does not depend on $\Lambda$ in the limit $\Lambda\to \infty$. Certainly, this condition does not uniquely define $\alpha$. In particular, it is possible to perform a finite redefinition of the coupling constant $\alpha = \alpha(\alpha')$. After this redefinition

$$\alpha_0(\alpha,\Lambda/\mu) \to
\alpha_0(\alpha(\alpha'),\Lambda/\mu).$$

Because the function $\alpha(\alpha')$ is finite, the expression

$$d^{-1}\Big(\alpha_0(\alpha(\alpha'),\Lambda/\mu),\Lambda/p\Big)$$

is also finite. Therefore, $\alpha'$ can be also chosen as a renormalized coupling constant.

Let us choose the renormalized coupling constant in an arbitrary way. It is defined by a function $\alpha_0(\alpha,\Lambda/\mu)$. Then we construct a new function

$$\label{B_Definition}
\alpha_0(\alpha,x=0) \equiv b(\alpha),$$

where

$$x \equiv \ln\frac{\Lambda}{\mu}.$$

Evidently, the function $b(\alpha)$ does not depend on $\Lambda$ and is finite in the limit $\Lambda\to \infty$. Then we define a new renormalized coupling constant according to the prescription

$$\alpha'(\alpha) \equiv b(\alpha).$$

Then, from Eq. (\[B\_Definition\]) we obtain

$$\alpha_0(\alpha(\alpha'),x=0) = \alpha'.$$

Therefore, choosing $\alpha'$ as a renormalized coupling constant we satisfy the first condition in Eq. (\[NSVZ\_Scheme\]) with $x_0=0$.

Similarly, by definition, the renormalization constant $Z$ is chosen in such a way that the expression

$$Z(\alpha,\Lambda/\mu)
G\Big(\alpha_0(\alpha,\Lambda/\mu),\Lambda/p\Big)$$

is finite. Again, this condition does not uniquely define the function $Z$. In particular, it is possible to multiply $Z$ by an arbitrary finite function of the renormalized coupling constant $\alpha$. This possibility can be used for constructing the renormalization constant which satisfies the second condition in Eq. (\[NSVZ\_Scheme\]). For this purpose we consider an arbitrary renormalization constant $Z$ and construct the function

$$\label{G_Definition}
Z(\alpha,x=0) \equiv g(\alpha).$$

Then the renormalization constant $Z'$ is defined according to the prescription

$$\label{Z_Prime_Definition}
Z'(\alpha',\Lambda/\mu) \equiv \frac{1}{g(\alpha(\alpha'))}
Z\Big(\alpha(\alpha'),\Lambda/\mu\Big).$$

It is evident that the expression

$$Z'(\alpha',\Lambda/\mu)
G\Big(\alpha_0(\alpha',\Lambda/\mu),\Lambda/p\Big)$$

does not depend on $\Lambda$ in the limit $\Lambda\to
\infty$. Therefore the function $Z'(\alpha',\Lambda/\mu)$ can be chosen as a renormalization constant. Moreover, due to Eqs. (\[G\_Definition\]) and (\[Z\_Prime\_Definition\]) this function satisfies the condition

$$Z'(\alpha',x=0) = 1.$$

Therefore, we have constructed the renormalization constants $Z_3'(\alpha',\Lambda/\mu)$ and $Z'(\alpha',\Lambda/\mu)$ which satisfy the conditions (\[NSVZ\_Scheme\_Z3\]) with $x_0=0$.

Let us verify these results by an explicit three-loop calculation. &gt;From Eq. (\[Alpha0\]) we obtain

$$\frac{1}{\alpha_0(\alpha,x=0)} = \frac{1}{\alpha} - \frac{1}{\pi}
b_1 - \frac{\alpha}{\pi^2} b_2 - \frac{\alpha^2}{\pi^3} b_3 +
O(\alpha^3).$$

As a consequence, the function $b(\alpha)$ defined by Eq. (\[B\_Definition\]) is given by

$$\label{Three_Loop_Alpha_Prime}
b(\alpha) = \alpha \Big(1 - \frac{\alpha}{\pi} b_1 -
\frac{\alpha^2}{\pi^2} b_2 - \frac{\alpha^3}{\pi^3} b_3\Big)^{-1}
+ O(\alpha^5)\equiv \alpha'(\alpha).$$

Rewriting the bare coupling constant $\alpha_0$ in terms of $\alpha'$ defined by this equation we obtain

$$\frac{1}{\alpha_0} = \frac{1}{\alpha'} - \frac{1}{\pi}
\ln\frac{\Lambda}{\mu} - \frac{\alpha'}{\pi^2}
\ln\frac{\Lambda}{\mu} -\frac{(\alpha')^2}{\pi^3}
\Big(\frac{1}{2}\ln{}^2 \frac{\Lambda}{\mu} -
\ln\frac{\Lambda}{\mu} \sum\limits_{I=1}^n c_I \ln
\frac{M_I}{\Lambda} - \frac{3}{2} \ln \frac{\Lambda}{\mu}\Big) +
O\left((\alpha')^3\right).$$

This choice is a particular case of Eq. (\[Alpha0\]) with $b_1=b_2=b_3=0$. Therefore, the function $d^{-1}\left(\alpha_0(\alpha',\Lambda/\mu),\Lambda/p\right)$ does not depend on $\Lambda$ in the limit $\Lambda\to \infty$. Moreover, it is evident that

$$\label{Three_Loop_NSVZ_Scheme_Alpha}
\alpha_0(\alpha',x=0) = \alpha' + O\left((\alpha')^5\right).$$

Similarly, one can construct the function $Z'(\alpha',\Lambda/\mu)$ in the two-loop approximation. As a starting point we consider the renormalization constant (\[Two\_Loop\_Z\]). Then we find

$$g(\alpha)= Z(\alpha,x=0) = 1 + \frac{\alpha}{\pi} g_1 +
\frac{\alpha^2}{\pi^2} g_2 + O(\alpha^3).$$

Writing this function in terms of $\alpha'$ defined by Eq. (\[Three\_Loop\_Alpha\_Prime\]) we obtain

$$g\left(\alpha(\alpha')\right) = 1 + \frac{\alpha'}{\pi} g_1 +
\frac{(\alpha')^2}{\pi^2} (g_2 - g_1 b_1) +
O\left((\alpha')^3\right).$$

As a consequence,

$$\begin{aligned}
&& Z'(\alpha',\Lambda/\mu) = g\left(\alpha(\alpha')\right)^{-1}
Z\Big(\alpha(\alpha'),\Lambda/\mu\Big)\nonumber\\
&& = 1 + \frac{\alpha'}{\pi}\ln\frac{\Lambda}{\mu} +
\frac{(\alpha')^2}{\pi^2} \ln{}^2\frac{\Lambda}{\mu} -
\frac{(\alpha')^2}{\pi^2}\Big( \sum\limits_{I=1}^n c_I \ln
\frac{M_I}{\Lambda} + \frac{3}{2}\Big) +
O\left((\alpha')^3\right).\end{aligned}$$

This function can be obtained from Eq. (\[Two\_Loop\_Z\]) if one sets $g_1=g_2 = 0$. (Earlier we have already set $b_1=0$.) Therefore, the finiteness of the renormalized two-point Green function of the matter superfields is guaranteed. Moreover, the function $Z'$ evidently satisfies the required condition

$$\label{Two_Loop_NSVZ_Scheme_Z}
Z'(\alpha',x=0) = 1 + O\left((\alpha')^3\right).$$

Thus, the NSVZ scheme is constructed in the three-loop approximation. Certainly, this procedure can be generalized to higher orders. For example, considering the next order of the perturbation theory we find a scheme in which Eq. (\[Three\_Loop\_NSVZ\_Scheme\_Alpha\]) is valid also for terms proportional to $\alpha_0^5$, and Eq. (\[Two\_Loop\_NSVZ\_Scheme\_Z\]) is valid for terms proportional to $\alpha_0^3$.

[100]{}

V. I. Ogievetsky and L. Mezincescu, Sov. Phys. Usp.  [**18**]{} (1975) 960 \[Usp. Fiz. Nauk [**117**]{} (1975) 637\]. P. C. West, “Introduction to supersymmetry and supergravity,” Singapore, Singapore: World Scientific (1990) 425 p.

I. L. Buchbinder and S. M. Kuzenko, “Ideas and methods of supersymmetry and supergravity: Or a walk through superspace,” Bristol, UK: IOP (1998) 656 p.

“Subtleties in quantum field theory: Lev Lipatov Festschrift,” Gatchina, Russia: St. Petersburg Nucl. Phys. Inst. (2010) 267 pp, ed. D.I.Diakonov.

M. T. Grisaru and W. Siegel, Nucl. Phys. B [**201**]{} (1982) 292 \[Erratum-ibid. B [**206**]{} (1982) 496\]. S. Mandelstam, Nucl. Phys. B [**213**]{} (1983) 149. L. Brink, O. Lindgren and B. E. W. Nilsson, Nucl. Phys. B [**212**]{} (1983) 401.

P. S. Howe, K. S. Stelle and P. K. Townsend, Nucl. Phys. B [**236**]{} (1984) 125. L. V. Avdeev, O. V. Tarasov and A. A. Vladimirov, Phys. Lett. B [**96**]{} (1980) 94. W. Siegel, Phys. Lett. B [**84**]{} (1979) 193. M. T. Grisaru, M. Rocek and W. Siegel, Phys. Rev. Lett.  [**45**]{} (1980) 1063. W. E. Caswell and D. Zanon, Nucl. Phys. B [**182**]{} (1981) 125. G. ’t Hooft and M. J. G. Veltman, Nucl. Phys. B [**44**]{} (1972) 189. C. G. Bollini and J. J. Giambiagi, Nuovo Cim. B [**12**]{} (1972) 20. J. F. Ashmore, Lett. Nuovo Cim.  [**4**]{} (1972) 289. G. M. Cicuta and E. Montaldi, Lett. Nuovo Cim.  [**4**]{} (1972) 329. R. Delbourgo and V. B. Prasad, J. Phys. G [**1**]{} (1975) 377.

W. A. Bardeen, A. J. Buras, D. W. Duke and T. Muta, Phys. Rev. D [**18**]{} (1978) 3998. A. L. Kataev, M. D. Vardiashvili, Phys. Lett. B [**221**]{} (1989) 377 \[Erratum-ibid. B [**241**]{} (1990) 644\]. W. Siegel, Phys. Lett. B [**94**]{} (1980) 37. L. V. Avdeev, G. A. Chochia and A. A. Vladimirov, Phys. Lett. B [**105**]{} (1981) 272. L. V. Avdeev, Phys. Lett. B [**117**]{} (1982) 317. L. V. Avdeev and A. A. Vladimirov, Nucl. Phys. B [**219**]{} (1983) 262. V. N. Velizhanin, Phys. Lett. B [**696**]{} (2011) 560. P. S. Howe, K. S. Stelle and P. C. West, Phys. Lett. B [**124**]{} (1983) 55. L. V. Avdeev and O. V. Tarasov, Phys. Lett. B [**112**]{} (1982) 356. I. Jack, D. R. T. Jones, P. Kant and L. Mihaila, JHEP [**0709**]{} (2007) 058. V. N. Velizhanin, Nucl. Phys. B [**818**]{} (2009) 95. M. T. Grisaru, W. Siegel and M. Rocek, Nucl. Phys. B [**159**]{} (1979) 429. V. A. Novikov, M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B [**229**]{} (1983) 381. V. A. Novikov, M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Phys. Lett. B [**166**]{} (1986) 329; Sov. J. Nucl. Phys.  [**43**]{} (1986) 294; \[Yad. Fiz.  [**43**]{} (1986) 459.\] M. A. Shifman and A. I. Vainshtein, Nucl. Phys. B [**277**]{} (1986) 456; Sov. Phys. JETP [**64**]{} (1986) 428; \[Zh. Eksp. Teor. Fiz.  [**91**]{} (1986) 723.\] A. I. Vainshtein, V. I. Zakharov and M. A. Shifman, JETP Lett.  [**42**]{} (1985) 224 \[Pisma Zh. Eksp. Teor. Fiz.  [**42**]{} (1985) 182\]. M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Phys. Lett. B [**166**]{} (1986) 334. I. Jack, D. R. T. Jones and C. G. North, Phys. Lett. B [**386**]{} (1996) 138 \[hep-ph/9606323\].

I. Jack, D. R. T. Jones and C. G. North, Nucl. Phys. B [**486**]{} (1997) 479. I. Jack, D. R. T. Jones and A. Pickering, Phys. Lett. B [**435**]{} (1998) 61.

S. Ferrara and B. Zumino, Nucl. Phys. B [**79**]{} (1974) 413.

D. R. T. Jones, Nucl. Phys. B [**87**]{} (1975) 127.

R. V. Harlander, D. R. T. Jones, P. Kant, L. Mihaila and M. Steinhauser, JHEP [**0612**]{} (2006) 024 \[hep-ph/0610206\].

M. A. Shifman and A. I. Vainshtein, Sov. J. Nucl. Phys.  [**44**]{} (1986) 321 \[Yad. Fiz.  [**44**]{} (1986) 498\]. J. Mas, M. Perez-Victoria and C. Seijas, JHEP [**0203**]{} (2002) 049 \[hep-th/0202082\]. A. A. Slavnov, Nucl. Phys. B [**31**]{} (1971) 301. A. A. Slavnov, Theor.Math.Phys. [**13**]{} (1972) 1064 \[Teor. Mat. Fiz.  [**13**]{} (1972) 174\]. V. K. Krivoshchekov, Theor. Math. Phys. [**36**]{} (1978) 745 \[Teor. Mat. Fiz.  [**36**]{} (1978) 291\]. P. C. West, Nucl. Phys. B [**268**]{} (1986) 113. A. A. Soloshenko and K. V. Stepanyantz, Theor. Math. Phys.  [**140**]{} (2004) 1264 \[Teor. Mat. Fiz.  [**140**]{} (2004) 430\]. A. B. Pimenov, E. S. Shevtsova and K. V. Stepanyantz, Phys. Lett. B [**686**]{} (2010) 293. K. V. Stepanyantz, Phys. Part. Nucl. Lett.  [**8**]{} (2011) 321. K. V. Stepanyantz, Int. J. Theor. Phys.  [**51**]{} (2012) 276. A. V. Smilga and A. Vainshtein, Nucl. Phys. B [**704**]{} (2005) 445.

K. V. Stepanyantz, arXiv:1108.1491 \[hep-th\]. K. V. Stepanyantz, J. Phys. Conf. Ser.  [**343**]{} (2012) 012115. K. V. Stepanyantz, J. Phys. Conf. Ser.  [**368**]{} (2012) 012052 \[arXiv:1203.5525 \[hep-th\]\]. K. V. Stepanyantz, Nucl. Phys. B [**852**]{} (2011) 71. L. D. Faddeev and A. A. Slavnov, “Gauge Fields. Introduction To Quantum Theory,” Nauka, Moscow, 1978 and Front. Phys.  [**50**]{} (1980) 1 \[Front. Phys.  [**83**]{} (1990) 1\]. A. A. Slavnov, Theor. Math. Phys. [**33**]{} (1977) 977 \[Teor. Mat. Fiz.  [**33**]{} (1977) 210\]. N. N. Bogolyubov and D. V. Shirkov, “Introduction To The Theory Of Quantized Fields,” Nauka, Moscow, 1984; and Intersci. Monogr. Phys. Astron.  [**3**]{} (1959) 1.

J. C. Collins, “Renormalization. An Introduction To Renormalization, The Renormalization Group, And The Operator Product Expansion,” Cambridge, Uk: Univ. Pr. ( 1984) 380p

A. A. Vladimirov, Theor. Math. Phys.  [**25**]{} (1976) 1170 \[Teor. Mat. Fiz.  [**25**]{} (1975) 335\]. N. Arkani-Hamed and H. Murayama, JHEP [**0006**]{} (2000) 030 \[hep-th/9707133\]. S. G. Gorishny, A. L. Kataev, S. A. Larin and L. R. Surguladze, Phys. Lett. B [**256**]{} (1991) 81. A. L. Kataev, arXiv:1305.4605 \[hep-th\]. A. L. Kataev, Phys. Lett. B [**691**]{} (2010) 82. A. A. Soloshenko and K. V. Stepanyants, Theor. Math. Phys.  [**134**]{} (2003) 377 \[Teor. Mat. Fiz.  [**134**]{} (2003) 430\].

[^1]: The detailed definition of the $\overline{\mbox{MS}}$ scheme is given in Ref.[@Kataev:1988sq].

[^2]: The limit $p\to 0$, in which the derivative of the function $d^{-1}$ vanishes, is obtained due to the factor $\delta^4(p)$.

[^3]: The expression presented in [@Soloshenko:2003sx] is slightly different, because it includes diagrams with one-loop counterterm insertions on the matter lines and, therefore, corresponds to $Z_{\mbox{\scriptsize 1-loop}} G_{\mbox{\scriptsize 2-loop}}$.

[^4]: In order to obtain the results of Ref. [@Jack:1996vg], it is necessary to set $\alpha =
    g^2/4\pi$, $\gamma(\alpha) = 2\gamma(g)$, $\beta(\alpha) = g
    \beta(g)/2\pi$.
