Infinite derivative gravity is a theory of gravity which attempts to remove cosmological and black hole singularities by adding extra terms to the Einstein–Hilbert action, which weaken gravity at short distances.

History

In 1987, Krasnikov considered an infinite set of higher derivative terms acting on the curvature terms and showed that by choosing the coefficients wisely, the propagator would be ghost-free and exponentially suppressed in the ultraviolet regime.^[1] Tomboulis (1997) later extended this work.^[2] By looking at an equivalent scalar-tensor theory, Biswas, Mazumdar and Siegel (2005) looked at bouncing FRW solutions.^[3] In 2011, Biswas, Gerwick, Koivisto and Mazumdar demonstrated that the most general infinite derivative action in 4 dimensions, around constant curvature backgrounds, parity invariant and torsion free, can be expressed by:^[4]

${\displaystyle S=\int \mathrm {d} ^{4}x{\sqrt {-g))\left(M_{P}^{2}R+RF_{1}(\Box )R+R^{\mu \nu }F_{2}(\Box )R_{\mu \nu }+C^{\mu \nu \lambda \sigma }F_{3}(\Box )C_{\mu \nu \lambda \sigma }\right)}$

where the ${\displaystyle F_{i}(\Box )=\sum _{n=0}^{\infty }f_{i_{n))\left(\Box /M^{2}\right)^{n))$ are functions of the D'Alembert operator ${\displaystyle \Box =g^{\mu \nu }\nabla _{\mu }\nabla _{\nu ))$ and a mass scale ${\displaystyle M}$, ${\displaystyle R}$ is the Ricci scalar, ${\displaystyle R_{\mu \nu ))$ is the Ricci tensor and ${\displaystyle C_{\mu \nu \lambda \sigma ))$ is the Weyl tensor.^[5] In order to avoid ghosts, the propagator (which is a combination of the ${\displaystyle F_{i}(\Box )}$s) must be the exponential of an entire function. A lower bound was obtained on the mass scale of IDG using experimental data on the strength of gravity at short distances,^[6] as well as by using data on inflation^[7] and on the bending of light around the Sun.^[8] The GHY boundary terms were found using the ADM 3+1 spacetime decomposition.^[9] One can show that the entropy for this theory is finite in various contexts.^[10][11]

The effect of IDG on black holes and the propagator was examined by Modesto.^[12][13][14] Modesto further looked at the renormalisability of the theory,^[15][16] as well as showing that it could generate "super-accelerated" bouncing solutions instead of a big bang singularity.^[17] Calcagni and Nardelli investigated the effect of IDG on the diffusion equation.^[18] IDG modifies the way gravitational waves are produced and how they propagate through space. The amount of power radiated away through gravitational waves by binary systems is reduced, although this effect is far smaller than the current observational precision.^[19] This theory is shown to be stable and propagates finite number of degrees of freedom.^[20]

Avoidance of singularities

This action can produce a bouncing cosmology, by taking a flat FRW metric with a scale factor ${\displaystyle a(t)=\cosh(\sigma t)}$ or ${\displaystyle a(t)=e^{\lambda t^{2))}$, thus avoiding the cosmological singularity problem.^{[3][21][22][23]} The propagator around a flat space background was obtained in 2013.^[24]

This action avoids a curvature singularity for a small perturbation to a flat background near the origin, while recovering the ${\displaystyle 1/r}$ fall of the GR potential at large distances. This is done using the linearised equations of motion which is a valid approximation because if the perturbation is small enough and the mass scale ${\displaystyle M}$ is large enough, then the perturbation will always be small enough that quadratic terms can be neglected.^[4] It also avoids the Hawking–Penrose singularity in this context.^[25][26]

Stability of black hole singularities

It was shown that in non-local gravity, Schwarzschild singularities are stable to small perturbations.^[27] Further stability analysis of black holes was carried out by Myung and Park.^[28]

Equations of motion

The equations of motion for this action are^[5]

${\displaystyle {\begin{aligned}T^{\alpha \beta }&=P^{\alpha \beta }\\&=G^{\alpha \beta }+4G^{\alpha \beta }F_{1}(\Box )R+g^{\alpha \beta }RF_{1}(\Box )R-4\left(\nabla ^{\alpha }\nabla ^{\beta }-g^{\alpha \beta }\Box \right)F_{1}(\Box )R-2\Omega _{1}^{\alpha \beta }+g^{\alpha \beta }\left(\Omega _{1\sigma }^{\sigma }\right)\\&\qquad +4{R^{\beta ))_{\mu }R^{\mu \alpha }-g^{\alpha \beta }R^{\mu \nu }F_{2}(\Box )R_{\mu \nu }-4\left(F_{2}(\Box )R^{\mu (\beta }\right)_{;\mu }^{;\alpha )}+2\Box \left(F_{2}(\Box )R^{\alpha \beta }\right)+2g^{\alpha \beta }\left(F_{2}(\Box )R^{\mu \nu }\right)_{;\mu ;\nu }-2\Omega _{2}^{\alpha \beta }+g^{\alpha \beta }\left(\Omega _{2\sigma }^{\sigma }+{\bar {\Omega ))_{2}\right)-4\Delta _{2}^{\alpha \beta }\\&\qquad -g^{\alpha \beta }C^{\mu \nu \lambda \sigma }F_{3}(\Box )C_{\mu \nu \lambda \sigma }+4{C^{\alpha ))_{\rho \theta \psi }F_{3}(\Box )C^{\beta \rho \theta \psi }-4\left[2\nabla _{\mu }\nabla _{\nu }+R_{\mu \nu }\right]F_{3}(\Box )C^{\beta \mu \nu \alpha }-2\Omega _{3}^{\alpha \beta }+g^{\alpha \beta }\left(\Omega _{3\gamma }^{\gamma }+{\bar {\Omega ))_{3}\right)-8\Delta _{3}^{\alpha \beta }\end{aligned))}$

where

${\displaystyle {\begin{aligned}\Omega _{1}^{\alpha \beta }&=\sum _{n=1}^{\infty }f_{1_{n))\sum _{m=0}^{n-1}\nabla ^{\alpha }\Box ^{m}R\nabla ^{\beta }\Box ^{n-m-1}R,\\{\bar {\Omega ))_{1}&=\sum _{n=1}^{\infty }f_{1_{n))\sum _{m=0}^{n-1}\Box ^{m}R\Box ^{n-m}R,\\\Omega _{2}^{\alpha \beta }&=\sum _{n=1}^{\infty }f_{1_{n))\sum _{m=0}^{n-1}\nabla ^{\alpha }\Box ^{m}{R^{\mu ))_{\nu }\nabla ^{\beta }\Box ^{n-m-1}{R^{\nu ))_{\mu },\\{\bar {\Omega ))_{2}&=\sum _{n=1}^{\infty }f_{1_{n))\sum _{m=0}^{n-1}\Box ^{m}{R^{\mu ))_{\nu }\Box ^{n-m}{R^{\nu ))_{\mu },\\\Delta _{2}^{\alpha \beta }&={\frac {1}{2))\sum _{n=1}^{\infty }f_{2_{n))\sum _{\ell =0}^{n-1}\nabla _{\nu }\left[\Box ^{\ell }{R^{\nu ))_{\sigma }\nabla ^{(\alpha }\Box ^{n-\ell -1}R^{\beta )\sigma }-\Box ^{\ell }\nabla ^{(\alpha }{R^{n}u}_{\sigma }\Box ^{n-\ell -1}R^{\beta )\sigma }\right],\\\Omega _{3}^{\alpha \beta }&=\sum _{n=1}^{\infty }f_{3_{n))\sum _{\ell =0}^{n-1}\nabla ^{\alpha }\Box ^{\ell }{C^{\mu ))_{\nu \lambda \sigma }\nabla ^{\beta }\Box ^{n-\ell -1}{C_{\mu ))^{\nu \lambda \sigma },\\{\bar {\Omega ))_{3}&=\sum _{n=1}^{\infty }f_{3_{n))\sum _{\ell =0}^{n-1}\Box ^{\ell }{C^{\mu ))_{\nu \lambda \sigma }\Box ^{n-\ell }{C_{\mu ))^{\nu \lambda \sigma },\\\Delta _{3}^{\alpha \beta }&={\frac {1}{2))\sum _{n=1}^{\infty }f_{3_{n))\sum _{\ell =0}^{n-1}\nabla _{\nu }\left[\Box ^{\ell }{C^{\lambda \nu ))_{\sigma \mu }\Box ^{n-\ell -1}{C_{\lambda ))^{(\beta |\sigma \mu |;\alpha )}-\Box ^{\ell }\nabla ^{(\alpha }C_{\sigma \mu }^{\lambda \nu }{C_{\lambda ))^{\beta )\sigma \mu }\right].\end{aligned))}$