State-Selective Signatures of Quantum and Classical Gravitational Environments

Abstract: A unified framework is developed for determining whether a gravitational-wave (GW) background behaves as a classical field or as a genuinely quantum environment. Unified here means that both descriptions originate from the same tidal coupling derived from geodesic deviation, which yields an identical quadratic interaction Hamiltonian for the detector; the only distinction lies in whether the GW degrees of freedom are modeled as classical phase-randomized coherent states or as quantized graviton modes. Within this common framework, the reduced dynamics of a quantum harmonic oscillator exhibit a sharp structural contrast: a quantized graviton bath preserves coherence within the lowest phonon-number manifold, forming a protected sector at leading order, whereas a classical stochastic GW field inevitably induces decoherence even inside this subspace. This difference provides an operational criterion for diagnosing the classical or quantum nature of gravitational waves using mesoscopic optomechanical systems. Our results establish decoherence structure - not merely its magnitude - as a sensitive probe of gravitational quantumness and delineate the experimental regimes under which such tests may become feasible.

• The paper develops a unified open-quantum-systems model showing that quadratic tidal coupling protects the {|0⟩,|1⟩} phonon subspace from vacuum-graviton decoherence through a Δn=±2 selection rule.
• It finds that thermal gravitons and phase-randomized classical backgrounds decohere all detector states, while deterministic classical waves produce oscillations without decoherence, making stochasticity essential to the classical signature.
• The proposed ratio R=Γ₀₂/(2Γ₀₁) equals 1 for generic linear Markovian noise but exceeds 1 for a quantum vacuum, while GHz HBAR measurements constrain graviton occupation to roughly n̄c≲8×10³³.

Overview and motivation

The paper develops a unified open-quantum-systems framework for distinguishing whether a gravitational-wave (GW) background behaves as a genuinely quantum environment or as a classical stochastic field. The unification is structural: both descriptions derive from the same tidal coupling obtained via geodesic deviation in the proper detector frame, yielding an identical quadratic interaction Hamiltonian between a single mechanical mode of a mesoscopic resonator and the transverse–traceless (TT) metric perturbation. The only difference lies in the statistical state assigned to the gravitational degrees of freedom: quantized graviton modes (vacuum or thermal Gibbs states) versus an ensemble of phase-randomized coherent states representing classical stochastic radiation.

The central claim is that the structure of decoherence—not its magnitude—provides an operational diagnostic of gravitational quantumness. Specifically, for a quantized graviton bath in the vacuum state, the reduced dynamics of the detector exhibits a protected Fock subspace within the lowest phonon manifold $\{|0\rangle,|1\rangle\}$ at leading order, whereas classical phase-randomized backgrounds and thermal graviton states induce decoherence even inside this subspace. This is a stronger statement than rate-based arguments: it is rooted in selection rules imposed by the quadratic detector–gravity coupling rather than in absolute decay rates.

Detector–gravity coupling

Starting from linearized gravity on Minkowski space in TT gauge, the authors work in Fermi normal coordinates along a freely falling detector worldline. In the long-wavelength approximation, geodesic deviation gives $\ddot{\xi}^i = \tfrac{1}{2}\ddot{h}^{TT}_{ij}\xi^j$, leading to the effective Hamiltonian

$$H_D = \frac{p_1^2}{2m_0} + \dot{h}^{TT}_{11}(t)\,\xi_1 p_1 + \frac{1}{2}m_0\omega_c^2\xi_1^2,$$

first introduced by Speliotopoulos and later used in quantum-gravity detector models. For a GW propagating along $z$ with plus polarization, this reduces to a time-dependent harmonic oscillator Hamiltonian with a bilinear $\xi_1 p_1 + p_1\xi_1$ coupling proportional to $\dot{\chi}(t)$.

A key quantitative motivation is scale separation. For micron-scale GHz devices (e.g., HBARs with $\omega_c/2\pi \simeq 5.96$ GHz), the direct GW-induced displacement $\Delta x_{\rm GW} = hL$ is roughly $10^{-29}$ m for $h\sim 10^{-23}$, giving $\ddot{\xi}^i = \tfrac{1}{2}\ddot{h}^{TT}_{ij}\xi^j$0. Displacement sensing is therefore hopeless; the relevant signature must be statistical—encoded in how the gravitational environment modifies coherence structure. LIGO-scale mirrors, by contrast, operate at high thermal occupation where any quantum signature is suppressed.

Quantum graviton environment

Quantizing the TT modes yields a total Hamiltonian with free graviton and detector terms plus an interaction $\ddot{\xi}^i = \tfrac{1}{2}\ddot{h}^{TT}_{ij}\xi^j$1, i.e., a two-phonon (quadratic) coupling with selection rule $\ddot{\xi}^i = \tfrac{1}{2}\ddot{h}^{TT}_{ij}\xi^j$2. Under Born–Markov and rotating-wave approximations, the reduced master equation is Lindblad-type with dissipators built from $\ddot{\xi}^i = \tfrac{1}{2}\ddot{h}^{TT}_{ij}\xi^j$3 and $\ddot{\xi}^i = \tfrac{1}{2}\ddot{h}^{TT}_{ij}\xi^j$4, weighted by $\ddot{\xi}^i = \tfrac{1}{2}\ddot{h}^{TT}_{ij}\xi^j$5 and $\ddot{\xi}^i = \tfrac{1}{2}\ddot{h}^{TT}_{ij}\xi^j$6 respectively, where $\ddot{\xi}^i = \tfrac{1}{2}\ddot{h}^{TT}_{ij}\xi^j$7 is the thermal occupation at the detector frequency.

Solving perturbatively to first order in the dissipative Liouvillian, the paper finds:

Initial state	Short-time evolution	Vacuum ($\ddot{\xi}^i = \tfrac{1}{2}\ddot{h}^{TT}_{ij}\xi^j$8) behavior
$\ddot{\xi}^i = \tfrac{1}{2}\ddot{h}^{TT}_{ij}\xi^j$9	$$H_D = \frac{p_1^2}{2m_0} + \dot{h}^{TT}_{11}(t)\,\xi_1 p_1 + \frac{1}{2}m_0\omega_c^2\xi_1^2,$$0	protected
$$H_D = \frac{p_1^2}{2m_0} + \dot{h}^{TT}_{11}(t)\,\xi_1 p_1 + \frac{1}{2}m_0\omega_c^2\xi_1^2,$$1	$$H_D = \frac{p_1^2}{2m_0} + \dot{h}^{TT}_{11}(t)\,\xi_1 p_1 + \frac{1}{2}m_0\omega_c^2\xi_1^2,$$2	protected
$$H_D = \frac{p_1^2}{2m_0} + \dot{h}^{TT}_{11}(t)\,\xi_1 p_1 + \frac{1}{2}m_0\omega_c^2\xi_1^2,$$3	coherences decay at $$H_D = \frac{p_1^2}{2m_0} + \dot{h}^{TT}_{11}(t)\,\xi_1 p_1 + \frac{1}{2}m_0\omega_c^2\xi_1^2,$$4	fully coherent
$$H_D = \frac{p_1^2}{2m_0} + \dot{h}^{TT}_{11}(t)\,\xi_1 p_1 + \frac{1}{2}m_0\omega_c^2\xi_1^2,$$5	$$H_D = \frac{p_1^2}{2m_0} + \dot{h}^{TT}_{11}(t)\,\xi_1 p_1 + \frac{1}{2}m_0\omega_c^2\xi_1^2,$$6	decoheres
$$H_D = \frac{p_1^2}{2m_0} + \dot{h}^{TT}_{11}(t)\,\xi_1 p_1 + \frac{1}{2}m_0\omega_c^2\xi_1^2,$$7	off-diagonal decays even at $$H_D = \frac{p_1^2}{2m_0} + \dot{h}^{TT}_{11}(t)\,\xi_1 p_1 + \frac{1}{2}m_0\omega_c^2\xi_1^2,$$8	decoheres

with $$H_D = \frac{p_1^2}{2m_0} + \dot{h}^{TT}_{11}(t)\,\xi_1 p_1 + \frac{1}{2}m_0\omega_c^2\xi_1^2,$$9. The decisive result is that vacuum fluctuations are blind to superpositions confined to $z$0 but not to those involving levels separated by two quanta—the latter couple via spontaneous emission of real on-shell gravitons through the $z$1 channel. Population redistribution from $z$2 to $z$3 at zero temperature is shown to preserve trace and complete positivity; the growth of ground-state population reflects unidirectional emission, not coherence restoration.

Thermal graviton states lift the protection: they interpolate between vacuum and classical Gaussian statistics, and can be physically motivated by pre-inflationary radiation eras or stimulated production, though the paper is careful to note "thermal" refers only to the density-matrix form, not dynamical graviton thermalization.

Classical stochastic limit

The classical background is modeled as a tensor product of phase-randomized coherent states, $z$4. A bookkeeping device separating $z$5 (graviton sector) from $z$6 (detector sector) tracks the classical limit $z$7 constant as $z$8, $z$9.

In this limit, decoherence occurs irrespective of detector state, including within $\xi_1 p_1 + p_1\xi_1$0, at rate $\xi_1 p_1 + p_1\xi_1$1. Crucially, Appendix D shows that a pure coherent state with fixed phase produces purely oscillatory $\xi_1 p_1 + p_1\xi_1$2-dependent contributions that vanish under secular averaging—so deterministic classical GWs induce no decoherence at all. Stochasticity, not classicality per se, is the essential ingredient behind classical gravitational decoherence. The paper also notes a technical subtlety: taking $\xi_1 p_1 + p_1\xi_1$3 prematurely destroys the oscillatory structure needed for the cancellation; the classical limit must be applied after obtaining the full reduced dynamics.

Decoherence timescales and spectral dependence

The general scaling is $\xi_1 p_1 + p_1\xi_1$4, since $\xi_1 p_1 + p_1\xi_1$5. For a power-law spectral density $\xi_1 p_1 + p_1\xi_1$6, the rate scales as $\xi_1 p_1 + p_1\xi_1$7; free gravitons in $\xi_1 p_1 + p_1\xi_1$8 dimensions give $\xi_1 p_1 + p_1\xi_1$9, so higher-frequency detectors experience stronger gravitational decoherence. The paper emphasizes that no absolute decoherence time can be predicted without specifying the environmental spectral profile—a point reinforcing the shift from magnitude-based to structure-based diagnostics.

State-selective experimental protocol

The proposed protocol exploits the ratio

$\dot{\chi}(t)$0

measured via $\dot{\chi}(t)$1/$\dot{\chi}(t)$2-type population and Ramsey measurements on superpositions $\dot{\chi}(t)$3 and $\dot{\chi}(t)$4. Generic Markovian environments with Lindblad operators linear in oscillator variables obey $\dot{\chi}(t)$5, hence $\dot{\chi}(t)$6 and $\dot{\chi}(t)$7—the same universal scaling holds for thermal graviton states and phase-randomized coherent backgrounds. Only the quantum vacuum yields $\dot{\chi}(t)$8, specifically $\dot{\chi}(t)$9, because $\omega_c/2\pi \simeq 5.96$0 while $\omega_c/2\pi \simeq 5.96$1. A statistically significant deviation of $\omega_c/2\pi \simeq 5.96$2 from unity therefore cannot be reproduced by any non-vacuum configuration derived from the same microscopic coupling.

Quantitatively, using HBAR parameters ($\omega_c/2\pi \simeq 5.96$3 GHz), the gravitational contribution to lowest-manifold decoherence evaluates to $\omega_c/2\pi \simeq 5.96$4. Demanding no excess over measured dephasing rates $\omega_c/2\pi \simeq 5.96$5 yields the bound

$\omega_c/2\pi \simeq 5.96$6

i.e., present GHz mechanical experiments already constrain non-vacuum graviton occupation at their operating frequency to below $\omega_c/2\pi \simeq 5.96$7. Although numerically large, this bound is model-independent within the perturbative regime and directly reflects Planck suppression ($\omega_c/2\pi \simeq 5.96$8). The paper notes that standard slow-roll inflation does not populate the GHz band (redshifted tensor modes cut off near $\omega_c/2\pi \simeq 5.96$9 Hz), so any nonzero occupation there would require alternative mechanisms such as preheating, phase transitions, or primordial black hole binaries.

Self-consistency of the Born–Markov factorization is verified in Appendix E: spontaneous graviton emission proceeds at rate $\Delta x_{\rm GW} = hL$0, so the induced bath occupation $\Delta x_{\rm GW} = hL$1 over all experimentally relevant times.

Limitations and open questions

Several caveats bear directly on the strength of the results. First, the entire analysis is perturbative (Born–Markov, first order in the dissipative Liouvillian); the protected-subspace statement is strictly a leading-order result, and the discrimination between $\Delta x_{\rm GW} = hL$2 and $\Delta x_{\rm GW} = hL$3 relies on the assumption that non-gravitational environmental channels are Markovian and linear in oscillator variables—an assumption that may fail for structured or nonlinear baths. Second, Planck suppression makes direct observation of vacuum-induced two-quantum processes unrealistic in the near term; the practically accessible output is the occupation bound $\Delta x_{\rm GW} = hL$4 and structural bounds on $\Delta x_{\rm GW} = hL$5, not detection of vacuum gravitons. Third, the treatment neglects squeezed graviton states, which time-dependent backgrounds generically produce; incorporating anomalous environmental correlations could modify the master-equation structure and is left as an open extension. Finally, the framework constrains only the occupation number at the detector frequency and assumes no specific cosmological origin for the background.

Conclusion

This work reframes gravitational decoherence diagnostics around the qualitative structure of reduced dynamics rather than decay-rate magnitudes. Within a single microscopic tidal coupling, the quantum vacuum preserves coherence in the $\Delta x_{\rm GW} = hL$6 manifold via a $\Delta x_{\rm GW} = hL$7 selection rule, while thermal, phase-randomized, and other non-vacuum backgrounds destroy this protection. The dimensionless selectivity ratio $\Delta x_{\rm GW} = hL$8 provides a falsifiable, noise-robust observable, and existing GHz optomechanical platforms already impose a quantum-regime bound $\Delta x_{\rm GW} = hL$9 on non-vacuum gravitational backgrounds. The approach elevates coherence protection from a radiation-theory curiosity to a structural probe of the quantum state of gravity, contingent on continued advances in mechanical coherence times and state-preparation fidelity.