The gravitational field of static branes in linearized ghostfree gravity
Abstract
We study the gravitational field of static branes in dimensional Minkowski space in the framework of linearized ghostfree (GF) gravity. The concrete models of GF gravity we consider are parametrized by the nonlocal form factors and , where is the scale of nonlocality. We show that the singular behavior of the gravitational field of branes in General Relativity is cured by shortrange modifications introduced by the nonlocalities, and we derive exact expressions of the regularized gravitational fields, whose geometry can be written as a warped metric. For large distances compared to the scale of nonlocality, , our solutions approach those found in linearized General Relativity.
pacs:
04.20.Dw, 04.20.Jb, 04.50.Kd, 11.27.+d AlbertaThy0218I Introduction
Einstein’s theory of General Relativity (GR) describes gravitational physics on the scale of our solar system remarkably well, and it has a welldefined Newtonian limit in the context of weak gravitational fields. However, one of the predictions of GR is the existence of spacetime singularities: spacetime loses geodesic completeness, and the curvature diverges. These problematic features in the ultraviolet (UV) regime are deemed physically unreasonable, and hence one may ask the question: what is the correct UV completion of gravity?
A fairly generic approach is to UV complete GR by adding terms that are quadratic or of higher order in the curvature, or contain more derivatives. This is a natural procedure since these terms are generated by the effective action during quantization anyway. For an example, see e.g. Stelle (1977, 1978). As it turns out, the gravitational potential of a point mass becomes regular for this class of theories (see Asorey et al. (1997); Modesto et al. (2015) for details). Unfortunately, the propagators of these theories contain ghost modes, reflecting an inherent instability of this class of theories Stelle (1977, 1978); Biswas et al. (2006); Barnaby and Kamran (2008). Hence, the addition of other higher derivative terms may mitigate the UV behavior but is usually accompanied by extra unphysical ghost or tachyon modes.
An interesting approach is to consider a theory with infinitely many derivatives, which is equivalent to a nonlocal modification of GR. Nonlocal field theories were considered a long time ago (see, e.g., Efimov (1967, 1968, 1972, 1974); Efimov et al. (1977)); see also Hehl and Mashhoon (2009); Blome et al. (2010) for a more recent approach. Also, they appear naturally in the context of noncommutative geometry deformation of GR Nicolini et al. (2006); Spallucci et al. (2006) (see the review Nicolini (2005) and references therein). The initial value problem in nonlocal theories was studied in Barnaby and Kamran (2008); Barnaby (2011).
In some cases, with a proper choice of nonlocal form factors, UV singularities may be avoided while no extra propagating degrees of freedom appear. Such modifications of GR are called ghostfree (GF) gravity (see, e.g., Tomboulis (1997); Modesto (2012a); Biswas et al. (2012); Modesto (2012b); Biswas et al. (2013, 2014); Modesto and Rachwal (2014); Tomboulis (2015a, b); Modesto and Rachwał (2017) and references therein). These theories are usually characterized by a mass parameter or a length scale at which nonlocalities become important. At large scales GF gravity is expected to reproduce GR. The theory may have another energy scale, , where quantum fluctuations of the metric become large. We assume that (or, equivalently, ). In other words, in the approach adopted in this paper the metric is treated as classical, that is, its quantum fluctuations are considered to be small. Therefore the length parameter is chosen to be larger than the Planck length, the string scale, or some other scale, depending on the fundamental quantum theory of gravity.
Recently, the entropy of black holes Conroy et al. (2015a), quantum effects like oneloop renormalization Modesto et al. (2017a), and the Unruh effect Kajuri (2017); Modesto et al. (2017b) have been studied in GF gravity in higher dimensions. The study of GF gravity as applied to the problem of cosmological singularities can be found in Biswas et al. (2010); Calcagni et al. (2014). GF gravity applications to the problem of black hole singularities were considered in Hossenfelder et al. (2010); Modesto et al. (2011); Zhang et al. (2015); Conroy et al. (2015b); Li et al. (2015); Calcagni and Modesto (2017); Cornell et al. (2017). It is well known that linearized GF gravity regularizes the gravitational field of pointlike sources Biswas et al. (2012); Modesto et al. (2011); Frolov et al. (2015), and recent studies indicate that this may remain true even in the full nonlinear GF theory Koshelev et al. (2018). Moreover, all fourdimensional curvature invariants are finite everywhere, and at the location of the point particle spacetime approaches conformal flatness Buoninfante et al. (2018). One can expect that for black holes of large mass and size the GF modification of gravity results only in small corrections. In particular, the properties of their horizons remain practically the same. For black holes of small masses (), GF gravity effects may become very important. For example, dynamical solutions, like collapsing matter or headon collision of particles in four and higher dimensions, were studied in Frolov (2015); Frolov et al. (2015), showing that there exists a mass gap for black hole formation in GF gravity.
One can expect that these properties (regularity of solutions and mass gap effect) might be valid in the complete GF gravity even in the strong field regime. Before studying this rather complicated problem it is instructive to demonstrate that linearized GR gravity regularizes the gravitational fields of not only pointlike particles, but also the field of other infinitely thin sources, e.g., cosmic strings, membranes, and branes in arbitrary dimensions. The aim of this paper is to demonstrate this explicitly and to obtain explicit solutions of the linearized GF gravity equations for such objects.
Let us first discuss the wellknown example of a cosmic string Vilenkin (1981) in the context of GR. Expressed in Cartesian coordinates, a straight cosmic string in Minkowski space located along the axis has the stressenergy tensor
(1) 
Defining , a teststring solution in the framework of linearized GR in the harmonic gauge is given by
(2) 
where and is an integration constant.
To linear order the curvature tensor vanishes everywhere except for the origin of the plane along the axis, where it has a like singularity producing an angle deficit. Accordingly, the geometry of a straight cosmic string corresponds to Minkowski space with a conical deficit angle .
Is this curvature singularity along the axis cured in GF gravity? In order to answer this and similar questions in quite some generality, in this work we study the static gravitational field of branes in dimensional Minkowski space, employing the framework of linearized GF gravity. In particular, we show that in the case of and (cosmic strings) the conical singularity is resolved. Moreover, we derive explicitly the regular potentials for any brane in any number of spacetime dimensions.
This paper is organized as follows: In Sec. II we briefly sketch the derivation of linearized GF gravity, before discussing the general gravitational field of branes in Sec. III. In Sec. IV we list a set of explicit examples and discuss their behavior in detail, before discussing the obtained results in Sec. V.
Ii Linearized ghostfree gravity
Consider a weak perturbation on Minkowski space in Cartesian coordinates,
(3) 
where is the Minkowski metric, and is assumed to be small. The dynamics of can be derived by variation of the GF gravity action. In order to study the action of linearized GF gravity, it is sufficient to consider only terms linear and quadratic in curvature, with generic nonlocal form factors. Using the symmetry properties of the Riemann tensor, the Bianchi identities, and the commutativity of the covariant derivations up to linear in the curvature terms one can show Barvinsky and Vilkovisky (1990); Modesto et al. (2015) that there are only two independent nonlocal form factors that characterize nonlocal linearized gravity in arbitrary dimensions. In dimensions, a generic action for linearized GF gravity written in Cartesian coordinates then takes the form Conroy et al. (2015a); Frolov and Zelnikov (2016a)
(4)  
where and are arbitrary, dimensionless form factors and is the gravitational constant in dimensions. This general action can be employed to describe various linearized gravitational theories Biswas et al. (2013):

GR is recovered for ,

L(R) gravity for the choice , ,

Weyl gravity for , , , where is a parameter of dimension mass.
At any rate, in order to recover GR in the infrared (IR) domain the form factors functions must satisfy the condition .
Let now be a stressenergy tensor of matter,
(5) 
then the linearized field equations for become
(6) 
The resulting nonlocal theory of gravity is ghostfree if it is described by form factors and that are entire functions of the operator . Such entire functions can be written as the exponential of any finite polynomial. According to the classification adopted in Frolov and Zelnikov (2016a), we call the simplest example of a GF theory with the choice
(7) 
This version of GF gravity has been extensively studied in the literature Tomboulis (1997); Biswas et al. (2012); Modesto (2012a, b); Biswas et al. (2014, 2013); Modesto and Rachwal (2014); Tomboulis (2015a, b); Modesto and Rachwał (2017), and this choice corresponds to the absence of propagating spin zero gravitational modes. The mass parameter describes the length scale below which nonlocalities start become relevant and GF gravity solutions deviate significantly from those obtained in GR. Recently Frolov and Zelnikov (2016a), more general theories of the type with
(8) 
have been studied. It was demonstrated that in the presence of timedependent sources gravities for odd suffer from instabilities, whereas they are stable in the case of even .
Iii Pbranes in Ddimensional Minkowski space
iii.1 brane ansatz
Let us consider a plane dimensional timelike surface (“brane”) embedded in dimensional Minkowski space. We use Cartesian coordinates such that are the coordinates on the brane (), while are the spatial coordinates in the bulk directions (with ). Thus we have and the brane is located at . From now on, we call the number the codimensionality of the brane.
The stressenergy tensor of the brane is
(9) 
where is the (positive) surface tension. In this approximation for the stressenergy tensor we assume that the thickness of the matter source is much smaller than the characteristic length parameter .
The presence of stressenergy (9) will lead to a deviation from Minkowski spacetime which we shall refer to as . As an ansatz for the perturbed spacetime we choose a warped geometry Dabholkar et al. (1990); Stelle (1998)
(10) 
We shall assume that the above metric in Cartesian coordinates deviates only slightly from Minkowski space. As such, there exists a gauge where has the form
(11) 
Here, the functions and depend on the distance from the brane defined as . We consider these functions as perturbations, such that and .
This geometry has the following isometries:

Poincaré symmetry in the sector,

rotation symmetry in the sector.
The full isometry group of (11) is hence . The discrete symmetries further guarantee that the surface described by (the brane) is geodesic and hence minimal.
Let us remark that the Newtonian limit can be read off from the component, and we define the Newtonian potential to be given by
(12) 
for later convenience. Incidentally, for the cosmic string in GR, Eq. (2), one has , implying that test particles do not feel any forces acting on them in the Newtonian limit.
iii.2 GF linearized equations and their solution
For the warped geometry (11) the linearized equations (6) take the following form:
(13)  
(14)  
where we defined such that is the dimensional Laplacian, and is the Euclidean metric in the sector. In a simple case when
(15) 
the homogeneous equation (14) can be solved by setting
(16) 
The remaining equation becomes
(17) 
where we introduced the constant prefactor
(18) 
Equation (17) can be solved by using the method of Green functions Frolov and Zelnikov (2016a, b). Given the function , the Green function for is
(19) 
where denotes the Bessel function of the first kind. The Green functions for can be determined via
(20) 
The above equality follows directly from the differential properties of Frolov and Zelnikov (2016a). The exact solution for a brane in dimensional spacetime is then given by
(21) 
Let us notice that the radial dependence of is universal and does not depend on the dimensionality of the brane: it only depends on the brane codimensionality . The parameter enters only in a constant prefactor. In this linear approximation to the full GF gravity, one can also consider a superposition of like branes: the obtained Green functions can be used to generate linearized GF gravity solutions corresponding to thick branes. It is clear that such a “smeared” solution remains regular and has the correct GR behavior at far distances. In a special case when the thickness of the brane is of the order of Planck length, the gravitational field at the brane position would be regular, but slightly different from the GF solution. However, at the scale it will coincide with the latter.
Iv and theory
Let us consider the potential (21) in the cases of and theory. That is, from now on we put as required by . In this case, the relation (16) simplifies and takes the form
(22) 
One of the immediate consequences of this relation is that for codimension 2, that is, when , the Newtonian potential vanishes. This is a generalization of the result for the cosmic string.
In what follows, we will first derive the explicit Green functions for and theory, and then discuss the relation to the Green functions of linearized GR.
iv.1 Green functions of and theory
In order to distinguish the Green functions for different versions of GF gravity, we use the corresponding index as a superscript and write the Green function in the form . The Green function for in theory is Frolov and Zelnikov (2016a)
(23) 
where denotes the lower incomplete gamma function. For we find the following expressions:
(24) 
where Ein(x) is the complementary exponential integral and erf(x) is the error function:
(25) 
In the above, denotes the exponential integral, and is the Euler–Mascheroni constant Olver et al. (2010). The expressions for have been calculated by using Eq. (20), and the above set (24) is sufficient to calculate for any .
In theory, the general Green function for takes the more complicated form Frolov and Zelnikov (2016a)
(26)  
where we defined , denotes the gamma function, and corresponds to the generalized hypergeometric function Olver et al. (2010). Again, the cases can be obtained by making use of Eq. (20). We find the following expressions:
(27)  
where again .
Lastly, the Green functions of linearized GR can be obtained from both the Green functions , or the Green functions , by considering the limit . Denoting the GR Green functions as , we find the expressions
(28) 
For , the logarithmic dependence on leads to a divergence; however, it is of no physical relevance for matter sources with compact support.
Note that both the Green functions for theory (24) and theory (IV.1) are manifestly regular at . In the cases their normalization has been chosen such that for ; in the cases their Green functions are unbounded for large radii, which is why we chose instead for . In contrast, the GR Green functions are singular at .
See the plots of these Green functions in Fig. 1.
iv.2 Linearized curvature
It is also useful to study the curvature of the brane solutions. In linear approximation, the Riemann tensor, the Ricci tensor, and the Ricci scalar are Conroy (2017)
(29)  
(30)  
(31) 
Besides the curvature tensors and invariants, it is interesting to consider the quantity
(32) 
where is the linearized Einstein tensor and is the timelike Killing vector. We see that in GR, corresponds to the energy density perceived by a static observer tangential to .
Substituting (11) and the interrelation of and , see Eq. (16), into Eqs. (29)–(31), we find for the Ricci scalar as well as for the energy density
(33) 
The above justifies the interpretation of the Ricci scalar as a rescaled energy density. Moreover, the quadratic curvature invariants take the form
(34)  
(35) 
where denotes the tracefree Ricci tensor (65). In , one may consider the Chern–Pontryagin pseudoscalar. It vanishes for warped geometries such as (11):
(36) 
For more details of the curvature components, see appendix A. Due to the radial dependence of the function , one has the following identities:
(37)  
(38) 
Using Eqs. (20) and (21), the above invariants can be rewritten in terms of higherdimensional Green functions. Moreover, by virtue of Eq. (21) and the above relations, for a qualitative study of the behavior of the curvature invariants it is sufficient to substitute into the expressions (37) and (38). Hence let us consider the dimensionless invariants
(39)  
(40) 
We may think of the invariant as a rescaled energy density as per Eq. (33). Using Eq. (20) as well as Eqs. (37) and (38), the above invariants can be recast into the form
(41)  
(42) 
where we suppressed the superscripts “” for sake of clarity. In Fig. 2 we visualize these invariants for some typical cases in and theory; for the other cases, see the plots in appendix B.
iv.3 Explicit examples
Let us now discuss various examples of branes in dimensional Minkowski space in more detail. To that end, one may consider branes for different values of , subject to . In dimensions, we then have the following special cases of branes:

point particle: (),

cosmic string: (),

domain wall: (),

angle deficit configurations: (, ).
Point particle in any
A point particle can be regarded as the somewhat degenerate case of a 0brane, and one can simply set such that . This case has been previously studied in Frolov and Zelnikov (2016a), and as follows from Eqs. (24) and (IV.1), their gravitational potential is regularized. The Newtonian potential takes the form
(43) 
where has now dimensions of a mass. The above correctly reproduces the wellknown trivial Newtonian limit in dimensions, as well as the modified Newtonian potential in, say, theory in :
(44) 
The limiting case reproduces the potential.
As was recently shown Buoninfante et al. (2018), the Weyl tensor vanishes at the location of pointlike particles in theory. Here, let us generalize this result in two directions by considering point particles in any dimension for a general theory. Assuming that the gravitational field of the point particle is regular Biswas et al. (2012); Modesto et al. (2011); Frolov et al. (2015), the metric function has the following expansion at small distances:
(45) 
where and are some constant factors. Then, to leading order, one has
(46) 
Upon this substitution it is easy to see that the Weyl tensor (appendix A) vanishes identically iff :
(47)  
(48)  
(49) 
Therefore, as long as point particles have a regular gravitational field, the Weyl tensor vanishes at their location in theories. No such theorem holds for branes with , except when , as we shall see below.
Cosmic string in
A cosmic string in , our example from the Introduction, is described by and such that
(50) 
The fact that vanishes corresponds to the absence of a gravitational field of the cosmic string in the Newtonian approximation and in linearized GR. In GF theories, however, one may calculate that the curvature is nonzero, even for . In the case of theory, the invariants take a simple form ():
(51)  
(52)  
(53) 
Let us notice that for all expressions for the curvature invariants can be expressed solely in terms of , and does not enter them.
Domain walls in any , conformal flatness
Domain walls are surfaces of dimension such that and . Therefore one has
(54) 
In this special case, the metric (11) is conformally flat. An easy way to see this is to introduce the new “radial” coordinate according to
(55) 
upon which the metric (11) becomes . Using the relations in appendix A, as well as the onedimensional identity , one can easily check that the Weyl tensor vanishes for , as it must.
Angle deficit configurations
We have seen that a cosmic string, in the GR limit , does not gravitate but creates an angle deficit in the twodimensional surface that it pierces. This triviality of the gravitational field is mirrored in the vanishing of the Newtonian potential for a cosmic string.
Lastly, let us generalize these solutions to the case for any such that and
(56) 
The shape of the warped geometry (11) then guarantees that one still has an angle deficit of in the plane, regardless of the dimensionality of the brane. However, as in the case of the cosmic string, the curvature invariants within GF theory do not vanish for this case.
V Conclusions
In this paper we obtained and studied the gravitational field of infinitely thin objects such as point particles, strings, and branes, in the context linearized GF gravity. The characteristic property of these objects is that their stressenergy tensor is described by a function distribution localized at a point, straight line, or flat planes, respectively. We considered two versions of linearized GF gravity, and , the form factors being
(57) 
In the above, the parameter is the energy scale at which the GF theories start to deviate significantly from GR. Equivalently, one may think of as the length scale of nonlocalities induced by the form factors and . Consequently, for one has , thereby recovering GR in the infrared.
In our calculations we assumed that the thickness of the branes is much smaller than , justifying the approximation of the stressenergy of the branes by simple distributions. In the linear approximation to the full GF gravity, the obtained Green functions can be used to generate solutions also in the case of thick branes.
The explicit brane solutions in linearized GF gravity in arbitrary number of spacetime dimensions possess two common properties: First, all of them are finite and regular at the position of the brane. Secondly, at far distances from the brane the GF gravity solutions asymptotically coincide with the solutions found in linearized GR (for the same, like stressenergy tensor). This second property is connected with the assumption that for , guaranteeing the proper infrared limit of the GF theory. The first property is a consequence of the adopted nonlocality of the theory. It has a rather simple interpretation. The equation
(58) 
for a static source can be rewritten into the equivalent
(59) 
Given a choice of nonlocal form factors, the effective current is smeared. This corresponds precisely to the quantity we defined in Sec. IV.2. As we have shown explicitly, this function is smooth for like sources. For pointlike particles this effect is well known, and has been generalized in this work to static branes in dimensional Minkowski space.
Let us emphasize that the solutions for the gravitational potential of the branes are quite similar for and models. This happens because the sources are static. In the timedependent case this is not true. In particular, this difference has been illustrated in case of the radiation of a timedependent like source in the connection with a GF massless scalar field Frolov and Zelnikov (2016b). In general, theory is sensitive to the difference of the spacelike and timelike directions as a result of the Lorentz signature of the metric. , on the other hand, has a better behavior. However, as argued above, for the static sources as considered in this work, this difference between and is irrelevant.
It might be interesting to study the gravitational field of branes in the complete, nonlinear GF gravity. One might hope that such a problem can be treated because of the high symmetry of the brane sources. It is also interesting to consider linearized GF gravity solutions on constant curvature backgrounds like (A)dS and to study how GF gravity might manifest itself in the AdS/CFT correspondence.
Acknowledgments
J.B. is grateful for a Vanier Canada Graduate Scholarship administered by the Natural Sciences and Engineering Research Council of Canada as well as for the Golden Bell Jar Graduate Scholarship in Physics by the University of Alberta. V.F. and A.Z. thank the Natural Sciences and Engineering Research Council of Canada and the Killam Trust for their financial support. V.F. also thanks the Yukawa Institute for Theoretical Physics at Kyoto University, where this work was completed during the workshop YITPT1702 “Gravity and Cosmology 2018.”
Appendix A Curvature expressions
Let us calculate the curvature tensors for the metric
(11) 
subject to the substitution
(22) 
We denote the flat background metric on the sector as , where now . The flat background metric on the sector is denoted by . Then, the Riemann tensor takes the form
(60)  
(61)  
(62) 
The Ricci tensor is
(63) 
and the Ricci scalar is simply
(64) 
We can construct the tracefree Ricci tensor:
(65)  
(66)  
(67) 
Let us define the Weyl tensor
(68)  
(69)  
For the components we obtain
(70)  
(71)  
(72) 
Using the above relations, the expressions in Sec. IV.2 can be readily derived.
Appendix B Curvature invariants visualized
Footnotes
 Because the form factors and are dimensionless, the d’Alembertian can only enter these in combination with at least one length scale