# Exploring gravity with the MIGA large scale atom interferometer

B. Canuel, A. Bertoldi, L. Amand, E. Pozzo Di Borgo, T. Chantrait, C. Danquigny, B. Fang, A. Freise, R. Geiger, J. Gillot, et al.

► To cite this version:

B. Canuel, A. Bertoldi, L. Amand, E. Pozzo Di Borgo, T. Chantrait, et al. Exploring gravity with the MIGA large scale atom interferometer. Scientific Reports, 2018, 8 (1), pp.14064. 10.1038/s41598-018-32165-z. hal-01874243

HAL Id: hal-01874243

https://hal.science/hal-01874243

Submitted on 14 Sep 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

# Exploring gravity with the MIGA large scale atom interferometer

B. Canuel<sup>1,2</sup>, A. Bertoldi<sup>1,2</sup>, L. Amand<sup>1,3</sup>, E. Pozzo di Borgo<sup>1,4</sup>, T. Chantrait<sup>1,3</sup>, C. Danquigny<sup>1,4</sup>, M. Dovale Álvarez<sup>1,5</sup>, B. Fang<sup>1,3</sup>, A. Freise<sup>5</sup>, R. Geiger<sup>1,3</sup>, J. Gillot<sup>1,2</sup>, S. Henry<sup>6</sup>, J. Hinderer<sup>1,7</sup>, D. Holleville<sup>1,3</sup>, J. Junca<sup>1,2</sup>, G. Lefèvre<sup>1,2</sup>, M. Merzougui<sup>1,8</sup>, N. Mielec<sup>1,3</sup>, T. Monfret<sup>1,9</sup>, S. Pelisson<sup>1,2</sup>, M. Prevedelli<sup>10</sup>, S. Reynaud<sup>1,11</sup>, I. Riou<sup>1,2</sup>, Y. Rogister<sup>1,7</sup>, S. Rosat<sup>1,7</sup>, E. Cormier<sup>1,12</sup>, A. Landragin<sup>1,3</sup>, W. Chaibi<sup>1,8</sup>, S. Gaffet<sup>1,9,13</sup> & P. Bouyer<sup>1,2</sup>

We present the MIGA experiment, an underground long baseline atom interferometer to study gravity at large scale. The hybrid atom-laser antenna will use several atom interferometers simultaneously interrogated by the resonant mode of an optical cavity. The instrument will be a demonstrator for gravitational wave detection in a frequency band (100 mHz–1 Hz) not explored by classical ground and space-based observatories, and interesting for potential astrophysical sources. In the initial instrument configuration, standard atom interferometry techniques will be adopted, which will bring to a peak strain sensitivity of $2 \cdot 10^{-33}/\sqrt{Hz}$ at 2 Hz. This demonstrator will enable to study the techniques to push further the sensitivity for the future development of gravitational wave detectors based on large scale atom interferometers. The experiment will be realized at the underground facility of the Laboratoire Souterrain à Bas Bruit (LSBB) in Rustrel–France, an exceptional site located away from major anthropogenic disturbances and showing very low background noise. In the following, we present the measurement principle of an in-cavity atom interferometer, derive the method for Gravitational Wave signal extraction from the antenna and determine the expected strain sensitivity. We then detail the functioning of the different systems of the antenna and describe the properties of the installation site.

After its demonstration in 1991<sup>1–4</sup>, atom interferometry (AI) rapidly revealed its huge potential for the precise measurement of inertial forces<sup>5–7</sup>. The following development of cold atom techniques lead to a large number of compact, stable and accurate experiments investigating both fundamental and applied aspects. Indeed, atom interferometry rapidly found a large range of applications such as measurement of fundamental constants<sup>8–11</sup>, gravimetry<sup>12–14</sup>, gradiometry<sup>15</sup>, underground survey<sup>16,17</sup> and inertial navigation<sup>18,19</sup>. The relentless progress in terms of sensitivity of AI based instruments now enables to investigate fundamental questions like those related to the interface of gravity and quantum mechanics<sup>20,21</sup>. The ability of AI to reveal extremely small changes in the inertial field can be used for General Relativity (GR) tests such as the validation of the Weak Equivalence Principle<sup>22–27</sup> and the local Lorentz invariance of post-Newtonian gravity<sup>28</sup>. Such experiments would particularly profit from micro gravity environments allowing extremely long interaction times, and therefore several GR

Received: 25 May 2018

Accepted: 8 August 2018

Published online: 14 September 2018

<sup>1</sup>MIGA Consortium, Talence, France. <sup>2</sup>LP2N, Laboratoire Photonique, Numérique et Nanosciences, Université Bordeaux-IOGS-CNRS:UMR 5298, rue F. Mitterrand, F-33400, Talence, France. <sup>3</sup>LNE-SYRTE, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, 61, avenue de l'Observatoire, F-75014, Paris, France. <sup>4</sup>UMR 1114 EMMAH, Université d'Avignon et des Pays de Vaucluse, INRA, BP 21239, F-84916, Avignon Cedex 9, France. <sup>5</sup>School of Physics and Astronomy and Institute of Gravitational Wave Astronomy, University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK. <sup>6</sup>Oxford University, Department of Physics, Denys Wilkinson Building, Keble Road, Oxford, OX1 3RH, UK. <sup>7</sup>Institut de Physique du Globe de Strasbourg, UMR 7516, Université de Strasbourg/EOST, CNRS, 5 rue Descartes, 67084, Strasbourg, France. <sup>8</sup>Laboratoire ARTEMIS, Université Côte d'Azur, CNRS, Observatoire Côte d'Azur, Bd de l'Observatoire, F-06304, Nice cedex 4, France. <sup>9</sup>Géoazur, Université Côte d'Azur, IRD, CNRS, Observatoire de la Côte d'Azur, 250 rue Albert Einstein, Sophia Antipolis, 06560, Valbonne, France. <sup>10</sup>Dipartimento di Fisica e Astronomia, Università di Bologna, Via Berti-Pichat 6/2, I-40126, Bologna, Italy. <sup>11</sup>Laboratoire Kastler Brossel, CNRS, Sorbonne Université, ENS-PSL Université, Collège de France, Campus Pierre et Marie Curie, F-75252, Paris, France. <sup>12</sup>CELIA, Centre Lasers Intenses et Applications, Université Bordeaux, CNRS, CEA, UMR 5107, F-33405, Talence, France. <sup>13</sup>LSBB, Laboratoire Souterrain à Bas Bruit, UNS, UAPV, CNRS:UMS 3538, AMU, La Grande Combe, F-84400, Rustrel, France. B. Canuel and A. Bertoldi contributed equally. Correspondence and requests for materials should be addressed to B.C. (email: benjamin.canuel@institutoptique.fr)
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Scheme of the MIGA antenna. A set of three atom interferometers at positions X₁,₂,₃ are interrogated by the resonant field of an optical cavity.</figcaption></figure>

space missions based on AI were also proposed to measure for example the Lens-Thirring effect²⁹ and the Weak Equivalence Principle³⁰⁻³². On the base of these developments, the idea also emerged to use AI for the detection of Gravitational Waves (GWs) at low frequency³³⁻³⁶.

The first direct GW detection by the two giant optical interferometers of Advanced LIGO³⁷ in September 2015 opened a new area for physics: GW detectors will reveal new information about massive astrophysical objects such as neutron stars, black holes, pulsars and their dynamics. The transient signal of this first observation, a chirp in frequency from 35 to 250 Hz lasting about 150 ms with a peak strain amplitude of 10⁻²¹, corresponds to the merging phase of a black hole binary system. Due to their bandwidth limited to the frequency range 10 Hz–10 kHz, only the last evolution phase of binary systems is observable with current GW detectors. Before their coalescence, the same sources emit at lower frequencies quasi-continuous GW signals in their “inspiral” phase. A new class of low frequency detectors would enable to observe the signal of such sources years before they enter in the bandwidth of ground-based optical detectors. For example, this first observed source, GW150914, was emitting at a frequency of 16 mHz 5 years before coalescence with a characteristic strain amplitude of the order of 10⁻²⁰³⁸. Low frequency GW detectors would therefore open the possibility of multi-band GW astronomy with long term observation of all evolution phases of binary systems³⁸. Such observatories would also enable to precisely predict the event of coalescence in time and space³⁹, which would ease coincident observations in the electromagnetic domain⁴⁰. Multi-band GW astronomy offers a great scientific payoff with the perspective of multi-messenger astronomy⁴¹, but also promises new gravity and cosmology tests. Low frequency detectors would therefore open a completely new area for GW astronomy⁴². In this context, an original concept for GW detection appeared that combines AI methods and laser ranging techniques used in optical GW detectors. It consists in measuring by AI the effect of GWs on a laser link correlating distant atomic sensors⁴³⁻⁴⁵. Such approach promises to overcome many limitations that affect purely optical interferometers and opens the way toward sub-Hz GW observation.

This new application gave birth to a new field, *large scale atom interferometry*. The use of AI techniques for the realization of large instruments will drastically increase applications of matter-wave interferometry by offering measurements of space-time fluctuations of the gravitational forces along one direction of space. Such instruments will open groundbreaking perspectives not only for GW detection but also for geosciences, allowing the precise mapping of mass distributions and transfers around the detector. This field will be explored by the Matter wave-laser based Interferometer Gravitation Antenna (MIGA)⁴⁶,⁴⁷, the first very long baseline AI instrument, now under construction. This demonstrator will consist of a network of AI sensors interrogated by the resonant field of an optical cavity, a configuration allowing the interferometric control of the phase front of the manipulation beams⁴⁸, as well as an increase of the instrument sensitivity thanks to the high momentum transferred to the atomic wavefunction⁴⁹. The instrument will be located in dedicated tunnels, 500 m underground at the “Laboratoire Souterrain à Bas Bruit” (LSBB), in Rustrel (France) in an environment with very low background noise and located away from major anthropogenic disturbances. MIGA will offer a multidisciplinary approach for the study of both fundamental and applied aspects of gravitation by correlating the main signal of the experiment with highly-precise instrumentation specific to geoscience. In this article, we will detail the measurement strategy of MIGA, present its general design and analyze the noise contributions and expected sensitivity, in relation to both strain measurements and detection of geophysical phenomena.

## MIGA Measurement Scheme

**Combining atom and laser interferometry.** The working principle of MIGA consists in correlating a set of three distant light pulse atom interferometers using the resonant field of a long ultra-stable optical cavity (see Fig. 1).

Simultaneous matter wave interferometers along the cavity are realized by time-modulation of the laser injected inside the resonator. The signal at the output of each interferometer will depend on the relative phase accumulated along the different interferometric path followed by the matter waves. The atom interferometric phase is linked to the phase of the cavity field, imprinted on the matter waves during each interrogation pulse.

The response of each atom interferometer therefore depends on local inertial forces sensed by the atoms, or any other effects modifying the optical phase seen by the matter waves during the interferometer sequence. Such effects may arise from strain variation inside the optical cavity induced by the effect of GWs or some residual motion of the cavity mirrors. MIGA is therefore an hybrid atom-laser interferometer using an array of atom sensors that will simultaneously measure inertial effects and strain variation inside the optical cavity.

**Geometry of MIGA atom interferometers.** MIGA will make use of 3 pulses Mach-Zehnder Interferometers, interrogated in the Bragg regime by two horizontal cavity fields (see Fig. 2).

After launching on a vertical trajectory, the atoms experience a set of π/2-π-π/2 cavity pulses that separate, deflect and recombine the matter waves. The energy-momentum conservation during the light-matter interaction
<figure><img src="image_4.png" /><figcaption>Figure 2. Principle of a Mach Zehnder interferometer for matter waves.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>Figure 3. An atom source placed at position X along the cavity is interrogated by the counter-propagating fields $E_c^+(X, t)$ and $E_c^-(X, t)$. $\delta x_1(t)$ and $\delta x_2(t)$ are respectively the position fluctuations of the left and right mirrors with respect to the cavity baseline at rest L.</figcaption></figure>

implies to couple atomic states of momenta $|+n\hbar k\rangle$ and $|-n\hbar k\rangle$ where $k$ is the wave vector of the interrogation field and $n$ the Bragg diffraction order. The transition probability $P$ between the different states at the interferometer output is given by a two wave interference formula $P = 1/2[1 - \cos(\Delta\phi_{AT})]$ where $\Delta\phi_{AT}$ is the atom phase shift accumulated inside the interferometer.

Using the sensitivity function formalism⁵⁰, the response of the atom interferometer at time $t$ can be expressed as a function of the variations of the local phase difference $\Delta\varphi(t)$ for the two counter-propagating interrogation fields:

$$ \Delta\phi_{AT}(t) = n \int_{-\infty}^{\infty} s(T-t) \frac{d\Delta\varphi(T)}{dT} dT = n \Delta\varphi(t) * \frac{ds(t)}{dt} \quad (1) $$

where $s(t)$ is the sensitivity function of a three pulse atom interferometer, detailed in⁵⁰, and “*” denotes the convolution.

**Response of an in-cavity atom interferometer.** In the following, we derive the response of an atom interferometer interrogated by a cavity at resonance in presence of strain variations induced by GWs, taking into account the major noise contributions, due to cavity mirror vibrations and frequency noise of the input optical field. We will consider that the GW propagates in the direction perpendicular to the cavity axis $x$. The notations used for the calculation are summarized in Fig. 3. We describe the electromagnetic field in the cavity as a superposition of two counter-propagating waves, $E_c^+(t)$, respectively propagating towards positive and negative $x$. As stated by Eq. 1, the atom interferometer response is determined by the relative phase $\Delta\varphi(t)$ between these fields, which is imprinted on the atoms at position X. We will determine $\Delta\varphi(t)$ as a function of monochromatic variations of the mirror positions $\delta x_{1,2}(t) = \delta x_{1,2} \cos(\Omega_{x_{1,2}}t)$, the frequency noise of the input laser $\delta\nu_{in}(t) = \delta\nu_{in} \cos(\Omega_{\nu_{in}}t)$, and the GW strain amplitude $h(t) = h \cos(\Omega_{gw}t)$.

At steady state, the intra-cavity field propagating in the forward direction at position X can be written in the general form:

$$ E_c^+(t) = E_{c0}(t)e^{i\varphi^+(t)}. $$

At time $t$, photons from the field propagating in the backwards direction at position X, $E_c^-(t)$ are emitted from the forward field at retarded time $t_r = t - \Delta t_r$, where $\Delta t_r$ is the light propagation delay on the round trip between the AI and the cavity end mirror. We then have:
$$E_c^-(t) = E_c^+(t_r) = E_c^+(t - \Delta t_r)$$

which means:

$$\varphi^-(t) = \varphi^+(t - \Delta t_r).$$

To first order in $\Delta t_r$, considering that the phase variations are slow with respect to the cavity round trip, one obtains:

$$\begin{aligned} \varphi^-(t) &= \varphi^+(t) - \Delta t_r \frac{d\varphi^+(t)}{dt} \\ \Delta\varphi(t) &= \Delta t_r \frac{d\varphi^+(t)}{dt} = 2\pi \Delta t_r \nu(t). \end{aligned} \qquad (2)$$

The relative phase imprinted on the atom can be expressed as the product of the propagation delay and the instantaneous frequency of the intracavity field at steady state $\nu(t)$; we calculate both terms thereafter.

**Calculation of $\nu(t)$**. The frequency of the intracavity field $\nu(t)$ can be expressed as the sum of the different contributions:

$$\nu(t) = \nu_0 + \delta\nu_{\text{cin}}(t) + \delta\nu_{\delta x_1}(t) + \delta\nu_{\delta x_2}(t) + \delta\nu_{\text{gw}}(t) \qquad (3)$$

where $\delta\nu_{\text{cin}}(t)$, $\delta\nu_{\delta x_1}(t)$ and $\delta\nu_{\text{gw}}(t)$ accounts for fluctuations of the intra-cavity field frequency induced by variations of the input frequency, vibration of cavity mirrors and GW effects. Calculation of $\delta\nu_{\delta x_1}(t)$ and $\delta\nu_{\text{gw}}(t)$ (to first order in $\frac{\Omega_{x1,2}}{\omega_p}$ and $\frac{\Omega_{gw}}{\omega_p}$) are carried out in annex A and B:

$$\begin{aligned} \delta\nu_{\text{gw}}(t) &= \frac{1}{2} \frac{\Omega_{\text{gw}}}{\omega_p} \nu_0 h \sin(\Omega_{\text{gw}} t) \\ \delta\nu_{\delta x_{1,2}}(t) &= \pm \frac{\Omega_{x_{1,2}}}{\omega_p} \nu_0 \frac{\delta x_{1,2}}{L} \sin(\Omega_{x_{1,2}} t) \end{aligned}$$

where $\omega_p/2\pi = c/(4LF)$ is the frequency pole of the cavity, and $F$ its finesse. Since the cavity acts as a first order low pass filter for fluctuations of the input frequency$^{51}$, we have:

$$\delta\nu_{\text{cin}}(t) = \frac{\delta\nu_{\text{in}} \cos\left[\Omega_{\nu_{\text{in}}} t + \arctan\left(\frac{\Omega_{\nu_{\text{in}}}}{\omega_p}\right)\right]}{\sqrt{1 + \left(\frac{\Omega_{\nu_{\text{in}}}}{\omega_p}\right)^2}}.$$

To first order in $\frac{\Omega_{\nu_{\text{in}}}}{\omega_p}$ we have therefore:

$$\delta\nu_{\text{cin}}(t) \approx \delta\nu_{\text{in}} \cos(\Omega_{\nu_{\text{in}}} t) - \frac{\Omega_{\nu_{\text{in}}}}{\omega_p} \delta\nu_{\text{in}} \sin(\Omega_{\nu_{\text{in}}} t).$$

The instantaneous frequency of the intracavity field can then be expressed:

$$\nu(t) = \nu_0 \left\{ 1 + \frac{\delta\nu_{\text{in}}}{\nu_0} \cos(\Omega_{\nu_{\text{in}}} t) - \frac{\Omega_{\nu_{\text{in}}}}{\omega_p} \frac{\delta\nu_{\text{in}}}{\nu_0} \sin(\Omega_{\nu_{\text{in}}} t) + \frac{\Omega_{x_1}}{\omega_p} \frac{\delta x_1}{L} \sin(\Omega_{x_1} t) - \frac{\Omega_{x_2}}{\omega_p} \frac{\delta x_2}{L} \sin(\Omega_{x_2} t) + \frac{1}{2} \frac{\Omega_{\text{gw}}}{\omega_p} h \sin(\Omega_{\text{gw}} t) \right\} \qquad (4)$$

**Calculation of the propagation delay $\Delta t_r$**. The relativistic invariant for the electromagnetic field propagating in the x direction can be written:

$$ds^2 = c^2 dt^2 - dx^2 + h(t) dx^2 = 0.$$

To first order in $h$:

$$dx = \pm \left[ 1 + \frac{1}{2} h(t) \right] dt \qquad (5)$$
where the plus (resp. minus) sign corresponds to the light propagating along the cavity from left to right (resp. right to left). Photons from the resonant field $E_c^+$ are emitted at position X at retarded time $t_r$ and then reflected on the cavity end mirror at position $L + \delta x_2(t_1)$ at time $t_{r'}$. Integrating Eq. 5 brings to:

$$ (L - X) + \delta x_2(t_{r'}) = c(t_{r'} - t_r) + \frac{c}{2} \int_{t_r}^{t_{r'}} h(\xi) d\xi. \qquad (6) $$

Photons are then reflected back and return to position X at time $t$:

$$ -(L - X) - \delta x_2(t_{r'}) = c(-t + t_{r'}) - \frac{c}{2} \int_{t_{r'}}^{t} h(\xi) d\xi. \qquad (7) $$

Subtracting Eqs 6 and 7, the emission time $t_r$ of the photon can be expressed:

$$ t_r = t - \frac{2(L - X)}{c} - \frac{2\delta x_2(t_{r'})}{c} + \frac{1}{2} \int_{t_r}^{t} h(\xi) d\xi. \qquad (8) $$

This implicit equation can be simplified to first order in $h$ and $\delta x_2$:

$$ t_r = t - \frac{2(L - X)}{c} - \frac{2\delta x_2\left(t - \frac{L-X}{c}\right)}{c} + \frac{1}{2} \int_{t-\frac{2(L-X)}{c}}^{t} h(\xi) d\xi. \qquad (9) $$

$\Delta t_r = t - t_r$ can then be expressed as:

$$ \Delta t_r = \frac{2(L - X)}{c} + \frac{2\delta x_2\left(t - \frac{L-X}{c}\right)}{c} - h\frac{L - X}{c}\text{sinc}\left(\Omega_{gw}\frac{L - X}{c}\right)\text{sinc}\left(\Omega_{gw}\left(t - \frac{L - X}{c}\right)\right). \qquad (10) $$

Expanding in series the cosine terms, we obtain to first order in $\Omega_{gw}L/c$ and $\Omega_{x_2}L/c$:

$$ \Delta t_r = \frac{2(L - X)}{c} + \frac{2\delta x_2}{c} \left( \cos(\Omega_{x_2}t) + \Omega_{x_2} \frac{L - X}{c} \sin(\Omega_{x_2}t) \right) \\ -h \frac{L - X}{c} \left( \cos(\Omega_{gw}t) + \Omega_{gw} \frac{L - X}{c} \sin(\Omega_{gw}t) \right). \qquad (11) $$

Calculation of $\Delta\phi(t)$. From Eq. 2, by keeping only first order terms in $\delta\nu$, $\delta x_{1,2}$ and $h$:

$$ \Delta\phi(t) = 4\pi\nu_0 \frac{(L-X)}{c} \left\{ 1 + \frac{\delta\nu_{in}}{\nu_0} \cos(\Omega_{\nu_{in}}t) - \frac{\Omega_{\nu_{in}}}{\omega_p} \frac{\delta\nu_{in}}{\nu_0} \sin(\Omega_{\nu_{in}}t) - \frac{h}{2} \cos(\Omega_{gw}t) \right. \\ - \frac{h}{2} \Omega_{gw} \left[ \frac{L-X}{c} - \frac{1}{\omega_p} \right] \sin(\Omega_{gw}t) + \frac{\delta x_1}{L} \frac{\Omega_{x_1}}{\omega_p} \sin(\Omega_{x_1}t) \\ \left. + \frac{\delta x_2}{L} \Omega_{x_2} \left[ \frac{L}{c} - \frac{1}{\omega_p} \right] \sin(\Omega_{x_2}t) \right\} + \frac{4\pi\nu_0}{c} \delta x_2 \cos(\Omega_{x_2}t). $$

The terms $\frac{1}{c}$ and $\frac{L}{c}$ are respectively the cavity photon lifetime and half of the cavity transit time; for $F \gg 1$ it holds $\left[\frac{L}{c} - \frac{1}{\omega_p}\right] \approx -\frac{1}{\omega_p}$. We thus obtain:

$$ \Delta\phi(t) = 4\pi\nu_0 \frac{(L-X)}{c} \left\{ 1 + \frac{\delta\nu_{in}}{\nu_0} \cos(\Omega_{\nu_{in}}t) - \frac{h}{2} \cos(\Omega_{gw}t) + \frac{1}{\omega_p} \left[ -\Omega_{\nu_{in}} \frac{\delta\nu_{in}}{\nu_0} \sin(\Omega_{\nu_{in}}t) \right. \right. \\ \left. + \frac{h}{2} \Omega_{gw} \sin(\Omega_{gw}t) + \frac{\delta x_1}{L} \Omega_{x_1} \sin(\Omega_{x_1}t) - \frac{\delta x_2}{L} \Omega_{x_2} \sin(\Omega_{x_2}t) \right] \\ \left. + \frac{4\pi\nu_0}{c} \delta x_2 \cos(\Omega_{x_2}t). \right. $$

Considering time-fluctuating effects with characteristic frequencies smaller than the frequency pole of cavity, we obtain:

$$ \Delta\phi(t) = 4\pi\nu_0 \frac{(L-X)}{c} \left\{ 1 + \frac{\delta\nu(t)}{\nu_0} - \frac{h(t)}{2} + \frac{1}{\omega_p} \left[ \frac{\delta\nu'(t)}{\nu_0} - \frac{h'(t)}{2} + \frac{\delta x'_2(t) - \delta x'_1(t)}{L} \right] \right. \\ \left. + \frac{4\pi\nu_0}{c} \delta x_2(t) \right\} \qquad (12) $$
where character “a” denotes the time derivative. The response of the AI can be expressed from Eq. 1:

$$ \Delta\phi_{AT}(X, t) = n\Delta\phi(X, t) * s'(t) + \Delta\phi_I(X, t) + \alpha(X, t). \quad (13) $$

In this expression, we introduced the detection noise $\alpha(X, t)$ and the inertial signal $\Delta\phi_I(X, t)$ associated with the local gravitational acceleration experienced by the atoms. The latter term reads:

$$ \Delta\phi_I(X, t) = \frac{4n\pi\nu_0}{c}x(X, t) * s'(t) \quad (14) $$

where $x(X, t)$ represents the motion of the atoms along the laser beam direction due to the fluctuations of the local gravity.

The differential signal between two atom interferometers placed in cavity at position $X_1$ and $X_2$ is therefore:

$$ \Delta\phi_{AT}(X_1, t) - \Delta\phi_{AT}(X_2, t) = \frac{4n\pi\nu_0}{c}(X_2 - X_1) \left( \frac{\delta\nu(t)}{\nu_0} - \frac{h(t)}{2} + \frac{x(X_1, t) - x(X_2, t)}{X_2 - X_1} \right. \\ \left. + \frac{1}{\omega_p} \left[ \frac{\delta\nu'(t)}{\nu_0} - \frac{h'(t)}{2} + \frac{\delta x'_2(t) - \delta x'_1(t)}{L} \right] \right) * s'(t) + \alpha(X_1, t) - \alpha(X_2, t). \quad (15) $$

Eq. 15 shows that an in-cavity gradiometer has a response similar to a free space one, with additional terms in $\frac{1}{\omega_p}$ taking into account the response of the cavity to the different noises sources. Moreover, in analogy with the free space configuration, the in-cavity gradiometer presents a strong reduction of the influence of end cavity mirror vibrations with respect to the signal of a single AI.

## Analysis of MIGA Strain Sensitivity

**Equivalent GW strain sensitivity.** In this section we derive the sensitivity of the MIGA instrument to GW strain variation. For sake of simplicity, we will only consider the signal $\Gamma(t)$ obtained from the largest atom gradiometer available from the antenna that will use atom sources close to the cavity mirrors (i.e. $X_1 = 0$ and $X_2 = L$):

$$ \Gamma(t) = \Delta\phi_{AT}(0, t) - \Delta\phi_{AT}(L, t). \quad (16) $$

We calculate the strain sensitivity that can be obtained from this configuration, defined as the minimum detectable GW Power Spectral Density (PSD). Considering that the fluctuations $\delta\nu(t)$, $h(t)$, $\delta x_1(t)$, $\delta x_2(t)$, $\alpha(X_1, t)$, $\alpha(X_2, t)$ and the strain Newtonian Noise $NN(t) = \frac{x(0, t) - x(L, t)}{L}$ are uncorrelated, the PSD $S_Γ$ of the gradiometer signal can be expressed:

$$ S_Γ(\omega) = \left[ \frac{4n\pi\nu_0}{c} L \right]^2 \left( 1 + \frac{\omega^2}{\omega_p^2} \left[ \frac{S_{\delta\nu}(\omega)}{\nu_0^2} + \frac{S_h(\omega)}{4} \right] + \frac{2\omega^2 S_x(\omega)}{\omega_p^2 L^2} + S_{NN}(\omega) \right) |\omega G(\omega)|^2 + 2S_\alpha(\omega) \quad (17) $$

where $G(\omega)$ is the Fourier transform of the AI sensitivity function $s(t)$, and $S_\alpha(\omega)$ denotes the PSD of the fluctuations of $u(t)$. We consider in Eq. 17 that $S_{\delta\nu}(\omega) = S_{\delta x_1}(\omega) = S_{\delta x_2}(\omega)$ and $S_{\alpha(X_1, \dots)}(\omega) = S_{\alpha(X_2, \dots)}(\omega) = S_\alpha(\omega)$.

Using the gradiometer signal, the signal to noise ratio SNR($\omega$) for a GW detection at a frequency $\omega$ is defined by the ratio between the GW term and all other terms of Eq. 17:

$$ \text{SNR}(\omega) = \frac{\left[ 1 + \frac{\omega^2}{\omega_p^2} \right] \frac{S_h(\omega)}{4}}{\left[ 1 + \frac{\omega^2}{\omega_p^2} \right] \frac{S_{\delta\nu}(\omega)}{\nu_0^2} + \frac{2\omega^2 S_x(\omega)}{\omega_p^2 L^2} + S_{NN}(\omega) + \frac{2S_\alpha(\omega)}{\left[ \frac{4n\pi\nu_0}{c} L \right]^2 |\omega G(\omega)|^2}}. \quad (18) $$

The strain sensitivity to GW is then defined by the GW PSD which corresponds to a SNR of 1:

$$ S_h(\omega) = \frac{4S_{\delta\nu}(\omega)}{\nu_0^2} + \frac{4}{\left[1 + \frac{\omega^2}{\omega_p^2}\right]} \left( \frac{2\omega^2 S_x(\omega)}{\omega_p^2 L^2} + S_{NN}(\omega) + \frac{2S_\alpha(\omega)}{\left[\frac{4n\pi\nu_0}{c}L\right]^2 |\omega G(\omega)|^2} \right). \quad (19) $$

**MIGA strain sensitivity curve and influence of instrumental noise.** We now determine the strain sensitivity curve of MIGA from the projection of the different noise sources identified in Eq. 19. The antenna will use a resonator of length $L = 200$ m with an initial Finesse of $F = 100$ that corresponds to a cavity frequency pole of $\frac{\omega_p}{2\pi} = 2.5$ kHz (see Sec. 4.2 for a discussion on cavity parameters). We will consider the use of Rb atom sources interrogated on the $D_2$ line, setting the wavelength of the resonating field to $\frac{c}{\nu_0} = 780$ nm. The total interrogation time of the atom will be $2T = 500$ ms and the Rabi frequency $\frac{\Omega_{eff}}{2\pi} = 10$ kHz.
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Top: average seismic acceleration PSD measured in the LSBB galleries with a STS-2 sensor during a typical quiet period. Bottom: projection of the different noise sources on strain sensitivity.</figcaption></figure>

We consider that the laser at the input of the interrogation cavity will be pre-stabilized with a state-of-the-art reference cavity. The relative frequency noise of the input laser becomes thus limited by the thermal noise of the ultra-stable reference cavity at the level of $S_{\delta\nu}(f) = 0.01 \times (1 \text{ Hz/f})\text{Hz}^2/\text{Hz}^{52-54}$. The projection of such noise on strain sensitivity (first term of Eq. 19) is plotted in blue on Fig. 4-bottom.

The Newtonian Noise is a gravity noise sensed by the atom interferometers of the antenna due to density fluctuations of the medium surrounding the experiment. Such noise was extensively studied in the field of optical GW detectors$^{55,56}$ and the two main sources at low frequencies are coming from seismic and atmospheric perturbations. The third term of Eq. 19 represents the sum of these contributions for the atom sources of MIGA. These contributions were previously studied in$^{47}$ at LSBB, showing a limit of the strain sensitivity smaller than $10^{-16} \text{ Hz}^{-1/2}$ above 0.1 Hz.

The displacement noise $S_x(\omega)$ of the cavity optics that will affect the measurements will depend on the system used for suspending and controlling the cavity. To obtain an upper limit of such noise, we consider the worst case of non-suspended mirrors. Figure 4-top shows the average seismic acceleration PSD measured in the LSBB galleries with a STS-2 sensor during a typical quiet period (see section 5.1 for a discussion on seismic properties at the LSBB laboratory). According to Eq. 19, displacement noise of cavity mirrors can contribute to limit strain sensitivity through phase modulation of the intra-cavity field (second term of Eq. 19). This contribution is plotted in dashed red in Fig. 4-bottom. It must be noted that seismic noise may also affect input optics and create an extra frequency noise on the input laser $S'_{\delta\nu}(\omega) = \frac{\nu_0^2}{c^2} \omega^2 S_{\text{xin}}(\omega)^{57}$, where $S_{\text{xin}}$ is the displacement noise PSD of the input optics. The projection of this extra input laser frequency noise is plotted in solid red in Fig. 4-bottom.

The last limitation to strain sensitivity in Eq. 19 comes from detection noise. At the output of each AI, the atom phase is measured by the fluorescence of the atom clouds. The noise in such process comes from the randomness associated with the quantum projection process, and will be $S_\alpha(\omega) = 1$ (mrad)$^2/\text{Hz}$ for an atom flux of $10^6$ atom/s. The projection of detection noise on strain measurement is plotted in Fig. 4-bottom considering Bragg transition of order $n=1$ and $n=10$ (respectively solid black and solid violet curves).

We observe that in the initial instrument configuration corresponding to $n=1$, the detection noise will be largely dominant on strain measurement in all detector bandwidth, setting an optimum strain sensitivity of $2.2 \cdot 10^{-13}$ at frequency $1/(2T) = 2$ Hz. MIGA will serve as a demonstrator to study techniques to push further the sensitivity by several orders of magnitude for the future realization of gravitational wave detectors based on large scale atom interferometry. Such developments will require sensitivities of the order of $10^{-20}$ at 1 Hz. The
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Overview of the MIGA instrument with the main sub-systems. Three atomic heads at positions X₁,₂,₃ launch atomic clouds in an almost vertical parabolic flight (dotted lines); the atoms are manipulated in the upper part of the parabola with a Bragg interferometric sequence by way of light pulses at 780 nm (red horizontal lines) resonant with two horizontal cavities. The ultra-high-vacuum system encompassing the optical cavities, the mirrors payloads and their stabilization system is represented in gray; the atomic heads are connected to its lower side. The control system of the experimental setup, the laser systems dedicated to each atomic head, and the μ-metal shield enclosing each interferometric region and the related atomic head are not represented in the figure.</figcaption></figure>

strain sensitivity of an atom-gradiometer such as MIGA depends on the phase sensitivity of each atom interferometer, and increases linearly with the baseline size. Among the possibilities to tune the sensitivity of each atom interferometer we list the use of: long interrogation times, by adopting ultra-cold atomic sources<sup>58</sup> or optically guided AI<sup>59</sup>; large momentum splitting techniques<sup>20,60</sup>; improved atom flux; non-classical input states to achieve sub-shot-noise sensitivity<sup>61,62</sup>. Adopting alkali-earth-like atoms coherently manipulated on optical transitions will mitigate the impact of the interrogation laser noise on long baseline interferometers<sup>63</sup>.

## MIGA Antenna Design

For the realization of the atom–laser antenna we rely on robust and well tested technology for the atom interferometers:

* rubidium atoms, the workhorse solution in AI;

* laser cooled atomic sources;

* free falling atom sensors;

* two photon Bragg transitions for the coherent manipulation of the matter waves;

* quantum-projection-noise limited sensitivity.

The unique exception to a standard atom interferometer is represented by the cavity enhanced interrogation (Par. 2.3), an essential component of the atom–laser antenna. For the design of the instrument to be installed underground at LSBB, we will adopt an interrogation time T = 250 ms, and a cavity length of 200 m; these choices translate into an expected strain sensitivity of 2.2·10<sup>-13</sup> at 2 Hz (Fig. 4). Several solutions that go beyond this set of choices as a baseline device will be later taken into account to upgrade the instrument performance.

This size envisaged for the instrument brings AI from a laboratory scale to that of a big infrastructure, with consequent technological and engineering issues in terms of geometry, vacuum, magnetic field etc. In this section we will describe the experimental scheme of the MIGA antenna, shown in Fig. 5. Three atomic clouds are prepared at the same height and horizontally separated by a distance L, using for each a 3D-MOT loaded with a 2D-MOT. The atomic clouds are successively launched in free fall along the vertical direction using a moving optical molasses. Before entering in the interferometric region, the atomic ensembles are prepared in relation to their internal and external degrees of freedom, by means of two sets of horizontal Raman beams and optical pumping techniques. The horizontal acceleration at the position of the three atomic clouds is then simultaneously measured using Bragg interferometry, implemented with a π/2–π–π/2 pulse sequence symmetric with respect to the apogees of the atomic trajectories. After the Bragg pulse sequence, the atomic phase of each atom interferometer encoded as population difference of two atomic states is measured by fluorescence detection, when the clouds pass in the detection regions during their free fall.

**Atomic source.** Each atom interferometer of the MIGA antenna (see Fig. 6) uses clouds of N ≈ 10⁶ atoms at a temperature of a few μK with a repetition rate of 1 Hz. Such sources are based on a 2D-MOT, which loads a 3D-MOT cooled to sub-Doppler temperature and launched vertically in a moving molasses, with a controlled velocity of ≈4 m/s. Before entering the interferometric region, the atomic clouds are prepared on the internal m<sub>F</sub>=0 state, sensitive to magnetic fields only at the second order, and with a horizontal velocity that fulfills the
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Description of the MIGA atom interferometers, from$^{47}$. (a) View of the AI including 2D-MOT, 3D-MOT, preparation and detection systems. A set of 3 different μ-metal shields (red, brown, green colors) are used to screen the system from magnetic fields. (b) Picture of the 2D-MOT, 3D-MOT, preparation and detection systems. (c) Picture of the interior μ-metal shield, placed along the vacuum system of the interrogation cavity (brown shield of (a)).</figcaption></figure>

Bragg condition for diffraction in the interferometric cavities. Two sets of horizontal Raman beams, tuned on the transition between the magnetically insensitive atomic sub-levels, are used to this scope, together with resonant optical beams to remove the atomic population on unwanted levels.

At the apex of their parabolic trajectory, the clouds experience a series of in-cavity $π/2$, $π$, $π/2$ pulses before returning in the detection region where the transition probability is measured. The Bragg interferometer uses the same internal state for both atomic paths, and the output ports are separated only for their horizontal velocity. A Raman pulse is then applied to reflect the horizontal velocity of one of the two output ports, and obtain internal state labeling for the two atomic populations. The counting on the two output ports uses fluorescence detection with atom shot noise sensitivity$^{64,65}$. The normalized population ratio gives the transition probability of the interferometer, and then the atom phase shift.

The preparation phases for the atomic ensembles - cooling, trapping, launch, optical pumping and velocity selection - and the successive fluorescence detection make use of a commercial laser system based on telecom technology and frequency doubling to obtain the required light beams at 780 nm$^{66}$.

**Cavity enhanced interrogation.** Two horizontal cavities of 200 m are used to coherently interrogate the three atomic ensembles simultaneously launched along the vertical direction. A $π/2-π-π/2$ Bragg pulse sequence is applied in a symmetric fashion with respect to the apogee of the atomic parabolas: the matter wave splitting and recombination are implemented through the lower cavity via $π/2$ pulses, the redirection with the upper cavity via a $π$ pulse. The laser radiation to interrogate the atoms is obtained by injecting the fundamental transverse mode TEM$_{00}$ of each cavity, and the Bragg pulses are shaped with acousto-optic modulators (AOMs). The resonance condition of each interrogation beam to the corresponding cavity is granted by generating the 780 nm light via frequency doubling of a master telecom laser at 1560 nm with a Periodically Poled Lithium Niobate (PPLN) crystal$^{67,68}$; the telecom laser is amplified, pre-stabilized to $δν/v_0 = 10^{-15}$ level at 1 Hz using an ultra-stable reference cavity, and finally injected in the two interrogation cavities, as shown in Fig. 7. In order to keep the resonance condition, the length of each cavity is locked at low frequency -outside the detector bandwidth- on the interrogation laser at 1560 nm. The frequency of the pulse shaping AOM at 780 nm is chosen so as to match the cavity resonance condition for the diffracted beam; such frequency is fixed by the differential optical length of the cavity at 1560 nm and 780 nm, because of the effect of the two reflective coatings.

Due to the cavities' long storage time with respect to the duration of the Bragg pulses, the interrogation fields will suffer some degree of deformation of their temporal amplitude profiles$^{69}$. This deformation scales with the cavity finesse, and it can have an adverse impact on the interferometers by increasing the minimum interaction time of the velocity-selective atomic transitions as well as their power requirements$^{70}$. On the other hand, the cavities offer spatial filtering of the interferometric beams due to their frequency-dependent resonance conditions. This filtering, which also scales with the cavity finesse, effectively reduces the sensitivity of the interferometers to laser wavefront distortions, which are a leading source of noise in current state-of-the-art detectors$^{71}$. The finesse of the cavities must therefore be chosen as a trade-off between these two cavity-induced effects on the atom optics pulses.
<figure><img src="image_9.png" /><figcaption>**Figure 7.** Scheme of the laser interrogation system, from$^{47}$ A master laser at 1560 nm (represented in green) is amplified and locked to a reference cavity, before being frequency doubled to obtain the interrogation radiation at 780 nm (represented in red). The master laser is frequency locked to one of the two interrogation cavities, whereas the component at 780 nm is pulsed by means of an AOM.</figcaption></figure>

For the 200 m MIGA cavities a finesse of 100 is chosen to strike a balance between the two effects for Bragg orders $n < 10^{70}$. The resulting cavity bandwidth yields a minimum interaction time 1.5 times larger than in the absence of the cavity for $n = 10$, but below $1/\omega_r$ for all $n$ (with $\omega_r = \hbar k^2/2M$ the 2-photon recoil frequency and $M$ the mass of the atom). In this configuration, most higher order spatial modes will be optically suppressed at the $10^{-2}$ level. A lower finesse would lead to less dilation of the interaction time and slightly improved power enhancement of the short beam splitter pulses, especially for the higher-order diffraction processes, but would incur a worse optical suppression of higher order modes, partially negating the benefit of the cavity-assisted atom optics.

The radii of curvature (ROC) of the cavity mirrors are chosen to yield a beam waist sufficient to interrogate the atomic clouds as they thermally expand during the measurement. Furthermore, the resulting cavity configuration provides suppression of mode degeneracy for Hermite-Gauss modes up to order 15, taking into account manufacturing tolerances of the ROC. For the chosen ROC of 555 m, the g-factor of the resonator is 0.64 and the beam radius is 7.28 mm at the waist and 8.04 mm at the mirrors, providing a good margin over the cloud radius at the last interferometric pulse, which is expected to be ≈5 mm for a 1 µK $^{87}$Rb ensemble at the end of the 2T = 500 ms pulse sequence. The resulting cavity configuration is robust to ROC deviations, although a reliable alignment sensing and control system will have to be implemented.

The two cavities share a common payload on each side to hold the mirrors, placed at a vertical distance of ≈30.6 cm to have an interrogation time of $T = 250$ ms. The impact of ground seismic noise on the position of the cavity mirrors will be reduced by means of an anti-vibration system. The interest of this system is two-fold: it limits the related phase noise contribution on each atom interferometer to a negligible level with respect to the atomic shot noise contribution and insures that the cavity remains close to resonance inside the detection bandwidth. We are working on two different approaches: a passive system of mechanical filters to suspend each payload, or an active stabilization of each mirror position using piezoelectric actuators. The main constraints on the anti-vibration system are set by the seismic noise level at the chosen location (see Par. 5.1), and the response function of the atom interferometers to mirror acceleration noise (see Par. 2.3). Furthermore, relative length fluctuations of the two cavities must be controlled, since they would determine a shift of the readout atomic phase.

The systematic shift on the interferometric measurement due to the Coriolis effect is differentially canceled in the gradiometric measurement, if the launch directions for the atomic ensembles do not change between the experimental runs. The effect of a change between two sensors in the relative launch direction of an angle $\Delta\theta$ in the horizontal direction perpendicular to $k$ will induce a Coriolis contribution to the differential phase$^{72}$

$$ \phi_C|_{n=1} = -2\Omega_E \sin\alpha_{lat} kT^2 v_l \Delta\theta \approx 470 \cdot \Delta\theta \quad (20) $$

where $\Omega_E = 72.9\mu rad/s$ is the Earth's rotation rate, $\alpha_{lat}$ the latitude equal to 43.9° at LSBB. For a systematic effect below the atomic shot noise, equal to 1 mrad in terms of the interferometric differential phase, a launch stability at 1 $\mu rad$ level is required for each sensor.

**Environment and control system.** The operation of each atom interferometer and of the cavity enhanced optical link implementing the atomic interrogation sets stringent environmental requirements. The cold atomic ensembles must be manipulated in regions evacuated to high vacuum in order to limit collisions with the thermal particles in the residual atmosphere. Differential vacuum techniques are adopted to obtain the vacuum levels required by different experimental regions: the typical pressure in the 2D-MOT cell is ≈10⁻⁸ mbar, optimized to achieve a good atomic flux to load the 3D-MOT; the 3D-MOT, the detection, and the interferometer regions are operated at a pressure below 10⁻⁹ mbar to have a long lifetime of the atomic samples and reduce systematic effects related to the presence of background pressure. The vacuum level along the optical link results less demanding than what required at the interferometric region, if the phase front distortions caused by the residual index of refraction are taken into account. Nevertheless, a similar requirement is finally obtained when one considers the difficulty to implement differential vacuum between the respective regions to limit optical diffraction of the Bragg beams into the cavity and the related decoherence effect. The interferometric interrogation time $T = 250$ ms translates into a vertical distance of ~30.6 cm between the optical axis of the two cavities; a single vacuum tube with a diameter $\Phi = 0.5$ m is chosen to host them, since a configuration with two separated tubes with a reduced diameter and conductance would make challenging the evacuation of the system to the required level. The target
vacuum pressure in the different experimental regions will be obtained using vibration free pumping systems, relying on a combination of ion pumps and non-evaporable getter pumps.

The magnetic susceptibility of Rb<sup>73</sup> requires magnetic shielding of a factor >10<sup>4</sup> for the regions where the atoms are prepared and coherently manipulated<sup>74–76</sup>. As can be seen in Fig. 5, 3 different shields will screen each atom interferometer: one is dedicated to the atomic source and detection region, and two to the interferometric region.

The experimental sequence, timing and data acquisition are computer controlled. The synchronization of the three atomic sensors is granted by the use of common interrogation beams, which is the same element on which relies the high rejection ratio of common mode noise in the gradiometric configuration. On the other hand, the optimization of the three atomic signals requires a good synchronization of the preparation and launching sequence for each of the three atomic cloud, which will be obtained by referencing and precisely phasing the operation of each atomic head to a common time signal.

## MIGA Installation Site: The LSBB Underground Science Platform

MIGA will be installed at the LSBB, an underground low-noise laboratory located in Rustrel, near the city of Apt in Vaucluse, France (see Fig. 8). The LSBB is a European interdisciplinary laboratory for science and technology created in 1997 from the decommissioning of a launching control system of the French strategic nuclear defense operative during the Cold War. The LSBB is now an underground scientific platform<sup>77</sup> characterized by an ultra-low noise environment, both seismic and electromagnetic, as a result of the distance of the site from heavy industrial and human activities. The LSBB fosters trans-disciplinary interactions and interdisciplinary approaches, pursuing both fundamental and applied research. The result is a broad scientific and industrial expertise<sup>78</sup>, besides the openness to European and international research programs and groups. The LSBB is an ideal facility for site studies of next generation GW detectors and more generally for the improvement of low-frequency sensitivity of existing and future GW antennas. In such environment, MIGA aims at studying Newtonian Noise and testing advanced detector geometries for its cancellation<sup>57</sup>.

The LSBB is embedded in the Fontaine-de-Vaucluse watershed, which is one of the world's largest karst aquifers<sup>79</sup>, and covers 54 ha in surface area. The underground facility includes 4 km-long horizontal drifts at a depth ranging from 0 m to 518 m, with north to south and north-east to south-west as its main orientations (Fig. 8). The galleries give access to underground wells, vaults and voids, where a passive temperature stability better than 0.1 °C results from the thermal blanketing supplied by the upper carbonate rock layer. Air pressure and circulation are controlled in the galleries by airlocks. At the deepest point, 518 m under the surface, is located the 1,250 m³ vault previously equipped as nuclear launch control room, hence designed to remain operational even in the case of a nearby nuclear blast. This place is nowadays used for the most demanding experiments in terms of seismic and electromagnetic residual noise. The whole system of underground drifts and surface areas is connected to energy, telephony, optical fiber Internet, and GPS synchronization. The LSBB's environment provides benefits for developing ultra-sensitive instrumentation, calibration and comparison of highly sensitive sensors, optimization of experimental protocols and the operation of industrial prototypes. The LSBB takes advantages of its remote location and its robust design and provides access to the karstified carbonate reservoir through galleries that may be assimilated to kilometric horizontal boreholes. The site provides an original access to the karst medium, complementary to usually investigated caves and springs<sup>80</sup>. Therefore, it is of broad interest for hydrogeologists and geophysicists who study karst systems<sup>81–83</sup>. Such a configuration helps also to develop research programs on poro-elastic dynamics of fractured media, induced seismicity and processes for internal damaging of rock massifs, and the handling of fluids and gas in reservoir<sup>84</sup>. A broad range high sensitivity instruments (e.g. optical strainmeters, a hydrostatic long baseline tiltmeter, seismometers, superconducting magnetometers and gravimeters, muon cameras for rock density, cold atom gravimeters) is studied and developed at LSBB; their metrological performances are characterized and cross-compared. The interdisciplinary ability of the LSBB allows to take into account the effects induced by the disturbances of the physical environment of the sensor on the measure it produces, like the rock mass tilts<sup>83</sup> or influence of gravity gradient variations on deformation measurements<sup>85,86</sup>.

**LSBB low noise properties.** The underground environment of LSBB is characterized by exceptional low-noise properties, which is the basic condition to realize high sensitivity experiments in diverse fields, like for example the first observation of simultaneously recorded seismic and seismo-magnetic signals caused by an earthquake<sup>87</sup>, the acquisition of extremely weak biological signals<sup>88</sup>, and the operation of a portable atomic gravimeter at high sensitivity without an anti-vibration platform<sup>89</sup>.

A network of 6 broadband seismometers is deployed at LSBB to monitor the seismic noise variations at the site<sup>78</sup>. Sample acceleration noise signals recorded at the station RUSF.01, located underground in a very quiet place, are presented in Fig. 9: their probability density function (PDF) is compared to the Peterson's models<sup>90</sup>, commonly used as a reference for the definition of the quality of a seismic recording site. The seismic data for three orthogonal components are reported for three different measurement intervals: the whole 2011 year, which includes the Mw 9.1 Tohoku-Oki March 11, 2011 mega-thrust Japanese earthquake; a 6 hour interval during which the same event took place; a 24 hours interval during 2011 without seismic events. The graphs show that the site has remarkable seismic properties, and, except during important transient signals, the highest probability of noise occurrence over a long interval is close to the new low-noise Peterson's model for all the three components and for the whole considered frequency band (Fig. 9-left). At frequencies below 0.2 Hz, the one-year PDF is broader and spreads between the high and low noise models of Peterson, since it includes the worldwide seismic activity and hence surface waves with common frequencies below 1/20 Hz and high amplitudes; among all events the strongest one has been the Tohoku-Oki March 11 earthquake, for which high seismic energy was measured at low frequency (Fig. 9-center). The noise PDF over a quiet day gives a sharp high probability that at low frequencies borders (horizontal components) or even drops below (vertical component) the low noise model of Peterson (Fig. 9-right).
<figure><img src="image_10.png" /><figcaption>**Figure 8.** (A) The LSBB is placed in south-east France. (B) Topographic map with the galleries (thick black lines) and permanent broadband seismic stations installed at the surface (red rectangle) and in the galleries (blue rectangles). (C) On the top image the 4-km LSBB galleries (cyan lines) and the length of its different branches (red). On the bottom images, from left to right: view of the top of the mountain with the shielded antenna dome intended to receive the launching order for the nuclear rockets; the launching control room in a EM shielded volume; the access to the galleries.</figcaption></figure>

The low-noise seismic properties of LSBB have been confirmed also with gravity measurements realized with a high sensitivity superconducting gravimeter (model iOSG from GWR Instruments Inc., Fig. 10-left), installed at the underground site in 2015 as a complementary instrument to the MIGA experiment. Figure 10-right shows a recent measurement of the noise level at LSBB, compared with the curves obtained in Strasbourg (France) and at the Black Forest Observatory (Germany). The instrument installed at LSBB has a noise performance among the best in a worldwide network of superconducting gravimeters⁹¹,⁹², with an amplitude spectral density of 1.8 nm·s⁻²·Hz⁻¹/² at 1 mHz. Notably, superconducting gravimeters have a characteristic noise lower than the Earth's background noise at frequencies below 1 mHz⁹³, and permit thus to study low-frequency seismic and sub-seismic modes.

The specific location of LSBB beneath a massif determines a unique sheltering effect as concerning electromagnetic noise and its effect on the measurement of the atom interferometric phase⁹⁴. The laboratory lies in a very quiet environment, the Regional Natural Park of Luberon: minimal human activities within two kilometers reduce the magnetic interference from railways and high-voltage power lines. Despite being in a moderately seismic area, the seismic noise spectra recorded at the site are close to the worldwide minima, hence there is no significant movement of magnetic mass, like the content of the nearby water reservoir, which could perturb the magnetic background noise.

An exceptionally low noise electromagnetic environment is provided by the former control room of the military facility: such infrastructure was designed to withstand the effects of a nearby nuclear blast, including the electromagnetic pulse (EMP) in its deepest zone. The room consists of an unconventional Faraday cage, where the electromagnetic shielding effect is provided both by karstic 500 m thick rocks loaded with water, which gives
<figure><img src="image_11.png" /><figcaption>**Figure 9.** Probability density function of seismic ambient noise (color scale) of RUSF.01 broadband station for year 2011, including quiet days and days with earthquakes (left); 6 hours of ground motion which include the Tohoku-Oki earthquake of March 11, 2011 (center); a quiet day (24 hours) of year 2011 without seismic events (right). The seismic noise PSD for three components (top, Z; middle, NS; bottom, EW) is compared to the Peterson's high and low noise models (black lines).</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 10.** (left) The superconducting gravimeter installed at LSBB in the frame of MIGA experiment. (right) Comparison of the noise power spectral density for superconducting gravimeters at different world sites: LSBB at Rustrel - France (“RU”, red); Strasbourg - France (“ST”, blue); Black Forest Observatory near Schiltach - Germany (“BFO”, green). The signals are obtained from daily power spectral densities on raw data sampled at 1 second. The sharp drop at high frequency (>0.1 Hz) is due to the anti-aliasing filters present in the superconducting gravimeters. The Peterson's low noise model (NLMN in solid gray) and the seismological GSN median noise model of Berger¹²³ (GSN median in dashed-dotted gray) are plotted for reference.</figcaption></figure>

a high-frequency cut-off around 200 Hz, and by the 2 m thick reinforced concrete capsule with a 1 cm thick steel inner coating, which reduces the cut-off to 10 Hz. The absence of μ-metal makes it a perfect low frequency pass-band filter, which is different from a zero Gauss chamber. The performance of the shielding was measured with a 3-axis Superconductive Quantum Interference Device (SQUID) magnetometer. The residual noise level is about 2 fT/√Hz above 40 Hz⁹⁵, and the attenuation in the DC domain leads to a residual field inside the Capsule of about 6 μT, compared to the 46 μT expected at the LSBB latitude. Inside the vault is placed a cabin measuring
<figure><img src="image_13.png" /><figcaption>**Figure 11.** From<sup>124</sup>. Frequency spectra of magnetic fluctuations measured at four different locations in the LSBB complex, indicated on the map in Fig. 12: the Capsule (CAP), the Anti-blast Gallery (GAS), the Secondary Gallery (GGB), and the Safety Gallery (C).</figcaption></figure>

20 × 6 m² and weighting 30 t, which hangs from the ceiling and rests on eight shock absorbers ensuring the decoupling from the ground movement. This location is used to perform experiments that require the most demanding quality in terms of seismic and electromagnetic environment.

Figure 11 shows frequency spectra recorded using a portable SQUID magnetometer in the Capsule, and at three other points inside the tunnel complex. In the range 0.01–10 Hz of interest to MIGA, all curves show a similar level till 60 mHz, except at the Capsule where the steel walls attenuate the geomagnetic fluctuations. Above 0.1 Hz, the magnitude of the fluctuations varies with the screening of the surrounding rock and any interference from nearby equipment. The spectra recorded at point C show broad peaks at the Schumann resonance frequencies. Other peaks can be attributed to mechanical vibrations or electromagnetic interference (50 Hz power lines and its harmonics are clear). At higher frequencies, the instrument noise starts to dominate. The signals are screened by the telluric currents induced in surround rock by the changing external field. This creates a significant gradient in these magnetic signals across the laboratory. A search for magnetic field changes correlated with large changes in groundwater flow has ruled out signals above 0.2 nT<sup>96</sup>.

**MIGA infrastructure at LSBB.** The initial approach of installing MIGA in the “Main Gallery” at LSBB (see Fig. 12) has been discarded, in relation to the 5% slope of the tunnel, which will introduce the projection of the gravitational acceleration into the measurement, and also to the requirement to have a quiet dedicated site and not a shared environment used to reach other experiments. Two horizontal galleries will be instead bored to host the MIGA instrument, as shown in Fig. 12. The choice to have two orthogonal tunnels is related to the possibility to mitigate the impact of the interrogation laser technical noise on strain measurements by adopting a Michelson-Morley configuration, in the same way as LIGO and VIRGO do. Nevertheless, mainly for budgetary reasons, only one tunnel will be initially equipped with an atom-laser antenna. The second tunnel will be equipped later to realize a 2D gravitational antenna. The orientation and position of the tunnels are chosen so as to exploit the existing galleries for the required access points at their extremities for safety reasons. This configuration forbids the use of drilling machines for the excavation, which will be realized using explosive charges. The explosions will be exploited by geoscientists at LSBB to study the response of the karstic mountain to man-made seismic activity. Boring the two 200 m long, 3.2 m wide tunnels will require a 12–18 months period of time, depending on the quality of the rocks, and an estimated cost of about 4 M€.

A preliminary design of the infrastructure is presented in Fig. 13; the size of the two tunnels and of the auxiliary rooms required for the cavity injection and for the atom interferometers is defined in terms of the instrument encumbrance, safety and environmental requirements. The tunnels will have to host the vacuum system used by the optical link, the AIs, the cavity injection optical setups and various electronic and data acquisition equipments. Apart from the large cavities at each end of the galleries, widenings will be present every 50 m along the tunnel to have the possibility to change the distance of the three initial AIs, and that of increasing the number of AIs so as to improve the spatial resolution of the experiment.

**Geophysics at LSBB with MIGA.** MIGA will extend the concept of correlated interferometry from the laboratory scale to that of a geological site, the LSBB at Rustrel, using an underground array of atom sensors distributed along a 200 m horizontal arm. Several techniques based on correlated atom interferometry will be implemented to characterize the gravitational field of the site, such as the simultaneous measurement of gravity acceleration and gradient<sup>97</sup> and the measurement of gravitational curvature<sup>20,98</sup>. It will be thus possible to investigate several geological phenomena, like the non-invasive detection of underground density anomalies<sup>17</sup>, the gravity perturbations due to local density changes caused by fault evolution as proposed in<sup>45</sup>, and the characterization of gravity-gradient noise, also called Newtonian Noise (NN)<sup>47,57</sup>.

The MIGA differential signal is sensitive to gravity gradient, related to the term $x(X_1) - x(X_2)$ of Eq. 15. The fluctuation of this term determines the NN, a tidal effect that affects any couple of test masses, both macroscopic and microscopic, and is considered to be a fundamental limit for ground based GW detectors. Nevertheless, GW and NN have different length scales over the detector's dimensions: GW signals have extremely long characteristic
<figure><img src="image_14.png" /><figcaption>**Figure 12.** (Above) Map of the underground galleries at LSBB, with the locations adopted for the characterization of magnetic fluctuations at the site indicated by red points. The place where MIGA will be installed is highlighted with a dashed blue rectangle. (Below) Zoom of the MIGA infrastructure: the two orthogonal boreholes will use the main gallery at their far ends for access, whereas the anti-blast gallery will be used as a safety exit. The three AIs will be located in the room used to inject the cavities (a), at the other end of one gallery (c) and at its mid position (b).</figcaption></figure>

lengths and are seen as pure gradients at the distances of interest, whereas NN has shorter characteristic lengths going from the meter to a few kilometers⁵⁵. Effects due to NN may become discernible using a network of correlated sensors distributed along the antenna's direction⁵⁷. The MIGA instrument will help to model and
<figure><img src="image_15.png" /><figcaption>**Figure 13.** Design of the galleries dedicated to the MIGA experiment at LSBB. The three atom interferometers of the antenna will be located at (a), (b) and (c). The optical setups for cavity injection will be hosted in room (a). The two MIGA galleries will be equipped with enlargements regularly spaced in order to add other AIs in the future.</figcaption></figure>

characterize the space-time correlation properties of the NN in an underground environment. Such study could pave the way to novel rejection methods for this kind of disturbance applicable to future sub-Hz ground based GW detectors.

MIGA will provide absolute gravity and gravity-gradient measurements, which will be used to obtain density maps for the surrounding volume via inversion algorithms. The limited resolution of the measurement, due to the reduced number of atomic sensors along the cavity arm, their fixed position and directionality, will be mitigated by correlating the AI measurements with other kinds of gravity measurements provided by the instrument network deployed at LSBB. This approach will also implement a long baseline hybrid gravity sensor, broadening the concept introduced in⁹⁹. The device will be calibrated on the signal produced by heavy source masses¹⁰⁰,¹⁰¹ placed at a variable distance up to 150 m, exploiting the nearby tunnels already present at the site. The instrument will provide a gravity gradient sensitivity of the order of $10^{-3}$ Eötvös at 1 s ($1 \text{ E} = 10^{-9} \text{ s}^{-2}$), an order of magnitude below the reported sensitivity of superconducting devices over much smaller probe distances¹⁰². The absolute readout of the AI-based device would finally allow to periodically calibrate the other geophysical instruments operated at LSBB.

In the following phase, the underground structure of the surrounding karstic environment will be studied, to detect subsurface cavities¹⁰³ or monitor groundwater dynamics. Indeed, better characterization of complex underground reservoirs is expected from recent and future developments of geophysical methods¹⁰⁴. Their application to karst is probably the most challenging¹⁰⁵ because karst heterogeneity is multi-scale, self and hierarchical organized and induces the most complex underground fluid dynamics. Peculiarly, the lack of non-invasive methods producing multi-scale 4D imaging of underground fluids remains a bottleneck for understanding and modeling this dynamics. One of the important questions is to have enough resolution and depth of investigation all at once to detail all the features controlling the groundwater circulation and storage from matrix porosity or micro-fracturing to major faults and karst conduits. Currently, only integrating methods directly or indirectly related to water content such as seismic, Electrical Resistivity Tomography (ERT), Magnetic Resonance Sounding (MRS) or gravimetry allow estimating the variation of water mass. Since no method presents at the same time the required resolution, depth of investigation and fluids sensitivity, coupling and comparisons appear to be a promising way for imaging underground structure and fluids dynamics (e.g.⁸²,¹⁰⁶). The place of gravimetry surveys is thus increasing in hydrological studies⁸⁵,¹⁰⁷⁻¹⁰⁹. Combining conventional instruments and methods from hydrogeology with cold atom based gravity measurements will allow better understanding and modeling of karst aquifers¹¹⁰, for which 4D data are currently lacking to constrain spatially distributed models¹¹¹.

With MIGA we also plan to apply aperture synthesis imaging techniques to geophysics, adopting well developed techniques from radio-astronomy¹¹² and optics¹¹³. To this aim, the distance between the atomic sensors installed on the optical link will be changed during successive measurement runs, so as to probe different Fourier components of the gravity field.

## Conclusions and Outlook

We presented the underground MIGA atom-laser antenna that is being built at LSBB in Rustrel to measure space-time strain in the infrasound with an expected peak sensitivity of $2 \cdot 10^{-13}/\sqrt{\text{Hz}}$ at 2 Hz. The experiment is presented in what will be its environment, characterized by an exceptionally low seismic and electromagnetic noise, and by the presence of a network of geophysical instruments monitoring the site. The setup will represent a demonstrator for gravitational wave detection using atom interferometry, in a frequency band not explored by classical ground and space-based observatories and will also permit to measure geophysical phenomena with unprecedented sensitivity. The underground location of the setup as well as its size pose several technological and scientific challenges, which are addressed in the article together with the devised solutions. However, this location offers a great opportunity to confront MIGA measurements with other gravity monitoring sensors and to explore innovative interdisciplinary applications of atom interferometry, especially in geoscience.

After the first phase, described in this article, the MIGA instrument will be used to test several advanced techniques to improve strain sensitivity of large baseline atom-gradiometers that could include the use of higher atom flux, large momentum beam splitters⁶⁰,¹¹⁴,¹¹⁵, non-classical input states⁶¹,⁶², asynchronous interrogation of
<figure><img src="image_16.png" /><figcaption>Figure 14. GW strain sensitivity curves of the “AI array” proposed in$^{57}$ and those of eLISA, aLIGO and ET (data from$^{125}$). We also plot the strength of different low frequency GW sources (from$^{125}$) and the one corresponding to the first GW detection, GW150914$^{37}$.</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>Figure 15. A laser field $A^+$ is emitted at position 0 and reflected at distance L. $A^-$ is the returning field.</figcaption></figure>

many ensembles$^{116}$, increase of the interrogation time via optically guided interferometry$^{59}$ or intermediate coherence preserving interrogations and feedback$^{117}$. These studies are needed for the future realization of sub-Hz GW detectors that requires sensitivity improvement by several orders of magnitude. In such perspective, specific metrological studies that can be done with MIGA are also required to determine for example the impact of wave-front errors$^{118}$ or diffused light, but also to realize improved NN rejection schemes$^{57}$.

The MIGA demonstrator could then pave the way towards the realization of a low frequency GW detector such as the large “AI array” described in$^{57}$. This configuration considers the use of large momentum splitting, improved detection noise and NN rejection scheme to obtain sensitivities of about $3 \cdot 10^{-23}/\sqrt{\text{Hz}}$ at 2 Hz. As seen in Fig. 14, the strain sensitivity of this configuration would complement those of space or ground detectors such as eLISA$^{119}$, aLIGO$^{120}$ or ET$^{121}$ and would offer observation of different low frequency GW sources.

## Annex A

**Calculation of the intracavity field frequency Modulation Under The Influence of a Monochromatic GW.** We first consider the effect of a monochromatic GW $h(t) = h\cos(\Omega_{gw}t)$ on a single round-trip of distance 2L. We consider that a laser field $A^+(t) = A_0e^{-i\omega t}$ is emitted at position 0 (see Fig. 15) and reflected at distance L. When it comes back to point with position 0, the field $A^-(t)$ can be expressed:

$$A^{-}(t) = A^{+}(t_{r}) \qquad (21)$$

where $t_r$ is the retarded time calculated in Sec. 2.3 (see Eq. 9) which can be expressed to first order in $\Omega_{gw}L/c$:

$$t_r = t - \frac{2L}{c} + h \frac{L}{c} \cos(\Omega_{gw}(t - L/c)). \qquad (22)$$

$A^-(t)$ results to be:

$$A^{-}(t) = A_{0}e^{-i\omega(t - \frac{2L}{c})}e^{-i\hbar\frac{\omega L}{c}\cos(\Omega_{gw}(t - L/c))}. \qquad (23)$$

Expanding the second exponential term in $h$, one obtains to first order:

$$A^{-}(t) = A_{0}e^{-i\omega(t - \frac{2L}{c})}\left(1 - i\hbar\frac{\omega L}{c}\cos(\Omega_{gw}(t - L/c))\right) \qquad (24)$$

and then:
$$A^{-}(t) = A_{0}e^{-i\omega(t-\frac{2L}{c})} - \frac{1}{2}i\hbar A_{0}\frac{\omega L}{c}e^{2i\frac{\omega L}{c}e^{-i\frac{\Omega_{gw}L}{c}e^{-i(\omega+\Omega_{gw})t}} - \frac{1}{2}i\hbar A_{0}\frac{\omega L}{c}e^{2i\frac{\omega L}{c}e^{-i\frac{\Omega_{gw}L}{c}e^{-i(\omega-\Omega_{gw})t}}. \quad (25)$$

After one round trip, the effect of the GW is to introduce two sidebands at frequencies $(\omega \pm \Omega_{gw})/2\pi$ of amplitude proportional to h. We calculate the effect of a GW on the circulating field of a cavity following closely$^{122}$ in Sec. 1.3. Going back and forth many times inside the cavity, the intracavity field will remain to first order in h of the form:

$$A(t) = \left( A_0 + \frac{1}{2}hA_1e^{-i\Omega_{gw}t} + \frac{1}{2}hA_2e^{i\Omega_{gw}t} \right) e^{-i\omega t} \quad (26)$$

which enables to defines a “generalized amplitude” for the cavity field as a rank 3 vector A:

$$A = (A_0; A_1; A_2). \quad (27)$$

We consider in the following a field $A_{in} = (A_0; 0; 0)$ incident on a cavity of length L formed with two mirrors with coefficients of reflection and transmission $(r, t_c)$ (see Fig. 16). The amplitude of the cavity field after one round trip A' can then be expressed as a matrix product:

$$A' = -r^2 X A \quad (28)$$

where X is the linear round trip operator:

$$X = e^{2iu} \begin{pmatrix} 1 & 0 & 0 \\ -i u \operatorname{sinc}(v) e^{i v} & e^{2i v} & 0 \\ -i u \operatorname{sinc}(v) e^{-i v} & 0 & e^{-2i v} \end{pmatrix} \quad (29)$$

where $u = \frac{\omega L}{c}$ and $v = \frac{\Omega_{gw}L}{c}$. Considering only the first order in v of sinc(v), one obtains:

$$X = e^{2iu} \begin{pmatrix} 1 & 0 & 0 \\ -i u e^{i v} & e^{2i v} & 0 \\ -i u e^{-i v} & 0 & e^{-2i v} \end{pmatrix}. \quad (30)$$

At steady state, the interference of the intracavity waves on the first mirror can be written (see Fig. 16):

$$A = t_c A_{in} - r^2 X A \quad (31)$$

which gives:

$$A = [1 + r^2 X]^{-1} t_c A_{in}. \quad (32)$$

After Inversion of the $1 + r^2 X$ Matrix, we calculate the generalized amplitude of the resonating field:

$$A = A_0 \left[ \frac{t_c}{1 + r^2 e^{2iu}}, \frac{ir^2 t_c u e^{i(2u+v)}}{(1 + r^2 e^{2iu})(1 + r^2 e^{2i(u+v)}),} \frac{ir^2 t_c u e^{i(2u+v)}}{(1 + r^2 e^{2iu})(r^2 e^{2iu} + e^{2iv})} \right]. \quad (33)$$

The cavity being at resonance, we have $e^{2iu} = -1$. From Eq. 26, A(t) can then be written:

$$A(t) = A_0 \left[ \frac{t_c}{1 - r^2} - \frac{ihr^2 t_c u e^{i t \Omega_{gw} + i v}}{2(1 - r^2)(-r^2 + e^{2iv})} - \frac{ihr^2 t_c u e^{i v - i t \Omega_{gw}}}{2(1 - r^2)(1 - r^2 e^{2iv})} \right] e^{i \omega t}. \quad (34)$$

The resonating field can then be expressed in the form $A(t) = A_0 e^{i(\omega t + \phi(t))}$ where:

$$\tan \phi(t) \approx \phi(t) = \frac{hr^2 u (r^2 \cos(v + t\Omega_{gw}) - \cos(v - t\Omega_{gw}))}{r^4 - 2r^2 \cos(2v) + 1}. \quad (35)$$

To first order in v:

$$\phi(t) = h u - \frac{r^2}{r^2 - 1} \left( \cos(t \Omega_{gw}) - \frac{(r^2 + 1)}{r^2 - 1} v \sin(t \Omega_{gw}) \right). \quad (36)$$

Substituting $\frac{r^2+1}{r^2-1} \approx -\frac{2F}{\pi}$ and $\frac{r^2}{r^2-1} \approx -\frac{F}{\pi}$ where F is the Finesse of the cavity, we obtain:

$$\phi(t) = -\frac{FL}{\pi c} \omega h \left( \cos(t\Omega_{gw}) + \frac{2F\Omega_{gw}L}{\pi c} \sin(t\Omega_{gw}) \right). \quad (37)$$

The frequency noise $\delta\nu_{gw}(t)$ induced by the GW on the resonating field can then be expressed as:

<footnote>$E_{in}$, $E_{ref}$ and $E_{trans}$ are respectively the incident, reflected and transmitted amplitudes on a cavity mirror; we obtain $E_{ref} = e^{in^2/2}rE_{in}$ and $E_{trans} = tE_{in}$.</footnote>
<figure><img src="image_18.png" /><figcaption>Figure 16. A laser field A<sub>in</sub> is incident on a cavity of length L formed with two mirrors with coefficients of reflection and transmission (r, t<sub>c</sub>). A is the circulating field and A' its copy after one round trip.</figcaption></figure>

<figure><img src="image_19.png" /><figcaption>Figure 17. A laser field A<sup>+</sup> is emitted at position 0 and reflected at distance L on a vibrating mirror. A<sup>-</sup> is the returning field.</figcaption></figure>

$$ \delta\nu_{gw}(t) = \frac{1}{2\pi} \frac{d\phi(t)}{dt} = \frac{F\Omega_{gw}L}{\pi c} \nu_0 h \left( \sin(t\Omega_{gw}) - \frac{2F\Omega_{gw}L}{\pi c} \cos(t\Omega_{gw}) \right) \quad (38) $$

$$ \delta\nu_{gw}(t) = \frac{1}{2} \frac{\Omega_{gw}}{\omega_p} \nu_0 h \left( \sin(t\Omega_{gw}) - \frac{\Omega_{gw}}{\omega_p} \cos(t\Omega_{gw}) \right) \quad (39) $$

where $\omega_p/2\pi = c/(4LF)$ is the frequency pole of the cavity. For $\Omega_{gw} \ll \omega_p$ we approximate to first order:

$$ \delta\nu_{gw}(t) = \frac{1}{2} \frac{\Omega_{gw}}{\omega_p} \nu_0 h \sin(t\Omega_{gw}). \quad (40) $$

## Annex B

**Calculation of the Intracavity Field Frequency Modulation Under The Influence of a Monochromatic Cavity Mirror Vibration.** As in the previous Annex, we first consider the effect of a monochromatic mirror vibration $\delta x(t) = \delta x \cos \Omega_x t$ on a single round-trip of distance $2L$: a laser field $A^+(t) = A_0 e^{-i\omega t}$ is emitted at position 0 (see Fig. 17) and reflected at distance L on a vibrating mirror. When its comes back to point position 0, the field $A^-(t)$ can be expressed:

$$ A^{-}(t) = A^{+}(t_{r}) \quad (41) $$

where $t_r$ is the retarded time introduced in Sec. 2.3 (see Eq. 9):

$$ t_r = t - \frac{2L}{c} - \frac{2\delta x}{c} \cos(\Omega_x(t - L/c)). \quad (42) $$

The effect of a mirror vibration of amplitude $\delta x$ on a round trip is therefore completely similar to the one of a GW of amplitude $h = -\frac{2\delta x}{c}$ (see Eqs 42 and 22).

If we consider a cavity formed by two mirrors placed at position $x_1, x_2$ (see Fig. 3), the effect of the vibration of a cavity mirror on the frequency noise $\delta\nu_x(t)$ of the resonating field can then be determined in analogy with Annex A. We obtain:
$$ \delta\nu_{\delta x_{1,2}}(t) = \pm \frac{\Omega_{x_{1,2}}}{\omega_p} \nu_0 \frac{\delta x_{1,2}}{L} \sin(\Omega_{x_{1,2}} t). \quad (43) $$

The effect of both cavity mirrors on the frequency noise of the circulating field is anti-symmetric. This can be simply explained by the fact that the effect on the accumulated phase on a single round trip is opposite for the two mirrors.

## Data Availability

The datasets studied during the current study are available from the corresponding author on reasonable request.

## References

1. Carnal, O. & Mlynek, J. Young's double-slit experiment with atoms: A simple atom interferometer. Phys. Rev. Lett. 66, 2689, https://doi.org/10.1103/PhysRevLett.66.2689 (1991).

2. Keith, D. W., Ekstrom, C. R., Turchette, Q. A. & Pritchard, D. E. An interferometer for atoms. Phys. Rev. Lett. 66, 2693, https://doi.org/10.1103/PhysRevLett.66.2693 (1991).

3. Riehle, F., Kisters, T., Witte, A., Helmcke, J. & Bordé, C. J. Optical Ramsey spectroscopy in a rotating frame: Sagnac effect in a matter-wave interferometer. Phys. Rev. Lett. 67, 177, https://doi.org/10.1103/PhysRevLett.67.177 (1991).

4. Kasevich, M. & Chu, S. Atomic interferometry using stimulated Raman transitions. Phys. Rev. Lett. 67, 181, https://doi.org/10.1103/PhysRevLett.67.181 (1991).

5. Gustavson, L., Landragin, A. & Kasevich, M. A. Rotation sensing with a dual atom-interferometer Sagnac gyroscope. Classical Quant. Grav. 17, 2385, https://doi.org/10.1088/0264-9381/17/12/311 (2000).

6. Peters, A., Chung, K. Y. & Chu, S. High-precision gravity measurements using atom interferometry. Metrologia 38, 25, https://doi.org/10.1088/0026-1394/38/1/4 (2001).

7. Canuel, B. et al. Six-axis inertial sensor using cold-atom interferometry. Phys. Rev. Lett. 97, 010402, https://doi.org/10.1103/PhysRevLett.97.010402 (2006).

8. Fixler, J. B., Foster, G. T., McGuirk, J. M. & Kasevich, M. A. Atom interferometer measurement of the Newtonian constant of gravity. Science 315, 74, https://doi.org/10.1126/science.1135459 (2007).

9. Rosi, G., Sorrentino, F., Cacciapuoti, L., Prevedelli, M. & Tino, G. Precision measurement of the Newtonian gravitational constant using cold atoms. Nature 510, 518, https://doi.org/10.1038/nature13433 (2014).

10. Bouchendira, R., Cladé, P., Guellati-Khélifa, S., Nez, F. & Biraben, F. New determination of the fine structure constant and test of the quantum electrodynamics. Phys. Rev. Lett. 106, 080801, https://doi.org/10.1103/PhysRevLett.106.080801 (2011).

11. Parker, R. H., Yu, C., Zhong, W., Estey, B. & Müller, H. Measurement of the fine-structure constant as a test of the standard model. Science 360, 191-195, https://doi.org/10.1126/science.aap7706 (2018).

12. Peters, A., Chung, K. Y. & Chu, S. Measurement of gravitational acceleration by dropping atoms. Nature 400, 849, https://doi.org/10.1038/23655 (1999).

13. Gillot, P., Francis, O., Landragin, A., Santos, F. P. D. & Merlet, S. Stability comparison of two absolute gravimeters: optical versus atomic interferometers. Metrologia 51, L15-L17, https://doi.org/10.1088/0026-1394/51/5/L15 (2014).

14. Freier, C. et al. Mobile quantum gravity sensor with unprecedented stability. J. Phys. Conf. Ser. 723, 012050, https://doi.org/10.1088/1742-6596/723/1/012050 (2016).

15. Snadden, M. J., McGuirk, J. M., Bouyer, P., Haritos, K. G. & Kasevich, M. A. Measurement of the earth's gravity gradient with an atom interferometer-based gravity gradiometer. Phys. Rev. Lett. 81, 971, https://doi.org/10.1103/PhysRevLett.81.971 (1998).

16. de Angelis, M. et al. Precision gravimetry with atomic sensors. Meas. Sci. Technol. 20, 022001, https://doi.org/10.1088/0957-0233/20/2/022001 (2009).

17. Metje, N., Chapman, D. N., Rogers, C. D. F. & Bongs, K. Seeing through the ground: the potential of gravity gradient as a complementary technology. Adv. Civ. Eng. 903758, https://doi.org/10.1155/2011/903758 (2011).

18. Kasevich, M. A. & Dubetsky, B. United States Patent 7317184 (2008).

19. Battelier, B. et al. Development of compact cold-atom sensors for inertial navigation. In Stuhler, J. & Shields, A. J. (eds) Quantum Optics, https://doi.org/10.1117/12.2228351 (SPIE, 2016).

20. Asenbaum, P. et al. Phase shift in an atom interferometer due to spacetime curvature across its wave function. Phys. Rev. Lett. 118, 183602, https://doi.org/10.1103/physrevlett.118.183602 (2017).

21. Geiger, R. & Trupke, M. Proposal for a quantum test of the weak equivalence principle with entangled atomic species. Phys. Rev. Lett. 120, 043602, https://doi.org/10.1103/PhysRevLett.120.043602 (2018).

22. Dimopoulos, S., Graham, P., Hogan, J. & Kasevich, M. Testing general relativity with atom interferometry. Phys. Rev. Lett. 98, 111102, https://doi.org/10.1103/PhysRevLett.98.111102 (2007).

23. Geiger, R. et al. Detecting inertial effects with airborne matter-wave interferometry. Nat. Commun. 2, 474, https://doi.org/10.1038/ncomms1479 (2011).

24. Schlippert, D. et al. Quantum test of the universality of free fall. Phys. Rev. Lett. 112, 203002, https://doi.org/10.1103/PhysRevLett.112.203002 (2014).

25. Barrett, B. et al. Correlative methods for dual-species quantum tests of the weak equivalence principle. New J. Phys. 17, 085010, https://doi.org/10.1088/1367-2630/17/8/085010 (2015).

26. Barrett, B., Bertoldi, A. & Bouyer, P. Inertial quantum sensors using light and matter. Phys. Scripta 91, 053006, https://doi.org/10.1088/0031-8949/91/5/053006 (2016).

27. Rosi, G. et al. Quantum test of the equivalence principle for atoms in coherent superposition of internal energy states. Nat. Commun. 8, 15529, https://doi.org/10.1038/ncomms15529 (2017).

28. Müller, H., Chiow, S.-W., Herrmann, S., Chu, S. & Chung, K.-Y. Atom-interferometry tests of the isotropy of post-Newtonian gravity. Phys. Rev. Lett. 100, 031101, https://doi.org/10.1103/PhysRevLett.100.031101 (2008).

29. Jentsch, C., Müller, T., Rasel, E. M. & Ertmer, W. HYPER: A satellite mission in fundamental physics based on high precision atom interferometry. Gen. Rel. Gravit. 36, 2197, https://doi.org/10.1023/b:gerg.0000046179.26175.fa (2004).

30. Aguilera, D. et al. STE-QUEST test of the universality of free fall using cold atom interferometry. Classical Quant. Grav. 31, 115010, https://doi.org/10.1088/0264-9381/31/11/115010 (2014).

31. Schubert, C. et al. Differential atom interferometry with $^{87}$Rb and $^{85}$Rb for testing the WEP in STE-QUEST. arXiv:1312.5963 [physics.atom-ph] (2013).

32. Tino, G. M. et al. Precision gravity tests with atom interferometry in space. Nucl. Phys. B - Proc. Suppl. 243-244, 203-217, https://doi.org/10.1016/j.nuclphysbps.2013.09.023 (2013).

33. Bordé, C. J., Sharma, J., Tourrenc, P. & Damour, T. Theoretical approaches to laser spectroscopy in the presence of gravitational fields. J. Physique 44, 983-990, https://doi.org/10.1051/jphyslet:019830044024098300 (1983).
34. Chiao, R. Y. & Speliotopoulos, A. D. Towards MIGO, the matter-wave interferometric gravitational-wave observatory, and the intersection of quantum mechanics with general relativity. *J. Mod. Opt.* **51**, 861, https://doi.org/10.1080/09500340408233603 (2004).

35. Foffa, S., Gasparini, A., Papucci, M. & Sturani, R. Sensitivity of a small matter-wave interferometer to gravitational waves. *Phys. Rev. D* **73**, 022001, https://doi.org/10.1103/PhysRevD.73.022001 (2006).

36. Tino, G. M. & Vetrano, F. Is it possible to detect gravitational waves with atom interferometers? *Classical Quant. Grav.* **24**, 2167, https://doi.org/10.1088/0264-9381/24/9/001 (2007).

37. LIGO Scientific Collaboration and Virgo Collaboration. Observation of gravitational waves from a binary black hole merger. *Phys. Rev. Lett.* **116**, 061102, https://doi.org/10.1103/PhysRevLett.116.061102 (2016).

38. Sesana, A. Prospects for multiband gravitational-wave astronomy after GW150914. *Phys. Rev. Lett.* **116**, 231102, https://doi.org/10.1103/PhysRevLett.116.231102 (2016).

39. Graham, P. W. & Jung, S. Localizing gravitational wave sources with single-baseline atom interferometers. *Phys. Rev. D* **97**, 024052 (2018).

40. Abbott, B. P. et al. Multi-messenger observations of a binary neutron star merger. *Astrophys. J* **848**, L12, https://doi.org/10.3847/2041-8213/aa91c9 (2017).

41. Christensen, N. L. for the LIGO Scientific Collaboration & the Virgo Collaboration. Multimessenger Astronomy. *ArXiv e-prints*, 1105.5843 (2011).

42. Mandel, I., Sesana, A. & Vecchio, A. The astrophysical science case for a decihertz gravitational-wave detector. *Classical Quant. Grav.* **35**, 054004, https://doi.org/10.1088/1361-6382/aaa7e0 (2018).

43. Dimopoulos, S., Graham, P. W., Hogan, J. M., Kasevich, M. A. & Rajendran, S. Gravitational wave detection with atom interferometry. *Phys. Lett. B* **678**, 37, https://doi.org/10.1016/j.physletb.2009.06.011 (2009).

44. Dimopoulos, S., Graham, P., Hogan, J., Kasevich, M. & Rajendran, S. Atomic gravitational wave interferometric sensor. *Phys. Rev. D* **78**, 122002, https://doi.org/10.1103/PhysRevD.78.122002 (2008).

45. Harms, J. et al. Low-frequency terrestrial gravitational-wave detectors. *Phys. Rev. D* **88**, 122003, https://doi.org/10.1103/PhysRevD.88.122003 (2013).

46. Canuel, B. et al. The matter-wave laser interferometer gravitation antenna (MIGA): New perspectives for fundamental physics and geosciences. *E3S Web of Conf.* **4**, 01004, https://doi.org/10.1051/e3sconf/20140401004 (2014).

47. Canuel, B. et al. MIGA: combining laser and matter wave interferometry for mass distribution monitoring and advanced geodesy. *Proc. SPIE* **9900**, 990008, https://doi.org/10.1117/12.22288251 (2016).

48. Hamilton, P. et al. Atom interferometry in an optical cavity. *Phys. Rev. Lett.* **114**, 100405, https://doi.org/10.1103/PhysRevLett.114.100405 (2015).

49. Riou, I. et al. A marginally stable optical resonator for enhanced atom interferometry. *J. Phys. B* **50**, 155002, https://doi.org/10.1088/1361-6455/aa7592 (2017).

50. Cheinet, P. et al. Measurement of the sensitivity function in a time-domain atomic interferometer. *IEEE T. Instrum. Meas.* **57**, 1141, https://doi.org/10.1109/TIM.2007.915148 (2008).

51. Zhu, M. & Hall, J. L. Stabilization of optical phase/frequency of a laser system: application to a commercial dye laser with an external stabilizer. *J. Opt. Soc. Am. B* **10**, 802, https://doi.org/10.1364/josab.10.000802 (1993).

52. Numata, K., Kemery, A. & Camp, J. Thermal-noise limit in the frequency stabilization of lasers with rigid cavities. *Phys. Rev. Lett.* **93**, 250602, https://doi.org/10.1103/PhysRevLett.93.250602 (2004).

53. Thorpe, M. J., Rippe, L., Fortier, T. M., Kirchner, M. S. & Rosenband, T. Frequency stabilization to $6 \times 10^{-16}$ via spectral-hole burning. *Nat. Photon.* **5**, 688–693, https://doi.org/10.1038/nphoton.2011.215 (2011).

54. Jiang, Y. Y. et al. Making optical atomic clocks more stable with $10^{-16}$-level laser stabilization. *Nat. Photon.* **5**, 158–161, https://doi.org/10.1038/nphoton.2010.313 (2011).

55. Saulson, P. Terrestrial gravitational noise on a gravitational wave antenna. *Phys. Rev. D* **30**, 732–736, https://doi.org/10.1103/PhysRevD.30.732 (1984).

56. Harms, J. Terrestrial gravity fluctuations. *Living Rev. Relativ.* **18**, 3, https://doi.org/10.1007/lrr-2015-3 (2015).

57. Chaibi, W. et al. Low frequency gravitational wave detection with ground-based atom interferometer arrays. *Phys. Rev. D* **93**, 021101, https://doi.org/10.1103/PhysRevD.93.021101 (2016).

58. Abend, S. et al. Atom-chip fountain gravimeter. *Phys. Rev. Lett.* **117**, https://doi.org/10.1103/physrevlett.117.203003 (2016).

59. McDonald, G. et al. Optically guided linear Mach-Zehnder atom interferometer. *Phys. Rev. A* **87**, 013632, https://doi.org/10.1103/PhysRevA.87.013632 (2013).

60. McDonald, G. D. et al. 80hk momentum separation with Bloch oscillations in an optically guided atom interferometer. *Phys. Rev. A* **88**, 053620, https://doi.org/10.1103/PhysRevA.88.053620 (2013).

61. Hosten, O., Engelsen, N. J., Krishnakumar, R. & Kasevich, M. A. Measurement noise 100 times lower than the quantum-projection limit using entangled atoms. *Nature* **529**, 505–508, https://doi.org/10.1038/nature16176 (2016).

62. Cox, K. C., Greve, G. P., Weiner, J. M. & Thompson, J. K. Deterministic squeezed states with collective measurements and feedback. *Phys. Rev. Lett.* **116**, https://doi.org/10.1103/physrevlett.116.093602 (2016).

63. Graham, P. W., Hogan, J. M., Kasevich, M. A. & Rajendran, S. New method for gravitational wave detection with atomic sensors. *Phys. Rev. Lett.* **110**, 171102, https://doi.org/10.1103/PhysRevLett.110.171102 (2013).

64. Santarelli, G. et al. Quantum projection noise in an atomic fountain: A high stability cesium frequency standard. *Phys. Rev. Lett.* **82**, 4619–4622, https://doi.org/10.1103/physrevlett.82.4619 (1999).

65. Rocco, E. et al. Fluorescence detection at the atom shot noise limit for atom interferometry. *New J. Phys.* **16**, 093046, https://doi.org/10.1088/1367-2630/16/9/093046 (2014).

66. http://www.muquans.com/ (2018).

67. Chiow, S.-w., Kovachy, T., Hogan, J. M. & Kasevich, M. A. Generation of 43 W of quasi-continuous 780 nm laser light via high-efficiency, single-pass frequency doubling in periodically poled lithium niobate crystals. *Opt. Lett.* **37**, 3861, https://doi.org/10.1364/OL.37.003861 (2012).

68. Sané, S. S. et al. 11 W narrow linewidth laser source at 780 nm for laser cooling and manipulation of rubidium. *Opt. Expr.* **20**, 8915, https://doi.org/10.1364/oe.20.008915 (2012).

69. Fang, B. et al. Improving the phase response of an atom interferometer by means of temporal pulse shaping. *New J. Phys.* **20**, 023020, https://doi.org/10.1088/1367-2630/aaa37c (2018).

70. Dovale-Alvarez, M. et al. Fundamental limitations of cavity-assisted atom interferometry. *Phys. Rev. A* **96**, 053820–10, https://doi.org/10.1103/PhysRevA.96.053820 (2017).

71. Kovachy, T. et al. Quantum superposition at the half-metre scale. *Nature* **528**, 530–533, https://doi.org/10.1038/nature16155 (2015).

72. Bertoldi, A. et al. Atom interferometry gravity-gradiometer for the determination of the Newtonian gravitational constant G. *Eur. Phys. J. D* **40**, 271–279, https://doi.org/10.1140/epjd/e2006-00212-2 (2006).

73. Steck, D. A. Rubidium 87 D line data, http://steck.us/alkalidata/rubidium87numbers.1.6.pdf (2001).

74. Dickerson, S. et al. A high-performance magnetic shield with large length-to-diameter ratio. *Rev. Sci. Instrum.* **83**, 065108, https://doi.org/10.1063/1.4720943 (2012).
75. Schuldt, T. et al. Design of a dual species atom interferometer for space. *Exp. Astronom.* 1–40, https://doi.org/10.1007/s10686-014-9433-y (2015).

76. Kubelka-Lange, A. et al. A three-layer magnetic shielding for the MAIUS-1 mission on a sounding rocket. *Rev. Sci. Instrum.* **87**, 063101, https://doi.org/10.1063/1.4952586 (2016).

77. Bettini, A. The world deep underground laboratories. *Eur. Phys. J.-Plus.* **127**, 114, https://doi.org/10.1140/epjp/i2012-12114-y (2012).

78. Gaffet, S. et al. A 3D broadband seismic array at LSBB. *IRIS Data Services Newsletter* **11**, 3 http://ds.iris.edu/ds/newsletter/vol11/no3/a-3d-broadband-seismic-array-at-lsbb/ (2009).

79. Ford, D. & Williams, P. *Karst Hydrogeology and Geomorphology* https://doi.org/10.1002/9781118684986 (John Wiley & Sons Ltd, 2007).

80. Garry, B. et al. Contribution of artificial galleries to the knowledge of karstic system behaviour in addition to natural cavern data. *Int. J. Spel.* **37**, 75, https://doi.org/10.5038/1827-806X.37.1.7 (2008).

81. Network of hydrogeological research sites http://hplus.ore.fr/ (2018).

82. Carrière, S. D. et al. The role of porous matrix in water flow regulation within a karst unsaturated zone: an integrated hydrogeophysical approach. *Hydrogeol. J.* **24**, 1905, https://doi.org/10.1007/s10040-016-1425-8 (2016).

83. Lesparre, N. et al. New insights on fractures deformation from tiltmeter data measured inside the fontaine de vaucluse karst system. *Geophys. J. Int.* **208**, 1389, https://doi.org/10.1093/gji/ggw446 (2016).

84. Wang, J., Guglielmi, Y. & Gaffet, S. Collaborative projects between two USA-France national subsurface laboratories to improve imaging of fractured-porous rocks properties and coupled THMCB processes. In Dudt & Mathier (eds) *Rock Mechanics in Civil and Environmental Engineering*, 857–860 (Taylor and Francis Group, London, 2010).

85. Deville, S., Jacob, T., Chery, J. & Champollion, C. On the impact of topography and building mask on time varying gravity due to local hydrology. *Geophys. J. Int.* **192**, 82, https://doi.org/10.1093/gji/ggs007 (2012).

86. Fores, B., Champollion, C., Moigne, N. L., Bayer, R. & Chéry, J. Assessing the precision of the iGrav superconducting gravimeter for hydrological models and karstic hydrological process identification. *Geophys. J. Int.* **208**, 269–280, https://doi.org/10.1093/gji/ggw396 (2016).

87. Gaffet, S. et al. Simultaneous seismic and magnetic measurements in the Low-Noise Underground Laboratory (LSBB) of Rustrel, France, during the 2001 january 26 Indian earthquake. *Geophys. J. Int.* **155**, 981, https://doi.org/10.1111/j.1365-246X.2003.02095.x (2003).

88. Zandi, A. S. et al. Scalp EEG acquisition in a low-noise environment: A quantitative assessment. *IEEE Trans. Biomed. Eng* **58**, 2407–2417, https://doi.org/10.1109/tbme.2011.2158647 (2011).

89. Farah, T. et al. Underground operation at best sensitivity of the mobile lne-syrrte cold atom gravimeter. *Gyroscopy Navig.* **5**, 266, https://doi.org/10.1134/S2075108714040051 (2014).

90. Peterson, J. Observations and modelling of seismic background noise. http://earthquake.usgs.gov/regional/asl/pubs/files/ofr93-322.pdf, U.S. Geol. Surv. Open-File Rept. 93-332, Albuquerque, New Mexico (1993).

91. Rosat, S. et al. First analyses of the iOSG-type superconducting gravimeter at the low noise underground laboratory (LSBB URL) of Rustrel, France. *E3S Web of Conf.* **12**, 06003, https://doi.org/10.1051/e3sconf/20161206003 (2016).

92. Rosat, S. et al. A two-year analysis of the iOSG-24 superconducting gravimeter at the low noise underground laboratory (LSBB URL) of Rustrel, France: Environmental noise estimate. *J. Geodyn.* **119**, 1–8, https://doi.org/10.1016/j.jog.2018.05.009 (2018).

93. Rosat, S. & Hinderer, J. Noise levels of superconducting gravimeters: Updated comparison and time stability. *Bull. Seism. Soc. Am.* **101**, 1233, https://doi.org/10.1785/0120100217 (2011).

94. Canuel, B. Etude d'un gyromètre à atomes froids, https://tel.archives-ouvertes.fr/tel-00193288. PhD Thesis, Université Paris XI (2007).

95. Waysand, G. et al. First characterization of the ultra-shielded chamber in the low-noise underground laboratory (LSBB) of Rustrel Pays d'Apt. *Instrum. Meth. A* **444**, 336, https://doi.org/10.1016/S0168-9002(99)01377-7 (2000).

96. Henry, S. et al. Monitoring geomagnetic signals of groundwater movement using multiple underground SQUID magnetometers. *E3S Web of Conferences* **4**, 02004, https://doi.org/10.1051/e3sconf/20140402004 (2014).

97. Sorrentino, F. et al. Simultaneous measurement of gravity acceleration and gravity gradient with an atom interferometer. *Appl. Phys. Lett.* **101**, 114106, https://doi.org/10.1063/1.4751112 (2012).

98. Rosi, G. et al. Measurement of the gravity field curvature by atom interferometry. *Phys. Rev. Lett.* **114**, 013001, https://doi.org/10.1103/PhysRevLett.114.013001 (2015).

99. Lautier, J. et al. Hybridizing matter-wave and classical accelerometers. *Appl. Phys. Lett.* **A05**, 144102, https://doi.org/10.1063/1.4897358 (2014).

100. Gray, S. D., Parmentola, J. A. & LeSchack, R. Estimating the weight of very heavy objects with a gravity gradiometer. *J. Phys. D Appl. Phys.* **28**, 2378, https://doi.org/10.1088/0022-3727/28/11/024 (1995).

101. Lamporesi, G. et al. Source mass and positioning system for an accurate measurement of g. *Rev. Sci. Instrum.* **78**, 075109, https://doi.org/10.1063/1.2751090 (2007).

102. Moody, M. V., Paik, H. J. & Canavan, E. R. Three-axis superconducting gravity gradiometer for sensitive gravity experiments. *Rev. Sci. Instrum.* **73**, 3957–3974, https://doi.org/10.1063/1.1511798 (2002).

103. Butler, D. K. Microgravimetric and gravity gradient techniques for detection of subsurface cavities. *Geophys.* **49**, 1084, https://doi.org/10.1190/1.1441723 (1984).

104. Berkowitz, B. Characterizing flow and transport in fractured geological media: A review. *Adv. Water Resour.* **25**, 861, https://doi.org/10.1016/s0309-1708(02)00042-8 (2002).

105. Chalikakis, K., Plagnes, V., Guerin, R., Valois, R. & Bosch, F. P. Contribution of geophysical methods to karst-system exploration: an overview. *Hydrogeol. J.* **19**, 1169–1180, https://doi.org/10.1007/s10040-011-0746-x (2011).

106. Carrière, S. D., Chalikakis, K., Sénéchal, G., Danquigny, C. & Emblanch, C. Combining electrical resistivity tomography and ground penetrating radar to study geological structuring of karst unsaturated zone. *J. Appl. Geophys.* **94**, 31, https://doi.org/10.1016/j.jappgeo.2013.03.014 (2013).

107. Hector, B. et al. Hydro-gravimetry in West-Africa: First results from the Djougou (Benin) superconducting gravimeter. *J. Geodyn.* **80**, 34, https://doi.org/10.1016/j.jog.2014.04.003 (2014).

108. Hasan, S., Troch, P. A., Bogaart, P. W. & Kroner, C. Evaluating catchment-scale hydrological modeling by means of terrestrial gravity observations. *Water Resour. Res* **44**, https://doi.org/10.1029/2007wr006321 (2008).

109. Pool, D. R. & Eychaner, J. H. Measurements of aquifer-storage change and specific yield using gravity surveys. *Ground Water* **33**, 425, https://doi.org/10.1111/j.1745-6584.1995.tb00299.x (1995).

110. Geiger, R. et al. Matter-wave laser interferometric gravitation antenna (MIGA): New perspectives for fundamental physics and geosciences. In Augé, E., Dumarchez, J. & Trân Thanh Văn, J. (eds.) *Proceedings of the 50th Rencontres de Moriond “100 years after GR”, La Thuile (Italy)*, 21–28 March 2015, http://arxiv.org/abs/1505.07137 (ARSIF, 2015).

111. Ghasemizadeh, R. et al. Review: Groundwater flow and transport modeling of karst aquifers, with particular reference to the North Coast Limestone aquifer system of Puerto Rico. *Hydrogeol. J.* **20**, 1441, https://doi.org/10.1007/s10040-012-0897-4 (2012).

112. Thompson, A. R., Moran, J. M. & Swenson, G. W. *Interferometry and Synthesis in Radio Astronomy, Second Edition* https://doi.org/10.1002/9783527617845 (John Wiley & Sons, inc., 2007).

113. Greenaway, A. H. Optical aperture synthesis. *Meas. Sci. Technol.* **2**, 1, https://doi.org/10.1088/0957-0233/2/1/001 (1991).
114. Chiow, S.-W., Kovachy, T., Chien, H.-C. & Kasevich, M. A. 102hk large area atom interferometers. *Phys. Rev. Lett.* **107**, 130403, https://doi.org/10.1103/PhysRevLett.107.130403 (2011).

115. Estey, B., Yu, C., Müller, H., Kuan, P.-C. & Lan, S.-Y. High-resolution atom interferometers with suppressed diffraction phases. *Phys. Rev. Lett.* **115**, 083002, https://doi.org/10.1103/PhysRevLett.115.083002 (2015).

116. Dutta, I. et al. Continuous cold-atom inertial sensor with 1 nrad/sec rotation stability. *Phys. Rev. Lett.* **116**, https://doi.org/10.1103/physrevlett.116.183003 (2016).

117. Kohlhaas, R. et al. Phase locking a clock oscillator to a coherent atomic ensemble. *Phys. Rev. X* **5**, 021011, https://doi.org/10.1103/PhysRevX.5.021011 (2015).

118. Louchet-Chauvet, A. et al. The influence of transverse motion within an atomic gravimeter. *New J. Phys.* **13**, 065025, https://doi.org/10.1088/1367-2630/13/6/065025 (2011).

119. Amaro-Seoane, P. et al. eLISA: Astrophysics and cosmology in the millihertz regime. *GW Notes* **6**, 4–110, 1201.3621 (2013).

120. The LIGO Scientific Collaboration. Advanced LIGO. *Classical Quant. Grav.* **32**, 074001, https://doi.org/10.1088/0264-9381/32/7/074001 (2015).

121. Punturo, M. et al. The Einstein Telescope: a third-generation gravitational wave observatory. *Classical Quant. Grav.* **27**, 194002, https://doi.org/10.1088/0264-9381/27/19/194002 (2010).

122. The Virgo Collaboration. The VIRGO Physics Book, vol. II, https://www.ego-gw.it/public/events/vesf/2010/Presentations/Interferometer_Materials-Vinet.pdf (2006).

123. Berger, J., Davis, P. & Ekström, G. Ambient Earth Noise: A survey of the Global Seismographic Network. *J. Geophys. Res.* **109**, B11307, https://doi.org/10.1029/2004JB003408 (2004).

124. Henry, S., Pozzo di Borgo, E., Danquigny, C. & Abi, B. Simultaneous geomagnetic monitoring with multiple squids and fluxgate sensors across underground laboratories. *E3S Web of Conferences* **12**, 02003, https://doi.org/10.1051/e3sconf/20161202003 (2016).

125. Moore, C. J., Cole, R. H. & Berry, C. P. L. Gravitational-wave sensitivity curves. *Classical Quant. Grav.* **32**, 015014, https://doi.org/10.1088/0264-9381/32/1/015014 (2014).

## Acknowledgements

This work was realized with the financial support of the French State through the "Agence Nationale de la Recherche" (ANR) in the frame of the "Investissement d'avenir" programs: Equipex MIGA (ANR-11-EQPX-0028), IdEx Bordeaux - LAPHIA (ANR-10-IDEX-03-02) and FIRST-TF (ANR-10-LABX-48-01). This work was also supported by the région d'Aquitaine (project IASIG-3D) and by the city of Paris (Emergence project HSENS-MWGRAV). We also acknowledge support from the CPER LSBB2020 project; funded by the "région PACA", the "département du Vaucluse", the MIGA Equipex and the "FEDER PA0000321 programmation 2014-2020". We also thank the "Pôle de compétitivité Route des lasers- Bordeaux" cluster for his support. G.L. thanks DGA for financial support. M.P. thanks LAPHIA-IdEx Bordeaux for partial financial support. M.D.A. and A.F. acknowledge financial support of the Defence Science and Technology Laboratory (DSTL) and the UK National Quantum Technology Hub in Sensors and Metrology with EPSRC Grant No. EP/M013294/1. A.F. was supported by the Science and Technology Facilities Council Consolidated Grant (No. ST/N000633/1).

## Author Contributions

B.C. and A.B. wrote the manuscript. B.C., A.B., R.G., M.P., S.R., W.C., S.P. and P.B. contributed to develop the model of the instrument response and the noise projections. L.A., T.C., J.G., D.H., M.M. contributed to the design of the atom source, vacuum system and mirror suspensions. B.C., L.A., B.F., R.G., J.G., D.H., J.J., G.L., N.M., I.R. and A.L. participated to assembly and preliminary experiments with the atom source. C.D., S.G. and T.M. contributed to geophysics applications of MIGA. M.D.A. and A.F participated to interrogation cavity design. E.C. contributed to interrogation laser design. S.G. contributed to seismic noise analysis and design of infrastructure works. J.H., S.R., Y.R. contributed to the analysis of site gravimetry using the SG instrument. E.P.d.B. and S.H. participated to the analysis of the electromagnetic noise of the site. P.B. supervised the work and the design of the antenna. All authors reviewed the manuscript.

## Additional Information

**Competing Interests:** The authors declare no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_20.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2018