Arxiv:1805.08566V1 [Hep-Th] 22 May 2018 Euaiainpoeue Hsrglrzditgaini Ultraviolet Is Integration Mom Regularized Are This Which the Parameter, As Procedure

Home , Regularization (physics)

arXiv:1805.08566v1 [hep-th] 22 May 2018 euaiainpoeue hsrglrzditgaini ultraviolet is integration mom regularized are This which the parameter, as procedure. regularization regularization such dimensional parameters the on masses, Villars depend integrals regularized These oeflto odsrb h yaiso ytm n processe and systems properties. of fractal and/or dynamics memory the equatio long-range integro-differential describe of to theory tool The powerful physics. p and by mechanics represented are derivatives these space momentum In type. 3,adi ensaqatmfil hoya o-rva ii ftheory of limit non-trivial a as theory field quantum a proces defines This cut-off. theory. it the and diverge of the [3], parameters d the ”renormalizable”, short of very called redefinition at are singular theories which sp field theories, the momentum quantum of in used structure of arise the are divergences i.e. which These moment, theories, high 2]. field [1, quantum divergence of ultraviolet feature characteristic A Introduction 1 iuul itiue ytm n lsia ed.Uigtefractio the Using fields. classical and systems f distributed [10, tinuously for temperature gravity positive and at theory theory field field quantum quantum fractional the the 8], [7, mechanics quantum eetytesailfatoa-re eiaie aebe activ been have derivatives fractional-order spatial the Recently h eomlzto eursterglrzto ftelo integrals loop the of regularization the requires renormalization The rcinlcluu n rcinldffrnileutos[,5 hav 5] [4, equations differential fractional and fractional-order calculus on Fractional based procedure regularization a suggest We rcinlcluu losu otk noacutfatoa power fractional account into take to us allows calculus Fractional euaie opitgasdpn nprmtrta steo the is that parameter non-i on of depend derivatives integrals fractional loop Riesz regularized the of application Regularizati on operators. fractional-order by replaced are eiaie h euaiainpoeuei demonstrated is procedure regularization ϕ The derivative. 4 ter on -theory nti ae,w rps e euaiain hr integer where regularization, new propose we paper, this In rcinlDrvtv euaiaini QFT in Regularization Derivative Fractional n dmninlped-ulda space-time. pseudo-Euclidean -dimensional kblsnIsiueo ula Physics, Nuclear of Institute Skobeltsyn oooo ocwSaeUniversity, State Moscow Lomonosov -al [email protected] E-mail: aiyE Tarasov E. Vasily ocw199,Russia 119991, Moscow Abstract 1 nfrqatmfil hoisbased theories field quantum for on tgrodr ssgetd The suggested. is orders nteger o clrmsls ed in fields massless scalar for rder odrdffrniloperators differential -order ihpwrlwnon-locality, power-law with s scle h renormalization the called is s wrfntoso momentum. of functions ower 11]. l sdi h space-fractional the in used ely so o-nee resi a is orders non-integer of ns a acls ecnconsider can we calculus, nal sacs ntenro class narrow the In istances. > α finite. cscnb eoe ya by removed be can nces nhg-nrypyis is physics, high-energy in sdi h corresponding the in used nu u-ff h Pauli- the cut-off, entum atoa pc-ie[9], space-time ractional eiaie fteRiesz the of derivatives ieapiainin application wide a e lwnnoaiyo con- of nonlocality -law c rmmdso very of modes from ace ftefractional the of 0 nmmnu space. momentum in ihaultraviolet a with space-time fractional differential equations in the quantum field theory. For the first time, the fractional-order Laplace and d’Alembert operators have been suggested by Riesz in [12]. Then non-integer powers of d’Alembertian are considered in different works (for example, see Section 28 in [4] and [13, 14, 15]). The fractional Laplace and d’Alembert operators of the Riesz type are the basis for the fractional field theory in multidimensional spaces. As it was shown in [16, 17, 18, 19, 20], the continuum equations with fractional derivatives of the Riesz type can be directly connected to lattice models with long-range properties. A connection between the dynamics of lattice system with long-range properties and the fractional continuum equations are proved by using the transform operation [16, 17, 18, 19]. There are different definitions of fractional derivatives such as Riemann-Liouville, Caputo, Grünwald-Letnikov, Marchaud, Weyl, Sonin-Letnikov, Riesz [4, 5]. Unfortunately all these fractional derivatives have a lot of unusual properties. For example, the well-known Leibniz rule does not hold for differentiation of non-integer orders [26]. It should be noted that the use of the fractional derivative of non-integer order is actually equivalent to using an infinite number of derivatives of all integer orders. For example, the RL α Riemann-Liouville derivative Da+ can be represented in the form of the infinite series ∞ Γ(α + 1) (x − a)k−α dkf(x) ( RLDα f)(x)= (1) a+ Γ(k + 1)Γ(α − k + 1) Γ(k − α + 1) dxk Xk=0 for analytic (expandable in a power series on the interval) functions on (a, b), (see Lemma 15.3 of [4]). Therefore, an application of fractional derivatives means that we take into account a contribution of derivatives of all integer orders with the power-law weight. In the suggested FD-regularization we use fractional derivatives of non-integer orders α close to integer values. The use of these fractional derivatives of non-integer orders α close to integer values in the regularization procedure means that we consider a weak nonlocality, i.e. a slight deviation from the locality.

2 Fractional Laplacian for Euclidean space Rn

2.1 Fractional integration in the Riesz form The Riesz fractional integral of order α for Rn is deﬁned [12] (see also Section 25 of [4]) by equation R α 1 α−n I f(P )= f(Q) rPQ dQ, (2) H (α) Rn n Z where P and Q are points of the space Rn, and πn/22αΓ(α/2) H (α)= . (3) n Γ((n − α)/2) Let us give some properties of the Riesz fractional integral that are proved in [12]:

2 1. The semi-group property of the Riesz fractional integration

RIα1 ( RIα2 f(P )) = RIα1+α2 f(P ), (4) where α1 > 0, α2 > 0, α1 + α2 < n. 2. The action of the Laplace operator on the Riesz fractional integral

∆ RIα+2f(P )= − RIαf(P ), (5) and we also have ∆k RIα+2kf(P )=(−1)k RIαf(P ), (k ∈ N). (6) 3. The Riesz fractional integration of exponential function does not change this function

RIα eixj = eixj , (7)

− n α/2 n n R α i P =1 aj xj i P =1 aj xj I e j = ajxj e j . (8) j=1 ! X There is the following important property

RIαf(P )=(−1)m RIα+2m ∆mf(P ). (9)

This property allows us to use an analytic continuation of RIαf(P ) for negative values of α > −2m, where m ∈ R (see [12]). In this case, the semigroup property (4) can be used for α1 > −2m, α2 > −2m, α1 +α2 > −2m, and m ∈ R. Analogously, we can consider the property (6) in the form ∆k RIαf(P )=(−1)k RIα−2kf(P ), (k ∈ N). (10)

2.2 Fractional Laplacian in the Riesz form For the ﬁrst time fractional Laplace operators have been suggested by Riesz in [12] (see also Section 25 of [4]). The fractional Laplacian (−∆)α/2 in the Riesz form for n-dimensional Euclidean space Rn can be considered [4] as an inverse Fourier’s integral transform F −1 of |k|α by ((−∆)α/2ϕ)(x)= F −1 |k|α(Fϕ)(k) , (11) where α> 0 and x ∈ Rn. For α> 0, the fractional Laplacian of the Riesz form can be deﬁned [4] as the hypersingular integral α/2 1 1 m 4 (−∆) ϕ (x)= α+n (∆z ϕ)(z) d z, (12) d (m, α) Rn |z| n Z

3 m where m>α, and (∆z ϕ)(z) is a finite difference of order m of a field ϕ(x) with a vector step z ∈ Rn and centered at the point x ∈ Rn:

m m! (∆mϕ)(z)= (−1)j ϕ(x − jz). z j!(m − j)! j=0 X

The constant dn(m, α) is defined by π1+n/2A (α) d (m, α)= m , n 2αΓ(1 + α/2)Γ((n + α)/2) sin(πα/2) where m m! A (α)= (−1)j−1 jα. m j!(m − j)! j=0 X Note that the hypersingular integral (12) does not depend on the choice of m>α. The Fourier transform F of the fractional Laplacian is given by F{(−∆)α/2ϕ}(k) = |k|α(Fϕ)(k). This equation is valid for the Lizorkin space [4] and the space C∞(Rn) of infinitely differentiable functions on Rn with compact support.

2.3 Fractional Laplacian in the Riesz-Trujillo form Using the property (9), we can deﬁne the fractional Laplace operator by the equation

RT ∆α/2f(P )=(−1)m RI2m−α ∆mf(P ), (13) where m > α/2 and m ∈ R. This form of deﬁnition of the fractional Laplacian has been suggested by J.J. Trujillo in 2012 (see also [25]). It is important that the semi-group property holds for these operators

RT ∆α1/2 RT ∆α2/2f(P )= RT ∆(α1+α2)/2f(P ), (14) where 0 <α1 < 2m, 0 <α2 < 2m, α1 + α2 < 2m. We can see that RT ∆kf(P )=(−1)k ∆kf(P ), (15) where k ∈ N. In the Riesz’s notation [12], this equation means that

RI−2kf(P )=(−1)k ∆kf(P ). (16)

We also have a generalization in the form

RIα−2kf(P )=(−1)k ∆k RIαf(P ). (17)

Note that the value of m can be chosen suﬃciently large to fulﬁll all the conditions on the fractional order α in the property (14).

4 3 Fractional d’Alembertian for pseudo-Euclidean space- Rn time 1,n−1

For the ﬁrst time the fractional d’Alembertian has been suggested by Riesz in [12] (see also [25]). Note that in the Riesz paper [12], the d’Alembertian is denoted by ∆ and it is called Rn for the Laplace operator for 1,n−1. We will use the generally accepted notation. Rn For the pseudo-Euclidean space-time 1,n−1, we use

n 2 2 2 2 rQP = rPQ =(x1 − y1) − (xj − yj) , (18) j=2 X and the operator ∂2 n ∂2 = − . (19) ∂x2 ∂x2 1 j=2 j X The Riesz fractional integal is deﬁned as

R α 1 α−n I f(P )= f(Q) rPQ dQ, (20) Hn(α) Rn Z 1,n−1 where Hn(α) is deﬁned by the equations (see Eq. 20 and 88 of [12])

Γ((α +2 − n)/2)Γ(α) H (α)= π(n−2)/22α−1 Γ(α/2) Γ((α +2 − n)/2) = π(n−1)/2 . (21) n Γ((α + 1)/2)

This expression of Hn(α) does not coincide with the Euclidean case (3). To ensure the conver- gence of the integral, we should assume that the parameter α satisﬁes the condition α−n> −2, i.e. α > n − 2. There is the following important property (see Eq. 67 of Chapter 3 in [12])

RIαf(P )= RIα+2m mf(P ), (22) that allows us to use an analytic continuation of RIαf(P ) for negative values of α if

α > n − 2 − 2m.

Using the property (22), we can deﬁne the fractional d’Alembertian for n-dimensional Rn pseudo-Euclidean space-time 1,n−1 by the equation

α/2f(P )= RI2m−α mf(P ), (23) where m> (α + n − 2)/2 and m ∈ R.

5 It is important that the semi-group property

α1/2 α2/2f(P )= (α1+α2)/2f(P ), (24) holds for operators (23) if 0 <α1 < 2m+2−n,0 <α2 < 2m+2−n, and α1 +α2 < 2m+2−n. α/2 Rn The fractional d’Alembertian for n-dimensional pseudo-Euclidean space-time 1,n−1 can be considered as an inverse Fourier’s integral transform F −1 of |k|α by

(α/2ϕ)(x)= F −1 |k|α(Fϕ)(k) , (25) Rn where α> 0 and x ∈ 1,n−1.

4 Scalar ﬁeld in pseudo-Euclidean space-time

Rn Let us consider a scalar ﬁeld ϕ(x) in the n-dimensional pseudo-Euclidean space-time 1,n−1 that is described by the ﬁeld equation

+ m2 ϕ(x)= J(x), (26)

Rn where is the d’Alembert operator, ϕ(x) is a real ﬁeld, and x ∈ 1,n−1 is the space-time vector with components xµ, where µ = 0, 1, 2, ..., n. Suppose that the scalar ﬁeld ϕ(x) has a source J(x), and we have put ~ = 1. Field equation (26) follows from stationary action principle, δS[ϕ] = 0, where the action S[ϕ] has the form

S[ϕ]= dnx L[ϕ(x)] (27) Rn Z 1,n−1 with the Lagrangian 1 L[ϕ(x)] = − ϕ(x) + m2 ϕ(x). (28) 2 The solution of (28) can be represented in the form

n ϕ(x)= − ΣF (x − y) J(y) d y, (29) Z where ΣF (x − y) is the so-called Feynman propagator, obeying

2 ( + m − iε) ΣF (x − y)(x)= −δ(x). (30)

Here δ(x) is the Dirac delta-function. It is easy to see that ΣF (x − y) has the Fourier representation dnk e−ikx Σ (x − y)(x)= . (31) F (2π)n k2 − m2 + iε Z 6 4 We can consider ϕ theory with an interaction constant g. It is known that ΣF (0) is a divergent quantity, which modiﬁes the free particle propagator and contributes to the self- energy. In momentum space it corresponds to the loop integral dnq 1 g Σ (0) = g . (32) F (2π)n q2 − m2 Z There are nth powers of q in the numerator and two in the denominator, so the integral diverges quadratically at large q, i.e. it is ultra-violet divergent. Expression (32) corresponds to a one- loop diagram Feynman that has the order g. Another divergent diagram is the O(g2) graph. The corresponding expression of the Feynman integral is

dnq dnq δ(q + q − p − p ) dnq 1 g2 1 2 1 2 1 2 = g2 . (33) (2π)n (2π)n (q2 − m2)(q2 − m2) (2π)n (q2 − m2)((p + p − q)2 − m2) Z 1 2 Z 1 2 Here there are nth powers of q in numerator and the four powers in denominator. For n = 4, we have four powers of q in both numerator and denominator, so we get a logarithmic divergence.

5 Regularization

Regularization is a method of isolating the divergences. There are several techniques of regularization. The most intuitive one is to introduce a lattice reguariazation [21]. 1. The Pauli-Villars regularization modiﬁes the quadratic terms by subtraction oﬀ the same quadratic terms with a much larger mass

1 1 1 1 ∂ ϕ ∂ ϕ − m2 ϕ2 − ∂ φ ∂ φ − M 2 φ2 . (34) 2 µ µ 2 2 µ µ 2 The Pauli-Villars regularization is a contribution of another field (the Pauli-Villars field φ) with the same quantum numbers as the original fieeld, but having the opposite statistics. For example, in the ϕ4-theory, the Pauli-Villars field φ is used, and the Pauli-Villars regularization changes the propagator as

2 2 1 1 1 e−p /Λ −→ − . (35) p2 + m2 p2 + m2 p2 + m2 p2 + M 2 2. The Gaussian cutoﬀ regularization is implemented by the modiﬁed kinetic term in the form 1 2 ∂ ϕ e/Λ ∂ ϕ, (36) 2 µ µ which changes the propagator 2 2 1 e−p /Λ −→ . (37) p2 + m2 p2 + m2

7 3. The higher derivative regularization also modiﬁes the quadratic terms such as

1 1 ∂ ϕ 1 − ∂ ϕ + m2 ϕ 1 − ϕ. (38) 2 µ Λ2 µ 2 Λ2 Then the propagator is modiﬁed

1 p2 1 −→ 1+ . (39) p2 + m2 Λ2 p2 + m2 6 Regularization by fractional derivatives

We propose a new regularization, where integer-order diﬀerential operators are replaced by fractional-order operators. For example, we can replace , which is the d’Alembert operator in the ﬁeld equations, by the fractional d’Alembert operator α:

−→ α/2. (40)

Here we use the dimensionless coordinates. The fractional derivative regularization modiﬁes the quadratic terms such as 1 1 − ϕ(x) + m2 ϕ(x) −→ − ϕ(x) α/2 + m2 ϕ(x). (41) 2 2 Rn Using the fractional d’Alembertian for pseudo-Euclidean space-time 1,n−1, we can replace ﬁeld equation (26) by the equation

α/2 + m2 ϕ(x)= J(x), (42)

α Rn where is the fractional d’Alembert operator, ϕ(x) is a real ﬁeld, and x ∈ 1,n−1 is the space-time vector with components xµ, where µ =0, 1, 2, ..., n. The solution of (42) is n ϕ(x)= − ΣF (x − y) J(y) d y, (43) Z where ΣF (x − y) is the Feynman propagator, obeying

α/2 2 ( + m − iε) ΣF (x − y)(x)= −δ(x). (44)

Using (25), it is easy to see that ΣF (x − y)(x) has the Fourier representation

dnk e−ikx Σ (x)= . (45) F (2π)n (k2)α/2 − m2 + iε Z

8 In this case, ΣF (0) is a divergent quantity, which modiﬁes the free particle propagator and contributes to the self-energy. In momentum space it corresponds to the loop integral dnq 1 g Σ (0) = g . (46) F (2π)n (q2)α/2 − m2 Z There are nth powers of q in the numerator and α in the denominator. Another divergent diagram is the O(g2) graph. The corresponding expression of the Feynman integral is

dnq dnq δ(q + q − p − p ) g2 1 2 1 2 1 2 = (2π)n (2π)n ((q2)α/2 − m2)((q2)α/2 − m2) Z 1 2 dnq 1 1 = g2 · . (47) (2π)n (q2)α/2 − m2 ((p + p − q)2)α/2 − m2 Z 1 2 Here there are nth powers of q in numerator and the 2α powers in denominator.

7 FD-regularization for massless theory

For simpliﬁcation, we will consider a massless scalar ﬁeld theory. For this case, expression (33) with m = 0 has the form dnq dnq δ(q + q − p − p ) dnq 1 1 g2 1 2 1 2 1 2 = g2 · . (48) (2π)n (2π)n ((q2)α/2)((q2)α/2) (2π)n (q2)α/2 ((p + p − q)2)α/2 Z 1 2 Z 1 2 The denominators in the integrand of (33) are combined by using the Feynman’s parametric integral formula 1 Γ(a + b) 1 za−1(1 − z)b−1 = dz . (49) Aa Bb Γ(a)Γ(b) [Az + B(1 − z)]a+b Z0 Changing the variables q → p and p1 + p2 → q, the integral (48) takes the form dnp 1 1 I(q)= g2 · . (50) (2π)n (p2)α/2 ((p − q)2)α/2 Z Using formula (49), we get

1 1 1 zα/2−1 (1 − z)α/2−1 · = dz . (51) (p2)α/2 ((p − q)2)α/2 [p2 − 2pq(1 − z)+ q2(1 − z)]α Z0 As a result, we have

1 dnp 1 I(q)= g2 zα/2−1 (1 − z)α/2−1 dz . (52) (2π)n [p2 − 2pq(1 − z)+ q2(1 − z)]α Z0 Z

9 Using p =(p0, r) and the polar coordinates, we get that momentum integral of (52) is dnp I[n, α, q]= = (p2 − 2pQ + M 2)α Z 2π(n−1)/2 +∞ ∞ rn−2 dr = dp0 2 2 2 α , (53) Γ((n − 1)/2) −∞ [p − r − 2pQ + M ] Z Z0 0 where Q = q(1 − z), M 2 = q2(1 − z). (54)

This integral is Lorentz invariant, so we evaluate it in the frame qµ = (µ, 0). Then 2pQ = ′ 2µ(1 − z)p0. Changing variables to pµ = pµ − qµ(1 − z), which implies that

2 2 ′ 2 2 2 p0 − 2q(1 − z)p0 = p0 − 2µ(1 − z)p0 =(p0) − q (1 − z) , (55) we have (n−1)/2 +∞ ∞ n−2 2π ′ r dr I[n, α, q]= dp0 ′ 2 2 2 α , (56) Γ((n − 1)/2) −∞ [(p ) − r + q z(1 − z)] Z Z0 0 where we use q2(1 − z) − q2(1 − z)2 = q2z(1 − z). Using the beta function ∞ t2x−1 Γ(x)Γ(y) 2 = B(x, y)= , (57) (1 + t2)x+y Γ(x + y) Z0 we get ∞ sβ ds Γ((1 + β)/2)Γ(α − (1 + β)/2) = , (58) (s2 + M 2)α 2(M 2)α−(1+β)/2 Γ(α) Z0 where we can use 2 ′ 2 2 M = −(p0) − q z(1 − z), β = n − 2. (59) Using (58) with (59), equation (56) gives

+∞ ′ (n−1)/2 Γ(α − (n − 1)/2) dp0 I[n, α, q]= π ′ 2 2 α−(n−1)/2 = Γ(α) −∞ [−(p ) − q z(1 − z)] Z 0 +∞ ′ −α+(n−1)/2 (n−1)/2 Γ(α − (n − 1)/2) dp0 =(−1) π ′ 2 2 α−(n−1)/2 . Γ(α) −∞ [(p ) + q z(1 − z)] Z 0 Using (58), we get Γ(α − n/2) 1 I[n, α, q]= i πn/2 . (60) Γ(α) [−q2z(1 − z)]α−n/2 As a result, we obtain dnp Γ(α − n/2) 1 I[n, α, q]= = iπn/2 . (61) [p2 − 2pq(1 − z)+ q2(1 − z)]α Γ(α) [−q2z(1 − z)]α−n/2 Z 10 Using (61), expression (52) has the form Γ(α − n/2) 1 zα/2−1 (1 − z)α/2−1 dz I(q)= ig2πn/2 . (62) Γ(α) [−q2z(1 − z)]α−n/2 Z0 Then we have Γ(α − n/2) 1 1 I(q)= ig2πn/2 z(n−α)/2−1 (1 − z)(n−α)/2−1 dz. (63) Γ(α) [−q2]α−n/2 Z0 For n = 4 and α =2+ ε, we get 1 =(−1)n/2−α 1+ ε ln(q2)+ O(ε2) . (64) [−q2]α−n/2 Here we can use (−1)n/2−α =(−1)−ε = eiπε =1+ iπε. (65) We can introduce the function 1 A(n, α)= z(n−α)/2−1 (1 − z)(n−α)/2−1 dz. (66) Z0 For n = 4 and α =2+ ε, expression (66) takes the form 1 ε 1 A(4, 2+ ε)= (z(1 − z))−ε/2 dz =1 − ln(z(1 − z)) dz + O(ε2). (67) 2 Z0 Z0 Using 1 ln(z(1 − z)) dz = −2, (68) Z0 equation (67) is written as A(4, 2+ ε)=1+ ε + O(ε2). (69) For 4-dimensional space time, Γ(α − n/2) 1 = . (70) Γ(α) (α − 1)(α − 2) For n = 4 and α =2+ ε, equation (70) gives Γ(α − n/2) 1 1 = = − 1. (71) Γ(α) (1 + ε) ε ε Then integral (63) with n = 4 and α =2+ ε, has the form 1 I(q)= ig2π2 − 1 (1 + ε)(1+ iπε) 1+ ε ln(q2) . (72) ε As a result, we get ig2π2 I(q)= + ig2π2 ln(q2) − π3g2 + O(ε). (73) ε

11 Remarks For the special cases, we have π A(2, 1) = π, A(4, 1) = , A(4, 2)=1, A(3, 1)=1. (74) 8 The well-known asymptotic expression for the gamma-function is

1 1 π2 Γ(ε)= − γ + γ2 + ε + O(ε2) (ε → 0). (75) ε 2 6 For 2-dimensional space-time,

Γ(α − n/2) 1 = . (76) Γ(α) α − 1

If α =1+ ε, then Γ(α − n/2) 1 = . (77) Γ(α) ε For 4-dimensional space time (n=4),

Γ(α − n/2) 1 = . (78) Γ(α) (α − 1)(α − 2)

If α =1+ ε, then Γ(α − n/2) 1 1 = = − − 1. (79) Γ(α) ε (ε − 1) ε If α =2+ ε, then Γ(α − n/2) 1 1 = = − 1. (80) Γ(α) (1 + ε) ε ε Remark. By changing variables to p′ = p − q(1 − z), we see that the denominator in the integrand of equation (51) is the square of (p′)2 −m2 +q2z(1−z). But dnp′ = dnp, so relabelling p′ → p, expression (48) becomes

1 dnp 1 I(q)= g2 zα/2−1 (1 − z)α/2−1 dz . (81) (2π)n [p2 + q2z(1 − z)]α Z0 Z Conclusion

In this paper, we proposed a new regularization of quantum ﬁeld theories, which is based on the Riesz fractional derivatives of non-integer order α > 0. The fractional derivatives of non- integer order are characterized by nonlocality, the proposed FD regularization is actually a regularization by space-time nonlocality.

12 Note that the Riesz fractional derivatives of non-integer order [12, 25] can be used not only for FD regularization. An application of fractional calculus in quantum field theory allows us to take into account space and time nonlocality [28, 29, 10, 11, 21]. It should be noted that we can use the Riesz fractional derivatives [12, 25], the fractional Laplacian and d’Alembertian in the lattice quantum field theories [21]. Using the exact discretization of the fractional Laplacian d’Alembertian [20], we can consider exact discrete analogs [27] of the quantum field theories. The fractional derivatives of the Riesz type [12, 25] are directly connected to the long-range properties and nonlocal interactions [18, 19, 27]. Because of this, this tool can be used not only for regularization, but also for constructing new quantum theories in high-energy physics and condensed matter physics.

Appendix

For the fractional case, we should consider the integrals ∞ sγ ds K[α,β,γ,M 2]= (82) (sα + M 2)β Z0 instead of the integrals of equation (58). Using new variable z such that z2 = sα, and

z = sα/2, s = z2/α, (83)

ds = (2/α) z2/α−1 dz, sγ = z2γ/α, (84) equation (82) gives 2 ∞ z2(γ+1)/α−1 K[α,β,γ,M 2]= dz . (85) α (z2 + M 2)β Z0 Then using equation (58), we obtain

2Γ((γ + 1)/α)Γ(β − (γ + 1)/α) K[α,β,γ,M 2]= . (86) 2 α (M 2)β−(γ+1)/α Γ(β) Conﬂict of Interests

The author declares that there is no conﬂict of interests regarding the publication of this paper.

Acknowledgments

This work was partially supported by grant NSh-7989.2016.2

13 References [1] N.N. Bogoliubov, D.V. Shirkov, Introduction to theory of Quantized Fields Third edition (John Wiley and Sons, New York, 1980) 638 pages. [2] L.H. Ryder, Quantum Field Theory (Cambridge University Press, Cambridge, 1985). 487 pages. [3] J.C. Collins, Renormalization: An Introduction to Renormalization, the Renormaliza- tion Group and the Operator-Product Expansion (Cambridge University Press, Cambridge, 1984). 392 pages. [4] S.G. Samko, A.A. Kilbas, O.I. Marichev, Fractional Integrals and Derivatives Theory and Applications (Gordon and Breach, New York, 1993) 1006 pages. [5] A.A. Kilbas, H.M. Srivastava, J.J. Trujillo, Theory and Applications of Fractional Differ- ential Equations (Elsevier, Amsterdam, 2006) 353 pages. [6] V.E. Tarasov, Fractional Dynamics: Applications of Fractional Calculus to Dynamics of Particles, Fields and Media (Springer, New York, 2011). 450 pages. [7] V.E. Tarasov, ”Weyl quantization of fractional derivatives”, Journal of Mathematical Physics. Vol.49. No.10. (2008) 102112. 6 pages. (arXiv:0907.2699) [8] V.E. Tarasov, ”Quantum Mechanics of Non-Hamiltonian and Dissipative Systems (Elsevier Science, 2008) 540 pages. [9] G. Calcagni, ”Quantum field theory, gravity and cosmology in a fractal universe”, Journal of High Energy Physics. Vol.2010. March (2010) 120. 38 pages. (arXiv:1001.0571) [10] S.C. Lim, ”Fractional derivative quantum fields at positive temperature”, Physica A. Vol.363. No.2. (2006) 269-281. [11] S.C. Lim, L.P. Teo, ”Casimir effect associated with fractional Klein-Gordon field”, in the book J. Klafter, S.C. Lim, R. Metzler (Eds.), Fractional Dynamics. Recent Advances (World Scientific, Singapore, 2011). pp.483-506. (arXiv:1103.1673) [12] M. Riesz, ”L’intégrale de Riemann-Liouville et le probléme de Cauchy”, Acta Mathematica. Vol.81. No.1. (1949) 1-222. in French. [13] C.G. Bollini, J.J. Giambiagi, ”Arbitrary powers of D’Alembertians and the Huygens’ principle”, Journal of Mathematical Physics. Vol.34. No.2. (1993) 610-621. [14] D.G. Barci, C.G. Bollini, L.E. Oxman, M.C. Rocca, ”Lorentz-invariant pseudo-differential wave equations”, International Journal of Theoretical Physics. Vol.37. No.12. (1998) 3015- 3030. [15] R.L.P.G. do Amaral, E.C. Marino, ”Canonical quantization of theories containing fractional powers of the d’Alembertian operator”, Journal of Physics A. Vol.25. No.19. (1992) 5183-5202. [16] V.E. Tarasov, ”Continuous limit of discrete systems with long-range interaction”, Journal of Physics A. Vol.39. No.48. (2006) 14895-14910. (arXiv:0711.0826) [17] V.E. Tarasov, ”Map of discrete system into continuous”, Journal of Mathematical Physics. Vol.47. No.9. (2006) 092901. 24 pages. (arXiv:0711.2612)

14 [18] V.E. Tarasov, ”Toward lattice fractional vector calculus”, Journal of Physics A. Vol.47. No.35. (2014) 355204. 51 pages. [19] V.E. Tarasov, ”Lattice fractional calculus”, Applied Mathematics and Computation. Vol.257. (2015) 12-33. [20] V.E. Tarasov, ”Exact discretization of fractional Laplacian”, Computers and Mathematics with Applications. 2017. Vol.73. No.5. P.855863. DOI: 10.1016/j.camwa.2017.01.012 [21] V.E. Tarasov, ”Fractional quantum field theory: From lattice to continuum”, Advances in High Energy Physics. Vol.2014. (2014) Article ID 957863. 14 pages. [22] G. ’t Hooft, M. Veltman, ”Regularization and renormalization of gauge fields”, Nuclear Physics B. Vol.44. No.1. (1972) 189-213. [23] G. Leibbrandt, ”Introduction to the technique of dimensional regularization”, Reviews of Modern Physics. Vol.47. No.4. (1975) 849-876. [24] K. Svozil, ”Quantum field theory on fractal spacetime: a new regularization method”, Journal of Physics A. Vol.20. No.12. (1987) 3861-3875. [25] H. Prado, M. Rivero, J.J. Trujillo, M.P. Velasco, ”New results from old investigation: A note on fractional m-dimensional differential operators. The fractional Laplacian”, Frac- tional Calculus and Applied Analysis. Vol.18. No.2. (2015) 290-306. DOI: 10.1515/fca- 2015-0020 [26] V.E. Tarasov, ”No violation of the Leibniz rule. No fractional derivative”, Communications in Nonlinear Science and Numerical Simulation. Vol.18. No.11. (2013) 2945-2948. [27] V.E. Tarasov, ”Exact discretization by Fourier transforms”, Communications in Nonlinear Science and Numerical Simulation. Vol.37. (2016) 3161. [28] G.V. Efimov, ”Non-local quantum theory of the scalar field”, Communications in Mathe- matical Physics. Vol.5 (1). (1967)4256. DOI: 10.1007/BF01646357 [29] G.V. Efimov, The problems of the quantum theory of nonlocal interactions. Moscow: Nauka, 1985. 216p.