We present an additive-manufactured microwave cavity for a Ramsey-type, double resonance, compact cold-atom clock. Atoms can be laser cooled inside the cavity using a grating magneto-optic trap with the cavity providing an excellent TE011-like mode while maintaining sufficient optical access for atomic detection. The cavity features a low Q-factor of 360 which conveniently reduces the cavity pulling of the future clock. Despite the potential porosity of the additive-manufacturing process, we demonstrate that the cavity is well-suited for vacuum. A preliminary clock setup using cold atoms allows for measuring the Zeeman spectrum and Rabi oscillations in the cavity which enables us to infer excellent field uniformity and homogeneity, respectively, across the volume accessed by the cold atoms. Ramsey spectroscopy is demonstrated, indicating that the cavity is suitable for clock applications. Finally, we discuss the limitations of the future clock.

In the field of atomic clocks and frequency standards, the application of laser cooling and trapping has brought significant breakthroughs, in particular, with cold atomic fountain clocks.1 It is also widely used in more portable or compact industrial clocks2,3 and gravimeters.4,5 In the meantime, Ramsey-type optical-microwave double resonance (DR) clocks based on hot vapor-cells have been studied extensively for future portable,6 compact and space-grade,7,8 and miniature cell9 clocks. Although they achieve state-of-the-art stabilities, they remain limited by a variety of effects such as strong coupling with fluctuating external parameters,6,10 vapor-cell frequency ageing,11 and high relaxation rates. Similar limitations apply to vapor-cell coherent population trapping (CPT) clocks. Cold-atom CPT and DR frequency standards are, hence, excellent candidates for the next generation of compact atomic clocks.

To reduce the cooling and trapping complexity for reduced size, weight, and power (SWaP) applications, grating magneto-optical traps (GMOTs) can be used to create a 3D MOT.12 In this case, a silicon grating is etched to diffract one incoming beam into three opposing beams. The four beams allow the creation of a 3D MOT. This technology has already been applied to CPT Ramsey clocks13–15 and a CPT mini-fountain.16 However, CPT clocks’ low contrast and strong light-shift sensitivities make them less stable in the short- and long-term, respectively. Unlike in miniature cell clocks, portable and high-performance CPT and DR Ramsey clocks based on vapor-cell or cold atoms share similar SWaPs, which eventually favors the latter at the price of developing a microwave cavity allowing for cooling and trapping the atoms.

To date, different microwave cavity geometries for compact cold-atom and vapor-cell clocks have been used: spherical,2 cylindrical,17,18 and loop-gap resonator (LGR).19–22 lGR techniques allow the creation of sub-wavelength, and hence, more compact cavities as small as a few mm wide.23 Additive manufacturing (AM) of such microwave cavities for compact vapor-cell clocks has already been demonstrated24,25 with a typical quality factor of Q c 70. While a high cavity quality factor favors low phase shift and, hence, more accurate clocks,19 a low-quality factor lowers the cavity-pulling effect26 which is of great interest for realizing more stable clocks.

Here, we present the conception, simulation, and realization of a vacuum-grade aluminum AM microwave cavity for a compact 87Rb cold-atom clock. The atom cooling and trapping are achieved using only one laser beam incident on a GMOT hosted inside the cavity.

In DR Ramsey 87Rb clocks, the clock transition with ν 0 = 6.835 GHz [see Fig. 1(c)] is interrogated with light and microwaves in a pulsed manner. After a trapping and cooling stage for a time t l o a d, a short optical state selection stage is utilized to prepare the atoms in the relevant initial clock state for interrogation. This state preparation beam propagates perpendicular to the cavity axis. The atoms are then driven between the atomic levels | F = 2 , m F = 0 and | F = 1 , m F = 0 using two π/2 microwave pulses of length t m separated by the Ramsey time T R. The atomic populations are subsequently readout using an absorption method, with an additional beam, also aligned perpendicular to the cavity axis co-parallel to the state selection beam, during a time t d e t. This concept is illustrated in Figs. 1(a) and 1(b).

FIG. 1.

(a) Optical intensity with respect to time. Pulses lengths are not to scale. The different pulses heights represent the differences in intensities. (b) π / 2 microwave pulses separated by Ramsey time T R. (c) Relevant 87Rb D2 line energy levels with corresponding excitation frequencies with respect to time.

FIG. 1.

(a) Optical intensity with respect to time. Pulses lengths are not to scale. The different pulses heights represent the differences in intensities. (b) π / 2 microwave pulses separated by Ramsey time T R. (c) Relevant 87Rb D2 line energy levels with corresponding excitation frequencies with respect to time.

Close modal
In these conditions, the so-called light-shift effects on the clock transition are minimized compared to CW or CPT operation where the light is on during the clock part of the cycle. In principle, the short-term stability of the clock linearly improves with t m and T R as the atomic transition quality factor Q a equals,27 
(1)
and the short-term stability of a signal-to-noise ratio (SNR)-limited clock is estimated as
(2)
with T C = t l o a d + t p r e p + 2 t m + T R + t d e t is the total cycle time and C is the fringe contrast.

Cold atoms allow for increased interrogation times T R due to low atomic speed (low atom temperature) and, hence, extended interrogation time in a finite interaction region. High-performance vapor-cell Ramsey operations typically show Ramsey times of the order of 3 ms6, whereas a GMOT CPT mini-fountain demonstrated T R as long as 100 ms.16 

Atom trapping and cooling is achieved using a diffraction grating, placed at the bottom of the cavity (see Fig. 2). In this geometry, a single input beam is sufficient. Two side holes in the cavity body are placed at the relevant height to prepare and probe the atom cloud. This simplifies the cavity design requirements as only three holes instead of the traditional six are necessary. Furthermore, the flat grating allows us to interrogate the atoms right after the MOT sequence which would not be possible in a pyramid structure.20 In principle, the atoms can be interrogated either with an absorption or fluorescence measurement.

FIG. 2.

Concept of a microwave cavity with integrated grating. The ensemble is placed in a vacuum chamber.

FIG. 2.

Concept of a microwave cavity with integrated grating. The ensemble is placed in a vacuum chamber.

Close modal
The LGR cavity geometry model is presented in Fig. 3. With this electrode structure, cavities much smaller than the wavelength of the mode can be realized which is of great interest for compact physics packages. The use of additive manufacturing allows for greater control and repeatability on such exotic electrode structures which cannot be easily attained with traditional manufacturing, especially when it comes to producing unique and small batches. A coarse estimate of the cavity resonance frequency ν c can be made using Eq. (3) with n being the number of electrodes, R and r the outer and inner diameters, and t and w the width and thickness of the gaps,28 
(3)
with ϵ ( F / m ) and μ ( N / A 2 ) being the permittivity and permeability inside of the cavity, respectively. Due to the absence of a vapor cell, the design is slightly simpler and the values for vacuum permittivity and permeability are used. Although Eq. (3) yields the correct order of magnitude for the mode frequency, it is necessary to take the full geometry into account such as the height of the cavity h, the electrode support angle α (see Fig. 1), the cavity material and the placement and geometry of the microwave coupling loop.
FIG. 3.

Left: planar view of the LGR geometry. The dashed line indicates the zero-degree reference for the microwave coupling loop, see Fig. 8. Right: AM cavity body prior to holes being drilled.

FIG. 3.

Left: planar view of the LGR geometry. The dashed line indicates the zero-degree reference for the microwave coupling loop, see Fig. 8. Right: AM cavity body prior to holes being drilled.

Close modal

The grating is excluded from the simulations for the sake of time saving and because it has shown little impact on the final mode frequency. The cavity geometry is slightly different from previous vapor-cell designs with the addition of newer constraints. First, the presence of two 4-mm diameter side holes for the purpose of MOT diagnostics, state preparation, and detection at the position of the atom cloud. Second, an inner cavity diameter that is sufficiently large to house the 2 × 2 cm2 grating chip and an upper aperture wide enough to allow the alignment of the trapping and cooling beam onto the grating surface.

These last requirements mainly put lower limit constraints on R and r.

Eigenfrequency analysis sweeping n, h, α, R, r, w, and t is first performed to find a suitable set of parameters yielding a TE011-like mode. The distance to the closest neighboring mode should be much larger than the full width at half maximum (FWHM) of the mode. An example of such a sweep is shown in Fig. 4, sweeping t and w which shows that a large set of ( t , w ) combinations is suitable for realizing the cavity indicated by the area between the two white contours. The figure shows the chosen combination which favors high t and w so that eventual manufacturing inaccuracies have a smaller impact. For all suitable configurations, the closest neighboring mode is at least 800 MHz away.

FIG. 4.

Best TE011-like mode frequency (GHz). The black line represents configurations with 6.8347 GHz mode frequency. The white lines are ±100 MHz boundaries that represent the range that can be recovered by adjusting the cavity height.

FIG. 4.

Best TE011-like mode frequency (GHz). The black line represents configurations with 6.8347 GHz mode frequency. The white lines are ±100 MHz boundaries that represent the range that can be recovered by adjusting the cavity height.

Close modal

The expected magnetic field Hz component and phase along the central axis of the cavity above the grating are presented in Fig. 5. In the current experimental design, the atoms are prepared above the grating and then fall toward it in the least uniform part of the cavity field. Future designs should make use of the most uniform part of the cavity by pushing the atoms away from the grating. In all cases, as the field spatially varies between the two points where the atoms are interrogated, the pulses might have to show different injected microwave power or duration to ensure π / 2 pulses.19 In a sphere of 1 mm diameter—which is of the order of the initial size of the atomic cloud—at its initial position, the Hz field and phase show deviations of 4% and 42 ppm, respectively.

FIG. 5.

Hz amplitude deviation and phase above the grating. The dashed line represents the initial atom cloud position. The deviation is normalized by the value of the Hz field at the initial position of the atoms. The gray area shows the region accessible to the first version of the clock setup.

FIG. 5.

Hz amplitude deviation and phase above the grating. The dashed line represents the initial atom cloud position. The deviation is normalized by the value of the Hz field at the initial position of the atoms. The gray area shows the region accessible to the first version of the clock setup.

Close modal

The final cavity model is presented in Fig. 6, to scale. The cavity body shown in Fig. 3 was manufactured by a commercial additive-manufacturing service, using a standard CL31AL material according to DIN EN 1706 AlSi10Mg(b), while the other simpler and less critical parts were manufactured in our workshop in 6082 Al alloy.

FIG. 6.

Final cavity model, to scale.

FIG. 6.

Final cavity model, to scale.

Close modal

Figure 7 compares the simulated and measured cavity S11 microwave reflection spectrum over a wide frequency range. The simulation and measurement agree well on the mode positions with a small discrepancy of 60 MHz for the mode of interest, which can be easily compensated by adjusting the height of the cavity. As predicted by the simulation, the left and right neighboring modes are 1.2 and 1.8 GHz away, respectively. They do not overlap with the mode of interest so that the field geometry is not degraded. The 8.5 GHz resonance is associated with a mode whose field energy is concentrated around the microwave coupling loop and low sensitivity to the other cavity parameters. The loop is realized as a small, rectangular copper wire loop of approximately 9 mm2 cross section with one end connected to the SMA core.29 According to our simulations, the precise loop geometry does not strongly influence the mode of interest, but rather the 8.5 GHz mode, which is of no major concern for our application due to the large distance between the two modes. However, the loop position relative to the electrodes is quite critical for the efficient coupling of microwave radiation to the cavity, as illustrated by Fig. 8.

FIG. 7.

Simulated and measured S11 spectrum of the cavity. The main figure and the inset share the same axis units.

FIG. 7.

Simulated and measured S11 spectrum of the cavity. The main figure and the inset share the same axis units.

Close modal
FIG. 8.

Measured microwave mode peak depth with respect to the rotation angle between the microwave coupling loop and the cavity. See Fig. 3 for the 0° reference position. Given the angle, coupling to the mode of interest can be either enhanced or depleted.

FIG. 8.

Measured microwave mode peak depth with respect to the rotation angle between the microwave coupling loop and the cavity. See Fig. 3 for the 0° reference position. Given the angle, coupling to the mode of interest can be either enhanced or depleted.

Close modal

The simulated and measured S11 of the 6.835 GHz mode (see the inset of Fig. 7 and Table I) differ in depth and FWHM. This is explained first because the simulation did not include any wall roughness which is not realistic, given the AM process. Second, for the sake of simulation time, the grating was not included in the simulation. With its etched nature, it exhibits increased impedance to surface currents which also participates in power losses and, hence, line broadening. Lastly, the loss of mode amplitude is attributed to the high sensitivity of the mode to the angle between the microwave coupling loop and the electrodes, see the measured data in Fig. 8. Although the angle was manually optimized to get a peak depth around −17 dB, corresponding to the loop placed approximately under one of the electrode wings, it was not perfectly preserved when securing the cavity, as evidenced by Fig. 7.

TABLE I.

Cavity mode FWHM and Q-factor simulation vs measurement comparison.

FWHM (MHz)QcDepth (dB)
Simulation 1709 −23 
Measurement 19 360 −9 
FWHM (MHz)QcDepth (dB)
Simulation 1709 −23 
Measurement 19 360 −9 
The rather low-quality factor minimizes the cavity-pulling effect, indeed, the cavity-pulling shift,6,26
(4)
with Q c and Q a being the cavity and atomic transition quality factors, respectively, c ( θ ) a function of the microwave pulse area that can be made < 0.005,26 and Δ ν c the cavity detuning which, under vacuum, only depends on the cavity temperature. The temperature coefficient is measured as Δ ν c / Δ T c a v = ( 156 ± 9 ) kHz / K which is similar to the other models of AM cavities.6,25

The cavity is mounted in a commercial vacuum chamber, with the grating placed at the bottom of the cavity. Resistively heated dispensers provide the necessary 87Rb vapor. The outgassing properties of the aluminum-based additively manufactured material were studied in a dedicated test chamber using the gas accumulation method.30 Its outgassing properties were compared to a sample of bulk aluminum. Both samples show similar outgassing rates estimated around 7 × 10−9 mbar m3 s−1 m−2 which suggests little to no degradation caused by the additive-manufacturing material compared to the traditional aluminum alloy. Furthermore, measurement of the cold-atom cloud 1/e lifetime31 gives an estimate of the pressure of the order of 1 × 10 8 mbar, which also indicates this material’s suitability for in-vacuum cold-atom applications.

The magnetic field is controlled with three orthogonal sets of Helmholtz coils placed outside of the chamber to minimize the field in the grating plane and apply a C-field along the principal axis of the cavity to separate the Zeeman sublevels. A pair of anti-Helmholtz coils are placed inside the vacuum chamber to allow trapping and cooling. The cavity loop is connected to a commercial microwave generator (Keysight E8257D).

A single 780 nm extended cavity diode laser (ECDL), frequency stabilized with saturated absorption spectroscopy is used throughout. This light is passed through an electro-optical modulator (EOM) to allow the generation of optical sidebands on carrier frequency for hyperfine repumping of the atoms. A tapered amplifier is then used to amplify the light before dividing it into three beam paths for trapping, state selection, and state-detection. Each beam path contains a double-pass acousto-optic modulator (AOM) for intensity and frequency control. Trapping light is coupled into an optical fiber, the output of which is incident on the grating through an aperture in the top of the cavity after expansion and collimation. State selection and detection light are coupled into a common additional fiber. From the fiber, the light is linearly polarized with respect to the magnetic bias field and split with a non-polarizing beam splitter (NPBS). This allows a portion of the light to be monitored on a reference photodiode while the rest of the light is aligned through the side holes in the cavity body and retro-reflected by a mirror back onto a second signal photodiode. This allows normalization of the absorption signal for fluctuations in laser intensity during state-detection.

The experimental cycle is presented in Fig. 1. The traditional Ramsey cycle is preceded by a trapping and cooling phase by optical molasses of duration t l o a d where the intensity is linearly ramped down from 5 I s a t to 0.1 I s a t,32 while simultaneously scanning the detuning from 12to 30 MHz. After this step, 1 × 10 6 atoms are in the states | F = 2 , m F = ± 2 , ± 1 , 0 with temperature 10 μ K. A state selection phase of duration t p r e p utilizing light of intensity 0.02 I s a t that is linearly polarized with respect to the C-field is then used to prepare the atoms in the clock state. This light is incident on the atoms through the cavity side-hole and is tuned to the | F = 2 | F = 2 state with repump light driving the | F = 1 | F = 2 transition. Due to selection rules, this pumps more than 95% of the atoms into the | F = 2 , m F = 0 level,33 increasing the signal contrast.

The atoms are then driven with two π / 2 microwave pulses of duration t m = 167 μ s separated by a free evolution (or Ramsey) time T R. The population P 2 in level | F = 2 is estimated by measuring the absorption of a side beam tuned on | F = 2 | F = 3 . This value is normalized by the absorption of the same beam with 5 % | F = 1 | F = 2 sideband. This second pulse recycles atoms initially shelved in the F = 1 ground-state back to the cycling transition and allows normalization for variation of the number of atoms present at each shot. The C-field is optimized for atom cooling and is kept constant during the clock cycle to avoid eddy currents.

Rabi oscillations (see Fig. 9) are measured by setting T R = 0 ms. To sweep the microwave pulse area, the microwave time is kept constant and the microwave power is swept. The opposite approach, constant power and sweeping microwave time is disregarded to avoid the influence of time-dependent relaxation processes and atom losses. The first minimum of the Rabi oscillations is used to calibrate the input microwave power to generate π / 2 pulses.

FIG. 9.

Rabi oscillations. The pulse area axis is rescaled assuming the first minimum is π.

FIG. 9.

Rabi oscillations. The pulse area axis is rescaled assuming the first minimum is π.

Close modal
The decay of the Rabi oscillation is due to the residual magnetic field inhomogeneities inside of the atomic volume.34 The number of visible oscillations, 9, is a sign of excellent H z homogeneity, which is comparable to other compact cold-atom clocks that show 14 or more oscillations,35 whereas high-performance vapor-cell clocks such as a 87Rb Ramsey maser show several oscillations as low as 3.36 In the case of vapor-cell clocks, the atoms cover a larger volume inside the cavity, and the field homogeneity is dominated by the TE011-like distribution. In contrast, the oscillation decay presented in Fig. 9 is compatible with a pulse area θ with a gaussian distribution. Indeed, if one assumes a gaussian distribution with a standard deviation σ, relative to the pulse area mean θ, one can show that the form of the textbook co s 2 ( θ / 2 ) Rabi oscillations is modified as follows:
(5)
with A and C as two fit parameters that account for the absorption measurement. With this model, from a fit to Fig. 9, we measure σ = ( 6.3 ± 0.1 ) %. The distribution of the field felt by the atoms is given by the convolution of the distribution of the atoms—which is assumed gaussian—and the distribution of the TE011-like field. The smaller the volume spanned by the atoms, the more uniform, the field distribution.

The Zeeman spectrum is presented in Fig. 10. It is intentionally acquired with repump light in the state selection beam extinguished, shelving atoms in the | F = 1 ground-state. The Ramsey cycle is set with T R = 0 ms and optimal microwave power deduced from the Rabi oscillations. Note that the amplitude of π transitions (i.e., transitions with m F = m F ) peaks at ± 1.1 MHz are slightly asymmetric which is attributed to a small imbalance in the σ ±-polarized beams used for trapping and cooling. The Zeeman spectrum shows peaks instead of dips because, in the absence of repump light, the majority of the atoms populate the sublevels | F = 1 , m F = ± 1 , 0 before the microwave cycle.

FIG. 10.

Zeeman spectrum without state selection. The peaks circled in red are the σ ± transitions.

FIG. 10.

Zeeman spectrum without state selection. The peaks circled in red are the σ ± transitions.

Close modal
Note that the σ ± transitions (i.e., transitions with m F = m F ± 1) peaks are only slightly above the noise level, which indicates excellent field homogeneity and alignment with the C-field. The Zeeman spectrum allows us to measure a field orientation factor ξ of 97 %, which is a measure of field uniformity that was defined in Ref. 29, namely,
(6)
with S π and S σ being the Zeeman spectrum π and σ components only, respectively. V is the volume spanned by the interrogated atoms.

Ramsey fringes were obtained for Ramsey times up to 20 ms, limited in the current setup by the atoms falling under gravity and leaving the transverse optical beam after this time.37  Figure 11 shows the optimal37 SNR Ramsey fringes with T R = 10 ms. The FWHM of the central fringe is 49 Hz and is compatible with Eq. (1). The contrast is 58% due to the state selection. Without the state selection, the maximum contrast is 33%. The noise on the fringes is primarily attributed to magnetic field noise.

FIG. 11.

Ramsey fringes for T R = 10 ms and t m = 167 μ s. The main figure and the inset share the same axes.

FIG. 11.

Ramsey fringes for T R = 10 ms and t m = 167 μ s. The main figure and the inset share the same axes.

Close modal

The nominal cycle time sequence is shown in Table II.

TABLE II.

Nominal time sequence.

PulseDuration
Tload 516 ms 
tprep 1 ms 
tm 167 μ
TR 10 ms 
tdet 1.4 ms 
TC (total) 528.7 ms 
PulseDuration
Tload 516 ms 
tprep 1 ms 
tm 167 μ
TR 10 ms 
tdet 1.4 ms 
TC (total) 528.7 ms 

At one second, the Allan deviation of the future clock should be σ ( τ = 1 s ) = 4.4 × 10 11. The SNR of Eq. (2) is estimated from a sine fit to the central fringe of Fig. 11.38 A preliminary estimation shows that the Dick effect contribution39 should be at the level of σ D i c k ( τ = 1 s ) = 1.2 × 10 12. A more detailed short-term analysis of the increased performances clock is presented in a more advanced version of the setup.37 Finally, we estimate the cavity pulling [see Eq. (4)] to be of the order of 2.1 × 10 15 for a reasonable relative cavity temperature fluctuation of 10 3, which is typical at 10 5 s.

We demonstrated the successful design, simulation, and implementation of a microwave cavity for a novel compact cold-atom DR clock based on intra-cavity atom cooling using a diffractive optical element. An efficient simulation workflow allows us to design the cavity with a small mismatch of the resonance frequency which is easily compensated via a small adaptation of the cavity height. Rabi oscillations and Zeeman spectrum show proof of good field uniformity and homogeneity, respectively.

Short-term performances of the future proof-of-concept clock are estimated to be closer than that of miniature vapor-cell clocks,9,40 with a Dick effect limitation, one order of magnitude lower. Accordant stability measurement has been demonstrated in an early-stage clock setup with an optimized signal.37 In order to improve on the short-term stability of future setups, the Ramsey time could be increased by either putting the cavity upside down and interrogating the atoms further down or operating the clock in a mini-fountain fashion with Ramsey times up to 50 and 100 ms,16 respectively. This will also allow us to interrogate the atoms in regions where the field is more uniform, as suggested by Fig. 5. To improve the Dick effect limit, T C must be reduced, which can be achieved by implementing atom recapturing at the end of the cycle.41,42

E. Batori, C. Affolderbach, and G. Mileti thank V. Dolgovskiy for his contributions in the early design phase, P. Scherler for technical assistance, and G. Bourban for his continued support. This work was financially supported by the Swiss Space Center (Swiss Confederation) and the European Space Agency (ESA, Contract No. 4000131046). The views expressed herein can in no way be taken to reflect the official opinion of the European Space Agency. A.B. was supported by a Ph.D. studentship from the Defence Science and Technology Laboratory (Dstl). We gratefully acknowledge support from the International Network for Microfabrication of Atomic Quantum Sensors (No. EPSRC EP/W026929/1).

The authors have no conflicts to disclose.

Etienne Batori: Conceptualization (equal); Data curation (equal); Formal analysis (equal); Investigation (lead); Methodology (lead); Software (lead); Validation (lead); Visualization (lead); Writing – original draft (lead); Writing – review & editing (lead). Alan Bregazzi: Conceptualization (equal); Data curation (equal); Formal analysis (equal); Investigation (equal); Methodology (equal); Software (equal); Validation (equal); Visualization (supporting); Writing – original draft (supporting); Writing – review & editing (equal). Ben Lewis: Conceptualization (equal); Investigation (equal); Methodology (equal); Software (equal). Paul Griffin: Conceptualization (equal); Data curation (equal); Funding acquisition (equal); Investigation (supporting); Methodology (supporting); Project administration (equal); Resources (equal); Supervision (equal); Validation (supporting); Writing – original draft (supporting); Writing – review & editing (equal). Erling Riis: Conceptualization (equal); Funding acquisition (equal); Project administration (equal); Resources (equal); Supervision (equal); Writing – review & editing (supporting). Gaetano Mileti: Conceptualization (equal); Funding acquisition (equal); Project administration (equal); Resources (lead); Supervision (lead); Writing – review & editing (supporting). Christoph Affolderbach: Conceptualization (equal); Data curation (equal); Formal analysis (equal); Funding acquisition (equal); Investigation (supporting); Methodology (supporting); Project administration (lead); Resources (equal); Supervision (equal); Validation (supporting); Visualization (supporting); Writing – original draft (equal); Writing – review & editing (equal).

The data that support the findings of this study are openly available in figshare at doi.org/10.6084/m9.figshare.23276357, Ref. 23276357.

1.
R.
Wynands
and
S.
Weyers
,
Metrologia
42
,
S64
(
2005
).
2.
F. X.
Esnault
,
N.
Rossetto
,
D.
Holleville
,
J.
Delporte
, and
N.
Dimarcq
,
Adv. Space Res.
47
,
854
(
2011
).
3.
SpectraDynamics
, “CRb–clock” (2019), available at https://spectradynamics.com/products/crb–clock/.
4.
IX Blue
, “Absolute quantum gravimeter” (2021), available at https://www.ixblue.com/quantum-gravimeter/.
5.
Q.
Bodart
,
S.
Merlet
,
N.
Malossi
,
F. P.
Dos Santos
,
P.
Bouyer
, and
A.
Landragin
,
Appl. Phys. Lett.
96
,
134101
(
2010
).
6.
N.
Almat
,
M.
Gharavipour
,
W.
Moreno
,
C.
Affolderbach
, and
G.
Mileti
, in
2019 Joint Conference of the IEEE International Frequency Control Symposium & European Frequency and Time Forum
(
IEEE
,
2019
), pp.
1
2
.
7.
E.
Batori
,
N.
Almat
,
C.
Affolderbach
, and
G.
Mileti
,
Adv. Space Res.
68
(12), 4723–4733 (
2020
).
8.
M.
Huang
,
A.
Little
, and
J.
Camparo
, in
2022 Joint Conference of the European Frequency and Time Forum & IEEE International Frequency Control Symposium
(
IEEE
,
2022
), pp.
1
4
.
9.
E.
Batori
,
C.
Affolderbach
,
M.
Pellaton
,
F.
Gruet
,
M.
Violetti
,
Y.
Su
,
A. K.
Skrivervik
, and
G.
Mileti
,
Phys. Rev. Appl.
18
,
054039
(
2022
).
10.
W.
Moreno
,
M.
Pellaton
,
C.
Affolderbach
, and
G.
Mileti
,
IEEE Trans. Ultrason. Ferroelectr. Freq. Control
65
,
1500
(
2018
).
11.
S.
Abdullah
,
C.
Affolderbach
,
F.
Gruet
, and
G.
Mileti
,
Appl. Phys. Lett.
106
,
163505
(
2015
).
12.
C. C.
Nshii
,
M.
Vangeleyn
,
J. P.
Cotter
,
P. F.
Griffin
,
E. A.
Hinds
,
C. N.
Ironside
,
P.
See
,
A. G.
Sinclair
,
E.
Riis
, and
A. S.
Arnold
,
Nat. Nanotechnol.
8
,
321
(
2013
).
13.
R.
Elvin
,
M. W.
Wright
,
B.
Lewis
,
B. L.
Keliehor
,
A.
Bregazzi
,
J. P.
McGilligan
,
A. S.
Arnold
,
P. F.
Griffin
, and
E.
Riis
,
Adv. Opt. Technol.
9
,
297
(
2020
).
14.
R.
Elvin
,
G. W.
Hoth
,
M.
Wright
,
B.
Lewis
,
J. P.
McGilligan
,
A. S.
Arnold
,
P. F.
Griffin
, and
E.
Riis
,
Opt. Express
27
,
38359
(
2019
).
15.
G. W.
Hoth
,
R.
Elvin
,
M.
Wright
,
B.
Lewis
,
A. S.
Arnold
,
P. F.
Griffin
, and
E.
Riis
,
Proc. SPIE
10934
,
109342E
(
2019
).
16.
B.
Lewis
,
R.
Elvin
,
A. S.
Arnold
,
E.
Riis
, and
P. F.
Griffin
,
Appl. Phys. Lett.
121
,
164001
(
2022
).
17.
F. G.
Ascarrunz
,
Y. O.
Dudin
,
M. C.
Delgado Aramburo
,
L. I.
Ascarrunz
,
J.
Savory
,
A.
Banducci
, and
S. R.
Jefferts
, in
2018 IEEE International Frequency Control Symposium
(
IEEE
,
2018
), pp.
1
3
.
18.
M.
Gozzelino
,
S.
Micalizio
,
C. E.
Calosso
,
A.
Godone
,
H.
Lin
, and
F.
Levi
,
IEEE Trans. Ultrason. Ferroelectr. Freq. Control
68
,
872
(
2021
).
19.
S.
Lee
,
G. W.
Choi
,
H.-G.
Hong
,
T. Y.
Kwon
,
S.-B.
Lee
,
M.-S.
Heo
, and
S. E.
Park
,
Appl. Phys. Lett.
119
,
064002
(
2021
).
20.
M.
Givon
,
L.
Habib
,
A.
Waxman
,
Y.
Bar-Haim
,
O.
Amit
,
Y.
Cina
,
G.
Boneh
,
D.
Groswasser
,
T.
David
,
B.
Levy
,
A.
Stern
, and
R.
Folman
, arXiv:2208.09038 (2022).
21.
Q.
Hao
,
W.
Xue
,
W.
Li
,
F.
Xu
,
X.
Wang
,
W.
Guo
,
P.
Yun
, and
S.
Zhang
,
IEEE Trans. Ultrason. Ferroelectr. Freq. Control
67
,
873
(
2020
).
22.
S.
Kang
,
M.
Gharavipour
,
F.
Gruet
,
C.
Affolderbach
, and
G.
Mileti
, in
2015 Joint Conference of the IEEE International Frequency Control Symposium & European Frequency and Time Forum
(
IEEE
,
2015
), pp.
800
803
.
23.
M.
Violetti
,
M.
Pellaton
,
C.
Affolderbach
,
F.
Merli
,
J.-F.
Zurcher
,
G.
Mileti
, and
A. K.
Skrivervik
,
IEEE Sens. J.
14
,
3193
(
2014
).
24.
M.
Pellaton
,
C.
Affolderbach
,
A. K.
Skrivervik
,
A. E.
Ivanov
,
T.
Debogovic
,
E.
Rijk
, and
G.
Mileti
,
Electron. Lett.
54
,
691
(
2018
).
25.
C.
Affolderbach
,
W.
Moreno
,
A. E.
Ivanov
,
T.
Debogovic
,
M.
Pellaton
,
A. K.
Skrivervik
,
E.
de Rijk
, and
G.
Mileti
,
Appl. Phys. Lett.
112
,
113502
(
2018
).
26.
S.
Micalizio
,
A.
Godone
,
F.
Levi
, and
C.
Calosso
,
IEEE Trans. Ultrason. Ferroelectr. Freq. Control
57
,
1524
(
2010
).
27.
S.
Micalizio
,
C. E.
Calosso
,
F.
Levi
, and
A.
Godone
,
Phys. Rev. A
88
,
33401
(
2013
).
28.
W.
Froncisz
and
J. S.
Hyde
,
J. Magn. Reson.
47
,
515
(
1982
).
29.
C.
Stefanucci
,
T.
Bandi
,
F.
Merli
,
M.
Pellaton
,
C.
Affolderbach
,
G.
Mileti
, and
A. K.
Skrivervik
,
Rev. Sci. Instrum.
83
,
104706
(
2012
).
30.
P. A.
Redhead
,
J. Vac. Sci. Technol. A
20
,
1667
(
2002
).
31.
C.
Monroe
,
W.
Swann
,
H.
Robinson
, and
C.
Wieman
,
Phys. Rev. Lett.
65
,
1571
(
1990
).
32.
D. A.
Steck
, “Rubidium 87 D Line Data,” available online at http://steck.us/alkalidata (revision 2.2.2, 9 July 2021).
33.
Y.-X.
Duan
,
B.
Wang
,
J.-F.
Xiang
,
Q.
Liu
,
Q.-Z.
Qu
,
D.-S.
, and
L.
Liu
,
Chin. Phys. Lett.
34
,
073201
(
2017
).
34.
W.
Moreno
,
M.
Pellaton
,
N.
Almat
,
M.
Gharavipour
,
C.
Affolderbach
, and
G.
Mileti
, in
2019 Joint Conference of the IEEE International Frequency Control Symposium & European Frequency and Time Forum
(
IEEE
,
2019
).
35.
P.
Liu
,
Y.
Meng
,
J.
Wan
,
X.
Wang
,
Y.
Wang
,
L.
Xiao
,
H.
Cheng
, and
L.
Liu
,
Phys. Rev. A
92
,
62101
(
2015
).
36.
A.
Godone
,
S.
Micalizio
,
F.
Levi
, and
C.
Calosso
,
Rev. Sci. Instrum.
82
,
074703
(
2011
).
37.
A.
Bregazzi
,
E.
Batori
,
B.
Lewis
,
C.
Affolderbach
,
G.
Mileti
,
E.
Riis
, and
P. F.
Griffin
, arXiv:2305.02944.
38.
S.
Negusse
and
P.
Zetterberg
,
IEEE Trans. Instrum. Measurement
63
,
1514
(
2014
).
39.
S.
Micalizio
,
A.
Godone
,
C.
Calosso
,
F.
Levi
,
C.
Affolderbach
, and
F.
Gruet
,
IEEE Trans. Ultrason. Ferroelectr. Freq. Control
59
,
457
(
2012
).
40.
C.
Carlé
,
M.
Petersen
,
N.
Passilly
,
M. A.
Hafiz
,
E.
De Clercq
, and
R.
Boudot
,
IEEE Trans. Ultrason. Ferroelectr. Freq. Control
68
,
3249
(
2021
).
41.
F.-X.
Esnault
,
E.
Blanshan
,
E. N.
Ivanov
,
R. E.
Scholten
,
J.
Kitching
, and
E. A.
Donley
,
Phys. Rev. A
88
,
42120
(
2013
).
42.
A. V.
Rakholia
,
H. J.
McGuinness
, and
G. W.
Biedermann
,
Phys. Rev. Appl.
2
,
054012
(
2014
).