Skip to content

Rydberg Electrometry

A Rydberg electrometer infers an electric-field property from the large dipole response of a highly excited atom. In the most common vapor-cell implementation, two optical fields create Rydberg electromagnetically induced transparency (EIT), a radio-frequency field couples the addressed Rydberg state to a neighboring state, and probe transmission supplies the data. The familiar Autler–Townes formula can then convert a resolved spectral splitting into a local field amplitude.

That one-line conversion is the middle of a measurement, not its beginning or end. A complete result must say which field component and amplitude convention are estimated, how a laser scan maps onto an atomic energy interval, how the glass cell changes the incident field, which atoms and positions contribute, and how the recorded photons determine an uncertainty. These distinctions are especially important when the splitting is unresolved, the field is detuned, or the desired measurand is the field outside the cell.

This page is the canonical home for the end-to-end Rydberg-electrometry inverse problem: operating regimes, field-to-spectrum forward models, likelihoods, sensitivity and bandwidth, polarization reconstruction, cell-transfer corrections, and evidence standards. Rydberg Atoms Basics owns the high-nn spectrum, scaling laws, lifetimes, polarizabilities, and transition dipoles. Electromagnetically Induced Transparency owns dark states and optical propagation, while Autler–Townes Splitting owns the general dressed-state line-shape derivation. Sensing Case Studies retains the historical account of the 2012 vapor-cell demonstration and the subsequent evidence record.

An electric field is a vector-valued function of position and time. A useful measurement therefore begins with a contract such as

SR=(θ,A,C,M,T,E).\mathcal S_{\mathrm R} = (\theta,\mathcal A,\mathcal C, \mathcal M,\mathcal T,\mathcal E).

Here

  • θ\theta is the estimand, such as a local peak amplitude, polarization, carrier detuning, modulation quadrature, or external-field parameter;
  • A\mathcal A specifies the atomic species, resolved sublevels, temperature, density, and Rydberg preparation;
  • C\mathcal C contains the optical and radio-frequency controls, scan axis, timing, and local oscillator;
  • M\mathcal M is the photon-count, voltage, camera, or fluorescence likelihood;
  • T\mathcal T maps the field outside the apparatus to the field sampled by the atoms; and
  • E\mathcal E is the estimator, calibration chain, and uncertainty model.

For a nearly monochromatic field, write the real field as

E(r,t)=Re⁡[E0(r)e−iωt].\mathbf E(\mathbf r,t) = \operatorname{Re} \left[ \mathbf E_0(\mathbf r) e^{-i\omega t} \right].

E0\mathbf E_0 is a peak complex phasor. For a linearly polarized sinusoid, Erms=E0/2E_{\mathrm{rms}}=E_0/\sqrt 2. Autler–Townes Rabi frequencies below use the peak convention. A result quoted in V m−1\mathrm{V\,m^{-1}} without saying peak, RMS, or spectral-amplitude convention is incomplete.

Four-level Rydberg electrometer and inference chain

The usual four-level transducer, three amplitude regimes, and the inference chain from an external field to optical data. Atomic structure calibrates the local projected field; an external-field claim additionally requires the cell transfer function and spatial weighting.

Use a ladder ∣g⟩↔∣e⟩↔∣r⟩|g\rangle\leftrightarrow|e\rangle\leftrightarrow|r\rangle with probe and coupling Rabi coefficients Ωp\Omega_p and Ωc\Omega_c. A field near angular frequency ωμ\omega_\mu couples ∣r⟩|r\rangle to a second Rydberg state ∣s⟩|s\rangle with coefficient Ωμ\Omega_\mu. Define atom-minus-field detunings

Δp=ωeg−ωp,Δc=ωre−ωc,Δμ=ωsr−ωμ,\begin{aligned} \Delta_p &= \omega_{eg}-\omega_p,\\ \Delta_c &= \omega_{re}-\omega_c,\\ \Delta_\mu &= \omega_{sr}-\omega_\mu, \end{aligned}

and cumulative detunings

δ2=Δp+Δc,δ3=Δp+Δc+Δμ.\delta_2 = \Delta_p+\Delta_c, \qquad \delta_3 = \Delta_p+\Delta_c+\Delta_\mu.

Angular frequencies are used unless 2π2\pi is written explicitly. In a rotating-wave frame, one convenient Hamiltonian is

Hℏ=(0Ωp∗/200Ωp/2ΔpΩc∗/200Ωc/2δ2Ωμ∗/200Ωμ/2δ3).\frac{H}{\hbar} = \begin{pmatrix} 0 & \Omega_p^*/2 & 0 & 0\\ \Omega_p/2 & \Delta_p & \Omega_c^*/2 & 0\\ 0 & \Omega_c/2 & \delta_2 & \Omega_\mu^*/2\\ 0 & 0 & \Omega_\mu/2 & \delta_3 \end{pmatrix}.

The model is useful only after specifying the states hidden inside each label. Hyperfine and Zeeman sublevels, unwanted polarization components, nearby Rydberg manifolds, and off-resonant couplings can turn one four-state diagram into many interfering pathways.

For peak phasor E0\mathbf E_0 and a selected transition,

Ωμ=−⟨s∣d∣r⟩⋅E0ℏ.\Omega_\mu = -\frac{ \langle s|\mathbf d|r\rangle \mathbin{\cdot} \mathbf E_0 }{\hbar}.

Its magnitude is

∣Ωμ∣=deffE0ℏ,deff=∣⟨s∣d∣r⟩⋅ϵ^∣.|\Omega_\mu| = \frac{d_{\mathrm{eff}}E_0}{\hbar}, \qquad d_{\mathrm{eff}} = \left| \langle s|\mathbf d|r\rangle \mathbin{\cdot} \widehat{\boldsymbol\epsilon} \right|.

The sign and complex phase matter in a closed-loop or coherent-mixing experiment, even though an isolated splitting depends only on ∣Ωμ∣|\Omega_\mu|. The rotating-wave approximation also requires the counter-rotating and neighboring-transition corrections to remain small on the required uncertainty scale.

In the weak-probe limit, ρgg≃1\rho_{gg}\simeq1 and the steady optical coherence can be written as a continued fraction. With coherence half-widths γge\gamma_{ge}, γgr\gamma_{gr}, and γgs\gamma_{gs},

2ρgeΩp∗=iγge−iΔp+∣Ωc∣2/4γgr−iδ2+∣Ωμ∣2/4γgs−iδ3.\frac{2\rho_{ge}}{\Omega_p^*} = \frac{i}{ \gamma_{ge}-i\Delta_p + \dfrac{|\Omega_c|^2/4}{ \gamma_{gr}-i\delta_2 + \dfrac{|\Omega_\mu|^2/4}{ \gamma_{gs}-i\delta_3 } } }.

This expression displays the nested transduction: the RF field changes a Rydberg coherence, the coupling field maps that change onto the optical coherence, and propagation maps the coherence onto transmitted light. It is not a universal vapor-cell fit function. Saturation, optical pumping, additional levels, velocity changes during transit, and propagation through appreciable optical depth require optical Bloch or master-equation modeling.

For a dilute homogeneous slice, the probe susceptibility is proportional to ρge/Ωp∗\rho_{ge}/\Omega_p^* and the intensity transmission has the form

Tp(Δp)=exp⁡ ⁣[−kpLIm⁡χ(Δp)].T_p(\Delta_p) = \exp\!\left[ -k_p L\operatorname{Im}\chi(\Delta_p) \right].

If spectral bin jj records a photon count njn_j, a basic detector model is

nj∼Poisson⁡[μj(θ,λ)],n_j \sim \operatorname{Poisson} \left[ \mu_j(\theta,\boldsymbol\lambda) \right],

where λ\boldsymbol\lambda collects laser power, detunings, linewidths, background, gain, and other nuisance parameters. Fitting peak positions is one possible estimator; fitting this physical likelihood is more reliable near the resolution boundary.

Ignore the weak optical probe momentarily and diagonalize the RF-coupled Rydberg pair. Relative to ∣r⟩|r\rangle, its rotating-frame block can be written

Hrsℏ=(0Ωμ∗/2Ωμ/2Δμ).\frac{H_{rs}}{\hbar} = \begin{pmatrix} 0 & \Omega_\mu^*/2\\ \Omega_\mu/2 & \Delta_\mu \end{pmatrix}.

The dressed shifts are

λ±=Δμ2±12Δμ2+∣Ωμ∣2,\lambda_\pm = \frac{\Delta_\mu}{2} \pm \frac12 \sqrt{ \Delta_\mu^2+|\Omega_\mu|^2 },

so the eigenvalue gap is

Sω=Δμ2+∣Ωμ∣2.S_\omega = \sqrt{ \Delta_\mu^2+|\Omega_\mu|^2 }.

On resonance, Sω=∣Ωμ∣S_\omega=|\Omega_\mu|. If a justified spectral estimator returns the corresponding splitting Sν=Sω/(2π)S_\nu=S_\omega/(2\pi) in hertz, then

E0,atom=hSνdeff.E_{0,\mathrm{atom}} = \frac{hS_\nu}{d_{\mathrm{eff}}}.

This direct conversion applies to the local, projected, peak field when the RF is resonant and the observed peak separation is a valid estimator of the dressed gap. Off resonance,

E0,atom=ℏdeffSω2−Δμ2.E_{0,\mathrm{atom}} = \frac{\hbar}{d_{\mathrm{eff}}} \sqrt{S_\omega^2-\Delta_\mu^2}.

Near the resolution threshold, raw local maxima need not coincide with the real parts of the response poles. The full line shape should then determine Ωμ\Omega_\mu and Δμ\Delta_\mu jointly.

The atomic gap and the separation on a plotted laser-frequency axis are not always equal. Let a moving atom have signed axial wave vectors kpk_p and kck_c. Its optical detunings are

Δp(v)=Δp+kpv,Δc(v)=Δc+kcv.\Delta_p(v) = \Delta_p+k_pv, \qquad \Delta_c(v) = \Delta_c+k_cv.

When the coupling laser is scanned while a resonant probe selects atoms near v=0v=0, the coupling-axis separation directly gives

Δνc=∣Ωμ∣2π\Delta\nu_c = \frac{|\Omega_\mu|}{2\pi}

in the ideal resonant limit. If instead the probe is scanned and the coupling frequency is held fixed, the strongly absorbing velocity class approximately satisfies Δp(v)=0\Delta_p(v)=0. For counter-propagating beams this gives

δ2≃−kckpΔp,\delta_2 \simeq -\frac{k_c}{k_p}\Delta_p,

and therefore

Δνp≃νpνc∣Ωμ∣2π=λcλp∣Ωμ∣2π.\Delta\nu_p \simeq \frac{\nu_p}{\nu_c} \frac{|\Omega_\mu|}{2\pi} = \frac{\lambda_c}{\lambda_p} \frac{|\Omega_\mu|}{2\pi}.

Thus a probe-axis splitting must be divided by νp/νc=λc/λp\nu_p/\nu_c=\lambda_c/\lambda_p before applying the atomic conversion. This result follows from the stated scan and propagation conventions; a memorized wavelength factor can be inverted by changing either one. Finite linewidth, optical pumping, temperature, and multilevel structure can also modify the simple selected-velocity argument, so precision work should fit the velocity-averaged response.

A field does not become undetectable merely because two peaks cannot be resolved. It can change the EIT depth, width, curvature, or phase. For independent Poisson spectral bins, the classical Fisher information for a field parameter EE is

FE=∑j1μj(E)(∂μj(E)∂E)2.F_E = \sum_j \frac{1}{\mu_j(E)} \left( \frac{\partial\mu_j(E)}{\partial E} \right)^2.

With nuisance parameters λ\boldsymbol\lambda, the relevant bound comes from the inverse of the full Fisher matrix,

Var⁡(E^)≥[F−1]EE,\operatorname{Var}(\widehat E) \geq \left[F^{-1}\right]_{EE},

not from 1/FE1/F_E while pretending that laser detuning and contrast are known. The best probe detuning is usually where the expected signal changes steeply and reproducibly, not necessarily at line center. Modulation and balanced or homodyne optical readout can move the estimate away from low-frequency laser noise, but their duty cycle and reference noise belong in the comparison.

If the observable is even in an unknown field, its small-field response may begin as

y(E)=y(0)+aE2+O(E4).y(E) = y(0)+aE^2+O(E^4).

Then ∂Ey=0\partial_Ey=0 at E=0E=0: the sign or phase is absent, and the local Fisher information can vanish even with a visibly perturbed spectrum at larger fields. A known bias or coherent local oscillator restores a linear term.

Let a strong local oscillator and a weak signal address the same effective Rydberg transition:

ELO(t)=ELcos⁡(ωLt),Es(t)=Escos⁡(ωst+ϕ).\begin{aligned} E_{\mathrm{LO}}(t) &= E_L\cos(\omega_Lt),\\ E_{\mathrm s}(t) &= E_s\cos(\omega_st+\phi). \end{aligned}

For Es≪ELE_s\ll E_L and ωs−ωL=ωIF\omega_s-\omega_L=\omega_{\mathrm{IF}} within the atomic and detector bandwidth, the slowly varying envelope is

Eenv(t)≃EL+Escos⁡(ωIFt+ϕ).E_{\mathrm{env}}(t) \simeq E_L + E_s\cos(\omega_{\mathrm{IF}}t+\phi).

Linearizing an optical observable y=f(Eenv)y=f(E_{\mathrm{env}}) gives

y(t)≃f(EL)+f′(EL)Escos⁡(ωIFt+ϕ).y(t) \simeq f(E_L) + f'(E_L)E_s \cos(\omega_{\mathrm{IF}}t+\phi).

The atoms now down-convert the signal to an intermediate-frequency optical beat note. Its amplitude estimates a signal quadrature and its phase is defined relative to the local oscillator. This extends detection below the resolved Autler–Townes regime and enables coherent reception. The result is not reference-free: local-oscillator amplitude, phase noise, spatial overlap, and detuning are resources and nuisance parameters.

Rydberg electrometry spans several physical responses. They should not be collapsed into one calibration formula.

RegimeDominant observableLeading field dependenceInformation lost without a reference
resonant RF, resolvedAutler–Townes gapSν∝E0S_\nu\propto E_0absolute carrier phase
resonant RF, unresolvedEIT line shape or mixer beatmodel dependent; linear with a coherent biasphase or sign in an intensity-only measurement
far-detuned RFdifferential AC Stark shiftΔν∝E02\Delta\nu\propto E_0^2sign and generally carrier phase
static or quasistaticdifferential DC Stark shiftΔν∝Edc2\Delta\nu\propto E_{\mathrm{dc}}^2 absent a biasfield sign

For a nondegenerate state with scalar static polarizability αi(0)\alpha_i(0),

ΔUi(dc)=−12αi(0)Edc2.\Delta U_i^{(\mathrm{dc})} = -\frac12\alpha_i(0)E_{\mathrm{dc}}^2.

The measured transition shift is set by the differential polarizability,

Δνsr(dc)=−αs(0)−αr(0)2hEdc2.\Delta\nu_{sr}^{(\mathrm{dc})} = -\frac{ \alpha_s(0)-\alpha_r(0) }{2h} E_{\mathrm{dc}}^2.

For a far-off-resonant sinusoid with peak amplitude E0E_0, the cycle average introduces another factor of 1/21/2:

ΔUi(ac)=−14αi(ω)E02.\Delta U_i^{(\mathrm{ac})} = -\frac14 \alpha_i(\omega)E_0^2.

Tensor polarizability, near-degenerate mixing, and strong fields require a Stark-map or Floquet calculation rather than these scalar formulas. The canonical derivations are on Stark Effect in Atoms and Dynamic Polarizability.

A known static bias EbE_b can linearize a small perturbation ee:

(Eb+e)2=Eb2+2Ebe+e2.(E_b+e)^2 = E_b^2+2E_be+e^2.

It also imports bias calibration and drift. In glass alkali-vapor cells, surface charge can screen low-frequency external fields or create inhomogeneous offsets, so the field at the atoms must be measured or modeled rather than inferred from electrode voltage alone.

Choose a quantization axis and decompose the RF phasor into spherical components,

Eq=eq∗⋅E0,q∈{−1,0,+1}.E_q = \mathbf e_q^* \mathbin{\cdot} \mathbf E_0, \qquad q\in\{-1,0,+1\}.

For one angular-momentum pathway,

Ωmm′(q)=−Eqℏ⟨F′m′∣dq∣Fm⟩,\Omega_{m m'}^{(q)} = -\frac{E_q}{\hbar} \langle F'm'|d_q|Fm\rangle,

with

⟨F′m′∣dq∣Fm⟩=(−1)F′−m′(F′1F−m′qm)⟨F′∥d∥F⟩.\langle F'm'|d_q|Fm\rangle = (-1)^{F'-m'} \begin{pmatrix} F'&1&F\\ -m'&q&m \end{pmatrix} \langle F'\|d\|F\rangle.

Selection rules therefore make the spectrum a polarization-dependent sum of pathways, not a direct meter of ∣E∣|\mathbf E|. Optical pumping changes their weights. Some populated states may be uncoupled and leave an unsplit central feature; others have different dipole moments and broaden or resolve the doublet.

A vector measurement needs enough independent settings to identify the desired phasor parameters. Schematically,

y=f(E0,B0,ϵp,ϵc,λ).\mathbf y = \mathbf f \left( \mathbf E_0, \mathbf B_0, \boldsymbol\epsilon_p, \boldsymbol\epsilon_c, \boldsymbol\lambda \right).

Rotating a known quantization axis, changing optical polarizations, resolving Zeeman components, or using multiple transitions can provide independent constraints. The Jacobian with respect to the estimated field components must have full rank after nuisance parameters are included. Measuring three peak separations is not automatically a three-dimensional vector measurement if the same unknown polarization impurity affects all three.

For velocity distribution f(v)f(v), the susceptibility of a dilute thermal ensemble begins with

χ‾(Δp,Δc)=∫−∞∞f(v)χ[Δp+kpv,Δc+kcv] dv.\overline\chi(\Delta_p,\Delta_c) = \int_{-\infty}^{\infty} f(v) \chi \left[ \Delta_p+k_pv, \Delta_c+k_cv \right] \,dv.

Counter-propagating probe and coupling beams reduce the residual two-photon Doppler wave vector kp+kck_p+k_c, but unequal wavelengths prevent perfect cancellation in the common two-photon ladder. Beam angle, transit time, collisions, laser linewidth, and velocity-changing processes further modify the response. Three-photon ladders can be arranged with a smaller net wave vector, trading additional lasers and levels for narrower Doppler response.

At carrier frequency ω\omega, write the local field as

Eatom(r,ω)=Tcell(r,ω)Eexternal(ω).\mathbf E_{\mathrm{atom}}(\mathbf r,\omega) = \mathsf T_{\mathrm{cell}}(\mathbf r,\omega) \mathbf E_{\mathrm{external}}(\omega).

Tcell\mathsf T_{\mathrm{cell}} can include dielectric boundary conditions, standing waves, nearby conductors, electrodes, supports, resonators, and polarization rotation. At low frequency it can include charge transport and screening. It is generally a complex, position-dependent tensor rather than a single scalar correction.

The detected spectrum is weighted by probe and coupling intensities and by atomic trajectories. A schematic forward model is

μj=Dj{∫d3r dv W(r,v)P[χj(Eatom,v,r)]},\mu_j = \mathcal D_j \left\{ \int d^3r\,dv\, W(\mathbf r,v) \mathcal P \left[ \chi_j (\mathbf E_{\mathrm{atom}},v,\mathbf r) \right] \right\},

where P\mathcal P denotes optical propagation and Dj\mathcal D_j the detector response. Because splitting, absorption, propagation, and detector averaging are nonlinear, one generally cannot replace the inhomogeneous field by its spatial mean. A field gradient may appear as broadening, unequal peaks, or a multicomponent line shape.

This distinction separates two legitimate measurands:

  1. Atomic local field: a weighted field experienced by atoms in the illuminated cell volume.
  2. External or incident field: a parameter before the cell, inferred only through a validated electromagnetic and spatial-transfer model.

A subwavelength optical spot can localize the readout more finely than the RF wavelength, but spatial resolution is set by the full weighting kernel, including beam waist, depth of field, atomic motion, camera pixels, and inverse regularization. It is not simply the cell size.

For total averaging time τ\tau in a stationary white-noise regime, define

σE(τ)≃ηEτ,\sigma_E(\tau) \simeq \frac{\eta_E}{\sqrt\tau},

where ηE\eta_E has units V m−1 Hz−1/2\mathrm{V\,m^{-1}\,Hz^{-1/2}}. The reported value should identify:

  • peak or RMS field;
  • local or external field;
  • carrier frequency, detuning, and polarization;
  • modulation or analysis frequency and resolution bandwidth;
  • active volume and spatial weighting;
  • optical powers, atomic density, and readout protocol;
  • whether dead time and initialization are included; and
  • the averaging interval over which 1/τ1/\sqrt\tau behavior was verified.

For a resolved resonant splitting,

σE=hdeffσSν,\sigma_E = \frac{h}{d_{\mathrm{eff}}} \sigma_{S_\nu},

before scan-axis, detuning, and transfer uncertainties. This relation converts a splitting-estimator uncertainty; it does not predict that uncertainty. Photon shot noise, laser amplitude and frequency noise, atomic projection or number fluctuations, thermal radiation, detector noise, and model mismatch all enter through the data likelihood.

The word bandwidth has at least three meanings here.

  1. Carrier coverage is the set of RF frequencies addressable by available Rydberg transitions or by Stark, Zeeman, or dressing-based tuning. Broad coverage across configurations is not simultaneous broadband reception.
  2. Instantaneous modulation bandwidth is the highest envelope or phase variation faithfully followed at one carrier setting. It depends on optical pumping, Rabi frequencies, coherence, transit replacement, and detector electronics.
  3. Measurement bandwidth or resolution bandwidth is the noise-equivalent bandwidth used to turn fluctuations into a spectral sensitivity.

For a first-order response with time constant τR\tau_R,

H(ωm)=H(0)1+iωmτR,f3dB=12πτR.H(\omega_m) = \frac{H(0)}{1+i\omega_m\tau_R}, \qquad f_{3\mathrm{dB}} = \frac{1}{2\pi\tau_R}.

In a thermal vapor, the transit estimate

τtr∼wv⊥\tau_{\mathrm{tr}} \sim \frac{w}{v_\perp}

can be comparable to or shorter than internal relaxation times. Smaller beams can increase response speed but reduce atom number and interaction length, so bandwidth and sensitivity are jointly engineered.

At low amplitude, an Autler–Townes estimator fails when the line shape no longer identifies two branches. Line-shape, modulation, or mixer estimators can continue below that point. At high amplitude, the selected pair is no longer isolated: nearby Rydberg levels, multiple RF photons, spatially varying splittings, counter-rotating terms, ionization, and detector scan range can limit the model. A trustworthy dynamic range is the interval over which linearity, identifiability, and coverage have all been validated.

The resonant local-field chain is attractive because hh is exact in the SI and a transition dipole can be calculated from atomic structure:

Sν⟶∣Ωμ∣⟶E0,atom.S_\nu \longrightarrow |\Omega_\mu| \longrightarrow E_{0,\mathrm{atom}}.

It is better described as atom-referenced than as automatically self-calibrating. The transition identity, polarization projection, scan-axis scale, RF detuning, and line-shape estimator must still be established. To claim an external field, append

E0,atom⟶Tcell−1⟶E0,external.E_{0,\mathrm{atom}} \longrightarrow \mathsf T_{\mathrm{cell}}^{-1} \longrightarrow E_{0,\mathrm{external}}.

For an approximately scalar transfer magnitude TT and independent small uncertainties, a useful first budget is

(u(Eexternal)Eexternal)2≃(u(Sν)Sν)2+(u(D)D)2+(u(deff)deff)2+(u(T)T)2+umodel2.\left( \frac{u(E_{\mathrm{external}})}{E_{\mathrm{external}}} \right)^2 \simeq \left(\frac{u(S_\nu)}{S_\nu}\right)^2 + \left(\frac{u(D)}{D}\right)^2 + \left(\frac{u(d_{\mathrm{eff}})}{d_{\mathrm{eff}}}\right)^2 + \left(\frac{u(T)}{T}\right)^2 +u_{\mathrm{model}}^2.

DD denotes the applicable scan-axis factor. Correlations require covariance terms, and line-shape model discrepancy is often not well represented by a single Gaussian contribution. A robust analysis varies beam position and power, reverses or rotates polarizations, changes scan direction, fits alternative plausible models, maps the cell electromagnetically, and reports which corrections were measured rather than assumed.

Precision and accuracy answer different questions. Repeated spectra can give a small σE\sigma_E while an incorrect cell transfer, pathway dipole, or frequency scale biases every result. Conversely, a traceable calibration does not guarantee high sensitivity or communications-ready bandwidth.

Suppose a coupling-laser scan resolves a resonant splitting

Sν=1.000 MHzS_\nu = 1.000\ \mathrm{MHz}

for an isolated pathway with

deff=1000 ea0=8.4784×10−27 C m.d_{\mathrm{eff}} = 1000\,e a_0 = 8.4784\times10^{-27}\ \mathrm{C\,m}.

The local peak field is

E0,atom=hSνdeff=7.815×10−2 V m−1=0.7815 mV cm−1.\begin{aligned} E_{0,\mathrm{atom}} &= \frac{hS_\nu}{d_{\mathrm{eff}}}\\ &= 7.815\times10^{-2}\ \mathrm{V\,m^{-1}}\\ &= 0.7815\ \mathrm{mV\,cm^{-1}}. \end{aligned}

For a linearly polarized sinusoid,

Erms,atom=5.526×10−2 V m−1.E_{\mathrm{rms,atom}} = 5.526\times10^{-2}\ \mathrm{V\,m^{-1}}.

If the same atomic gap were read by scanning a 780 nm780\ \mathrm{nm} probe while a counter-propagating 480 nm480\ \mathrm{nm} coupling laser remained fixed, the ideal probe-axis separation would instead be

Δνp=480780(1.000 MHz)≃0.615 MHz.\Delta\nu_p = \frac{480}{780} (1.000\ \mathrm{MHz}) \simeq 0.615\ \mathrm{MHz}.

Applying hΔνp/deffh\Delta\nu_p/d_{\mathrm{eff}} directly would underestimate the field. Finally, if a validated scalar cell model gives T=0.80T=0.80, the inferred external peak field is E0,external=0.0977 V m−1E_{0,\mathrm{external}}=0.0977\ \mathrm{V\,m^{-1}}. The local atomic result does not require that last correction; the external result does.

GoalUseful modeMain calibration objectDominant failure mode
absolute local RF amplituderesonant resolved splittingdeffd_{\mathrm{eff}}, detuning, scan axismultilevel or inhomogeneous peaks
very weak coherent signalRF homodyne or heterodyneLO amplitude, phase, transfer gainreference noise and drift
polarizationresolved sublevels and varied axesangular pathways and optical polarizationrank-deficient reconstruction
off-resonant carrierAC Stark or dressed-state modeldynamic polarizability or Floquet spectrumquadratic ambiguity and nearby levels
low-frequency or DC fieldDC Stark shift with bias or electrodesdifferential polarizability and screeningwall charge and field inhomogeneity
field imagecamera or scanned optical volumespatial kernel and cell mapinterpreting pixel size as resolution
communication waveformlocked optical readout and mixertransfer function and linear rangeconfusing data rate with carrier coverage

No mode is universally best. Resolved splitting provides a transparent frequency-to-field conversion but has a linewidth-imposed lower boundary. Mixing provides weak-field and phase sensitivity but consumes a reference. Stark sensing accesses detuned or low-frequency fields but is quadratic unless biased. The estimand decides which compromise is appropriate.

The strongest claim supported by a result depends on what was actually validated.

  • Observed response: a reproducible optical change correlated with an applied field.
  • Calibrated local amplitude: a state-resolved model, frequency scale, polarization, and estimator connect that response to the field at the atoms.
  • External-field measurement: the cell and spatial-transfer model are independently validated over the relevant frequency and geometry.
  • Sensitivity: noise spectra or repeated estimates establish a stated amplitude uncertainty per square root bandwidth under fixed conditions.
  • Receiver performance: linearity, phase or modulation fidelity, dynamic range, instantaneous bandwidth, and error rate are measured.
  • Quantum advantage: a matched classical receiver with the same aperture, power, bandwidth, prior information, and estimation loss is outperformed.

Rydberg atoms are quantum systems, but using them as transducers is not by itself evidence of quantum-enhanced metrology. Large atomic dipoles, subwavelength optical readout, broad reconfigurable carrier coverage, and an atom-referenced local calibration are meaningful capabilities without that extra claim.

Calling every peak distance a Rabi frequency

Section titled “Calling every peak distance a Rabi frequency”

Near threshold, detuning, unequal strengths, overlapping sublevels, and inhomogeneous fields move local maxima away from the dressed-state gap. Fit and test the physical response.

The standard Rabi coefficient uses the peak sinusoidal field. Quoting an RMS result without the factor 2\sqrt2, or mixing field and power conventions, creates a systematic calibration error.

Applying a wavelength factor without naming the scan

Section titled “Applying a wavelength factor without naming the scan”

The probe-scan correction and coupling-scan calibration differ. Derive the mapping from signed Doppler detunings for the actual beam geometry.

Treating one radial matrix element as the full dipole

Section titled “Treating one radial matrix element as the full dipole”

Angular, hyperfine, Zeeman, and polarization factors determine deffd_{\mathrm{eff}}. Unresolved pathways can make a single effective value model dependent.

Calling the incident field the atomic field

Section titled “Calling the incident field the atomic field”

Glass, standing waves, charges, supports, and resonators transform the field. An atomic measurement is local until that transfer has been included.

Averaging the field before calculating the spectrum

Section titled “Averaging the field before calculating the spectrum”

The line shape is nonlinear in local coupling and optical depth. Spatially average the physical response in the correct propagation order.

Reporting carrier range as instantaneous bandwidth

Section titled “Reporting carrier range as instantaneous bandwidth”

Switching among transitions can cover many decades of carrier frequency while one operating point follows only a finite modulation band.

Equating atomic operation with quantum advantage

Section titled “Equating atomic operation with quantum advantage”

An atom-referenced transducer may be excellent without beating a resource-matched non-atomic receiver. State the comparison and loss function.

For a resonant, isolated RF transition and peak field convention,

∣Ωμ∣=deffE0,atomℏ,E0,atom=hSνdeff.|\Omega_\mu| = \frac{d_{\mathrm{eff}}E_{0,\mathrm{atom}}}{\hbar}, \qquad E_{0,\mathrm{atom}} = \frac{hS_\nu}{d_{\mathrm{eff}}}.

For RF detuning Δμ\Delta_\mu, the dressed angular-frequency gap is

Sω=Δμ2+∣Ωμ∣2.S_\omega = \sqrt{\Delta_\mu^2+|\Omega_\mu|^2}.

In the ideal counter-propagating probe-scan geometry,

Δνp≃νpνc∣Ωμ∣2π,\Delta\nu_p \simeq \frac{\nu_p}{\nu_c} \frac{|\Omega_\mu|}{2\pi},

whereas a coupling-laser scan reads the atomic gap directly. The local and external fields are related by

Eatom(r,ω)=Tcell(r,ω)Eexternal(ω).\mathbf E_{\mathrm{atom}}(\mathbf r,\omega) = \mathsf T_{\mathrm{cell}}(\mathbf r,\omega) \mathbf E_{\mathrm{external}}(\omega).

These equations become a metrological result only when the scan, polarization, spatial weighting, likelihood, bandwidth, and uncertainty budget are stated.

  1. T. F. Gallagher, Rydberg Atoms (Cambridge University Press, 1994).
  2. A. K. Mohapatra, T. R. Jackson, and C. S. Adams, “Coherent optical detection of highly excited Rydberg states using electromagnetically induced transparency,” Physical Review Letters 98, 113003 (2007).
  3. J. A. Sedlacek, A. Schwettmann, H. Kübler, R. Löw, T. Pfau, and J. P. Shaffer, “Microwave electrometry with Rydberg atoms in a vapour cell using bright atomic resonances,” Nature Physics 8, 819–824 (2012).
  4. J. A. Sedlacek, A. Schwettmann, H. Kübler, and J. P. Shaffer, “Atom-based vector microwave electrometry using rubidium Rydberg atoms in a vapor cell,” Physical Review Letters 111, 063001 (2013).
  5. C. L. Holloway et al., “Broadband Rydberg atom-based electric-field probe for SI-traceable, self-calibrated measurements,” IEEE Transactions on Antennas and Propagation 62, 6169–6182 (2014).
  6. H. Q. Fan et al., “Subwavelength microwave electric-field imaging using Rydberg atoms inside atomic vapor cells,” Optics Letters 39, 3030–3033 (2014).
  7. H. Fan et al., “Effect of vapor-cell geometry on Rydberg-atom-based measurements of radio-frequency electric fields,” Physical Review Applied 4, 044015 (2015).
  8. S. Kumar, H. Fan, H. Kübler, J. Sheng, and J. P. Shaffer, “Atom-based sensing of weak radio frequency electric fields using homodyne readout,” Scientific Reports 7, 42981 (2017).
  9. K. C. Cox, D. H. Meyer, F. K. Fatemi, and P. D. Kunz, “Quantum-limited atomic receiver in the electrically small regime,” Physical Review Letters 121, 110502 (2018).
  10. M. T. Simons, A. H. Haddab, J. A. Gordon, and C. L. Holloway, “A Rydberg atom-based mixer: Measuring the phase of a radio frequency wave,” Applied Physics Letters 114, 114101 (2019).
  11. J. A. Gordon, C. L. Holloway, M. T. Simons, and A. H. Haddab, “Weak electric-field detection with sub-1 Hz resolution at radio frequencies using a Rydberg atom-based mixer,” AIP Advances 9, 045030 (2019).
  12. M. Jing et al., “Atomic superheterodyne receiver based on microwave-dressed Rydberg spectroscopy,” Nature Physics 16, 911–915 (2020).
  13. D. H. Meyer, Z. A. Castillo, K. C. Cox, and P. D. Kunz, “Assessment of Rydberg atoms for wideband electric field sensing,” Journal of Physics B 53, 034001 (2020).
  14. Y.-Y. Jau and T. Carter, “Vapor-cell-based atomic electrometry for detection frequencies below 1 kHz,” Physical Review Applied 13, 054034 (2020).
  15. S. M. Bohaichuk et al., “Origins of Rydberg-atom electrometer transient response and its impact on radio-frequency pulse sensing,” Physical Review Applied 18, 034030 (2022).
  16. N. Schlossberger et al., “Rydberg states of alkali atoms in atomic vapour as SI-traceable field probes and communications receivers,” Nature Reviews Physics 6, 606–620 (2024).
  17. N. Schlossberger et al., “Zeeman-resolved Autler–Townes splitting in Rydberg atoms with tunable resonances and a single transition dipole moment,” Physical Review A 109, L021702 (2024).
  18. B. N. Miller, D. H. Meyer, T. Virtanen, C. M. O’Brien, and K. C. Cox, “RydIQule: A graph-based paradigm for modeling Rydberg and atomic sensors,” Computer Physics Communications 294, 108952 (2024).

A coupling-laser scan resolves a resonant splitting of 2.40 MHz2.40\ \mathrm{MHz}. The effective peak-field dipole is 750 ea0750\,e a_0. Find the local peak and RMS fields.

Solution

Using ea0=8.47835×10−30 C me a_0=8.47835\times10^{-30}\ \mathrm{C\,m},

deff=6.3588×10−27 C m.d_{\mathrm{eff}} = 6.3588\times10^{-27}\ \mathrm{C\,m}.

Therefore

E0=h(2.40×106 Hz)deff=0.2501 V m−1.E_0 = \frac{h(2.40\times10^6\ \mathrm{Hz})}{d_{\mathrm{eff}}} = 0.2501\ \mathrm{V\,m^{-1}}.

For a linear sinusoid,

Erms=E02=0.1768 V m−1.E_{\mathrm{rms}} = \frac{E_0}{\sqrt2} = 0.1768\ \mathrm{V\,m^{-1}}.

An RF-dressed spectrum gives Sω/(2π)=8.0 MHzS_\omega/(2\pi)=8.0\ \mathrm{MHz} while an independent frequency measurement gives ∣Δμ∣/(2π)=4.8 MHz|\Delta_\mu|/(2\pi)=4.8\ \mathrm{MHz}. Find ∣Ωμ∣/(2π)|\Omega_\mu|/(2\pi).

Solution

The generalized gap obeys

Sω2=Δμ2+∣Ωμ∣2.S_\omega^2 = \Delta_\mu^2+|\Omega_\mu|^2.

The common factor 2π2\pi cancels, so

∣Ωμ∣2π=(8.0 MHz)2−(4.8 MHz)2=6.4 MHz.\frac{|\Omega_\mu|}{2\pi} = \sqrt{(8.0\ \mathrm{MHz})^2-(4.8\ \mathrm{MHz})^2} = 6.4\ \mathrm{MHz}.

Using 8.0 MHz8.0\ \mathrm{MHz} directly would overestimate the field by 25%25\%.

Take counter-propagating beams with signed wave vectors kp>0k_p>0 and kc<0k_c<0. The coupling laser is fixed on resonance. Use the velocity-class condition Δp+kpv=0\Delta_p+k_pv=0 to express the two-photon detuning in terms of the scanned probe detuning.

Solution

The selected velocity is

v=−Δpkp.v = -\frac{\Delta_p}{k_p}.

With Δc=0\Delta_c=0,

δ2(v)=Δp+(kp+kc)v=Δp−kp+kckpΔp=−kckpΔp.\begin{aligned} \delta_2(v) &= \Delta_p+(k_p+k_c)v\\ &= \Delta_p - \frac{k_p+k_c}{k_p}\Delta_p\\ &= -\frac{k_c}{k_p}\Delta_p. \end{aligned}

Because −kc/kp=∣kc∣/kp=νc/νp-k_c/k_p=|k_c|/k_p=\nu_c/\nu_p, a two-photon separation SωS_\omega appears on the probe axis as

Δωp=νpνcSω.\Delta\omega_p = \frac{\nu_p}{\nu_c}S_\omega.

A prepared state has only a π\pi-allowed RF transition relative to the bias axis. What does an Autler–Townes measurement reveal about a purely transverse field?

Solution

A π\pi transition couples to the spherical component E0E_0, parallel to the quantization axis. A purely transverse field has only q=±1q=\pm1 components, so the chosen transition has Ωμ=0\Omega_\mu=0 in the ideal selection-rule model. The absence of splitting is a polarization blind spot, not evidence that the total electric field vanishes. Additional σ+\sigma^+ and σ−\sigma^- pathways, or a rotated quantization axis, are needed for reconstruction.

At one operating point the mean detected count is μ(E)=μ0(1+aE)\mu(E)=\mu_0(1+aE) over the local linear range. Assuming Poisson counts, find the single-shot Fisher information and its small-field limit.

Solution

For a Poisson mean,

FE=1μ(E)(dμdE)2.F_E = \frac{1}{\mu(E)} \left( \frac{d\mu}{dE} \right)^2.

Since dμ/dE=μ0ad\mu/dE=\mu_0a,

FE=μ0a21+aE.F_E = \frac{\mu_0a^2}{1+aE}.

At ∣aE∣≪1|aE|\ll1,

FE≃μ0a2,σE≥1∣a∣μ0.F_E \simeq \mu_0a^2, \qquad \sigma_E \geq \frac{1}{|a|\sqrt{\mu_0}}.

This is a photon-shot-noise bound for the assumed likelihood, not a complete instrument sensitivity.

6. Show how a local oscillator restores linear response

Section titled “6. Show how a local oscillator restores linear response”

Suppose the optical observable is y(E)=y0+aE2y(E)=y_0+aE^2. Set E(t)=EL+Escos⁡(ωIFt+ϕ)E(t)=E_L+E_s\cos(\omega_{\mathrm{IF}}t+\phi) with Es≪ELE_s\ll E_L. Find the component at ωIF\omega_{\mathrm{IF}}.

Solution

Expanding to first order in EsE_s gives

y(t)=y0+a[EL+Escos⁡(ωIFt+ϕ)]2≃y0+aEL2+2aELEscos⁡(ωIFt+ϕ).\begin{aligned} y(t) &= y_0+a \left[ E_L+E_s\cos(\omega_{\mathrm{IF}}t+\phi) \right]^2\\ &\simeq y_0+aE_L^2 +2aE_LE_s \cos(\omega_{\mathrm{IF}}t+\phi). \end{aligned}

The intermediate-frequency amplitude is 2∣a∣ELEs2|a|E_LE_s, linear in the weak signal. Its phase is measured relative to the local oscillator.

7. Propagate an external-field uncertainty

Section titled “7. Propagate an external-field uncertainty”

A resonant measurement has independent relative standard uncertainties of 0.8%0.8\% from splitting, 0.4%0.4\% from the scan factor, 0.6%0.6\% from the dipole, and 2.5%2.5\% from the cell transfer. Ignore model discrepancy. Find the combined relative uncertainty.

Solution

Add the independent relative variances:

u(E)E=(0.008)2+(0.004)2+(0.006)2+(0.025)2=0.0272.\begin{aligned} \frac{u(E)}{E} &= \sqrt{ (0.008)^2+(0.004)^2 +(0.006)^2+(0.025)^2 }\\ &= 0.0272. \end{aligned}

The combined standard uncertainty is about 2.7%2.7\%, dominated by the cell transfer. Improving the spectral fit alone would have little effect.

An experiment reports a precisely measured Autler–Townes splitting and a calculated radial Rydberg matrix element, then calls the incident microwave field self-calibrated. List four additional demonstrations needed to support that wording.

Solution

A defensible incident-field claim should at least demonstrate:

  1. the populated hyperfine and Zeeman pathways and their angular dipole factors;
  2. polarization at the atoms and the peak-versus-RMS convention;
  3. the scan-axis frequency calibration, RF detuning, and validity of the splitting estimator; and
  4. the position-dependent transfer from the incident field through the cell, including standing waves and optical spatial weighting.

It should also test line-shape alternatives, quantify laser and detector noise, and propagate all terms into an uncertainty budget. Without the cell transfer, the result can still be a well-calibrated local atomic field.