Experimental Parameters for Detection of Spin Entanglement via Quantum Gravity: Implications and Proposed Models

1. Introduction

In the realm of theoretical physics, quantum mechanics and quantum field theory have been able to describe three out of four fundamental known forces till now which are the electromagnetic interaction, the strong force, and the weak force. Gravity, however, still does not have a complete theoretical model. Our current model of gravity is based on Einstein’s General Theory of Relativity which has been remarkably accurate in describing many phenomena from the bending of light around mass to gravitational waves. However a fundamental problem we face is that space-time interacts with energy-momentum which are properties of fundamental particles which are currently best described by the quantum theory. General Relativity and Quantum Theory are incompatible models due to a variety of reasons [1], thus we must look for a way to introduce gravity to quantum physics.

In this paper we would look at an experimental proposal that aims to show that gravity can influence matter at quantum scales through the double Stern Gerlach Interferometer experiment devised by Bose et al in their paper [2]. This experiment might be able to show how gravitational field could cause entanglement between two particles. The first portion would look at the experimental parameters that would maximize spin entanglement in the experiment suggested by Bose et al under optimal conditions. Later we discuss what the implications of this experiment are and exploring further possible models.

2. Double Stern Gerlach Spin Entanglement

As suggested in Ref. [2], we keep two test masses in superposition of spacially localized states $|L\rangle$ and $|R\rangle$. Their evolution happens under gravitational interaction and all other forces are minimized. These two states are localized Gaussian wavepackets and due to their relatively small diameters with respect to their distances $\delta x$ and $\Delta x$, we can assume $\langle L|R\rangle =0$.

In this setup, we take two masses $m_1$ and $m_2$ with the centers of the left and right states of both the particles separated by the distance $\Delta x$ as shown in the figure.

Closest distance between ${|R\rangle }_1$ and ${|L\rangle }_2$ is $(\Delta x-\delta x)$ and farthest distance between ${|L\rangle }_1$ and ${|R\rangle }_2$ is $(\Delta x+\delta x)$.

The Classical gravitational potential energy is $E=-\frac{G{m}_1{m}_2}{d}$.

The initial state at $t=0$ is:

$${|\Psi (t=0)\rangle }_{12}=\frac{1}{\sqrt{2}}{\left(\right|L\rangle }_1+{|R\rangle }_1)\otimes \frac{1}{\sqrt{2}}{\left(\right|L\rangle }_2+{|R\rangle }_2)$$

(1)

$$=\frac{1}{2}{\left(\right|L\rangle }_1\otimes {|L\rangle }_2+{|L\rangle }_1\otimes {|R\rangle }_2+{|R\rangle }_1\otimes {|L\rangle }_2+{|R\rangle }_1\otimes {|R\rangle }_2)$$

Figure 1. Schematic depiction of the double Stern-Gerlach experiment as proposed by Bose et al. (a) Two mesoscopic particles of equal mass in this case in superposition of two spatially localized states $|L\rangle +|R\rangle$. (b) Evolution of states in SG interferometer into spin entangled states ${|L,\uparrow \rangle }_j+{|R,\downarrow \rangle }_j(j=1,2)$ Only the $T_{00}$ component of energy-momentum tensor is non-negligible due to speeds much slower than the speed of light c[2], hence only $h_{00}$ component evolves (where $h_{\mu \nu }$ is the small perturbation in linearized gravity $g_{\mu \nu }={\eta }_{\mu \nu }+h_{\mu \nu }$, ${\eta }_{\mu \nu }$ being the Minkowskian metric.)

Figure 1. Schematic depiction of the double Stern-Gerlach experiment as proposed by Bose et al. (a) Two mesoscopic particles of equal mass in this case in superposition of two spatially localized states $|L\rangle +|R\rangle$. (b) Evolution of states in SG interferometer into spin entangled states ${|L,\uparrow \rangle }_j+{|R,\downarrow \rangle }_j(j=1,2)$ Only the $T_{00}$ component of energy-momentum tensor is non-negligible due to speeds much slower than the speed of light c[2], hence only $h_{00}$ component evolves (where $h_{\mu \nu }$ is the small perturbation in linearized gravity $g_{\mu \nu }={\eta }_{\mu \nu }+h_{\mu \nu }$, ${\eta }_{\mu \nu }$ being the Minkowskian metric.)

Using the Schrödinger’s equation :

$$i\hslash \frac{\partial }{\partial t}|\Psi \rangle =\widehat{H}|\Psi \rangle$$

(2)

The general solution is found to be $|\Psi \left(t\right)\rangle =e^{-iHt/\hslash }|\Psi \left(0\right)\rangle$ when the Hamiltonian is time-independent and $\widehat{U}\left(t\right)=e^{-iHt/\hslash }$ is the unitary time-evolution operator.

As E is an eigenvalue of the Hamiltonian, we can write the final state as $|\Psi \left(t\right)\rangle =e^{-iEt/\hslash }|\Psi \left(0\right)\rangle$ At $t=\tau$, the state evolves to:

$${|\Psi (t=\tau )\rangle }_{12}=\frac{e^{i\varphi }}{2}{\left(\right|L\rangle }_1\otimes {|L\rangle }_2+e^{i\Delta {\varphi }_{LR}}{|L\rangle }_1\otimes {|R\rangle }_2+e^{i\Delta {\varphi }_{RL}}{|R\rangle }_1\otimes {|L\rangle }_2+{|R\rangle }_1\otimes {|R\rangle }_2)$$

(3)

Where ${\varphi }_{RL}\sim \frac{G{m}_1{m}_2\tau }{\hslash (\Delta x-\delta x)},{\varphi }_{LR}\sim \frac{G{m}_1{m}_2\tau }{\hslash (\Delta x+\delta x)},\varphi \sim \frac{G{m}_1{m}_2\tau }{\hslash \Delta x}$

$\Delta {\varphi }_{LR}={\varphi }_{LR}-\varphi$ and $\Delta {\varphi }_{RL}={\varphi }_{RL}-\varphi$.

A state is said to be entangled if it cannot be factorized (non-separable) to be written as a product of states such as $|\Psi \rangle ={|\psi \rangle }_1\otimes {|\psi \rangle }_2$ (a separable state).

If we try to factorize Eq(3), we obtain:

$$\frac{e^{i\varphi }}{2}{\left(\right|L\rangle }_1\otimes {\left(\right|L\rangle }_2+e^{i\Delta {\varphi }_{LR}}{|R\rangle }_2{)+|R\rangle }_1\otimes {\left(e^{i\Delta {\varphi }_{RL}}\right|L\rangle }_2+{|R\rangle }_2\left)\right).$$

The condition for this to be completely factorizable would be if ${\left(\right|L\rangle }_2+e^{i\Delta {\varphi }_{LR}}{|R\rangle }_2)\sim {\left(e^{i\Delta {\varphi }_{RL}}\right|L\rangle }_2+{|R\rangle }_2)e^{i\Delta {\varphi }_{RL}}{|L\rangle }_2+{|R\rangle }_2=e^{i\Delta {\varphi }_{RL}}{\left(\right|L\rangle }_2+e^{-i\Delta {\varphi }_{RL}}{|R\rangle }_2)$ means $e^{-i\Delta {\varphi }_{RL}}=e^{i(2n\pi -\Delta {\varphi }_{RL})}$ where n is an arbitrary integer.

Therefore for the state to be entangled, $2n\pi -\Delta {\varphi }_{RL}\ne \Delta {\varphi }_{LR}$

This gives us the entanglement condition:

$$\Delta {\varphi }_{LR}+\Delta {\varphi }_{RL}\ne 2n\pi$$

(4)

The SG interferometer accelerates the particles through a magnetic gradient for a time ${\tau }_{\mathrm{acc}}$ and then the magnetic field is switched off for a greater amount of time $\tau$ enough to entangle both the superposition states.

In this process, the two joint states of the test masses evolve as shown in Eq.(1) and Eq.(3) the orbital states ${|L\rangle }_j$ and ${|R\rangle }_j$ get replaced by the spin-orbital qubit states ${|L,\uparrow \rangle }_j$ and ${|R,\uparrow \rangle }_j$.(Ref.[2]).

Replacing the states in Eq.(3) with the spin-orbital states gives us:

$$|\Psi (t=\tau +{\tau }_{\mathrm{acc}}){\rangle }_{12}=\frac{e^{i\varphi }}{2}{\left(\right|L,\uparrow \rangle }_1\otimes {|L,\uparrow \rangle }_2+e^{i\Delta {\varphi }_{LR}}{|L,\uparrow \rangle }_1\otimes {|R,\downarrow \rangle }_2+e^{i\Delta {\varphi }_{RL}}{|R,\downarrow \rangle }_1\otimes {|L,\uparrow \rangle }_2+{|R,\downarrow \rangle }_1\otimes {|R,\downarrow \rangle }_2)$$

(5)

In the final step, the superposition is brought back through unitary transformations ${\left(\right|L,\uparrow \rangle }_j\to {|C,\uparrow \rangle }_j{,|R\downarrow \rangle }_j\to {|C,\downarrow \rangle }_j$:

$$|\Psi {(t=\tau +2{\tau }_{\mathrm{acc}}\rangle }_{12}=\frac{e^{i\varphi }}{2}{\left(\right|\uparrow \uparrow \rangle }_{12}+e^{i\Delta {\varphi }_{LR}}{|\uparrow \downarrow \rangle }_{12}+e^{i\Delta {\varphi }_{RL}}{|\downarrow \uparrow \rangle }_{12}+{|\downarrow \downarrow \rangle }_{12}{\left)\right|C\rangle }_1{|C\rangle }_2$$

(6)

This entanglement can be measured using an entanglement witness $\mathcal{W}$ which is discussed later in Appendix A. If $\mathcal{W}$ satisfies a certain condition for example exceeding unity in the one Bose et al. suggested, then the state is proven to be entangled.

3. Minimization of Casimir-Polder Interaction

While performing the experiment, it is important to reduce extra noise that may present itself due to other forces between the two particles we are using. One of them is the Casimir force which we will try to reduce.

We will be taking a linear setup for the experiment[3] and the particles as diamond nanoparticles of mass $\sim {10}^{-14}kg$.

The formula for the Casimir force between two neutral dielectric particles is

$$V_{CP}=-\frac{23\hslash c}{4\pi }\frac{R^6}{{(\Delta x-\delta x)}^7}{\left(\frac{\epsilon -1}{\epsilon +2}\right)}^2$$

(7)

Where R is the radius of the particles, $\Delta x-\delta x$ is the minimum distance between them, c is the speed of light in vacuum ($\sim 3\times {10}^8$m/s), ℏ is reduced Planck’s constant ($1.05\times {10}^{-34}$Js)and $\epsilon$ is the dielectric constant (∼ 5.7 for diamonds).

For a minimum difference of one order of magnitude[3], we obtain $0.1|V_G|\ge |V_{CP}|$

$$\begin{array}{c}\hfill \frac{G{m}^2}{\Delta x-\delta x}\ge \frac{230\hslash c}{4\pi }\frac{R^6}{{(\Delta x-\delta x)}^7}{\left(\frac{\epsilon -1}{\epsilon +2}\right)}^2\\ \hfill \Rightarrow \frac{G{(4\pi \rho R^3/3)}^2}{\Delta x-\delta x}\ge \frac{230\hslash c}{4\pi }\frac{R^6}{{(\Delta x-\delta x)}^7}{\left(\frac{\epsilon -1}{\epsilon +2}\right)}^2\\ \hfill \Rightarrow \frac{{(\Delta x-\delta x)}^7}{\Delta x-\delta x}\ge \frac{9}{G{\left(4\pi \rho R^3\right)}^2}\frac{23\hslash c{R}^6}{40\pi }{\left(\frac{\epsilon -1}{\epsilon +2}\right)}^2\\ \hfill \Rightarrow {(\Delta x-\delta x)}^6\ge \frac{230\hslash c}{4\pi G}{\left(\frac{3}{4\pi \rho }\frac{\epsilon -1}{\epsilon +2}\right)}^2\end{array}$$

$$\Delta x-\delta x\ge {\left[\frac{230\hslash c}{4\pi G}{\left(\frac{3}{4\pi \rho }\frac{\epsilon -1}{\epsilon +2}\right)}^2\right]}^{1/6}\equiv {(\Delta x-\delta x)}_{\min }\approx 157\mu m$$

(8)

Thus we need the minimum distance between the entangled particles to be $157\mu m$ for the Casimir force to be negligible (one order of magnitude lower than gravitational force in this case).

4. Conditions for Maximal Entanglement

We will be converting the quantum state to a density operator ($\widehat{\rho }$) in order to perform certain operations on it.

The state can be written as $|\Psi \rangle =\frac{e^{i\varphi }}{2}\left(\right|\uparrow \uparrow \rangle +e^{i\Delta {\varphi }_{LR}}|\uparrow \downarrow \rangle +e^{i\Delta {\varphi }_{RL}}|\downarrow \uparrow \rangle +|\downarrow \downarrow \rangle )$

Here $|\uparrow \rangle =\left(\begin{array}{c}1\\ 0\end{array}\right),|\downarrow \rangle =\left(\begin{array}{c}0\\ 1\end{array}\right)$ and $|\alpha \beta \rangle =|\alpha \rangle \otimes |\beta \rangle$

The density matrix is evaluated to be

$$\widehat{\rho }=\frac{1}{4}\left(\begin{array}{cccc}1& e^{-i\Delta {\varphi }_{LR}}& e^{-i\Delta {\varphi }_{RL}}& 1\\ e^{i\Delta {\varphi }_{LR}}& 1& e^{i(\Delta {\varphi }_{LR}-\Delta {\varphi }_{RL})}& e^{i\Delta {\varphi }_{LR}}\\ e^{i\Delta {\varphi }_{RL}}& e^{i(\Delta {\varphi }_{RL}-\Delta {\varphi }_{LR})}& 1& e^{i\Delta {\varphi }_{RL}}\\ 1& e^{-i\Delta {\varphi }_{LR}}& e^{-i\Delta {\varphi }_{RL}}& 1\end{array}\right)$$

The bipartite von Neumann entropy can help us find the maximally entangled state of a quantum system and is represented by

$$\mathcal{S}\left({\rho }_A\right)=-\mathrm{Tr}[{\rho }_{A}log{\rho }_A]=-\mathrm{Tr}[{\rho }_{B}log{\rho }_B]=\mathcal{S}\left({\rho }_B\right)$$

(9)

Where ${\rho }_A={\mathrm{Tr}}_B\left(\rho \right)$ and ${\rho }_B={\mathrm{Tr}}_A\left(\rho \right)$ are the reduced density matrices.

The reduced density matrix was calculated to be:

$${\widehat{\rho }}_1=\frac{1}{4}\left(\begin{array}{cc}2& e^{-i\Delta {\varphi }_{RL}}+e^{i\Delta {\varphi }_{LR}}\\ e^{i\Delta {\varphi }_{RL}}+e^{-i\Delta {\varphi }_{LR}}& 2\end{array}\right)$$

$$\mathcal{S}\left({\rho }_1\right)=-\mathrm{Tr}({\rho }_{1}log\left({\rho }_1\right))$$

The system is maximally entangled when ${\rho }_1$ and ${\rho }_2$ are maximally mixed[4] which is achieved when they are diagonal matrices with diagonal elements as $1/d$ where d is the dimension of the Hilbert space of the basis (which is 2 in our case).

To get these conditions, we set the non-diagonal entries to 0 first:

$e^{-i\Delta {\varphi }_{RL}}+e^{i\Delta {\varphi }_{LR}}=0$ and $e^{i\Delta {\varphi }_{RL}}+e^{-i\Delta {\varphi }_{LR}}=0$

This gives us two equations $\left\{\begin{array}{c}cos(-\Delta {\varphi }_{RL})+isin(-\Delta {\varphi }_{RL})+cos(\Delta {\varphi }_{LR})+isin(\Delta {\varphi }_{LR})=0\hfill \\ cos(\Delta {\varphi }_{RL})+isin(\Delta {\varphi }_{RL})+cos(-\Delta {\varphi }_{LR})+isin(-\Delta {\varphi }_{LR})=0\hfill \end{array}\right.$

They can be written as

$$\left\{\begin{array}{c}(cos\Delta {\varphi }_{RL}+cos\Delta {\varphi }_{LR})+i(sin\Delta {\varphi }_{LR}-sin\Delta {\varphi }_{RL})=0\hfill \\ (cos\Delta {\varphi }_{RL}+cos\Delta {\varphi }_{LR})+i(sin\Delta {\varphi }_{RL}-sin\Delta {\varphi }_{LR})=0\hfill \end{array}\right.$$

By collecting the real and imaginary parts of the expressions.

For $z\in \mathbb{C}$ to be 0, $\mathrm{Re}\left(z\right)=0$ and $\mathrm{Im}\left(z\right)=0$ (where Re and Im are the real an imaginary components of the complex number respectively)

Looking at the imaginary part, $sin\Delta {\varphi }_{LR}-sin\Delta {\varphi }_{RL}=sin\Delta {\varphi }_{RL}-sin\Delta {\varphi }_{LR}=0$

This implies that $sin\Delta {\varphi }_{LR}=sin\Delta {\varphi }_{RL}$ Which means either

$$\Delta {\varphi }_{LR}=\Delta {\varphi }_{RL}$$

(10)

or

$$\Delta {\varphi }_{LR}=(2k+1)\pi -\Delta {\varphi }_{RL}$$

(11)

where k is an arbitrary integer.

Substituting Eq(10) in the real part, $cos\Delta {\varphi }_{LR}+cos\Delta {\varphi }_{RL}=cos\Delta {\varphi }_{LR}+\Delta {\varphi }_{LR}=2cos\Delta {\varphi }_{LR}$

This is not a satisfactory condition as $2cos\Delta {\varphi }_{LR}=0$ for only some values of $\Delta {\varphi }_{LR}$

Thus we can assume, $\Delta {\varphi }_{LR}\ne \Delta {\varphi }_{RL}$ This leads us to using Eq(11) as follows:

$cos\Delta {\varphi }_{LR}+cos\Delta {\varphi }_{RL}=cos\Delta {\varphi }_{LR}+cos((2k+1)\pi -\Delta {\varphi }_{LR})=cos\Delta {\varphi }_{LR}-cos\Delta {\varphi }_{LR}=0$

And that satisfies the given conditions, therefore

$$\Delta {\varphi }_{LR}+\Delta {\varphi }_{RL}=(2k+1)\pi$$

(12)

This is our condition for non-diagonal elements being 0. Now to verify for maximally mixed state, we can substitute these values in the reduced density matrix.

$${\rho }_1=\frac{1}{4}\left(\begin{array}{cc}2& 0\\ 0& 2\end{array}\right)=\left(\begin{array}{cc}{2}^{-1}& 0\\ 0& 2^{-1}\end{array}\right)$$

So $\begin{array}{l}S\left({\rho }_1\right)=-\mathrm{Tr}(\frac{1}{2}\left(\begin{array}{cc}1& 0\\ 0& 1\end{array}\right)log\left(\begin{array}{cc}{2}^{-1}& 0\\ 0& 2^{-1}\end{array}\right))=-\mathrm{Tr}\left(\frac{1}{2}\left(\begin{array}{cc}1& 0\\ 0& 1\end{array}\right)\left(\begin{array}{cc}log\left(2^{-1}\right)& 0\\ 0& log\left(2^{-1}\right)\end{array}\right)\right)\\ =-\mathrm{Tr}\left(\frac{1}{2}\left(\begin{array}{cc}-log2& 0\\ 0& -log2\end{array}\right)\right)=-\mathrm{Tr}\left(\left(\begin{array}{cc}(-log2)/2& 0\\ 0& (-log2)/2\end{array}\right)\right)\\ =-(-\frac{log2}{2}-\frac{log2}{2})=log2=log(dim\left({\mathcal{H}}_1\right))\end{array}$

By analogy with the classical entropy formula, the entanglement entropy has the following bounds:

$$0\le S\left({\rho }_1\right)\le log\left(d_{{\mathcal{H}}_1}\right)$$

As the system saturates the upper pound, the state is maximally entangled.

Figure 2. Plot of Entanglement Entropy as a function of $\Delta {\varphi }_{LR}$(from 0 to $\pi /2$) and $\Delta {\varphi }_{RL}$(from $\pi /2$ to $\pi$). It can be seen that the entropy achieves its maximum value of $log2$ when the sum of the phases approach $\pi$.

Figure 2. Plot of Entanglement Entropy as a function of $\Delta {\varphi }_{LR}$(from 0 to $\pi /2$) and $\Delta {\varphi }_{RL}$(from $\pi /2$ to $\pi$). It can be seen that the entropy achieves its maximum value of $log2$ when the sum of the phases approach $\pi$.

Now we can substitute the values of the phases to get an expression for the distances. For simplicity, we will take $k=0$ in Eq(12).

$$\Delta {\varphi }_{LR}+\Delta {\varphi }_{RL}=\pi \Rightarrow {\varphi }_{LR}-\varphi +{\varphi }_{RL}-\varphi =\pi \frac{G{m}^2\tau }{\hslash (\Delta x+\delta x)}-\frac{G{m}^2\tau }{\hslash \Delta x}+\frac{G{m}^2\tau }{\hslash (\Delta x-\delta x)}-\frac{G{m}^2\tau }{\hslash \Delta x}=\pi$$

gives

$$\frac{G{m}^2\tau }{\hslash }(\frac{1}{\Delta x+\delta x}-\frac{2}{\Delta x}+\frac{1}{\Delta x-\delta x})=\pi$$

And finally,

$$\frac{2{\left(\delta x\right)}^2}{\Delta x({(\Delta x)}^2-{\left(\delta x\right)}^2)}=\frac{\pi \hslash }{G{m}^2\tau }$$

(13)

The expression for the distance between the left and right states is[2]:

$$\delta x\sim \frac{1}{2}\frac{g{\mu }_B{\partial }_{x}B}{m}{\tau }_{\mathrm{acc}}^2$$

(14)

Where $g\sim 2$ is the electronic g-factor, ${\mu }_B\approx 9.274·{10}^{-24}J{T}^{-1}$ is the Bohr magneton, ${\partial }_{x}B$ is the magnetic field gradient in the x-direction. The mass of the particles we are taking is $={10}^{-14}$kg, the magnetic field gradient is $={10}^{7}T{m}^{-1}$, ${\tau }_{\mathrm{acc}}$ is the time the masses are subjected to the magnetic field gradient and is a variable.

Substituting the distance for minimized CP-effect, $\Delta x-\delta x=157\mu m\Rightarrow \Delta x=\delta x+157\mu m$, we get:

$$\frac{2{\left(\frac{g{\mu }_B{\partial }_{x}B}{2m}{\tau }_{\mathrm{acc}}^2\right)}^2}{(\frac{g{\mu }_B{\partial }_{x}B}{2m}{\tau }_{\mathrm{acc}}^2+157\times {10}^{-6})({(\frac{g{\mu }_B{\partial }_{x}B}{2m}{\tau }_{\mathrm{acc}}^2+157\times {10}^{-6})}^2-{\left(\frac{g{\mu }_B{\partial }_{x}B}{2m}{\tau }_{\mathrm{acc}}^2\right)}^2)}=\frac{\pi \hslash }{G{m}^2\tau }$$

(15)

As the experiment cannot be completely free from outside interaction, there would always be some sort of decoherence in the system. The collisional decoherence time is approximately of the same magnitude as the total fall time of the particle ${\tau }_{\mathrm{decoherence}}\sim {\tau }_{\mathrm{total}}=\tau +2{\tau }_{\mathrm{acc}}$[2]. This will lead to a strong loss of coherence which will interfere with the gravitational entanglement. To minimize this, we can minimize the total time taken for the particle to fall (${\tau }_{\mathrm{total}})$

Minimizing ${\tau }_{\mathrm{total}}$ by a simple python program mentioned in Appendix C gives $\tau \approx 8.34349s$ and ${\tau }_{\mathrm{acc}}\approx 0.58829s$ Using these values, we can obtain the distances between the entangled states as

$$\delta x=3209.59\mu m\mathrm{and}\Delta x=3366.59\mu m$$

Calculating the phase values:

$${\varphi }_{LR}=\frac{G{m}^2\tau }{\hslash (\Delta x+\delta x)}=0.08059546$$

$${\varphi }_{RL}=\frac{G{m}^2\tau }{\hslash (\Delta x-\delta x)}3.37586158$$

$$\varphi =\frac{G{m}^2\tau }{\hslash \Delta x}=0.157432377$$

This gives $\Delta {\varphi }_{RL}=3.218429203,\Delta {\varphi }_{LR}=-0.076836917$

These values can be substituted in Eq(12) to verify that the maximal entanglement conditions are met:

$$\Delta {\varphi }_{RL}+\Delta {\varphi }_{LR}=3.141592286\approx \pi$$

Some rounding error of $1.17·{10}^{-5}$% persists but it is almost negligible.

5. Discussion

In order to make inferences from the proposed experiment, we will first need to know its assumptions and limitations.

The main theoretical assumptions of the Bose et al [2] paper is that the gravitational interaction between two masses is mediated by a gravitational field (not a direct action at a distance) and that entanglement between two systems cannot be created by Local Operations and Classical Communication (LOCC).

Another assumption made by them is using the Newtonian potential for mediating the gravitational entanglement which is not the true degree of freedom of the gravitational field and hence might have affected the implications of the experiment. (Discussed thoroughly in Ref. [5])

Based on the given assumptions, the entanglement generated by gravity would be non-classical. However, even if it eliminates the possibility of a completely classical model of gravity, it does not strictly imply that gravity is quantum and a semi-classical model of a gravity is still a valid possibility.

Due to the success of the general theory of relativity in explaining gravity, a good start to uncovering the mystery would be to look at the semi-classical models of gravity. The classical Einstein’s Field Equations are

$$G_{\mu \nu }=8\pi G{T}_{\mu \nu }$$

(16)

Where $G_{\mu \nu }$ is the Einstein Tensor and $T_{\mu \nu }$ is the stress-energy tensor. However in a quantum scenario, the stress-energy tensor has to be an operator ${\widehat{T}}_{\mu \nu }$ which would be inconsistent with the left half of the equation which is why we take its expectation value according to Møller-Rosenfeld [6,7], giving us a scalar value of the operator.

$$G_{\mu \nu }=8\pi G\langle {\widehat{T}}_{\mu \nu }\rangle$$

(17)

Here $\langle {\widehat{T}}_{\mu \nu }\rangle =\langle \psi |{\widehat{T}}_{\mu \nu }|\psi \rangle$. We do run into technical difficulties in this equation one of them being that the expectation value of the energy-momentum tensor is not conserved while the Einstein tensor is (${\nabla }^{\mu }\langle {\widehat{T}}_{\mu \nu }\rangle \ne 0,{\nabla }^{\mu }{G}_{\mu \nu }=0$). Thus we need to renormalize $T_{\mu \nu }$ which would give us ${\nabla }^{\mu }\langle T_{\mu \nu }\rangle =0$ and other conditions discussed by Wald [8].

Another flaw with the field equations semi-classical gravity is that the expectation value in some cases implies that the net curvature of the superposition states would be averaged out in the center and we would expect a test particle to move towards that center. This is not what experiments have shown so it is a massive downside of this model. ([9,10]) Due to the limitations of this model, other semi-classical models have also been proposed, an important one being the stochastic model [11].

Figure 3. a) Expected behaviour in Møller-Rosenfeld semi-classical gravity vs b) Observed behaviour of particle in superposition $|L\rangle +|R\rangle$

Figure 3. a) Expected behaviour in Møller-Rosenfeld semi-classical gravity vs b) Observed behaviour of particle in superposition $|L\rangle +|R\rangle$

Two of the current leading models for quantum gravity are the Loop Quantum Gravity model [12] and the String Theory [13] which would not be discussed here being beyond the scope of this paper.

In conclusion, new experiments might help us eliminate more erroneous models and get closer to the right one.