Quantum spreading of a self-gravitating wave-packet in singularity free gravity

In this paper we will study for the first time how the wave-packet of a self-gravitating meso-scopic system spreads in theories beyond Einstein's general relativity. In particular, we will consider a ghost-free infinite derivative gravity, which resolves the $1/r$ singularity in the potential - such that the gradient of the potential vanishes within the scale of non-locality. We will show that a quantum wave-packet spreads faster for a ghost-free and singularity-free gravity as compared to the Newtonian case, therefore providing us a unique scenario for testing classical and quantum properties of short-distance gravity in a laboratory in the near future.


I. INTRODUCTION
On large distances and late times the gravitational interaction is well described by the theory of general relativity (GR) that, indeed, has been very successful since Einstein's initial work, being tested to a very high precision in the infrared (IR) [1]. The most recent success of GR comes from the observation of gravitational waves from merging of binary blackholes which gave a further confirmation of its predictions [2]. Despite these great achievements, our knowledge of the gravitational interaction in the ultraviolet (UV) is still very limited: suffice to say that the inverse-square law of the Newtonian potential has been tested only up to 5.6×10 −5 m in torsionbalance experiments so far [3]. This means that any modification from the Newtonian 1/r-fall is expected to happen in the large range of values going from the lower bound 0.004 eV to the Planck scale M p ∼ 10 19 GeV. This is the place where nature should manifest a different behaviour compared to GR and where either quantum or classical modification from GR should appear in order to solve problems that still remain unsolved as, for example, blackhole and cosmological singularities that make Einstein's theory incomplete in the UV. There have been many theoretical attempts that try to modify GR in the UV regime but none of them have been sufficiently satisfactory so far. In fact, only the experiment will be able to tell us whether the gravitational interaction is really quantum or not, and, in both cases, whether the classical properties are also modified.
Recently, a new scenario has been proposed in which by studying the quantum spread of the solitonic wavepacket for a self-gravitating meso-scopic system one can test and constrain modified theories of gravity in the near future [4]. In this framework, the so called infinite derivative gravity (IDG) [5] was considered as example of alternative theory: it belongs to the class of nonlocal ghost-and singularity-free theories of gravity. It was shown that in the non-relativistic and in the weak-field regimes the dynamics of the matter-sector is governed by a Schrödinger equation with a non-linear self-interaction term [4]: where G = 1/M 2 p and M s represents the scale of new physics at which non-locality-effects should become relevant, i.e. 0.004eV ≤ M s ≤ 10 19 GeV.
Such an integro-differential equation can have two completely different physical interpretations, correspondingly the wave-function ψ( x, t) can assume two different meanings 1 : 1. It can appear when gravity is quantized and directly coupled to the stress-energy tensor operator. In this case, it is derived as a Hartree equation in a mean-field approximation; ψ( x, t) has the meaning of wave-function associated to an N -particle state, with large number of particles (N → ∞), i.e. a condensate [4].
2. Moreover, Eq. (1) can be seen as a fundamental equation describing the dynamics of a selfgravitating one-particle system, when considering a semi-classical approach where gravity is coupled to the expectation value of the stress-energy tensor; in this case ψ( x, t) represents a one-particle wavefunction. In such semi-classical framework gravity is treated as a classical interaction, while matter is quantized.
Note that, in case 1. the non-linearity emerges when considering the limit of large number of particles, i.e. in the mean-field regime; while in 2. one has non-linearity even for a one-particle state, bringing to a modified Schrödinger equation.
In this respect, a semi-classical approach to IDG would seem more speculative as, not only we would modify GR, but also quantum mechanics. While considering the case of quantized gravity, and studying the dynamics of a selfgravitating condensate, could be particularly more interesting as the main motivations of IDG concern problems emerging when one tries to quantize gravity. However, in this manuscript we wish to make a general treatment by considering both cases of semi-classical and quantized gravitational interaction, and discuss the experimental feasibility of the model in both cases.
It is also worth emphasizing that the non-linear potential term in Eq. (1) can be split in two parts as follow [4]: From the last decomposition one can notice that the first term contains all information about the non-local nature of the gravitational interaction, while the second one has the same form of the usual Newtonian self-potential that also appears in the well-known Schrödinger-Newton equation, see [7][8][9]. In Ref. [4] it was shown that Eq. (1) admits stationary solitonic-like solutions for the ground-state and it was found that in the case of IDG the energy E and the spread σ of the solitonic wave-packet turn out to be larger compared to the respective ones in Newtonian gravity, i.e. E IDG ≥ E N and σ IDG > σ N ; these are effects induced by the non-local nature of the gravitational interaction [4].
The expression of the ground-state energy is given by that in the limit M s σ > 2 recovers the energy of New- The above Eq. (3) shows the action of two kinds of forces that are completely different in nature: a quantum-mechanical kinetic contribution that tends to spread the wave-packet and the gravitational potential which takes into account the attractiveness of gravity coming from the non-linear term of Eq. (1). In a stationary scenario the two contributions balance each other and the soliton-like solution above can be found.
In this paper, unlike Ref. [4], we are more interested in studying non-stationary solutions of Eq. (1). In particular we want to understand how the spreading of the wavepacket is affected by the presence of a non-local gravitational self-interaction. The analysis that we will present will apply to both cases of semi-classical and quantized gravity as the main equation is mathematically the same.
As pointed out in Refs. [10][11][12] 2 , where numerical studies of the Schrödinger-Newton equation were made, there should exist a threshold mass µ, such that the collapse of the wave-function induced by gravity will take place for any m > µ. In Ref. [10] it was noticed that the collapsing behavior appears only if the initial state of the quantum system has negative energy, such that the attractive contribution of self-gravity dominates. From this last observation, we understand that a possible way to find an analytical estimation for the threshold mass µ is to equate kinetic and gravitational contributions in Eq. (3); thus we obtain (4) Note that in the limit when M s σ > 2, Eq. (4) gives the threshold mass similar to the case of Newtonian theory: Eqs. (4) and (5) clearly show that non-locality implies a larger value of the threshold mass, i.e. µ IDG > µ N for any values of M s and σ. If we choose the current lower bound on the scale of non-locality, M s = 0.004 eV [25], and σ = 500 nm, 3 the values of the threshold masses are µ IDG 3.5 × 10 −17 kg and µ N 6.7 × 10 −18 kg. From an analytical estimation one expects that for masses, m > µ IDG , the self-gravitating quantum wavepacket collapses, while for masses, m < µ IDG , one would expect no collapse of the wave-packet, but only a slowdown of the spreading compared to that of the freeparticle case.
We now wish to study the quantum spreading of a selfgravitating wave-packet and understand how it is affected by singularity-free gravity, without taking into account 2 In Refs. [10][11][12] the authors mainly focused on the semi-classical approach, where gravity is treated classically. However, the following treatment will also apply to the case of quantized gravity as one has to consider the same integro-differential equation (1). Let us keep in mind that in the semi-classical approach the quantum wave-packet represents a one-particle wave-function; while, when gravity is quantized, it is associated to the dynamics of a many-particle system, i.e. a condensate. 3 It is worthwhile to note that the special choice we have made for the value of the spread, σ = 500 nm, corresponds to the actual slit separation d in a Talbot-Lau interferometry setup [13]. The slit separation is related to both length L of the device and de Broglie wave-length λ = h/mv, where h is the Planck's constant, m the mass of the particle and v is its velocity, through the relation L = d 2 /λ. The interference with larger masses requires smaller wavelengths, which in turn means shorter slitseparations. Moreover, the most massive quantum systems which have been seen showing interference are organic molecules with a mass of the order of 10 −22 -10 −21 kg [14].
the collapsing phase. It means we will work in a regime in which we can assume that the non-linear contribution in Eq. (1) is smaller compared to the kinetic term: In this regime non-linearity-effects are sufficiently small and it allows us to find non-stationary solutions by applying the Fourier analysis to Eq. (1). From Eq. (6), we can also define the dimensionless parameter that quantifies the degree of non-linearity due to selfgravity. It depends on the initial data through the mass m and σ, that can represent the initial spread of the selfgravitating quantum wave-packet. When ξ < 1, we can assume that the non-linear effects are sufficiently small. In the Newtonian limit, M s σ 0 > 2, we will obtain ξ ∼ Gm 3 σ.
A comparison between modified theories of gravity, IDG in our case, and Newton's gravity, will provide us a new and unique framework to test short-distance gravity beyond Einstein's GR. This paper is organized as follows: first of all we will briefly introduce ghost-free and singularity-free IDG; then we will study the spreading solutions of Eq. (1) with the aim of comparing free, Newton and IDG cases; finally there will be a summary and a discussion on current and near future experimental scenarios in both cases of semi-classical and quantized gravity.

II. INFINITE DERIVATIVE GHOST-FREE AND SINGULARITY-FREE GRAVITY
There have been many attempts to modify GR by introducing higher order derivative contribution in the action, especially a conformal gravity containing quadratic terms in the curvature like R 2 , R µν R µν , R µνρσ R µνρσ . Such quadratic theory of gravity turns out to be conformal as well as renormalizable, but it suffers from the presence of a massive spin-2 ghost field that makes the theory classically unstable and non-unitary at the quantum level [15].
Recently, it has been noticed that by considering an infinite number of derivatives in the quadratic curvature gravitational action one can prevent the presence of ghost [5,16]. At the same time, such a ghost-free action also improves the behaviour of the gravitational interaction in the UV regime showing a non-singular potential and a vanishing gravitational force: Φ → const and F g → 0 as r → 0, where F g represents the mutual force between two particles separated by the distance r [5].
The most general torsion-free, parity-invariant and quadratic covariant action that contains an infinite number of derivatives has been constructed around constant curvature backgrounds, and reads [5,17]: where α is a dimensionful coupling, s ≡ /M 2 s and ≡ g µν ∇ µ ∇ ν , where µ, ν = 0, 1, 2, 3, and the mostly positive metric signature, (−, +, +, +), is chosen. The information about the presence of infinite derivatives is contained in the three gravitational form factors F i ( s ) which have to be analytic functions of , n , thus we can smoothly recover GR when we take the limit → 0. These form factors can be further constrained by requiring general covariance, that no additional dynamical degrees of freedom propagate other than the massless graviton and the ghost-free condition, that around Minkowski background is given by [5]. Note that around a constant curvature spacetime we can set F 3 = 0 [5], without loss of generality.
As we have already mentioned above, the parameter M s represents the scale of non-locality where gravity, described by this class of ghost-and singularity-free theories, shows a non-local nature [16,[18][19][20][21][22][23][24]. The current constraints on M s comes from torsion-balance experiments which have seen no departure from the Newtonian 1/r-fall up to a distance of 5.6×10 −5 meters, that implies M s ≥ 0.004 eV [25,26].
Furthermore, it is worth noting that IDG-theory can also resolve cosmological singularity, see [5,16,27], and the non-local nature of gravity can possibly even play a crucial role in the resolution of blackhole singularity as pointed out in Ref. [28]; while at the quantum level it is believed that the action in Eq. (8) describes a gravitational theory that is UV-finite beyond 1-loop [18,19,22,23].
The ghost-free condition of IDG demands a special choice for the gravitational form factors, see [5,16] Generally, a( s ) should be exponential of an entire function [5,16,19,22], in order to avoid additional dynamical degrees of freedom other than the massless spin-2 graviton, therefore no propagating ghost-like states. In fact, any generalised form of exponential of an entire function yields a similar behavior in the UV and IR regimes, namely a similar non-singular modified gravitational potential that for large distances recovers the Newtonian 1/r-fall [25,32]. By linearizing the action in Eq. (8) and going to momentum space one can easily show that the choice in Eq. (9) does not introduce any extra degrees of freedom in the gravity sector. Indeed, as shown in Ref. [5,33,34], the gauge-independent part of the propagator corresponding to the linearized action around Minkowski spacetime is given by where P 2 and P 0 are the well known spin projector operators that project any symmetric two-rank tensor along the spin-2 and spin-0 components, respectively; Π GR = P 2 /k 2 − P 0 s /2k 2 is the GR propagator. For the special choice a(−k 2 ) = e k 2 /M 2 s , there are no additional poles in the complex plane and thus only the massless graviton propagates.
We are interested in the non-relativistic, weak-field and static spacetime approximations, such that we can compute the gravitational potential from which, in turn, one can write down the Hamiltonian interaction coupling gravity and matter sectors in both cases of semi-classical and quantized gravitational interaction, as it has been done in Ref. [4].
As shown in Ref. [4], in the semi-classical approach gravity is coupled to the expectation value of the quantum energy-stress tensor, and the field equation for the potential, with the choice Eq. (9), reads 5 : whose solution is given by 5 Since terms with derivatives of order higher than four are usually neglected at small curvature, one could be brought to think that the Taylor expansion of the exponential e −∇ 2 /M 2 s can be truncated. However, it is not the case here: in fact, in Eq. (8) we are considered the most general quadratic-curvature action, and the infinite-order in derivatives comes from the form factors F i ( ) and not from higher order curvature-invariants. Moreover, even if we wanted to truncate the series, we would suffer from the ghost problem again.
i.e. one has a classical gravitational potential generated by the probability density |ψ( x, t)| 2 that plays the role of a semi-classical source. In a Hamiltonian formulation the potential in Eq. (12) contributes to the self-interaction of matter as described by the modified Schrödinger equation in Eq. (1), see [4].
In the case of quantized gravity instead, the graviton field is directly coupled to the quantum energy-stress tensor so that the analog of Eq. (11) reads whose solution is given bŷ (14) Note that we have usedτ 00 ≡ρ = mψ † ψ. By calculating also in this case the Hamiltonian interaction, it becomes clear that the quantum gravitational potential in Eq. (14) does not introduce any non-linearity in a Nparticle Schrödinger equation, but the non-linear integrodifferential equation (1) would emerge when considering mean-field regime for the many-body system (N → ∞) [4].
In the following section we will study the spreading solutions of Eq. (1) and, by comparing to the case of Newtonian gravity, we will be able to see which is the effect of non-locality on a self-gravitating wave-packet in IDG-theory. The analysis will hold for both cases of semiclassical and quantized gravity as the main dynamical equation is mathematically the same, i.e. Eq. (1).

III. SPREADING SOLUTIONS FOR A SELF-GRAVITATING WAVE-PACKET
We now wish to study non-stationary solutions of the non-linear integro-differential equation in Eq. (1) by working in a regime where non-linearity-effects can be considered sufficiently small such that there will be no gravity-induced collapse. Such a regime is the one described by the inequality in Eq. (6) which gives the range of masses (see Eqs. (4) and (5)) for which the self-gravitating wave-packet will not collapse, as the attractive contribution due to gravity is not dominating, but it can allow only the spread of the wave-packet. In this scenario, we are allowed to study Eq. (1) in the Fourier space.
Let us suppose we start with an initial Gaussian wavepacket: ψ( x, 0) = 1 where σ 0 is the initial spread. A formal expression of the wave-packet at a generic time t > 0 can be found in terms of its Fourier transform: where we have used φ( k, ω) = 2πφ( k)δ(ω − ω( k)), and φ( k) can be obtained by calculating the anti-Fourier transform at the initial time t = 0: where k ≡ | k|. By using Eq. (16), and the expression of the IDG potential in the momentum space, and by acting with´d 3 xdte −i( k· x−ωt) on both sides of the modified Schrödinger equation in Eq. (1), we obtain the dispersion-frequency ω IDG as a function of k: where (20) In order to solve the integral in Eq. (20) we can make the following change of integration variables: X := k − k , Y := k + k , thus k · k = (Y 2 − X 2 )/4 and the integral turns out to be decoupled in two other integrals that can be easily calculated by using polar coordinates: 6   note that in the regime M s σ 0 > 2, gives the corresponding Newtonian limit: In Fig. 1 it is shown the behavior of the dispersion relation ω(k) in the free, Newtonian and IDG cases, and one can immediately notice that as the parameter M s increases the frequency ω IDG tends to ω N . The dispersion relation in Eq. (22) is crucial in order to determine the time evolution of the wave-packet, indeed from Eq. (16) the solution ψ IDG ( x, t) is expressed in terms of the frequency ω IDG : (24) The integral in Eq. (24) cannot be solved analytically, but we will be able to find numerical solutions.
First of all, note that in the free-particle case there exists a well know analytical solution that describes a quantum-mechanical spreading of the wave-packet, and the corresponding probability density reads: from which one can see that there exist a time-scale for the spreading, i.e. a time after which the particle turns out to be de-localized, and it is given by: Note that the time-scale in Eq. (26) can be also obtained by imposing the equality τ free = 1/2ω free (1/σ 0 ). We argue that in the same way we can also obtain an analytical estimation for the spreading time-scale of a self-gravitating system, thus by using Eq. (22) and imposing τ IDG = 1/2ω IDG (1/σ 0 ) we obtain 7 that in the case of Newtonian self-interaction reduces to In the opposite regime M s σ 0 < 2, when non-localityeffects become dominant, the time-scale in Eq. (27) assume the following form: that is the same time-scale that we would have if there was a constant gravitational potential.
Since we are considering non-linear effects sufficiently small, in order to be consistent with the inequality in Eq. (6) we need to require that can be seen as a quantifier of non-linearity, as the one in Eqs. (6)- (8), and from which we can determine again a threshold value for the mass: that in the regime M s σ 0 > 2 reduces to the one of Newton's gravity: If we choose the values M s = 0.004 eV and σ 0 = 500 nm for the scale of non-locality and the initial spread, respectively, we obtain χ IDG 3.1 × 10 −17 kg and χ N 6.1 × 10 −18 kg. These values are of the same order of the threshold masses µ IDG and µ N that we have found above by equating kinetic and gravitational energy contributions. Thus, Eq. (30) is consistent with the inequality in Eq. (6).
It is very clear by comparing Eqs. (26), (27) and (28) that namely the gravitational self-interaction causes a slowdown of the spreading of the wave-packet compared to the free-particle case. Moreover, since at short distances IDG interaction is weaker than the Newtonian one, a selfgravitating wake-packet in the IDG case will spread more quickly compared to Newton's theory.
We have solved numerically the integral in Eq. (24) for both IDG and Newton's theory, and in Fig. 2 we have plotted the radial probability density ρ( x, t) = 4π| x| 2 |ψ( x, t)| 2 in the free, Newtonian and IDG cases at two fixed values of time as a function of the radial coordinate | x|. We can immediately notice that the results in    the two plots are in agreement with the analytical estimation made in Eq. (33). Such gravitational inhibitions of the spreading, not only would offer a way to explore classical and/or quantum properties of the gravitational interaction, but would also provide a new framework to test short-distance gravity beyond GR, by investigating the real nature of the gravitational potential.

IV. DISCUSSION
In the previous section we have found quantum spreading solutions for a self-gravitating wave-packet both in the case of IDG and Newtonian self-interactions. The results we have obtained, especially the analytical estimation in Eq. (33) and the numerical solutions in Fig. 2, provide us a unique window of opportunity to test modified theories of gravity in a laboratory, in particular IDGas an example of singularity-free theory of gravity.
Although the previous analysis holds for both cases of semi-classical and quantized gravity, in order to discuss the experimental testability of the model we need to distinguish between the two cases. As we have already mentioned above: 1. In the case of quantized gravitational interaction, Eq. (1) describes the dynamics of a condensate, where the mutual gravitational interaction of all components would give an effective self-potentialcontribution as result of the mean-field approximation; 2. In the semi-classical approach Eq. (1) can describe the dynamics of a self-gravitating one particle system, as for instance elementary particles, or molecules whose center-of-mass' dynamics would be taken into account.
The first case might be more interesting to explore: let us remind that the main motivations of IDG concern problems arising when trying to quantize the gravitational interaction, for example unitarity and renormalizability as already mentioned in Section II. Moreover, a semi-classical approach to IDG, would also imply a modification of quantum mechanics such that the fundamental equation governing the dynamics of a single-particle state would be non-linear. In this respect, case 1. seems less speculative than case 2.. However, in this section we wish to study the current and future experimental testability of our predictions about both classical and quantum aspects of gravity.
As far as case 1. is concern, there is a very promising experiment that is aimed to test quantum mechanics of weakly coupled Bose-Einstein condensates (BEC) in a freely falling system, where the spread of the quantum wave-packet would be tested in microgravity, see Ref. [40]. The used BEC was made of about 10 4 atoms, but technology is progressing and in the near future it will be possible to consider BEC with 10 6 , or even more, atoms, allowing us to compare our predictions with the experimental data in such a way to constrain the scale of non-locality M s , that so far has been only bounded in torsion-balance experiments [25].
Regarding case 2., molecule-interferometry [35] seems to be one of the most favorable scenario to verify predictions of the semi-classical approach. Although it is not the configuration analyzed in this paper, it is also worth mentioning that another promising scenario aimed to test semi-classical gravity is given by optomechanics as explained in Refs. [36][37][38], where one considers manybody systems with a well-localized wave-function for the center of mass.
All these kinds of experiments are very sensitive and, unfortunately, there are several sources of noise that need to be taken into account, as for example decoherence effects, see Ref. [39] for a review, and they represent a big challenge to overcome.
The choices of the initial spread σ 0 and the mass m are very crucial in order to determine the time-scale for the spreading. First of all, to have appreciable gravityinduced effects we need values of the mass that are not much smaller than the threshold mass. Moreover, if we take σ 0 = 500 nm and m ∼ χ N 6.0×10 −18 kg, the timescales are of the order of τ IDG , τ N ∼ O(10 5 )-O(10 6 ) sec and τ free ∼ O(10 4 ) sec. In an interferometric experimental setup, for example, these values would not be suitable to test short-distance gravity, because the time-scales are too large compared to the coherence-times achieved by a modern matter-wave interferometer (1-3 seconds) [41]. This means that we need smaller initial spreads so that the values of the time-scales will also decrease. For in-stance, if we choose σ 0 = 1 nm, m = 4.8 × 10 −17 kg and M s = 300 eV, the time-scales turn out to be of the order of τ free ∼ O(10 −1 ) sec, τ IDG ∼ O(1) sec and τ N ∼ O(10) sec, which are more suitable values of time to test and compare modified theories of gravity.
In the table I we have shown some values of the three time-scales for fixed values of initial spreads and masses. It is very clear that by decreasing the initial spread, the time-scales also decrease. Those numerical values apply to both cases 1. and 2.; of course the setup can be different depending on the kind of experiments, and so the preparation of the initial state will also differ. In tests with BEC [40] the experimental configuration is given by an asymmetric Mach-Zehnder interferometer, while molecular-interferometry can be performed, for instance, with a Talbot-Lau interferometer [35]. In both cases the spread σ 0 of the initial wave-packet is related in some way to the size of the slit separation d in the interferometric setup.
For example, as we have already mentioned in the footnote 3, in a Talbot-Lau interferometry setup the initial spread σ 0 is of the same order of the slit separation d, that in turn is expressed in terms of the length L of the interferometric device and of the the wave-length λ of the particle: L = d 2 /λ. The smallest slit separation that has been achieved so far is d = 500 nm [14], which also implies σ 0 ∼ 500 nm. However, we have seen that in order to build a very suitable experimental scenario in which we can test modified theories of gravity we need at least a value σ 0 = 1 nm for the initial spread, which means that technology should decrease at least of two orders of magnitude the slit-separation d. Moreover, we also need masses of the order of O(10 −18 -10 −17 ) kg, and this is one of the biggest challenge to overcome: in fact, the most massive systems for which interference patterns have been observed are organic molecules with a mass of 10 −22 -10 −21 kg, see Ref. [14]. In view of this last observation, it is worthwhile to mention that in Ref. [42] the authors present a new exciting proposal in which one might be able to perform quantum-interference with superconducting spheres with masses of the order of 10 −14 kg, and in such a regime both gravitational and quantum effects should be not negligible.
In this manuscript we have extended the window of opportunity to test classical and quantum properties of modified gravity provided in Ref. [4], where stationary properties of a quantum wave-packet were taken into account. Indeed, here we have made additional predictions that might allow us to test and constrain gravitational theories by studying non-stationary properties, in particular the quantum spreading of the wave-packet associated to a self-gravitating meso-scopic system in alternative theories beyond Einstein's GR. We have found a unique feature of IDG, as example of singularity-free theory of gravity, that predicts a faster spreading of the wave-packet compared to Newton's theory. A future observation of these predictions, even in a table-top experiment, might allow us a deeper and clearer understanding of short-distance gravity and let us learn more about its real nature, whether it is classical or quantum.