Landau-Khalatnikov-Fradkin transformation and the mystery of even $\zeta$-values in Euclidean massless correlators

TL;DR

This paper uses the Landau-Khalatnikov-Fradkin transformation to explain why even zeta-values appear in multi-loop Euclidean massless correlators, revealing a new basis that removes these values from perturbative series.

Contribution

It introduces a non-perturbative identity linking massless propagators in different gauges and constructs a basis that eliminates even zeta-values in perturbative expansions.

Findings

Elimination of even zeta-values in perturbative series.

Exact relation between hatted and standard zeta-functions.

Insight into the gauge dependence of massless correlators.

Abstract

The Landau-Khalatnikov-Fradkin (LKF) transformation is a powerful and elegant transformation allowing to study the gauge dependence of the propagator of charged particles interacting with gauge fields. With the help of this transformation, we derive a non-perturbative identity between massless propagators in two different gauges. From this identity, we find that the corresponding perturbative series can be exactly expressed in terms of a hatted transcendental basis that eliminates all even Euler $\zeta$-functions. This explains the mystery of even $\zeta$-values observed in multi-loop calculations of Euclidean massless correlators for almost three decades now. Our construction further allows us to derive an exact formula relating hatted and standard $\zeta$-functions to all orders of perturbation theory.

Equations

$${\hat{\zeta }}_3\equiv {\zeta }_3+\frac{3\epsilon }{2}{\zeta }_4-\frac{5{\epsilon }^3}{2}{\zeta }_6,{\hat{\zeta }}_5\equiv {\zeta }_5+\frac{5\epsilon }{2}{\zeta }_6,{\hat{\zeta }}_7\equiv {\zeta }_7,$$

$$S_F(p,\xi )=\frac{1}{i\hat{p}}P(p,\xi ),$$

$$S_F(x,\xi )=\hat{x}X(x,\xi ).$$

$$S_F(p,\xi )=\int \frac{{\mathrm{d}}^{d}x}{(2\pi {)}^{d/2}}{e}^{\mathrm{i}px}{S}_F(x,\xi ),$$

$$S_F(x,\xi )=\int \frac{{\mathrm{d}}^{d}p}{(2\pi {)}^{d/2}}{e}^{-\mathrm{i}px}{S}_F(p,\xi ).$$

$$S_F(x,\xi )=S_F(x,\eta )e^{\mathrm{i}(D(x)-D(0))},$$

$$D(x)=-\mathrm{i}\Delta e^2{\mu }^{4-d}\int \frac{{\mathrm{d}}^{d}p}{(2\pi {)}^d}\frac{e^{-\mathrm{i}px}}{p^4},\Delta =\xi -\eta .$$

$$T_{\alpha }(m^2)=\int \frac{{\mathrm{d}}^{d}p}{(2\pi {)}^d}\frac{e^{-\mathrm{i}px}}{(p^2+m^2{)}^{\alpha }}.$$

$$S_F(x,\xi )=S_F(x,\eta )e^{\mathrm{i}D(x)}.$$

$$\int {\mathrm{d}}^{d}x\frac{e^{\mathrm{i}px}}{x^{2\alpha }}=\frac{2^{2\tilde{\alpha }}{\pi }^{d/2}a(\alpha )}{p^{2\tilde{\alpha }}},a(\alpha )=\frac{\Gamma (\tilde{\alpha })}{\Gamma (\alpha )},\tilde{\alpha }=\frac{d}{2}-\alpha ,$$

$$\int {\mathrm{d}}^{d}p\frac{e^{-\mathrm{i}px}}{p^{2\alpha }}=\frac{2^{2\tilde{\alpha }}{\pi }^{d/2}a(\alpha )}{x^{2\tilde{\alpha }}}.$$

$$D(x)=-\mathrm{i}\Delta e^2({\mu }^2{x}^2{)}^{2-d/2}\frac{\Gamma (d/2-2)}{2^4(\pi {)}^{d/2}},$$

$$D(x)=\frac{\mathrm{i}\Delta A}{\epsilon }\Gamma (1-\epsilon )(\pi {\mu }^2{x}^2{)}^{\epsilon },A=\frac{{\alpha }_{\mathrm{em}}}{4\pi }=\frac{e^2}{(4\pi {)}^2}.$$

$$S_C(x,\xi )=S_C(x,\eta )e^{\mathrm{i}D(x)}.$$

$$P(p,\eta )=\sum _{m=0}^{\infty }{a}_m(\eta )A^m{\left(\frac{{\tilde{\mu }}^2}{p^2}\right)}^{m\epsilon }.$$

$${\tilde{\mu }}^2=4\pi {\mu }^2,$$

$$P(p,\xi )=\sum _{m=0}^{\infty }{a}_m(\xi )A^m{\left(\frac{{\tilde{\mu }}^2}{p^2}\right)}^{m\epsilon },$$

$$a_m(\xi )=a_m(\eta )\frac{\Gamma (2-(m+1)\epsilon )}{\Gamma (1+m\epsilon )}\times$$

$$\times \sum _{l=0}^{\infty }\frac{\Gamma (1+(m+l)\epsilon ){\Gamma }^l(1-\epsilon )}{l!\Gamma (2-(m+l+1)\epsilon )}\frac{(\Delta A{)}^l}{(-\epsilon {)}^l}{\left(\frac{{\tilde{\mu }}^2}{p^2}\right)}^{l\epsilon }.$$

$${\overline{\mu }}^2={\tilde{\mu }}^2{e}^{-\gamma },$$

$${\mu }_G^{2\epsilon }={\tilde{\mu }}^{2\epsilon }\frac{{\Gamma }^2(1-\epsilon )\Gamma (1+\epsilon )}{\Gamma (2-2\epsilon )}.$$

$${\mu }_g^{2\epsilon }={\tilde{\mu }}^{2\epsilon }\frac{{\Gamma }^2(1-\epsilon )\Gamma (1+\epsilon )}{\Gamma (1-2\epsilon )}.$$

$${\mu }_{\text{MV}}^{2\epsilon }=\frac{{\tilde{\mu }}^{2\epsilon }}{\Gamma (1-\epsilon )}.$$

$$a_m(\xi )$$

$$\times {\Phi }_p(m,l,\epsilon )\frac{(\Delta A{)}^l}{(-\epsilon {)}^{l}l!}{\left(\frac{{\mu }_p^2}{p^2}\right)}^{l\epsilon },$$

$${\Phi }_{\text{MV}}(m,l,\epsilon )=\frac{\Gamma (1-(m+1)\epsilon )\Gamma (1+(m+l)\epsilon ){\Gamma }^{2l}(1-\epsilon )}{\Gamma (1+m\epsilon )\Gamma (1-(m+l+1)\epsilon )},$$

$${\Phi }_g(m,l,\epsilon )={\Phi }_{\text{MV}}(m,l,\epsilon )\frac{{\Gamma }^l(1-2\epsilon )}{{\Gamma }^{3l}(1-\epsilon ){\Gamma }^l(1+\epsilon )},$$

$$\displaystyle\Gamma(1+\beta\varepsilon)=\exp\Big{[}-\gamma\beta\varepsilon+\sum_{s=2}^{\infty}\,(-1)^{s}\,\eta_{s}\beta^{s}\varepsilon^{s}\Bigr{]},~{}~{}\eta_{s}=\frac{\zeta_{s}}{s}\,.$$

$$\Phi_{\text{MV}}(m,l,\varepsilon)=\exp\Big{[}\sum_{s=2}^{\infty}\,\eta_{s}\,p_{s}(m,l)\,\varepsilon^{s}\Bigr{]}\,,$$

$$p_s(m,l)=$$

$$\displaystyle+(-1)^{s}\Bigl{\{}(m+l)^{s}-m^{s}\Bigr{\}}\,,$$

$$p_1(m,l)=0,p_2(m,l)=0.$$

$$S_C(p,\eta )=\frac{1}{p^2}\sum _{m=0}^{\infty }{a}_m^c(\eta )A^m{\left(\frac{{\tilde{\mu }}^2}{p^2}\right)}^{m\epsilon },$$

Full text

I Introduction

Gauge invariance governs the dynamics of systems of charged particles with deep consequences in elementary particle physics and beyond. Through the gauge principle, it gives rise to gauge field theories the prototype of which is quantum electrodynamics (QED). In the latter, gauge freedom can be expressed via a covariant gauge fixing procedure that introduces an explicit dependence of the Lagrangian on a gauge fixing parameter $\xi$. While physical quantities should not depend on this parameter, precious information can be obtained by studying the $\xi$-dependence of various correlation functions.

Such a task can be carried out with the help of the Landau-Khalatnikov-Fradkin (LKF) transformation Landau:1955zz (see also [Johnson:1959zz, ; Sonoda:2000kn, ]) that elegantly relates the QED fermion propagator in two different $\xi$-gauges (and similarly for the fermion-photon vertex). Its most important applications Curtis:1990zs ; Burden:1998gr ; Jia:2016udu are related to the study of the gauge covariance of QED Schwinger-Dyson equations and their solutions. This allows, e.g., to construct a charged-particle-photon vertex ansatz both in scalar Fernandez-Rangel:2016zac and spinor QED Kizilersu:2009kg . Other applications Bashir:2002sp ; Jia:2016wyu are focused on estimating large orders of perturbation theory. Indeed, and this will play a crucial role in what follows, the non-perturbative nature of the LKF transformation allows to fix some of the coefficients of the all-order expansion of the fermion propagator. Starting with a perturbative propagator in some fixed gauge, say $\eta$, all the coefficients depending on the difference between the gauge fixing parameters of the two propagators, $\xi-\eta$, get fixed by a weak coupling expansion of the LKF-transformed initial one. Such estimations have been carried out for QED in various dimensions (see [Bashir:2002sp, ; Jia:2016wyu, ]), for generalizations to brane worlds Ahmad:2016dsb and for more general SU(N) gauge theories DeMeerleer:2018txc .

A seemingly unrelated topic is focused on the multi-loop structure of propagator-type functions (p-functions 111Following Baikov:2018gap , by p-functions we understand ($\overline{\text{MS}}$-renormalized) Euclidean 2-point functions (that can also be obtained from 3-point functions by setting one external momentum to zero with the help of infra-red rearrangement) expressible in terms of massless propagator-type Feynman integrals also known as p-integrals.). About three decades ago, it was noticed that all contributions proportional to $\zeta_{4}=\pi^{4}/90$ mysteriously cancel out in the Adler function at three-loops Gorishnii:1990vf . Two decades later, it was shown that the four-loop contribution is also $\pi$-free and that a similar fact holds for the coefficient function of the Bjorken sum rule Baikov:2010je . There is by now mounting evidence, see, e.g., Baikov:2016tgj ; Herzog:2017ohr ; Davies:2017hyl ; Moch:2018wjh , that various massless Euclidean physical quantities demonstrate striking regularities in terms proportional to even $\zeta$-function values, $\zeta_{2n}$, e.g., to $\pi^{2n}$ with $n$ being a positive integer 222Notice also that, within a Schwinger-Dyson equation approach in fixed dimension, renormalized Euclidean massless correlators were shown to be expressed only in terms of odd zeta-values Kreimer:2006ua .. Additional cancellations of $\pi^{2n}$ terms have been observed in the so-called C-scheme Jamin:2017mul . Such puzzling facts have recently given rise to the “no-$\pi$ theorem”. The latter is based on the observation Broadhurst:1999xk ; Baikov:2010hf that the $\varepsilon$-dependent transformation of the $\zeta$-values:

[TABLE]

eliminates even zetas from the expansion of four-loop p-integrals. A generalization of (1) to 5- and 6-loops is available in Refs. [Baikov:2018wgs, ; Georgoudis:2018olj, ] and [Baikov:2018gap, ], respectively (the latter used the results obtained in Lee:2011jt ) 333After the submission of our paper, a new paper Baikov:2019zmy appeared on the ArXiv, where the results (1) have been extended to 7 loop order.. The results (1) and their extensions in Georgoudis:2018olj ; Baikov:2018gap ; Baikov:2019zmy give a possibility to predict the terms $\sim\pi^{2n}$ in higher orders of perturbation theory (see their evaluation in Baikov:2018wgs ; Baikov:2018gap ; Baikov:2019zmy ). Note that the results Baikov:2018wgs ; Georgoudis:2018olj ; Baikov:2018gap ; Baikov:2019zmy also contain multi-zeta values the consideration of which is beyond the scope of the present study.

In the present paper, we shall use the LKF transformation in order to study general properties of the coefficients of the propagator. We will show how the transformation naturally reveals the existence of the hatted transcendental basis. Moreover, it will allow us to extend the results of Eq. (1) to any order in $\varepsilon$.

The appearance of the hatted transcendental basis from the LKF transformation can be naturally understood in the following way. The LKF transformation produces all-loop results for very restricted objects: the difference of fermion propagators in two gauges. So, at every order of the $\varepsilon$-expansion these all-loop results should contain (at least, a part of) the basic properties of the corresponding master integrals, i.e., the all-loop results should be expressed in the form of (at least, a part of) the corresponding hatted $\zeta$-values. In a sense, it is not the full set of the hatted $\zeta$-values but only the one-fold ones. This comes from the fact that the results produced by the LKF transformation contain only products of $\Gamma$-functions and, thus, their expansions contain only the simple one-fold $\zeta$-values.

The paper is organized as follows. Sec. II introduces the LKF transformation in coordinate space. The corresponding results in momentum space are presented in Secs. III and IV for spinor QED and its scalar counter-part, respectively. Sec. V then provides the evaluation of the basic elements of the LKF transformation in the form of hatted $\zeta$-values. The results are summarized in Sec. VI and some of the first coefficients of the expansion of the hatted $\zeta$-values in terms of the usual ones are displayed in App. A.

II LKF transformation

In the following, we shall consider QED in an Euclidean space of dimension $d$ ($d=4-2\varepsilon$). The general form of the fermion propagator $S_{F}(p,\xi)$ in some gauge $\xi$ reads:

[TABLE]

where the tensorial structure, e.g., the factor $\hat{p}$ containing Dirac $\gamma$-matrices, has been extracted. It is also convenient to introduce the $x$-space representation $S_{F}(x,\xi)$ of the fermion propagator as:

[TABLE]

The two representations, $S_{F}(x,\xi)$ and $S_{F}(p,\xi)$, are related by the Fourier transform which is defined as:

[TABLE]

The famous LKF transformation connects in a very simple way the fermion propagator in two different gauges, e.g., $\xi$ and $\eta$. In dimensional regularization, it reads:

[TABLE]

where

[TABLE]

Note that, in dimensional regularization, the term $D(0)$ is proportional to the massless tadpole $T_{2}$, the massive counterpart of which is defined as:

[TABLE]

The tadpole $T_{\alpha}(m^{2})\sim\delta(\alpha-d/2)$ in the massless limit and, thus, $D(0)=0$ in the framework of dimensional regularization. So, Eq. (5) can be simplified as follows:

[TABLE]

We may now proceed in calculating $D(x)$. In order to do so, it is possible to use the following simple formulas for the Fourier transform of massless propagators (see, for example, Kotikov:2018wxe ):

[TABLE]

This yields:

[TABLE]

or, equivalently, with the parameter $\varepsilon$ made explicit:

[TABLE]

From Eq. (11), we see that $D(x)$ contributes with a common factor $\Delta A$ accompanied by the singularity $\varepsilon^{-1}$.

Note that, in the case of the scalar QED, the corresponding scalar propagator $S_{C}(x,\xi)$ has a similar transformation between $\xi$ and $\eta$ gauges, i.e.,

[TABLE]

III LKF transformation in momentum space

Let’s assume that, for some gauge fixing parameter $\eta$, the fermion propagator $S_{F}(p,\eta)$ with external momentum $p$ has the form (2) with $P(p,\eta)$ reading:

[TABLE]

In Eq. (13), the $a_{m}(\eta)$ are coefficients of the loop expansion of the propagator and $\tilde{\mu}$ is the renormalization scale:

[TABLE]

which lies somehow between the MS-scale $\mu$ and the $\overline{\text{MS}}$-scale $\overline{\mu}$. Then, the LKF transformation shows that, for another gauge parameter $\xi$, the fermion propagator can be expressed as:

[TABLE]

where

[TABLE]

In order to derive (16), we used the fermion propagator $S_{F}(p,\eta)$ with $P(p,\eta)$ given by (13), did the Fourier transform to $S_{F}(x,\eta)$ and applied the LKF transformation (8). As a final step, we took the inverse Fourier transform and obtained $S_{F}(p,\xi)$ with $P(p,\xi)$ given by (15).

III.1 Scale fixing

In our present study, we consider only the case of the so-called MS-like schemes (see, for example, $\tilde{\mu}$ above in (14)). In such schemes, we need to fix specific terms coming from the application of dimensional regularization. Such a procedure will be called scale fixing and will play a crucial role in our analysis.

Let’s first recall that the $\overline{\text{MS}}$-scale $\overline{\mu}$ is related to the previously defined scale $\tilde{\mu}$ with the help of:

[TABLE]

where $\gamma$ is the Euler constant. An advantage of the $\overline{\text{MS}}$-scale is that it subtracts the Euler constant $\gamma$ from the $\varepsilon$-expansion. Moreover, it is well known that, in calculations of two-point massless diagrams, the final results do not display any $\zeta_{2}$  444Strictly speaking, $\zeta_{2}$ can appear in some formulas such as sum rules in deep-inelastic scattering. They originate from an analytic continuation Kotikov:2005gr of certain special forms of p-integrals. We will not consider this case in the present study.. So it is convenient to choose some scale which also subtracts $\zeta_{2}$ in intermediate steps of the calculation. For this purpose, we shall consider two different scales.

The first one is the popular $G$-scale Chetyrkin:1980pr , which subtracts the coefficient in factor of the singularity $1/\varepsilon$ in the one-loop scalar p-type integral, i.e.,

[TABLE]

Following Broadhurst:1999xk , we shall use a slight modification of this scale that we will refer to as the $g$-scale and in which an additional factor $1/(1-2\varepsilon)$ is subtracted from the one-loop result, i.e.,

[TABLE]

The advantage of the $g$-scale (over the $G$-scale) will reveal itself in discussions below related to the so-called transcendental weight of various contributions.

We shall also introduce a new scale which is based on old calculations of massless diagrams performed by Vladimirov who added Vladimirov:1979zm an additional factor $\Gamma(1-\varepsilon)$ to each loop contribution. The latter corresponds to adding the factor $\Gamma^{-1}(1-\varepsilon)$ to the corresponding scale. We shall refer to this scale as the minimal Vladimirov-scale, or MV-scale, and define 555Notice that the form (20) has been used once to define the $\overline{\text{MS}}$ scheme (see Errata to Ref. Kataev:1988sq ).:

[TABLE]

As we will show below, the use of the MV-scale leads to simpler results in comparison with the $g$ one. Hence, the MV-scale is more appropriate to our analysis and all our basic results will be given in the MV-scale. After that we will discuss the differences coming from the use of the $g$-scale.

In both the MV-scale and $g$-scale, we can rewrite the result (16) in the following general form:

[TABLE]

where $p=$MV$,g$ see also Sec. IV for very similar results in the case of scalar QED. In Eq. (21), the factor $(1-(m+1)\varepsilon)/(1-(m+l+1)\varepsilon)$ has been specially extracted from $\Phi_{p}(m,l,\varepsilon)$ in order to insure equal transcendental level, i.e., the same value of $s$ for $\zeta_{s}$ at every order of the $\varepsilon$-expansion of $\Phi_{p}(m,l,\varepsilon)$ (see below). Central to the present work, the factors $\Phi_{\text{MV}}(m,l,\varepsilon)$ and $\Phi_{g}(m,l,\varepsilon)$ read:

[TABLE]

and may be expressed as expansions in $\zeta_{i}$ $(i\geq 3)$ that we now proceed on studying.

III.2 MV-scale

The $\Gamma$-function $\Gamma(1+\beta\varepsilon)$ has the following expansion:

[TABLE]

Substituting Eq. (23) in Eq. (22a), yields for the factor $\Phi_{\text{MV}}(m,l,\varepsilon)$:

[TABLE]

where

[TABLE]

and, as expected from the MV-scale, we do have:

[TABLE]

As can be see from Eq. (24), $\Phi_{\text{MV}}(m,l,\varepsilon)$ contains $\zeta_{s}$-function values of a given weight (or transcendental level) $s$ in factor of $\varepsilon^{s}$. Such a property strongly constrains the coefficients of the $\varepsilon$-series thereby simplifying our analysis. It is reminiscent of the one earlier found in Ref. Kotikov:2000pm . When judiciously used, it sometimes allows to derive results without any calculations (as in Ref. Kotikov:2002ab ). In other cases, it simplifies the structure of the results which can then be predicted as an ansatz in a very simple way (see Refs. Fleischer:1998nb ; Kotikov:2007cy ). For a recent application of such property, see the recent papers Dixon:2019uzg and references and discussions therein.

IV Scalar QED

In the case of scalar (spin-[math]) QED, the LKF transformation leads to expressions which are very similar to spinor QED that we have considered so far.

Indeed, let $S_{C}(p,\eta)$ be a scalar propagator with external momentum $p$ and gauge fixing parameter $\eta$; its general expression reads:

[TABLE]

where $a^{c}_{m}(\eta)$ are coefficients of the loop expansion of the propagator and $\tilde{\mu}$ the renormalization scale (14). Then, the LKF transformation shows that, for another gauge parameter $\xi$, the scalar propagator can be expressed as:

[TABLE]

where

[TABLE]

In order to derive (29), we used the scalar propagator $S_{C}(p,\eta)$ of Eq. (27), did the Fourier transform to $S_{C}(x,\eta)$ and applied the LKF transformation (12) which is identical to (8). As a final step, we took the inverse Fourier transform and obtained $S_{C}(p,\xi)$ in (28).

For the MV- and $g$-scales, we can rewrite Eq. (29) in the following form:

[TABLE]

where $p=$MV$,g$. We therefore see that the difference between Eqs. (30) and (21) is only in the factor $(1-(m+1)\varepsilon)/(1-(m+l+1)\varepsilon)$ which did not play any role in the above analysis. Hence, for scalar QED, we can repeat all evaluations done for spinor QED and we will end up with the same set of hatted $\zeta$-values as before in Eqs. (45) and (46).

V Solution of the recurrence relations

We now focus on the polynomial $p_{s}(m,l)$ of Eq. (25) that is conveniently separated in even and odd $s$ values. Then, we see that the following recursion relations hold:

[TABLE]

Specific to the MV-scheme, these relations only depend on $L$ which leads to strong simplifications. Nevertheless, they are difficult to solve for arbitrary $k$. It is simpler to proceed by explicitly considering the first values of $k$:

[TABLE]

showing that $p_{s}$ takes the form of a polynomial in $L$ in factor of $p_{3}$. Then, taking $Lp_{3}$ from (32b) and substituting it in (32c) yields:

[TABLE]

which reveals that the even polynomial $p_{6}$ can be entirely expressed in terms of the lower order odd ones, $p_{3}$ and $p_{5}$.

We may automate this procedure for higher values of $k$. The general expression of $p_{s}$ is given by:

[TABLE]

see Eqs. (52) of App. A for values of the first 20 coefficients and Eqs. (53) for values of $p_{s}$ with $7\leq s\leq 20$. Taking $L^{k}p_{3}$ from the equations for $p_{2k-1}$ and substituting them in the equations for $p_{2k}$ yields:

[TABLE]

where the 20 first coefficients, $C_{2k,2s-1}$, are presented in App. A, Eq. (54). From these results, it is possible to determine the exact $k$-dependence of $C_{2k,2s-1}$, which has the following structure:

[TABLE]

with the first coefficients $b_{2m-1}$ taking the values:

[TABLE]

Examining the numerators of $b_{2m-1}$, one can see that they are proportional to the numerators of Bernoulli numbers. Indeed, a closer inspection reveals that, accurate to a sign, the coefficients $b_{2m-1}$ coincide with the zero values of Euler polynomials $E_{n}(x)$:

[TABLE]

and therefore to Bernoulli and Genocchi numbers, $B_{m}$ and $G_{m}$, respectively, because

[TABLE]

Hence, the compact formula for the coefficients $b_{2m-1}$, expressed through the well known Bernoulli numbers $B_{m}$, reads:

[TABLE]

Together with (36), Eq. (40) provides an exact analytic expression for $p_{2k}$, Eq. (35), for arbitrary values of $k$.

V.1 Hatted $\zeta$-values

At this point, it is convenient to represent the argument of the exponential in the r.h.s. of (24) as follows:

[TABLE]

With the help of Eq. (35), the first term in the r.h.s. of Eq. (41) may be expressed as:

[TABLE]

Then, Eq. (41) can be written as $\sum_{s=2}^{\infty}\,\hat{\eta}_{2s-1}\,p_{2s-1}\,\varepsilon^{2s-1}$ where

[TABLE]

Thus, Eq. (24) can be represented as:

[TABLE]

where

[TABLE]

with

[TABLE]

Together with (46) and (40), Eq. (45) provides an exact expression for the hatted $\zeta$-values in terms of the standard ones valid for all $\varepsilon$, see Eqs. (55) for explicit expressions of $\hat{\zeta}_{2s-1}$ for $s\leq 10$.

V.2 $g$-scale

We may proceed in a similar way for the factor $\Phi_{g}(m,l,\varepsilon)$. Substituting Eq. (23) in Eq. (22b), yields:

[TABLE]

where the new polynomial $p^{g}_{s}(m,l)$ can be expressed in terms of $p_{s}(m,l)$, Eq. (25), as:

[TABLE]

Eq. (48) is such that $\delta_{s}(m,l)=0$ for $s=1$ and $s=2$ and, thus,

[TABLE]

similarly to the Vladimirov case, Eq. (26).

We may then consider the even and odd values of $s$ separately as in (31) leading to the following recursion relations:

[TABLE]

These recurrence relations depend on the variable $l$ but not on the product $L=l(l+1)$ as it was for the MV-scale. So, the $g$-scale recursion relations (50) are essentially more complicated than the MV-scale ones, (31). Fortunately, it is very simple to see that in the relations:

[TABLE]

the coefficients $C_{2k,2s-1}$ are exactly the same as in Eq. (35) because the corrections $\delta_{2k}$ and $\delta_{2k-1}$ exactly cancel each other. So, the hatted $\zeta$-values for the $g$-scale are identical to the ones of the MV-scale, (45) and (46).

VI Summary

From the result (21) corresponding to the LKF transformation of the fermion propagator (and similarly, from (29), corresponding to the LKF transformation of the scalar propagator) we have found peculiar recursion relations (31a) and (31b) between even and odd values of the polynomial associated to the uniformly transcendental factor $\Phi_{\text{MV}}(m,l,\varepsilon)$ (22a). These relations are simple in the new MV-scheme that we have introduced in Eq. (20). They relate the even and odd parts in a rather simple way (see Eq. (35)) which reveals the possibility (44) to express all results for $\Phi_{\text{MV}}(m,l,\varepsilon)$ in terms of hatted $\zeta$-values. As can be seen from Eqs. (50a) and (50b) in the more popular $g$-scheme, the corresponding recursion relations are slightly more complicated but lead to the same relations (35) between even and odd parts of the polynomial associated to $\Phi_{g}(m,l,\varepsilon)$ (22b) and, correspondingly, to the same hatted $\zeta$-values. Our careful study of the recursion relations (31a) and (31b) allowed us to derive exact formulas, Eqs. (43) and (45), relating hatted and standard $\zeta$-values to all orders of perturbation theory. The coefficients of the relations are expressed trough the well-known Bernoulli numbers, $B_{2m}$, as can be see from Eqs. (46) and (40). The numerical values of some of these coefficients and some explicit relations between hatted and standard $\zeta$-values are presented in App. A. Our results provide stringent constraints on multi-loop calculations at any order in perturbation theory.

Acknowledgements.

We are grateful to M. Bellon, A. Kataev, R. Lee and S. Mikhailov for valuable comments and discussions.

Appendix A Explicit expressions for $A_{k,m}$, $p_{s}$, $C_{2k,2s-1}$ and hatted $\zeta$-values

A.1 Coefficients $A_{k,m}$ of Eq. (34) for $k\leq 20$

[TABLE]

A.2 Polynomials $p_{s}$ of Eq. (34) for $7\leq s\leq 20$

[TABLE]

A.3 Coefficients $C_{2k,2s-1}$ of Eq. (35) for $k\leq 10$

[TABLE]

A.4 Hatted $\zeta$-values, $\hat{\zeta}_{2s-1}$, of Eq. (45) for $s\leq 10$

We display below the first terms of the $\varepsilon$-series generated by our exact result Eq. (45). For the sake of clarity, we display in a box the terms which were known up to the present work from the latest publication Baikov:2018gap . Moreover, we also display in a dashed box the terms calculated in the paper Baikov:2019zmy which appeared in ArXiv two months after the first version of our paper.

[TABLE]