Introduction

The question of whether gravity is quantum has remained, despite a century of effort, without a direct experimental answer. Tests with matter-wave interferometers and levitated mirrors  have pushed the regime of demonstrated superposition into the mesoscopic but not yet into the regime where gravitational self-energy effects are accessible. The Bose–Marletto–Vedral (BMV) proposal  sharpened the question by identifying a tabletop test: place two massive particles each in spatial superposition, near enough that their mutual gravitational interaction can entangle their internal degrees of freedom, and detect that entanglement with a witness. By the LOCC argument, mediation of entanglement by the gravitational field implies that the field cannot be modelled as classical ; this conclusion has been further scrutinised by analysing the causality constraints that the same setup places on the gravitational field . The QGEM team  have developed a concrete experimental design using nanodiamonds with embedded NV centres in Stern–Gerlach interferometers, with target masses of \(\sim 10^{-17}\) kg and superposition arms of \(\sim 100\) \(\mu\)m.

A standard quantum-gravity treatment of the BMV setup predicts that the entanglement reaches the witness threshold after a time \[\begin{equation} \tau_{\mathrm{BMV}}\;\sim\;\frac{\hbar\, d^3}{G m^2 \Delta^2}\,, \label{eq:tau-BMV-intro} \end{equation}\] where \(d\) is the inter-particle separation and \(\Delta\) the superposition arm. Throughout this time, the gravitational field is treated as if it neither decoheres the source masses nor itself fluctuates appreciably – gravity is a passive mediator of phase.

Quantum-Geometric Duality (QGD) takes a different view of the gravitational field in this scenario. In the constrained Feynman–Vernon influence functional formulation , the linearised Wheeler–DeWitt constraint forces the field to remain in a dressed-coherent state for each matter configuration, and the result is a non-vanishing decoherence rate at first order in \(G\): \[\begin{equation} \Gamma\;=\;\frac{G m^2}{\hbar\,\Delta}\quad\text{per particle}\,, \label{eq:gamma-intro} \end{equation}\] the Diósi–Penrose rate  with \(d\to\Delta\) (the relevant length scale for a single mass in superposition is its own arm separation, not the inter-particle distance). A companion paper  establishes the same rate from the Hamiltonian-constraint perspective in linearised gravity.

The purpose of this paper is to combine these two ingredients – the standard BMV phase formula and the QGD per-particle decoherence rate – and ask: what does QGD predict for the visibility of the BMV witness?

The answer is sharp. At the BMV witness criterion (\(\phi_{\mathrm{BMV}} = \pi/2\), when the entanglement is half-formed), the off-diagonal element of the two-particle reduced density matrix is suppressed by a visibility factor \[\begin{equation} V\bigl(\tau_{\mathrm{BMV}}\bigr)\;=\;\exp\!\Bigl[\,-\,\mathcal{O}(1)\,(d/\Delta)^3\,\Bigr]\,. \label{eq:V-cube-intro} \end{equation}\] The cube law arises algebraically: substituting \(\Gamma\propto 1/\Delta\) into \(V = \exp(-2\Gamma\tau_{\mathrm{BMV}})\) with \(\tau_{\mathrm{BMV}}\propto d^3/\Delta^2\) yields the exponent \(\propto d^3/\Delta^3\). The mass dependence of \(\Gamma\) and \(\tau_{\mathrm{BMV}}\) cancels exactly, so the suppression depends only on the geometric ratio \(\Delta/d\). The numerical prefactor depends only on the arm geometry: \(\pi/2\) for parallel-axis arms and \(\pi\) for perpendicular (QGEM “diamond”) arms.

Three features make Eq. [eq:V-cube-intro] an unusually clean experimental discriminator. First, the mass-independence is unique among gravitational decoherence proposals: the Diósi–Penrose rate , continuous spontaneous localisation , and Penrose’s reduction-time hypothesis  all carry explicit mass scaling. Second, the prediction is parameter-free within QGD’s constrained influence functional – there is no analogue of the CSL collapse rate \(\lambda\) to fit. Third, the angular dependence of the leading BMV phase is a Legendre \(P_{2}(\cos\theta)\) in the arm-to-axis angle, with the parallel-axis (\(\theta=0\)) and perpendicular (\(\theta=\pi/2\)) geometries differing by a factor of two in the leading exponent and by the sign of the next-to-leading correction. Between these two cardinal orientations sits the magic angle \(\cos^{2}\theta_{m} = 1/3\) where the dipole term vanishes identically and a higher-multipole quintic law takes over. An experiment that can rotate its arms relative to the inter-particle axis therefore sweeps three qualitatively distinct QGD signatures with the same hardware: cubic suppression with prefactor \(\pi/2\) at \(\theta=0\), cubic suppression with prefactor \(\pi\) at \(\theta=\pi/2\), and quintic suppression with prefactor \(9\pi/7\) at \(\theta = \theta_{m}\).

Section 2 sets up the two-mass dynamics in QGD’s constrained-IF picture. Section 3 derives the leading-order cube law for both cardinal arm geometries. Section 4 computes the exact-phase corrections, identifies a sign-flip in the sub-leading term between them, and notes that the perpendicular formula has no singularity at \(\Delta\to d\) where the parallel formula breaks down. Section 5 extends the analysis to arbitrary arm angle \(\theta\), identifies the magic angle, and derives both the leading-order quintic visibility and its exact-phase form. Section 6 gives explicit numerical predictions for QGEM’s proposed parameters in all three geometries and compares with competing theories. Section 7 discusses limitations, falsification routes, and the place of this prediction among other near-term tests.

The Two-Mass Dynamics in QGD’s Constrained Influence Functional

State and field

Two equal masses \(m\), labelled \(A\) and \(B\), are each prepared in a balanced spatial superposition, \[\begin{equation} \ket{\psi_A}\;=\;\frac{1}{\sqrt 2}\bigl(\ket{L}_A + \ket{R}_A\bigr)\,, \qquad \ket{\psi_B}\;=\;\frac{1}{\sqrt 2}\bigl(\ket{L}_B + \ket{R}_B\bigr)\,, \end{equation}\] with arm separation \(\Delta\) and centre-to-centre separation \(d\) between the two particles. Two arm geometries are of experimental interest:

• Parallel-axis – the two arms of each particle lie along the inter-particle axis; the four branch distances are \(r_{LL}=r_{RR}=d\) and \(r_{LR}=d-\Delta\), \(r_{RL}=d+\Delta\) (so the minimum branch distance is \(d-\Delta\), requiring \(\Delta < d\)).

• Perpendicular-arm (QGEM “diamond”) – both arms of each particle lie in a plane transverse to the inter-particle axis, displaced by \(\pm\Delta/2\) from the central line; the branch distances are \(r_{LL}=r_{RR}=d\) and \(r_{LR}=r_{RL}=\sqrt{d^2+\Delta^2}\) (no constraint between \(\Delta\) and \(d\)).

The combined matter state has four computational-basis amplitudes, \[\begin{equation} \ket{\Psi^{\mathrm{matter}}}\;=\;\frac{1}{2}\sum_{i,j\in\{L,R\}}\ket{i}_A\ket{j}_B\,, \end{equation}\] each of which sources its own gravitational configuration. In the constrained Feynman–Vernon picture , integrating out the linearised gravitational degrees of freedom subject to the Wheeler–DeWitt constraint replaces this state with a dressed total state \[\begin{equation} \ket{\Psi(t)}\;=\;\frac{1}{2}\sum_{i,j} \ket{i}_A\ket{j}_B\,e^{-i\phi_{ij}(t)}\, \ket{\alpha_{ij}}\,, \label{eq:dressed-state} \end{equation}\] where \(\ket{\alpha_{ij}}\) is the constraint-enforced dressed coherent state of the gravitational field for matter configuration \((i,j)\) and the phase \(\phi_{ij}(t)\) is the Newtonian gravitational interaction phase for that configuration: \[\begin{equation} \phi_{ij}(t)\;=\;-\,\frac{E_{ij}\,t}{\hbar}\,, \qquad E_{ij}\;=\;-\,\frac{G m^2}{r_{ij}}\,. \label{eq:phi-ij} \end{equation}\]

Reduced matter density matrix

Tracing over the gravitational field, \[\begin{equation} \rho^{\mathrm{matter}}_{ij,kl}(t) \;=\;\frac{1}{4}\,\bra{\alpha_{kl}}\alpha_{ij}\rangle\, e^{-i(\phi_{ij}-\phi_{kl})\,t/\hbar}\,. \label{eq:rho-matter} \end{equation}\] Two distinct objects appear:

• the phase \(\phi_{ij}-\phi_{kl}\) is the standard BMV entanglement phase, generating coherent two-particle correlations as \(t\) grows;

• the dressed-state overlap \(\bra{\alpha_{kl}}\alpha_{ij}\rangle\) is the QGD-specific decoherence factor. In the conventional unconstrained Feynman–Vernon analysis this overlap equals unity (the gravitational field is pure noise that adds no off-diagonal damping at first order in \(G\)). In the constrained analysis it is a non-trivial Gaussian.

The QGD decoherence factor

Two coherent states satisfy \(|\bra{\alpha_a}\alpha_b\rangle| = \exp(-\|\delta\alpha_{ab}\|^2/2)\). In the constrained-IF picture the squared overlap exponent for two configurations differing by displacement of a single mass \(m\) over a distance \(\Delta\) saturates  at \[\begin{equation} \|\delta\alpha\|^2_{\mathrm{sat}} \;=\;\xi\,\ln(d/\ell_P)\,, \qquad \xi\;=\;\frac{G m^2}{\hbar c}\,, \label{eq:saturation} \end{equation}\] with \(\ell_P\) the Planck length. This logarithmic saturation describes the equilibrium of the constraint-dressed field after sound-crossing time \(d/c\). Beyond saturation, the constraint enforces a continued linear growth of the overlap exponent at the per-particle rate \[\begin{equation} \Gamma\;=\;\frac{G m^2}{\hbar\,\Delta} \label{eq:gamma-per-particle} \end{equation}\] for each particle whose configuration differs between the two branches. This linear growth is the post-saturation decoherence rate; it reduces to the Diósi–Penrose rate  for a single mass with the single-particle separation reinterpreted as the arm separation \(\Delta\) (the relevant length scale for one mass in superposition). The companion paper  derives the same rate from the linearised Hamiltonian-constraint perspective.

Throughout this work we are interested in the regime \(\tau \gg d/c\) where saturation has been reached and the linear-rate formula applies; for QGEM-class parameters \(d/c \sim 10^{-12}\) s while \(\tau_{\mathrm{BMV}} \sim 1\) s, so the saturated regime is overwhelmingly the relevant one.

The competition

The BMV proposal asks the experimenter to wait until the off-diagonal phase \(\phi_{ij}-\phi_{kl}\) has grown large enough to be witnessed. QGD adds a second clock: the dressed-state overlap \(|\bra{\alpha_{kl}}\alpha_{ij}\rangle|\) decays exponentially during this wait, at a rate set by \(2\Gamma\) for the channels in which both particles change configuration (e.g., \(LL \leftrightarrow RR\)). The remainder of the paper computes the visibility of the off-diagonal element at the time the BMV phase witness criterion is reached, in the saturated regime.

The Cube Law: Leading-Order Derivation for Both Geometries

Parallel-axis arms

For arms aligned with the inter-particle axis, the Newtonian interaction energies of the four configurations are \(E_{LL}=E_{RR}=-Gm^2/d\) and \(E_{LR}=-Gm^2/(d-\Delta)\), \(E_{RL}=-Gm^2/(d+\Delta)\). The off-diagonal phase between \(LL\) and \(RR\) is \(\phi_{LL}-\phi_{RR}=0\) (both diagonal channels share the same energy); the relevant phase for the witness is the one between the two “crossed” channels and the diagonal: \[\begin{equation} \Delta\phi_{\mathrm{par}}^{\mathrm{lead}}(t) \;=\;\frac{Gm^2 t}{\hbar} \,\biggl[\,\frac{2}{d}\,-\,\frac{1}{d-\Delta}\,-\,\frac{1}{d+\Delta}\,\biggr] \;=\;-\,\frac{Gm^2 t}{\hbar}\,\frac{2\Delta^2}{d(d^2-\Delta^2)}\,. \label{eq:phi-par-exact} \end{equation}\] The leading-order expansion in \(\Delta/d\) gives the standard BMV formula: \[\begin{equation} \Delta\phi_{\mathrm{par}}^{\mathrm{lead}}(t) \;\simeq\;\frac{2\,Gm^2\,\Delta^2\,t}{\hbar\,d^3}\,. \label{eq:phi-par-lead} \end{equation}\] The BMV witness criterion is reached when this phase equals \(\pi/2\): \[\begin{equation} \tau_{\mathrm{BMV}}^{\mathrm{par,lead}} \;=\;\frac{\pi\hbar\,d^3}{4\,Gm^2\,\Delta^2}\,. \label{eq:tau-par-lead} \end{equation}\]

The QGD visibility at this time, using the per-particle rate \(\Gamma=Gm^2/(\hbar\Delta)\) from Eq. [eq:gamma-per-particle], is the dressed-state-overlap modulus \[\begin{equation} V_{\mathrm{lead}}^{\mathrm{par}}\bigl(\tau_{\mathrm{BMV}}^{\mathrm{par,lead}}\bigr) \;=\;\exp\!\Bigl(-\,2\Gamma\,\tau_{\mathrm{BMV}}^{\mathrm{par,lead}}\Bigr) \;=\;\exp\!\Bigl(-\,\frac{\pi}{2}\,\bigl(d/\Delta\bigr)^3\Bigr)\,. \label{eq:V-par-lead} \end{equation}\] The factor of two in \(2\Gamma\) counts the two particles that each contribute coherence loss to the off-diagonal element. The mass dependence cancels: \(\Gamma\propto m^2\) and \(\tau_{\mathrm{BMV}}\propto 1/m^2\), so the exponent is geometric.

Perpendicular-arm (QGEM “diamond”) geometry

For arms transverse to the inter-particle axis, the four configurations have \(E_{LL}=E_{RR}=-Gm^2/d\) and \(E_{LR}=E_{RL}=-Gm^2/\sqrt{d^2+\Delta^2}\). The witness phase difference is \[\begin{equation} \Delta\phi_{\mathrm{perp}}(t) \;=\;\frac{2\,Gm^2 t}{\hbar} \,\biggl[\,\frac{1}{d}\,-\,\frac{1}{\sqrt{d^2+\Delta^2}}\,\biggr]\,. \label{eq:phi-perp-exact} \end{equation}\] At leading order in \(\Delta/d\), \[\begin{equation} \Delta\phi_{\mathrm{perp}}^{\mathrm{lead}}(t) \;\simeq\;\frac{Gm^2\,\Delta^2\,t}{\hbar\,d^3}\,, \label{eq:phi-perp-lead} \end{equation}\] which is exactly half the parallel-axis leading-order formula [eq:phi-par-lead]. The phase therefore takes twice as long to reach the \(\pi/2\) witness criterion: \[\begin{equation} \tau_{\mathrm{BMV}}^{\mathrm{perp,lead}} \;=\;\frac{\pi\hbar\,d^3}{2\,Gm^2\,\Delta^2} \;=\;2\,\tau_{\mathrm{BMV}}^{\mathrm{par,lead}}\,. \label{eq:tau-perp-lead} \end{equation}\] Substituting into the visibility formula gives \[\begin{equation} V_{\mathrm{lead}}^{\mathrm{perp}}\bigl(\tau_{\mathrm{BMV}}^{\mathrm{perp,lead}}\bigr) \;=\;\exp\!\Bigl(-\,\pi\,\bigl(d/\Delta\bigr)^3\Bigr)\,, \label{eq:V-perp-lead} \end{equation}\] twice the suppression coefficient of the parallel-axis case at the same \(\Delta/d\).

The cube law

Both leading-order results take the form \[\begin{equation} V_{\mathrm{lead}}^{(g)}\;=\;\exp\!\Bigl(-\,c_g\,(d/\Delta)^3\Bigr)\,, \qquad c_{\mathrm{par}}=\frac{\pi}{2}\,,\qquad c_{\mathrm{perp}}=\pi\,, \label{eq:cube-law} \end{equation}\] where the prefactor \(c_g\) depends only on the arm geometry \(g\). Observation that the suppression follows a cube law in \(d/\Delta\) would already be a strong qualitative signature of the underlying mechanism (per-particle rate \(\Gamma\propto 1/\Delta\) competing against \(\tau_{\mathrm{BMV}}\propto d^3/\Delta^2\)). Three features deserve emphasis:

1. Mass independence. The mass enters \(\Gamma\propto m^2\) and \(\tau_{\mathrm{BMV}}\propto 1/m^2\) symmetrically; the exponent in \(V\) is purely geometric. This contrasts sharply with the Diósi–Penrose collapse rate  (which scales as \(m^2\)), continuous spontaneous localisation  (which scales as the squared local mass density and thus is parameter-tuned), and Penrose’s reduction-time formula (also \(\propto m^2\)).

2. Parameter freedom. The prefactor \(c_g\) is fixed by the geometry alone; QGD has no analogue of the CSL collapse rate \(\lambda\) to fit. Up to the small ambiguity of choosing the witness criterion at \(\phi = \pi/2\) versus \(\phi = \pi\) (which doubles \(c_g\)), the prediction is rigid.

3. Geometry sensitivity. Choosing parallel versus perpendicular arms doubles the exponent at the same \(\Delta/d\). An experiment that can run with both arm orientations therefore provides an internal cross-check of the mechanism.

The next section examines the corrections that arise when the leading-order BMV phase is not adequate – specifically, when \(\Delta\) is comparable to \(d\) rather than much smaller, which is precisely the QGEM operating regime.

Exact-Phase Corrections and the Sign-Flip Between Geometries

The leading-order cube law of Section 3 assumes \(\Delta\ll d\). QGEM-class proposals operate at \(\Delta \sim d\) (sometimes \(\Delta > d\) in the parallel-axis layout, where the formula breaks down), so the sub-leading corrections matter for any quantitative comparison. The corrections arise from using the exact branch energies in Eqs. [eq:phi-par-exact] and [eq:phi-perp-exact] rather than their leading-order expansions. Both geometries develop sub-leading \((\Delta/d)\) corrections to the visibility exponent, but with opposite signs.

Parallel-axis exact-phase visibility

Setting the exact phase [eq:phi-par-exact] equal to \(\pi/2\) gives \[\begin{equation} \tau_{\mathrm{BMV}}^{\mathrm{par,exact}} \;=\;\frac{\pi\hbar\,d\,(d^2-\Delta^2)}{4\,Gm^2\,\Delta^2} \;=\;\tau_{\mathrm{BMV}}^{\mathrm{par,lead}}\,\bigl(1-(\Delta/d)^2\bigr)\,, \label{eq:tau-par-exact} \end{equation}\] i.e. the exact phase reaches the witness threshold faster than the leading-order formula predicts (because the singularity at \(\Delta\to d\) enhances the rate of phase accumulation). QGD therefore has less time to decohere the superposition, and the visibility prediction is correspondingly more permissive: \[\begin{equation} V_{\mathrm{exact}}^{\mathrm{par}}\bigl(d/\Delta\bigr) \;=\;\exp\!\biggl[\,-\,\frac{\pi}{2}\,\Bigl(\bigl(d/\Delta\bigr)^3-\bigl(d/\Delta\bigr)\Bigr)\biggr]\,. \label{eq:V-par-exact} \end{equation}\] The ratio to the leading-order prediction is \[\begin{equation} \frac{V_{\mathrm{exact}}^{\mathrm{par}}}{V_{\mathrm{lead}}^{\mathrm{par}}} \;=\;\exp\!\Bigl(\,+\frac{\pi}{2}\,(d/\Delta)\Bigr)\,, \label{eq:ratio-par} \end{equation}\] which at QGEM-relevant \(d/\Delta \sim 2\) is \(\exp(\pi)\approx 23\) – the leading-order formula over-suppresses the QGD visibility prediction by more than an order of magnitude in this regime.

Equation [eq:V-par-exact] reduces to the leading cube law in the wide-arm limit \(\Delta/d\to 0\) and diverges to \(V\to 1\) as \(\Delta\to d^-\), where the parallel-axis arms touch the partner particle and the linearised Newtonian formula breaks down. The expression is well-defined only for \(\Delta < d\).

Perpendicular-arm exact-phase visibility

Setting Eq. [eq:phi-perp-exact] equal to \(\pi/2\) and writing \(u\equiv\Delta/d\), one finds \[\begin{equation} \tau_{\mathrm{BMV}}^{\mathrm{perp,exact}} \;=\;\frac{\pi\hbar\,d}{4\,Gm^2}\,\frac{\sqrt{1+u^2}}{\sqrt{1+u^2}-1}\,. \label{eq:tau-perp-exact} \end{equation}\] After rationalising and substituting into the visibility formula (see Appendix 8), \[\begin{equation} V_{\mathrm{exact}}^{\mathrm{perp}}\bigl(u\bigr) \;=\;\exp\!\biggl[\,-\,\frac{\pi}{2}\,\frac{1+u^2+\sqrt{1+u^2}}{u^3}\,\biggr]\,. \label{eq:V-perp-exact} \end{equation}\] Two structural differences from the parallel-axis case:

• Sign of the sub-leading correction. Expanding Eq. [eq:V-perp-exact] in \(u\), \[\begin{equation} -\ln V_{\mathrm{exact}}^{\mathrm{perp}} \;=\;\frac{\pi}{u^3}\,+\,\frac{3\pi}{4}\,u^{-1}\,+\,\mathcal{O}(u)\,, \end{equation}\] which differs from the parallel-axis exact result \(-\ln V_{\mathrm{exact}}^{\mathrm{par}} = (\pi/2)\,u^{-3} - (\pi/2)\,u^{-1} + \cdots\) by the sign of the \(u^{-1}\) correction: parallel-axis is more permissive than its leading-order limit, perpendicular is more suppressive.

• No singularity at \(\Delta = d\). The perpendicular formula [eq:V-perp-exact] is regular for all positive \(u\), since the perpendicular arms never coincide with the partner particle. At \(u=1\) the suppression is \(V_{\mathrm{exact}}^{\mathrm{perp}}(1)=\exp\bigl(-\tfrac{\pi}{2}(2+\sqrt 2)\bigr)\approx 4.7\times 10^{-3}\), still well within the noise but far from the \(V\to 1\) blow-up of the parallel-axis formula.

Discrimination-feasible windows

Defining the experimentally interesting window as \(0.05 < V < 0.95\) (the regime where QGD can be distinguished from standard quantum gravity at modest signal-to-noise), the leading cube laws give \[\begin{equation} \begin{aligned} \text{parallel-axis (}V_{\mathrm{lead}}^{\mathrm{par}}\text{)}:\quad &\Delta/d \in [0.807,\,3.13]\,,\\ \text{perpendicular (}V_{\mathrm{lead}}^{\mathrm{perp}}\text{)}:\quad &\Delta/d \in [1.015,\,3.94]\,. \end{aligned} \label{eq:windows} \end{equation}\] The perpendicular window is shifted to larger \(\Delta/d\) by a factor \(\sim 1.26\) (the cube root of the prefactor ratio \(2\)). The parallel-axis window straddles the singularity at \(\Delta = d\); the realistic parallel-axis discrimination range is therefore \(\Delta/d \in [0.807,\,\sim 0.99]\) – a narrow but accessible regime. The perpendicular window is wider and entirely free of singularities.

Implication for QGEM design

At the QGEM design point \(\Delta = d/2\) (so \(u = 0.5\), \(d/\Delta = 2\)), the four visibility predictions are: \[\begin{equation} \begin{aligned} V_{\mathrm{lead}}^{\mathrm{par}}\bigl(2\bigr)\;&=\;\exp(-4\pi)\;\approx\;3.5\times 10^{-6}\,,\\ V_{\mathrm{exact}}^{\mathrm{par}}\bigl(2\bigr)\;&=\;\exp(-3\pi)\;\approx\;8.0\times 10^{-5}\,,\\ V_{\mathrm{lead}}^{\mathrm{perp}}\bigl(2\bigr)\;&=\;\exp(-8\pi)\;\approx\;1.3\times 10^{-11}\,,\\ V_{\mathrm{exact}}^{\mathrm{perp}}\bigl(0.5\bigr)\;&=\;\exp(-29.76)\;\approx\;1.2\times 10^{-13}\,. \end{aligned} \label{eq:V-QGEM-design} \end{equation}\] All four values are well below detection threshold: at the nominal QGEM operating point the QGD prediction is essentially zero residual visibility. To leave detectable residual visibility, the experiment must move to \(\Delta/d \to 1\) from below in the parallel-axis layout (within \(\sim 15\%\) of the singularity) or to \(\Delta/d \sim 2.5\) in the perpendicular layout.

Experimental design takeaway: the choice of arm orientation alters the QGD-vs-standard-QG discrimination by approximately five orders of magnitude in \(V\) at the QGEM design point. This sensitivity is not present in any standard quantum-gravity treatment of BMV and is therefore a directly testable prediction. The QGEM design literature should specify the intended arm orientation – QGD-vs-QG discrimination depends on it.

Oblique-Arm Geometry and the Magic Angle

Sections 3 and 4 treated the two limiting arm orientations \(\theta = 0\) (parallel-axis) and \(\theta = \pi/2\) (perpendicular). Any physical BMV apparatus can in principle orient its arms at an arbitrary angle \(\theta\) to the inter-particle axis, with each particle’s two arms displaced along the common unit vector \(\hat n = (\cos\theta, \sin\theta, 0)\). The full one-parameter family interpolating between Eqs. [eq:V-par-lead] and [eq:V-perp-lead] contains an isolated angle at which the leading-order BMV phase vanishes identically; the QGD prediction at that angle therefore probes a higher multipole of the source distribution and is qualitatively distinct from the cubic suppression at the two limits.

Branch distances and witness phase at arbitrary \(\theta\)

For arms tilted by \(\theta\) relative to the inter-particle axis, the four branch separations are \[\begin{equation} r_{LL} = r_{RR} = d\,, \qquad r_{\pm}(\theta) \;=\; \sqrt{d^{2} \pm 2\,d\,\Delta\,\cos\theta + \Delta^{2}}\,. \label{eq:r-angular} \end{equation}\] At \(\theta = 0\) these reduce to \(r_{\pm} = d\pm\Delta\) (the parallel-axis case of Section 3); at \(\theta = \pi/2\) they collapse to \(r_{+} = r_{-} = \sqrt{d^{2}+\Delta^{2}}\) (the perpendicular case). The witness phase magnitude is \[\begin{equation} |\Delta\phi(\theta)| \;=\;\frac{G m^{2}\tau}{\hbar}\, \biggl|\,\frac{2}{d} - \frac{1}{r_{+}(\theta)} - \frac{1}{r_{-}(\theta)}\,\biggr|\,, \label{eq:phi-angular} \end{equation}\] the natural one-parameter extension of Eqs. [eq:phi-par-exact] and [eq:phi-perp-exact].

Multipole expansion and the \(P_{2}(\cos\theta)\) structure

Writing \(u\equiv\Delta/d\) and expanding the two inverse branch distances to fourth order in \(u\) (the odd-in-\(u\) terms cancel between \(r_{+}\) and \(r_{-}\) because of the \(\pm\) displacement) yields \[\begin{equation} |\Delta\phi(\theta)| \;=\;\frac{G m^{2}\tau}{\hbar\,d}\, \Bigl[\, \bigl|\,1 - 3\cos^{2}\theta\,\bigr|\,u^{2} \,+\,\tfrac{1}{4}\,\bigl|\,3 - 30\cos^{2}\theta + 35\cos^{4}\theta\,\bigr|\,u^{4} \,+\,\mathcal{O}(u^{6})\, \Bigr]\,. \label{eq:phi-angular-expanded} \end{equation}\] The angular factor in the leading \(u^{2}\) term is \(-2 P_{2}(\cos\theta)\) where \(P_{2}\) is the second Legendre polynomial: the BMV witness phase is a Newtonian dipole-dipole interaction at large \(d/\Delta\). Equation [eq:phi-angular-expanded] reproduces \(|\Delta\phi| = (2 G m^{2}\tau/\hbar d)\,u^{2}\) at \(\theta = 0\) (parallel-axis leading-order, Eq. [eq:phi-par-lead]) and \(|\Delta\phi| = (G m^{2}\tau/\hbar d)\,u^{2}\) at \(\theta = \pi/2\) (perpendicular, Eq. [eq:phi-perp-lead]).

The magic angle: \(\cos^{2}\theta_{m} = 1/3\)

The leading \(u^{2}\) coefficient \(|1-3\cos^{2}\theta|\) in Eq. [eq:phi-angular-expanded] vanishes at the unique angle \[\begin{equation} \cos^{2}\theta_{m} = \tfrac{1}{3} \quad\Longleftrightarrow\quad \theta_{m} = \arccos\!\bigl(1/\sqrt{3}\bigr) \approx 54.7356^{\circ}\,. \label{eq:magic-angle} \end{equation}\] At \(\theta_{m}\) the dipole contribution to the BMV phase is identically zero and the next-order (quadrupole) term takes over. Its coefficient evaluates to \[\begin{equation} \tfrac{1}{4}\,\bigl|\,3 - 30/3 + 35/9\,\bigr| \;=\;\tfrac{1}{4}\cdot\tfrac{28}{9} \;=\;\tfrac{7}{9}\,, \end{equation}\] so \[\begin{equation} |\Delta\phi(\theta_{m})|_{\mathrm{lead}} \;=\;\frac{7}{9}\,\frac{G m^{2}\tau}{\hbar\,d}\,(\Delta/d)^{4}\,. \label{eq:phi-magic-lead} \end{equation}\] The witness criterion \(|\Delta\phi|=\pi/2\) now fixes \[\begin{equation} \tau_{\mathrm{BMV}}^{\mathrm{magic,lead}} \;=\;\frac{9\pi\hbar\,d^{5}}{14\,G m^{2}\,\Delta^{4}}\,, \end{equation}\] which is parametrically longer than the cubic-geometry witness times of Eqs. [eq:tau-par-lead] and [eq:tau-perp-lead] by a factor \(\propto (d/\Delta)^{2}\).

Combining with the per-particle QGD rate \(\Gamma = G m^{2}/(\hbar\Delta)\) from Eq. [eq:gamma-per-particle] and the two-particle exponent \(2\Gamma\tau\), the visibility at the magic angle is \[\begin{equation} V_{\mathrm{lead}}^{\mathrm{magic}}\;=\;\exp\!\Bigl[\,-\,\tfrac{9\pi}{7}\,(d/\Delta)^{5}\,\Bigr]\,, \label{eq:V-magic-lead} \end{equation}\] a quintic suppression law, qualitatively distinct from the cubic \((\Delta/d)^{3}\) scaling of the two cardinal geometries. No CSL or Diósi–Penrose model can reproduce this scaling because in those models the decoherence rate depends only on the local mass distribution and is blind to the arm orientation; the \((d/\Delta)^{5}\) law is a direct fingerprint of the geometric structure of the Newtonian phase that QGD competes against.

The general one-parameter family follows from Eq. [eq:phi-angular-expanded] at leading order: \[\begin{equation} V_{\mathrm{lead}}(\theta) \;=\;\exp\!\biggl[\,-\,\frac{\pi\,(d/\Delta)^{3}}{|1-3\cos^{2}\theta|}\,\biggr]\,, \label{eq:V-angular-lead} \end{equation}\] which reproduces \(V_{\mathrm{lead}}^{\mathrm{par}}\) at \(\theta = 0\) (denominator \(|1-3|=2\)), \(V_{\mathrm{lead}}^{\mathrm{perp}}\) at \(\theta = \pi/2\) (denominator \(|1-0|=1\)), and diverges as \(\theta\to\theta_{m}\) where the dipole vanishes and the leading-order approximation must be replaced by Eq. [eq:V-magic-lead]. The exponent of Eq. [eq:V-angular-lead] is strictly monotone in \(\theta\) on each of the two intervals \([0,\theta_{m})\) and \((\theta_{m},\pi/2]\): the magic angle is the worst-case orientation for residual visibility (best-case for QGD discrimination), with the cardinal geometries straddling it.

Exact-phase magic-angle visibility

Equation [eq:V-magic-lead] retains only the leading \(u^{4}\) term of the phase. The exact-phase calculation at \(\theta_{m}\) proceeds without expansion. Defining the dimensionless function \[\begin{equation} g_{m}(u)\;\equiv\;\Bigl|\,2 - \bigl(1 + (2/\sqrt 3)\,u + u^{2}\bigr)^{-1/2} - \bigl(1 - (2/\sqrt 3)\,u + u^{2}\bigr)^{-1/2}\,\Bigr|\,, \label{eq:gm-def} \end{equation}\] the exact magic-angle witness phase is \(|\Delta\phi(\theta_{m})| = (G m^{2}\tau/\hbar d)\,g_{m}(u)\). Expanding \(g_{m}(u)\) analytically gives the closed series (see Appendix 9) \[\begin{equation} g_{m}(u) \;=\;\frac{7}{9}\,u^{4}\,\Bigl[\,1 - \tfrac{4}{7}\,u^{2} + \mathcal{O}(u^{4})\,\Bigr]\,, \label{eq:gm-series} \end{equation}\] in which the \(u^{2}\) (dipole) coefficient is identically zero — the defining property of the magic angle — and the \(u^{4}\) coefficient is exactly \(7/9\), recovering Eq. [eq:phi-magic-lead]. The first non-trivial correction enters at \(u^{6}\) with relative coefficient \(-4/7\).

Setting \(|\Delta\phi(\theta_{m})| = \pi/2\) and combining with \(\Gamma = G m^{2}/(\hbar\Delta)\) as before, the exact-phase visibility exponent is \[\begin{equation} -\ln V_{\mathrm{exact}}^{\mathrm{magic}} \;=\;\frac{\pi}{u\,g_{m}(u)}\,, \qquad V_{\mathrm{exact}}^{\mathrm{magic}} \;=\;\exp\!\Bigl[\,-\,\pi/(u\,g_{m}(u))\,\Bigr]\,. \label{eq:V-magic-exact} \end{equation}\] Substituting Eq. [eq:gm-series] yields the analytic expansion \[\begin{equation} -\ln V_{\mathrm{exact}}^{\mathrm{magic}} \;=\;\frac{9\pi}{7}\,(d/\Delta)^{5} \,+\,\frac{36\pi}{49}\,(d/\Delta)^{3} \,+\,\mathcal{O}\!\bigl(d/\Delta\bigr)\,. \label{eq:V-magic-exact-expand} \end{equation}\] The leading term reproduces Eq. [eq:V-magic-lead]. The decisive qualitative feature is the sign of the sub-leading correction: the \(36\pi/49\) coefficient is positive, so \[\begin{equation} V_{\mathrm{exact}}^{\mathrm{magic}} < V_{\mathrm{lead}}^{\mathrm{magic}} \quad\text{for all}\ u > 0\,. \label{eq:magic-sign} \end{equation}\] The exact phase makes the magic-angle suppression stronger than the leading quintic. This is the same direction as the perpendicular geometry (Eq. [eq:V-perp-exact] versus Eq. [eq:V-perp-lead]) and the opposite direction to the parallel-axis geometry (Eq. [eq:V-par-exact] versus Eq. [eq:V-par-lead]), completing the exact-vs-leading trilogy across the three computed geometries.

Numerical regime and discrimination window

Representative values of the magic-angle exact-phase exponent, mass-independent: \[\begin{equation} \begin{aligned} -\ln V_{\mathrm{exact}}^{\mathrm{magic}}\bigl(u=0.5\bigr)\;&=\;153.4\quad(V\approx 2.4\times 10^{-67})\,,\\ -\ln V_{\mathrm{exact}}^{\mathrm{magic}}\bigl(u=1.0\bigr)\;&=\;8.99\quad(V\approx 1.24\times 10^{-4})\,,\\ -\ln V_{\mathrm{exact}}^{\mathrm{magic}}\bigl(u=2.0\bigr)\;&=\;1.54\quad(V\approx 0.21)\,. \end{aligned} \label{eq:magic-numerics} \end{equation}\] The corresponding leading-order quintic exponents are \(129.3\), \(4.04\), and \(0.126\); the exact-versus-leading gap is \(24.2\), \(4.95\), and \(1.41\) respectively, with the same sign at all three points (exact more suppressive than leading). The two-term cubic correction \((36\pi/49)(d/\Delta)^{3}\) of Eq. [eq:V-magic-exact-expand] reproduces the gap qualitatively (\(18.5\), \(2.31\), \(0.29\) at the same \(u\)) but underestimates it: higher-order corrections in \(d/\Delta\) contribute at the \(\sim 20\)–\(80\%\) level once \(u\lesssim 1\). The closed-form Eq. [eq:V-magic-exact] is therefore the only reliable evaluation in that regime.

The discrimination-feasible window \(0.05 < V < 0.95\), computed numerically from Eq. [eq:V-magic-exact], is \[\begin{equation} \text{magic-angle (exact): }\quad \Delta/d \in [1.465,\,31.6]\,. \label{eq:window-magic} \end{equation}\] The exact-phase perpendicular window is the comparable \([1.322,\,31.6]\), sharing the same upper edge \(\approx 31.6\) because both exponents scale as \(\pi/(2u)\) for \(u\to\infty\) and intersect the \(V=0.95\) threshold at the same point. The parallel exact-phase window \([0.664,\,0.984]\) is roughly two orders of magnitude narrower than either of the other two (constrained from above by the formula’s singularity at \(\Delta = d\), where the parallel-axis arms touch the partner particle). The magic-angle and perpendicular geometries therefore both retain QGD-distinguishable residual visibility across a span of arm separations roughly twenty times wider than the parallel-axis layout permits.

Experimental and theoretical takeaways

Equations [eq:V-angular-lead] and [eq:V-magic-exact] extend the QGD prediction of Section 3 from a binary choice between two cardinal arm orientations to a continuous one-parameter family with three qualitatively distinct regimes:

• Cubic regime, \(|\cos^{2}\theta - 1/3|\) bounded away from zero: visibility scales as \(\exp(-c(\theta)(d/\Delta)^{3})\) with \(c(\theta) = \pi/|1-3\cos^{2}\theta|\), interpolating smoothly between \(\pi/2\) (parallel) and \(\pi\) (perpendicular).

• Quintic regime, \(\theta = \theta_{m}\): visibility scales as \(\exp(-(9\pi/7)(d/\Delta)^{5})\), with the dipole BMV phase contribution exactly cancelled by the geometric condition \(\cos^{2}\theta = 1/3\) and the quadrupole taking over.

• Transition region, \(|\cos^{2}\theta - 1/3| \lesssim u^{2}\): the leading-order expansion breaks down and the exact formula \(V = \exp(-\pi/(u\,g(u,\theta)))\) must be used; the discrimination window widens monotonically as \(\theta\to\theta_{m}\).

For experimental design the implications are direct. An apparatus that can rotate its arms relative to the inter-particle axis can sweep three distinct QGD signatures in one platform. Demonstrating the \(\cos^{2}\theta_{m} = 1/3\) phase null is itself a non-trivial test of Newtonian gravity at mesoscopic separations: any theory whose phase response retains a non-zero dipole component at the magic angle (any modification that breaks \(P_{2}\) structure) would manifest as residual BMV signal where QGD and Newtonian gravity predict the quintic-suppressed null. A combined measurement at \(\theta = 0\), \(\theta = \pi/2\), and \(\theta = \theta_{m}\) tests three independent QGD predictions with the same hardware: cubic suppression with prefactor \(\pi/2\), cubic suppression with prefactor \(\pi\), and quintic suppression with prefactor \(9\pi/7\).

Mass independence is preserved at every angle. At each \(\theta\) the QGD rate \(\Gamma\propto m^{2}\) and the witness time \(\tau_{\mathrm{BMV}}(\theta)\propto 1/m^{2}\) cancel exactly in the exponent, so \(V(\theta)\) depends only on \(\Delta/d\) and \(\theta\) for any choice of arm orientation.

Numerical Predictions and Comparison with Competing Theories

Forecast for QGEM-class parameters

The QGEM proposal  considers two equal nanodiamond masses (\(m\sim 10^{-17}\) kg each) in Stern–Gerlach interferometers, with arm separations \(\Delta \sim 100\)–\(250\) \(\mu\)m and inter-particle separations \(d \sim 100\)–\(250\) \(\mu\)m. The QGD prediction depends only on the dimensionless ratio \(\Delta/d\) and on the arm-to-axis angle, so we tabulate the visibility forecasts as a function of \(\Delta/d\) in Table 1, separately for the three computed geometries. All entries are computed analytically from Eqs. [eq:V-par-lead], [eq:V-perp-lead], [eq:V-par-exact], [eq:V-perp-exact], [eq:V-magic-lead], and [eq:V-magic-exact]; the same numbers are returned by the companion code referenced in Section 6.4.

The QGEM nominal design \(\Delta \approx d\) is therefore in a regime where the parallel-axis QGD prediction breaks down formally and both the perpendicular and magic-angle predictions give strong but still observable suppression (\(V_{\mathrm{exact}}^{\mathrm{perp}}\sim 5\times 10^{-3}\), \(V_{\mathrm{exact}}^{\mathrm{magic}}\sim 10^{-4}\)). To work in the discrimination-feasible window the experiment can either:

• adopt the parallel-axis layout with \(\Delta/d \in [0.85,\,0.95]\) – a narrow regime within \(\sim 15\%\) of the formula singularity, where QGD predicts \(V_{\mathrm{exact}}^{\mathrm{par}} \in [0.49,\,0.84]\) (i.e. \(16\)–\(51\%\) visibility loss);

• adopt the perpendicular-axis layout with \(\Delta/d \in [1.32,\,31.6]\) from Eq. [eq:V-perp-exact], where QGD predicts \(V_{\mathrm{exact}}^{\mathrm{perp}}\) varying smoothly between \(0.05\) and \(0.95\) with a clear cube-law signature in the suppression as \(\Delta/d\) varies; or

• adopt the magic-angle layout with \(\Delta/d \in [1.46,\,31.6]\) from Eq. [eq:V-magic-exact], where QGD predicts \(V_{\mathrm{exact}}^{\mathrm{magic}}\) following the quintic law \(\propto \exp(-(9\pi/7)(d/\Delta)^{5} - (36\pi/49)(d/\Delta)^{3})\) – a different functional form whose suppression decays more slowly with \(\Delta/d\).

Comparison with competing theories

Five candidate predictions for the BMV/QGEM visibility, evaluated at the same parameters:

The QGD prediction is the only one in this list that is both mass-independent and parameter-free. The Diósi–Penrose rate is set by \(m^2\) and depends on a regularisation length scale \(\sigma_0\) associated with the mass smearing. CSL has the free collapse rate \(\lambda\) and the localisation length \(r_C\). The Pikovski et al. time-dilation decoherence depends on the composition (number of internal degrees of freedom) and ambient temperature. Only QGD predicts a sharp \((d/\Delta)^{3}\) cubic law in the cardinal arm orientations and a \((d/\Delta)^{5}\) quintic law at the magic angle, with no fitted parameters and no mass dependence. None of the collapse models has any analogue of the magic-angle quintic law: the decoherence rate in CSL/D-P/Pikovski depends on the local matter distribution and is therefore blind to the orientation of the superposition arms relative to the inter-particle axis.

Falsification routes

The experimental discrimination is binary at the QGEM design point:

• If QGEM observes the standard BMV entanglement amplitude (\(V \approx 1\)) at \(\Delta/d \sim 1\) in any of the three computed arm orientations, QGD’s constrained-IF picture is falsified.

• If QGEM observes a cube-law suppression of the visibility, with the geometry-dependent prefactor (parallel: \(\pi/2\), perpendicular: \(\pi\)) and no mass dependence across a parameter sweep, QGD is selected over Diósi–Penrose, CSL, and Pikovski-class models in a single experiment.

• If QGEM observes a quintic-law suppression at the magic angle \(\theta_{m} = \arccos(1/\sqrt 3)\) with prefactor \(9\pi/7\), this in addition rules out any decoherence mechanism whose rate is independent of the arm orientation relative to the inter-particle axis – a constraint that no collapse model can satisfy.

Failure modes that would not falsify QGD specifically: any visibility loss with explicit \(m^2\) scaling (consistent with any collapse model), or any temperature dependence (consistent with Pikovski et al.). Both can be controlled by the standard experimental sweeps that QGEM is already planning.

Companion code

The numerical values in Tables 1 and 2 are computed by an open-source Python script (released alongside this paper) implementing the six QGD visibility formulae (cardinal-cubic in both arm orientations, exact-phase corrections in both, magic-angle quintic, exact magic-angle) and the standard BMV phase formula. The unit-test suite covers: the cube-law scaling, mass-independence, the leading-to-exact crossover in both cardinal geometries, the \(P_{2}(\cos\theta)\) angular dependence, the magic-angle dipole-phase null, the quintic-law scaling with prefactor \(9\pi/7\), the exact magic-angle closed form, and the analytic-series accuracy below the numerical seam \(u = 10^{-2}\). All tests pass to machine precision.

Discussion and Conclusions

Summary of results

We have derived the QGD prediction for the BMV/QGEM gravitationally-induced-entanglement visibility, in the full one-parameter family of arm orientations parametrised by the angle \(\theta\) between the superposition axis and the inter-particle axis. In the two cardinal geometries (\(\theta = 0\) parallel, \(\theta = \pi/2\) perpendicular) and at the magic angle (\(\cos^{2}\theta_{m} = 1/3\)), the visibility formulae are \[\begin{equation} \begin{aligned} V_{\mathrm{lead}}^{\mathrm{par}}\;&=\;\exp\!\bigl[-\tfrac{\pi}{2}\bigl(d/\Delta\bigr)^3\bigr]\,,\\ V_{\mathrm{exact}}^{\mathrm{par}}\;&=\;\exp\!\bigl[-\tfrac{\pi}{2}\bigl((d/\Delta)^3-(d/\Delta)\bigr)\bigr]\,,\\ V_{\mathrm{lead}}^{\mathrm{perp}}\;&=\;\exp\!\bigl[-\pi\bigl(d/\Delta\bigr)^3\bigr]\,,\\ V_{\mathrm{exact}}^{\mathrm{perp}}\;&=\;\exp\!\Bigl[-\tfrac{\pi}{2}\bigl(1+u^2+\sqrt{1+u^2}\bigr)/u^3\Bigr]\,, \quad u\equiv\Delta/d\,,\\ V_{\mathrm{lead}}^{\mathrm{magic}}\;&=\;\exp\!\bigl[-\tfrac{9\pi}{7}\bigl(d/\Delta\bigr)^5\bigr]\,,\\ V_{\mathrm{exact}}^{\mathrm{magic}}\;&=\;\exp\!\bigl[-\pi/(u\,g_{m}(u))\bigr]\,,\quad g_{m}(u) \text{ from Eq.~\eqref{eq:gm-def}}\,, \end{aligned} \end{equation}\] together with the leading-order angular interpolation \(V_{\mathrm{lead}}(\theta) = \exp[-\pi(d/\Delta)^{3}/|1-3\cos^{2}\theta|]\) between the cardinal geometries (Eq. [eq:V-angular-lead]), all of them mass-independent and free of fitted parameters. The leading prefactor doubles between the cardinal cubic geometries; the sign of the next-to-leading correction flips between them; and at the magic angle the dipole BMV phase vanishes identically and the suppression becomes quintic with prefactor \(9\pi/7\), exact-phase strictly stronger than leading.

Origin of the cube law

Two physical scales control the BMV experiment. The phase scale is the gravitational interaction energy difference between branches, \(\Delta E \sim G m^2 \Delta^2 / d^3\) (for \(\Delta \ll d\) in either geometry, up to the geometry-dependent prefactor). The decoherence scale is the per-particle dressed-state-overlap rate \(\Gamma = G m^2/(\hbar\Delta)\) derived from the constrained Feynman–Vernon influence functional . The witness time is \(\tau \sim \hbar/\Delta E \sim \hbar d^3/(G m^2 \Delta^2)\). Their product \[\begin{equation} \Gamma\,\tau\;\sim\;\frac{G m^2}{\hbar \Delta}\cdot\frac{\hbar d^3}{G m^2 \Delta^2} \;=\;\bigl(d/\Delta\bigr)^3 \end{equation}\] gives the cube law and the cancellation of the mass. The cube power is therefore an algebraic consequence of the two timescales QGD identifies in the problem; it is not a free choice.

What is, and is not, derived

The cube-law and quintic-law forms of the visibility, the leading and exact-phase prefactors, the magic-angle identification, and the \(P_{2}(\cos\theta)\) angular interpolation are derived from QGD’s constrained Feynman–Vernon influence functional plus the standard Newtonian phase formula expanded in the Legendre generating function for the branch distances (Appendix 9). The dressed-coherent-state overlap structure of the gravitational field, the saturation of the overlap exponent at \(\xi\,\ln(d/\ell_P)\), and the post-saturation linear-rate \(\Gamma = G m^2/(\hbar\Delta)\) are inputs from the companion paper  on the constrained influence functional.

Two results are not derived but conjectured:

• Validity beyond linearised gravity. The derivation uses the Newtonian limit of linearised general relativity throughout. At QGEM-class parameters this is well within the regime of validity (graviton wavelength \(\gg\) particle separation), but extension to higher mass scales would require post-Newtonian corrections to both the BMV phase formula and the constrained-IF dressed-state overlap.

• Effect of imperfect coherent-state alignment. We have assumed that the gravitational field for each matter configuration is exactly in the constraint-enforced dressed coherent state. Departures from exact alignment (due to, e.g., decoherence of the gravitational field by external sources) would weaken the QGD prediction; the magnitude of this effect has not been estimated.

Limitations and caveats

Saturation regime: the linear-rate formula \(\Gamma\propto 1/\Delta\) applies for \(\tau > d/c\) (the gravitational sound-crossing time). At QGEM scales \(d/c \sim 10^{-12}\) s while \(\tau_{\mathrm{BMV}}\sim 1\) s, so the saturated regime dominates by twelve orders of magnitude. For interferometers with \(\tau \lesssim d/c\) the rate is suppressed and the prediction does not apply.

Visibility versus concurrence: we have computed the visibility of the off-diagonal element of the two-particle reduced density matrix, not the full two-qubit concurrence. For the BMV witness criterion, the off-diagonal modulus dominates the concurrence monotonically; the cube law therefore translates directly to a concurrence suppression. A full treatment of the four-by-four matter density matrix is straightforward and would refine the numerical prefactors by \(\mathcal{O}(1)\) but does not alter the cube law or the mass-independence.

Environmental decoherence: any non-gravitational decoherence (residual gas, blackbody radiation, internal degrees of freedom of the nanodiamond)  acts in parallel and adds to the QGD suppression. The QGEM background-decoherence budget  of \(\Gamma_{\mathrm{env}}\,\tau\lesssim 1\) is comparable to the QGD suppression at \(\Delta/d\sim 1\) in the parallel-axis layout, so distinguishing the gravitational signature from the background requires the parameter sweep across \(\Delta/d\) that the cube law predicts.

Relation to other master-equation derivations

Several master equations for gravitationally induced decoherence have been derived from quantising linearised gravity perturbatively. Anastopoulos and Hu  and Blencowe  obtain a Lindblad equation by tracing over the transverse-traceless graviton bath; the resulting rate is \(\mathcal{O}(G^2)\) in the matter–graviton coupling and depends on the temperature of the graviton bath. The QGD rate [eq:gamma-per-particle] is \(\mathcal{O}(G)\) in the same coupling because the constrained influence functional retains the Hamiltonian-constraint sector of the gravitational field, which the perturbative trace integrates out . At QGEM-class parameters the QGD rate exceeds the Anastopoulos–Hu/Blencowe rate by many orders of magnitude (the BMV cube-law prediction would be invisible in either of the perturbative treatments), so an observation of cube-law suppression at the QGEM design point would also discriminate between QGD and these perturbative master equations.

Place among other near-term tests of QGD

This prediction sits alongside two other QGD-specific experimental signatures, summarised in Table 3.

The single-mass tests  discriminate QGD from theories that predict no gravitational decoherence at all; they cannot distinguish QGD from collapse models that also predict mass-dependent decoherence. The BMV test is the only near-term proposal whose mass-independent cube law is incompatible with all extant collapse models simultaneously.

Conclusions

QGD predicts that gravitational decoherence of the source masses suppresses the BMV/QGEM entanglement visibility by a mass-independent factor whose form depends on the arm-to-axis angle \(\theta\). In the two cardinal cubic geometries the suppression is \(V \sim \exp\bigl(-(d/\Delta)^{3}\bigr)\) with prefactor \(\pi/2\) (parallel) or \(\pi\) (perpendicular); at the magic angle \(\cos^{2}\theta_{m} = 1/3\) the dipole BMV phase vanishes identically and the suppression becomes quintic \(V \sim \exp\bigl(-(9\pi/7)(d/\Delta)^{5}\bigr)\). Including the exact phase, the parallel-axis formula is more permissive than its leading-order limit while the perpendicular and magic-angle formulae are more suppressive; the magic-angle exact-phase exponent acquires a positive \((36\pi/49)(d/\Delta)^{3}\) cubic correction, completing the exact-versus-leading trilogy across the three computed geometries. At the QGEM design point \(\Delta = d/2\) all three predictions are essentially zero residual visibility; the discrimination-feasible regime requires either parallel-axis \(\Delta \sim 0.85d\)–\(0.95d\), perpendicular \(\Delta\in [1.3d, 31.6d]\), or magic-angle \(\Delta\in [1.5d, 31.6d]\). Observation of standard BMV entanglement at QGEM-class parameters falsifies QGD’s constrained-IF picture; observation of mass-independent cubic suppression in the cardinal geometries together with quintic suppression at the magic angle selects QGD over all extant collapse models in a single experiment, since no collapse mechanism depends on the orientation of the superposition arms relative to the inter-particle axis.

Appendix

Algebra of the Exact-Phase Visibility Formulae

This appendix records the algebraic steps from the exact branch energies of Section 2 to the exact-phase visibility formulae of Section 4. Throughout, write \(u\equiv\Delta/d\).

Parallel-axis arms: from energy to visibility

The four-branch interaction energies for the parallel-axis layout are \[\begin{equation} E_{LL} = E_{RR} = -\frac{G m^2}{d}\,, \qquad E_{LR} = -\frac{G m^2}{d - \Delta}\,, \qquad E_{RL} = -\frac{G m^2}{d + \Delta}\,. \end{equation}\] The witness phase difference at time \(\tau\) is \[\begin{equation} \Delta\phi_{\mathrm{par}}(\tau) = \frac{\tau}{\hbar}\bigl(E_{LL}+E_{RR}-E_{LR}-E_{RL}\bigr) = -\,\frac{G m^2 \tau}{\hbar}\Biggl[\frac{2}{d}-\frac{1}{d-\Delta}-\frac{1}{d+\Delta}\Biggr]\,. \end{equation}\] Algebraically, \[\begin{equation} \frac{2}{d}-\frac{1}{d-\Delta}-\frac{1}{d+\Delta} \;=\;\frac{2(d^2-\Delta^2) - d(d+\Delta) - d(d-\Delta)}{d(d^2-\Delta^2)} \;=\;-\,\frac{2\Delta^2}{d(d^2-\Delta^2)}\,, \end{equation}\] so \[\begin{equation} |\Delta\phi_{\mathrm{par}}(\tau)| \;=\;\frac{2\,G m^2\,\Delta^2\,\tau}{\hbar\,d(d^2-\Delta^2)}\,. \end{equation}\] Setting this equal to \(\pi/2\) (the BMV witness criterion) gives \[\begin{equation} \tau_{\mathrm{BMV}}^{\mathrm{par,exact}} \;=\;\frac{\pi\hbar\,d(d^2-\Delta^2)}{4\,G m^2\,\Delta^2} \;=\;\frac{\pi\hbar}{4 G m^2}\,\frac{d^3}{\Delta^2}\,(1-u^2)\,. \label{eq:tau-par-exact-app} \end{equation}\] The leading-order limit (\(u\to 0\)) gives \(\tau_{\mathrm{BMV}}^{\mathrm{par,lead}} = \pi\hbar d^3/(4 G m^2 \Delta^2)\) as in Eq. [eq:tau-par-lead].

The QGD visibility uses the per-particle rate \(\Gamma = G m^2/(\hbar\Delta)\) and the two-particle exponent \(2\Gamma\tau\): \[\begin{equation} \begin{aligned} -\ln V_{\mathrm{exact}}^{\mathrm{par}} &\;=\;2\Gamma\,\tau_{\mathrm{BMV}}^{\mathrm{par,exact}} \;=\;2\,\frac{G m^2}{\hbar\Delta}\,\cdot\,\frac{\pi\hbar\,d(d^2-\Delta^2)}{4\,G m^2\,\Delta^2}\\ &\;=\;\frac{\pi}{2}\,\frac{d(d^2-\Delta^2)}{\Delta^3} \;=\;\frac{\pi}{2}\,\Bigl[(d/\Delta)^3 - (d/\Delta)\Bigr]\,, \end{aligned} \end{equation}\] which is Eq. [eq:V-par-exact]. The leading cube-law term recovers Eq. [eq:V-par-lead]; the sub-leading \(-(d/\Delta)\) term carries a positive sign in the visibility (i.e. less suppression than the leading-order formula), because the exact phase reaches the witness threshold faster than its leading-order approximation predicts.

Perpendicular-arm geometry: from energy to visibility

The four-branch interaction energies for the perpendicular-arm layout are \[\begin{equation} E_{LL} = E_{RR} = -\frac{G m^2}{d}\,, \qquad E_{LR} = E_{RL} = -\frac{G m^2}{\sqrt{d^2+\Delta^2}}\,. \end{equation}\] The witness phase difference is \[\begin{equation} |\Delta\phi_{\mathrm{perp}}(\tau)| \;=\;\frac{2\,G m^2\,\tau}{\hbar}\,\Biggl[\frac{1}{d}-\frac{1}{\sqrt{d^2+\Delta^2}}\Biggr] \;=\;\frac{2\,G m^2\,\tau}{\hbar\,d}\,\Biggl[1-\frac{1}{\sqrt{1+u^2}}\Biggr]\,. \end{equation}\] Setting this equal to \(\pi/2\) gives \[\begin{equation} \tau_{\mathrm{BMV}}^{\mathrm{perp,exact}} \;=\;\frac{\pi\hbar\,d}{4 G m^2}\,\frac{1}{\,1-1/\sqrt{1+u^2}\,} \;=\;\frac{\pi\hbar\,d}{4 G m^2}\,\frac{\sqrt{1+u^2}}{\sqrt{1+u^2}-1}\,. \label{eq:tau-perp-exact-app} \end{equation}\] Rationalising the denominator, \[\begin{equation} \frac{\sqrt{1+u^2}}{\sqrt{1+u^2}-1} \;=\;\frac{\sqrt{1+u^2}\,(\sqrt{1+u^2}+1)}{(1+u^2)-1} \;=\;\frac{1+u^2+\sqrt{1+u^2}}{u^2}\,. \end{equation}\] Hence \[\begin{equation} \tau_{\mathrm{BMV}}^{\mathrm{perp,exact}} \;=\;\frac{\pi\hbar\,d}{4 G m^2}\,\frac{1+u^2+\sqrt{1+u^2}}{u^2}\,. \end{equation}\]

The two-particle visibility exponent is then \[\begin{equation} \begin{aligned} -\ln V_{\mathrm{exact}}^{\mathrm{perp}} &\;=\;2\Gamma\,\tau_{\mathrm{BMV}}^{\mathrm{perp,exact}} \;=\;2\,\frac{G m^2}{\hbar\Delta}\,\cdot\,\frac{\pi\hbar\,d}{4 G m^2}\,\frac{1+u^2+\sqrt{1+u^2}}{u^2}\\ &\;=\;\frac{\pi}{2}\,\frac{d}{\Delta}\,\frac{1+u^2+\sqrt{1+u^2}}{u^2} \;=\;\frac{\pi}{2}\,\frac{1+u^2+\sqrt{1+u^2}}{u^3}\,, \end{aligned} \end{equation}\] which is Eq. [eq:V-perp-exact]. In the wide-arm limit \(u\to 0\), \(1+u^2+\sqrt{1+u^2}\to 2+\mathcal{O}(u^2)\), so \(-\ln V_{\mathrm{exact}}^{\mathrm{perp}}\to \pi/u^3\), recovering the leading cube law \(V_{\mathrm{lead}}^{\mathrm{perp}} = \exp(-\pi(d/\Delta)^3)\) of Eq. [eq:V-perp-lead].

Sub-leading expansion

Expanding both exact-phase exponents in \(u\): \[\begin{equation} \begin{aligned} -\ln V_{\mathrm{exact}}^{\mathrm{par}} &\;=\;\frac{\pi}{2}\Bigl[u^{-3}\,-\,u^{-1}\Bigr] \;=\;\frac{\pi}{2}\,u^{-3}\,-\,\frac{\pi}{2}\,u^{-1}\,,\\ -\ln V_{\mathrm{exact}}^{\mathrm{perp}} &\;=\;\frac{\pi}{2}\,\frac{1+u^2+\sqrt{1+u^2}}{u^3}\\ &\;=\;\frac{\pi}{2}\,u^{-3}\bigl(2+\tfrac{3}{2}u^2+\mathcal{O}(u^4)\bigr) \;=\;\pi\,u^{-3}\,+\,\frac{3\pi}{4}\,u^{-1}\,+\,\mathcal{O}(u)\,. \end{aligned} \end{equation}\] The leading-order coefficients are \(\pi/2\) (parallel) and \(\pi\) (perpendicular), giving the factor-of-two prefactor difference of Eq. [eq:cube-law]. The sub-leading \(u^{-1}\) coefficients are \(-\pi/2\) (parallel) and \(+3\pi/4\) (perpendicular), with opposite signs: the parallel-axis exact formula is more permissive than its leading-order limit by \(\exp(\pi u/2)\), the perpendicular-axis exact formula is more suppressive by \(\exp(3\pi u/4)\).

Algebra of the Oblique-Arm and Magic-Angle Formulae

This appendix records the algebraic steps from the angular branch distances of Eq. [eq:r-angular] to the multipole expansion of Eq. [eq:phi-angular-expanded], and from there to the exact magic-angle series [eq:gm-series]. Throughout, \(u\equiv\Delta/d\) and the displacement vector is \(\hat n = (\cos\theta,\sin\theta,0)\).

Legendre generating function for the branch distances

The two non-trivial branch distances at general \(\theta\) are \(r_{\pm}(\theta) = d\sqrt{1 \pm 2 u\cos\theta + u^{2}}\). Their reciprocals admit the Legendre generating function expansion \[\begin{equation} \frac{1}{r_{\pm}(\theta)}\;=\;\frac{1}{d}\sum_{\ell\ge 0}\bigl(\mp 1\bigr)^{\ell}\,P_{\ell}(\cos\theta)\,u^{\ell}\,, \label{eq:legendre-generating} \end{equation}\] valid for \(u < 1\). Summing the two branches, the odd-\(\ell\) terms cancel, \[\begin{equation} \frac{1}{r_{+}(\theta)} + \frac{1}{r_{-}(\theta)} \;=\;\frac{2}{d}\,\sum_{\ell\ge 0}P_{2\ell}(\cos\theta)\,u^{2\ell}\,. \end{equation}\] The witness phase combination of Eq. [eq:phi-angular] is therefore \[\begin{equation} \frac{2}{d} - \frac{1}{r_{+}} - \frac{1}{r_{-}} \;=\;-\,\frac{2}{d}\,\sum_{\ell\ge 1}P_{2\ell}(\cos\theta)\,u^{2\ell}\,. \label{eq:phi-multipole} \end{equation}\] Inserting the standard Legendre values \[\begin{equation} P_{2}(\cos\theta) = \tfrac{1}{2}(3\cos^{2}\theta - 1)\,, \qquad P_{4}(\cos\theta) = \tfrac{1}{8}(35\cos^{4}\theta - 30\cos^{2}\theta + 3)\,, \end{equation}\] and taking the modulus gives, with \(|\Delta\phi| = (G m^{2}\tau/\hbar)\bigl|\,2/d - 1/r_{+} - 1/r_{-}\,\bigr|\), \[\begin{equation} \begin{aligned} |\Delta\phi(\theta)| &\;=\;\frac{G m^{2}\tau}{\hbar d}\,\Bigl[\,2\,|P_{2}|\,u^{2} + 2\,|P_{4}|\,u^{4} + \mathcal{O}(u^{6})\,\Bigr]\\ &\;=\;\frac{G m^{2}\tau}{\hbar d}\, \Bigl[\,|1-3\cos^{2}\theta|\,u^{2} + \tfrac{1}{4}|3 - 30\cos^{2}\theta + 35\cos^{4}\theta|\,u^{4} + \mathcal{O}(u^{6})\,\Bigr]\,, \end{aligned} \end{equation}\] recovering Eq. [eq:phi-angular-expanded]. The leading \(u^{2}\) angular factor is exactly \(-2\,P_{2}(\cos\theta)\) in modulus; its zero \(\cos^{2}\theta_{m} = 1/3\) defines the magic angle of Eq. [eq:magic-angle].

Magic-angle phase function \(g_{m}(u)\)

At \(\theta = \theta_{m}\), \(\cos\theta_{m} = 1/\sqrt 3\), so \(2\cos\theta_{m} = 2/\sqrt 3\). Define \[\begin{equation} x_{\pm}(u)\;\equiv\;\bigl(1 \pm \tfrac{2}{\sqrt 3}u + u^{2}\bigr)^{-1/2} \;=\;\sum_{\ell\ge 0}(\mp 1)^{\ell}\,P_{\ell}(1/\sqrt 3)\,u^{\ell}\,, \end{equation}\] so that \(g_{m}(u) = |2 - x_{+}(u) - x_{-}(u)|\) as in Eq. [eq:gm-def]. From the parity-summed series, \[\begin{equation} g_{m}(u) \;=\;\Bigl|\,-\,2\sum_{\ell\ge 1}P_{2\ell}(1/\sqrt 3)\,u^{2\ell}\,\Bigr|\,. \label{eq:gm-legendre} \end{equation}\] The relevant Legendre values are \[\begin{equation} \begin{aligned} P_{2}(1/\sqrt 3)\;&=\;\tfrac{1}{2}\,(3\cdot\tfrac{1}{3} - 1)\;=\;0\,,\\ P_{4}(1/\sqrt 3)\;&=\;\tfrac{1}{8}\,\Bigl(\tfrac{35}{9} - \tfrac{30}{3} + 3\Bigr)\;=\;\tfrac{1}{8}\cdot\tfrac{35 - 90 + 27}{9}\;=\;-\tfrac{7}{18}\,,\\ P_{6}(1/\sqrt 3)\;&=\;\tfrac{1}{16}\,\Bigl(\tfrac{231}{27} - \tfrac{315}{9} + \tfrac{105}{3} - 5\Bigr) \;=\;\tfrac{1}{16}\cdot\tfrac{231 - 945 + 945 - 135}{27}\;=\;\tfrac{2}{9}\,. \end{aligned} \label{eq:legendre-magic} \end{equation}\] The \(u^{2}\) coefficient of \(g_{m}\) vanishes (this is the magic-angle defining property \(P_{2}(1/\sqrt 3) = 0\)); the \(u^{4}\) coefficient is \(-2 P_{4}(1/\sqrt 3) = 7/9\); the \(u^{6}\) coefficient is \(-2 P_{6}(1/\sqrt 3) = -4/9\). Hence \[\begin{equation} g_{m}(u) \;=\;\tfrac{7}{9}\,u^{4}\,-\,\tfrac{4}{9}\,u^{6}\,+\,\mathcal{O}(u^{8}) \;=\;\tfrac{7}{9}\,u^{4}\,\Bigl[\,1 - \tfrac{4}{7}\,u^{2} + \mathcal{O}(u^{4})\,\Bigr]\,, \label{eq:gm-series-derived} \end{equation}\] which is Eq. [eq:gm-series]. The leading coefficient is exactly \(7/9\), reproducing Eq. [eq:phi-magic-lead] and hence \(V_{\mathrm{lead}}^{\mathrm{magic}} = \exp[-(9\pi/7)(d/\Delta)^{5}]\) of Eq. [eq:V-magic-lead]. The first correction enters at relative order \(u^{2}\) with negative coefficient \(-4/7\).

The choice of absolute value in Eq. [eq:gm-def] requires the two-term series to be sign-checked for the range \(u\in[0,1]\). Equation [eq:gm-series-derived] is positive in this range (\(1 - (4/7) u^{2} > 1 - 4/7 = 3/7 > 0\)), and direct numerical evaluation of \(|2 - x_{+}(u) - x_{-}(u)|\) confirms positivity at the few percent level for all \(u\in[0,1]\) checked. Numerically at \(u = 0.5\): the two-term truncation gives \(g_{m} \approx (7/9)(0.0625)(1 - 1/7) = 0.0417\); the exact closed form Eq. [eq:gm-def] gives \(g_{m}(0.5) = 0.04097\), agreeing to better than \(2\times 10^{-2}\) relative accuracy. At \(u = 0.1\) the agreement is sub-permille; at \(u = 0.01\) it is sub-microsecond — the series is the appropriate evaluation below the seam \(u_{\mathrm{seam}} = 10^{-2}\) where the closed form Eq. [eq:gm-def] loses thirteen significant digits to catastrophic cancellation \((x_{+} + x_{-}\to 2)\).

Exact-phase visibility expansion

Substituting Eq. [eq:gm-series-derived] into Eq. [eq:V-magic-exact], \[\begin{equation} -\ln V_{\mathrm{exact}}^{\mathrm{magic}} \;=\;\frac{\pi}{u\,g_{m}(u)} \;=\;\frac{\pi}{u}\cdot\frac{9}{7 u^{4}}\cdot\bigl[\,1 - \tfrac{4}{7}u^{2}\,\bigr]^{-1}\,. \end{equation}\] Expanding the bracket as a geometric series in \((4/7)u^{2}\) and keeping the first two terms, \[\begin{equation} -\ln V_{\mathrm{exact}}^{\mathrm{magic}} \;=\;\frac{9\pi}{7\,u^{5}}\,\Bigl[\,1 + \tfrac{4}{7}u^{2} + \mathcal{O}(u^{4})\,\Bigr] \;=\;\frac{9\pi}{7}\,(d/\Delta)^{5} + \frac{36\pi}{49}\,(d/\Delta)^{3} + \mathcal{O}(d/\Delta)\,, \end{equation}\] which is Eq. [eq:V-magic-exact-expand]. The sub-leading cubic correction has the positive sign \((+36\pi/49)\) and therefore strengthens the suppression relative to the pure-quintic leading term, as stated in Eq. [eq:magic-sign].

Cardinal-geometry consistency check

The Legendre generating-function form [eq:phi-multipole] reproduces the parallel-axis and perpendicular results of Sections 3 and 4 at \(\theta = 0\) and \(\theta = \pi/2\) respectively. At \(\theta = 0\), \(\cos\theta = 1\), and \(1/r_{\pm}(0) = (1/d)(1\mp u)^{-1} = (1/d)\sum_{\ell}(\mp 1)^{\ell} u^{\ell}\), recovering Eq. [eq:phi-par-exact]. At \(\theta = \pi/2\), \(\cos\theta = 0\) and only the \(u^{2\ell}\) terms with \(P_{2\ell}(0)\) non-zero contribute. Eq. [eq:legendre-generating] gives \(\sum_{\ell}P_{2\ell}(0)u^{2\ell} = (1+u^{2})^{-1/2}\) (the standard Taylor series at \(x = 0\)), so \[\begin{equation} \frac{1}{r_{+}(\pi/2)} + \frac{1}{r_{-}(\pi/2)} \;=\;\frac{2}{d}\,(1+u^{2})^{-1/2}\,, \end{equation}\] recovering Eq. [eq:phi-perp-exact]. Both cardinal geometries are therefore reproduced from the same Legendre identity, with the magic angle representing the unique zero of the leading multipole.

← Back to Research