Signatures of correlation of spacetime fluctuations in laser interferometers

TL;DR

This paper explores how laser interferometers can detect spacetime fluctuations predicted by gravity models by analyzing their signal patterns.

Contribution

The paper introduces a method to link interferometer signals to specific gravity models based on spacetime fluctuation correlations.

Findings

Three classes of spacetime fluctuations are characterized by their correlation functions and frequency behaviors.

Interferometers with and without arm cavities can detect these fluctuations using different sensitivity ranges.

Comparing modeled signals with experimental data helps constrain the strength and correlation scale of spacetime fluctuations.

Abstract

Spacetime fluctuations (SFs), a common feature of proposed gravity models, could be detected using laser interferometers. To advance this effort, we provide the correspondence between expected interferometer output signals and gravity models. We consider three classes of SFs, characterised by the decay behaviours and symmetries of their two-point correlation functions. For each, we identify the low- and high-frequency behaviour of the outputs and their dependence on the interferometer’s length. Capturing these requires sensitivity over a broad frequency range that spans the light-round-trip frequency, as provided by laboratory-scale setups, whereas detecting the presence or absence of SFs can be done with high narrowband sensitivity at the light-round-trip frequency. Our approach applies to interferometers with arm cavities, such as the km-long LIGO detectors, and those without, like laboratory-scale setups QUEST and GQuEST. Finally, we constrain the strength and correlation scale of SFs by comparing our modelled signals with experimental data.

Spacetime fluctuations have been predicted by a plethora of semiclassical and quantum models of gravity. Here, the authors discuss how characteristic signatures of classes of such spacetime fluctuations can be identified in laser interferometers.

Figures

MLI output for three correlation classes.Output signal power spectral densities (PSDs) of an MLI due to SFs. For correlation class (a), *factorised ρ*_F_, the scaled PSD $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${S}_{{{{\rm{NC}}}}}(\nu )=\frac{cS(f)}{{\Gamma }_{{{{\rm{S}}}}}{{{{\mathcal{L}}}}}^{3}}$$\end{document}$ vs. $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfo

Geometry of the MLI setup.a A schematic diagram of the Michelson laser interferometer (MLI). Each of the two interferometers in the Holometer^[28](#CR28)^, QUEST^[29](#CR29)^, and GQuEST^[25](#CR25)^ is an MLI. b Fabry-Pérot arm cavity design used in LIGO-type interferometers. Mirror M is introduced in each of the two arms of the MLI to create optical resonators (illustrated for arm D) with power transmissivity *T*_M_ = 1 − *R*_M_.

Excluding parameter values using data from Holometer and QUEST.Excluded region in the parameter space of *Γ*_S_ and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa={\ell }_{r}/{{{\mathcal{L}}}}$$\end{document}$ for correlation classes: (a) *factorised ρ*_F_, (b) *inverse ρ*_I*m*_, and (c) *exponential ρ*_E*m*_, where *m* = S (*m* = ST) denotes correlations depending on spatial (spacetime) separation. Data f

Full text

Introduction

A wide spectrum of ideas has been considered^1,2^ to understand the fundamental nature of gravity. Some, such as the idea of spacetime fluctuations (SFs), form a leitmotif in this effort. Since its first proposal by Wheeler^3^, SFs have been extensively examined^4^ in the context of a quantum description of gravity, as well as in different semiclassical models of gravity^5–11^ and the study of stochastic gravitational waves^12,13^. Spacetime has thus been hypothesised to be, for instance, classical but stochastic^14^, or classical but emerging from underlying quantum entanglement^15^, or holographic and having quantum perturbations^16^. These hypotheses suggest different mathematical forms for the correlation functions of SFs.

The scale of the correlations in SFs differs widely across models, ranging from the Planck length scale in effective field theories^17^ to long-range correlations in holographic models^18,19^. The latter hypothesise observable effects such as the violation of Lorentz invariance^20^, gravitational decoherence^2,21^, blurring of astronomical objects^22^, and interferometric noise^23–25^. Of these, the last has garnered attention through the development of gravitational-wave detectors^26^, and the Holometer^27,28^ experiment. Recently, sophisticated laboratory-scale Michelson laser interferometers (MLIs), such as QUEST^29^ and GQuEST^25^, which incorporate new quantum technologies, aim to search for SFs.

Unambiguous detection of SFs would constitute a breakthrough in understanding the fundamental nature of gravity. To enable such a detection, experimental designs require estimates of the strength and bandwidth of the expected interferometric output signal so they can be suitably optimised. Further, computing the interferometric output signal can help in understanding how SF signals, which are typically broadband, can be distinguished from instrumental noise. Moreover, quantitative theoretical predictions of the output signal are essential either to rule out the presence of SFs in interferometric data or, if the underlying gravity theory includes free parameters, to constrain those parameters. While many theoretical works have modelled the interferometric output signal to identify possible SF signatures^10,16,17,30–32^, they only offer model-specific predictions.

We take a more expansive approach in studying SFs. Irrespective of the classical or quantum description of a phenomenon, correlations in physical processes tend to decay either exponentially or polynomially with increasing separation \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r$$\end{document} between two spacetime points^33,34^. Exponential decay of correlations typically emanates from underlying physics that is short-ranged. In the context of gravity, this captures both quantum^15^ and semiclassical^11^ models. On the other hand, polynomial decay of correlations of the form r^−η^ ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta \in {{\mathbb{R}}}^{+}$$\end{document} ), emanate from long-range interactions. In the context of gravity, such correlations of SFs are expected to decay as a reciprocal of distance, i.e., 1/r^16,17,22,35^ or possibly with a higher inverse power-law such as r^−4/3^^36^. Furthermore, internal symmetries of the spacetime metric are reflected in induced symmetries of the correlation function, such as factorisation into spatial and temporal parts. This factorisation occurs in SF models such as the Oppenheim model^10^, a generalised Károlyházy model^36^, the continuous spontaneous localisation model^37,38^, and the Diósi-Penrose model^39–41^.

Thus motivated, we consider three possible classes of two-point correlation functions of the SFs: (a) factorised into spatial and temporal correlations, both decreasing with increasing spatial and temporal separation (for instance, a semiclassical model^10^); (b) an inverse of the separation between the two spacetime points (for instance, an SF model obeying the wave-equation in 3 + 1 dimensions^16,22,35^), with subclasses refining the definition of the separation as either a spatial or a spacetime separation; and (c) an exponential decay with the separation between the two spacetime points (for instance, due to quantum entanglement^15^, and in semiclassical^11^ models), with subclasses as in (b). Finding experimental evidence for or against one of the above classes of correlation functions describing SFs would constitute a breakthrough in elucidating the fundamental nature of gravity.

In this work, corresponding to each class of SFs, we identify three characteristic signatures: the low-frequency behaviour, the high-frequency behaviour and the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal{L}}}}$$\end{document} dependence of the interferometric output signal, where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal{L}}}}$$\end{document} is the arm length of the MLI. In order to observe all three signatures, it is sufficient if the MLIs are sensitive over two decades spanning the light round-trip frequency \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${f}_{{{{\rm{LRT}}}}}=c/(2{{{\mathcal{L}}}})$$\end{document} . However, even MLIs that are sensitive in a shorter span that includes fLRT, could allow observation of the three signatures to a limited extent. The projected sensitivities of laboratory-scale MLIs such as QUEST and GQuEST span the fLRT, unlike LIGO. Therefore, we find that these laboratory-scale MLIs allow observation of more signatures than LIGO in principle, thus providing more information on the class of the underlying correlation function and aiding in distinguishing between the classes (a)–(c).

We also show that MLIs with arm cavities have a significant advantage in detecting the bare presence or absence of SFs. This is due to the peak in the interferometric output signal for MLIs with arm cavities at their fLRT. We note that for LIGO, fLRT ≈ 37.5 kHz, which is outside the frequency span of the publicly available data. GWTC-3 run^42^ in LIGO reports data only up to 5 kHz. We await new results from experimental runs that report at a higher frequency range.

Our approach is agnostic to the microscopic origins of the SFs. It is important to emphasise that the microscopic origins of the SFs differ widely between the different gravity models, ranging from quantum entanglement of degrees of freedom resulting in suitable moments of the metric^15^ to geontropic^16^ or stochastic^11^ fluctuations of the metric itself. Our approach requires as input only the two-point correlation of the spacetime metric from any model of gravity and the geometry of the MLI. This enables us to compute the output signal power spectral density (PSD) of the MLI on which all our conclusions are based. Our approach also accommodates computing the signal PSD for MLIs with arm cavities, such as LIGO. Whether the presence of arm cavities provides any advantage towards detection has been much debated lately^16,18,43^. This debate even led to the exclusion of arm cavities from experimental designs of recent MLIs^25,28,29^. We now settle the debate for any SF from classes (a)−(c).

In addition, we compare the computed output signal PSD with the measured output signal from the Holometer and QUEST MLI experiments to constrain the free parameters in classes (a)−(c). We also compare our constraints with prior ones^16^.

Results

In this paper, we present three main results, namely, distinguishing between correlation classes (a)−(c) using three characteristic signatures identified in the output signal PSD of the MLI without arm cavities, the advantage enjoyed by the MLI with arm cavities in detecting the bare presence or absence of spacetime fluctuations, and constraining the free parameters in the correlation functions from different gravity models using the measured output signal from the Holometer and QUEST MLI experiments. To this end, we use the approach described in the “Methods” section to obtain the PSD of the optical path difference between the two arms of the MLI in the presence of random fluctuations w(r) at r ≡ (t, x, y, z) in the spacetime metric (Eq. (4)). Thus obtained, the interferometric output signal PSD for MLI without (Eq. (15)) and with (Eq. (22)) arm cavities involves the correlation function \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{w({t}_{1},{\vec{r}}_{1})w({t}_{2},{\vec{r}}_{2})}={\Gamma }_{{{{\rm{S}}}}}\,\rho \left(c{t}_{12},{\vec{r}}_{12}\right)$$\end{document} (Eq. (8)) where X_12_ = X_1_ − X_2_ ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\rm{X}}}}=t,\vec{r}$$\end{document} ) with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\vec{r}}_{i}\equiv \{{x}_{i},{y}_{i},{z}_{i}\}$$\end{document} (i = 1, 2) and ΓS is the strength of the SFs. Subsection “Correlation classes” in the “Methods” section lists the correlation functions defining the three correlation classes (a)−(c). Further, for the ease of comparison, we consider dimensionless frequency and PSDs defined in the “MLI output signal PSDs” subsection : ν (Eq. (24)), SNC (Eq. (25)), and SC (Eq. (26)). Note that we have two types of dimensionless PSDs: SNC independent of the correlation length ℓr, catering to correlation class (a) which is also independent of ℓr, and SC which depends on ℓr for correlation classes (b) and (c). From Eqs. (15, 17 and 22), it is clear that both SNC(ν) and SC(ν) are independent of ΓS.

We present the output signal PSDs for an MLI without (see Fig. 1(a)−(c)) and with arm cavities (see Fig. 2) for ρ corresponding to correlation classes: (a) factorised ρF (Eq. (10)), (b) inverse ρIm (Eqs. (11, 12)), and (c) exponential ρEm (Eqs. (13, 14)) where m = S, ST, with m = S (m = ST) denoting dependence on spatial (spacetime) separation.Fig. 1MLI output for three correlation classes.Output signal power spectral densities (PSDs) of an MLI due to SFs. For correlation class (a), factorised* ρF, the scaled PSD \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${S}_{{{{\rm{NC}}}}}(\nu )=\frac{cS(f)}{{\Gamma }_{{{{\rm{S}}}}}{{{{\mathcal{L}}}}}^{3}}$$\end{document} vs. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nu=\frac{\pi f}{2{f}_{{{{\rm{LRT}}}}}}$$\end{document} is plotted. For classes (b), inverse **ρ_Im_, and (c) exponential **ρ_Em_ with m = S, ST, the scaled PSD \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${S}_{{{{\rm{C}}}}}(\nu )=\frac{cS(f)}{{\Gamma }_{{{{\rm{S}}}}}{\ell }_{r}{{{{\mathcal{L}}}}}^{2}}$$\end{document} vs ν is plotted. Here, m = S (m = ST) denotes correlations depending on spatial (spacetime) separation. In (c), PSDs corresponding to ρES for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa={\ell }_{r}/{{{\mathcal{L}}}}=0.025$$\end{document} (red solid), κ = 0.01 (green dashed), κ = 0.005 (brown dot-dashed), and κ = 0.0025 (pink dotted), demonstrate dependence on κ. The PSD corresponding to ρEST for κ = 0.025 (blue points) is also plotted in (c). Small and large ν trends are as indicated by black dashed/dotted lines in (a)−(c). The black vertical line marks ν = 1.Fig. 2. Output of MLI with arm cavities.Scaled PSD SC(ν) corresponding to LIGO (red) with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal{L}}}}=4000$$\end{document}  m, QUEST (blue) with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal{L}}}}=3$$\end{document}  m and a hypothetical MLI without arm cavities (black dashed) with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal{L}}}}=4000$$\end{document}  m vs scaled frequency ν is plotted for (top) inverse ρIS, and (bottom) exponential ρES correlation functions with ℓ_r*_ = 0.03 m.

Distinguishing correlation functions

To list the three characteristic signatures of the interferometric output signal PSD corresponding to each correlation class, we first consider the case of a simple MLI without arm cavities, as in Fig. 3(a). We present the analytical expressions for the PSDs corresponding to each class (For more details, see Sec. V of the Supplementary Information for class (a) and Sec. VI of the Supplementary Information for classes (b) and (c)).Fig. 3. Geometry of the MLI setup.a A schematic diagram of the Michelson laser interferometer (MLI). Each of the two interferometers in the Holometer^28^, QUEST^29^, and GQuEST^25^ is an MLI. b Fabry-Pérot arm cavity design used in LIGO-type interferometers. Mirror M is introduced in each of the two arms of the MLI to create optical resonators (illustrated for arm D) with power transmissivity TM = 1 − RM.

Low-frequency limit* ν* ≪ 1: For class (b) (Fig. 1 (b)), the PSDs are proportional to ν^2^ with SC(ν = 0) = 0 in this limit. We analytically find the constants of proportionality as 2/(3π) and 1/2 for ρIS and ρIST respectively (see, Sec. VI of the Supplementary Information for details). This can also be understood intuitively, especially in the case of ρIS. In this particular case, the fluctuations satisfy the wave equation and therefore, the 4-d Fourier transform (Eq. (18)) of ρIS involves a Dirac delta function that forces \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2\pi f=c| {\vec{k}}_{1}|$$\end{document} (see Eq. (81) of the Supplementary Information for the exact form). This leaves the interferometer response function (Eq. (19)) as the only frequency-dependent contribution to the PSD, which, up to a constant factor, is \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{{\mathcal{L}}}}}^{2}{f}^{2}| {T}_{+}^{(x)}-{T}_{+}^{(z)}{| }^{2}/4$$\end{document} in this limit. The f^2^ in this term produces the quadratic trend of PSD in class (b).

For classes (a) and (c), the PSDs are almost flat in a log-log plot in the limit ν ≪ 1 in Fig. 1 (a), (c). We find the numerical fit in this limit for both cases to be of the form \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${e}^{-\alpha {\nu }^{2}}$$\end{document} where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha \in {\mathbb{R}}$$\end{document} is found numerically. In class (a), the exact numerical fit is \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${S}_{{{{\rm{NC}}}}}(\nu \, \ll \, 1)\approx 0.275{e}^{-\frac{5{\nu }^{2}}{2\pi }}$$\end{document} . This fit is independent of the value of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal{L}}}}$$\end{document} (see Sec. V of the Supplementary Information for details). In class (c), the value α in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${e}^{-\alpha {\nu }^{2}}$$\end{document} depends on the specific value of the ratio \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa={\ell }_{r}/{{{\mathcal{L}}}}$$\end{document} . For instance, α = 4/π for κ = 0.01 (see Sec. VI of the Supplementary Information for details).

High-frequency limit* ν* ≫ 1: As is evident from Fig. 1 (a)−(c), the PSDs for all correlation classes decrease with increasing ν. However, the rates at which they decay are starkly different. In Sec. VI of the Supplementary Information, we analytically justify why we find the following decay rates of SC(ν): ∝ 1/ν for ρIS, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\propto 1/\sqrt{\nu }$$\end{document} for ρIST, and ∝ 1/ν^2^ for ρES. We also present SNC(ν) vs ν for ρES in Fig. 4. Though the PSD in this case is evidently dependent on κ, and therefore ℓr, we use SNC(ν) to showcase the high-frequency behaviour, which is always ∝ 1/ν^2^ irrespective of the value of κ (as this is not very apparent from Fig. 1(c)). We also note here that the onset of this decay is delayed as the ratio κ decreases.Fig. 4κ-dependence of MLI output in class (c).Scaled output signal PSDs SNC(ν) = κ **SC(ν) vs ν are plotted for exponential* ρ*ES corresponding to κ = 0.025 (red solid), κ = 0.01 (green dashed), κ = 0.005 (brown dot-dot-dashed), and κ = 0.0025 (pink dotted), demonstrating the 1/ν^2^ behaviour at the high-frequency limit.

The other decay rates, such as SNC ∝ 1/ν^2^ for class (a) and SC(ν) ≈ 0.03 ν^−0.3^ e^−0.01ν^ for ρEST, are obtained through numerical fits. Of these, the functional form of the latter is obtained by considering the possible form of the cosine transform of ρEST (see Sec. VI of the Supplementary Information for details).

Dependence on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal{L}}}}$$\end{document} and κ: For class (a), we analytically show that SNC(ν) is independent of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal{L}}}}$$\end{document} in Eq. (68) of the Supplementary Information, and that SC(ν) for class (b) is independent of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal{L}}}}$$\end{document} in Eqs. (69, 70) of the Supplementary Information. For class (c), we analytically show that SC(ν) depends on κ in Eqs. (74-76) of the Supplementary Information. Figure 1(c) illustrates this κ-dependence of the PSDs.

As shown above, the different classes of spacetime fluctuations produce three characteristic signatures in their corresponding output signal PSDs. These can be used to identify the nature of the underlying SFs from interferometric data. Here, it is worth pointing out another feature that is present in the scaled PSDs of the three classes is the oscillatory behaviour at ν > 1. However, this feature originates from the interferometer response (See Eqs. (49-50) of the Supplementary Information for details) and is not a characteristic signature of the different classes.

To observe all three signatures characteristic of each correlation class, the interferometer should ideally be sensitive in the range of ν ≈ 0.1 to ν ≈ 10. In QUEST with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal{L}}}}=3$$\end{document}  m (resp., GQuEST with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal{L}}}}=5$$\end{document}  m), the sensitive bandwidth is designed to span from 1 MHz to 250 MHz^29^ (resp., 8 MHz to 40 MHz^25^), with a corresponding span of 0.03 ⩽ ν ⩽ 78 (resp., 0.42 ⩽ ν ⩽ 2.1). This illustrates that the bandwidths of both QUEST and GQuEST would allow observation of all three signatures, although the narrower bandwidth of GQuEST could limit the observation of the low- and high-frequency signatures to some extent.

On the other hand, experimental data from LIGO covers the frequency range only from about 10 Hz to 10 kHz, corresponding roughly to 0.0004 < ν < 0.4. However, for completeness, we discuss the low- and high-frequency behaviour of the PSD for LIGO in Sec. VIII of the Supplementary Information.

To summarise, QUEST and GQuEST, with their broader bandwidths, allow observation of all the characteristic signatures that could help in distinguishing between correlation functions using their interferometric PSD data.

Detecting SFs

To highlight the advantage LIGO enjoys in detecting the SFs, we list three key features of the interferometric output signal of the MLIs with arm cavities (see Fig. 3(b) for the geometry, and see Fig. 2 for the PSDs) for the different correlation classes. For further details, see Sec. VIII of the Supplementary Information.

1. For any fluctuation described by a Gaussian random process, the signal PSD (Eq. (22)) of MLIs with arm cavities has peaks at ν = m**π/2, for m = 1, 2, 3, ⋯  , of magnitude \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${T}_{{{{\rm{M}}}}}^{4}/{\left(1-\sqrt{{R}_{{{{\rm{M}}}}}}\right)}^{6}$$\end{document} . For LIGO with TM = 1 − RM = 0.014 of the input mirror of the arm cavity, the magnitude of the peak is ≈ 3.2 × 10^5^.
2. For class (b), SC for an MLI without arm cavities is independent of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal{L}}}}$$\end{document} , with a global maximum at ν = 1. Therefore, in this class, the Fabry-Pérot cavity response (Eq. (23)) enhances the signal strength (Eq. (22)) at every frequency ν. The strongest signal is at ν = π/2 (equivalently, f ≈ 37.5 kHz for LIGO), as expected. This is illustrated in the top panel of Fig. 2. This is also consistent with prior work^16^. A similar argument holds for class (a) as its SNC for an MLI without arm cavities is also independent of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal{L}}}}$$\end{document} .
3. For class (c), SC is directly proportional to κ. We compare SC for a given ℓr in two different setups, with and without arm cavities and with different arm lengths: LIGO with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal{L}}}}=4$$\end{document}  km, and QUEST with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal{L}}}}=3$$\end{document}  m. It is evident from the values of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal{L}}}}$$\end{document} that the ratio κ in QUEST is far greater than that in LIGO for a given ℓr. Thus, SC ∝ κ implies that the LIGO signal is reduced with respect to that of QUEST by the factor ≈ 10^−3^ (i.e., the ratio \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal{L}}}}$$\end{document} of QUEST to that of LIGO). This should be considered in conjunction with Feature 1 in this list (the presence of a peak at ν = π/2 in the LIGO signal is always enhanced by the Fabry-Pérot cavity gain ≈ 10^5^). Therefore, when measuring SF with any ℓr in these setups, SC at ν = π/2 of LIGO has a peak that exceeds the SC of QUEST. This is illustrated in the bottom panel of Fig. 2.

It is thus evident that for the classes (a)−(c) considered, LIGO has a clear advantage over QUEST and GQuEST in detecting the presence of SFs.

It is also evident from our arguments that this advantage of LIGO is not guaranteed for all correlation functions. For instance, let us consider a class of correlation functions for which SC ∝ κ^2^ for an MLI without arm cavities. For such a class of correlation functions, using arguments similar to those used in discussing Feature 3, we can see that the peak in the LIGO signal does not exceed the PSD of QUEST. Therefore, LIGO does not enjoy an advantage in detection for such a class. This line of thinking might help in understanding some gravity models that predict LIGO should not enjoy any advantage in detecting SFs^18,43^.

Comparison with experimental data

We compare our results with experimental data from the Holometer [^44^, Fig. 12], a now-retired MLI experiment with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal{L}}}}=40$$\end{document}  m, and QUEST [^45^, Fig. 4], a table-top MLI experiment ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal{L}}}}=1.8$$\end{document} m) that has just produced its first data. Note that we use QUEST data from the full measurement bandwidth, only part of which is published in Fig. 4 of^ 45^. While both the Holometer and QUEST experiments comprise twin MLIs designed to be cross-correlated, we restrict the present comparison to single-MLI PSDs for generality, avoiding additional assumptions on the correlation length of the SFs. We also do not consider data from LIGO or other gravitational-wave (GW) detectors, as these instruments operate and record data in a bandwidth that spans a relatively small range of ν and does not include the highest magnitude of the signal PSD, which is at ν = π/2.

The data from the Holometer and QUEST is obtained using the measurement technique in ref. ^46^, so the output signal sensitivity is improved by removing quantum shot noise, and potentially revealing stationary signals and classical noise produced in the MLI. The data is dominated by classical noise, and no signal was detected (i.e., the three characteristic signatures mentioned earlier cannot be discerned in the data). Thus, we can constrain the two parameters ΓS and κ (equivalently, ℓr). To obtain the constraints, we consider the magnitude of the PSD S(ν = νi) of the SFs described by each correlation class to be smaller than the experimentally measured spectral density Smeas(ν = νi), at frequencies {νi}. This experimentally measured Smeas(ν = νi) contains negligible measurement uncertainty for the data considered. For instance, in the Holometer (respectively, QUEST) experiment, the constraint is obtained by S(ν = νi) ⩽ Smeas(ν = νi) at every νi in the range 0.4 ⩽ ν ⩽ 4 (resp., 0.38 ⩽ ν ⩽ 1.8). The constraints are listed in Table 1 and plotted in Fig. 5 for classes (a)−(c). It is to be noted that the plot (Fig. 5) is a log-log plot of the parameter space {ΓS, κ}. Each grid line on the ΓS-axis marks an increase by a factor of 10, while that on the κ-axis corresponds to an increase by a factor of 10^0.2^ ≈ 1.6.Table 1. Constraint on ΓS and κ using experimental data from the Holometer and QUEST experiments for different correlation classesCorrelation classExperimentPrior boundsHolometerQUEST(a) Factorised ρFΓS ⩽ 1.628 × 10^−34^ΓS ⩽ 5.154 × 10^−33^None(b) Inverse ρIm(b1) Spatial separation m = SΓSκ ⩽ 6.075 × 10^−35^or \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha \equiv \frac{2{\Gamma }_{{{{\rm{S}}}}}\kappa {{{\mathcal{L}}}}}{3\sqrt{2\pi }{\ell }_{{{{\rm{P}}}}}}\leqslant 40$$\end{document} ΓSκ ⩽ 7.022 × 10^−33^or α ⩽ 208α ⩽ 0.6^16^(b2) Spacetime separation m = STΓSκ ⩽ 2.264 × 10^−35^ΓSκ ⩽ 3.537 × 10^−33^None(c) Exponential ρEm(c1) Spatial separation m = SΓSκ^2^ ⩽ 1.25 × 10^−35^ΓSκ^2^ ⩽ 6.173 × 10^−32^None(c2) Spacetime separation m = STΓSκ^2^ ⩽ 2.812 × 10^−36^ΓSκ^2^ ⩽ 1.234 × 10^−32^NoneSee discussion in the main text on prior bounds.Fig. 5. Excluding parameter values using data from Holometer and QUEST.Excluded region in the parameter space of ΓS and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa={\ell }_{r}/{{{\mathcal{L}}}}$$\end{document} for correlation classes: (a) factorised **ρF, (b) inverse **ρIm, and (c) exponential* ρ_Em_, where m = S (m = ST) denotes correlations depending on spatial (spacetime) separation. Data from the Holometer excludes the regions shaded in blue and red; data from QUEST excludes the blue region only. The black dashed (ΓSκ*^2^ ⩽ 6.173 × 10^−32^), dot-dashed (ΓSκ^2^ ⩽ 1.25 × 10^−35^), dotted (ΓSκ^2^ ⩽ 1.234 × 10^−32^), and dot-dot-dashed (ΓSκ^2^ ⩽ 2.812 × 10^−36^) lines denote curves fitted to the boundary of the excluded regions.

(a) Factorised correlation function: Considering ρF (Eq. (10)) is independent of ℓr, the constraint obtained involves only ΓS. We know that SNC is independent of ΓS and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal{L}}}}$$\end{document} from Eq. (68) of the Supplementary Information. Therefore, the constraint is given by

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Gamma }_{{{{\rm{S}}}}}\leqslant \,{{{\rm{min}}}}\,\left\{\frac{c\,{S}_{{{{\rm{meas}}}}}({\nu }_{i})}{{{{{\mathcal{L}}}}}^{3}\,{S}_{{{{\rm{NC}}}}}({\nu }_{i})}\right\}.$$\end{document}

The constraints obtained using the above expression are quoted in Table 1 and plotted in Fig. 5 (a). We note that we consider 10^−4^ ⩽ κ ⩽ 5 × 10^−2^. This range is chosen such that the two following conditions are satisfied: (1) κ ≪ 1 to ensure that the path difference remains a stationary random process, and (2) \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa \gg \lambda /{{{\mathcal{L}}}}$$\end{document} , where λ is the wavelength of light ( ~10^−6^ m for Holometer and QUEST), to ensure that the eikonal approximation remains valid. There are no constraints on the fluctuation strength considering ρF from prior work.

(b) Inverse correlation function: The constraint in each subclass is given by

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{r}{\Gamma }_{{{{\rm{S}}}}}\kappa \leqslant \,{{{\rm{min}}}}\,\left\{\frac{c\,{S}_{{{{\rm{meas}}}}}({\nu }_{i})}{{{{{\mathcal{L}}}}}^{3}\,{S}_{{{{\rm{C}}}}}({\nu }_{i})}\right\},\end{array}$$\end{document}

where SC(ν) is given by Eqs. (69,70) of the  Supplementary Information with m = S (respectively, m = ST) for the subclass (b1) (resp., (b2)). Equation (2) yields the constraints listed in Table 1 and plotted in Fig. 5(b1), (b2).

In the case of the subclass (b1), the constraints are compared with the constraint α ⩽ 0.6 in the Pixellon model^16^. This is because

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha \equiv \frac{2{\Gamma }_{{{{\rm{S}}}}}\kappa {{{\mathcal{L}}}}}{3\sqrt{2\pi }{\ell }_{{{{\rm{P}}}}}},$$\end{document}

where ℓP is the Planck length. In^16^, the constraint on α is obtained by comparing the expected signal PSD with the cross-spectral density (CSD) data from the twin interferometers in the Holometer (~10^−41^ m^2^ Hz^−1^). This is in contrast to our comparison with the PSD data from a single interferometer (~10^−39^ m^2^ Hz^−1^ for the Holometer), yielding the constraint α ⩽ 40. The two orders of magnitude difference between the constraints on α obtained using data from the Holometer can be explained by a corresponding difference between the magnitude of the CSD and PSD from the same experiment. In order to use the CSD data, an implicit assumption that the correlation length is larger than the distance between the two interferometers is made. We avoid such an assumption by restricting our comparison to the PSD data of the single interferometer. In subclass (b2), there are no constraints from prior work.

(c) Exponential correlation function: The constraints are obtained by setting S(ν = νi) ⩽ Smeas(ν = νi) at every νi. Here we numerically compute \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S(\nu={\nu }_{i})={\Gamma }_{{{{\rm{S}}}}}{\ell }_{r}{{{{\mathcal{L}}}}}^{2}{S}_{{{{\rm{C}}}}}(\nu={\nu }_{i})/c$$\end{document} using Eqs. (74-76) of the Supplementary Information. A numerical fit to such constraints yields the expressions listed in Table 1. In this class, there are no constraints from prior work.

Discussion

We have developed a methodology to systematically compute interferometric output signal power spectral densities, produced by statistically defined spacetime fluctuations (SFs), under an explicit list of assumptions. Using this methodology, we have compiled the interferometric output signals due to SFs in three correlation classes, for Michelson laser interferometers with and without arm cavities. This allows us to identify characteristic signatures in spectral densities for the different classes of correlation functions of SFs. We also find that (1) the laboratory-scale QUEST and GQuEST will have the broad bandwidth needed to observe all the characteristic signatures, while (2) LIGO is better suited for detecting the bare presence or absence of SFs. By comparing with the experimental data from the Holometer and QUEST, we have also constrained the strength and correlation scale of the SFs in the classes (a)−(c), which tallies with constraints obtained in contemporary work^16^.

Moreover, our methodology enables unambiguous computation of interferometric signals for other (current or future) theories of gravity just from the correlation function of the SFs and the geometry of the interferometer. It can also be applied to compute interferometric signals to search for stochastic gravitational waves^12,13^ or dark matter^47^. Lastly, our methodology may be applied in instrumental ‘noise hunting’ or calibration efforts for interferometers, for cases where the noise or calibration signal can be described as metric or phase fluctuations^48–50^ along the light path.

Future work could apply our methodology to other gravity models, such as continuous spontaneous localisation and Diósi-Penrose models, to constrain their free parameters using interferometric experiments. Further research could also extend our methodology to include more complicated interferometer configurations and multi-interferometer cross-correlations, which will then be used to further constrain the parameters of SFs.

Methods

Modelling light propagation

To investigate the effects of a fluctuating spacetime, we consider an isotropic spacetime metric g^β**γ^ (β, γ = 0, 1, 2, 3) of the form

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${g}^{00}=- \! 1+2w({{{\boldsymbol{r}}}}),\quad {g}^{ij}={\delta }_{ij},\quad {g}^{0i}={g}^{i0}=0,$$\end{document}

where i, j = 1, 2, 3 and w(r) is a random process in r ≡ (t, x, y, z), with w ≪ 1. Here 0 (resp., i ∈ {1, 2, 3}) corresponds to the timelike (resp., spacelike) component. While we consider this specific form in this work, our approach can encompass general fluctuations in every g^β**γ^ (see Sec. I of the Supplementary Information for details). To model the propagation of light of frequency Ω and wavelength λ = 2π**c/Ω in the given spacetime manifold, we solve for the electromagnetic tensor in the relativistic wave equation (RWE) (see Sec. I of the Supplementary Information with ref. ^51^ therein for details). This is subject to the following assumptions.

**Assumption (i) **Setting length and time scales: The correlation scales in length and time of the metric fluctuations w need to be longer than λ and 2π/Ω respectively. This effectively sets the wavelength as the smallest length scale in the system, i.e., the eikonal approximation. This allows us to neglect diffraction due to the fluctuations.

Applying the eikonal approximation (k = 2π/λ → ∞), the RWE reduces to the light propagating along the null geodesic for any general spacetime metric^52^. For light propagating along, say, the z-axis, the electric field solution for the RWE (see Sec. I of the Supplementary Information for details) is

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\vec{E}({{{\boldsymbol{r}}}}(t))={\vec{E}}_{in}(x,y){e}^{ik\Phi ({{{\boldsymbol{r}}}}(t))},$$\end{document}

where

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Phi ({{{\boldsymbol{r}}}}(t))=ct-z+{\Phi }_{{{{\rm{F}}}}}({{{\boldsymbol{r}}}}(t)),$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\Phi }_{{{{\rm{F}}}}}({{{\boldsymbol{r}}}}(t))=c\int_{0}^{t}\,{{{\rm{d}}}}\,{\!t}^{{\prime} }\,w({{{\boldsymbol{r}}}}({t}^{{\prime} })),$$\end{document}

and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\vec{E}}_{in}(x,y)$$\end{document} is the input transverse profile of the beam.

**Assumption (ii) **Slowly varying envelope approximation (SVEA): The above solution (5) also uses the SVEA, which is consistent with the eikonal approximation. It assumes a very small rate of metric-fluctuation-induced phase fluctuations ∂tΦF ≪ c and ∂iΦF ≪ 1 (i = x, y, z).

We now introduce the assumptions on the random SF process w.

Assumption (iii)* Stationarity:*w(r) is a stationary Gaussian random process with the expectation values,

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{w}=0,\quad \,{{\rm{and}}}\,$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{w({t}_{1},{\vec{r}}_{1})w({t}_{2},{\vec{r}}_{2})}={\Gamma }_{{{{\rm{S}}}}}\,\rho \left(c{t}_{12},{\vec{r}}_{12}\right).$$\end{document}

Here X_12_ = X_1_ − X_2_ ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\rm{X}}}}=t,\vec{r}$$\end{document} ) with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\vec{r}}_{i}\equiv \{{x}_{i},{y}_{i},{z}_{i}\}$$\end{document} (i = 1, 2) and ΓS is the strength of the SFs. Both ΓS and ρ are dimensionless, consistent with w being dimensionless. Subsequently, we consider the correlation function \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho \left(c{t}_{12},{\vec{r}}_{12}\right)$$\end{document} to be from the classes (a)−(c) to obtain the corresponding MLI output as illustrated in Fig. 1.

Assumption (iv)* Isotropy:* The two-point correlation function ρ is isotropic in space, i.e.,

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{array}{l}\rho \left({\delta }_{0},\{{\delta }_{1},{\delta }_{2},{\delta }_{3}\}\right)=\rho \left({\delta }_{0},\{{\delta }_{3},{\delta }_{1},{\delta }_{2}\}\right)\\=\rho \left({\delta }_{0},\{{\delta }_{3},{\delta }_{2},{\delta }_{1}\}\right)=\ldots,\end{array}$$\end{document}

for any separation δμ (μ = 0, 1, 2, 3). We consider the correlation length of ρ to be ℓr for all three spatial dimensions, and the correlation time to be ℓr/c.

Correlation classes

We now list the different correlation classes, which completes the description of the spacetime fluctuations.

(a) Factorised correlation function: This class of correlation functions of the form

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\rho }_{{{{\rm{F}}}}}\left(c{\Delta }_{{{{\rm{t}}}}},{\vec{\Delta }}_{{{{\rm{r}}}}}\right)=-\left(\frac{\parallel {\vec{\Delta }}_{{{{\rm{r}}}}}\parallel }{{\ell }_{r}}\right)\,\left(\frac{{\ell }_{r}\,\delta ({\Delta }_{{{{\rm{t}}}}})}{c}\right),$$\end{document}

is motivated by the Oppenheim model^10^, where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\vec{\Delta }}_{{{{\rm{r}}}}}$$\end{document} is any 3-vector in space and Δt is any time interval. We note that multiplicative factors consistent with Assumption (iv) are suitably introduced such that the delta-correlated ρF remains dimensionless. Also, it is evident that ρF is independent of ℓr. Thus, the only parameter in this correlation class is the fluctuation strength ΓS that scales the correlation function (Eq. (8)).

(b) Inverse correlation functions: This class, a subset of the polynomial decay of the correlations, is motivated by models such as those of Karolyhazy^35^, and Zurek^16^, as well as effective field theories^17^. We consider two sub-classes where the correlations decay as a function of spatial separation and spacetime separation, respectively, with the specific form of the correlation functions defining the subclass listed below.

(b1) Spatial separation:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\rho }_{{{{\rm{IS}}}}}\left(c{\Delta }_{{{{\rm{t}}}}},{\vec{\Delta }}_{{{{\rm{r}}}}}\right)=\frac{{\ell }_{r}}{\parallel {\vec{\Delta }}_{{{{\rm{r}}}}}\parallel }\Theta (\parallel {\vec{\Delta }}_{{{{\rm{r}}}}}\parallel -c| {\Delta }_{{{{\rm{t}}}}}| ),$$\end{document}

where the Heaviside theta function \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Theta (\parallel {\vec{\Delta }}_{{{{\rm{r}}}}}\parallel -c| {\Delta }_{{{{\rm{t}}}}}| )$$\end{document} vanishes for all pairs of time-like separated spacetime points. Such correlation is found in models that assume the fluctuations satisfy the wave equation^16,17,35^.

(b2) Spacetime separation:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\rho }_{{{{\rm{IST}}}}}\left(c{\Delta }_{{{{\rm{t}}}}},{\vec{\Delta }}_{{{{\rm{r}}}}}\right)=\frac{{\ell }_{r}\,\Theta (\parallel {\vec{\Delta }}_{{{{\rm{r}}}}}\parallel -c| {\Delta }_{{{{\rm{t}}}}}| )}{\sqrt{\parallel {\vec{\Delta }}_{{{{\rm{r}}}}}{\parallel }^{2}-{c}^{2}{\Delta }_{{{\rm{t}}}\,}^{2}}}.$$\end{document}

This is a generalisation of (b1). Carrying forward the step function from (b1), the sign convention in the spacetime separation is chosen such that the correlation remains a real-valued function.

While we use two parameters to characterise strength, ΓS and the scale ℓr, in this class, the function depends only on the product ΓSℓr, rendering the parameters degenerate.

(c) Exponential correlation functions: This class of correlation functions covers models motivated by entanglement between holographic degrees of freedom^15^ or a mesoscopic interpretation of gravity^11^. We again list the correlation functions for two sub-classes based on

(c1) Spatial separation:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\rho }_{{{{\rm{ES}}}}}\left(c{\Delta }_{{{{\rm{t}}}}},{\vec{\Delta }}_{{{{\rm{r}}}}}\right)={e}^{-\frac{\parallel {\vec{\Delta }}_{{{{\rm{r}}}}}\parallel }{{\ell }_{r}}}\Theta (\parallel {\vec{\Delta }}_{{{{\rm{r}}}}}\parallel -c| {\Delta }_{{{{\rm{t}}}}}| ).$$\end{document}

(c2) Spacetime separation:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\rho }_{{{{\rm{EST}}}}}\left(c{\Delta }_{{{{\rm{t}}}}},{\vec{\Delta }}_{{{{\rm{r}}}}}\right)={e}^{-\frac{\sqrt{\parallel {\vec{\Delta }}_{{{{\rm{r}}}}}{\parallel }^{2}-{c}^{2}{\Delta }_{{{{\rm{t}}}}}^{2}}}{{\ell }_{r}}}\,\Theta (\parallel {\vec{\Delta }}_{{{{\rm{r}}}}}\parallel -c| {\Delta }_{{{{\rm{t}}}}}| ).$$\end{document}

Class (c) is a true two-parameter model, unlike classes (a) and (b). Classes (b) and (c) cannot be factorised into spatial and temporal functions.

MLI output signal PSDs

In an MLI (Fig. 3 (a)), light propagates from a laser source at the input port A to a detector at the output port B via the two perpendicular arms, denoted by C and D. The 50/50 lossless beamsplitter is denoted by BS and is taken as the origin of the reference frame in our computation. The beamsplitter and the end mirrors are assumed to be perfect due to the very small optical losses that are typical in these high-precision interferometers. We can effectively assume the detector to be at the origin, as the effect of any phase fluctuations suffered by the light after interference at the BS is negligible compared to that of phase fluctuations incurred in the arms in practice. The arm length of the MLI without fluctuations is \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\mathcal{L}}}}$$\end{document} and light-round-trip time τ0 = 1/fLRT (we recall that \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${({f}_{{{{\rm{LRT}}}}})}^{-1}=2{{{\mathcal{L}}}}/c$$\end{document} ). The PSD of the optical path difference between the two arms (see Sec. II of the Supplementary Information for details) at the detector is then written as a cosine transform from time separation Δτ to frequency f.

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S(f)=\frac{{c}^{2}{\Gamma }_{{{{\rm{S}}}}}}{2\pi }\int_{0}^{\infty }\,{{\rm{d}}}\,{\Delta }_{\tau }\left[\sigma ({\Delta }_{\tau })-\xi ({\Delta }_{\tau })\right]\cos 2\pi f{\Delta }_{\tau },$$\end{document}

where

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma ({\Delta }_{\tau })= \int\limits_{0}^{{\tau }_{0}}\,{{{\rm{d}}}}{t}_{1}\int\limits_{0}^{{\tau }_{0}}{{{\rm{d}}}}\,{\!t}_{2}\,\rho (c({t}_{1}+{\Delta }_{\tau }-{t}_{2}),0,0,s({t}_{1})-s({t}_{2}))\\ = \int\limits_{0}^{{\tau }_{0}}\,{{{\rm{d}}}}{t}_{1}\int\limits_{0}^{{\tau }_{0}}{{{\rm{d}}}}\,{\!t}_{2}\,\rho (c({t}_{1}+{\Delta }_{\tau }-{t}_{2}),s({t}_{1})-s({t}_{2}),0,0),$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\xi ({\Delta }_{\tau })=\int_{0}^{{\tau }_{0}}\,{{{\rm{d}}}}\,{t}_{1}\,\int_{0}^{{\tau }_{0}}\,{{{\rm{d}}}}\,{t}_{2} \: \rho (c({t}_{1}+{\Delta }_{\tau }-{t}_{2}),s({t}_{1}),0,-s({t}_{2})).$$\end{document}

Also, s(t) = *c *t if t ⩽ τ0/2 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$s(t)=2{{{\mathcal{L}}}}-c t$$\end{document} if t > τ0/2. Assumptions (i)−(iv) are used to obtain Eq. (15), with Eq. (16a) using Assumption (iv). We reiterate that the expressions obtained in this subsection are true irrespective of the form of the correlation function ρ and are not limited to classes (a)−(c). σ arises from correlations in the spacetime metric fluctuations within an arm of the interferometer, and ξ corresponds to correlations of the metric fluctuations between the two arms. Here, we have also assumed that the width of the light beams are of the order of the wavelength and negligible (see Sec. II of the Supplementary Information for details).

A response-function-based approach allows straightforward extension of our methodology to different MLI geometries and detection schemes. These include for instance, the cross spectral density (CSD) of the output signal from two co-located, co-aligned MLIs (see Sec. II and III of the Supplementary Information for details) and the signal PSD from an MLI with Fabry-Pérot arm cavities (see Sec. IV of the Supplementary Information for details) in the presence of SFs. Here it is important to note that both our approaches of obtaining the spectral densities, namely, based on correlation integral in Eq. (15) and the response function in Eq. (17) assume weak stationarity of the optical path difference. However, the existence of such a stationarity is not guaranteed, even when assuming the underlying SFs to be stationary and Gaussian. Therefore, when computing the PSDs, we check if the covariance of the optical path difference is positive-definite, to ensure weak stationarity of the optical path difference in our investigation. We find in class (c), this holds only when \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa={\ell }_{r}/{{{\mathcal{L}}}} \, \ll \, 1$$\end{document} .

We thus rewrite the PSD in terms of the corresponding *interferometer response function * \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\widetilde{\chi }}_{{{{\rm{I}}}}}(f,{\vec{k}}_{1})$$\end{document} as an integral over a 3-dimensional reciprocal space (see Sec. III of the Supplementary Information for details).

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S(f)=\int{{{\rm{d}}}}^{3}{\vec{k}}_{1}\,{\Gamma }_{{{{\rm{S}}}}}\,\widetilde{\rho }(2\pi f,{\vec{k}}_{1})\,{\widetilde{\chi }}_{{{{\rm{I}}}}}(f,{\vec{k}}_{1}),$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\widetilde{\rho }(2\pi f,{\vec{k}}_{1})$$\end{document} is given by

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\widetilde{\rho }(2\pi f,{\vec{k}}_{1})=\frac{1}{{(2\pi )}^{4}}\,\int{{{\rm{d}}}}^{3}{\vec{r}}_{12}\int_{-\infty }^{\infty }\,{{{\rm{d}}}}\,{t}_{12}\,\rho \left(c{t}_{12},{\vec{r}}_{12}\right)\,{e}^{-i\left(2\pi f{t}_{12}+{\vec{k}}_{1}\cdot {\vec{r}}_{12}\right)},$$\end{document}

and the response function \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\widetilde{\chi }}_{{{{\rm{I}}}}}(f,{\vec{k}}_{1})$$\end{document} of the MLI (Fig. 3 (a)) is

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\widetilde{\chi }}_{{{{\rm{I}}}}}(\,f,{\vec{k}}_{1})={\left(\frac{{{{\mathcal{L}}}}}{2}\right)}^{2}\,{\left\vert {C}_{x}(\,f,{\vec{k}}_{1})-{C}_{z}(\,f,{\vec{k}}_{1})\right\vert }^{2},$$\end{document}

with

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${C}_{j}(\,f,{\vec{k}}_{1})={e}^{if{T}_{+}^{(j)}}\left\{\,{{{\rm{Sinc}}}}\,\left(f{T}_{+}^{(j)}\right)+{e}^{\frac{2\pi if{{{\mathcal{L}}}}}{c}}\,{{{\rm{Sinc}}}}\,\left(f{T}_{-}^{(j)}\right)\right\},$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${T}_{\pm }^{(j)}(\,f,{\vec{k}}_{1})=\frac{\pi {{{\mathcal{L}}}}}{c}\left(1\pm \frac{c}{2\pi f}\,{\vec{k}}_{1}\cdot {\hat{e}}_{j}\right),\quad (\,j=x,z).$$\end{document}

It is interesting to contrast the general 4-dimensional Fourier transform considered in Eq. (18) with the plane wave expansions of ‘metric perturbations’ considered in contemporary investigations of stochastic gravitational wave backgrounds^53^. These investigations consider stochastic gravitational waves as perturbations of the metric in the transverse-traceless gauge, and they require such perturbations to be a sum of plane waves that satisfy the wave equation. In contrast, we do not expect the stationary Gaussian SFs in the non-relativistic limit to satisfy the wave equation.

We also find the response function corresponding to an MLI with Fabry-Pérot arm cavities (see Fig. 3 (b) for the interferometer geometry; also see Sec. IV of the Supplementary Information for details). The signal PSD in terms of this response function is given by

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S(\,f)={\widetilde{\chi }}_{{{{\rm{FP}}}}}(\,f)\,\int{{{\rm{d}}}}^{3}{\vec{k}}_{1}\,{\Gamma }_{{{{\rm{S}}}}}\,\widetilde{\rho }(2\pi f,{\vec{k}}_{1})\,{\widetilde{\chi }}_{{{{\rm{I}}}}}(\,f,{\vec{k}}_{1}),$$\end{document}

where the Fabry-Pérot cavity response is given by

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\widetilde{\chi }}_{{{{\rm{FP}}}}}(\,\,f)={T}_{{{{\rm{M}}}}}^{4}{\left(\frac{1}{1-\sqrt{{R}_{{{{\rm{M}}}}}}}\right)}^{4} \left(\frac{1}{1+{R}_{{{{\rm{M}}}}}-2\sqrt{{R}_{{{{\rm{M}}}}}}\,\cos (2\pi f/{f}_{{{{\rm{LRT}}}}})}\right),$$\end{document}

and TM = 1 − RM is the power transmissivity, and RM is the power reflectivity of the mirror M rendering the arm cavity in either arm.

To assess and compare the behaviour of the computed signal PSDs for the different correlation function classes, we consider the dimensionless frequency and PSDs.

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nu \; \doteq \, \pi f/(2{f}_{{{{\rm{LRT}}}}})$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${S}_{{{{\rm{NC}}}}}(\nu ) \, \doteq \,\left(\frac{c}{{\Gamma }_{{{{\rm{S}}}}}\,{{{{\mathcal{L}}}}}^{3}}\right)\,S(\,\,f),\,{{\rm{and}}}\,$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${S}_{{{{\rm{C}}}}}(\nu ) \, \doteq \,\left(\frac{c}{{\Gamma }_{{{{\rm{S}}}}}\,{\ell }_{r}\,{{{{\mathcal{L}}}}}^{2}}\right)\,S(\,\,f).$$\end{document}

Note that we have two types of dimensionless PSDs: SNC independent of ℓr, catering to correlation class (a), which is also independent of ℓr, and SC, which depends on ℓr for correlation classes (b) and (c). From Eqs. (15, 17 and 22), it is clear that both SNC(ν) and SC(ν) are independent of ΓS.

Supplementary information

Supplementary Information Transparent Peer Review file