Diósi–Penrose model

The Diósi–Penrose model was introduced as a possible solution to the measurement problem, where the wave function collapse is related to gravity. The model was first suggested by L. Diósi when studying how possible gravitational fluctuations may affect the dynamics of quantum systems.[1][2] Later, following a different line of reasoning, R. Penrose arrived at an estimation for the collapse time of a superposition due to gravitational effects, which is the same (within an unimportant numerical factor) as that found by Diósi, hence the name Diósi–Penrose model. However, it should be pointed out that while Diósi gave a precise dynamical equation for the collapse,[2] Penrose took a more conservative approach, estimating only the collapse time of a superposition.[3]

The Diósi model

In the Diósi model, the wave-function collapse is induced by the interaction of the system with a classical noise field, where the spatial correlation function of this noise is related to the Newtonian potential. The evolution of the state vector ${\displaystyle |\psi _{t}\rangle }$ deviates from the Schrödinger equation and has the typical structure of the collapse models equations:

${\displaystyle {\begin{aligned}d|\psi _{t}\rangle &={\Bigg [}-{\frac {i}{\hbar }}H\,dt+{\sqrt {\frac {G}{\hbar }}}\int d\mathbf {x} \,{\big (}{\mathcal {M}}(\mathbf {x} )-\langle {\mathcal {M}}(\mathbf {x} )\rangle _{t}{\big )}\,dW(\mathbf {x} ,t)\\&\qquad {}-{\frac {G}{2\hbar }}\int d\mathbf {x} \int d\mathbf {y} \,{\frac {{\big (}{\mathcal {M}}(\mathbf {x} )-\langle {\mathcal {M}}(\mathbf {x} )\rangle _{t}{\big )}{\big (}{\mathcal {M}}(\mathbf {y} )-\langle {\mathcal {M}}(\mathbf {y} )\rangle _{t}{\big )}}{|\mathbf {x} -\mathbf {y} |}}\,dt{\Bigg ]}|\psi _{t}\rangle ,\end{aligned}}}$

(1)

where

${\displaystyle {\mathcal {M}}(\mathbf {x} )=\sum _{j}m_{j}\mu _{R_{0}}(\mathbf {x} -{\hat {\mathbf {x} }}_{j})}$

(2)

is the mass density function, with ${\displaystyle m_{j}}$, ${\displaystyle {\hat {\mathbf {x} }}_{j}}$ and ${\displaystyle \mu _{R_{0}}(\mathbf {x} )}$ respectively the mass, the position operator and the mass density function of the ${\displaystyle j}$-th particle of the system. ${\displaystyle R_{0}}$ is a parameter introduced to smear the mass density function, required since taking a point-like mass distribution

${\displaystyle {\mathcal {M}}_{\text{point}}(\mathbf {x} )=\sum _{j=1}^{N}m_{j}\delta (\mathbf {x} -{\hat {\mathbf {x} }}_{j})}$

would lead to divergences in the predictions of the model, e.g. an infinite collapse rate[4][5] or increase of energy.[6][7] Typically, two different distributions for the mass density ${\displaystyle \mu _{R_{0}}(\mathbf {x} -{\hat {\mathbf {x} }}_{j})}$ have been considered in the literature: a spherical or a Gaussian mass density profile, given respectively by

${\displaystyle \mu _{R_{0}}^{\text{s}}(\mathbf {x} -{\hat {\mathbf {x} }}_{j})={\frac {3}{4\pi R_{0}^{3}}}\theta {\big (}|\mathbf {x} -{\hat {\mathbf {x} }}_{j}|-R_{0}{\big )}}$

and

${\displaystyle \mu _{R_{0}}^{\text{g}}(\mathbf {x} -{\hat {\mathbf {x} }}_{j})={\frac {1}{(2\pi R_{0}^{2})^{3/2}}}\,\exp \left(-{\frac {(\mathbf {x} -{\hat {\mathbf {x} }}_{j})^{2}}{2R_{0}^{2}}}\right).}$

Choosing one or another distribution ${\displaystyle \mu _{R_{0}}(\mathbf {x} -{\hat {\mathbf {x} }}_{j})}$ does not affect significantly the predictions of the model, as long as the same value for ${\displaystyle R_{0}}$ is considered. The noise field ${\displaystyle w(\mathbf {x} ,t):={\frac {dW(\mathbf {x} ,t)}{dt}}}$ in Eq. (1) has zero average and correlation given by

${\displaystyle \mathbb {E} \left[w(\mathbf {x} ,t)w(\mathbf {y} ,s)\right]={\frac {1}{|\mathbf {x} -\mathbf {y} |}}\delta (t-s),}$

(3)

where “${\displaystyle \mathbb {E} }$” denotes the average over the noise. One can then understand from Eq. (1) and (3) in which sense the model is gravity-related: the coupling constant between the system and the noise is proportional to the gravitational constant ${\displaystyle G}$, and the spatial correlation of the noise field ${\displaystyle w(\mathbf {x} ,t)}$ has the typical form of a Newtonian potential. Similarly to other collapse models, the Diósi–Penrose model shares the following two features:

• The model describes a collapse in position.
• There is an amplification mechanism, which guarantees that more massive objects localize more effectively.

In order to show these features, it is convenient to write the master equation for the statistical operator ${\displaystyle \rho (t)=\mathbb {E} {\big [}|\psi _{t}\rangle \langle \psi _{t}|{\big ]}}$ corresponding to Eq. (1):

${\displaystyle {\frac {d\rho (t)}{dt}}=-{\frac {i}{\hbar }}[H,\rho (t)]+{\frac {G}{\hbar }}\iint {\frac {d\mathbf {x} \,d\mathbf {y} }{|\mathbf {x} -\mathbf {y} |}}\left({\mathcal {M}}(\mathbf {x} )\rho (t){\mathcal {M}}(\mathbf {y} )-{\frac {1}{2}}\{{\mathcal {M}}(\mathbf {x} ){\mathcal {M}}(\mathbf {y} ),\rho (t)\}\right).}$

(4)

It is interesting to point out that this master equation has more recently re-derived by L. Diósi using an hybrid approach where quantized massive particles interact with classical gravitational fields.[8]

If one considers the master equation in the position basis, introducing ${\displaystyle \rho ({\vec {\boldsymbol {a}}},{\vec {\boldsymbol {b}}},t):=\langle {\vec {\boldsymbol {a}}}|\rho (t)|{\vec {\boldsymbol {b}}}\rangle }$ with ${\displaystyle |{\vec {\boldsymbol {a}}}\rangle$ :=|{\boldsymbol {a}}_{1}\rangle \otimes \dots \otimes |{\boldsymbol {a}}_{N}\rangle } , where ${\displaystyle |{\boldsymbol {a}}_{j}\rangle }$ is a position eigenstate of the ${\displaystyle j}$-th particle, neglecting the free evolution, one finds

${\displaystyle {\frac {d\rho ({\vec {\boldsymbol {a}}},{\vec {\boldsymbol {b}}},t)}{dt}}=-\Lambda ({\vec {\boldsymbol {a}}},{\vec {\boldsymbol {b}}})\rho ({\vec {\boldsymbol {a}}},{\vec {\boldsymbol {b}}},t)}$

(5)

with

${\displaystyle \Lambda ({\vec {\boldsymbol {a}}},{\vec {\boldsymbol {b}}})={\frac {G}{2\hbar }}\iint d\mathbf {x} \,d\mathbf {y} \,{\frac {{\big (}{\mathcal {M}}(\mathbf {x} ,{\vec {\boldsymbol {a}}})-{\mathcal {M}}(\mathbf {x} ,{\vec {\boldsymbol {b}}}){\big )}{\big (}{\mathcal {M}}(\mathbf {y} ,{\vec {\boldsymbol {a}}})-{\mathcal {M}}(\mathbf {y} ,{\vec {\boldsymbol {b}}}){\big )}}{|\mathbf {x} -\mathbf {y} |}},}$

(6)

where

${\displaystyle {\mathcal {M}}(\mathbf {x} ,{\vec {\boldsymbol {a}}}):=\sum _{j}m_{j}\mu _{R_{0}}(\mathbf {x} -{\boldsymbol {a}}_{j})}$

is the mass density when the particles of the system are centered at the points ${\displaystyle {\boldsymbol {a}}_{1}}$, ..., ${\displaystyle {\boldsymbol {a}}_{N}}$. Eq. (5) can be solved exactly, and one gets

${\displaystyle \rho ({\vec {\boldsymbol {a}}},{\vec {\boldsymbol {b}}},t)=\rho ({\vec {\boldsymbol {a}}},{\vec {\boldsymbol {b}}},0)e^{-t/\tau _{\text{DP}}({\vec {\boldsymbol {a}}},{\vec {\boldsymbol {b}}})},}$

(7)

where

${\displaystyle \tau _{\text{DP}}({\vec {\boldsymbol {a}}},{\vec {\boldsymbol {b}}})={\frac {1}{\Lambda ({\vec {\boldsymbol {a}}},{\vec {\boldsymbol {b}}})}}.}$

(8)

As expected, for the diagonal terms of the density matrix, when ${\displaystyle {\vec {\boldsymbol {a}}}={\vec {\boldsymbol {b}}}}$, one has ${\displaystyle \Lambda ({\vec {\boldsymbol {a}}},{\vec {\boldsymbol {a}}})=0}$, i.e. the time of decay goes to infinity, implying that states with well-localized position are not affected by the collapse. On the contrary, the off-diagonal terms ${\displaystyle {\vec {\boldsymbol {a}}}\neq {\vec {\boldsymbol {b}}}}$, which are different from zero when a spatial superposition is involved, will decay with a time of decay given by Eq. (8).

To get an idea of the scale at which the gravitationally induced collapse becomes relevant, one can compute the time of decay in Eq. (8) for the case of a sphere with radius ${\displaystyle R_{0}}$ and mass ${\displaystyle m}$ in a spatial superposition at a distance ${\displaystyle d:=|{\boldsymbol {a}}-{\boldsymbol {b}}|}$. Then the time of decay can be computed[9]) using Eq. (8) with

${\displaystyle \Lambda (d)={\begin{cases}{\dfrac {6GM^{2}}{5R_{0}\hbar }}\left({\dfrac {5}{3}}\lambda ^{2}-{\dfrac {5}{4}}\lambda ^{3}+{\dfrac {1}{6}}\lambda ^{5}\right)&{\text{when }}0\leq \lambda \leq 1,\\{\dfrac {6GM^{2}}{5R_{0}\hbar }}\left(1-{\dfrac {5}{12\lambda }}\right)&{\text{when }}\lambda \geq 1,\end{cases}}}$

(8)

where ${\displaystyle \lambda =d/(2R_{0})}$. To give some examples, if one considers a proton, for which ${\displaystyle m\simeq 1.67\times 10^{-27}}$ kg and ${\displaystyle R_{0}\simeq 10^{-15}}$ m, in a superposition with ${\displaystyle d\gg R_{0}}$, one gets ${\displaystyle \tau _{\text{DP}}\simeq 10^{6}}$ years. On the contrary, for a dust grain with ${\displaystyle m\simeq 6\times 10^{-12}}$ kg and ${\displaystyle R_{0}\simeq 10^{-5}}$ m, one gets one gets ${\displaystyle \tau _{\text{DP}}\simeq 10^{-8}}$ s. Therefore, contrary to what might be expected considering the weaknesses of gravitational force, the effects of the gravity-related collapse become relevant already at the mesoscopic scale.

Recently, the model have been generalized by including dissipative[7] and non-Markovian[10] effects.

Penrose's proposal

It is well known that general relativity and quantum mechanics, our most fundamental theories for describing the universe, are not compatible, and the unification of the two is still missing. The standard approach to overcome this situation is to try to modify general relativity by quantizing gravity. Penrose suggests an opposite approach, what he calls “gravitization of quantum mechanics”, where quantum mechanics gets modified when gravitational effects become relevant.[3][4][9][11][12][13] The reasoning underlying this approach is the following one: take a massive system well-localized states in space. In this case, being the state well-localized, the induced space–time curvature is well defined. According to quantum mechanics, because of the superposition principle, the system can be placed (at least in principle) in a superposition of two well-localized states, which would lead to a superposition of two different space–times. The key idea is that since space–time metric should be well defined, nature “dislikes” these space–time superpositions and suppresses them by collapsing the wave function to one of the two localized states.

To set these ideas on a more quantitative ground, Penrose suggested that a way for measuring the difference between two space–times, in the Newtonian limit, is

${\displaystyle \Delta E_{G}={\frac {1}{G}}\int d{\boldsymbol {r}}\,{\big (}g_{a}({\boldsymbol {r}})-g_{b}({\boldsymbol {r}}){\big )}^{2},}$

(9)

where ${\displaystyle g_{i}({\boldsymbol {r}})}$ is the Newtoninan gravitational acceleration at the point where the system is localized around ${\displaystyle i}$. The acceleration ${\displaystyle g_{i}({\boldsymbol {r}})}$ can be written in terms of the corresponding gravitational potentials ${\displaystyle \Phi _{i}({\boldsymbol {r}})}$, i.e. ${\displaystyle g_{i}({\boldsymbol {r}})=-\nabla \Phi _{i}({\boldsymbol {r}})}$. Using this relation in Eq. (9), together with the Poisson equation ${\displaystyle \nabla ^{2}\Phi _{i}({\boldsymbol {r}})=4\pi G\mu _{i}({\boldsymbol {r}})}$, with ${\displaystyle \mu _{i}({\boldsymbol {r}})}$ giving the mass density when the state is localized around ${\displaystyle i}$, and its solution, one arrives at

${\displaystyle \Delta E_{G}=4\pi G\int d{\boldsymbol {r}}\int d{\boldsymbol {r}}'\,{\frac {[\mu _{a}({\boldsymbol {r}})-\mu _{b}({\boldsymbol {r}})][\mu _{a}({\boldsymbol {r}}')-\mu _{b}({\boldsymbol {r}}')]}{|{\boldsymbol {r}}-{\boldsymbol {r}}'|}}.}$

(10)

The corresponding decay time can be obtained by the Heisenberg time–energy uncertainty:

${\displaystyle \tau _{G}={\frac {\hbar }{\Delta E_{G}}},}$

(11)

which, apart for a factor ${\displaystyle 8\pi }$ simply due to the use of different conventions, is exactly the same as the time decay ${\displaystyle \tau _{\text{DP}}}$ derived by Diósi's model. This is the reason why the two proposals are named together as the Diósi–Penrose model.

More recently, Penrose suggested a new and quite elegant way to justify the need for a gravity-induced collapse, based on avoiding tensions between the superposition principle and the equivalence principle, the cornerstones of quantum mechanics and general relativity. In order to explain it, let us start by comparing the evolution of a generic state in the presence of uniform gravitational acceleration ${\displaystyle {\boldsymbol {g}}}$. One way to perform the calculation, what Penrose calls “Newtonian perspective”,[4][9] consists in working in an inertial frame, with space–time coordinates ${\displaystyle ({\boldsymbol {r}},t)}$ and solve the Schrödinger equation in presence of the potential ${\displaystyle V({\boldsymbol {x}})=m{\boldsymbol {g}}\cdot {\boldsymbol {x}}}$ (typically, one chooses the coordinates in such a way that the acceleration ${\displaystyle {\boldsymbol {g}}}$ is directed along the ${\displaystyle z}$ axis, in which case ${\displaystyle V(z)=mgz}$). Alternatively, because of the equivalence principle, one can choose to go in the free-fall reference frame, with coordinates ${\displaystyle ({\boldsymbol {R}},T)}$ related to ${\displaystyle ({\boldsymbol {r}},t)}$ by ${\displaystyle {\boldsymbol {R}}={\boldsymbol {r}}-{\frac {1}{2}}{\boldsymbol {g}}t^{2}}$ and ${\displaystyle T=t}$, solve the free Schrödinger equation in that reference frame, and then write the results in terms of the inertial coordinates ${\displaystyle ({\boldsymbol {r}},t)}$. This is what Penrose calls “Einsteinian perspective”. The solution ${\displaystyle \Psi ({\boldsymbol {r}},t)}$ obtained in the Einsteinian perspective and the one ${\displaystyle \psi ({\boldsymbol {r}},t)}$ obtained in the Newtonian perspective are related to each other by

${\displaystyle \Psi ({\boldsymbol {r}},t)=e^{i{\frac {m}{\hbar }}\left({\frac {1}{6}}g^{2}t^{3}-{\boldsymbol {g}}\cdot {\boldsymbol {r}}\,t\right)}\psi ({\boldsymbol {r}},t).}$

(12)

Being the two wave functions equivalent apart for an overall phase, they lead to the same physical predictions, which implies that there are no problems in this situation, when the gravitational field has always a well-defined value. However, if the space–time metric is not well defined, then we will be in a situation where there is a superposition of a gravitational field corresponding to the acceleration ${\displaystyle {\boldsymbol {g}}_{a}}$ and one corresponding to the acceleration ${\displaystyle {\boldsymbol {g}}_{b}}$. This does not create problems as far as one sticks to the Newtonian perspective. However, when using the Einstenian perspective, it will imply a phase difference between the two branches of the superposition given by ${\displaystyle e^{i{\frac {m}{\hbar }}\left({\frac {1}{6}}({\boldsymbol {g}}_{a}-{\boldsymbol {g}}_{b})^{2}t^{3}-({\boldsymbol {g}}_{a}-{\boldsymbol {g}}_{b})\cdot {\boldsymbol {r}}\,t\right)}}$. While the term in the exponent linear in the time ${\displaystyle t}$ does not lead to any conceptual difficulty, the first term, proportional to ${\displaystyle t^{3}}$, is problematic, since it is a non-relativistic residue of the so-called Unruh effect: in other words, the two terms in the superposition belong to different Hilbert spaces and, strictly speaking, cannot be superposed. Here is where the gravity-induced collapse plays a role, collapsing the superposition when the first term of the phase ${\displaystyle {\frac {1}{6}}(g_{a}-g_{b})^{2}t^{3}}$ becomes too large.

Further information on Penrose's idea for the gravity-induced collapse can be also found in the Penrose interpretation.

Experimental tests and theoretical bounds

Since the Diósi–Penrose model predicts deviations from standard quantum mechanics, the model can be tested. The only free parameter of the model is the size of the mass density distribution, given by ${\displaystyle R_{0}}$. All bounds present in the literature are based on an indirect effect of the gravitational-related collapse: a Brownian-like diffusion induced by the collapse on the motion of the particles. This Brownian-like diffusion is a common feature of all objective-collapse theories and, typically, allows to set the strongest bounds on the parameters of these models. The first bound on ${\displaystyle R_{0}}$ was set by Ghirardi et al.,[6] where it was shown that ${\displaystyle R_{0}>10^{-15}}$ m to avoid unrealistic heating due to this Brownian-like induced diffusion. Then the bound have been further restricted to ${\displaystyle R_{0}>4\times 10^{-14}}$ m by the analysis of the data from gravitational wave detectors.[14] and later to ${\displaystyle R_{0}\gtrsim 10^{-13}}$ m by studying the heating of neutron stars.[15]

Regarding direct interferometric tests of the model, where a system is prepared in a spatial superposition, there are two proposals currently considered: an optomechanical setup with a mesoscopic mirror to be placed in a superposition by a laser,[16] and experiments involving superpositions of Bose–Einstein condensates.[9]

See also

• Measurement problem
• Interpretation of quantum mechanics
• Penrose interpretation
• Wave function collapse
• Objective-collapse theory
• Ghirardi–Rimini–Weber theory

References

1. Diósi, L. (1987-03-16). "A universal master equation for the gravitational violation of quantum mechanics". Physics Letters A. 120 (8): 377–381. doi:10.1016/0375-9601(87)90681-5. ISSN 0375-9601.
2. Diósi, L. (1989-08-01). "Models for universal reduction of macroscopic quantum fluctuations". Physical Review A. 40 (3): 1165–1174. doi:10.1103/PhysRevA.40.1165. PMID 9902248.
3. Penrose, Roger (1996-05-01). "On Gravity's role in Quantum State Reduction". General Relativity and Gravitation. 28 (5): 581–600. doi:10.1007/BF02105068. ISSN 1572-9532. S2CID 44038399.
4. Penrose, Roger (2014-05-01). "On the Gravitization of Quantum Mechanics 1: Quantum State Reduction". Foundations of Physics. 44 (5): 557–575. doi:10.1007/s10701-013-9770-0. ISSN 1572-9516.
5. Diósi, Lajos (2013-06-10). "Gravity-related wave function collapse: mass density resolution". Journal of Physics: Conference Series. 442: 012001. doi:10.1088/1742-6596/442/1/012001. ISSN 1742-6596.
6. Ghirardi, GianCarlo; Grassi, Renata; Rimini, Alberto (1990-08-01). "Continuous-spontaneous-reduction model involving gravity". Physical Review A. 42 (3): 1057–1064. doi:10.1103/PhysRevA.42.1057. PMID 9904128.
7. Bahrami, M.; Smirne, A.; Bassi, A. (2014-12-01). "Role of gravity in the collapse of a wave function: A probe into the Diósi-Penrose model". Physical Review A. 90 (6): 062105. arXiv:1408.6460. doi:10.1103/PhysRevA.90.062105. S2CID 118601408.
8. Diósi, Lajos (2011-07-08). "The gravity-related decoherence master equation from hybrid dynamics". Journal of Physics: Conference Series. 306: 012006. doi:10.1088/1742-6596/306/1/012006. ISSN 1742-6596.
9. Howl, Richard; Penrose, Roger; Fuentes, Ivette (2019-04-25). "Exploring the unification of quantum theory and general relativity with a Bose–Einstein condensate". New Journal of Physics. 21 (4): 043047. doi:10.1088/1367-2630/ab104a. ISSN 1367-2630. Eq. (23).
10. Gasbarri, G.; Toroš, M.; Donadi, S.; Bassi, A. (2017-11-13). "Gravity induced wave function collapse". Physical Review D. 96 (10): 104013. arXiv:1701.02236. doi:10.1103/PhysRevD.96.104013. S2CID 118881695.
11. Penrose, Roger; Mermin, N. David (December 1990). "The Emperor's New Mind: Concerning Computers, Minds, and the Laws of Physics". American Journal of Physics. 58 (12): 1214–1216. doi:10.1119/1.16207. ISSN 0002-9505.
12. Penrose, Roger, author. (31 March 2016). The Road To Reality : a Complete Guide to the Laws of the Universe. ISBN 978-1-4464-1820-8. OCLC 959367695.CS1 maint: multiple names: authors list (link)
13. Penrose, Roger (1998-08-15). "Quantum computation, entanglement and state reduction". Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences. 356 (1743): 1927–1939. doi:10.1098/rsta.1998.0256. ISSN 1364-503X. S2CID 83378847.
14. Helou, Bassam; Slagmolen, B. J. J.; McClelland, David E.; Chen, Yanbei (2017-04-28). "LISA pathfinder appreciably constrains collapse models". Physical Review D. 95 (8): 084054. doi:10.1103/PhysRevD.95.084054.
15. Tilloy, Antoine; Stace, Thomas M. (2019-08-21). "Neutron Star Heating Constraints on Wave-Function Collapse Models" (PDF). Physical Review Letters. 123 (8): 080402. doi:10.1103/PhysRevLett.123.080402. PMID 31491197. S2CID 119272121.
16. Marshall, William; Simon, Christoph; Penrose, Roger; Bouwmeester, Dik (2003-09-23). "Towards Quantum Superpositions of a Mirror". Physical Review Letters. 91 (13): 130401. arXiv:quant-ph/0210001. doi:10.1103/PhysRevLett.91.130401. PMID 14525288. S2CID 16651036.