Stark Effect on Rubidium Rydberg Atoms

• Stark effects in rubidium Rydberg atoms are defined by electric field-induced energy shifts that transition from quadratic behavior in weak fields to complex nonperturbative dynamics at high fields.
• Numerical diagonalization of the Stark Hamiltonian, incorporating quantum defects, fine structure, and hyperfine splittings, enables precise mapping of energy levels and ionization thresholds.
• These effects underpin practical applications such as high-resolution field sensing, quantum gate implementations, and the investigation of electron dynamics under extreme conditions.

Rubidium Rydberg atoms—alkali atoms excited to high principal quantum number $n$—exhibit pronounced Stark effects across a variety of field regimes due to their large size and strong polarizability. The Stark effect encompasses both the behavior of energy levels under static (dc) and oscillating (ac) electric fields, with a crossover from perturbative quadratic shifts to complex nonperturbative dynamics including manifold mixing, avoided crossings, ionization, and resonance phenomena. In rubidium, fine structure, quantum defects, and hyperfine structure further enrich the spectra relative to the hydrogenic case. Stark effects on Rydberg atoms underpin high-resolution field sensing, quantum gate protocols, and studies of electron dynamics under extreme conditions.

1. Theoretical Formulation of the Stark Hamiltonian

In the presence of a static electric field $E$ oriented along $z$, the single-electron Hamiltonian for rubidium Rydberg atoms is

$H = H_0 + eEz$

where $H_0$ is the field-free Hamiltonian, including quantum defects $\delta_{\ell j}$ via

$E_{n\ell j} = -\frac{1}{2 (n - \delta_{\ell j})^2}$

for principal quantum number $n$, azimuthal quantum number $\ell$, and total angular momentum $j$ ($\hbar = 1$, atomic units). For non-degenerate states or weak fields, the energy shift is given by

$$\Delta E(E) = \Delta E^{(1)} + \Delta E^{(2)} + \cdots$$

with

$$\Delta E^{(1)} = eE \langle n\ell m | z | n\ell m \rangle = 0\;\;\text{for S states}$$

and the quadratic Stark effect

$$\Delta E^{(2)} = -\frac{1}{2} \alpha(n, \ell, m) E^2$$

where $\alpha$ is the polarizability. For hydrogen-like states, $\alpha \sim n^7$, and near-hydrogenic high-$\ell$ states, linear Stark shifts can appear due to manifold mixing.

At higher fields, nonperturbative diagonalization of $H$ is required. The classical ionization threshold is set by the saddle point in the Coulomb-plus-field potential:

$$V_C(r,z) = -\frac{1}{r} - Ez,\quad \text{with}\quad E_\mathrm{ion} = -2\sqrt{E}$$

defining the field above which levels enter the continuum (Grimmel et al., 2017, Grimmel et al., 2015).

2. Stark Spectra: Numerical Methods and Experimental Probes

High-precision mapping of Stark-shifted energy levels in rubidium Rydberg manifolds leverages both theoretical diagonalization and state-of-the-art experimental techniques.

2.1 Numerical Stark Map Calculation

Effective modeling requires constructing the Stark Hamiltonian in the basis $|n\ell j m_j\rangle$ with quantum defects and fine structure included. Matrix elements for the $z$ operator are determined via radial and angular factors, often following the technique of Zimmerman et al. (PRA 20, 2251), and the Hamiltonian is diagonalized separately for each $m_j$. Basis truncation at several thousand states is required for high $n$ (e.g., $n=43$ requires $\sim4000$ states; $n=70$ requires $\sim 10^4$ states) (Grimmel et al., 2017, Grimmel et al., 2015).

2.2 Complex Absorbing Potential for Ionizing Levels

To treat ionization, a non-Hermitian term $-i\eta W(r; E)$ is added (Grimmel et al., 2017), starting at the saddle point radius $r_c(E)=1/\sqrt{E}$ and scaling as $(r - r_c)^6$ for $r>r_c$. Diagonalization of $H_\mathrm{CAP}$ yields complex eigenvalues $E_c=E_r - i\Gamma/2$, with $\Gamma$ the ionization rate.

2.3 Stark Spectroscopy Measurements

Experimental determination of Stark maps typically uses electromagnetically induced transparency (EIT) spectroscopy in cold atom traps or vapor cells. Key implementations include:

• Two-photon ladder excitation (e.g., $$5S_{1/2}\rightarrow5P_{3/2} \rightarrow nS_{1/2}/nD_{5/2}$$) probed via transmission of a weak probe laser and scanned coupling laser.
• Direct measurement of induced Stark shifts by EIT in the presence of controlled external fields up to $>500$ V/cm (Grimmel et al., 2015).
• Ionization detection using extracted ion currents in strong-field regimes (Grimmel et al., 2017).

Agreement between measurement and calculation is routinely within a few MHz for $n$ up to 70 and fields up to (and exceeding) the classical ionization threshold.

3. Field Regimes and Stark Shift Behavior

3.1 Weak Field Regime: Quadratic Stark Effect

For low fields, the energy shift is strictly quadratic, $\Delta E=-\frac{1}{2}\alpha_n F^2$, with experimental polarizabilities (in MHz/$(\mathrm{V}/\mathrm{cm})^2$) verified to better than 1% precision (e.g., $\alpha_{20s}=0.0720(8)$ for $n=20$) (Tauschinsky et al., 2013). The scaling $\alpha_n\sim n^7$ holds well for S and D Rydberg states.

3.2 Nonperturbative and Ionization Regime

At higher fields, quadratic approximation fails due to manifold mixing and increasing coupling to the continuum. Stark maps display fan-like avoided crossings, level mixing, and ultimately line broadening and disappearance above the ionization field (e.g., $F_\mathrm{ion}=127$ V/cm for $n=43$ and $16.1$ V/cm for $n=70$) (Grimmel et al., 2017, Grimmel et al., 2015). Notably, long-lived “trapped” resonances can persist well above threshold ($\Gamma\lesssim25$ MHz).

3.3 High-$\ell$ States and Linear Stark Effect

Hydrogen-like high-$\ell$ states ($\ell \geq 4$) display first-order (linear) Stark effects due to near-degeneracy, with maximal dipole moments $\mu_\mathrm{max}\sim1.5 n^2 e a_0$ (e.g., for $n=57$, $|\mu| = (4.14\pm0.10)\times10^{-26}$ C·m, corresponding to $620\pm15\,\mathrm{MHz}/(\mathrm{V}/\mathrm{cm})$ (Duspayev et al., 2023)).

4. AC Stark Effect and Light Shifts

Interaction with oscillating (laser or microwave) fields induces energy shifts (“AC Stark shifts” or “light shifts”) in Rydberg levels. For a Rydberg state $|i\rangle$ under an off-resonant field $\mathbf{E}(t)=\mathbf{E}_0\cos(\omega t)$:

$$\Delta E_i = \sum_{j\neq i} \frac{|\langle i|\boldsymbol{\mu}\cdot\mathbf{E}_0|j\rangle|^2}{\hbar \Delta_{ij}}$$

where $\Delta_{ij} = \omega-\omega_{ij}$. For a two-level system, $$\Delta E_i \approx \hbar \Omega_{ij}^2/(4\Delta_{ij})$$, with Rabi frequency $\Omega_{ij}$ and $\langle nS|\mu|nP\rangle \sim n^2 e a_0$. The AC Stark shift thus scales steeply as $n^7 E_0^2$ (Brekke et al., 2023). Measurements in an optical dipole trap confirm that the shift is dominated by the free-electron ponderomotive potential, with corrections below 10% (Markert et al., 2010). At special “magic wavelengths” (e.g., $\lambda_\mathrm{magic}=1063.529(4)$ nm for $5s\!-\!18s$), the differential AC Stark shift vanishes (Goldschmidt et al., 2015).

5. Motional and Microfield Stark Effects

Moving atoms in a magnetic field experience a motional Stark effect via the Lorentz transformation: an atom at velocity $\mathbf{v}$ in $B$ sees an electric field $\mathbf{E}_L = \mathbf{v}\times\mathbf{B}$. The resulting quadratic energy shift is $\Delta E = -\frac{1}{2}\alpha_n (vB)^2$, with the familiar $n^7$ scaling (Kaiser et al., 2017). Precision studies reveal shifts on the order of 10 MHz for $n=100$ at typical thermal velocities ($v\sim400$ m/s) and 100 G fields.

In cold-atom ion sources, local Coulomb interactions generate stochastic microfields, leading to broadened and asymmetric Stark features. For zero applied field, the distribution converges to the Holtsmark law, while applied fields narrow the distribution and permit directional extraction of ion streams (Duspayev et al., 2023).

6. Stark-Tuning of Förster Resonances

The dc Stark effect enables tuning of Förster-type (resonant dipole-dipole) energy defects in few-atom Rydberg ensembles. The energy mismatch between initial and final product states,

$\Delta(E) = \frac{2E_2(E)-E_1(E)-E_3(E)}{\hbar}$

is tuned to zero by adjusting $E$, maximizing the dipole-dipole interaction (e.g., $E_\mathrm{res} = 1.79$ V/cm tunes the Rb $$37P_{3/2} + 37P_{3/2} \leftrightarrow 37S_{1/2} + 38S_{1/2}$$ channel) (Ryabtsev et al., 2010). Monte Carlo simulations show that such tuning enables coherent Rabi-like populations oscillations and high-fidelity blockade operation, contingent on spatial localization ($\lesssim 1\,\mu$m) and narrow laser linewidth ($\lesssim1$ MHz).

7. Applications and Precision Field Sensing

The exceptional field sensitivity of Rydberg Stark shifts underpins a range of field-imaging and sensing applications:

• Electric field mapping via EIT in the presence of externally applied or stray fields, resolving fields at the 0.1 V/cm scale (Behary et al., 5 Jan 2026).
• Real-time, non-invasive two-dimensional electron beam profiling using DC Stark shifts of Rydberg manifolds, exploiting $n^7$ polarizability scaling (Behary et al., 5 Jan 2026).
• Calibration of local microfield distributions in cold-atom ion sources, with direct applications in focused-ion-beam diagnostics and cold plasma studies (Duspayev et al., 2023).
• Precise knowledge of Stark (and hyperfine) shifts is critical for implementation of Rydberg-mediated quantum gates, as light and field shifts set protocol fidelity limits (Tauschinsky et al., 2013).

Table: Representative Stark Shift Parameters for Rubidium Rydberg States

State	Shift Type	Parameter, Value (Uncertainty)
$20s$	Quadratic	$\alpha_{20s} = 0.0720(8)$ MHz/$(\mathrm{V}/\mathrm{cm})^2$ (Tauschinsky et al., 2013)
$57F_{5/2}$	Linear	$|\mu| = (4.14\pm0.10)\times10^{-26}$ C·m; $620\pm15$ MHz/$(\mathrm{V}/\mathrm{cm})$ (Duspayev et al., 2023)
$60P_{1/2}$	Quadratic	$\alpha = 1122\pm10$ MHz/$(\mathrm{V}/\mathrm{cm})^2$ (Duspayev et al., 2023)

These benchmarks anchor quantitative Stark shift analyses, facilitate cross-comparison between experiment and theory, and set performance limitations for Rydberg-based measurement and control methodologies.

---

Stark effects in rubidium Rydberg atoms combine strong, tunable atomic response with rich nonperturbative structure, making them foundational to quantum measurement, field mapping, and strongly correlated few- and many-body quantum phenomena.