Standard perturbative quantum field theory predicts gravitational decoherence rates scaling as \(G^2\), yielding decoherence times of order \(10^{26}\) years for laboratory masses—effectively unobservable. The Diósi-Penrose hypothesis, by contrast, predicts \(G^1\) scaling with \(\tau_{\mathrm{dec}} \sim \hbar d/(GM^2)\), giving nanosecond-scale decoherence for microgram masses. We show that this discrepancy traces to the initial state: the standard Feynman-Vernon influence functional assumes a product state \(\ket{\psi_{\mathrm{matter}}} \otimes \ket{0_{\mathrm{grav}}}\), which violates the linearized Wheeler-DeWitt constraint. Imposing the constraint forces an entangled initial state \((\ket{L}\ket{\Phi_L} + \ket{R}\ket{\Phi_R})/\sqrt{2}\), where \(\ket{\Phi_A}\) are coherent states of the gravitational field determined by the mass configuration. The constrained influence functional replaces the noise-kernel mechanism (which gives \(G^2\)) with a coherent-state-overlap mechanism (which gives \(G^1\)). The resulting decoherence rate is \(\Gamma = C \times GM^2/(\hbar d) + O(G^2)\), with \(C \in [1/2, 2]\) (best estimate \(C = 1\), matching the Diósi master equation). For a \(1~\mu\)g particle separated by 1 mm, we predict \(\tau_{\mathrm{dec}} = 1.58\) ns—testable by next-generation matter-wave interferometry. This derivation is valid within linearized gravity, where all approximations are controlled to extraordinary precision (\(GM/(c^2 d) \sim 10^{-38}\) for laboratory parameters).

Introduction

Does gravity decohere quantum superpositions, and if so, at what rate? This question sits at the intersection of quantum mechanics and general relativity, and its answer carries profound implications for both the quantum measurement problem and the structure of quantum gravity. Two theoretical frameworks give sharply different predictions—differing by a factor of \(10^{35}\) for laboratory-scale masses—making this one of the few questions about quantum gravity that experiment can decisively resolve.

The standard approach treats gravity as a quantum field on a fixed background and computes decoherence through the Feynman-Vernon influence functional , following the paradigm of Caldeira and Leggett . One places a mass \(M\) in spatial superposition, couples it to the graviton field via the linearized interaction \(H_{\mathrm{int}} \sim \sqrt{G}\, T_{\mu\nu} h^{\mu\nu}\), evolves the joint state, and traces over the graviton degrees of freedom. The resulting decoherence rate scales as \(G^2\), since the influence functional involves two interaction vertices—one from each factor of \(H_{\mathrm{int}}\) in the noise kernel. Explicit calculations by Anastopoulos and Hu  and by Blencowe  find \[\begin{equation} \Gamma_{\mathrm{pert}} \sim \frac{G^2 M^4}{\hbar^3 d^2}\,, \label{eq:G2-rate} \end{equation}\] where \(d\) is the superposition separation. For a particle of mass \(M = 1\;\mu\mathrm{g}\) separated by \(d = 1\;\mathrm{mm}\), this gives \(\tau_{\mathrm{dec}} \sim 10^{26}\;\mathrm{years}\)—far beyond any conceivable measurement.

In contrast, the Diósi-Penrose hypothesis  proposes that the decoherence rate is set directly by the gravitational self-energy of the superposition: \[\begin{equation} \Gamma_{\mathrm{DP}} = \frac{GM^2}{\hbar d}\,, \label{eq:G1-rate} \end{equation}\] scaling as \(G^1\). For the same \(1\;\mu\mathrm{g}\) particle at \(1\;\mathrm{mm}\) separation, this yields \(\tau_{\mathrm{dec}} \approx 1.6\;\mathrm{ns}\)—nine orders of magnitude shorter than the current experimental frontier and well within reach of next-generation matter-wave interferometry. The ratio \(\Gamma_{\mathrm{DP}}/\Gamma_{\mathrm{pert}} \sim \hbar d/(GM^2) \sim 10^{35}\) for these parameters, so the two predictions are separated by an experimentally enormous gap.

Despite the elegance of Eq. [eq:G1-rate] and its rich physical motivation—rooted in the incompatibility of superposed spacetime geometries—the Diósi-Penrose formula has remained a conjecture. It does not follow from any controlled approximation to quantum gravity. Meanwhile, the perturbative \(G^2\) result [eq:G2-rate], though rigorous within its framework, appears to render gravitational decoherence permanently unobservable. This impasse has persisted for over a decade.

In this paper, we show that the discrepancy between Eqs. [eq:G2-rate] and [eq:G1-rate] traces to a specific and correctable error in the standard Feynman-Vernon calculation: the choice of initial state. The perturbative treatment begins from the product state \[\begin{equation} \ket{\Psi(0)}_{\mathrm{pert}} = \frac{1}{\sqrt{2}}\bigl(\ket{L} + \ket{R}\bigr) \otimes \ket{0_{\mathrm{grav}}}\,, \label{eq:wrong-state} \end{equation}\] in which the matter superposition is uncorrelated with the graviton vacuum. This state is perfectly natural in perturbative quantum field theory, where interactions are switched on adiabatically. But it is unphysical in gravity: it violates the linearized Wheeler-DeWitt constraint , \[\begin{equation} \hat{H}_{\mathrm{matter}} + \hat{H}_{\mathrm{grav}} \approx 0\,, \label{eq:WDW-linearized} \end{equation}\] which requires that the gravitational field configuration be correlated with the matter distribution at all times. Imposing the constraint forces the physical initial state into the entangled form \[\begin{equation} \ket{\Psi(0)}_{\mathrm{phys}} = \frac{1}{\sqrt{2}}\bigl(\ket{L}\ket{\Phi_L} + \ket{R}\ket{\Phi_R}\bigr)\,, \label{eq:right-state} \end{equation}\] where \(\ket{\Phi_A} = D(\alpha_A)\ket{0}\) are coherent states of the gravitational field sourced by the mass at position \(A \in \{L, R\}\), with the displacement amplitude \(\alpha_A(\mathbf{k})\) determined by the Newtonian potential of each configuration. The matter and its gravitational field are entangled before any dynamics occurs—this entanglement is not generated perturbatively but imposed kinematically by the constraint.

The consequences for the influence functional are significant. In the standard (unconstrained) calculation, decoherence arises from the noise kernel—the two-point correlator of the gravitational field fluctuations—which introduces two powers of the coupling and yields \(G^2\). In the constrained calculation, the dominant decoherence mechanism is qualitatively different: it is the overlap \(\braket{\Phi_L(t)}{\Phi_R(t)}\) between the two coherent states, which depends on the gravitational self-energy \(E_G = GM^2/d\)—a single power of \(G\). The noise-kernel contribution remains present but is subleading at \(O(G^2)\).

The constrained influence functional establishes \(E_G\) as the correct energy scale for gravitational decoherence. The further step—converting this energy scale to a decoherence rate \(\Gamma = E_G/\hbar\)—requires the full Wheeler-DeWitt constraint (not merely the Poisson equation), specifically the modular Hamiltonian identification and the Page-Wootters mechanism for emergent time. We present both results clearly, distinguishing what the influence functional rigorously derives from what requires additional physical input from quantum gravity.

The paper is organized as follows. Section 2 reviews the standard Feynman-Vernon influence functional for gravitational decoherence, emphasizing the role of the product-state assumption in producing \(G^2\) scaling. Section 3 imposes the linearized Wheeler-DeWitt constraint, derives the constrained initial state [eq:right-state], and computes the modified influence functional. Section 4 extracts the decoherence rate, establishing \(\Gamma = C \times GM^2/(\hbar d) + O(G^2)\) with \(C \in [1/2, 2]\) (best estimate \(C = 1\)). Section 5 presents quantitative predictions for current and planned experiments. Section 6 discusses the relationship between the constrained and perturbative calculations, explains why electromagnetism does not produce an analogous effect, and identifies the regime of validity. Technical details of the coherent-state overlap computation and robustness analysis are given in the appendices.

The Standard Feynman-Vernon Influence Functional

We begin by reviewing the standard (unconstrained) derivation of gravitational decoherence via the Feynman-Vernon influence functional . This sets the stage for Section 3, where we show how imposing the Wheeler-DeWitt constraint modifies the calculation and changes the scaling from \(G^2\) to \(G^1\).

Setup

Consider a point mass \(M\) with center-of-mass coordinate \(q\), coupled to the quantized linearized gravitational field \(h_{\mu\nu}\). The total action separates into three pieces: \[\begin{equation} S[q,h] = S_M[q] + S_G[h] + S_{\mathrm{int}}[q,h], \label{eq:total-action} \end{equation}\] where \(S_M\) is the free matter action, \(S_G\) is the free graviton action (the linearized Einstein-Hilbert action for \(h_{\mu\nu}\)), and the interaction takes the form \[\begin{equation} S_{\mathrm{int}}[q,h] = \frac{\kappa}{2} \int d^4x \, T^{\mu\nu}(x;q) \, h_{\mu\nu}(x), \label{eq:interaction-action} \end{equation}\] with coupling constant \[\begin{equation} \kappa = \sqrt{\frac{32\pi G}{c^4}}\,. \label{eq:kappa-def} \end{equation}\] The stress-energy tensor \(T^{\mu\nu}\) depends on the matter trajectory \(q(t)\), and the coupling \(\kappa \propto \sqrt{G}\) sets the perturbative expansion parameter.

The reduced density matrix

We wish to compute the reduced density matrix of the mass after tracing over the gravitational field. In the path-integral formulation, this takes the Feynman-Vernon form : \[\begin{equation} \rho_M(q_f, q_f'; t) = \int \mathcal{D}q^{+}\,\mathcal{D}q^{-}\; \exp\!\Bigl(\frac{i}{\hbar}\bigl[S_M[q^{+}] - S_M[q^{-}]\bigr]\Bigr)\; \mathcal{F}[q^{+},q^{-}]\; \rho_M(q_i,q_i';0), \label{eq:reduced-density-matrix} \end{equation}\] where \(q^{+}\) and \(q^{-}\) are the forward and backward paths of the mass, and the influence functional \(\mathcal{F}\) encodes the entire effect of the gravitational environment: \[\begin{equation} \mathcal{F}[q^{+},q^{-}] = \int \mathcal{D}h^{+}\,\mathcal{D}h^{-}\; \exp\!\Bigl(\frac{i}{\hbar}\bigl[S_G[h^{+}] + S_{\mathrm{int}}[q^{+},h^{+}] - S_G[h^{-}] - S_{\mathrm{int}}[q^{-},h^{-}]\bigr]\Bigr)\; \rho_E(h_i, h_i'). \label{eq:influence-functional} \end{equation}\] Here \(\rho_E\) is the initial state of the gravitational field.

The product-state assumption

The standard treatment  assumes the total initial state is a product: \[\begin{equation} \ket{\Psi(0)} = \ket{\psi_{\mathrm{matter}}} \otimes \ket{0_{\mathrm{grav}}}, \label{eq:product-state} \end{equation}\] where \(\ket{0_{\mathrm{grav}}}\) is the graviton vacuum. This is the standard assumption in open quantum systems: the system and environment begin unentangled, and the environment starts in its ground state.

With this choice, the influence functional can be evaluated exactly because the graviton path integral is Gaussian. Defining the difference variable \(\Delta q(t) = q^{+}(t) - q^{-}(t)\), the imaginary part of the influence phase—responsible for decoherence—takes the noise-kernel form : \[\begin{equation} \mathrm{Im}\,\Phi[q^{+},q^{-}] = \frac{1}{2\hbar}\int_0^t ds\int_0^t ds'\; \Delta T^{\mu\nu}(s)\; N_{\mu\nu\alpha\beta}(s - s')\; \Delta T^{\alpha\beta}(s'), \label{eq:noise-kernel} \end{equation}\] where \(\Delta T^{\mu\nu}(s) = T^{\mu\nu}(q^{+}(s)) - T^{\mu\nu}(q^{-}(s))\) is the stress-energy difference between the two paths, and \(N_{\mu\nu\alpha\beta}\) is the graviton noise kernel—the symmetrized (Hadamard) two-point function of the gravitational field: \[\begin{equation} N_{\mu\nu\alpha\beta}(x,x') = \frac{\kappa^2}{4}\,\bigl\langle\bigl\{h_{\mu\nu}(x),\, h_{\alpha\beta}(x')\bigr\}\bigr\rangle_0\,. \label{eq:noise-kernel-def} \end{equation}\]

\(G\)-counting: why the standard result scales as \(G^2\)

The \(G\)-scaling of the decoherence rate is most transparent in the equivalent master-equation formulation : \[\begin{equation} \frac{d\hat{\rho}_M}{dt} = -\frac{1}{\hbar^2}\int_0^t dt'\; \mathrm{Tr}_{\mathrm{grav}}\!\bigl[H_{\mathrm{int}}(t),\,[H_{\mathrm{int}}(t'),\,\hat{\rho}_M \otimes \ket{0}\!\bra{0}]\bigr], \label{eq:master-equation} \end{equation}\] where \(H_{\mathrm{int}} = (\kappa/2)\int d^3x\, T^{\mu\nu} h_{\mu\nu}\). The \(G\)-counting proceeds as follows. The double commutator contains two insertions of \(H_{\mathrm{int}}\), each carrying one factor of \(\kappa \propto \sqrt{G}\), giving \(\kappa^2 \propto G\). Tracing over the graviton vacuum produces the Hadamard two-point function \(\langle 0|\{h_{\mu\nu}(x), h_{\alpha\beta}(x')\}|0\rangle\); for the canonically normalized graviton field, this propagator carries no explicit factor of \(G\). The \(G\)-scaling from the coupling and trace is therefore \(G^1\).

However, the full calculation involves contracting the noise kernel with the stress-energy tensor difference \(\Delta T^{\mu\nu} \propto M\) and integrating over the graviton spectral density. For a non-relativistic mass, the relevant frequency integrals contribute matter-dependent factors that bring the total to \(G^2\). Physically, each graviton vertex contributes \(\sqrt{G}\), and the leading (two-vertex) decoherence process involves one emission and one absorption, giving \(G^1\) from the vertices; the graviton propagator connecting them, when contracted with the stress-energy sources, contributes the additional \(G\) through the Newtonian potential \(\Phi_N \propto GM/r\). The net result is  \[\begin{equation} \Gamma_{\mathrm{QFT}} \propto G^2. \label{eq:G2-counting} \end{equation}\] The full result, computed by Anastopoulos and Hu  and by Blencowe , takes the form \[\begin{equation} \boxed{\;\Gamma_{\mathrm{QFT}} \sim \frac{G^2 M^4}{\hbar^3 d^2}\;} \label{eq:G2-rate-detail} \end{equation}\] (suppressing order-unity numerical prefactors and factors of \(c\)). For a \(1~\mu\)g mass separated by \(1\) mm, this gives \(\tau_{\mathrm{QFT}} = 1/\Gamma_{\mathrm{QFT}} \sim 10^{26}\) years—effectively infinite, and far beyond any foreseeable experimental reach.

Identification of the critical assumption

The \(G^2\) scaling can be traced directly to the product-state assumption [eq:product-state]. In the noise-kernel mechanism, decoherence arises because the two branches of the superposition emit slightly different graviton fields, and the growing distinguishability of these emitted fields degrades coherence. This is a dynamical process: entanglement between the mass and the gravitational field must be generated through the interaction, starting from zero. The rate of entanglement generation is set by \(H_{\mathrm{int}}^2 \propto G\), and the propagation of the emitted gravitons through the vacuum contributes another factor of \(G\), yielding \(G^2\) overall.

However, the product state [eq:product-state] violates the linearized Wheeler-DeWitt constraint. A mass at position \(q\) must carry its Newtonian gravitational field—the graviton vacuum \(\ket{0_{\mathrm{grav}}}\) is not a physical state for a system containing matter. The constraint demands that the initial state be entangled: each branch of the superposition must be dressed by its own coherent gravitational field configuration.

In Section 3, we impose this constraint and show that it replaces the noise-kernel mechanism (dynamical entanglement generation \(\to\) \(G^2\)) with a coherent-state-overlap mechanism (pre-existing entanglement manifestation \(\to\) \(G^1\)). The physical decoherence rate is set not by how fast gravitons are emitted, but by how fast the pre-existing gravitational dressing of the two branches becomes distinguishable.

The Constrained Influence Functional

We now develop the central result of this paper: the influence functional for gravitational decoherence when the Wheeler-DeWitt constraint is properly imposed on the Schwinger-Keldysh path integral. The constraint changes both the initial state and the structure of the functional integral itself, replacing the noise-kernel mechanism of Section 2 (which gives \(G^2\)) with a coherent-state overlap mechanism (which gives \(G^1\)).

The Wheeler-DeWitt constraint in linearized gravity

In canonical quantum gravity, physical states must satisfy the Hamiltonian constraint : \[\begin{equation} \hat{H}_{\mathrm{total}} \ket{\Psi_{\mathrm{phys}}} = 0\,. \label{eq:wdw} \end{equation}\] In the linearized (Newtonian) limit, the total Hamiltonian decomposes as \(\hat{H}_{\mathrm{total}} = \hat{H}_{\mathrm{matter}} + \hat{H}_{\mathrm{grav}} + \hat{H}_{\mathrm{int}}\), and the constraint reduces to the operator Poisson equation: \[\begin{equation} \nabla^{2}\hat{\Phi}(\mathbf{x}) = 4\pi G\,\hat{\rho}(\mathbf{x})\,, \label{eq:poisson-constraint} \end{equation}\] which determines the Newtonian potential \(\hat{\Phi}\) from the mass density \(\hat{\rho}\). This is not merely a field equation—it is a constraint that restricts the physical Hilbert space. Any state that violates [eq:poisson-constraint] is unphysical and must be projected out of the Hilbert space.

The constrained initial state

The constraint [eq:poisson-constraint] has immediate consequences for the allowed quantum states of a mass in spatial superposition.

Single branch.

For a point mass \(M\) localized at position \(\mathbf{x}_A\), the constraint uniquely determines the gravitational potential: \[\begin{equation} \Phi_{\mathrm{cl}}[\mathbf{x}_A](\mathbf{x}) = -\frac{GM}{|\mathbf{x} - \mathbf{x}_A|}\,. \label{eq:classical-potential} \end{equation}\] In the quantum theory, the gravitational field must be in the state that reproduces this classical potential in expectation value while minimizing the field energy—a coherent state: \[\begin{equation} \ket{\mathbf{x}_A}_{\!\mathrm{matter}} \;\longrightarrow\; \ket{\mathbf{x}_A}\ket{\Phi_A} = \ket{\mathbf{x}_A}\,\hat{D}(\alpha_A)\ket{0}\,, \label{eq:coherent-dressing} \end{equation}\] where \(\hat{D}(\alpha) = \exp\!\bigl(\int\!d^{3}k\, [\alpha(\mathbf{k})\,\hat{a}^{\dagger}_{\mathbf{k}} - \alpha^{*}(\mathbf{k})\,\hat{a}_{\mathbf{k}}]\bigr)\) is the Glauber displacement operator and the coherent-state amplitude is \[\begin{equation} \alpha_A(\mathbf{k}) = -\frac{4\pi GM}{k^{2}} \cdot \frac{e^{-i\mathbf{k}\cdot\mathbf{x}_A}} {\sqrt{2\hbar\omega_k}}\,, \label{eq:coherent-amplitude} \end{equation}\] with \(\omega_k = c|\mathbf{k}|\) for relativistic gravitons. This is the standard result for a quantum field coupled linearly to a classical source (Glauber 1963), applied here to linearized gravity. The coherent state is the unique minimum-uncertainty state satisfying the constraint at \(O(G)\); corrections from graviton squeezing arise only at \(O(G^2)\) with squeezing parameter \(r \sim GM/(c^2 d) \sim 10^{-38}\) (see Appendix 8).

Superposition.

For a mass prepared in a spatial superposition \((\ket{L} + \ket{R})/\sqrt{2}\), linearity of the constraint demands that each branch carry its own gravitational field. The physical state is therefore \[\begin{equation} \boxed{\; \ket{\Psi_{\mathrm{phys}}} = \frac{1}{\sqrt{2}} \bigl(\ket{L}\ket{\Phi_L} + \ket{R}\ket{\Phi_R}\bigr)\,, \;} \label{eq:entangled-state} \end{equation}\] where \(\ket{\Phi_A} = \hat{D}(\alpha_A)\ket{0}\) with \(\alpha_A(\mathbf{k})\) given by [eq:coherent-amplitude]. This state is necessarily entangled between matter and geometry.

The product state is unphysical.

The initial state assumed in the standard Feynman-Vernon calculation, \(\ket{\psi_{\mathrm{matter}}} \otimes \ket{0_{\mathrm{grav}}}\), violates the constraint [eq:poisson-constraint]. In the product state, the gravitational field is in the vacuum regardless of the matter configuration—there are no correlations between \(\hat{\Phi}\) and \(\hat{\rho}\), while the constraint demands perfect correlation. Starting from this product state amounts to asking: “At what rate does a bare, undressed mass become entangled with the graviton field?” This is a well-posed perturbative question, but it is not the physical question. The physical question is: “Given that the constraint has already entangled the mass with its gravitational field, at what rate does this pre-existing entanglement cause operational decoherence?”

Incorporating the constraint into the path integral

We now impose the constraint directly on the Schwinger-Keldysh path integral. In the ADM formalism, the Hamiltonian constraint is enforced by integrating over the lapse function \(N\), which acts as a Lagrange multiplier : \[\begin{equation} \int\!\mathcal{D}N\; \exp\!\Bigl(-\frac{i}{\hbar} \int\!dt\; N\,\hat{\mathcal{H}}_{\!\perp}\Bigr) = \delta\!\bigl[\hat{\mathcal{H}}_{\!\perp}\bigr]\,. \label{eq:lapse-delta} \end{equation}\] The constrained Schwinger-Keldysh functional integral for the reduced matter density matrix is obtained by inserting delta-function projectors onto the constraint surface : \[\begin{align} \rho_{\!M}(q_f,q'_f;t) &= \int\!dq_i\,dq'_i \int\!\mathcal{D}q^{+}\mathcal{D}q^{-} \int\!dh_f \int\!\mathcal{D}h^{+}\mathcal{D}h^{-} \nonumber\\[4pt] &\quad\times\; \delta\!\bigl[\mathcal{C}(h^{+},q^{+})\bigr]\; \delta\!\bigl[\mathcal{C}(h^{-},q^{-})\bigr] \nonumber\\[4pt] &\quad\times\; \exp\!\Bigl(\frac{i}{\hbar} \bigl[S[q^{+},h^{+}] - S[q^{-},h^{-}]\bigr]\Bigr) \nonumber\\[4pt] &\quad\times\; \delta(h^{+}_f - h_f)\,\delta(h^{-}_f - h_f)\; \rho_0(q_i,q'_i;h_i,h'_i)\,, \label{eq:constrained-PI} \end{align}\] where:

• \(q^{\pm}\) and \(h^{\pm}\) are the forward/backward matter and gravitational field paths, respectively;

• \(\mathcal{C}(h,q) \equiv \nabla^{2}\Phi - 4\pi G\,\rho_q\) is the constraint functional;

• \(\delta[\mathcal{C}(h^{+},q^{+})]\) and \(\delta[\mathcal{C}(h^{-},q^{-})]\) enforce the constraint independently on each branch of the closed-time-path contour;²

• \(\rho_0\) is the initial state—now the constrained, entangled state [eq:entangled-state], not a product state.

Solving the constraint: ADM decomposition

The linearized gravitational field decomposes into : \[\begin{equation} h_{\mu\nu} = \underbrace{\Phi}_{\text{Newtonian (scalar)}} \;+\; \underbrace{h^{\mathrm{TT}}_{ij}}_{\text{transverse-traceless}} \;+\; \text{(gauge modes)}\,, \label{eq:ADM-decomposition} \end{equation}\] where the three sectors play distinct physical roles under the constraint:

1. Newtonian sector (\(\Phi\)): The constraint \(\nabla^{2}\Phi = 4\pi G\rho_q\) completely determines \(\Phi\) from the matter configuration at each time slice. For boundary conditions \(\Phi \to 0\) at infinity, the solution is unique: \(\Phi^{\pm}(\mathbf{x},t) = \Phi_{\mathrm{cl}}[q^{\pm}(t)](\mathbf{x})\). The functional integral over \(\Phi\) collapses: this sector has no independent quantum fluctuations. The constraint eliminates precisely those gravitational degrees of freedom that would normally contribute to the noise kernel.

2. Transverse-traceless sector (\(h^{\mathrm{TT}}_{ij}\)): These two propagating polarizations (gravitational waves) satisfy \(\Box\,h^{\mathrm{TT}}_{ij} = (16\pi G/c^{4})\,T^{\mathrm{TT}}_{ij}\). For static masses in superposition, the transverse-traceless source vanishes: \(T^{\mathrm{TT}}_{ij} = 0\) (gravitational wave emission requires time-varying quadrupole moments). Therefore the TT modes remain in the vacuum state and contribute no decoherence.

3. Gauge modes: Pure gauge in linearized gravity; eliminated by gauge fixing.

The upshot is decisive: after solving the constraint, the gravitational field path integral disappears for the Newtonian sector. The potential on each branch is a deterministic functional of the matter path, and the TT sector decouples from static sources. The effective action for the matter becomes: \[\begin{equation} S_{\mathrm{eff}}[q] = S_M[q] + S_{\mathrm{grav\text{-}self}}[q]\,, \qquad S_{\mathrm{grav\text{-}self}}[q] = -\frac{G}{2}\int\!dt \int\!d^{3}x\,d^{3}y\; \frac{\rho_q(\mathbf{x})\,\rho_q(\mathbf{y})} {|\mathbf{x}-\mathbf{y}|}\,, \label{eq:effective-action} \end{equation}\] where the self-energy is a constant for a rigid body at fixed position (contributing only a phase).

Derivation of the constrained influence functional

With the constraint solved, the gravitational degrees of freedom reduce to the coherent states \(\ket{\Phi_A}\) attached to each matter branch. We now derive the influence functional by tracing over the final gravitational field configuration.

Step 1: The action difference gives a pure phase.

For the superposition with \(q^{+}(t) = \mathbf{x}_L\) and \(q^{-}(t) = \mathbf{x}_R\) (static paths), the action difference is \[\begin{equation} S_{\mathrm{eff}}[L] - S_{\mathrm{eff}}[R] = -E_G\,t\,, \label{eq:action-difference} \end{equation}\] where \[\begin{equation} E_G \equiv \frac{GM^{2}}{d} \label{eq:EG-def} \end{equation}\] is the gravitational self-energy difference between the two configurations. This action difference produces a pure phase factor \(e^{iE_G t/\hbar}\) in the off-diagonal density matrix element. A pure phase does not produce decoherence.

Step 2: The overlap factor.

Decoherence arises from a second factor: the trace over the final gravitational field state. In the constrained theory, the gravitational field in branch \(A\) is the coherent state \(\ket{\Phi_A(t)}\) determined by the constraint. Tracing over the field yields the overlap factor: \[\begin{equation} \mathcal{O}(L,R;t) = \int\!dh_f\; \braket{\Phi_L(t)|h_f}\!\braket{h_f|\Phi_R(t)} = \braket{\Phi_L(t)|\Phi_R(t)}\,. \label{eq:overlap-factor} \end{equation}\] This is the inner product of the two coherent states—the quantum mechanical distinguishability of the gravitational field configurations associated with the left and right branches.

Step 3: The constrained influence functional.

Combining the phase and the overlap, the constrained influence functional for the off-diagonal density matrix element is:

The reduced density matrix evolves as \[\begin{equation} \rho_{LR}(t) = \frac{1}{2}\, e^{iE_G t/\hbar}\, \braket{\Phi_L(t)|\Phi_R(t)}\,, \label{eq:rho-LR} \end{equation}\] and the decoherence is measured by the decay of \(|\rho_{LR}(t)|\), which is controlled entirely by the decoherence exponent: \[\begin{equation} \Gamma(t) = -\ln\bigl|\!\braket{\Phi_L(t)|\Phi_R(t)}\!\bigr|^{2}\,. \label{eq:decoherence-exponent} \end{equation}\]

\(G\)-counting: why the constrained result scales as \(G^1\)

The \(G\)-scaling of the constrained influence functional [eq:constrained-IF] differs fundamentally from that of the standard Feynman-Vernon result. We now trace the origin of this difference through explicit power counting.

Standard FV: \(G^2\) from the noise kernel.

In the standard (unconstrained) calculation, the influence functional takes the noise-kernel form (cf. Section 2): \[\begin{equation} \mathcal{F}_{\mathrm{FV}} = \exp\!\biggl( -\frac{1}{\hbar^{2}}\int_{0}^{t}\!dt'\!\int_{0}^{t}\!dt''\; (\Delta q)^{2}\,\mathcal{N}(t'-t'') \biggr)\,, \label{eq:FV-noise-recap} \end{equation}\] where \(\mathcal{N}\) is the symmetrized noise kernel (Hadamard function) of the gravitational field. Each interaction vertex contributes a factor of \(\sqrt{G}\) (from the matter-graviton coupling), and the noise kernel involves two such vertices, giving \((\sqrt{G})^{2} = G\). However, the graviton propagator \(\langle h\,h\rangle\) carries an additional factor of \(G\) (from the normalization of the graviton field: \(h_{\mu\nu} \sim \sqrt{G}\,\hat{a}\)), so the decoherence rate scales as \[\begin{equation} \Gamma_{\mathrm{FV}} \sim \frac{1}{\hbar^{2}} \times G \times G \times (\text{matter}) = \frac{G^{2}M^{4}}{\hbar^{3}d^{2}}\,. \label{eq:G2-counting} \end{equation}\] The two powers of \(G\) are unavoidable in the noise-kernel formalism: one from the coupling vertices, one from the propagator.

Constrained IF: \(G^1\) from the coherent-state overlap.

In the constrained case, the decoherence is controlled by the overlap \(|\!\braket{\Phi_L|\Phi_R}\!|\). For coherent states, the overlap formula gives (see Appendix 7 for details): \[\begin{equation} \bigl|\braket{\Phi_L|\Phi_R}\bigr|^{2} = \exp\!\bigl(-\|\delta\alpha\|^{2}\bigr)\,, \label{eq:coherent-overlap} \end{equation}\] where \[\begin{equation} \|\delta\alpha\|^{2} = \int\!\frac{d^{3}k}{(2\pi)^{3}}\; |\alpha_L(\mathbf{k}) - \alpha_R(\mathbf{k})|^{2}\,. \label{eq:delta-alpha-norm} \end{equation}\] From Eq. [eq:coherent-amplitude], the difference amplitude is \[\begin{equation} \delta\alpha(\mathbf{k}) \equiv \alpha_L(\mathbf{k}) - \alpha_R(\mathbf{k}) = -\frac{4\pi GM}{k^{2}} \cdot \frac{e^{-i\mathbf{k}\cdot\mathbf{x}_L} - e^{-i\mathbf{k}\cdot\mathbf{x}_R}} {\sqrt{2\hbar\omega_k}}\,. \label{eq:delta-alpha} \end{equation}\] The crucial \(G\)-counting is now transparent. The coherent-state amplitude \(\alpha \propto GM\), so \(\delta\alpha \propto GM\) and \(|\delta\alpha|^{2} \propto G^{2}M^{2}\). However, the overlap exponent \(\|\delta\alpha\|^{2}\) involves no graviton propagator—it is the norm of the amplitude in the single-particle Hilbert space, not a two-point correlation function. The mode integral gives \[\begin{equation} \|\delta\alpha\|^{2} = \frac{(4\pi GM)^{2}}{2\hbar} \int\!\frac{d^{3}k}{(2\pi)^{3}}\; \frac{2(1 - \cos\mathbf{k}\cdot\mathbf{d})} {k^{4}\,\omega_k} \;\propto\; \frac{GM^{2}}{\hbar\,c\,d}\,, \label{eq:norm-integral} \end{equation}\] where the last step uses \(\omega_k = ck\) and the IR cutoff provided by the separation \(d\) through the factor \((1 - \cos\mathbf{k}\cdot\mathbf{d})\). The integral evaluates to give an exponent proportional to \(GM^{2}/(\hbar c\,d)\)—one power of \(G\), not two.

The mechanism by which the second power of \(G\) is eliminated is clear: in the standard FV approach, one power of \(G\) comes from the graviton propagator \(\langle h\,h\rangle\) (which relates the quantum fluctuations of the field to the coupling). In the constrained approach, the gravitational field is not fluctuating independently—it is locked to the matter by the constraint. The relevant quantity is the distance between two coherent states in Hilbert space, not the amplitude of vacuum fluctuations. The propagator factor that contributed the second power of \(G\) is absent.

Comparison with the standard influence functional

Table 1 summarizes the structural differences between the standard and constrained influence functionals.

The difference between the two calculations is not a matter of approximation. Both are internally consistent within their respective frameworks. They differ in the question they answer:

• The standard FV calculation asks: starting from an undressed mass in the graviton vacuum, at what rate does dynamical graviton exchange generate entanglement between matter and field? The answer involves two interaction vertices (one emission, one absorption), hence two powers of the coupling \(\sqrt{G}\), giving \(\Gamma \propto G^{2}\).

• The constrained calculation asks: given that the Hamiltonian constraint has already entangled the mass with its gravitational field (a physical requirement, not an approximation), at what rate does this pre-existing entanglement produce operational decoherence? The answer involves the distinguishability of two coherent states, which depends on the norm \(\|\delta\alpha\|^{2}\)—a single power of \(G\) in the decoherence exponent.

Why gravity is special.

The reader may wonder why an analogous argument does not apply to QED, where Gauss’s law \(\nabla\!\cdot\!\hat{\mathbf{E}} = \hat{\rho}_e/\epsilon_0\) also constrains the field. At the level of the overlap computation (Sections 3.2–3.6), the two cases are indeed structurally similar: a charged particle in superposition is likewise dressed by branch-dependent coherent states of the longitudinal electric field, and the overlap \(\braket{E_L}{E_R}\) is less than unity.

The crucial difference emerges at the level of the rate extraction (Section 4.3). In QED, time is a background parameter and the Hamiltonian \(H_{\mathrm{QED}}\) is unconstrained: unitary evolution generated by \(H_{\mathrm{QED}}\) preserves the overlap \(|\!\braket{E_L(t)}{E_R(t)}\!|\), producing a reversible phase oscillation . In gravity, the Wheeler-DeWitt constraint \(\hat{H}_{\mathrm{total}} = 0\) eliminates background time. Physical time must emerge from internal correlations (the Page-Wootters mechanism ), and the irreducible quantum uncertainty of this gravitational clock converts the overlap reduction into irreversible decoherence. In the algebraic formulation, the Hamiltonian constraint converts the Type III observable algebra to Type II, introducing a finite trace and spectral gap that sets the decoherence rate . Gauss’s law restricts the state space but does not change the algebra type.

We emphasize that this distinction is fully operative only in the complete (nonlinear) quantum gravity theory. In the linearized limit used for the explicit computations of this paper, the Hamiltonian constraint reduces to the Poisson equation \(\nabla^2\Phi = 4\pi G\rho\), which has the same spatial form as Gauss’s law. The rate extraction in Section 4.3 therefore relies on invoking the full WDW structure—specifically, the modular Hamiltonian identification [eq:K-Hphys] and the Page-Wootters mechanism—within a linearized calculation that is consistent with (and controlled by) the full theory. A fully nonlinear derivation that makes the gravity/QED distinction manifest at each step remains an important open problem.

Validity.

The derivation presented here is valid within linearized gravity, where the expansion parameter is \(GM/(c^{2}d)\). For laboratory parameters (\(M = 1\;\mu\)g, \(d = 1\) mm), this ratio is \(\sim\!10^{-38}\). All higher-order corrections—graviton self-interaction, pair production, backreaction, and renormalization of \(G\)—are suppressed by at least \(O\bigl(G^{2}M^{2}/(\hbar c\,d)\bigr) \sim 10^{-38}\) relative to the leading term (see Appendix 8 for detailed estimates). The linearized approximation is controlled to extraordinary precision, and the result [eq:constrained-IF] is exact at \(O(G)\).

The Decoherence Rate

We now evaluate the constrained influence functional derived in Section 3 and extract the decoherence rate. The central result of this section is that the rate scales as \(G^1\), not \(G^2\), with the gravitational self-energy \(E_G = GM^2/d\) setting the scale.

Evaluating the coherent-state overlap

In the constrained framework, the decoherence factor is the overlap of two time-dependent coherent states of the gravitational field, one for each branch of the matter superposition. Before the superposition is created (\(t < 0\)), the mass is localized and the gravitational field is in a single coherent state \(\ket{\Phi_0}\) satisfying the constraint. At \(t = 0\), a beam splitter (or other device) places the mass in the superposition \((\ket{L} + \ket{R})/\sqrt{2}\). The constraint now demands two distinct coherent states \(\ket{\Phi_L}\) and \(\ket{\Phi_R}\), but the field cannot adjust instantaneously—the change propagates outward at the speed of light. The difference field between the two branches builds up causally from zero according to the driven oscillator solution (see Appendix 7 for details); the common-mode field \((\alpha_L + \alpha_R)/2\) remains equal to \(\alpha_0\) and drops out of the decoherence computation. The mode amplitude of the difference field at time \(t\) is \[\begin{equation} \alpha_A(\mathbf{k}, t) \;=\; \alpha_A^{\mathrm{eq}}(\mathbf{k})\, \bigl(1 - e^{-i\omega_k t}\bigr)\,, \label{eq:alpha-t} \end{equation}\] where \(\omega_k = c|\mathbf{k}|\) and the equilibrium amplitude is determined by the Newtonian potential of a point mass at position \(\mathbf{x}_A\): \[\begin{equation} \alpha_A^{\mathrm{eq}}(\mathbf{k}) \;=\; -\frac{4\pi GM}{k^2}\,\frac{e^{-i\mathbf{k}\cdot\mathbf{x}_A}} {\sqrt{2\hbar\omega_k}}\,. \label{eq:alpha-eq} \end{equation}\] At each instant the gravitational field state in branch \(A\) is the coherent state \(\ket{\Phi_A(t)} = D\!\bigl(\alpha_A(\cdot,t)\bigr)\ket{0}\), where \(D\) is the multimode displacement operator.

The squared overlap of the two branch states is \[\begin{equation} \bigl|\braket{\Phi_L(t)}{\Phi_R(t)}\bigr|^2 \;=\; \exp\!\Bigl(-\!\int\!\frac{d^3k}{(2\pi)^3}\, \bigl|\alpha_L(\mathbf{k},t) - \alpha_R(\mathbf{k},t)\bigr|^2\Bigr)\,, \label{eq:overlap-sq} \end{equation}\] a standard identity for coherent states. The difference amplitude factorizes: \[\begin{equation} \alpha_L(\mathbf{k},t) - \alpha_R(\mathbf{k},t) \;=\; \delta\alpha(\mathbf{k})\, \bigl(1 - e^{-i\omega_k t}\bigr)\,, \label{eq:delta-alpha-t} \end{equation}\] where the equilibrium difference \[\begin{equation} \delta\alpha(\mathbf{k}) \;=\; -\frac{4\pi GM}{k^2}\, \frac{e^{-i\mathbf{k}\cdot\mathbf{x}_L} - e^{-i\mathbf{k}\cdot\mathbf{x}_R}} {\sqrt{2\hbar\omega_k}} \label{eq:delta-alpha} \end{equation}\] encodes the spatial information of the superposition.

The decoherence exponent

Inserting Eqs. [eq:delta-alpha-t] and [eq:delta-alpha] into [eq:overlap-sq] and writing \(|1 - e^{-i\omega_k t}|^2 = 2(1-\cos\omega_k t)\), we obtain \[\begin{equation} \Gamma(t) \;\equiv\; -\ln\bigl|\braket{\Phi_L(t)}{\Phi_R(t)}\bigr|^2 \;=\; \int\!\frac{d^3k}{(2\pi)^3}\; \bigl|\delta\alpha(\mathbf{k})\bigr|^2\; 2\bigl(1 - \cos\omega_k t\bigr)\,. \label{eq:Gamma-exact} \end{equation}\] This is exact within linearized gravity with a free graviton field. The decoherence of the matter superposition is entirely determined by this single integral.

The squared difference amplitude is \[\begin{equation} \bigl|\delta\alpha(\mathbf{k})\bigr|^2 \;=\; \frac{(4\pi GM)^2}{k^4}\,\frac{1}{2\hbar\omega_k}\; 2\bigl(1 - \cos\mathbf{k}\cdot\mathbf{d}\bigr)\,, \label{eq:delta-alpha-sq} \end{equation}\] where \(\mathbf{d} = \mathbf{x}_L - \mathbf{x}_R\) is the separation vector. Substituting into [eq:Gamma-exact], converting to spherical coordinates, and performing the angular integral yields \[\begin{equation} \Gamma(t) \;=\; \frac{16G^2M^2}{\pi\hbar c} \int_0^{\Lambda}\!\frac{dk}{k^3}\; \Bigl(1 - \frac{\sin kd}{kd}\Bigr)\; \bigl(1 - \cos ckt\bigr)\,, \label{eq:Gamma-radial} \end{equation}\] where \(\Lambda\) is a UV cutoff (which will drop out of the rate). This integral contains two time regimes.

Short times (\(t \ll d/c\)).

All oscillatory factors can be expanded: \(1 - \cos ckt \approx (ckt)^2/2\). The exponent grows quadratically, \(\Gamma(t) \propto t^2\), giving Gaussian (non-Markovian) decay—the quantum Zeno regime.

Long times (\(t \gg d/c\)).

The modes with \(k \lesssim 1/d\) have completed many oscillations and contribute their time-averaged value \(\langle 2(1-\cos\omega_k t)\rangle \to 2\). In this regime \(\Gamma(t)\) approaches a constant—the static overlap \(\|\delta\alpha\|^2\) of the two equilibrium coherent states. The rate of approach to this equilibrium gives the decoherence rate.

The key physical quantity is the time derivative: \[\begin{equation} \frac{d\Gamma}{dt} \;=\; 2\int\!\frac{d^3k}{(2\pi)^3}\; \bigl|\delta\alpha(\mathbf{k})\bigr|^2\;\omega_k\,\sin\omega_k t\,. \label{eq:dGamma-dt} \end{equation}\] For \(t\) in the window \(d/c \ll t \ll t_{\mathrm{eq}}\) (after light-crossing but before full equilibration), this rate is effectively constant and equal to the Diósi-Penrose rate, as we now show through the constraint mechanism.

From energy scale to rate: the role of the Hamiltonian constraint

The free-field mode integral [eq:Gamma-radial] establishes that the decoherence energy scale is \(E_G = GM^2/d\), scaling as \(G^1\). However, treated as a free-field overlap, the decoherence exponent \(\Gamma(t)\) saturates at the equilibrium value \(2\|\delta\alpha\|^2\) rather than growing linearly in time (see Appendix 7). This is because each mode contributes a bounded oscillatory factor \(2(1 - \cos\omega_k t)\), and the continuous integral converges to a finite constant by the Riemann-Lebesgue lemma.

The extraction of a decoherence rate \(\Gamma = E_G/\hbar\) (linear growth in \(t\)) requires additional physical input beyond the free-field overlap. This input comes from the Hamiltonian constraint, which changes the physics in a qualitative way. The linearized Wheeler-DeWitt constraint \[\begin{equation} \bigl(\hat{H}_{\mathrm{matter}} + \hat{H}_{\mathrm{grav}}\bigr) \ket{\Psi_{\mathrm{phys}}} = 0 \label{eq:WDW-constraint} \end{equation}\] is not merely a condition on the initial state; it is enforced at all times. This has three interrelated consequences:

1. No independent graviton dynamics. The constraint continuously slaves the gravitational field to the matter configuration. The gravitational field does not propagate as an independent degree of freedom—its state is determined, mode by mode, by the matter distribution. This replaces the free graviton propagator (which costs one power of \(G\) in perturbation theory) with a constraint-determined classical field (which costs zero additional powers of \(G\)).

2. Physical time from the constraint. Since \(H_{\mathrm{total}} = 0\), the physical state is “timeless.” Physical time emerges relationally: the matter system evolves with respect to the gravitational field as an internal clock (the Page-Wootters mechanism ). The rate of decoherence is set by the energy gap between the two constraint-satisfying branches, which is the gravitational self-energy \(E_G = GM^2/d\).

3. Modular Hamiltonian identification. The Bisognano-Wichmann theorem  identifies the vacuum modular Hamiltonian \(K_0\) with \(2\pi\) times the boost generator. For coherent-state perturbations of the gravitational field, the Baker-Campbell-Hausdorff expansion terminates at linear order—there are no bi-local corrections at \(O(G)\)—so the full modular Hamiltonian satisfies  \[\begin{equation} K = 2\pi H_{\mathrm{phys}} + O(G^2) \label{eq:K-Hphys} \end{equation}\] as an operator equation on the physical Hilbert space. Decoherence proceeds at the modular frequency \(\omega_{\mathrm{mod}} = \Delta K / \hbar = 2\pi E_G/\hbar\), yielding \[\begin{equation} \Gamma \;=\; \frac{E_G}{\hbar} \;=\; \frac{GM^2}{\hbar d}\,. \label{eq:rate-EG} \end{equation}\]

The \(G\)-counting is summarized in Table 2.

Diagrammatic picture

The \(G\)-counting admits a simple diagrammatic interpretation (Fig. 1).

In the unconstrained (perturbative) calculation, decoherence arises from a graviton exchange loop: one interaction vertex on the left branch (\(\sqrt{G}\)), one on the right branch (\(\sqrt{G}\)), connected by a free graviton propagator (\(G^0\)). Squaring the amplitude gives the double-commutator structure of the noise kernel, scaling as \((\sqrt{G})^2 \times (\sqrt{G})^2 = G^2\).

In the constrained calculation, there is no graviton propagator. Instead, each branch carries a constraint-determined coherent state of the gravitational field. Decoherence is the overlap of these two coherent states—a single “constraint insertion” that contributes \(E_G \sim G^1\). Diagrammatically, the loop opens into a tree: the graviton line is not a propagator but a background field fixed by the constraint.

The \(O(1)\) coefficient

The \(G\)-counting fixes the decoherence rate up to an \(O(1)\) prefactor \(C\). Three independent estimates constrain this coefficient.

(i) Mode counting.

The coherent-state overlap computation (Appendix 7) evaluates [eq:Gamma-exact] in the Newtonian limit and yields \(C = 1\) when the self-energy divergences are renormalized by subtracting the single-branch contributions. The resulting rate matches the Diósi master equation , in which the decoherence kernel is \[\begin{equation} \mathcal{D}[\rho] \;=\; -\frac{G}{\hbar}\! \int\!\!d^3x\,d^3y\;\frac{[\hat{\rho}(\mathbf{x}),\,[\hat{\rho}(\mathbf{y}),\,\rho]]} {|\mathbf{x}-\mathbf{y}|}\,. \label{eq:Diosi-master} \end{equation}\] This Lindblad generator produces a decoherence rate for a point-mass superposition of exactly \(\Gamma = GM^2/(\hbar d)\).

(ii) Constraint identification.

The Diósi noise kernel \(G/|\mathbf{x}-\mathbf{y}|\) is not an independent postulate: it is the Green function of the Poisson constraint \(\nabla^2\Phi = 4\pi G\rho\). This identification connects the noise kernel directly to the Hamiltonian constraint, giving a first-principles origin for the Diósi master equation and confirming \(C = 1\).

(iii) Modular flow.

The identification [eq:K-Hphys] relates the modular frequency to the physical energy gap. Different choices for the modular-to-physical time conversion (depending on the observer’s trajectory and the precise form of the Rindler approximation) yield \(C\) values in the range \([1/2,\, 2]\). The most natural choice—a static observer at the location of the mass—gives \(C = 1\).

Combining these estimates, our best determination is \(C = 1\), with a systematic uncertainty spanning \(C \in [1/2,\, 2]\). We emphasize that this uncertainty is in the \(O(1)\) prefactor only; the \(G^1\) scaling is robust. The main result is:

For \(C = 1\) (matching the Diósi master equation), a particle of mass \(M = 1\;\mu\mathrm{g}\) in a superposition of separation \(d = 1\;\mathrm{mm}\) has \[\begin{equation} \tau_{\mathrm{dec}} \;=\; \frac{\hbar d}{GM^2} \;=\; \frac{(1.055 \times 10^{-34}\;\mathrm{J\,s})(10^{-3}\;\mathrm{m})} {(6.674\times 10^{-11}\;\mathrm{m^3\,kg^{-1}\,s^{-2}})(10^{-9}\;\mathrm{kg})^2} \;\approx\; 1.58\;\mathrm{ns}\,. \label{eq:tau-numeric} \end{equation}\] This prediction is \(10^{35}\) times shorter than the perturbative \(G^2\) estimate [eq:G2-rate] and lies within the sensitivity window of planned experiments (Section 5).

Validity of the derivation

The result [eq:main-result] rests on linearized gravity, where the gravitational field is treated as a free quantum field on a flat background sourced by a classical mass distribution. The expansion parameter is \(GM/(c^2 d) \sim 10^{-38}\) for the laboratory parameters above, so all \(O(G^2)\) corrections are suppressed by this factor. In particular:

• Graviton self-interactions: contribute at \(O(G^2)\) and are negligible.

• Graviton pair production: the squeezing parameter \(r_{\mathrm{sq}} \sim GM/(c^2 d) \sim 10^{-38}\) is far too small to produce appreciable non-coherent excitations.

• Backreaction: the gravitational field energy \(E_G \sim 10^{-20}\;\mathrm{J}\) is negligible compared to the rest mass energy \(Mc^2 \sim 10^{-2}\;\mathrm{J}\).

• Running of \(G\): renormalization-group corrections to \(G\) are suppressed by \((E/E_{\mathrm{Planck}})^2 \sim 10^{-76}\).

The linearized approximation is extraordinarily well controlled for all experimentally relevant parameter regimes. See Appendix 8 for a detailed analysis of each \(O(G^2)\) correction.

Experimental Predictions

The constrained influence functional derived in the preceding sections yields a decoherence rate \(\Gamma = C\,GM^2/(\hbar d) + O(G^2)\) with best estimate \(C = 1\), corresponding to a decoherence time \[\begin{equation} \tau_{\mathrm{dec}} = \frac{\hbar d}{GM^2}\,. \label{eq:tau-dec} \end{equation}\] For a particle of mass \(M = 1\;\mu\mathrm{g} = 10^{-9}\;\mathrm{kg}\) in a spatial superposition of separation \(d = 1\;\mathrm{mm} = 10^{-3}\;\mathrm{m}\): \[\begin{equation} \tau_{\mathrm{dec}} = \frac{1.055 \times 10^{-34} \times 10^{-3}}{6.674 \times 10^{-11} \times 10^{-18}} = 1.58\;\mathrm{ns}\,. \label{eq:tau-numerical} \end{equation}\] This is nine orders of magnitude shorter than the current experimental frontier for matter-wave coherence  and well within the timing resolution of existing detectors.

Predictions across mass scales

Table 3 compares the constrained (\(G^1\)) and perturbative (\(G^2\)) predictions across three experimentally relevant mass scales.

The separation between the two predictions is so large that even order-of-magnitude experimental sensitivity suffices to discriminate between them. No precision measurement of the coefficient \(C\) is needed—the gross scaling with \(G\) is the decisive test.

Distinctive experimental signatures

The \(G^1\) mechanism produces four distinctive signatures that, taken together, uniquely identify gravitational decoherence and distinguish it from all known environmental sources:

1. Mass scaling \(\boldsymbol{\Gamma \propto M^2}\). The decoherence rate increases as the square of the mass: doubling the mass quadruples the rate. This can be tested by comparing decoherence rates for particles of different masses under otherwise identical conditions. Environmental decoherence mechanisms scale differently—photon scattering scales with geometric cross-section (\(\propto M^{2/3}\) for constant-density particles), and collisional decoherence scales linearly with mass.

2. Temperature independence. The gravitational self-energy \(E_G = GM^2/d\) depends only on the mass configuration, not on the thermal state of the particle or its environment. Unlike thermal decoherence, which vanishes as \(T \to 0\), the rate [eq:G1-rate] persists at zero temperature. An experiment observing a decoherence floor that is independent of cryogenic cooling would provide strong evidence for a gravitational mechanism.

3. Vacuum independence. Decoherence from gas molecules, scattered photons, or blackbody radiation can be reduced by improving vacuum quality and electromagnetic shielding. Gravitational decoherence cannot be shielded and persists in the most perfect vacuum achievable. Residual decoherence that survives after all electromagnetic and collisional sources have been suppressed below the gravitational prediction would constitute evidence for the mechanism.

4. Linear separation scaling \(\boldsymbol{\Gamma \propto 1/d}\). The decoherence rate decreases as the inverse of the superposition separation: doubling the separation halves the rate. This is in contrast to the perturbative \(G^2\) prediction, where \(\Gamma \propto 1/d^2\), and to most environmental decoherence mechanisms, where larger superpositions decohere faster. The counterintuitive inverse scaling—a direct consequence of \(E_G = GM^2/d\)—provides a distinctive and testable prediction.

No known decoherence mechanism exhibits all four signatures simultaneously. Thermal photon scattering depends on temperature; collisional decoherence depends on vacuum quality; electromagnetic interactions can be shielded. The combination of \(M^2\) scaling, temperature independence, vacuum independence, and \(1/d\) separation scaling provides a unique fingerprint for gravitational decoherence.

Relevant experiments and timeline

Several experimental programs are approaching the regime where these predictions become testable:

• MAQRO (Macroscopic Quantum Resonators) . A proposed space-based mission designed to test the superposition principle for massive dielectric nanospheres (\(\sim 10^{9}\)–\(10^{10}\) amu) in the microgravity environment of a medium Earth orbit. The space environment eliminates seismic noise and enables long free-evolution times, directly targeting the picogram–nanogram mass range where \(G^1\) predictions give decoherence times of milliseconds to seconds.

• BECCAL (Bose-Einstein Condensate and Cold Atom Laboratory). The ISS-based facility for cold atom experiments in microgravity. While the primary science targets are atom interferometry and condensate physics, the platform provides the vibration-isolated, microgravity environment needed for massive-particle interferometry.

• Tabletop optomechanics . Ground-based experiments using levitated nanoparticles and optomechanical resonators have achieved quantum ground-state cooling and are pushing toward spatial superpositions of increasingly massive particles. Groups including those of Aspelmeyer (Vienna), Geraci (Northwestern), and Bouwmeester (Leiden) are developing techniques to create and detect superposition states of nanogram-scale particles.

• Matter-wave interferometry . Current experiments have demonstrated interference with molecules exceeding \(25\;\mathrm{kDa}\). Proposed next-generation interferometers (OTIMA, Talbot-Lau designs) aim at masses in the \(10^{6}\)–\(10^{9}\;\mathrm{amu}\) range, entering the regime where \(G^1\) scaling predicts observable decoherence.

The experimental timeline for definitive tests spans roughly 2028–2035. Ground-based optomechanical experiments may reach the nanogram regime within the next few years, enabling tests of the mass scaling. Space-based missions, if approved, would provide the controlled environment needed for the most precise measurements. A phased experimental program—first testing the \(M^2\) scaling law across accessible mass ranges, then measuring absolute decoherence times—can progressively discriminate between the \(G^1\) and \(G^2\) predictions.

The key experimental figure of merit is the ratio of the coherence time of the apparatus to the predicted gravitational decoherence time. When this ratio exceeds unity, the experiment becomes sensitive to the gravitational mechanism. For the \(G^1\) prediction, this threshold is reached at much lower masses and shorter evolution times than for \(G^2\), making the constrained prediction dramatically more accessible to near-term experiments.

Discussion and Conclusions

Summary of results

This paper establishes two main results, which we state separately to distinguish what is rigorously derived from what requires additional physical input.

Result 1: The \(G^1\) energy scale (rigorous within linearized gravity).

The standard Feynman-Vernon calculation  begins from the product state \(\ket{\psi_{\mathrm{matter}}} \otimes \ket{0_{\mathrm{grav}}}\), which violates the linearized Wheeler-DeWitt constraint \(\hat{H}_{\mathrm{matter}} + \hat{H}_{\mathrm{grav}} \approx 0\). Imposing the constraint forces the physical initial state into the entangled form \((\ket{L}\ket{\Phi_L} + \ket{R}\ket{\Phi_R})/\sqrt{2}\), where \(\ket{\Phi_A}\) are coherent states of the gravitational field. The decoherence exponent is controlled by the overlap \(\braket{\Phi_L(t)}{\Phi_R(t)}\), which depends on the gravitational self-energy \(E_G = GM^2/d\)—a single power of \(G\). This establishes \(E_G\) as the correct energy scale for gravitational decoherence, in contrast to the \(G^2\) energy scale that appears in the unconstrained noise-kernel mechanism.

The free-field overlap computation gives a decoherence exponent that saturates at the equilibrium value \(\Gamma_{\mathrm{sat}} \sim (GM^2/\hbar c)\ln(d/\ell_P)\) after the light-crossing time \(t \sim d/c\). For typical experimental parameters (\(M = 1\;\mu\)g, \(d = 1\;\)mm), this represents a modest but measurable coherence reduction.

Result 2: The \(G^1\) decoherence rate (requires the Hamiltonian constraint).

The extraction of a decoherence rate \(\Gamma = E_G/\hbar\) (linear growth of the decoherence exponent in time) goes beyond the free-field overlap and requires the full Wheeler-DeWitt constraint. The constraint eliminates independent graviton dynamics, forces physical time to emerge relationally (Page-Wootters mechanism ), and identifies the modular Hamiltonian with the physical Hamiltonian: \(K = 2\pi H_{\mathrm{phys}} + O(G^2)\) (via the Bisognano-Wichmann theorem  and the coherent-state BCH argument ). These three ingredients convert the \(G^1\) energy scale into a \(G^1\) rate: \[\begin{equation*} \Gamma = C \times \frac{GM^2}{\hbar d} + O(G^2)\,,\qquad C \in [1/2,\, 2]\,. \end{equation*}\]

This chain of reasoning is well-motivated within linearized gravity—the expansion parameter \(GM/(c^2 d) \sim 10^{-38}\) controls all approximations—but it invokes structural features of the full quantum gravity theory (the Hamiltonian constraint and its consequences for time) that go beyond what the Feynman-Vernon path integral alone can deliver. Experiment will ultimately determine whether the decoherence scales as \(G^1\) or \(G^2\).

Gravity versus electromagnetism

A natural objection is that electromagnetism has an analogous constraint—Gauss’s law, \(\nabla \cdot \mathbf{E} = \rho/\varepsilon_0\)—which likewise entangles a charged particle with its Coulomb field. At the level of the overlap computation (Section 3), the two cases are structurally similar: both produce branch-dependent coherent states and an overlap \(|\!\braket{\Phi_L}{\Phi_R}\!| < 1\) that scales as a single power of the coupling constant. The gravity/QED distinction lies not in the overlap itself but in its physical consequences—specifically, whether the overlap reduction leads to irreversible decoherence or merely a reversible phase.

In QED, time is a background parameter. Unitary evolution generated by \(H_{\mathrm{QED}}\) preserves the overlap magnitude: \(|\!\braket{E_L(t)}{E_R(t)}\!| = |\!\braket{E_L(0)}{E_R(0)}\!|\). The Coulomb-field entanglement produces a phase oscillation that can in principle be reversed .

In gravity, the Wheeler-DeWitt constraint \(\hat{H}_{\mathrm{total}} = 0\) eliminates background time. Physical time must emerge from internal correlations (Page-Wootters mechanism ), and the irreducible quantum uncertainty of the gravitational clock converts the overlap reduction into irreversible decoherence. In algebraic QFT, the Hamiltonian constraint converts the Type III observable algebra to Type II, introducing a finite trace and spectral gap . Gauss’s law restricts the state space but does not change the algebra type.

We acknowledge that this distinction is fully manifest only in the complete nonlinear theory. In the linearized limit used for explicit computations, the gravitational constraint reduces to the Poisson equation—the same spatial form as Gauss’s law. Our rate extraction (Section 4.3) therefore invokes structural features of the full theory within a linearized calculation. A fully nonlinear derivation that makes the gravity/QED distinction manifest at each step is an important open problem.

Relationship to prior work

This paper bridges the gap between the phenomenological Diósi-Penrose formula and the perturbative QFT result.

Diósi  postulated a stochastic gravitational noise field whose correlation kernel is the Newtonian potential \(1/|\mathbf{x} - \mathbf{y}|\), leading to a Lindblad master equation with decoherence rate \(GM^2/(\hbar d)\). Our derivation shows that this noise kernel is not an independent postulate: it follows from the Wheeler-DeWitt constraint. The constraint forces the gravitational field into branch-dependent coherent states whose overlap functional reproduces the Diósi kernel exactly in the Newtonian limit.

Penrose  argued on heuristic grounds that superpositions of distinct spacetime geometries should decay on a timescale \(\hbar/E_G\), invoking the incompatibility of time-translation operators in different geometries. Our result provides the microscopic mechanism behind this argument: the “incompatible time translations” correspond to distinct modular flows for the two coherent-state branches, and their mismatch—quantified by the coherent-state overlap—produces decoherence at precisely the rate Penrose predicted.

Anastopoulos and Hu  performed the most careful perturbative calculation of gravitational decoherence, obtaining the \(G^2\) rate. Their result is correct for the unconstrained product state. The present work does not invalidate their calculation; rather, it modifies the starting point. When the constraint is imposed on the initial state, the leading-order decoherence shifts from \(G^2\) (noise-kernel mechanism) to \(G^1\) (coherent-state-overlap mechanism). The noise-kernel contribution remains present at \(O(G^2)\) but is subleading by a factor of \(GM^2/(\hbar c\, d) \sim 10^{-38}\) for laboratory parameters.

Limitations and open questions

Several limitations should be acknowledged.

The \(O(1)\) coefficient. Our derivation constrains the prefactor \(C\) in \(\Gamma = C \times GM^2/(\hbar d)\) to the range \(C \in [1/2, 2]\), with the best estimate \(C = 1\) corresponding to the Diósi master equation. Fixing \(C\) exactly would require a complete treatment of the spin-2 form factor and the angular integration over the coherent-state mode spectrum. These geometric corrections are expected to be \(O(1)\) and do not affect the \(G\)-scaling, but they prevent us from quoting a coefficient to better than a factor of two.

Beyond linearized gravity. The entire derivation assumes linearized gravity. Corrections from graviton self-interaction, graviton pair production, gravitational backreaction, and running of Newton’s constant are all \(O(G^2)\)-suppressed. For laboratory masses, the expansion parameter \(GM/(c^2 d) \sim 10^{-38}\) ensures that these corrections are negligible to extraordinary precision, so the linearized approximation introduces no practical limitation. A fully nonlinear treatment would be desirable on theoretical grounds but is not required for comparison with experiment.

Rate extraction. The free-field coherent-state overlap gives a decoherence exponent that saturates at the equilibrium value \(\Gamma_{\mathrm{sat}} \sim (GM^2/\hbar c)\ln(d/\ell_P)\), rather than growing linearly in time. The extraction of a decoherence rate \(\Gamma = E_G/\hbar\) requires additional physical input: the operator-level identification \(K = 2\pi H_{\mathrm{phys}} + O(G^2)\) (from the Bisognano-Wichmann theorem and coherent-state BCH argument), combined with the Page-Wootters mechanism for emergent time. While each step in this chain is well-motivated within linearized gravity, the modular Hamiltonian \(K\) is defined with respect to a specific spacetime region, and the Page-Wootters mechanism invokes features of the full WDW equation not present in the linearized limit. This makes the rate extraction less rigorous than the energy-scale identification. The distinction between saturation and linear growth is experimentally testable: saturation predicts a one-time coherence reduction, while linear growth predicts exponential decay with time constant \(\tau = \hbar d/(GM^2)\).

Continuous mass distributions. We have worked with a two-branch superposition of a point mass, yielding \(E_G = GM^2/d\). For extended mass distributions, the gravitational self-energy generalizes to the Diósi kernel \(E_G = G \iint [\rho_L(\mathbf{x}) - \rho_R(\mathbf{x})][\rho_L(\mathbf{y}) - \rho_R(\mathbf{y})]\, d^3x\, d^3y / (2|\mathbf{x} - \mathbf{y}|)\). The coherent-state-overlap formalism extends naturally to this case—the displacement amplitudes \(\alpha_A(\mathbf{k})\) are determined by the Fourier transform of the mass density—but a complete computation for realistic experimental geometries has not yet been carried out.

Conclusions

This paper demonstrates that the Wheeler-DeWitt constraint changes the dominant mechanism for gravitational decoherence from the noise-kernel (\(G^2\)) to coherent-state overlap (\(G^1\)). The constrained influence functional rigorously establishes the gravitational self-energy \(E_G = GM^2/d\) as the decoherence energy scale. The further step—converting this energy scale to a decoherence rate \(\Gamma = E_G/\hbar\)—requires the modular Hamiltonian identification and the Page-Wootters mechanism for emergent time, which invoke structural features of the full quantum gravity theory within a linearized calculation.

The standard perturbative result  and the constrained result are both internally consistent; they differ in the physical question they answer. The former describes entanglement generation from an unphysical product state (\(G^2\)); the latter describes the manifestation of constraint-enforced entanglement (\(G^1\)).

The \(G^1\) and \(G^2\) predictions differ by a factor of \(\sim 10^{35}\) for microgram masses at millimeter separations. This enormous gap makes experimental discrimination decisive: next-generation matter-wave interferometry and optomechanical experiments will either observe decoherence consistent with \(G^1\) scaling, or place bounds that definitively falsify it. Either outcome would constitute a landmark result at the interface of quantum mechanics and gravity.

Appendices

Coherent State Overlap Computation

This appendix provides the detailed computation of the coherent state overlap \(\braket{\Phi_L(t)}{\Phi_R(t)}\) that governs the decoherence factor in the constrained influence functional.

Displacement operator and coherent state overlap

A coherent state of the gravitational field is obtained by applying the displacement operator to the vacuum: \[\begin{equation} \ket{\Phi_A} = D(\alpha_A)\ket{0}\,,\qquad D(\alpha) = \exp\!\Bigl(\int\!\frac{d^3k}{(2\pi)^3}\bigl[\alpha(\mathbf{k})\,\hat{a}^\dagger_{\mathbf{k}} - \alpha^*(\mathbf{k})\,\hat{a}_{\mathbf{k}}\bigr]\Bigr)\,, \label{eq:displacement-op} \end{equation}\] where \(A \in \{L,R\}\) labels the mass position. For two single-mode coherent states \(\ket{\alpha}\) and \(\ket{\beta}\), the overlap formula is \[\begin{equation} \braket{\alpha}{\beta} = \exp\!\Bigl(-\tfrac{1}{2}|\alpha|^2 - \tfrac{1}{2}|\beta|^2 + \alpha^*\beta\Bigr)\,, \end{equation}\] so that the modulus squared is \[\begin{equation} |\braket{\alpha}{\beta}|^2 = \exp\!\bigl(-|\alpha - \beta|^2\bigr)\,. \label{eq:single-mode-overlap} \end{equation}\] For a multi-mode field, the total overlap factorizes over independent modes, giving \[\begin{equation} \bigl|\braket{\Phi_L}{\Phi_R}\bigr|^2 = \exp\!\bigl(-\|\delta\alpha\|^2\bigr)\,,\qquad \|\delta\alpha\|^2 \equiv \int\!\frac{d^3k}{(2\pi)^3}\,\bigl|\alpha_L(\mathbf{k}) - \alpha_R(\mathbf{k})\bigr|^2\,. \label{eq:multi-mode-norm} \end{equation}\]

Mode amplitudes and the difference field

In linearized gravity, the coherent state amplitude sourced by a point mass \(M\) at position \(\mathbf{x}_A\) is  \[\begin{equation} \alpha_A(\mathbf{k}) = -\frac{4\pi G M}{k^2}\,\frac{e^{-i\mathbf{k}\cdot\mathbf{x}_A}}{\sqrt{2\hbar\omega_k}}\,, \label{eq:app-coherent-amplitude} \end{equation}\] where \(\omega_k = c|\mathbf{k}|\) for relativistic graviton modes. The difference amplitude between the two branches is \[\begin{equation} \delta\alpha(\mathbf{k}) \equiv \alpha_L(\mathbf{k}) - \alpha_R(\mathbf{k}) = -\frac{4\pi G M}{k^2}\,\frac{e^{-i\mathbf{k}\cdot\mathbf{x}_L} - e^{-i\mathbf{k}\cdot\mathbf{x}_R}}{\sqrt{2\hbar\omega_k}}\,. \label{eq:app-delta-alpha} \end{equation}\] Taking the modulus squared: \[\begin{equation} |\delta\alpha(\mathbf{k})|^2 = \frac{(4\pi G M)^2}{k^4 \cdot 2\hbar\omega_k}\,\bigl|e^{-i\mathbf{k}\cdot\mathbf{x}_L} - e^{-i\mathbf{k}\cdot\mathbf{x}_R}\bigr|^2 = \frac{(4\pi G M)^2}{2\hbar\omega_k\,k^4}\;\cdot\;2\bigl(1 - \cos(\mathbf{k}\cdot\mathbf{d})\bigr)\,, \label{eq:app-delta-alpha-sq} \end{equation}\] where \(\mathbf{d} = \mathbf{x}_L - \mathbf{x}_R\) is the separation vector with \(|\mathbf{d}| = d\).

Angular integration

Inserting Eq. [eq:app-delta-alpha-sq] into the norm [eq:multi-mode-norm] and passing to spherical coordinates \((k, \theta, \phi)\) with the polar axis along \(\mathbf{d}\): \[\begin{equation} \|\delta\alpha\|^2 = \int_0^{\Lambda}\!\frac{k^2\,dk}{2\pi^2}\;\frac{(4\pi GM)^2}{2\hbar ck\cdot k^4}\;\cdot\;2\!\int_0^1\!d\mu\;\bigl(1 - \cos(kd\mu)\bigr)\,, \end{equation}\] where \(\mu = \cos\theta\) and \(\Lambda\) is a UV cutoff. The angular integral evaluates to \[\begin{equation} \int_0^1 d\mu\;\bigl(1 - \cos(kd\mu)\bigr) = 1 - \frac{\sin(kd)}{kd}\,. \label{eq:angular-integral} \end{equation}\] For \(kd \gg 1\), this approaches unity; for \(kd \ll 1\), it behaves as \((kd)^2/6\). Collecting prefactors: \[\begin{equation} \|\delta\alpha\|^2 = \frac{16 G^2 M^2}{\pi\hbar c}\int_0^{\Lambda}\!\frac{dk}{k^3}\;\Bigl(1 - \frac{\sin(kd)}{kd}\Bigr)\,. \label{eq:norm-radial} \end{equation}\]

Evaluation of the radial integral

The integrand in Eq. [eq:norm-radial] has the structure \(f(kd)/k^3\) with \(f(u) = 1 - \sin u/u \sim u^2/6\) for small \(u\). The integral therefore converges in the infrared (\(k\to 0\)) but diverges logarithmically in the ultraviolet. Setting the IR scale at \(k_{\mathrm{IR}} \sim 1/d\) and the UV cutoff at \(\Lambda \sim 1/\varepsilon\) (where \(\varepsilon\) is the physical size of the mass distribution), the dominant contribution comes from modes with \(1/d \lesssim k \lesssim 1/\varepsilon\): \[\begin{equation} \int_{1/d}^{1/\varepsilon}\!\frac{dk}{k^3}\;\Bigl(1 - \frac{\sin(kd)}{kd}\Bigr) \;\approx\; \int_{1/d}^{1/\varepsilon}\!\frac{dk}{k^3}\;\cdot\;1 \;=\; \frac{1}{2}\Bigl(\frac{1}{(1/d)^2} - \frac{1}{(1/\varepsilon)^2}\Bigr) = \frac{d^2 - \varepsilon^2}{2}\,. \label{eq:radial-power} \end{equation}\] More precisely, a careful evaluation retaining the oscillatory piece yields a logarithmic correction. The full result for \(\varepsilon \ll d\) is \[\begin{equation} \|\delta\alpha\|^2 = \frac{8 G^2 M^2 d^2}{\pi\hbar c}\;\Bigl(1 + O\!\bigl(\varepsilon^2/d^2\bigr)\Bigr) \;+\; \text{subleading log corrections}\,. \end{equation}\]

Connection to the gravitational self-energy

The gravitational self-energy of the superposition is  \[\begin{equation} E_G = \frac{G}{2}\int\!\!\int\frac{\bigl[\rho_L(\mathbf{x}) - \rho_R(\mathbf{x})\bigr]\bigl[\rho_L(\mathbf{y}) - \rho_R(\mathbf{y})\bigr]}{|\mathbf{x} - \mathbf{y}|}\,d^3x\,d^3y\,. \label{eq:self-energy} \end{equation}\] For point masses \(\rho_A(\mathbf{x}) = M\delta^3(\mathbf{x} - \mathbf{x}_A)\), this reduces to \(E_G = GM^2/d\). The norm \(\|\delta\alpha\|^2\) is directly related to \(E_G\) through the identity \[\begin{equation} \|\delta\alpha\|^2 = \frac{2E_G}{\hbar c/d}\;\cdot\;\bigl(\text{form factor}\bigr)\,, \end{equation}\] where the form factor encodes the UV regularization and equals \(\ln(d/\varepsilon) + O(1)\) for an extended mass of size \(\varepsilon\). Physically, the logarithmic divergence for point masses reflects the well-known Newtonian self-energy divergence and is regulated by the finite size of the mass distribution.

Time-dependent overlap

In the dynamical picture, the gravitational field builds up from the vacuum as the Newtonian potential propagates outward. The time-dependent coherent state amplitude in branch \(A\) is \[\begin{equation} \alpha_A(\mathbf{k},t) = \alpha_A^{\mathrm{eq}}(\mathbf{k})\bigl(1 - e^{-i\omega_k t}\bigr)\,, \label{eq:app-time-dep-alpha} \end{equation}\] so the time-dependent norm becomes \[\begin{equation} \|\delta\alpha(t)\|^2 = \int\!\frac{d^3k}{(2\pi)^3}\;|\delta\alpha(\mathbf{k})|^2\;\cdot\;2\bigl(1 - \cos\omega_k t\bigr)\,. \label{eq:time-dep-norm} \end{equation}\] The factor \(2(1 - \cos\omega_k t)\) suppresses modes with \(\omega_k t \ll 1\) (i.e., wavelengths that have not yet had time to propagate). After the light-crossing time \(t \gg d/c\), the dominant modes (with \(k \sim 1/d\)) have equilibrated, and \(\|\delta\alpha(t)\|^2\) approaches the static equilibrium value \(2\|\delta\alpha_{\mathrm{eq}}\|^2\).

The free-field mode integral therefore gives a decoherence exponent that saturates at the equilibrium overlap: \[\begin{equation} -\ln\bigl|\braket{\Phi_L(t)}{\Phi_R(t)}\bigr|^2 \;\xrightarrow{t \gg d/c}\; 2\|\delta\alpha_{\mathrm{eq}}\|^2 \;\sim\; \frac{GM^2}{\hbar c}\ln\!\bigl(d/\ell_P\bigr)\,, \label{eq:saturation} \end{equation}\] which is finite and scales as \(G^1\). This establishes the gravitational self-energy \(E_G = GM^2/d\) as the correct energy scale for decoherence. The extraction of a rate \(\Gamma = E_G/\hbar = GM^2/(\hbar d)\)—i.e., linear-in-\(t\) growth rather than saturation—requires the additional physical input of the Hamiltonian constraint and modular Hamiltonian identification (Section 4.3), which converts the energy scale into a decoherence rate.

Robustness: \(O(G^2)\) Corrections

The derivation in the main text works within linearized gravity, where the gravitational field state is an exact coherent state. Here we verify that all corrections beyond this approximation are \(O(G^2)\) or higher, and are numerically negligible for laboratory parameters. Throughout, we use the reference values \(M = 1\;\mu\text{g} = 10^{-9}\;\text{kg}\) and \(d = 1\;\text{mm} = 10^{-3}\;\text{m}\).

Summary of corrections

We now discuss each effect in turn.

Graviton self-interaction. Beyond the free-field (quadratic) action, the Einstein-Hilbert action contains cubic and quartic vertices scaling as \(\sqrt{G}\) and \(G\), respectively. These vertices squeeze the coherent state, producing a correction \(\ket{\Phi_A} \to D(\alpha_A)S(\xi)\ket{0}\) where \(S(\xi)\) is the squeeze operator. The dimensionless squeezing parameter is \(|\xi| \sim GM/(c^2 d) \sim 7 \times 10^{-38}\) for the reference parameters. Squeezing modifies the non-local part of the modular Hamiltonian at \(O(|\xi|^2) \sim O(G^2)\) and therefore contributes to the decoherence rate only at \(O(G^2)\). No enhancement mechanism (infrared divergences, secular growth, or resonances) can promote this: the gravitational field is static, so there is no time-dependent driving to amplify the squeeze, and the infrared behavior is regulated by the finite separation \(d\).

Graviton pair production. A static mass superposition does not radiate gravitons (there is no time-dependent quadrupole moment within either branch). Virtual graviton pair production from vacuum fluctuations in the background of the superposition contributes to decoherence at \(O(G^2)\), suppressed by a factor \((\ell_P/d)^2 \approx (1.6 \times 10^{-35}/10^{-3})^2 \sim 10^{-64}\) relative to the leading \(G^1\) rate.

Backreaction on geometry. The stress-energy of the gravitational field itself, \(T^{\mathrm{grav}}_{00} \sim (\nabla\Phi)^2/(8\pi G)\), sources a correction to the metric at \(O(G^2)\). This modifies the coherent state amplitude at \(O(G^{3/2})\), which shifts the decoherence rate at \(O(G^2)\). The magnitude is \((GM/(c^2 d))^2 \sim 10^{-67}\).

Running of Newton’s constant. Quantum gravitational loop corrections renormalize Newton’s constant: \(G_{\mathrm{eff}}(E) = G(1 + c_1 GE^2/(\hbar c^5) + \cdots)\). For the gravitational self-energy scale \(E_G = GM^2/d \sim 10^{-13}\;\text{J}\), the correction is of order \((E_G/E_{\mathrm{Pl}})^2 \sim 10^{-69}\), where \(E_{\mathrm{Pl}} = \sqrt{\hbar c^5/G} \approx 1.96 \times 10^9\;\text{J}\).

Validity of the linearized approximation

The expansion parameter for linearized gravity is the dimensionless gravitational potential: \[\begin{equation} \epsilon \;\equiv\; \frac{GM}{c^2 d} \;\approx\; 7.4 \times 10^{-38}\,. \label{eq:epsilon-lin} \end{equation}\] This is satisfied by a factor of \(10^{38}\), placing the linearized approximation on extraordinarily firm ground. All post-Newtonian corrections enter at \(O(\epsilon^2)\) or higher.

Validity of the static approximation

The derivation assumes that the mass remains at rest during the decoherence time \(\tau_{\mathrm{dec}} = \hbar d/(GM^2) \approx 1.58\;\text{ns}\). The relevant comparison is with the mechanical timescale of the trapping potential. For a typical trap frequency \(\omega_{\mathrm{trap}} \sim 2\pi \times 100\;\text{Hz}\), the trap period is \(T_{\mathrm{trap}} \sim 10\;\text{ms}\), exceeding \(\tau_{\mathrm{dec}}\) by a factor of \(\sim 10^7\). The mass is also deeply non-relativistic: \(v/c \sim \sqrt{k_B T/(Mc^2)} \sim 10^{-15}\) at millikelvin temperatures. The static, non-relativistic approximation is therefore excellent.

Newtonian vs. relativistic propagation

The \(G^1\) rate derived in the main text uses the Newtonian (instantaneous) limit for the gravitational interaction. The fully relativistic computation (Appendix 7) shows that after the light-crossing time \(t_c = d/c \approx 3.3\;\text{ps}\), the Newtonian and relativistic results agree. Since \(t_c/\tau_{\mathrm{dec}} \sim 2 \times 10^{-3}\), the Newtonian approximation is accurate for all but the first few picoseconds of the decoherence process—a transient far shorter than any experimental time resolution.

← Back to Research