Skip to content

Precision AMO experiments translate small changes in a field, frequency, or acceleration into quantum phases and measured populations. Their performance depends on more than the smallest reported noise: interrogation time, contrast, bandwidth, repetition rate, calibration, and the interpretation of the signal all matter. This page explains these established relationships and how to compare sensors under the same conditions.

For dated measurements, unresolved questions, and possible research projects, see Precision AMO tests. The calculations below are idealized examples, not forecasts for a particular instrument.

Required background. Precision metrology introduces uncertainty and measurement models. Ramsey interferometry explains phase-to-population conversion, and Fisher information provides the local estimation bound used below.

For small perturbations about a calibrated operating point, a linear sensor can be described by a signal response, unwanted responses, and noise. If the response is also time invariant, channel ii has the frequency-domain form

y~i(f)=Hiθ(f)θ~(f)+∑jHij(f)x~j(f)+n~i(f).\widetilde y_i(f)=H_{i\theta}(f)\widetilde\theta(f) +\sum_jH_{ij}(f)\widetilde x_j(f)+\widetilde n_i(f).

Here θ\theta is the quantity to be estimated, xjx_j are environmental or control inputs, and nin_i is residual noise. The transfer function HH tells you both how strongly and with what phase the sensor responds. For pulsed, time-varying, or nonlinear measurements, use the corresponding two-time response kernel or full measurement model instead of assuming a convolution. Degen, Reinhard, and Cappellaro (2017) develop this way of comparing quantum sensors.

Device Simple response and its conditions
Atomic magnetometer ωL=γB\omega_L=\gamma B for the chosen Zeeman manifold and field regime; the response also depends on pumping, relaxation, and measurement bandwidth.
Clock Detuning accumulates a phase during interrogation; the final population depends on that phase, contrast, and readout.
Light-pulse accelerometer ϕa=keff⋅a T2\phi_a=\mathbf k_{\mathrm{eff}}\cdot\mathbf a\,T^2 for the ideal three-pulse sequence with constant acceleration and negligible pulse duration.
Frequency-ratio comparison Different transitions respond differently to dimensionless constants and environmental shifts.

These scale factors are not complete instruments. Spatial averaging, field gradients, propagation delays, state preparation, and detection can modify the measured response.

Suppose a weak parameter θ(t)\theta(t) shifts an angular frequency by gθ(t)g\theta(t), with constant coupling gg. Neglect pulse duration and decay for the moment. A Ramsey free-evolution interval from 00 to TT accumulates

ϕ(T)=g∫0Tθ(t) dt.\phi(T)=g\int_0^T\theta(t)\,dt.

For a complex harmonic input θ(t)=θ0ei2πft\theta(t)=\theta_0e^{i2\pi ft}, whose real part is the physical signal, the phase amplitude is HT(f)θ0H_T(f)\theta_0, with

HT(f)=gTeiπfTsinc⁡(πfT),sinc⁡(x)=sin⁡xx.H_T(f)=gT e^{i\pi fT}\operatorname{sinc}(\pi fT), \qquad \operatorname{sinc}(x)=\frac{\sin x}{x}.

At zero frequency ∣HT∣=∣g∣T|H_T|=|g|T: longer interrogation amplifies a static signal. But the first zero moves to f=1/Tf=1/T. A signal that completes one full cycle during the interval contributes zero net phase, regardless of its instantaneous amplitude.

Magnitude of ideal Ramsey response for interrogation times T0 and twice T0, showing a larger static response and a lower first zero for the longer interval.

The exact ideal-filter magnitudes ∣HT(f)∣/(∣g∣T0)|H_T(f)|/(|g|T_0) for T=T0T=T_0 and T=2T0T=2T_0. Doubling the interval doubles the static response and halves the frequency of the first zero. The curves describe phase accumulation before decoherence, sampling between shots, or detection noise is included.

An echo or a sequence of control pulses changes the sign or weighting of the integrand and can select an oscillatory signal instead. The corresponding sensitivity function is the natural starting point for comparing protocols at the frequency of interest.

For a regular one-parameter statistical model and a locally unbiased estimator, ν\nu independent repetitions obey

Var⁡(θ^)≥1νFC≥1νFQ.\operatorname{Var}(\hat\theta)\geq\frac{1}{\nu F_C} \geq\frac{1}{\nu F_Q}.

FCF_C is the Fisher information in the chosen measurement, while FQF_Q is the maximum permitted by the prepared state. A state with large FQF_Q is helpful only when the readout accesses that information. Pezzè et al. (2018) explain this connection for entangled atomic ensembles.

If each cycle lasts TcT_c, independent repeated measurements have asymptotic information rate FC/TcF_C/T_c. Preparation, cooling, readout, and calibration time count alongside interrogation. Unknown systematic bias is not removed by this statistical bound; correlated cycles require a joint noise model.

A longer interrogation can reduce the information rate

Section titled “A longer interrogation can reduce the information rate”

For NN independent ideal two-state probes measured at the steepest point of a Ramsey fringe, let the contrast be C(T)=e−T/T2C(T)=e^{-T/T_2}. The binary-outcome Fisher information for a static θ\theta is

FC(T)=Ng2T2e−2T/T2.F_C(T)=N g^2T^2e^{-2T/T_2}.

With dead time TdT_d, the rate is

F˙C(T)=Ng2T2e−2T/T2T+Td.\dot F_C(T)=\frac{N g^2T^2e^{-2T/T_2}}{T+T_d}.

If Td=0T_d=0, differentiating its logarithm gives 1/T−2/T2=01/T-2/T_2=0, so the optimum is T=T2/2T=T_2/2. For nonzero TdT_d, the stationary condition is 2/T−2/T2−1/(T+Td)=02/T-2/T_2-1/(T+T_d)=0. This is an example under fixed NN and exponential contrast decay, not a universal optimum. Technical noise and a different signal waveform change it.

Accepted trials are only part of the experiment

Section titled “Accepted trials are only part of the experiment”

If a fraction pp of attempts is accepted independently of θ\theta, and rejected trials contain no useful information, the rate from accepted trials is pFC,acc/TcpF_{C,\mathrm{acc}}/T_c. Quoting only the information per accepted trial can hide a low success probability.

When acceptance itself depends on θ\theta, the accept/reject record carries information. For a resolved binary selection flag, the full classical information per attempted trial decomposes as

FC,full=[∂θp]2p(1−p)+pFC,acc+(1−p)FC,rej,0<p<1.F_{C,\mathrm{full}}= \frac{[\partial_\theta p]^2}{p(1-p)} +pF_{C,\mathrm{acc}}+(1-p)F_{C,\mathrm{rej}}, \qquad 0<p<1.

This follows by factoring the joint probability into the flag probability and the conditional outcome distribution. Conditional scores have zero mean, so their cross terms vanish. A valid comparison uses the same total resources and specifies which records are retained.

Separate fluctuations of repeated data from uncertainty in corrections and in the physical response model. When contributions are uncorrelated, their standard uncertainties can be added in quadrature. Shared references, constants, or environmental corrections instead produce covariance.

For example, two measurements with equal variance u2u^2 and correlation ρ\rho have

umean2=u22(1+ρ),udifference2=2u2(1−ρ).u_{\mathrm{mean}}^2=\frac{u^2}{2}(1+\rho), \qquad u_{\mathrm{difference}}^2=2u^2(1-\rho).

Positive correlation weakens the improvement from averaging but can help a difference reject common noise. Whether that also rejects the desired signal depends on its response in the two channels.

For three positive frequency ratios with directions RAB=νA/νBR_{AB}=\nu_A/\nu_B, RBC=νB/νCR_{BC}=\nu_B/\nu_C, and RCA=νC/νAR_{CA}=\nu_C/\nu_A, consistency requires

ln⁡RAB+ln⁡RBC+ln⁡RCA=0.\ln R_{AB}+\ln R_{BC}+\ln R_{CA}=0.

This is a useful experimental test only when at least one ratio is measured independently. Defining the third ratio from the other two makes closure an algebraic identity. Include shared-oscillator and transfer-link covariances when assigning an uncertainty to the closure residual. Measurements taken at different times also require a model for actual frequency changes.

Optical clocks and Ludlow et al. (2015) explain the clock-specific corrections. Near Earth’s surface, the weak-field relation Δν/ν≃gΔh/c2\Delta\nu/\nu\simeq g\Delta h/c^2 gives a fractional shift near 1.1×10−191.1\times10^{-19} per millimetre. That conversion alone does not establish millimetre height resolution: averaging time, transfer uncertainty, and the geopotential model still matter.

The atom-interferometer scale factor grows as T2T^2, but the phase is periodic. For a one-axis accelerometer the acceleration change producing 2π2\pi phase is 2π/(keffT2)2\pi/(k_{\mathrm{eff}}T^2). A population measurement of a cosine fringe has additional within-period ambiguities unless a second quadrature, a prior, or a lock identifies the operating branch.

Increasing TT can therefore improve local sensitivity while making it harder to follow a rapidly changing acceleration. An auxiliary accelerometer, different interrogation times, or adaptive phase control can help identify the fringe. Cronin, Schmiedmayer, and Pritchard (2009) develop the interferometer response; atom-interferometric sensors connects it to practical noise and corrections.

Many precision tests use a linearized response matrix

y=Ac+ϵ,\mathbf y=A\mathbf c+\boldsymbol\epsilon,

where c\mathbf c contains unknown interaction coefficients and AA combines atomic, molecular, nuclear, and instrument responses. The rank of AA determines which combinations are identifiable. Almost parallel response columns can leave large uncertainty in individual coefficients even when the measured combinations are precise. Setting all but one coefficient to zero is an additional model assumption.

A spin-independent Yukawa modification is one illustrative parameterization:

V(r)=−Gm1m2r(1+αYe−r/λ),λ=ℏmbc.V(r)=-\frac{Gm_1m_2}{r}\left(1+\alpha_Ye^{-r/\lambda}\right), \qquad \lambda=\frac{\hbar}{m_bc}.

αY\alpha_Y describes a strength relative to gravity for the specified sources. Composition-dependent or spin-dependent forces require their own charges and operators. The range alone does not determine a detector’s response: source geometry, gradients, and shielding or screening assumptions can matter.

Likewise, a bound on an oscillation of a dimensionless constant is not yet a bound on every dark-matter model. The conversion to field coupling may require the field density and coherence. Those assumptions are distinct from the measured frequency-ratio limit itself. Safronova et al. (2018) develop these interpretations; fundamental symmetry tests and variation-of-constants searches give the relevant AMO response models.

For a useful comparison, hold the physical task fixed. A static magnetic field, a known-frequency oscillation, and an unknown transient require different estimators. State the observation time, accepted frequency band, total cycle time, particle or photon resources, and environmental conditions.

Then compare the quantity that answers the question: an uncertainty after a specified time, a frequency-dependent detection probability, or a model-bound curve with a stated coverage procedure. A low amplitude spectral density in one band does not establish a low uncertainty for all waveforms. Reversing a field or comparing two sensors can isolate a response symmetry, but an unwanted effect can have that same symmetry.

For research directions and measurements with their publication dates, use Precision AMO tests. The experimental implementations are developing; the transfer functions, statistical assumptions, and limiting calculations above explain what those developments must be compared against.

Protocol A uses uncorrelated atoms, has Fisher information FA=1.0×104F_A=1.0\times10^4 per shot, and cycles every 1.0 s1.0\ \mathrm{s}. Protocol B uses a squeezed state, has FB=2.5×104F_B=2.5\times10^4 per accepted shot, takes 2.0 s2.0\ \mathrm{s} per attempted cycle, and accepts 80%80\% of attempts. Assume the acceptance probability is parameter independent, rejected shots contain no information, and the measurements attain their classical Cramér–Rao bounds in the repeated-trial limit.

  1. Compute the Fisher-information rate for each protocol.
  2. Which has the smaller ideal statistical variance after one hour?
  3. By what factor?
Solution

For A,

F˙A=1041.0 s=1.0×104 s−1.\dot F_A = \frac{10^4}{1.0\ \mathrm{s}} = 1.0\times10^4\ \mathrm{s}^{-1}.

For B, the accepted information per attempt is 0.80FB=2.0×1040.80F_B=2.0\times10^4, so

F˙B=2.0×1042.0 s=1.0×104 s−1.\dot F_B = \frac{2.0\times10^4}{2.0\ \mathrm{s}} = 1.0\times10^4\ \mathrm{s}^{-1}.

The protocols have equal ideal information rate and therefore equal Cramér–Rao variance after one hour. Protocol B has a single-shot resource gain but no wall-clock statistical gain under the stated overhead. Their systematic errors and dynamic ranges would still need comparison.

2. Atom-interferometer scale and phase ambiguity

Section titled “2. Atom-interferometer scale and phase ambiguity”

An ideal light-pulse accelerometer has keff=1.61×107 m−1k_{\mathrm{eff}}=1.61\times10^7\ \mathrm{m}^{-1} and T=0.20 sT=0.20\ \mathrm{s}.

  1. Find the phase change caused by δa=1.0×10−8 m s−2\delta a=1.0\times10^{-8}\ \mathrm{m\,s^{-2}}.
  2. Find the acceleration interval corresponding to a 2π2\pi phase change.
  3. Explain why sensitivity and dynamic range must be reported together.
Solution

Using δϕ=keffδaT2\delta\phi=k_{\mathrm{eff}}\delta aT^2,

δϕ=(1.61×107)(1.0×10−8)(0.20)2=6.44×10−3 rad.\delta\phi = (1.61\times10^7) (1.0\times10^{-8}) (0.20)^2 = 6.44\times10^{-3}\ \mathrm{rad}.

One 2π2\pi period corresponds to

Δa2π=2πkeffT2≈9.76×10−6 m s−2.\Delta a_{2\pi} = \frac{2\pi} {k_{\mathrm{eff}}T^2} \approx 9.76\times10^{-6}\ \mathrm{m\,s^{-2}}.

Increasing TT amplifies a small acceleration but shrinks the phase-period acceleration interval as T−2T^{-2}. A cosine population measurement also needs a known branch or phase quadrature. A high-sensitivity interferometer may need a chirp, midfringe lock, auxiliary accelerometer, multiple scale factors, or phase-unwrapping model in a dynamic environment.

Two clock ratios r1r_1 and r2r_2 each have standard uncertainty u=3.0×10−18u=3.0\times10^{-18} and correlation coefficient ρ=0.60\rho=0.60 because they share an oscillator and link. Find the standard uncertainty of:

  1. the average (r1+r2)/2(r_1+r_2)/2;
  2. the difference r1−r2r_1-r_2.

Treat the ratios as small fractional deviations around nominal values.

Solution

The covariance is Cov⁡(r1,r2)=ρu2\operatorname{Cov}(r_1,r_2)=\rho u^2. Therefore

Var⁡(r1+r22)=u2+u2+2ρu24=u22(1+ρ).\operatorname{Var} \left( \frac{r_1+r_2}{2} \right) = \frac{u^2+u^2+2\rho u^2}{4} = \frac{u^2}{2}(1+\rho).

Thus

uavg=3.0×10−180.8≈2.68×10−18.u_{\mathrm{avg}} = 3.0\times10^{-18} \sqrt{0.8} \approx 2.68\times10^{-18}.

For the difference,

Var⁡(r1−r2)=2u2(1−ρ),\operatorname{Var}(r_1-r_2) = 2u^2(1-\rho),

so

udiff=3.0×10−180.8≈2.68×10−18.u_{\mathrm{diff}} = 3.0\times10^{-18} \sqrt{0.8} \approx 2.68\times10^{-18}.

The equal numerical values are specific to ρ=0.60\rho=0.60. Positive common-mode correlation weakens averaging but improves a difference relative to the uncorrelated case.

A magnetometer has one-sided white magnetic-noise amplitude SB1/2=15 fT/HzS_B^{1/2}=15\ \mathrm{fT}/\sqrt{\mathrm{Hz}} over the relevant band. For a known-phase sinusoidal signal spanning many periods, with unity calibrated response and an ideal matched estimator, estimate the statistical amplitude uncertainty after τ=400 s\tau=400\ \mathrm{s} using uB≃SB1/2/τu_B\simeq S_B^{1/2}/\sqrt{\tau}. Why is this not automatically a detection limit for an unknown-frequency signal?

Solution

The ideal uncertainty is

uB≃15 fT/Hz400 s=0.75 fT.u_B \simeq \frac{15\ \mathrm{fT}/\sqrt{\mathrm{Hz}}} {\sqrt{400\ \mathrm{s}}} = 0.75\ \mathrm{fT}.

An unknown-frequency search scans multiple frequencies, phases, and possibly time windows. It must account for the spectral window, nonwhite noise, calibration, line broadening, trials factor, and detection efficiency. Environmental lines can also occupy the same band. The one-template uncertainty is therefore not a global discovery threshold.