Eur. Phys. J. C (2016) 76:180 THE EUROPEAN , _, CrossMark

D0110.1140/epjc/s10052-016-4022-z PHYSICAL JOURNAL C

Regular Article - Theoretical Physics

Renormalization on noncommutative torus

D. D'Ascanio1a, P. Pisani1b, D. V. Vassilevich23c

1 Instituto de Física La Plata-CONICET, Universidad Nacional de La Plata, C.C. 67 (1900), La Plata, Argentina

2 CMCC, Universidade Federal do ABC, Santo André, SP CEP 09210-180, Brazil

3 Department of Physics, Tomsk State University, Tomsk, Russia

Received: 13 February 2016 / Accepted: 14 March 2016

© The Author(s) 2016. This article is published with open access at Springerlink.com

Abstract We study a self-interacting scalar y4 theory on the d-dimensional noncommutative torus. We determine, for the particular cases d = 2 and d = 4, the counterterms required by one-loop renormalization. We discuss higher loops in two dimensions and two-loop contributions to the self-energy in four dimensions. Our analysis points toward the absence of any problems related to the ultraviolet/infrared mixing and thus to renormalizability of the theory. However, we find another potentially troubling phenomenon which is a wild behavior of the two-point amplitude as a function of the noncommutativity matrix 0.

1 Introduction

One of the motivations for considering quantum field theories on noncommutative spaces was the hope that they may be ultraviolet (UV) finite. It was shown, however, that UV divergences persist on the noncommutative (NC) Moyal plane [1,2]. Moreover, though certain Feynman diagrams are less UV divergent than in the commutative case, they develop singularities at some special, typically zero, value of the external momenta. When such diagrams appear as subgraphs of higher-order diagrams, the latter diagrams become divergent in a nonrenormalizable manner. This phenomenon [3-5], called the UV/IR mixing [4], is the main obstacle to renormalization of NC field theories.

It was believed for some time that the UV/IR mixing appears exclusively in Euclidean signature spaces. However, it was demonstrated [6] that similar problems exist in Minkowski spacetime as well.

Various methods were proposed to deal with this problem. Of course, the supersymmetry helps to achieve renormalizability of noncommutative theories [7,8]. Grosse and

a e-mail: dascanio@fisica.unlp.edu.ar b e-mail: pisani@fisica.unlp.edu.ar c e-mail: dvassil@gmail.com

Wulkenhaar [9,10] motivated by the Langmann-Szabo duality [11] proposed to add to the action an oscillator term which breaks translation invariance but ensures renormalizability. Modifications of the momentum dependence of the kinetic term were considered in [12]. Taking the noncommutativity parameter nilpotent [13] also improves renormalization. It has been shown [14] that spontaneous symmetry breaking softens the UV/IR mixing. A fairly recent review is Ref. [15].

In this work, we take a different path. We consider a noncommutative y4 theory on a torus. Sensitivity of UV divergences in NC theories to the presence of compact dimension (and even eventual disappearance of such divergences) has been stressed already in [2]; see also [16]. Note, however, that, due to a different implementation of noncommutativity, the existence of a compact dimension in the two-dimensional case considered in [2] guarantees the finiteness of tadpole contributions. This is not the case in the model considered in the present article, where quantum corrections are UV-divergent and must be properly renormalized.

One may get an idea on the structure of counterterms, singularities of the propagators etc. by looking at the heat kernel expansion (see, e.g. [17]). Roughly speaking, the relevant operators1 on noncommutative spaces are generalized Laplacians that contain gauge fields and potentials (as usual Laplacians), but these gauge fields and potentials act by left or right Moyal multiplications on the fluctuations Sy. If the generalized Laplacian contains only left or only right Moyal multiplications, the structure of corresponding heat kernel coefficients is very simple on both NC torus [18] and NC plane [19]—they look almost as their commutative counterparts with star products instead of usual products. For interesting theories, however, the relevant operators contain both left and right multiplications. For such operators on the NC plane the structure of heat kernel coefficients is very compli-

1 For bosonic theories, these are the operators L appearing in the second variation of classical action, S2 = J(Sy)L (y)(Sy), with Sy being a fluctuation and y—a background field.

Published online: 01 April 2016

Springer

cated [20,21] thus reflecting the presence of the UV/IR mixing. The situation changes drastically on NC torus [22]. If the noncommutativity parameter satisfies the so-called Dio-phantine condition or is rational, the heat kernel coefficients (and thus the one-loop counterterms) assume a very simple form if written in terms of a suitably defined trace operation on the algebra of smooth functions on the torus. We shall use this observation to formulate our proposal for a (presumably) renormalizable <4 theory on NC torus.

Let us stress that the notion of locality does not make much sense in noncommutative theories since the star product itself is nonlocal. Instead of local polynomial actions one has to use traces of polynomials constructed from the fields and their derivatives. There are more different traces on Td than on Rd. This observation will be crucial for our construction of admissible counterterms.

Here we like to mention several papers that considered quantum field theories on NC torus. In Ref. [7] it was demonstrated that supersymmetric Yang-Mills theory on T3 with rational 6 is one-loop renormalizable. Pure Yang-Mills theories were considered in [23] at one loop. Some arguments regarding the higher-order behavior were also presented. Relations between NC theories on Td with rational 6 and matrix models were studied in [24,25].

The purpose of this paper is to set up the stage for renor-malization on NC torus and to discuss basic features of this process. First, we write down the model and introduce new counterterms for Diophantine and rational 6. We analyze in detail two- and four-point functions at one loop. In d = 2, the only superficially divergent diagrams are the one-loop two-point functions. We demonstrate that the insertion of these diagrams (together with counterterms) into internal lines of other diagrams does not lead to any divergences, so that there is no UV/IR mixing (at least in its classical formulation [4]) on T2. In d = 4, we analyze the two-loop self-energy diagrams. All our findings, though do not contain a complete proof, strongly suggest that the introduction of new coun-terterms does make the < 4 theory on NC torus in d = 2 and d = 4 renormalizable.

The counterterms depend in a very essential way on the number theory nature of 6. But not only this, we show that also renormalized two-point functions (too) strongly depend on 6. More precisely, we compare the two-point functions in d = 2 for two close values 61 and 62 of the noncommutativ-ity matrix, one being rational, and the other irrational (Dio-phantine). We find that the typical variation of the two-point function is ~ln ||61 - 621|. However, this does not necessarily mean that the theory has no prediction power. We discuss the implications and possible ways out in the Conclusions of the article.

The paper is organized as follows. The next section contains the definitions that will be used throughout the text. In Sect. 3 we consider the two-point functions at one-loop order

and analyze their renormalization and variation with 0. Section 4 is dedicated to four-point functions at one loop. Higher loops in d = 2 are considered in Sect. 5 and two-loop two-point functions in d = 4 in Sect. 6. Our results are discussed in Sect. 7. The behavior of some double sums is analyzed in Appendices A and B.

2 The model

As a base manifold we take the d-dimensional noncommu-tative (NC) torus Td with unit radii; see [26]. The algebra A of smooth functions on Td is formed by the Fourier-type series

V = X V P UP ■

where the Fourier coefficients vP e C vanish at | p| faster than any power of p. The unitary Up satisfy

up uq = en pe q U.

where 6 is a constant and non-degenerate skew-symmetric d x d matrix. Expressions such as p6q represent the quadratic form 6^v p^qv. One may think of Up's as of plane waves eipx. Then the well-known Moyal product

(p * f)(x) = exp (-Ín 31) p(x)f(y)\y=x reproduces (2),

eípx * eíqx = eíx peq ei (P+q)x

One should keep in mind, however, that the Moyal star product has to be understood as a formal expansion in the noncommutativity parameter, i.e., it is not convergent.

There is a trace on the algebra C ) defined through

T(p) =

which can be implemented in A by t(Up ) = 8p .2

To proceed, we need some number theory preliminaries concerning the matrix elements of 6. Let us define the set

Ze = {q e Zd/eq e Zd}.

Clearly, Z0 is a Z-linear space whose dimension is the rank of the rational part of 0 .As we will see, field modes vP with momentum p e Z0 present a distinct renormalization

We use &p to denote 1 if p = 0, and 0 otherwise.

behavior in the sense that they are affected differently by quantum corrections.

On the other hand, it has been demonstrated in [22] that the heat kernel expansion and thus the one-loop divergences in a wide range of quantum field theories on the NC torus are well under control if the "irrational" part of 6 satisfies a certain Diophantine condition; namely, there should be two positive constants, C and j, such that

inf \вq - к\ >

for all q e Zd\Z6. In the last section of this article, we will see that this Diophantine condition becomes crucial for the determination of the divergences of the double sums corresponding to some two-loop Feynman diagrams.

Our starting point is then the following action for a self-interacting scalar particle on the NC torus:

S[p] = 2 т(др dp) + 1 m2 r(p2) + X т(р4).

All products in (8) are in the noncommutative algebra, i.e., they are star-products. Since these will be the only products used in this work, we shall never write the symbol * explicitly. The scalar field < undergoes a self-interaction given by the four-point vertex which in Fourier space can be written as

X X! s*1+-"+*4e

i л(к,1в к2+кзв k4)

Рк1 Рк2 Рк3 <Рк4.

к1,...,к4

The free propagator is given by

< p P p')free =

The heat kernel analysis of [22] suggests that the theory may be renormalized by adding the counterterm action3

Sc.t. =

^ \~2 (p)p (p)-p + x1 (p)p (p3)-p

+ X2 (P2)p (p2)-p^

(in addition to counterterms for the couplings in (8) and eventual renormalization of the field p).

< p P p') =

S p + p '

p2 + m2 + p)

where £ (p)—the self-energy of the scalar particle—is given by the contributions of one-particle irreducible (1PI) two-point functions,

1PI diagrams = -

p + p '

( p2 + m 2)2

Ъ( p).

We will analyze the perturbative structure of the self-energy,

Ъ(p) = h E1(p) + h2 E2(p) + ■■■ ,

with particular emphasis in d = 2 and d = 4, in order to determine the kind of counterterms required by renormalization.

One-loop contributions £i(p) to the self-energy ^(p) arise from all (connected) contractions between two external fields and the fields in the vertex (9). In the commutative case all such contractions would give the same contribution because the vertex is invariant under any permutation of the internal momenta k1,..., k4. However, this invariance is lost in the presence of the twisting factor exp iж(к1вk2 + k30k4), which is only invariant under cyclic permutations of the internal momenta so there are three sets of four equivalent contractions. Since there are only two external fields, two of these sets of contractions give the same contribution due to momentum conservation. There are thus eight contractions which give the same contribution and a different type of contribution from the other four contractions. In terms of Feynman diagrams, the former are related to planar diagrams (Fig. 1) whereas the latter correspond to nonplanar ones (Fig. 2). It is well known in noncommutative theories that the distinction between planar and nonplanar contributions plays a crucial role in the description of quantum corrections to any n-point function; the NC torus is not an exception to this fact.

As a consequence, E1 (p) can be written as

E1(p) = 8X S1(0) + 4X S1(p),

3 One-loop renormalization of self-energy diagrams

In this section we analyze the one-loop two-point functions. Quantum corrections generate a full propagator

Fig. 1 One-loop planar contribution to the self-energy

3 These terms are certain (Dixmier-type) traces on the NC torus; see

[22]. In these sense, they generalize the trace terms in (8). They may also be interpreted as usual traces after projecting the fields to a subalgebra [27].

Fig. 2 One-loop nonplanar contribution to the self-energy

where S1 represents the sum

e2nik 8p

S1( P) = X

k2 + m2

One can easily see that S1(p) is divergent for certain values of p—determined by the numerical character of 0—so we need an appropriate definition of this series that provides a regularization of its divergences.

In this article we regularize the divergencies of Feynman diagrams by introducing an arbitrary complex power e of the free propagators [with Re(e) large enough],

k2 + m2 ^ (k2 + m2)1+e'

and then performing the analytic extension to e = 0; eventual divergencies then emerge as poles of this analytic extension. At some point, this technique can be related to dimensional regularization.

Let us study, in general, the sum

e2nik 0p

Sn (p,e) = £

{(k2 + m 2 )n }1+e'

whose analytic extension to e = 0 for n = 1 defines the expression 51 (p) given in (16). If we introduce the Schwinger proper time representation we obtain

Sn (p,€)

T(n + ne) Jo

n 2 fC

r(n + ne) jo

n+ne-1 e-tm2 ^^ e-tk2e2ni k8P

keZd -tm2

dttn(1+e)-1-§ e-tm2 ^ e-2- |k+8p|2 .

In the last line of this expression we have used Poisson resum-mation,

£ f (k) = (2n)dY^ f(2nk),

keZd keZd

for f (k) = exp (-tk2 + 2nik0p) and f its Fourier transform. It is convenient to consider separately the case in which p e Z0, as defined in (6); recall that for rational 0 this set is infinite, whereas for irrational 0 the set Z0 is trivial. For p e Z0 each term in the sum of expression (19) decreases exponentially for t ^ 0 so the integration can be performed in the vicinity of e = 0 and the result reads

2nn d n ^ K2-n(2nmlk + 0pI) Sn(p,e) = ^^ m2-n £ 2

(n - 1)! + O (e),

Ik + 8 p| 2 -n

where K represents the modified Bessel function. The sum S1(p), originally defined in (16), is then given—for n = 1

and any p e Z0 —by the convergent series in the r.h.s. of (21).

On the contrary, if the external momentum p belongs to the set Z0, then the term in the series (19) with k = —0p does not present the exponential decrease for small t so the integration must be performed for Re(e) > —1 + d/2n. If we separate this term we get, after integration in t,

Sn ( p,e)

T(n + ne - f )

d ^„d-2n-2ne

r(n + ne)

2n n+ne m 2

S - n -ne

T(n + ne)

k=-8 p

Kn+ne-2 (2nm |k + 8p|)

|k + 8 p| 2

— f —ii

This expression shows that, for p e Z0, the analytic extension of Sn(p, e) has a simple pole at e = 0 if n < d/2. In particular, for n = 1 we obtain

S1(p, e) = -- (-1)2 Vd md-2- + - (-1)2 Vd m 2 e 2

x I log m2 - f (I) - y J

d , „ Kd__ 1(2nm|k + 8p|) + 2n m 2-1 °

k = -8p |k + 8P| 2

d _ i n 1

+ o (e),

where Vd is the volume of the sphere Sd 1. Therefore, the original sum (16) can be written, for p e Z0, as

S1(p) = -- (-1)2 Vdm

- + (finite terms). (24) e

In conclusion, S1(p) is conditionally convergent for p / Z8

(19) but diverges as ~ md-2/e otherwise, in particular for p = 0.

The divergent contribution of 51(0) to Ei(p) (see (15)) can be removed by an ordinary mass redefinition (see (12)),

22 m ^ m

(-n)2 m

However, due to the term S1 (p) in (15), S1 (p) might still be divergent if the external momentum p belongs to Z0 so we need to introduce new mass terms in the action

X \vp I

M2 = m

r (2) «

for those field components yp such that p e Z0. This is one of the counterterms present in expression (11).

In consequence, after appropriate O (k) mass renormal-izations, S1(p) is finite for any value of p. Note that, upon

quantum corrections, the mass of the field takes a different value for field components with momentum in Z6. In particular, for irrational 6 the new term (26) in the action can be written as

1 ß2 т(р)т(р),

discussed below in Sect. 7. Note that since 61 = 62 both two-point functions, Sf (p)6l and Sf (p)62, are always finite.

One may find some similarities between this situation and the one in the matrix model approach to noncommutivity, where the effective action behaves quite irregularly for some relations between parameters of the theory [28].

so only the zero-momentum component щ of the field gets a different mass.

Although the one-loop correction to the self-energy for any value of the external momentum is rendered finite by the mass renormalizations, the sum of all these contributions— implicit in the effective action—can be seen to be convergent, for irrational в, only under the Diophantine condition [22].

Having computed these corrections, we want to analyze the dependence of two-point functions with the numerical character of в. Let us then consider two noncommutativ-ity matrices, вх and в2, one being rational while the other -irrational. Even though the difference ||вх - в21| may be arbitrarily small, the counterterms vary drastically from вх to в2. This large variation is pretty harmless if it can be removed from the amplitudes by a finite renormalization of couplings. Let us see if this is the case at the example of the one-loop two-point function in d = 2. Let вх be rational, and в2 be a Diophantine noncommutativity matrix very close to вх. The planar diagram does no depend on в, so that we shall consider nonplanar contributions only. Let us take p e Zei \{0}. If apart from the pole term in (27) one allows for a finite renormalization of д2, the finite part of Sx(p, e) may be shifted to an arbitrary p-independent value. Therefore, the renormalized two-point function reads

SR (pM = S1 + 2n £ K0 (2л m\к + в1 p\), к=-в1 p

where s1 has to be fixed by a suitable normalization condition. Since 62 is Diophantine, S1(p, e)62 is not divergent, and its renormalized value is just the e ^ 0 limit of (21),

SR(p)e2 = 2nK0(2nm\(e1 - в2)p\)

+ 2n ^ K0(2nm\к + в2p\), к=-в1 p

where we separated one of the terms in the infinite sum in (21). Consider Sf (p)61 —Sf (p)62 inthelimit ||6i —621| ^ 0. The contributions of k = —6\p cancel in this limit, as one can easily see. sx may depend on 61, but definitely not on 62. Therefore,

Sf (p)61 — Sf (p)62 = —2n K o(2n m | (61 — 62) p\) + O (1) = 2n ln (61 — 62)p^ O(1). (31)

Hence, the variation of two-point function grows indefinitely as 62 approaches 61. Some implications of this result will be

4 One-loop renormalization of four-point functions

In order to complete the analysis of one-loop divergencies we consider the four-point function with external momenta pb p2, p3, p4. The contributions of the different Feynman diagrams to the s -channel ( p1; p2 entering the same vertex) are given by

64 X2 eni(plвp2+pзвp4) L(p1 + p2, 0), (Fig. 3) (32)

64 X2 eni(plвp2+pзвp4) L(p1 + p2, px + pi), (Fig. 4) (33) 32 X2 eni(plвp2+pз[L(p1 + p2, p2)+L(p1 + p2, p4)],

(Fig. 5) (34)

32 X2 eni(plвp2+pз[L(px + p2, pi)+L(px + p2, ps)],

(Fig. 6) (35)

16 X2 eni(plвp2+pзвp4) L(px + p2, px + pa), (Fig. 7)

16 X2 eni(pl()p2-pзвp4) L(px + p2, px + p3), (Fig. 8).

In these expressions L (p, q) is defined as the analytic extension to e = 0 of

L ( p, q,e) = X

e2niквq

кеТА

{[(к + p)2 + m2](к2 + m2)}1+e '

The two terms in (34) and (35) correspond to the cases where the incoming momentum p1 + p2 enters the diagrams from the left or from the right.

In order to study the analytic extension of the sum (38) we introduce Feynman parameters u,v to collect both propagators into a single denominator, we use the Schwinger proper time representation and then the Poisson resummation formula; the result reads

e2ni квq Г(2 + 2e)

L ( p-q'e) = £ " Г2(1 + eW„A

S(u + v - 1) (uv)e

{к2 + 2ukp + up2 + m2}2+2e

n 2 r1 e rœ

r2n , , du [u(1 - u)]e dt Г2(1 + e) J0 J0

x ti_2 x

1-2+2e

-t [m2+u(1-u) p2 ] 2

e-Zf- \к+вq\2-2лiup(k+вq)

кеТ

As before, if q / Z0 then each term in the series is exponentially decreasing for small t so it can be integrated in some neighborhood of e = 0; L(p, q) is thus finite for q / Z0. On the other hand, for q e Z0, integration in t gives

L(p, q,e) =

(-n)2 f

r d -1) L 1

2n d-2

du [m2 + u (1 - u )p2]2

- log [m2 + u(1 - u) p2]

+ log Vu(1 - u) + y + f (2 - 1)

+ 2n2 I du [m2 + u(1 — u)p2]4 J0

_^ e—2niup(k+0 q)

k=—0q |k + 0q |2 2 Kd/2—|k+0 q m 2 + u(1—u) p2j + O (e).

The sums L( p, q) that determine the contributions of the diagrams displayed in Figs. 3, 4, 5, 6, 7, and 8 can then be written, for q e Z0, as

2n d-1

L( p, q) =

(-n)2 f

r (2 -1) L

du (m2 + u(1 - u)p2)2

x--+ (finite terms).

This expression is finite for d = 2 (with a branch cut at p2 = 4m2) and diverges as n2/2e (independently of p) for d = 4. In higher dimensions the residue depends on p.

Fig. 3 Planar contribution to the four-point function

Fig. 4 Nonplanar contribution to the four-point function

Fig. 5 Nonplanar contribution to the four-point function

Fig. 7 Nonplanar contribution to the four-point function

Fig. 8 Nonplanar contribution to the four-point function

Let us therefore analyze the counterterms that are needed in four dimensions in order to remove the resulting divergencies of the four-point functions. The contribution (32)— corresponding to the planar diagram in Fig. 3—contains an UV divergence, which can be removed by a renormalization of the self-coupling constant,

X ^ X ^ 1 + 4n2X

Besides, contributions (34) and (35)—corresponding to the diagrams in Figs. 5 and 6—together with the t- and u-channels are also divergent if any of the incoming momenta belongs to the set Ze. This type of divergence can be removed by introducing the following self-interaction, corresponding to the second term in (11):

X1 <k 5k+k1+k2+k3 <k1 <Pk2 <Pk3

keZ8 k1,k2 ,k3eZd ein k18k2+in k38k

X1 = 8n 2X2 —.

For irrational 8, this new interaction reads

X1 r(<p)r(<p ).

Lastly, contributions (33), (36) and (37)—corresponding to the diagrams in Figs. 4,7 and 8—present a divergence whenever the sum of two incoming momenta belongs to Ze whose cancellation requires the following self-interaction, corresponding to the third term in (11):

X2 X 2 Äk1+k2 -k $k1 +

keZ8 k1 ,...,k4

••+k4 <k1 Pk2 Pk3 Pk4e

in k1<8k2+in k3&k

Fig. 6 Nonplanar contribution to the four-point function

X2 = 6n2 X2 —.

For irrational e, the counterterm (46) reads X2 x(cp2)x(cp2).

After the introduction of these counterterms four-point functions in T4 are rendered finite for any value of the external momenta. Note that all j -functions associated with the coupling constants X,X1,X2 are positive.

5 Higher loops at two dimensions

Before analyzing higher order of perturbation series on T2 let us briefly recall the UV/IR mixing problem on noncom-mutative plane Rd. The nonplanar diagrams on Rd behave better in the ultraviolet than their commutative counterparts since (p6)—1 (with p being an external momentum) serves as an effective ultraviolet cutoff. However, the divergences are restored in the commutative limit, 6 ^ 0, implying also a singularity at p ^ 0. According to Ref. [4], these singularities cause troubles with the convergence of loop integrals at zero momenta if nonplanar diagrams are inserted into internal lines of other diagrams. Note that in two dimensions 1PI diagrams are at most logarithmically divergent, so that the IR singularity may also be at most ln | p| . This singularity is rather mild. Thus one does not expect much troubles with the UV/IR mixing in d = 2. For this reason, our consideration of the two-torus will also be rather sketchy.

Turning to T2, we first note that there are only two diagrams, Figs. 1 and 2, which are superficially divergent. By using the expression (23) and basic properties of K0, one can easily show that after adding the counterterm from (26) the nonplanar diagram Fig. 2 with p e Z6 becomes a bounded function of p. For p e Z6, there may be a growing contribution to S1(p), which comes from the momentum kp in (22) that minimizes ^ + 6p|. It reads

2nKo(2nmlkp + 0pi) ~ -2n ln Ikp + 0pi.

By the Diophantine condition (7), this contribution is restricted by 2n(1 + j) ln |p| at large |pi Therefore, the renormalized 2-point function on T2 has a logarithmic singularity, but in contrast to R2 this singularity is UV rather than IR. The UV singularities on the quantum plane, discussed in [29], were found more severe than on the commutative plane. However, the singularities on T2 are very mild. Indeed, if one inserts the renormalized diagram of Fig. 2 into an internal line with the momentum p of some other diagram one gets (at large | p| ) a multiplier of ln | p| from the diagram itself and (p2 + m2)—1 from an extra propagator. The overall contribution behaves as ln |p| ■ (p2 + m2)—1 and does even improve the convergence of larger diagram.

We saw that in the < 4 theory on T2 (i) all superficially divergent diagrams can be renormalized by the counterterms that we have proposed, and (ii) the insertion of renormalized superficially divergent diagrams does not lead to any problems with convergence. Hence, there is no UV/IR mixing in this model, and it is likely renormalizable.

6 Two-loop self-energy at four dimensions

In this last section we describe the diagrams that contribute to S2(p)—the second order correction to the self-energy— in the four-dimensional torus. For corresponding analysis on R4 one may consult Ref. [30]. In this section we restrict ourselves to the case of a pure irrational (Diophantine) noncom-mutativity parameter. Therefore, Z6 = {0}. We analyze the divergences of two-loop diagrams and point out the main difficulties one finds in computing the remaining double sums; in the course of this analysis we will see the importance of the Diophantine condition on the matrix 6. Let us also remark that, since we are interested in the renormalizability of two-point functions, we will neglect divergent contributions which are either independent or quadratic in the external momentum p for they can be removed by mass and field renormalizations of order O (X2).

Before considering two-loop diagrams, we analyze the O(X2) contributions of the counterterms already introduced in the previous sections, i.e. one-loop diagrams built with the leading quantum corrections to the parameters m2, X, as well as with the new parameters /2,X1,X2.In the first place, the nonplanar tadpole in Fig. 2 gives an O (X2) contribution,

- 32 X2m2 S2( p)—.

from the insertion of the mass correction (25) into its internal propagator, as well as another O (X2) contribution,

16 X2 Sx(p)

from the insertion of the X correction (42) into its vertex. In (50) S2(p) is defined as the analytic extension to e = 0 of (18) (for n = 2 and d = 4). Second, a planar tadpole of the type shown in Fig. 1 at the vertex X2 gives

p2 + m2 e

Equations (50)-(52) represent nonlocal (and not tracelike) divergencies introduced in the self-energy by the renor-malization of the parameters in the action; contributions of /j2 and X1 are not taken into account for they are either p-independent or quadratic in p.

The two-loop diagrams can be built in two different ways: either by inserting the one-loop self-energy S1(p) into the internal propagator of a planar or a nonplanar tadpole (see Figs. 9, 10, respectively), so that both external legs enter the same vertex, or by attaching each external leg to a different vertex so that the two loops share a common internal momentum, as in Figs. 11, 12, 13, and 14.

For the first type of diagrams, the insertion into a planar tadpole, Fig. 9, gives a contribution which, though divergent, does not depend on the external momentum. As we have

Fig. 9 Planar diagram with a £1(p) insertion

Fig. 10 Nonplanar diagram with a E1 (p) insertion

Fig. 11 Planar diagram contributing with -16Л2 U(p, 0) to E2(p)

Fig. 12 Nonplanar diagram contributing with -32Л2 U(p, p) to E2(p)

Fig. 13 Nonplanar diagram contributing with -32Л2 V(p, 0) to E2(p)

Fig. 14 Nonplanar diagram contributing with -16Л2 V(p, p) to E2(p)

already explained above, such contributions are harmless and will be discarded. The contributions corresponding into the insertion into a nonplanar tadpole, Fig. 10, are given, up to O (I2), by

e2nikep

- 4Я X ^ ■ E1(k)•

(k2 + m 2)2

Replacing (15) into this expression we obtain

- 32Л2 S1(0)S2(p) - 16Л2 T(p),

where T(p) is defined as the analytic extension to e = 0 of the series

T (p,e) =

2n ik 0l

k,lsZ4

{[(k + p)2 + m2](l2 + m2)2}1+e '

Note that the first term in (54), though p-dependent, is exactly canceled by (50).

Lastly, the contributions of diagrams which contain overlapping divergencies read (see Figs. 11, 12, 13, 14),

-16Л2 U(p, 0) - 32Л2 U(p, p) - 32Л2 V(p, 0) - 16À2 V(p, p),

where the sums U, V are defined as the analytic extensions to e = 0 of the series

U(p, q,€)

k,l sZ4

V(p, q,€)

k,l sZ4

e2n i l в q

{[(k+p)2+m2](l2+m 2 )[(l+k )2+m 2)]}1+^

e2ni ke q e2n i kel

{[(k+p)2+m2](l2+m 2 )[(l+k )2+m 2)]}1+e'

Since U (p, 0) corresponds to the diagram of the ordinary commutative case, its divergence is a quadratic polynomial in p. According to Appendix A, the sum U(p, p, e) behaves as

- 32Л2 U(p, p, e) = -16 Л2 S1(p) —

+ quad. pol. + O (e),

where we have represented by "quad.pol." terms which, though eventually divergent, are quadratic expressions4 in p. Therefore, the nonlocal non-trace-like divergence introduced by U(p, p) completely cancels (51).

In consequence, two-loop renormalization of the self-energy demands that the remaining nonlocal divergencies (given by (52), the second term in (54) and the last two terms in (56)) cancel among each other. In other words, the remaining potentially divergent terms read

(53) 48 Л2

1 n2 2

-2--16Л2 T(p,e)

p2 + m2 e

-32Л2 V(p, 0, e) - 16Л2 V(p, p, e).

The divergences of this expression at e = 0 must repeat the structure of quadratic (in p) counterterms. I.e., they have to be of the form of a quadratic polynomial in p plus a term proportional to Sp.

The double sums in (60) can all be treated in a unified way. The divergences of these sums in Z8 at e = 0 arise from the fact that the denominator increases only with a sixth power at infinity. Nevertheless, the twisting factor e2nikel contributes to regularize the series. This certainly happens in the continuum case where the corresponding integration in R8 is finite. However, in the discrete case there exist four-dimensional subspaces—with null measure in R8—for which the twisting factor vanishes.

4 More precisely, "quad. pol." should have the form a + bp2, i.e. the

coefficient in front of pM pv has to be proportional to the unit matrix.

Therefore, the divergencies of T(p. e) can be attributed to the subseries in the subspaces for which k01 = 0 and, simultaneously, the denominator increases at infinity with a power which is less or equal than four. Such subseries correspond to l = 0 and k = 0. Thus,

T(p, €) = X

{(l2 + m 2 )2( p2 + m 2 )}1+e

{((k + p)2 + m2)2m4}1+e

+ O(1). (61)

The second sum in (61) does not depend on p, while the first one is proportional to S2(0, e); see (18). We conclude that

T (p,e) = —.-2 — + quad. pol. + O (1). (62)

p2 + m2 2e

The same formula is reobtained in Appendix B, by using mathematically rigorous methods.

To analyze V (p, q, e), let us first change the summation variables, so that

V(p, q,e)

J2ni k61

= Y _e_

, {[(k+p)2 +m2][(l — q)2 +m2][(l +k — q)2 +m2)]}^

Using the same argument as above, potentially divergent terms come from the subsets k = 0, l = 0, and k + l = 0, which contribute as

( p2 + m2)1+e

52(0, e) +

(q 2 + m 2 )1+e

L(p + q, 0, e)

(q2 + m 2)1+e

L(p - q, 0,e).

The remaining terms in the series (63) are expected to decrease for large k, l as long as \k — в11 does not decrease too fast, which is guaranteed by the Diophantine condition on в. This implies in particular

У (p, 0,e) = -2-2— + quad. pol. + O(1), (65)

p2 + m2 2e

У (p, p,€) = 3 -2 — + quad. pol. + O(1). (66)

p2 + m2 2e

Unfortunately, we cannot reconfirm (65) and (66) by more rigorous methods.

From Eqs. (62), (65), and (66) one concludes that (60) does not contain any divergences except for the ones that can be removed by a renormalization of couplings in the action

SM = 1 т(дф дф) + 1 m2 r(cp2) + X x(cpA) 1

+ 2 m2 t((p) t(<p) + т((р) r((p3)

+ X2 Т(р( )т(р2).

Let us recall that in this section we assumed that в is irrational.

Although our analysis is far from being exhaustive, we believe it strongly suggest that the p4 theory on T^ is renor-malizable.

7 Conclusions

In this paper, we analyzed the renormalization of a scalar field theory on Tde with a quartic self-interaction after the introduction of a new type of nonlocal (but trace-like) interactions suggested by the previous heat kernel calculations [22]. At one loop our analysis is complete. We also argued that in two dimensions no problems appear at higher orders as well. In four dimensions, we checked the renormalization of self-energy at two-loop order relying on our understanding of the behavior of double sums, which we were able to reconfirm by rigorous methods for all diagrams but one. Our findings strongly suggest that the p4 theory on T2 and T4 is renormalizable. The renormalization always strongly depends on the Diophantine character of the noncommutativ-ity matrix в. We cannot exclude completely that some more elaborate multiple-trace counterterms will be needed, though their algebraic nature is less clear than that of the ones listed (11). To check this, one has to calculate the two-loop four-point functions.

On the technical side, it is important to develop the theory of regularized multiple sums with twisting factors. To the best of our knowledge, such sums have not been considered in the mathematics literature so far (see, e.g. [31]).

While calculating the renormalized two-point functions we encountered a potentially troubling phenomenon: these functions depend too strongly on the noncommutativity matrix. In other words, an arbitrarily small error in в may cause an arbitrarily large variation in the two-point functions. Or, the value of two-point functions cannot be predicted unless we know в with an infinite precision. This does not necessarily mean, however, that the theory is meaningless. We can suggest the following explanations and ways to overcome the difficulty.

1. Since the plane waves Up do not commute even classically, see (2), they probably do not form a good basis. Therefore, the correlation functions of plane waves may be of little physical relevance by themselves. The problem is then to find a physically motivated basis of states that will ensure a kind of "smooth" dependence of the correlation functions on в.

2. One can try to achieve a meaningful answer by smearing the correlation functions over a small vicinity of a given в. The key issue is to find an appropriate measure.

3. Finally, perhaps one can extend the model to fix в sharply to certain value, e.g. by some topological considerations.

Although we cannot offer much details of any of the items above, we believe that these directions deserve further study.

Note that for p e/ Z0 the only nonlocal divergence in this expression is contained in the first term, so that

Acknowledgments DVV was supported in part by FAPESP, Project 2012/00333-7, CNPq, Project 306208/2013-0 and by the Tomsk State University Competitiveness Improvement Program. DD acknowledges support from CONICET and UNLP (Proj. 11/X615), Argentina. PP acknowledges support from CONICET (PIP 1787/681), ANPCyT (PICT-2011-0605) and UNLP (Proj. 11/X615), Argentina.

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP3.

Appendix A: The double sum U(p, p)

Let us consider the analytic extension to e = 0 of the sum

U(p, p, e)

e2ni lep

k,l sZ

{[(k+p)2+m2](l2+m 2 )[(l+k )2+m 2)]}1+e

_^ e2n i l ep

= Z r(l_L p)2 . m2-11+e L(l, 0' e)'

[(l + p)2 + m2]

Using (39) at d = 4, this expression reads

r2 e2n ile p

U(p, p, e) =

T2(1 + e) [(l + p)2 + m2]1+e

x / du [u(1 - u)]e J0

x i dtt-1+2e e-t(m2+«(1-«)l2)

x Ze '

n- k2-2niulk

If we separate the term corresponding to k = 0, integrate in t and expand about e = 0 we obtain

U(p, p,e) = n2 T(2e) £

e2ni l ep

lsZ4 {(l + p)2 + m2}1+e

f ' (u(1 - u))e x du ■

/o (m2 + u (1 - u )l2)2e

+ 2n 2Z

e2ni l ep

—2n iulk

l (l+p)2+m2t=0J0

x K0(2n|kWm2 + u(1 - u)l2) + O(e)

U(p, p,e) = — S1 (p) + quad. pol. + O(e). 2e

In particular, for irrational 0 this expression holds for any p = 0 and its divergent part is canceled by (51); the contribution U(0, 0) is, of course, regularized by the renormalization of the parameter /x2.

Appendix B: The double sum T(p)

The double sum

T (p, e) = X

k,lsZ4

can be written as T ( p,e) =

2n ik el

{[(k + p)2 + m2](l2 + m2)2}1+e

( p2 + m 2 )1+e

S2(0, e)

k =0 l

2nike l

{[(k + p)2 + m2](/2 + m2)2}1+e ■

We will show that the divergent part of the second term in this expression does not depend on p or, equivalently, that the expression

a2n ik el

(l2 + m2)2(1+e)

J(k + p)2 + m2]1+e

(k2 + m2)

m2) 1+e J

is finite at e = 0. In order to do that, we expand (B.7) in Taylor series in p keeping potentially divergent terms only,

a2n ikel

(l2 + m2)2(1+e)

x . --(¿±4fi- + 2(1 + ,e)(2 + py\ . (B.8)

(k2 + m2)2+e

(k2 + m2)3+e

Assuming 0 is proportional to the standard four-dimensional symplectic matrix times a constant5 we can replace (k ■ p)2 by 1 p2k2, which makes (B.8) convergent.

Finally, we show that the difference between expressions (B.7) and (B.8) is finite at e = 0. In fact, after Poisson inversion this difference reads

5 In other words, we assume that the NC torus has two sets of noncom-(A.3) muting coordinates with the same noncommutativity parameter.

2n 2 XX K C (2n |l + 0 k |)

(k + p)2 + m2 k2 + m2 (k2 + m2)2

4(k ■ p)2

(k2 + m2)3

where the expression in square brackets is O (|k | 6) for large |k|. Now, for any given k one can separate in the l-sum the term corresponding to the smallest value of Ц+вk |. Under the Diophantine condition, this term decreases as ^|-6 log ^l whereas the remaining terms are bounded by e-2n\l+ek|, due to the behavior of K0 for large arguments. The subsequent sum in k is therefore absolutely convergent.

References

1. T. Filk, Divergencies in a field theory on quantum space. Phys. Lett. B 376, 53 (1996)

2. M. Chaichian, A. Demichev, P. Presnajder, Quantum field theory on noncommutative space-times and the persistence of ultraviolet divergences. Nucl. Phys. B 567,360 (2000). arXiv:hep-th/9812180

3. I. Chepelev, R. Roiban, Renormalization of quantum field theories on noncommutative Rd. 1. Scalars. JHEP 0005, 037 (2000). arXiv:hep-th/9911098

4. S. Minwalla, M. Van Raamsdonk, N. Seiberg, Noncommutative perturbative dynamics. JHEP 0002, 020 (2000). arXiv:hep-th/9912072

5. I.Y. Aref'eva, D.M. Belov, A.S. Koshelev, Two loop diagrams in noncommutative phi**4(4) theory. Phys. Lett. B 476, 431 (2000). arXiv:hep-th/9912075

6. D. Bahns, The ultraviolet infrared mixing problem on the noncommutative Moyal space. arXiv:1012.3707 [hep-th]

7. M.M. Sheikh-Jabbari, Renormalizability of the supersymmetric Yang-Mills theories on the noncommutative torus. JHEP 9906, 015 (1999). arXiv:hep-th/9903107

8. H.O. Girotti, M. Gomes, V.O. Rivelles, A.J. da Silva, A consistent noncommutative field theory: the Wess-Zumino model. Nucl. Phys. B 587, 299 (2000). arXiv:hep-th/0005272

9. H. Grosse, R. Wulkenhaar, Renormalization of phi**4 theory on noncommutative R**2 in the matrix base. JHEP 0312, 019 (2003). arXiv:hep-th/0307017

10. H. Grosse, R. Wulkenhaar, Renormalization of phi**4 theory on noncommutative R**4 in the matrix base. Commun. Math. Phys. 256, 305 (2005). arXiv:hep-th/0401128

11. E. Langmann, R.J. Szabo, Duality in scalar field theory on noncommutative phase spaces. Phys. Lett. B 533, 168 (2002). arXiv:hep-th/0202039

12. R. Gurau, J. Magnen, V. Rivasseau, A. Tanasa, A translationinvariant renormalizable non-commutative scalar model. Commun. Math. Phys. 287, 275 (2009). arXiv:0802.0791 [math-ph]

13. R. Fresneda, D.M. Gitman, D.V. Vassilevich, Nilpotent noncom-mutativity and renormalization. Phys. Rev. D 78, C25CC4 (2CCS). arXiv:CSC4.i566 [hep-th]

14. F. Ruiz Ruiz, UV/IR mixing and the Goldstone theorem in noncommutative field theory. Nucl. Phys. B 637, i43 (2CC2). arXiv:hep-th/02020ll

15. D.N. Blaschke, E. Kronberger, R.I.P. Sedmik, M. Wohlgenannt, Gauge theories on deformed spaces. SIGMA 6, C62 (2CiC). arXiv:iCC4.2i27 [hep-th]

16. M. Chaichian, A. Demichev, P. Presnajder, M.M. Sheikh-Jabbari, A. Tureanu, Quantum theories on noncommutative spaces with nontrivial topology: Aharonov-Bohm and Casimir effects. Nucl. Phys. B 611, 3S3 (2CCi). arXiv:hep-th/0l0l209

17. D.V. Vassilevich, Heat kernel expansion: user's manual. Phys. Rep. 388, 279 (2CC3). arXiv:hep-th/C3C6i3S

iS. D.V. Vassilevich, Noncommutative heat kernel. Lett. Math. Phys. 67, iS5 (2CC4). arXiv:hep-th/C3iCi44

i9. V. Gayral, B. Iochum, The spectral action for Moyal planes. J. Math. Phys. 46, C435C3 (2CC5). arXiv:hep-th/C4C2i47

2C. D.V. Vassilevich, Heat kernel, effective action and anomalies in noncommutative theories. JHEP 0508, CS5 (2CC5). arXiv:hep-th/C5C7i23

21. R. Bonezzi, O. Corradini, S.A. Franchino Viñas, P.A.G. Pisani, Worldline approach to noncommutative field theory. J. Phys. A 45, 4C54Ci (2Ci2). arXiv:i2C4.iCi3 [hep-th]

22. V. Gayral, B. Iochum, D.V. Vassilevich, Heat kernel and number theory on NC-torus. Commun. Math. Phys. 273, 4i5 (2CC7). arXiv:hep-th/C6C7C7S

23. T. Krajewski, R. Wulkenhaar, Perturbative quantum gauge fields on the noncommutative torus. Int. J. Mod. Phys. A 15, iCii (2CCC). arXiv:hep-th/99C3iS7

24. Z. Guralnik, R.C. Helling, K. Landsteiner, E. Lopez, Pertur-bative instabilities on the noncommutative torus, Morita duality and twisted boundary conditions. JHEP 0205, C25 (2CC2). arXiv:hep-th/C2C4C37

25. G. Landi, F. Lizzi, R.J. Szabo, Matrix quantum mechanics and soliton regularization of noncommutative field theory. Adv. Theor. Math. Phys. 8(i), i (2CC4). arXiv:hep-th/C4CiC72

26. A. Connes, C * algèbres et géométrie différentielle. C. R. Acad. Sci. Paris 290, 599-6C4 (i9SC). arXiv:hep-th/CiCiC93

27. D.V. Vassilevich, Induced Chern-Simons action on noncommutative torus. Mod. Phys. Lett. A 22, i255 (2CC7). arXiv:hep-th/C7CiCi7

2S. D.N. Blaschke, H. Steinacker, M. Wohlgenannt, Heat kernel expansion and induced action for the matrix model Dirac operator. JHEP 1103, CC2 (2Cii). doi:iC.iCC7/JHEPC3(2Cii)CC2. arXiv:iCi2.4344 [hep-th]

29. M. Chaichian, A. Demichev, P. Presnajder, Quantum field theory on the noncommutative plane with E(q)(2) symmetry. J. Math. Phys. 41, i647 (2CCC). arXiv:hep-th/99C4i32

3C. A. Micu, Noncommutative phi4 theory at two loops. JHEP 025, CiCi (2CCi). arXiv:hep-th/CCCSC57

3i. S. Paycha, Regularized Integrals, Sums and Traces (AMS, Providence, 2Ci2)