Nonlocal generalization of Galilean theories and gravity

In this paper we propose a wider class of symmetries including the Galilean shift symmetry as a subclass. We will show how to construct ghost-free nonlocal actions, consisting of infinite derivative operators, which are invariant under such symmetries, but whose functional form is not simply given by exponentials of entire functions. Motivated by this, we will consider the case of a scalar field and discuss the pole structure of the propagator which has infinitely many complex conjugate poles, but satisfies the tree-level unitarity. We will also show how the ultraviolet behavior of loop integrals can be ameliorated. Moreover, we will consider kinetic operators respecting the same symmetries in the context of linearized gravity. In such a scenario, the graviton propagator turns out to be ghost-free and the spacetime metric generated by a point-like source is nonsingular. These new nonlocal models can be seen as an infinite derivative generalization of Lee-Wick theories and open a new branch of nonlocal theories.


Introduction
Back in the fifties, when the renormalizability of quantum electrodynamics was still not totally established, physicists were also trying to improve the ultraviolet (UV) behavior of the loop integrals by introducing nonlocal form factors in the Lagrangians [1]. It was observed that integrals could be made more convergent by using non-polynomial differential operators, and infinite order derivative 1 actions with such form factors were also studied from a pure mathematical and axiomatic point of view [2].
Infinite derivative gauge theories and their UV behavior were first studied in Refs. [3,5,6], where the authors computed the propagator and noticed that the presence of exponentials of entire functions does not introduce any new poles other than the local ones. Subsequently, the same results were also obtained in the context of gravity [7][8][9], where it is possible to formulate an infinite derivative theory whose action is of higher order in the curvature invariants but still free from ghost fields, therefore unitary at the quantum level [21,[23][24][25]. Such theories are also often called ghost-free infinite derivative field theories. An important point to note here is that the infinite order derivatives make the interaction nonlocal, indeed the interaction vertex is smeared out in the momentum space.
In ghost-free infinite derivative gravity the propagator has been computed around maximally symmetric backgrounds [9,22]. It was shown that no other dynamical degree of freedom is introduced besides the massless graviton of general relativity (GR) [3,8,9,22]. At the classical level, the nonlocal gravitational interaction can solve blackhole [8,9,22,[27][28][29][30][31][32][33][34][35][36][37][38][39][40] and cosmological singularities [7,[41][42][43][44]. Infinite derivative gravity in 3 dimensions has been recently studied, in both massless and massive cases; see Ref. [45]. At a quantum level there are hints to show that the UV behavior of the theory is ameliorated [4,6,8,10,46,47]. Furthermore, in Ref. [48] it was noticed that in an infinite derivative extension of the Abelian Higgs model has no instabilities emerge as the β-function vanishes in the UV regime, i.e. for energies beyond the scale of nonlocality, p 2 > M 2 s . Infinite derivative Lagrangians were also studied in the context of thermal field theory [49][50][51], inflationary cosmology [52] and supersymmetry [53,54]. Recently, it has also been pointed out that in presence of multi particle interactions, such as mutiparticle scattering, the nonlocal scale can be transmuted from the UV to the IR depending on the number of particles involved in the scattering process. The emergence of such a new scale in the IR is an interesting result which demonstrates the existence of some complementarity principle in infinite derivative theories [55].
So far we have only mentioned nonlocal theories whose Lagrangians are constructed in terms of analytic operators. It is worthwhile to mention that nonanalytic operators like 1/ and ln( ) naturally emerge from one-loop quantum corrections to the effective action of quantum gravity [56][57][58][59], and also in causal-set theory [60].
In this paper, our aim is to understand if there exist nonlocal actions which exhibit some special symmetries. In particular, we will ask whether we can enlarge the Galilean symmetry [61], ∂ µ φ −→ ∂ µ φ + b µ , which is typical of theories with some specific derivative self-interactions and is often taken as a fundamental guiding principle to construct and constrain the interaction Lagrangians in both flat [61,62] and curved backgrounds [63][64][65][66][67][68][69]. A very interesting aspect of higher derivative theories with Galilean symmetry is that the field equations are second order in the derivatives, thus no instabilities arise [61]. However, we now wish to generalize to actions with derivatives of infinite order.
The paper is organized as follows.
• In Section 2 we will find nonlocal operators which are invariant under a larger class of transformations containing the Galilean transformations as a subclass.
• In Section 3, we will introduce a new nonlocal field theoretical model motivated by such a new symmetry; we will compute the propagator and discuss the pole structure, by showing tree-level unitarity. We will show that the UV properties of 1-loop integrals are ameliorated by the presence of derivatives of infinite order. Furthermore, we will notice that such nonlocal models can be seen as an infinite derivative generalization of the Lee-Wick theories, since the propagator possesses infinite complex conjugate poles.
• In Section 4, we will consider the same class of nonlocal operators in the context of linearized quadratic gravity, and also in this case we will compute the propagator by showing tree-level unitarity. Very interestingly, we will notice that it is possible to enlarge the class of nonlocal form factors which make the graviton propagator ghostfree at the tree-level, without restricting ourself to exponentials of entire functions. Moreover, we will show that the linearized spacetime metric generated by a point-like source is nonsingular.
• In Section 5, we will discuss our results and draw the conclusions.
• In Appendix A, we will discuss another larger class of transformations containing the shift transformations as a subclass, and consider a natural nonlocal extension of the Dirac action, which can be gauge invariant.
Throughout the paper we will work with the mostly positive metric signature (− + ++), and we use Natural Units = 1 = c.
2 Enlarging the Galilean symmetry: from local to nonlocal operators In standard local field theory, Lagrangians are constructed in terms of polynomials of fields and polynomials of derivatives of fields, therefore, the order of derivatives is always finite: where n is a positive finite integer and φ(x), in principle, can be any kind of tensorial field. It is often very important to ask what kind of transformations leave a Lagrangian (action) invariant, especially in order to find symmetries and conserved quantities. In particular, we are interested here in the Galilean shift symmetry in the Minkowski spacetime, defined by the following transformation: where {b µ , a} are five constant parameters which generate the whole Galilean symmetry.
Lagrangians invariant under the above transformations can be constructed, in terms of a massless scalar field φ(x), where only specific derivative self-interactions are present. For example, the following Lagrangian exhibits such a symmetry [61]: where Λ can be, for example, a cut-off energy scale above which an effective field theory description may break down. A more general form of the Lagrangian in Eq.(2.4) can be found in Ref. [61]. We now wish to understand whether we can enlarge the Galilean symmetry while working with nonlocal differential operators. First of all, for the sake of clarity, let us mention that a nonlocal Lagrangian is a function which can be also made up of nonpolynomial differential operators, like for instance where the non-polynomial operators contain infinite order covariant derivatives; M s is the energy scale of nonlocality at which new physics should manifest, and it is mathematically needed to make the arguments of logarithm and exponential dimensionless. For example, in the case of the exponential of the d'Alembertian, we can write the operator: where the derivative order n goes up to infinity. As also mentioned above in the Introduction, exponentials of entire functions, like the one in Eq.(2.6), have been used previously in the context of infinite derivative field theory and gravity [3,4,[6][7][8][9]14]. In order to satisfy the Galilean shift symmetry in Eqs.(2.2,2.3), we need to construct other nonlocal operators which are different from the ones already known in the literature. Further note that second order derivative operators, containing terms like the ones in Eq.(2.4), are invariant under the Galilean transformation in Eqs. (2.2,2.3). It is easy to understand that if φ is invariant under the Galilean transformation, then any power n φ will be also invariant, which in turn implies that the following nonlocal operator is also invariant under the Galilean transformation: which is slightly different from the one in Eq.(2.6).
We now wish to ask the following question -is there any function ψ(x), such that the transformation φ(x) −→ φ(x)+ψ(x), leaves invariant the nonlocal operator in Eq.(2.7), but not the local one? Or, in other words, can we find a function ψ(x) such that e − /M 2 s − 1 ψ = 0, but ψ = 0?

1D case
Let us start considering a 1-dimensional case as a warm-up exercise, namely let us find solutions to the following nonlocal differential equation: which also means to find the eigenfunctions ψ of the operator e −∂ 2 x /M 2 s of unit eigenvalue. First of all, note that the following property holds where k is an integer number and we have defined the differential operator Thus, we need to find solutions for the equation θ 2 ψ(x) = 0 : (2.11) where C k and D k are two integration constants, which need to be fixed by imposing that for k = 0, then we recover the Galilean shift symmetry in Eqs.(2.2,2.3). Indeed, by choosing we obtain which in the limit k → 0 gives ψ 0 (x) = bx + a : where b ≡ 2(1 + i)B √ πM s . The solutions in Eq.(2.13) are valid for positive integer, k > 0, while in the case of k < 0, the solutions read:

4D case
We can now apply the same procedure to the 4-dimensional case. Also in this case we will find a larger class of solutions which are eigenfunctions of the differential operator e − /M 2 s with eigenvalue equal to one, 16) and such that ψ(x) = 0. We can write where we have now defined We can show that the solutions of the differential equation θ 2 ψ(x) = 0 are given by: where c 2 ≡ c µ c µ = 1; the analog solution for k < 0 can be easily found too, as done above for the 1-dimensional case. In the limit k → 0 we obtain: which recovers the Galilean shift symmetry in Eqs.
Note that we can now rewrite the function ψ k in Eq. (2.19) in a more compact form as follows: Hence, we have explicitly shown that by working with the nonlocal, infinite derivative operator introduced in Eq.(2.7), we can enlarge the Galilean shift symmetry, which now becomes a subclass of a wider class described by the following family of parameters: Thus, the Galilean transformations correspond to the subfamily {a, b µ , 0} . In other words, the Galilean shift symmetry turns out to be a subclass (k = 0) of this larger symmetry which is expressed in terms of the following field transformation 2 :

Generic powers of
So far we have only performed our study for the nonlocal operator where the exponent is simply given by . However, we can also consider more generic entire functions in the exponent, as for example (− /M 2 s ) n or even non-polynomial entire functions. In the former general scenario, a new symmetry can be still found, but the solutions for ψ k (x) become more complicated as the power n increases.
For instance, we can consider the following more general operator: First of all, we can write In this more general case the solutions of the field equation is formally given by: where d l are constant parameters depending on the order 2n of the differential equation in Eq.(2.27) and satisfy the algebraic equation (−d 2 l ) n + 2πi = 0, while C k,l are integration constants. Anyway, in this paper we will mainly work with the nonlocal operator O 1 .

A nonlocal model with infinite complex conjugate poles
In this Section we wish to construct actions, or in other words Lagrangians, which respect the new symmetry in Eq.(2.23). By working with the nonlocal operator O 1 , defined in Eq.(2.7), possible action terms are given by with λ i being dimensionless coupling constants and the iteration is meant for n > 2. For instance, up to cubic vertex interaction, the action in Eq.(3.1) explicitly reads: is invariant under the transformation in Eq.(2.23) and, the first nonvanishing term in the expansion is given by

Nonlocal propagator
In this subsection, we will consider the simplest action possessing the symmetry in Eq.(2.23) and only work with a kinetic term, as we are mainly interested in computing the propagator. Such a kinetic action is which, in the local limit, recovers the Klein-Gordon action for a massless scalar field. The field equation can be easily found by variating the action and reads: The bare propagator is defined as the inverse of the kinetic operator, and in momentum space it is given by Note that the pole structure of the propagator is more complicated than the local one; besides the usual real massless pole at p 2 = 0, we now have an infinite number of complex conjugate poles. Indeed, the denominator in Eq.(3.9) vanishes for where k is an integer number: k = 0 corresponds to the only real massless pole, while each value of k = 0 is associated with two complex poles, whose conjugates are the ones corresponding to the opposite integer −k. In fact, the square root in Eq.(3.10) can be also decomposed in real and imaginary parts as follows: , with θ being the Heaviside step function, or equivalently ε(k) ≡ sign(k). See Fig. 1 for the graphic location of the poles. We can isolate the real massless pole in the propagator by writing: where the zeros of the function f (p) correspond to the complex poles. Indeed f (p) = 0, if and only if p 2 = i2πM 2 s k, with k = 0. Furthermore, the propagator has no essential singularities at infinity which, in principle allow us to define a Wick rotation from Minkowski to Euclidean space and vice-versa.
Note that the presence of complex poles, in general, may spoil perturbative unitarity [70], in such a way that predictability would be lost at the quantum level. For this reason, further investigations are needed in order to understand whether the unitarity condition, and so the optical theorem, are satisfied.

Unitarity with infinite complex conjugate poles
The unitarity condition of the S-matrix reads: which, by writing S = 1 + iT, can also be expressed as follows where the last equation represents the so called optical theorem. One of the implication of the optical theorem is that the imaginary part of any amplitude T can not be negative. For instance, in the case of a tree-level amplitude with constant three-vertex interaction, it requires the imaginary part of the propagator to satisfy the inequality: Im {Π(p)} ≥ 0. It is very important to understand whether the imaginary part of the propagator in Eq.(3.9) satisfies the optical theorem, i.e. whether the presence of infinite complex conjugate poles spoils unitarity or not. Let us introduce a new variable z := −p 2 /M 2 s , so that the nonlocal propagator in Eq.(3.9) will read , (3.15) whose poles are given by z = i2πk, where k is an integer number. We now wish to partially decompose the propagator in an infinite number of fractions, each corresponding to a single pole. Indeed, we can write then by using the following identity [71] 1 we obtain 3 By going back to the momentum variable p 2 , we obtain the nonlocal propagator in Eq.(3.9) expressed in partial decomposition: Therefore, we have been able to explicitly isolate all the poles of the propagators. The first term corresponds to the real massless pole, p 2 = 0, while the second term takes into account of all the infinite complex conjugate poles, p 2 = ±i2πM 2 s k with k = 0. By using this form of the propagator it is now easy to compute its imaginary part: where in the last step we have taken the limit → 0, and we have used the fact that Thus, the infinite complex poles do not contribute to the imaginary part of the propagator, and it is due to the fact that they appear in conjugate pairs. This implies that the propagator in Eqs. (3.9,3.19) satisfies the tree-level unitarity condition. It is worthwhile to mention that in Refs. [72][73][74][75][76], the authors have considered the unitarity issue with a finite number of complex conjugate poles; in particular in Refs. [73,74] it was rigorously shown that the optical theorem is preserved at all orders in pertubation theory. They have considered the so called Lee-Wick theories, whose Lagrangians are sixth order in the derivatives and the propagator is made up of the usual real pole plus a pair of complex conjugate poles. Our propagator in Eqs. (3.9,3.19) can be seen as a nonlocal extension of the Lee-Wick propagator, where we now have infinite pairs of complex conjugate poles whose infinite number is related to the presence of infinite order derivatives in the Lagrangian. In this respect, the new nonlocal model we have introduced can be seen as an infinite derivative Lee-Wick model.
Since the nonlocal propagators whose form given by the first term in Eq. (3.19) leads to unitary theories [21,[23][24][25], and knowing that the presence of a finite number of complex conjugate poles does not break the unitary condition either, we would expect also in our case the optical theorem to hold at all orders in perturbation theory. However, further investigations on perturbative unitarity with infinite complex conjugate poles will be subject of future works.

Ultraviolet behavior
The tree-level Lagrangians in local Galilean theories are free from any kind of instabilities as the corresponding field equations turn out to be of second order in the derivatives [61]. However, as shown in Ref. [77], quantum corrections can introduce new higher derivative terms which cause classical instabilities and unitarity violation, since ghost modes are introduced.
Very interestingly, nonlocal generalizations of Galilean theories can avoid any ghost degrees of freedom, not only at the tree-level but also when quantum loop corrections are taken into account, due to their improved UV behavior. As an example, we can consider the following nonlocal model: where λ is a dimensionless coupling constant; in the low energy regime, M 2 s , the action in Eq.(3.22) reduces to some version of local Galilean actions including, for instance, the term in Eq.(3.6).
The corresponding Feynman rules in Euclidean space are given by the nonlocal propagator Π(k) in Eq.(3.9) and by the following three-vertex: The presence of the exponentials in both propagator and interaction vertex ameliorates the UV behavior of the theory not only at the tree-level but also at higher loop orders. From power counting arguments we can easily understand how loop integrals will behave: since both propagator and vertices are exponentially suppressed at high energy, the UV behavior of the loop integrals is generally governed by e (−I−3V )k 2 /M 2 s , where I is the number of internal propagators and V counts the number of vertices. However, there can be cases in which one of the three momenta in the vertex is external and two are integrated as part of the loop, so in this case we get the upper bound for the superficial degree of divergence which is given by Note that D is always negative, which is a good hint in favor of the finiteness of loop integrals for the model in Eq. (3.22). We can explicitly show this feature, and make it more clear, by computing the self-energy Σ(p) at 1-loop for the action in Eq.(3.22): (3.25) The integral for the self-energy is finite, indeed by taking the UV regime (k 2 M 2 s ) of the integrand we obtain where the factor 6 in the exponent is in agreement with the power counting in Eq.(3.24): D = −2 − 2(2) = −6. Note that the integrals in Eq.(3.25) can be computed in Euclidean space and then analytically continued back to Minkowski space by Wick rotating. From Eq.(3.26) we can notice that the loop integral is exponentially suppressed at high energy and does not need to be renormalized. Therefore, for the model in Eq.(3.22) we would expect no ghost degrees of freedom to emerge when quantum loop corrections are taken into account, unlike the local case studied in Ref. [77]. However, more general studies regarding the UV behavior of these nonlocal models are needed and will be subjected to future investigations.

Comparison with infinite derivative field theory (IDT)
In the context of IDT, the simplest action we can consider is [11,14] which, by making the field redefinitionφ = e − /2M 2 s φ, can be equivalently written as which has only one massless real pole at p 2 = 0, and has the same form of the first term of the nonlocal propagator in Eq.(3.9). Therefore, the main difference between the propagators in Eq.

Comparison with p-adic string
The action of p-adic string is given by [16] where φ(x) is a scalar field and represents the open string tachyon, the scale m s is the string mass and the parameter g s is the open string coupling constant. By defining M 2 s ≡ m 2 s /lnp we can also write the action in Eq.(3.31) as follows The first derivation of the p-adic action assumed p to be a prime number, but it was realized that it can be defined for any positive integer, and even the case p = 1 can be treated in some way [80]. At first sight, the nonlocal action in Eq.(3.7) seems a particular case of the p-adic action in Eq.(3.32), indeed for p = 1 we obtain the term φ 2 . However, the p-adic string for p = 1 has nothing to do with the action in Eq.(3.7): the 1-adic action also has a potential term as shown in Ref. [80], while the action in Eq.(3.7) is purely kinetic. Moreover, unlike the propagator in Eq.(3.9), the p-adic propagator is given by which possesses no poles. As for the high energy regime, both propagators in Eqs.(3.9,3.33) have a similar ameliorated UV behavior.

A nonlocal gravity with infinite complex conjugate poles
In this Section we will construct a linearized gravitational action around the Minkowski background by using the non-polynomial form-factor introduced in the previous sections, in such a way that the graviton propagator has a similar structure as the one in Eq.(3.9). We will also find the corresponding non-linear action up to quadratic order in the curvature invariants, around Minkowski background. Let us consider the following gravitational action, which is the most general quadratic in the curvature, torsion free and parity invariant [3, 4, 7-9, 22] 5 : (4.1) where F i ( ) are three form factors which we can be uniquely determined around Minkowski background by fixing the form of the graviton propagator [8,9]. Since we are interested in the linearized regime we can always neglect the Riemann squared term R µνρσ F 3 ( )R µνρσ ; indeed the following identity holds where O(R 3 ) takes into account of terms of the order O(h 3 ) and div includes total derivatives. By linearizing around the Minkowski, where κ := √ 8πG and h µν is the linearized metric perturbation, we obtain the following action up to O(h 2 µν ) [9]: where h ≡ η µν h µν is the trace, = η µν ∂ µ ∂ ν is the flat d'Alembertian and The gauge independent part of the saturated graviton propagator around the Minkowski background has the following general expression [3,8,9]: where we have suppressed the tensorial indices for simplicity; P 2 and P 0 s are the so called spin projection operators, which project along the spin-2 and spin-0 components, respectively; see Refs. [9,26,32] for more details. Note that for f = 1 = g we recover the graviton propagator, Π GR (p) = P 2 /p 2 − P 0 s /2p 2 , of the Einstein general relativity.

Gravitational action and propagator
The gravitational analog of the nonlocal scalar propagator defined in Eq.(3.9) can be found by making the following choice for the functions f (p) and g(p) : which fixes the form of the gauge independent part of the saturated propagator as follows Note that the linearized action in Eq.(4.4) with the nonlocal operator in Eq.(4.7) is invariant under the following enlarged Galilean transformation 6 : where ψ k has been introduced in Eqs. (2.19,2.21). Also in this case, the propagator possesses a massless real pole, p 2 = 0, and an infinite number of massive complex conjugate poles. However, as we have shown in Subsection 3.2 tree-level unitarity is maintained. From the choice made in Eq.(4.7), we can also find the expressions for the form factors in the full non-linear gravitational action, indeed we can easily check that Eq.(4.7) implies Thus, the nonlocal quadratic gravitational action is given by the following expression: where G µν = R µν −1/2g µν R is the Einstein tensor. Note that the two quadratic terms with 1/ and 1/ 2 are individually non-analytic, but their particular combination in Eq.(4.10), and so the gravitational action in Eq. (4.11), are analytic functions of . Such an analiticity is also respected at the linearized level, where the nonlocal operator in the kinetic term is given by e − /M 2 s − 1 / ; see Eq.(4.12) of the next subsection. In relation with the discussion in Subsection 3.2, the nonlocal gravitational model in Eq.(4.11) can be seen as an infinite derivative generalization of (local) Lee-Wick theories of gravity [75,76,82].

Linearized metric solution
We now wish to compute the gravitational potential generated by a point-like source in the nonlocal theory of gravity introduced above. By varying the linearized action in Eq.(4.4), and imposing the choice in Eq.(4.7), we can obtain the following linearized field equations: 12) where T µν is the stress-energy tensor of the matter sector.
By working in the transverse gauge, we can write a static and spherically symmetric linearized metric as follows [83] ds 2 = −(1 + 2Φ)dt 2 + (1 − 2Ψ)(dr 2 + r 2 dΩ 2 ), (4.13) so that h 00 = −2Φ < 1, h ij = −2Ψδ ij < 1, h = 2(Φ − 3Ψ), with Φ and Ψ being the two unknown metric potentials. Then, by considering a static point-like source, of mass m, described by the stress-energy tensor T µν = mδ 0 µ δ 0 ν δ (3) ( r), one can show that the metric potentials solve the following modified Poisson equation: We have not been able to compute the integral in Eq.(4.15) analytically, however we can do it numerically and its behavior is shown in comparison with the Newtonian potential in Fig. 2.
We can notice that the metric potential is nonsingular at the origin, unlike the Newtonian one, indeed nonlocality regularizes its behavior. We can exactly compute the value of the metric potential at r = 0, and is given by: n −x is the Riemann zeta function. Note that the linearity of the metric potential holds as long as the following inequality is satisfied:

Comparison with infinite derivative gravity (IDG)
Now we briefly make a comparison between the nonlocal gravity model we have introduced above and the infinite derivative gravity (IDG) studied in Refs. [7][8][9]. The simplest form of the gravitational action in IDG is given by the following choice for the form factors [8,9,27]: which fixes the form of the gravitational action as follows The gauge independent part of the graviton propagator around the Minkowski background takes the following form [3,[6][7][8]: As in the scalar case analyzed in Section 2, we note that there is a net difference between the two propagators in Eqs. (4.8, 4.20). The main difference is contained in the function f (p). Indeed, in IDG the function f IDG (p) does not introduce any new degrees of freedom besides the massless spin-2 graviton. While, for the gravitational model proposed in this paper, the function f (p) possesses zeros which introduce extra poles in the graviton propagator in Eq.(4.8), and they are complex and infinite in number. Nevertheless, both satisfy the tree-level unitarity.

Conclusions
In this paper we have introduced a new nonlocal model motivated by a symmetry principle. We have noticed that by working with nonlocal actions it is possible to enlarge the Galilean shift symmetry. Indeed, we have found a new class of transformations, φ −→ φ + ψ k (see Eq.(2.23)), under which our nonlocal actions are invariant, and such a wider class includes the Galilean shift symmetry as a subclass, for k = 0. For a scalar field we have constructed some possible actions which respect such a symmetry, and studied the corresponding field equation and propagator. The pole structure of the propagator turns out to be more complicated; indeed, besides the massless pole of the local theory, we have also an infinite number of complex conjugate poles which, in principle, may spoil unitarity. However, we have noticed that all the infinite complex poles do not contribute to the imaginary part of the propagator, so that unitarity turns out to be preserved at the tree-level. We have good hints that the optical theorem holds at all orders in perturbation theory, but further investigations will be subject of future works. Moreover, we have made some comparisons with other two already well known nonlocal theories, i.e. infinite derivative field theory and p-adic string. Our new nonlocal model can be seen as an infinite derivative generalization of Lee-Wick theories, as we have propagators with pairs of conjugate poles which are not finite in number, but infinite.
We have analyzed these particular form factors also in the context of linearized gravity. We have constructed its nonlocal action quadratic in the curvature, which in the linear regime around Minkowski manifestly exhibits the enlarged Galilean symmetry. Also in this case we have computed the propagator which, besides the massless graviton pole, possesses an infinite number of complex conjugate poles, but no real ghost-modes. We have also computed the gravitational potential generated by a point-like source, which turns out to be nonsingular at the origin, unlike the Newtonian one.
In the framework of nonlocal theories, the requirement of using exponentials of entire functions seemed to be fundamental in order to avoid ghost-modes in higher derivative gravity 7 [7,9] and, thus, preserving unitarity at the perturbative level [21,[23][24][25]. Very interestingly, we have been able to enlarge the class of nonlocal form factors which make the propagator ghost-free, without restricting ourself to functions f (p) which are equal to exponential of entire functions. However, further studies to show that the optical theorem holds at each order in perturbation theory will be subject of future works.
Analogously to the scalar field case, this new model of nonlocal gravity can be seen as an infinite derivative generalization of Lee-Wick theories of gravity. However, we have to be careful in the extension to full non-linear gravity because Galilean shift symmetry must be broken when we couple it to gravity even in finite order derivative theories. We need to consider further extension [85,86], which will be investigated in future publication. These new form factors might be useful to be explored also in the context of string field theory.
Such a nonlocal extension of the Galilean symmetry might be useful to formulate nonlocal effective field theories, analogously to the local case with the Galilean symmetry, but now by demanding the enlarged symmetry in Eq.(2.23) to be satisfied. This might be crucial for discriminating local and nonlocal theories. We leave such a possibility for future investigations.
Let us consider a 1-dimensional case and find the solution for the following differential equation e −∂x/Ms − 1 χ(x) = 0 ⇐⇒ e −∂x/Ms χ(x) = χ(x), (A.2) which also means to find the eigenfunctions χ of the operator e −∂x/Ms of unit eigenvalue. First of all, note that the following property holds e −∂x/Ms = e −∂x/Ms+i2πk ≡ e θ 1 , where k is an integer number and we have defined the differential operator Thus, we need to find solutions for the equation θ 1 χ(x) = 0 : where c is an integration constant, which can be uniquely fixed by imposing that for k = 0 we recover the shift symmetry in Eq.(A.1). Hence, we have explicitly shown that by working with the non-polynomial operator introduced in Eq.(A.2), we can enlarge the shift symmetry, which now becomes a subclass of a wider class described by the following family of parameters: {c, k} , (A. 6) which is expressed in terms of the following field transformation: Thus, the shift transformations correspond to the subfamily {c, 0} . The generalization to a 4-dimensional case is given by replacing the operator ∂ x by d µ ∂ µ with d µ being a constant vector: such a non-polynomial operator might break the Lorentz invariance explicitly.
where k is an integer number, 1 is the 4 × 4 identity matrix, and we have defined the differential matrix operatorθ Thus, we need to find solutions for the equationθ 1 λ(x) = 0 : where λ 0 is a constant Dirac spinor, which can be uniquely fixed by imposing that for k = 0 we recover the standard shift symmetry, while l µ is a vector satisfying −l µ γ µ /M s +2πik1 = 0. Hence, we have explicitly shown that by working with the non-polynomial operator introduced in Eq.(A.8), we can enlarge the shift symmetry, which now becomes a subclass of a wider class described by the following family of parameters: which is a natural nonlocal extension of the Dirac action. Such an action is invariant under the global U (1). Furthermore, if we introduce a gauge field A µ and covariantize the derivatives, the following operator e iγ µ Dµ , with D µ = ∂ µ + ieA µ , one can show that also the local gauge symmetry is preserved. Indeed, if ψ −→ e iα(x) ψ and A µ −→ A µ + 1/e∂ µ α(x), then (iγ µ D µ ) n ψ −→ e iα(x) (iγ µ D µ ) n ψ. The fermionic nonlocal action in Eq.(A.14) will be further investigated in future publication.