Rydberg Atoms

A Rydberg platform begins with an ordinary atom and temporarily promotes one electron to a highly excited state. The promotion is brief, but it changes the relevant scales dramatically: electric response and interatomic interactions become large, nearby levels become dense, and weak stray fields become experimentally important. The useful resource is therefore not merely a large principal quantum number. It is a repeatable cycle of state preparation, coherent excitation, interaction, de-excitation, and calibrated measurement.

The central calculation is a driven two-state reduction of a detuned three-level ladder. It predicts the useful coupling together with the light shifts and intermediate-state scattering that limit it. Motion, state-dependent trapping, and imperfect detection then determine how that Hamiltonian appears in an experiment.

Required background. Rydberg atoms basics supplies quantum defects and pair-interaction scaling, optical tweezers supplies confinement and preparation, Rabi oscillations supplies coherent transfer, and adiabatic elimination supplies the controlled reduction used below. Rydberg blockade develops interaction-suppressed double excitation and its collective dynamics.

From an Excited State to a Platform

The same large response that makes a Rydberg state useful makes it fragile. Choosing a state is therefore a multi-objective decision rather than a search for the largest available $n$.

Design question	Why it matters	Typical diagnostic
Which isotope and internal ground state?	Fixes hyperfine structure, available optical paths, and qubit encoding	microwave and optical spectroscopy
Which $n\ell j m_j$ or hyperfine-resolved Rydberg state?	Sets lifetime, polarizability, angular interaction channels, and spectral isolation	Stark map and polarization-resolved spectrum
One-photon or multiphoton excitation?	Sets laser wavelength, intermediate-state scattering, Doppler sensitivity, and phase noise	single-atom Rabi and Ramsey data
Is the trap left on?	Determines differential light shifts, photoionization, and mechanical forces	spectroscopy versus trap depth and pulse timing
How is the Rydberg population detected?	Determines whether absence means Rydberg excitation, ordinary loss, or another dark state	independently calibrated readout matrix
Which pair potential is used?	Determines whether a scalar $C_6/R^6$ model is valid	pair spectroscopy versus distance and angle

High $n$ often strengthens interactions, but it also compresses the spectrum, increases sensitivity to electric fields, and can bring unwanted pair channels near resonance. Low orbital angular momentum is usually convenient for optical access, while high-angular-momentum and circular states can have very different lifetime and control requirements. Alkali atoms offer relatively simple one-valence-electron spectra. Alkaline-earth-like atoms add metastable clock states and autoionization-based readout, but also add electronic structure that must be modeled explicitly.

Optical Excitation

A three-level ladder

Many alkali-atom experiments use a ladder

$$|g\rangle \xleftrightarrow{\ \Omega_1\ } |p\rangle \xleftrightarrow{\ \Omega_2\ } |r\rangle,$$

where $|p\rangle$ is a short-lived intermediate state and $|r\rangle$ is the selected Rydberg state. Let

$$\begin{aligned} \Delta &= \omega_1-\omega_{pg}, \\ \delta &= \omega_1+\omega_2-\omega_{rg} \end{aligned}$$

be the one-photon and two-photon detunings. With the rotating-wave approximation and basis $(|g\rangle,|p\rangle,|r\rangle)$, one useful convention is

$$\frac{H}{\hbar} = \begin{pmatrix} 0 & \Omega_1^*/2 & 0 \\ \Omega_1/2 & -\Delta & \Omega_2^*/2 \\ 0 & \Omega_2/2 & -\delta \end{pmatrix}.$$

The complex phases of $\Omega_1$ and $\Omega_2$ contain laser phases and spatial phase factors. A different sign convention for detunings changes signs in the effective Hamiltonian, so every quoted light shift must be read with its convention.

Adiabatic elimination

If

$$|\Delta| \gg |\Omega_1|,\ |\Omega_2|,\ \Gamma_p,\ |\delta|,$$

and the controls vary slowly enough, the intermediate amplitude follows the long-lived amplitudes approximately:

$$c_p \simeq \frac{ \Omega_1 c_g+\Omega_2^*c_r }{ 2\Delta }.$$

Substitution gives an effective two-state Hamiltonian

$$\frac{H_{\mathrm{eff}}}{\hbar} \simeq \begin{pmatrix} \Delta_g^{\mathrm{AC}} & \Omega_{\mathrm{eff}}^*/2 \\ \Omega_{\mathrm{eff}}/2 & -\delta+\Delta_r^{\mathrm{AC}} \end{pmatrix},$$

with

$$\begin{aligned} \Omega_{\mathrm{eff}} &\simeq \frac{\Omega_1\Omega_2}{2\Delta}, \\ \Delta_g^{\mathrm{AC}} &\simeq \frac{|\Omega_1|^2}{4\Delta}, \\ \Delta_r^{\mathrm{AC}} &\simeq \frac{|\Omega_2|^2}{4\Delta}. \end{aligned}$$

Thus the experimentally relevant two-photon detuning is not just $\delta$. In this convention,

$$\delta_{\mathrm{eff}} = \delta +\Delta_g^{\mathrm{AC}} -\Delta_r^{\mathrm{AC}},$$

up to additional Zeeman, trap, many-body, and laser-frequency shifts. The differential light shift can exceed the effective Rabi frequency even when the intermediate state is only weakly populated.

Adiabatic elimination is not an identity. It can fail because the detuning is too small, because a pulse edge has high-frequency content, because nearby hyperfine levels contribute with different signs, or because a pulse fails to follow the relevant instantaneous eigenstates. A far-detuned dark state can be retained correctly by elimination; using a dark-state protocol is not itself a reason for failure. Sudden preparation of a bare state can additionally excite fast intermediate-state transients absent from the adiabatic solution. The mathematical treatment is in Adiabatic Elimination.

Intermediate-state scattering

If $\Gamma_p$ is the population decay rate of $|p\rangle$, then

$$\Gamma_{\mathrm{sc}}(t) = \Gamma_p |c_p(t)|^2.$$

For a normalized pure state in the retained subspace, Cauchy–Schwarz gives the leading adiabatic upper bound

$$\Gamma_{\mathrm{sc}} \lesssim \Gamma_p \frac{ |\Omega_1|^2+|\Omega_2|^2 }{ 4\Delta^2 }.$$

This bound is on the leading eliminated model; it does not cover a suddenly excited fast component or omitted intermediate levels. The actual population in $|p\rangle$ depends on the instantaneous state, interference between excitation paths, pulse shape, branching ratios, and all resolved intermediate levels. The notation also hides a common $2\pi$ error: if a natural linewidth is quoted as $\Gamma_p/(2\pi)$ in megahertz, the decay rate in the equation is $2\pi$ times that number.

One-photon and Raman alternatives

A direct one-photon transition removes intermediate-state scattering and the associated differential light shift, but may require ultraviolet light or a metastable initial state. A two-photon path provides convenient wavelengths and selection flexibility, but inherits two optical phases and an intermediate state. More elaborate Raman or stimulated-adiabatic protocols can suppress intermediate population, although they add control parameters and do not remove laser phase noise or off-resonant couplings automatically.

The available route must be evaluated for the actual isotope. Selection rules fix the intended coupling only after specifying quantization axis, polarizations, hyperfine basis, and magnetic-field regime. Atomic Selection Rules gives the angular-momentum treatment.

Top: the light-shift-compensated effective Rabi trajectory in the worked example gives $P_r=\sin^2(\pi t/2t_\pi)$ and a changing intermediate population, plotted as $P_p/(2.5\times10^{-3})$. The scattering probability depends on the time integral of $P_p$. Bottom: ballistic thermal averaging at $10\ \mu\mathrm K$ gives Gaussian Ramsey contrast for absorbed $780$ and $480\ \mathrm{nm}$ photons. A small Doppler width relative to a fast Rabi rate can still limit a longer Ramsey interrogation.

Phase and Doppler Geometry

For two absorbed photons, the two-photon coupling carries the spatial phase

$$\Omega_{\mathrm{eff}}(\mathbf r) \propto \exp\!\left( i\mathbf k_{\mathrm{eff}}\cdot\mathbf r \right),$$

where $\mathbf k_{\mathrm{eff}}$ is the signed sum of the participating wavevectors. Define the Doppler shift

$$\delta_D = \mathbf k_{\mathrm{eff}}\cdot\mathbf v.$$

The detuning seen by a moving atom is $\delta_{\mathrm{eff}}-\delta_D$, with the sign fixed by the absorption phase above. For a ladder driven by absorption of both photons,

$$\mathbf k_{\mathrm{eff}} = \mathbf k_1+\mathbf k_2.$$

Counterpropagating beams with unequal wavelengths reduce $|\mathbf k_{\mathrm{eff}}|$ but do not make it vanish. Copropagating beams increase the Doppler sensitivity. The geometry must also be compared with the trap axes: a cold atom can have different velocity widths along different directions, and local addressing beams can imprint site-dependent phases.

For a one-dimensional Maxwell distribution,

$$\sigma_v = \sqrt{\frac{k_{\mathrm B}T}{m}}, \qquad \frac{\sigma_D}{2\pi} = \frac{ |\mathbf k_{\mathrm{eff}}|\sigma_v }{ 2\pi }.$$

The absorption also changes center-of-mass kinetic energy. For a free atom with initial momentum $m\mathbf v$, resonance includes the recoil frequency

$$\omega_{\mathrm{rec}}=\frac{\hbar|\mathbf k_{\mathrm{eff}}|^2}{2m}, \qquad \delta_{\mathrm{atom}}=\delta_{\mathrm{eff}} -\mathbf k_{\mathrm{eff}}\cdot\mathbf v-\omega_{\mathrm{rec}}.$$

A trapped, resolved-sideband spectrum requires the quantized motional model instead of this free-particle resonance condition.

Cooling helps in two distinct ways. It narrows the Doppler distribution and reduces position-dependent sampling of the laser intensity and pair interaction. Neither effect is captured by quoting temperature alone; the relevant quantities are the motional distribution, trap frequencies, beam geometry, and pulse duration.

Worked Excitation and Scattering Estimate

Consider a simplified $^{87}\mathrm{Rb}$ ladder with

$$\begin{aligned} \frac{\Omega_1}{2\pi}&=100\ \mathrm{MHz}, & \frac{\Omega_2}{2\pi}&=20\ \mathrm{MHz}, \\ \frac{\Delta}{2\pi}&=1.0\ \mathrm{GHz}, & \frac{\Gamma_p}{2\pi}&=6.0\ \mathrm{MHz}. \end{aligned}$$

The effective Rabi frequency is

$$\frac{\Omega_{\mathrm{eff}}}{2\pi} \simeq \frac{ (100\ \mathrm{MHz})(20\ \mathrm{MHz}) }{ 2(1000\ \mathrm{MHz}) } = 1.0\ \mathrm{MHz},$$

so an ideal resonant pi pulse lasts

$$t_\pi = \frac{\pi}{\Omega_{\mathrm{eff}}} = 0.50\ \mu\mathrm s.$$

The light-shift scales are

$$\frac{\Delta_g^{\mathrm{AC}}}{2\pi} \simeq 2.5\ \mathrm{MHz}, \qquad \frac{\Delta_r^{\mathrm{AC}}}{2\pi} \simeq 0.10\ \mathrm{MHz}.$$

Their differential value is therefore about $2.4\ \mathrm{MHz}$ in magnitude, larger than $\Omega_{\mathrm{eff}}/(2\pi)$. A laser pair placed at the bare two-photon resonance would be substantially off resonance.

Choose real positive couplings and compensate the light shift by setting $\delta=\Delta_r^{\mathrm{AC}}-\Delta_g^{\mathrm{AC}}$, here $\delta/(2\pi)=-2.4\ \mathrm{MHz}$. Neglect scattering in computing the leading Rabi trajectory and assume the fast intermediate-state transient has been suppressed by turn-on. After removing a common phase,

$$c_g=\cos(\Omega_{\mathrm{eff}}t/2),\qquad c_r=-i\sin(\Omega_{\mathrm{eff}}t/2).$$

The two contributions to $c_p$ are in quadrature, giving

$$P_p(t)\simeq \frac{\Omega_1^2\cos^2(\Omega_{\mathrm{eff}}t/2) +\Omega_2^2\sin^2(\Omega_{\mathrm{eff}}t/2)} {4\Delta^2}.$$

It falls from $2.5\times10^{-3}$ to $1.0\times10^{-4}$ over the ideal pulse; its average is $1.3\times10^{-3}$. Thus the leading mean number of intermediate-state scattering events is

$$N_{\mathrm{sc}}=\Gamma_p\int_0^{t_\pi}P_p(t)\,dt \simeq0.0245.$$

For rare events this is also the first-order scattering probability, about $2.45\%$. Freezing the initial population throughout the pulse would instead give $4.7\%$. Neither estimate replaces the full multilevel open dynamics when turn-on transients, branching, or repeated scattering matter. Small intermediate population alone does not ensure a negligible integrated error.

Now take counterpropagating $780\ \mathrm{nm}$ and $480\ \mathrm{nm}$ beams. For $T=10\ \mu\mathrm K$,

$$\sigma_v \simeq 3.09\times10^{-2}\ \mathrm{m\,s^{-1}},$$

and

$$\frac{\sigma_D}{2\pi} \simeq 24.8\ \mathrm{kHz}.$$

The copropagating geometry gives approximately $104\ \mathrm{kHz}$. Both are below the nominal $1\ \mathrm{MHz}$ Rabi frequency. Nevertheless, a frozen Gaussian velocity ensemble during free flight gives the Ramsey envelope

$$C(t)=\left|\langle e^{i\mathbf k_{\mathrm{eff}}\cdot\mathbf v t}\rangle\right| =e^{-\sigma_D^2t^2/2}.$$

Its $1/e$ times are about $9.1\ \mu\mathrm s$ and $2.2\ \mu\mathrm s$, respectively. These are reversible inhomogeneous-dephasing scales, not radiative lifetimes.

This simplified calculation illustrates the relevant scales. Real $^{87}\mathrm{Rb}$ requires the intermediate hyperfine manifold, actual Clebsch–Gordan coefficients, beam waists, laser spectra, polarization impurities, and branching to states outside the chosen basis.

Electric Fields and Stark Maps

An isolated nondegenerate parity eigenstate with no permanent laboratory-frame dipole has a leading quadratic Stark shift. In a scalar-response model,

$$\Delta E_r = -\frac{1}{2}\alpha_r(0)E^2.$$

Because the Rydberg polarizability can be enormous, a field too small to matter for a ground-state transition can shift the Rydberg resonance by many linewidths. Near-degenerate manifolds require diagonalizing the Stark Hamiltonian; their response can become linear over an experimentally relevant range. A scalar response is also inadequate when tensor polarizability or field-induced sublevel mixing matters; one must specify the field direction and atomic state.

Compensation by spectroscopy

Suppose an electrode voltage $V$ produces field

$$E(V) = E_{\mathrm{stray}}+\beta(V-V_0).$$

In the quadratic regime, the transition frequency follows

$$\nu(V) = \nu_0 -\frac{\alpha_{\mathrm{diff}}}{2h} \left[ E_{\mathrm{stray}}+\beta(V-V_0) \right]^2.$$

Here $\nu_0$ is the frequency at zero field in the scalar model; for $\alpha_{\mathrm{diff}}>0$ the vertex is a maximum, not a minimum. A transverse residual field adds a constant quadratic shift. The vertex identifies compensation only along the electrode’s field direction. Three-dimensional compensation requires independent field directions or a model of the electrode response matrix. Cross-couplings, voltage offsets, and patch charges make this a matrix calibration, not three unrelated knob scans.

An experimental Stark map serves several purposes:

• identifies the intended Rydberg state and nearby avoided crossings;
• measures the residual electric field and electrode response;
• reveals whether a quadratic approximation is valid;
• locates interaction resonances and unwanted pair channels; and
• tracks slow drift from surface charging, adsorbates, or photoelectrons.

Compensation is not permanent. Ultraviolet or blue excitation light can alter surface charge, and moving optics or changing trap power can change the local environment. Repeating the Stark spectrum can reveal such drift.

Magnetic fields matter as well. They define the quantization axis and lift Zeeman degeneracies, but field gradients and calibration errors create site-dependent detunings. Electric and magnetic shifts should be separated by their symmetry under field reversal and by spectroscopy of multiple states.

Population Lifetime and Coherence

Radiative and blackbody processes

For ordinary low-angular-momentum Rydberg series at zero temperature, the radiative lifetime often scales approximately as $(n^*)^3$. At finite temperature, ambient blackbody photons drive transitions to nearby Rydberg states and can contribute photoionization. If these are independent losses from the selected state with negligible repopulation, the exponential lifetime satisfies

$$\frac{1}{T_1} = \Gamma_{\mathrm{rad}} +\Gamma_{\mathrm{BBR}} +\Gamma_{\mathrm{photo}} +\Gamma_{\mathrm{coll}} +\cdots.$$

The terms depend on state, temperature, electromagnetic environment, atomic density, and applied light. A tabulated free-space lifetime is therefore not automatically the in-apparatus lifetime. Blackbody transfer can populate other interacting Rydberg states rather than remove the atom. If those states feed back, population survival requires coupled rate equations, not a single exponential.

Lifetime is not coherence time

The coherence between $|g\rangle$ and $|r\rangle$ can decay faster than the Rydberg population:

$$\frac{1}{T_2} = \frac{1}{2T_1} +\Gamma_\phi.$$

Here $\Gamma_\phi$ collects pure dephasing from laser phase noise, Doppler shifts, electric-field noise, magnetic shifts, intensity fluctuations, interaction inhomogeneity, and unresolved state mixing. The equality assumes stable $|g\rangle$, exponential population loss from $|r\rangle$, and Markovian pure dephasing. Quasistatic noise can be stationary and still give nonexponential Ramsey contrast, as the Doppler example shows; a single constant $\Gamma_\phi$ then does not describe it.

Three measurements answer different questions:

Measurement	Main information	Important caveat
delayed population survival	$T_1$ and transfer to dark channels	detector may not distinguish decay products
Ramsey contrast	inhomogeneous and low-frequency dephasing	includes shot-to-shot detuning drift
spin echo or dynamical decoupling	refocusable noise versus irreversible loss	pulse errors can mimic decoherence

A fitted exponential time and a Gaussian $1/e$ time describe different envelopes even if their numerical values coincide.

Trapping During Rydberg Excitation

For a field varying little over the electronic orbit, the weak-field scalar optical potential uses the real dynamic polarizability and peak complex field amplitude:

$$U_a(\mathbf r) \simeq -\frac{1}{4} \operatorname{Re}\alpha_a(\omega_L) |\mathcal E(\mathbf r)|^2.$$

A Rydberg orbit can be large enough to violate this local-field assumption. In the off-resonant nearly-free-electron approximation its ponderomotive contribution is instead

$$U_{\mathrm{Ry}}(\mathbf R)\simeq \frac{e^2}{4m_e\omega_L^2} \int d^3r\,|\psi_r(\mathbf r)|^2 |\mathcal E(\mathbf R+\mathbf r)|^2+U_{\mathrm{core}}(\mathbf R).$$

Here $\mathbf R$ locates the ionic core and $\mathbf r$ is the electron’s coordinate relative to it. Spatial averaging can change trap curvature; the electronic state must remain well defined for this averaged-potential description. Topcu and Derevianko analyze this intensity-landscape effect.

Ground and Rydberg states generally have very different polarizabilities at a tweezer wavelength. The consequences include

• a differential transition shift;
• different trapping or antitrapping forces;
• position-dependent detuning across the atomic wavepacket;
• motional excitation during state changes; and
• possible photoionization or coupling to nearby Rydberg levels.

Dynamic Polarizability develops the frequency-dependent response theory. There are three common strategies.

Turn the trap off

The tweezer is extinguished during the Rydberg pulse and restored before imaging. This suppresses trap-induced differential shifts, but the atom evolves approximately freely when interaction forces, residual trapping, and gravity are negligible over the pulse:

$$\mathbf r(t) = \mathbf r(0)+\mathbf v(0)t.$$

Long dark windows reduce recapture and change pair separations. Trap switching can also impart parametric excitation if the timing or extinction is imperfect.

Leave the trap on

Leaving the light on can improve recapture if both relevant states are confined; an antitrapped Rydberg state can instead accelerate outward. This requires determining the Rydberg-state potential, differential light shift, and any photoionization. A single transition frequency measured at one trap depth does not establish that the potential is benign.

Engineer state-insensitive trapping

At selected wavelengths, polarizations, or dressed operating points, the differential polarizability can be reduced. Such “magic” conditions are state-specific and can be spoiled by tensor structure, nearby resonances, polarization impurity, or intensity variation. They must be demonstrated by spectroscopy and motional measurements for the states actually used.

Motion changes both drive and interaction

In a thermal wavepacket, the Rabi frequency, detuning, and interaction are random variables:

$$\Omega_i=\Omega(\mathbf r_i), \qquad \Delta_i=\Delta(\mathbf r_i,\mathbf v_i), \qquad V_{ij}=V(\mathbf r_i-\mathbf r_j).$$

Replacing them by mean values can overpredict contrast because coherent averaging is nonlinear. A reliable model samples the measured position and velocity distributions or propagates the motional degrees of freedom explicitly.

Long-Range Interactions as a Platform Resource

The atomic origin of resonant dipole–dipole and off-resonant van der Waals interactions is developed in Rydberg Atoms Basics. At the platform level, the essential question is whether a chosen pair state can be represented by a scalar potential

$$V_{ij} \simeq \frac{C_6}{R_{ij}^6}$$

over the experiment’s range of separation, angle, fields, and detuning.

That approximation can fail because

• several Zeeman pair states are coupled;
• the interaction is anisotropic;
• a Förster defect becomes comparable to a dipole coupling;
• external fields mix pair channels;
• higher multipoles matter at short range;
• atomic motion samples a broad range of $R$; or
• the laser couples to more than one molecular eigenpotential.

The remedy is not to fit an arbitrary effective $C_6$ to every dataset. First calculate or spectroscopically map the relevant pair eigenstates, then state the interval over which a reduced model is valid.

Sensitivity to position

For a van der Waals interaction,

$$\frac{\delta V}{V} \simeq -6\frac{\delta R}{R}.$$

At $R=5.0\ \mu\mathrm m$, a relative-separation uncertainty of $\sigma_R=50\ \mathrm{nm}$ gives

$$\frac{\sigma_V}{|V|} \simeq 6\frac{0.050}{5.0} = 0.060.$$

A seemingly small one-percent position uncertainty therefore produces a six-percent interaction uncertainty. Angular noise can add another first-order contribution when the potential is anisotropic.

Effective Array Hamiltonian

After single-atom and pair calibration, a commonly used rotating-frame model is

$$H ={} \frac{\hbar}{2} \sum_i \left[ \Omega_i e^{i\phi_i} |r_i\rangle\langle g_i| + \Omega_i^* e^{-i\phi_i} |g_i\rangle\langle r_i| \right] - \hbar\sum_i \Delta_i n_i + \sum_{i<j}V_{ij}n_i n_j,$$

Choose $\Omega_i\geq0$ and put the whole drive phase in $\phi_i$. Here

$$n_i = |r_i\rangle\langle r_i|.$$

This Hamiltonian is experimentally meaningful only with a parameter map:

• $\Omega_i$ includes beam intensity, polarization, and atomic position;
• $\phi_i$ includes optical phase and site geometry;
• $\Delta_i$ includes programmed detuning plus Stark, Zeeman, Doppler, and light shifts;
• $V_{ij}$ comes from a stated pair-state model or measurement; and
• omitted loss and dephasing channels have measured bounds.

Static site phases can always be removed by local basis rephasings for this pure density-density model. Time-dependent rephasings add detunings; exchange or hopping terms can make relative phases physically significant. If $\Omega_i=\Omega$ and $\Delta_i=\Delta$, the two internal states map to a spin-$1/2$ model. But “realizes an Ising model” does not mean every term beyond that ideal model is absent. It means the residual terms are calculated, measured, or bounded over the reported operating region.

Lattice models overview develops the many-body interpretation. Exercise 8 makes the interaction-induced longitudinal spin fields explicit.

Excitation and Detection

Coherent excitation

A resonant Rabi trace tests coherent transfer but does not identify every source of contrast loss. For an ideal two-level atom initially in $|g\rangle$,

$$P_r(t) = \frac{\Omega^2}{\Omega^2+\Delta^2} \sin^2\!\left[ \frac{t}{2} \sqrt{\Omega^2+\Delta^2} \right].$$

Observed contrast loss can arise from decay, dephasing, detuning inhomogeneity, Rabi-frequency inhomogeneity, imperfect preparation, or readout. A damped sine fit may summarize the data but generally does not identify which mechanism dominates.

Loss-based detection

In many alkali tweezer experiments, ground-state atoms are recaptured and imaged, whereas Rydberg excitation tends to produce absence. Absence is not a projector onto $|r\rangle$. It can also result from ordinary background loss, heating, failed recapture, off-resonant transfer, ionization, or an imaging error.

Let $a$ denote the measured outcome “absent.” A binary confusion model is

$$\begin{pmatrix} P(a) \\ P(\bar a) \end{pmatrix} = \begin{pmatrix} P(a|g) & P(a|r) \\ P(\bar a|g) & P(\bar a|r) \end{pmatrix} \begin{pmatrix} P_g \\ P_r \end{pmatrix}.$$

This binary equation assumes $P_g+P_r=1$ at readout and a measurement diagonal in that basis. Initial vacancies and leakage require extra columns or an explicitly conditional sample; two outcomes cannot identify several unknown populations. Even in the two-state model, inversion requires $P(a|r)\ne P(a|g)$. The probabilities must be calibrated under conditions close to those of the science sequence. Inverting it without propagating uncertainty can turn a small calibration error into a large population error.

More state-selective methods

Destructive field ionization directly detects charged fragments and can resolve arrival-time windows associated with different states. Rydberg de-excitation followed by recapture maps $|r\rangle$ back into a trapped state. Two-electron atoms can permit autoionization of the Rydberg electron, enabling state-selective detection while retaining additional internal information. Shelving and fluorescence protocols provide other mappings.

Worked Readout Correction

Suppose a ground-state atom survives and is detected with probability

$$s_g=0.98,$$

so $P(a|g)=0.02$. Suppose a prepared Rydberg atom is absent with probability

$$\ell_r=P(a|r)=0.95.$$

If the measured absent fraction is $P(a)=0.40$, then

$$P(a) = \ell_r P_r + (1-s_g)(1-P_r).$$

Solving,

$$P_r = \frac{ P(a)-(1-s_g) }{ \ell_r-(1-s_g) } = \frac{0.40-0.02}{0.95-0.02} \simeq 0.409.$$

Calling the raw absent fraction “$40\%$ Rydberg population” would be close in this example, but only accidentally. The correction grows when ordinary loss increases or Rydberg detection becomes less efficient. Correlated loss also invalidates a product of independent single-site confusion matrices.

Exercises

Exercise 1: Eliminate the intermediate state

Use the three-level Hamiltonian on this page. Assume real $\Delta$ and $|\Delta|\gg|\Omega_1|,|\Omega_2|,|\delta|$. Derive the leading effective Hamiltonian in the $(|g\rangle,|r\rangle)$ subspace. Identify the two-photon Rabi frequency and the two diagonal light shifts.

Solution

The amplitude equations are

$$\begin{aligned} i\dot c_g &= \frac{\Omega_1^*}{2}c_p, \\ i\dot c_p &= \frac{\Omega_1}{2}c_g -\Delta c_p +\frac{\Omega_2^*}{2}c_r, \\ i\dot c_r &= \frac{\Omega_2}{2}c_p -\delta c_r. \end{aligned}$$

At leading adiabatic order, set $\dot c_p\simeq0$:

$$c_p \simeq \frac{ \Omega_1 c_g+\Omega_2^*c_r }{ 2\Delta }.$$

Substitution gives

$$i\frac{d}{dt} \begin{pmatrix} c_g\\c_r \end{pmatrix} = \begin{pmatrix} \dfrac{|\Omega_1|^2}{4\Delta} & \dfrac{\Omega_1^*\Omega_2^*}{4\Delta} \\ \dfrac{\Omega_1\Omega_2}{4\Delta} & -\delta+\dfrac{|\Omega_2|^2}{4\Delta} \end{pmatrix} \begin{pmatrix} c_g\\c_r \end{pmatrix}.$$

Writing the off-diagonal element as $\Omega_{\mathrm{eff}}/2$ gives

$$\Omega_{\mathrm{eff}} = \frac{\Omega_1\Omega_2}{2\Delta},$$

with phases assigned according to the definitions of $\Omega_1$ and $\Omega_2$. The diagonal shifts are

$$\Delta_g^{\mathrm{AC}} = \frac{|\Omega_1|^2}{4\Delta}, \qquad \Delta_r^{\mathrm{AC}} = \frac{|\Omega_2|^2}{4\Delta}.$$

Changing the detuning convention changes their displayed signs but not the need to retain the differential shift.

Exercise 2: Balance the excitation legs

For the compensated ideal Rabi trajectory above, keep $\Delta>0$ and $\Omega_{\mathrm{eff}}=\Omega_1\Omega_2/(2\Delta)$ fixed while varying the positive ratio $r=\Omega_1/\Omega_2$. Minimize the integrated intermediate-state scattering. What assumptions make this optimization useful?

Solution

The mean intermediate population during a pi pulse is $(\Omega_1^2+\Omega_2^2)/(8\Delta^2)$ and $t_\pi=2\pi\Delta/(\Omega_1\Omega_2)$. Consequently

$$N_{\mathrm{sc}}=\frac{\pi\Gamma_p}{4\Delta} \left(r+\frac1r\right).$$

The inequality $r+1/r\geq2$ gives equal single-photon Rabi amplitudes as the optimum. The worked ratio $r=5$ gives $N_{\mathrm{sc}}\simeq0.0245$; equal legs at the same effective rate give $N_{\mathrm{sc}}\simeq0.00942$, or $0.94\%$. Their Rabi amplitudes are each $\sqrt{100\times20}\ \mathrm{MHz}\simeq44.7\ \mathrm{MHz}$ after division by $2\pi$.

The result neglects turn-on transients, other intermediate states, scattering backaction, and unequal laser-power costs. It also requires reapplying the differential-shift compensation as the ratio changes. Equal optical powers need not give equal Rabi amplitudes because the dipole matrix elements differ.

Exercise 3: A trapped Ramsey revival

In a classical thermal harmonic trap of frequency $\omega_t$, suppose both internal states experience the same potential and recoil is negligible. Ideal short Ramsey pulses sample the phase difference $k_{\mathrm{eff}}[x(t)-x(0)]$. Find its contrast and recover the ballistic short-time limit.

Solution

The trajectory is $x(t)=x_0\cos\omega_t t+(v_0/\omega_t)\sin\omega_t t$. Thermal equilibrium has independent Gaussian $x_0,v_0$, with variances $k_BT/(m\omega_t^2)$ and $k_BT/m$. Therefore

$$\operatorname{Var}[x(t)-x_0] =\frac{2k_BT}{m\omega_t^2}(1-\cos\omega_t t),$$

and Gaussian averaging gives

$$C(t)=\exp\!\left[ -\frac{k_{\mathrm{eff}}^2k_BT}{m\omega_t^2} (1-\cos\omega_t t)\right].$$

At short times this is $\exp[-k_{\mathrm{eff}}^2k_BTt^2/(2m)]$, the ballistic result. The ideal model revives at $2\pi/\omega_t$: initial velocity uncertainty has not caused irreversible decoherence. Different internal-state potentials, anharmonicity, decay, recoil, and pulse errors modify or suppress the revival.

Exercise 4: Locate a Stark-compensation point

A Rydberg transition has a quadratic shift

$$\Delta\nu(V) = -A(E_s+\beta V)^2,$$

with $A>0$. Show that the vertex occurs at $V_c=-E_s/\beta$. If spectra along one electrode direction find $V_c=0.36\ \mathrm V$ and $\beta=0.25\ \mathrm{V\,cm^{-1}\,V^{-1}}$, what component of stray electric field is compensated? Why does this not establish full three-dimensional compensation?

Solution

Differentiate:

$$\frac{d\Delta\nu}{dV} = -2A\beta(E_s+\beta V).$$

The derivative vanishes at

$$V_c = -\frac{E_s}{\beta}.$$

Therefore

$$E_s = -\beta V_c = -(0.25)(0.36)\ \mathrm{V\,cm^{-1}} = -0.090\ \mathrm{V\,cm^{-1}}.$$

The electrode compensates only the projection of the stray field along its response direction. Transverse components remain, and real electrodes can produce nonorthogonal fields. A multidimensional voltage scan or calibrated electrode-response matrix is required.

Exercise 5: Distinguish T1 from T2

A Rydberg state has independent characteristic times

$$T_{\mathrm{rad}}=150\ \mu\mathrm s, \qquad T_{\mathrm{BBR}}=300\ \mu\mathrm s, \qquad T_{\mathrm{photo}}=1000\ \mu\mathrm s.$$

Calculate $T_1$. If the pure-dephasing time is $T_\phi=80\ \mu\mathrm s$, calculate $T_2$ using $1/T_2=1/(2T_1)+1/T_\phi$.

Solution

Independent population-loss rates add:

$$\frac{1}{T_1} = \frac{1}{150} +\frac{1}{300} +\frac{1}{1000} \quad \mu\mathrm s^{-1}.$$

Therefore

$$T_1 \simeq 90.9\ \mu\mathrm s.$$

The coherence time is

$$T_2 = \left[ \frac{1}{2(90.9)} +\frac{1}{80} \right]^{-1} \mu\mathrm s \simeq 55.6\ \mu\mathrm s.$$

The result is shorter than both $2T_1$ and $T_\phi$ because both mechanisms degrade coherence.

Exercise 6: Convert position noise into interaction noise

Two atoms interact through $V=C_6/R^6$. Linearize the interaction around $R_0$ and find the fractional rms interaction noise for $R_0=5.0\ \mu\mathrm m$ and $\sigma_R=50\ \mathrm{nm}$. If $|V(R_0)|/h=20\ \mathrm{MHz}$, estimate the rms frequency variation.

Solution

Differentiating,

$$\frac{dV}{dR} = -\frac{6C_6}{R^7} = -\frac{6V}{R}.$$

For small fluctuations,

$$\frac{\sigma_V}{|V|} \simeq 6\frac{\sigma_R}{R_0} = 6\frac{0.050}{5.0} = 0.060.$$

The interaction-frequency variation is therefore

$$\frac{\sigma_V}{h} \simeq 0.060(20\ \mathrm{MHz}) = 1.2\ \mathrm{MHz}.$$

This linear rms estimate needs a narrow radial distribution with finite variance, not specifically a Gaussian one; it ignores angular dependence.

Exercise 7: Correct loss-based detection

The probability that a ground-state atom appears absent is $0.02$, and the probability that a Rydberg atom appears absent is $0.95$. An experiment measures an absent fraction of $0.40$.

1. Infer the Rydberg population.
2. Repeat for an absent fraction of $0.04$.
3. Explain why the second inference is especially sensitive to calibration uncertainty.

Solution

The readout equation is

$$p_a = 0.02(1-p_r)+0.95p_r = 0.02+0.93p_r.$$

For $p_a=0.40$,

$$p_r = \frac{0.40-0.02}{0.93} \simeq 0.409.$$

For $p_a=0.04$,

$$p_r = \frac{0.04-0.02}{0.93} \simeq 0.0215.$$

In the second case, the signal above the ground-state loss floor is only $0.02$. A small error in that floor is comparable to the inferred Rydberg population. Confidence intervals must therefore propagate calibration counts and any drift, rather than treating $0.02$ and $0.95$ as exact constants.

Exercise 8: The longitudinal field in the Ising model

Set all static drive phases to zero and define $\sigma_i^z=|r_i\rangle\langle r_i|-|g_i\rangle\langle g_i|$. Rewrite the array Hamiltonian in Pauli operators. Does uniform laser detuning imply a uniform longitudinal spin field in an open finite array?

Solution

Since $n_i=(1+\sigma_i^z)/2$,

$$H=\frac{\hbar}{2}\sum_i\Omega_i\sigma_i^x +\sum_i\left(-\frac{\hbar\Delta_i}{2} +\frac14\sum_{j\ne i}V_{ij}\right)\sigma_i^z +\frac14\sum_{i<j}V_{ij}\sigma_i^z\sigma_j^z +\text{constant}.$$

The interaction contributes to both the pair term and the longitudinal field. Edge atoms have a different interaction sum even if $\Delta_i$ is uniform. Canceling the longitudinal field requires $\Delta_i=(2\hbar)^{-1}\sum_{j\ne i}V_{ij}$, generally a site-dependent detuning in a finite array. Omitting these single-spin terms changes the physical resonance condition.

References

• C. S. Adams, J. D. Pritchard, and J. P. Shaffer, “Rydberg atom quantum technologies,” Journal of Physics B 53, 012002 (2020), doi:10.1088/1361-6455/ab52ef.
• L. Béguin, A. Vernier, R. Chicireanu, T. Lahaye, and A. Browaeys, “Direct measurement of the van der Waals interaction between two Rydberg atoms,” Physical Review Letters 110, 263201 (2013), doi:10.1103/PhysRevLett.110.263201.
• I. I. Beterov, I. I. Ryabtsev, D. B. Tretyakov, and V. M. Entin, “Quasiclassical calculations of blackbody-radiation-induced depopulation rates and effective lifetimes of Rydberg $nS$, $nP$, and $nD$ alkali-metal atoms with $n\leq80$,” Physical Review A 79, 052504 (2009); erratum 80, 059902 (2009), doi:10.1103/PhysRevA.79.052504.
• A. Browaeys and T. Lahaye, “Many-body physics with individually controlled Rydberg atoms,” Nature Physics 16, 132–142 (2020), doi:10.1038/s41567-019-0733-z.
• T. F. Gallagher, “Rydberg atoms,” Reports on Progress in Physics 51, 143–188 (1988), doi:10.1088/0034-4885/51/2/001.
• T. F. Gallagher, Rydberg Atoms (Cambridge University Press, 1994), doi:10.1017/CBO9780511524530.
• D. Jaksch, J. I. Cirac, P. Zoller, S. L. Rolston, R. Côté, and M. D. Lukin, “Fast quantum gates for neutral atoms,” Physical Review Letters 85, 2208–2211 (2000), doi:10.1103/PhysRevLett.85.2208.
• M. Morgado and S. Whitlock, “Quantum simulation and computing with Rydberg-interacting qubits,” AVS Quantum Science 3, 023501 (2021), doi:10.1116/5.0036562.
• S. Ravets, H. Labuhn, D. Barredo, L. Béguin, T. Lahaye, and A. Browaeys, “Coherent dipole–dipole coupling between two single Rydberg atoms at an electrically-tuned Förster resonance,” Nature Physics 10, 914–917 (2014), doi:10.1038/nphys3119.
• M. Saffman, “Quantum computing with atomic qubits and Rydberg interactions: progress and challenges,” Journal of Physics B 49, 202001 (2016), doi:10.1088/0953-4075/49/20/202001.
• M. Saffman, T. G. Walker, and K. Mølmer, “Quantum information with Rydberg atoms,” Reviews of Modern Physics 82, 2313–2363 (2010), doi:10.1103/RevModPhys.82.2313.
• T. Topcu and A. Derevianko, “Intensity landscape and the possibility of magic trapping of alkali-metal Rydberg atoms in infrared optical lattices,” Physical Review A 88, 043407 (2013), doi:10.1103/PhysRevA.88.043407.