Quantum gravitational decoherence from fluctuating minimal length and deformation parameter at the Planck scale

Schemes of gravitationally induced decoherence are being actively investigated as possible mechanisms for the quantum-to-classical transition. Here, we introduce a decoherence process due to quantum gravity effects. We assume a foamy quantum spacetime with a fluctuating minimal length coinciding on average with the Planck scale. Considering deformed canonical commutation relations with a fluctuating deformation parameter, we derive a Lindblad master equation that yields localization in energy space and decoherence times consistent with the currently available observational evidence. Compared to other schemes of gravitational decoherence, we find that the decoherence rate predicted by our model is extremal, being minimal in the deep quantum regime below the Planck scale and maximal in the mesoscopic regime beyond it. We discuss possible experimental tests of our model based on cavity optomechanics setups with ultracold massive molecular oscillators and we provide preliminary estimates on the values of the physical parameters needed for actual laboratory implementations.

F ollowing early pioneering studies [1][2][3][4] , the investigation of the quantum-to-classical transition via the mechanism of decoherence has become a very active area of research, both experimentally and theoretically, playing an increasingly central role in the research area on the foundations of quantum mechanics (QM) and the appearance of a classical world at the macroscopic scale, as one may gather, e.g., from the many excellent existing reviews on the subject (i.e., see for instance refs. [5][6][7][8][9][10][11][12][13] and references therein). In broad terms, decoherence appears to be due to the inevitable interaction and the ensuing creation of entanglement between a given quantum system and the environment in which it is embedded. Indeed, no quantum system can truly be regarded as isolated, and the entanglement shared with the environment significantly affects the outcome of local measurements, even in the circumstance in which "classical" disturbances (such as dissipation and noise) can be neglected. Decoherence has been corroborated over the years by several experimental observations, among which it is worth mentioning its first detection obtained via cavity QED 14,15 , as well as the various laboratory tests involving superconductive devices 16 , trapped ions 17 and matter-wave interferometers 18 .
Apart from the domains of condensed matter and atomic physics, decoherence has been studied in more exotic and extreme contexts, ranging from particle physics to cosmological large-scale structures. For instance, in the framework of particle physics, decoherence has been investigated at the quantum field theory level in the regime of nontrivial interactions, both at zero and finite temperature 6,19,20 , and including lepton mixing and oscillations 21 . These studies and related ones have led to the opportunity of predicting new physics phenomenology, especially in connection with entanglement degradation in non-inertial frames and in the presence of strong gravitational fields 22 .
In principle, as in the above examples, the main features of each decoherence mechanism depend significantly on the type of interaction with the environment. On the other hand, in an attempt to explain wave function collapse and reduction to pointer states in general terms, Ghirardi, Rimini and Weber (GRW) proposed the existence of a universal decoherence mechanism and spontaneous localization with the introduction of just two model-dependent parameters, namely the collapse rate and the localization distance 23 (for important developments in GRW-like models, we refer the reader to refs. [24][25][26][27][28][29][30]. Inspired by the original GRW suggestion, Diosi 31 and Penrose 32 singled out classical Einstein gravity as the possible overarching interaction responsible for the loss of coherence in GRW theory, thus refining the previous intuition and fixing the free parameters in the GRW scheme of the quantum-to-classical transition. Since the innovative Diosi-Penrose insight, much effort has been devoted to experimentally test the gravitational realization of the GRW mechanism. In this regard, the attention has been mainly focused on gravitational time dilation 33 , optomechanics 34 and quantum clocks 35 . Very recently, a crucial experiment has demonstrated that the parameter-free Diosi-Penrose model based on classical Einstein gravity fails to account for the correct rates of radiation emission due to quantum state reduction, thus ruling it out as a feasible candidate for the implementation of the GRW decoherence mechanism 36 . The experimental falsification of the Diosi-Penrose version of GRW theory can be considered an important milestone as well as a new starting point. As a matter of fact, among the various ideas that revolve around the concept of gravitational decoherence (i.e., see for example refs. [37][38][39][40][41][42][43][44] ; for a comprehensive review and a more detailed list of references, see e.g., ref. 45 ), a significant portion deals with the high-energy regime in which quantum and gravitational effects are deemed to be comparably important. Although a complete and consistent theory of quantum gravity is still lacking, all the different current candidates including string theory, loop quantum gravity, noncommutative geometry and doubly special relativity, predict the existence of a minimal length at the Planck scale. An immediate consequence of this common aspect is the breakdown of the Heisenberg uncertainty principle (HUP), whose most famous formulation conveys that spatial resolution can be made arbitrarily small with a proper energetic probe. Therefore, as firstly derived in the framework of string theory 46 and strongly supported by thought experiments on micro black holes afterwards 47 , at the Planck scale the HUP must be superseded by a generalized uncertainty principle (GUP), whose minimal one-dimensional expression reads where β denotes the so-called deformation parameter and ℓ p is the Planck length. Starting from the above basic picture, the GUP has been discussed in the most disparate settings, from QM and quantum field theory [48][49][50][51][52][53][54] to noncommutative geometry and black hole physics [55][56][57][58][59][60] .
In this work, we introduce a model accounting for a universal quantum-gravitational decoherence process. The model assumes deformed canonical commutation relations (DCCRs) leading to the deformed uncertainty relation (1) that accounts for the presence of a minimal length scale. Furthermore, the model assumes the minimal length scale to be a fluctuating quantity, induced by a conjectured foamy character of quantum spacetime, whose average magnitude is fixed at the Planck length scale. As a consequence, the deformation parameter β that enters in the modified canonical momentum operator and in the DCCRs is promoted to a random quantity as well. As it is customary in the analysis of signatures of quantum gravitational effects at the classical macroscopic scale 61 , in the following we treat the contribution associated to the GUP as a small perturbation. We then show how to derive a master equation of the Lindblad-Gorini-Kossakowski-Sudarshan form 62,63 for the averaged quantum density matrix out of the corresponding βdeformed stochastic Schrödinger equation for the quantum state vector, and we study the physical consequences of such open quantum state dynamics. In light of the apparent experimental confutation of the Diosi-Penrose model, our proposal may be regarded as a possible alternative universal decoherence mechanism, provided that a set of reasonable assumptions holds. At the same time, a laboratory test of our model and its predictions would also yield an indirect evidence for the quantum nature of the gravitational interaction. Along these lines, building on a universal decoherence measure introduced in ref. 34 , we investigate a scheme for laboratory tests involving a cavity optomechanical system equipped with heavy molecular structures. We then provide numerical estimates and show that such a setup can probe our predictions with a high degree of accuracy.

Results
Deformed commutation relations, generalized uncertainty principle and modified Schrödinger dynamics. In order to define a consistent quantum mechanical theory in Hilbert space from Eq. (1), one can start from the β-deformed canonical commutation relation, as shown in ref. 64 . In compliance with the above discussions, we can then consider the deformation of the commutator (and, in turn, the modification of the uncertainty relations) as the signature of the quantum gravitational environment in which the system under examination is embedded. Notice that this kind of approach departs from refs. 42 and 43 which treat gravity as a classical decoherence source; at the same time, it also differs from ref. 41 that is based on perturbative quantum gravity.
For mirror-symmetric states, it is easy to prove that the uncertainty relation (1) can be straightforwardly derived from the following commutator: where the capital letters label the high-energy position and momentum operators, which are basically different from the usual (low-energy) ones of ordinary QM 51 . In the threedimensional case, the most general counterpart of the previous expression that preserves rotational isotropy is given by 65 : As frequently done in the literature 61,64 , in the following we will consider the case β 0 ¼ 2β; such a choice guarantees that the spatial geometrical structure is left commutative up to Oðβ; β 0 Þ. As a matter of fact, should a different selection be made for β 0 , we would havê which is subtler to handle, but eventually does not affect the final outcome, since the position operator does not play a relevant rôle in the upcoming discussions. It is worth remarking that from Eq.
(3) we can obtain an isotropic minimal uncertainty in position which is proportional to the Planck length, i.e., ΔX ' ffiffi ffi β p ' p . As already pointed out above, in light of Eq. (3) it is not possible to simply recover the position and momentum operators of ordinary QM, since they need to be properly modified so as to account for a non-vanishing value of the deformation parameter β. A suitable choice for these operators that complies with the DCCRs requires the identification: wherex j andp k satisfy the original HUP. The full expansion ofX j lies beyond the scope of the present paper; a more detailed treatment of these higher-order terms can be found in refs. 51 . For the sake of completeness, it must be noticed that, due to the existence of a minimal uncertainty in position, the operatorX j is no longer self-adjoint, but only symmetric 64 . Consequently, it is inappropriate to speak of position eigenstates, since the ensuing uncertainty is always non-vanishing. However, by relying on ΔX ≠ 0, one can build the so-called "quasi-position" space, which is made up of non-orthogonal states that represent the counterpart of the Dirac distribution δðx À x 0 Þ appearing in the standard QM framework. By thoroughly investigating the above setting, interesting phenomenological implications may be deduced even for simple problems, such as the particle in a box and the potential barrier 66 .
Collecting the above results, we can now seek the proper generalization of the Schrödinger time evolution equation in order to include the effects due to the DCCRs and the GUP. Starting from (5), it is immediate to achieve the modified quantum mechanical evolution equation for pure quantum states ψ in the form where we have omitted the hats on the operators to streamline the notation. By introducing H 0 = p 2 /2m and H = H 0 + V, the previous equation can be cast in the form where The magnitude of the deformation parameter β is commonly estimated to be of order unity 46,49,55,56,58,59 . From a more general perspective, β may be regarded as a dynamical variable 67 whose sign in several significant works is taken as either positive 46,47,49,56 or negative 55,57,59,60 . Note that, for the case β < 0, one can come across a "classical" regime at the Planck scale, that is, ΔXΔP ≥ 0. This possibility is predicted not only by the GUP, but also by doubly special relativity 68 as well as the cellular automaton formulation of QM due to 't Hooft 69 . Furthermore, the functional analysis of the position and momentum operators can be rigorously carried out also in this context, as shown in ref. 57 .
The above wide spectrum of estimates is compatible with the possibility that spacetime fluctuates as the Planck scale is approached, as predicted by major theoretical frameworks, such as the quantum foam scenario [70][71][72] and loop quantum gravity [73][74][75] . In turn, this suggests that space-time fluctuations fix the minimmal length scale only on average, thereby making the associated deformation parameter itself a fluctuating quantity. Equipped with a random β, Eq. (7) is promoted to a (linear) stochastic Schrödinger equation. Since the fluctuations in β should be related to the fluctuations of the metric tensor near the Planck threshold 71,72 , within the framework of non-relativistic QM it is natural to assume β to be a Gaussian white noise with a fixed mean and with a sharp auto-correlation so that where 〈 … 〉 denotes an average over fluctuations whose intensity is provided by the Planck time t p and β is the fixed mean. The consistency of our ansatz is supported by the observation that a white-noise process implies a large disproportion between the typical timescale of free motion and the auto-correlation time of the stochastic fluctuations. In addition, since quantum gravitational effects are expected to become maximally relevant at the Planck scale, the identification of t p as the correct time reference is tailor-made. This is a reasonable choice to describe Planck-scale effects but it is not the only one available, as one could always set the intensity of the fluctuations at the value α t p , where α would then be the remaining overall free parameter of the theory. In the following, we assume the previous parameter to be of order unity: α ≃ 1; by this choice, we assume the intensity of the fluctuations to be roughly given by the Planck time. Bearing this in mind, in order to make our model parameter-free, we need to fix the value of the constant mean appearing in Eq. (9). The original finding in the context of string theory concerning the magnitude of β being of order unity has since been confirmed by a large body of works starting from different frameworks; therefore, it is legitimate to set 〈β〉 = 1. Consequently, we obtain , by virtue of which we unambiguously identify our stochastic quantity.
Finally, from the β-deformed stochastic Schrödinger equation (7) for the state vector ψ , it is straightforward to derive the ensuing β-deformed stochastic Liouville-von Neumann equation for the dynamics of the quantum density matrix Equation (10) is the starting point for the analysis of the decoherence dynamics.
Fluctuating deformation parameter and induced master equation. Recalling the above hypothesis on the stochastic nature of the deformation parameter β entering in the DCCRs and the associated GUP, in the following we obtain a master equation in Lindblad form for the quantum density matrix averaged over the fluctuations of β. To this aim, we begin by working in the interaction picture via the standard unitary transformation which allows us to rewrite Eq. (7) as Consequently, after having restored the hidden dependence on the time t, the Liouville-von Neumann equation transforms accordingly as which has the formal exact solutioñ By inserting Eq. (14) into (13), one finds where in the integral in the r.h.s. we have made use of the Born-Markov approximation to let ϱ depend on t and not on t 0 . In order to provide the correct description in the mean of the quantum stochastic process under examination, we introduce the density matrix ρ = 〈ϱ〉 averaged over the fluctuations of the deformation parameter β. We remark that, being a convex combination of projectors, the average density matrix ρ is in general a mixed state subject to a non-unitary time evolution and hence to dissipation and decoherence. Bearing this in mind, averaging Eq. (15) by means of Eqs. (9) and keeping only the lowest-order contributions that do not exceed Oðβ 2 Þ, one has ∂ tρ ðtÞ ¼ Àσ It is worth noting that the intermediate result in Eq. (16) coincides with the one that is obtained exploiting the cumulant expansion method 76 , as applied in different contexts 38 . Finally, the Schrödinger representation for the average density matrix ρ is recovered by applying the transformation Assuming free motion (V = 0), the Lindblad-type master equation for the averaged density matrix ρ reads In the previous expression, the so-called dissipator is handily identified as the second term in the r.h.s. The Lindblad form of the master equation assures that the dynamical map is completely positive and trace preserving, thereby describing a bona fide quantum open system dynamics 77,78 . As required, in the limit σ ≈ 0 we recover the standard Liouville-von Neumann equation and the corresponding unitary dynamics. Furthermore, looking at Eq. (19), we see that one does not have to be concerned with the effects due to spatial non-commutativity; indeed, the position operator does not appear at all in the above equation, since the GUP correction only affects the momentum, which is the sole physical quantity that enters in H 0 .
From Eq. (19), we can also draw another interesting aspect by comparing the aforesaid expression with the outcome of previous works centered around an analogous issue. Indeed, in the dissipator the unperturbed Hamiltonian appears as H 2 0 , whereas in related treatments on gravitational decoherence (i.e., refs. 39,79 ) such dependence is linear. As a matter of fact, the present framework revolves around the existence of a minimal length ℓ p whilst in refs. 39,79 the original ansatz addresses the presence of a fundamental minimal time lapse. In principle, the two approaches might be reconciled provided that one resorts to a fully relativistic generalization of the DCCRs and of the GUP 51 to treat space and time on the same footing. However, in that case we would have to deal with two free parameters. Moreover, if we want to leave spacetime isotropy as well as the Poincaré algebra untouched, the choice of the aforementioned parameters is such that the nonrelativistic limit of the new DCCRS and of the new GUP significantly differs from Eq. (3) in the values of β and β 0 80 . In addition to that, under these circumstances we would have to work with fields instead of single-particle states, and it is not straightforward or even well-defined how to proceed in such an extremely complex scenario.
From the above analysis, we can conclude that, in order for quantum gravity to affect quantum coherence, one needs two basic ingredients: the existence of a minimal length at the Planck scale and a fluctuating GUP deformation parameter. In the following, we will study some of the most relevant quantities related to quantum gravitational decoherence as described by Eq. (19), specifically the entropy production and the decoherence time. In addition, we will discuss how the latter depends on the size and the characteristic energy scales for a variety of macroscopic and microscopic systems.
Physical results: entropy variation and decoherence time. In this section, we study the time evolution of the linear entropy and the quantum state purity, and we estimate the decoherence time for different physical systems. In particular, we find that the quantum gravitational decoherence associated with the GUP (1) entails a strong localization in energy space.
In what follows, we consider the linear entropy, defined as 81 : with trðρ 2 ðtÞÞ being the trace of the squared density matrix, that is the quantum state purity. Multiplying Eq. (20) by the constant factor d/(d − 1), with d being the Hilbert space dimension, the linear entropy can be normalized in such a way that SðtÞ 2 0; 1 ½ . Besides corresponding to the state mixedness, the linear entropy is also strictly connected to the von Neumann entropy, since the former represents the leading term of the latter in the Mercator logarithmic series expansion 82 around the pure state condition ρ 2 = ρ 83 .
The time evolution of the linear entropy reads Inserting Eq. (19) into the above expression and omitting the explicit time dependence, we have The first term in the r.h.s. identically vanishes due to the cyclic property of the trace, whereas the second contribution arising from the GUP must be manipulated a bit further. After some straightforward but tedious algebra, Eq. (22) can be recast in the form Now, by introducing the operator the function appearing in square brackets in Eq. (23) can be identified with which is the trace of a positive operator. Therefore, from Eq. (17) we have σ > 0, and hence As a result, the linear entropy is a monotonically increasing function of time and, viceversa, the quantum state purity is monotonically decreasing. Asymptotically, the averaged density matrix will tend to the maximally mixed state proportional to the identity, and the off-diagonal coherences will be completely washed out. A zero growth rate for the linear entropy can occur only in the instance [ρ, H 0 ] = 0, which yields the trivial solution corresponding to the time-invariant density matrix.
As hinted above, in order to estimate the decoherence time it is convenient to work in the momentum representation. By means of this choice, we have to work with ρ p;p 0 ðtÞ ¼ p ρðtÞ p 0 : ð27Þ By using the notation E(p) = p 2 /2m, Eq. (19) becomes whose solution is with ΔE 2 ¼ E 2 ðpÞ À E 2 ðp 0 Þ À Á . From Eq. (29), we observe that the time evolution conserves the diagonal elements whereas, as long as ΔE 2 ≠ 0 (i.e., no degeneracy), the off-diagonal terms of the density matrix decay exponentially, thereby realizing an effective localization in energy. The decay rate depends on both σ and the fourth power of the characteristic energy scale of the system. From Eq. (29), we can then identify the decoherence time τ D : The decoherence time is strongly influenced by the actual energy regime: the larger the deviation from the Planck scale is, the longer τ D becomes. We remark that Eq. (30) accounts for the size of the physical apparatus under consideration, as the squared mass appears explicitly in the denominator of the expression for τ D . In Fig. 1, we draw the contour plot of the decoherence time as a function of the mass and energy scales.
Our investigation suggests that the effects stemming from quantum gravity induce an efficient bridge between the quantum and the classical domain at currently reachable energies. In order to verify this claim, it is important to quantify the quantumgravitational decoherence times for typical classical and quantum systems at different scales. The physical dimensions and the corresponding τ D for the examples considered are reported in Table 1.
As long as classical systems are concerned, the quantumgravitational decoherence time is extremely short, and in one of the examples it is even close to the Planck time. This implies that quantum gravity could in principle be regarded as one of the main sources of the quantum-to-classical transition. Viceversa, when focusing on microscopic quantum systems, the results collected in Table 1 demonstrate that no coherence loss can be ascribed to the quantum nature of the gravitational interaction: the typical values for τ D in the last three cases are far longer than the currently estimated age of the Universe. This occurrence leaves no room for a gravity-induced classical behavior in microscopic quantum systems. This is to be expected, as it would be paradoxical to observe a quantum-to-classical transition at the microscopic level induced by quantum effects.
As a final step, to make Eq. (30) more transparent, we can write the decoherence time as a function of the Planck energy E p = m p c 2 and the energy of the system in the rest frame E = mc 2 . In so doing, we obtain In a similar fashion, also other decoherence times stemming from different physical settings 38,39,41,42 can be cast in the same form, which thus allows us to compare the predictions of our framework with those of the relevant gravitational decoherence models mentioned above. By focusing on the leading-order contributions only, we can build a table in which the main features of the different decoherence mechanisms are summarized.
By relying on the information contained in Table 2, we can pinpoint some interesting remarks. According to the above dependence of the various τ D on the energy regime E, it appears that our model provides the shortest decoherence time above a given energetic threshold (corresponding to a given mesoscopic scale), below which the scenario is completely reversed. As it might be expected, such an inversion occurs precisely at the Planck energy E = E p , with the ensuing decoherence time given by τ D = t p . This behavior is illustrated in Fig. 2, where we report τ D as a function of E. For the time being, we can observe that, since our model predicts the longest decoherence time in the quantum regime (below the mesoscopic threshold fixed by the Planck scale), it appears to be the one best suited for experimental verifications with quantum coherent systems.
Experimental verification. Concerning the experimental verification by laboratory tests of any of the gravitational decoherence mechanisms proposed thus far, the canonical route in order to detect such effects is to perform matter-wave interferometry with massive particles in superposition states 45 . However, since the seminal experiment with fullerene 84 and other relatively light molecules, as of today no sufficiently heavy mesoscopic system has been successfully prepared in tabletop laboratory tests. On the other hand, despite this persistent limitation, some significant improvements have been achieved in recent years; indeed, in less than a decade the largest molecule featuring quantum properties has moved from aggregates of <500 atoms 85 to structures made up of roughly 2000 atoms 86 , thus paving the way for further progress in the hopefully not too distant future.
Other practically viable avenues to test gravitational decoherence models rely on quantum simulations with atom-optical platforms 87,88 and analogue models 89 , which have been able to reproduce the spacetime curvature in the vicinity of black holes and wormholes.
Further potentially very efficient and promising schemes for the search of gravitational signatures in the decoherence processes are built upon quantum optical setups. The first proposal of a feasible experiment with a quantum optical system dates back to more then a decade ago 90,91 . According to these works, gravitational decoherence can be detected in near-Earth satellite tests involving entangled photons. The degradation of entanglement between the two photons is interpreted as a transfer of information from the bipartite quantum system to the environment via its coupling to each single photon, and the pervasive presence of gravity in empty interstellar space is regarded as the sole responsible for a   41 we have required the thermal background of gravitons to yield a contribution compatible with the energy of the analyzed quantum system, that is k B T ≃ E. similar occurrence. For more details on this research field, see the comprehensive review 92 ; for experimental developments and preliminary results, see ref. 93 . The same philosophy behind the quantum optical satellite probes is shared by laboratory tests of deformed commutation relations and other quantum gravity effects relying on systems of cavity optomechanics 94,95 , with the further advantage that optomechanical setups allow for the implementation of universal schemes for decoherence detection. Specifically, in ref. 34 the authors have introduced a very general method to verify any gravitational decoherence model starting from a readout of the loss of entanglement between two prepared subsystems. This phenomenon can be revealed by means of either quantum tomography or (as we will see below) Clauser-Horne-Shimony-Holt (CHSH) correlation measurements 96 . According to ref. 34 , a universal measure of decoherence is provided by means of the socalled min-entropy 97 . In a nutshell, from a physical point of view such an entropy accounts for the transfer (loss) of bipartite entanglement to a third party (an environment).
Formally, given a tripartite system made of distinguishable subsystems A, B and E, the degree of decoherence of A with respect to E is defined as 34 where R E!B is the set of all quantum operations from E to B, ρ AE is the reduced state of the bipartite system AE, i.e., the partial trace of the total density matrix ρ ABE with respect to B (i.e., ρ AE ¼ Tr B ðρ ABE Þ), Φ AB is a maximally entangled state on the bipartite system AB, 1 A is the unnormalized, fully decohered, maximally mixed state of subsystem A, and F is the Uhlmann fidelity, that is 81 Systems A and B are the two parties which have been initially prepared in an entangled state, whilst E is the environment. If there is some mechanism that determines a loss of coherence, the entanglement between A and B is constantly degraded and at the same time the entanglement between A and E is enhanced. In the absence of decoherence the maximum of the state fidelity is realized for A and B in a maximally entangled state, while for a fully decohered system the maximum of the fidelity is realized for A in the fully decohered state 1 A /d A , where d A is the dimension of the Hilbert space associated to subsystem A. This quantification of decoherence via the fidelity interplay between maximally entangled and maximally decohered states is physically very transparent. Actually, it was shown in ref. 34 that the decoherence measure Dec AjE ð Þ enjoys a remarkable relation to the minentropy H min AjE ð Þ ρ , that is the smallest Rényi conditional entropy providing a lower bound to the Shannon entropy of a statistical distribution ρ. Indeed, it turns out that 34,97 The detailed proof of this remarkable relation between entropy and decoherence is thoroughly discussed in ref. 34 .
Building on the above framework for the quantification of decoherence, we now discuss an optomechanical scheme for putting our model of quantum gravitational decoherence to the test. Let us consider two optomechanical systems, as shown in Fig. 3, where entangled photonic qubits are created in different conditions. In particular, one of the cavities possesses two fixed mirrors, while the other one is prepared so as to be subject to gravitational decoherence by allowing one of the mirrors to be movable 34 . In both cavities, the mechanical oscillator is an atom or a molecular structure trapped in a harmonic potential and the two systems are prepared initially in an entangled state. By applying an external laser field, these oscillators jump from the ground state to an excited level; they then decay back to the ground state by emitting the cavity radiation that allows for an accurate study of their vibrational modes.
In a realistic environment, besides gravitational decoherence, we must necessarily take into account the phenomenon of mechanical heating. Hence, when estimating the mean phonon number n of the mechanical oscillator, for the cavity with a movable mirror we have to include the simultaneous concurrence of the two effects, which yields where Λ grav = 1/τ D , Λ heat = k B T/ℏQ and γ m = ω m /Q expresses the ratio between the frequency ω m of the harmonic potential and the cavity quality factor Q. It is then necessary to reach very low temperatures and very high quality factors in order for the second term to be suppressed and achieve the sensitivity required to detect the conjectured gravitational effect. The required best conditions then go as follows: T = 10 nK, ω m = 1 s −1 , γ m = 10 −10 s −1 (and consequently Q = 10 10 ). These are still very demanding conditions; in particular, the required operating temperature is still some orders of magnitude below the lowest temperature T ≃ 0.1 mK achievable with current cooling technologies of optomechanical systems 98,99 .
Following the methods introduced in ref. 34 , a rigorous evaluation of the decoherence measure Dec(A|E) yields where g 0 = 1 s −1 is the single photon optomechanical coupling rate. Recalling the expression of the decoherence time, Eq. (31), we see that in order to achieve a gravitational decoherence rate comparable to or larger than the one due to the mechanical heating process one needs large molecular structures with mass of the order 10 −16 Kg (It is interesting to observe in passing that exactly this mass order of magnitude is the one required, in recent proposals for near-future tests of the quantum nature of the gravitational interaction 100,101 ). It is foreseen that the use of large molecular structures with masses of the order 10 −16 Kg will no longer represent an experimental limit in the next few years 100,101 .
In the above experimental conditions, quantum gravitational decoherence effects become clearly detectable and distinguishable from other environmental decoherence mechanisms. In Fig. 4 we report the behavior of the decoherence rate Dec(A|E) as a function of time both in the presence and in the absence of quantum gravitational decoherence. As expected, in the former situation the loss of coherence is faster, as there are two concurring decoherence rates (gravitational and due to mechanical heating).
The full-fledged experimental verification of our quantum gravitational decoherence model lies in the possibility of detecting the quantum correlations of the photon sources inside the cavities via CHSH inequality-type measurements. The interplay between Bell's theorem, the CHSH inequalities and gravity has an interesting story per se and has been addressed in several different contexts, see for instance refs. [102][103][104] . For the present goal, the Bell-type CHSH correlation to be evaluated reads For the effective two-level systems under consideration the standard CHSH observables are with σ x , σ y , σ z being the Pauli matrices.
As it can be deduced from Fig. 5, where we report the CHSH correlation ξ(t) using the same parameter values adopted in Fig. 4, the presence or the absence of our quantum gravitational decoherence mechanism strongly discriminates between different experimental outcomes. Indeed, the function ξ(t) for small (but accessible) times significantly departs from the value it would have if mechanical heating were the only source of decoherence. Therefore, the experimental scheme that we have illustrated appears suitable for the verification of our theoretical model.

Discussion
We have investigated the decoherence process associated with the existence of a minimal length at the Planck scale and the corresponding deformed quantum uncertainty relations. Assuming that the minimal spatial scale of length and the ensuing deformation parameter β are fixed only on average by space-time quantum fluctuations and are thus fluctuating random quantities, we have derived a master equation for the averaged quantum density matrix, thus showing that such quantum open dynamics can be cast in a Lindblad form, which guarantees complete positivity and trace preserving. The dissipator in the master equation depends on β in such a way to assure that the standard quantum mechanical dynamics is recovered in the limit β → 0. The evolutions of the linear entropy and quantum state purity are monotonically increasing and monotonically decreasing in time, respectively, yielding diagonal output states asymptotically. By resorting to the momentum representation, we have estimated the decoherence time and we have evaluated typical values in order of magnitude for some physical systems of varying mass and energy scales.
We have also proposed an experimental setup with which to test our predictions via cavity optomechanics at low temperatures performing CHSH correlation measurements. In order to discriminate the quantum gravitational decoherence effects from those due to mechanical heating, we need heavy molecular oscillators and ultracold temperatures that are beyond the reach of current technologies but could become available in the near future 100,101 .
In light of the above findings, we can conclude that the presence of a minimal length and a stochastic deformation parameter (as predicted by several leading candidates for a quantum theory of gravity) provides a valid decoherence mechanism capable of explaining some universal aspects of the quantum-to-classical Fig. 3 Optomechanical apparatus for the experimental test of quantum gravitational decoherence. The left cavity has two fixed mirrors, whereas the right cavity contains a movable mirror. The same external laser excites two atomic (or molecular) systems that act as sources of entangled photons when they decay back to their respective ground states emitting in the process entangled photon radiation into the cavity.  transition. By reasonably setting the mean deformation parameter at the constant value β ¼ 1= ffiffiffiffi t p p , we have cleared our model from any dependence on free parameters that would require to be constrained a posteriori by experimental data, as long as the would-be free parameter α controlling the fluctuation intensity is set to be of order one.
The working hypothesis of a stochastic contribution to the deformation parameter β is rooted in the dynamical scenarios for quantum gravitational effects 67 and the necessity to preserve the black hole complementarity principle. This ansatz leads to a transparent quantum-gravitational decoherence process that is quite immediate to discern from our formalism. For physical systems with mass and energy varying over a large spectrum, the values of the decoherence time τ D as summarized in Table 1 are entirely self-consistent and further corroborate the stochastic deformation picture. Furthermore, as reported in Fig. 2 and Table 2, comparison of our model with other relevant decoherence mechanisms proposed in the literature in recent years shows that our decoherence time is extremal, as it is the largest below the Planck scale and the smallest above it.
In summary, by assuming an emergent space-time equipped with a fluctuating fundamental length and a corresponding random deformation of the canonical commutation relations, our model yields a consistent and experimentally testable decoherence mechanism that explains the quantum-to-classical transition and provides a route for probing specific macroscopic consequences of the conjectured quantum nature of the gravitational interaction.
Concerning possible complementary and future developments, we recall that our quantum gravitational decoherence process occurs in energy space. On the other hand, the very same mechanism might arise also in connection with spatial localization. Indeed, the line of reasoning that we have carried out for the momentum-dependent GUP could be extended to include other features induced by quantum gravity, such as the so-called extended uncertainty principle, which accounts for the existence of a minimal momentum. Such a possibility was originally proposed in order to provide a foundation for the existence of the space-time curvature that becomes progressively more important at large distances 105 . The hypothesis of a minimal momentum scale is further supported by arguments based on the existence of a non-vanishing cosmological constant in the (anti)-de Sitter universe 106,107 . Following the same procedure developed in the present work, we plan to investigate the realization of a quantumgravitationally induced spatial localization mechanism.
On a final note, we remark that we have studied the implications of a GUP that is in fact the lowest order in a series of corrections to the Heisenberg uncertainty relation. Consequently, our results hold true as long as the quantum gravitational effects are not exceedingly strong. Such an instance is always realized at the energy scales that can be probed in the current laboratory tests, and hence our predictions are suited for the foreseeable experimental verifications of a universal decoherence process.
On the other hand, it is conceptually worthwhile to work out a complete theory capable of taking into account all the possible higher-order extensions of the uncertainty principle 60,108-110 . The introduction of additional corrections to the Heisenberg uncertainty relation would plausibly modify and refine various quantitative traits; however, such corrections are likely to preserve the general qualitative features outlined in the present work. We plan to tackle these issues in upcoming follow-on investigations.

Data availability
Data sharing not applicable to this article as no datasets were generated or analyzed during the current study.