Ultracold Rydberg molecules.

Ultracold molecules formed from association of a single Rydberg atom with surrounding atoms or molecules and those from double Rydberg excitations are discussed in this review. Ultralong-range Rydberg molecules possess a novel molecular bond resulting from scattering of the Rydberg electron from the perturber atoms or molecules. The strong interactions between Rydberg atoms in ultracold gases may lead to formation of macroscopic Rydberg macrodimers. The exquisite control over the properties of the Rydberg electron means that interesting and unusual few-body and quantum many-body features can be realized in such systems.

Introduction

The spectroscopy of excited atoms was instrumental in the development of early atomic physics as well as the evolution of quantum physics. The Balmer series in hydrogen[1], and the spectra of other atoms, notably the alkali-metal atoms by Rydberg[2], provided some of the most important clues to understanding nature on a microscopic scale[3].

The modern picture of the exotic molecules described in this brief review emerged from the investigation of spectral line broadening of highly excited atoms. What happens to an atom with a weakly bound electron in the presence of surrounding atoms or molecules? The answer to this question leads to an intriguing line of inquiry that culminates in our current understanding of ultracold Rydberg molecules, including ultralong-range Rydberg molecules and Rydberg macrodimers, where, respectively, single- and double-Rydberg excitations are involved in their formation.

The review is organized to first provide a historical footing for much of what will follow. This perspective is additionally relevant for a proper accounting of how two of the most important topics in low-energy physics, the scattering length and the pseudopotential, emerged. We then describe the ultralong-range Rydberg molecule and Rydberg macrodimer Hamiltonians. Finally, we report on the spectroscopic properties and laboratory realizations of both types of Rydberg molecules. The few-body and many-body implications of these molecules in atomic quantum gases are then discussed.

Formation of ultracold Rydberg and macrodimer molecules

Collisional interactions between a Rydberg atom and a ground-state atom

The first measurements of the collisional properties of Rydberg atoms were made by Amaldi and Segre[4]. In their ground-breaking experiment, they observed that the shift of atomic Rydberg lines depended on the type of gas in which the sodium or potassium vapors they studied were immersed. They discovered that in different gases (H2, N2, He, and Ar), the alkali-metal Rydberg lines sometimes shifted to the red, and at other times shifted to the blue of the unperturbed atomic lines observed at relatively low gas densities. It was left to Fermi to develop the first theory of collisional shift and broadening of Rydberg lines. In the process of explaining the experimental observations of Amaldi and Segre, Fermi introduced two important concepts to low-energy physics: the scattering length and the pseudopotential[5].

To appreciate Fermi’s insight, consider how a weakly bound electron in a gas interacts with nearby atoms or molecules. If the neighboring species which interact with the Rydberg electron do not possess a permanent electric dipole moment (PEDM), then the interaction of the electron with the surrounding ground-state atoms is in the first-order, the polarization potential. An immediate expectation is that all spectral lines shift to the red, see Box 1.

Fermi recognized that scattering mattered: In the low-energy scattering of the Rydberg electron from perturber atoms, the wave function will be slowly varying, except near the perturbers. He, therefore, invented another wave function as the mean value of ψ around the perturbers, whose extent was assumed to be small compared with the Rydberg electron de Broglie wavelength, but still large enough to contain a sizable number of perturber atoms. The wave function, regular at the origin, can be written u(r) = rψ(r) = (r − a(0)), with a(0) = − the s-wave electron-scattering length, obtained in the zero-momentum limit of the s-wave phase shift, δ(k). In the above, r is the distance from the Rydberg electron to the perturber atoms.

With this approach, the spectral energy shift in the long-wavelength limit becomes, e.g., Eq. (17) in ref. [5] in units of cm−1with me the mass of the electron. The total energy shift is Δ + , wherein units of cm−1, is the shift due to dielectric screening of the Rydberg electron by the gas atoms at density ρ and dielectric constant , see Box 1. In the above equations, e is the elementary charge on the electron, ħ is the reduced Planck constant, and c is the speed of light. In atomic units, e = ħ = 1. The crucial point that unveils itself is that when a(0) > 0 there can be a spectral line shift to the blue if Δ > , but when a(0) < 0 the spectral line shift is always to the red.

To understand the mechanism for the formation of ultralong-range Rydberg molecules, one should consider the interaction that leads to Eq. (1). This interaction between the Rydberg electron at position r and the perturber atom at a distance R = Rz from the Rydberg ion, is the Fermi pseudopotential, henceforth in atomic unitsHere we have extended the concept of zero-energy scattering length to include a momentum dependence, k. It is now possible, without resorting to solving the full molecular Schrödinger equation, to obtain the energy shift due to V(r − R). By employing the Rydberg atom wave function at energy E, the total energy[6-8],In this equation, n and l, are respectively, the principal and orbital angular momentum quantum numbers of the Rydberg electron at energy , while m is the projection quantum number of l onto the internuclear axis. Equation (4) includes the energy shift of the Rydberg electron due to the presence of the perturber atom. We note here that any contributions due to higher-order partial-wave interactions, and relativistic spin corrections are omitted, and discussed later in this review.

A comparison of Eqs. (1) and (4) offers a reinterpretation of density. Whereas the dielectric broadening and shift of the Rydberg lines in Eq. (1) depend on the density of perturbing gas particles, the energy shift in Eq. (4) is proportional to the quantum mechanical electron probability density, i.e., the modulus of wave function squared at the position of the perturber atom.

The form of Eq. (4) suggests that the Born-Oppenheimer (BO) potential energy curves, UBO(R), will be oscillatory functions of R, since the Rydberg wave function oscillates, see Fig. 1. When a(0) is negative, the potentials will possess localized wells at extremely large separations, which scale with n2, and can support bound vibrational levels. Such Rydberg molecules were predicted by Greene et al.[9]. The properties of these bound states can be predicted from the near threshold continuum scattering. The undulations in the BO potential energy curves have been confirmed with detailed quantum chemical calculations, e.g., in ref.[10].

Fig. 1

Ultracold Rydberg molecules and macrodimers. Schematics of Rydberg molecular electronic density and potential energy curves; a depicts an ultralong-range Rydberg molecule while b shows a Rydberg macrodimer. The ultralong-range Rydberg molecule consists of a ground state atom or a polar molecule embedded in the electronic orbital of the Rydberg atom. The potential energy curves exhibit oscillations. A Rydberg macrodimer consists of two Rydberg atoms interacting via macroscopic electrostatic multipolar interactions. The Rydberg electron orbitals (~ 0.1 μm) in this schematic are not shown in proportion to the separation distances

Box 1 Lines redshift in a dielectric medium

The interaction between a Rydberg electron and a ground state atom is, in the first-order perturbation theory, a polarization potential W(q) = −, where α is the polarizability of the ground state particles, and q is the electron position from the ith perturber core. All spectral lines are hence expected to be redshifted from the unperturbed levels, as the polarization potential for a perturber in a non-degenerate ground state is always attractive. The shift of the Rydberg lines due to dielectric screening is obtained in cm−1 as in Eq. (23) in ref. [5],where is the dielectric constant, normalized to vacuum permittivity and ρ is the gas density. Fermi determined that this simple picture was flawed because the Rydberg electron-perturber interaction is fundamentally a scattering process.

Interactions between two Rydberg atoms

When the Rydberg atom density is sufficiently high, Rydberg atoms can interact electrostatically at large distances via van der Waals (vdW)-type forces[11-13]. In a typical ultracold atomic gas, Rydberg–Rydberg interactions can still be substantial at internuclear separations exceeding 10 μm[14-16]. The long-range interactions between Rydberg atoms induce Rydberg blockade[17], due to the large multipolar transition moments between Rydberg states. These interactions can be expressed in an expansion in powers of R−1,where exist in closed form in refs. [15,18,19].

In the above expression, (r, θ, ϕ) are the electronic coordinates of the Rydberg electron on each atom, and are spherical harmonics, with l < the smaller of l1 or l2. Note that l1 and l2 are the moments of the charge distribution and not the orbital angular momenta of the electrons. For example, the set l1 = l2 = 1 gives the dipole–dipole interaction. To obtain the BO potential energy curves, the Hamiltonian, consisting of the individual atomic Hamiltonians and the Rydberg–Rydberg interaction is diagonalized[15]. The full Hamiltonian may also include interactions with external electric and magnetic fields, which are present in experiments. The field dressed potential energy surfaces are dependent on the angle between the internuclear axis and a lab-frame quantization axis[20-22] and R. The lab-frame axis can be determined by an external field or laser polarization.

There exist avoided crossings in the doubly-excited Rydberg BO potential curves due to the large density of states. Each term in the multipole expansion is R-dependent and can be either repulsive or attractive. The R-dependence of the multipolar terms in the expansion in Eq. (5) hybridizes the atomic states as the Rydberg atoms approach, resulting in complex couplings between different molecular states. The hybridization occurs due to the electric fields produced by the Rydberg electron from one atom acting on the other. The interplay between different terms in the multipolar series can conspire to produce attractive wells at large R, in which Rydberg macrodimers may form, Fig. 1. The name macrodimers is apt because of the enormous internuclear separations, >1 μm[19,23-27]. The potential wells can be controlled with small electric fields which shift the Rydberg energy levels and induce avoided crossings of the potential energy surfaces[19,24].

Experiments on Rydberg macrodimers are strong tests of Rydberg atom interactions and serve as a direct test for the accuracy of calculations of the doubly-excited Rydberg BO potentials. These Rydberg molecules can regulate dephasing and decay in thermal vapors which affects applications such as electric field sensing[28,29], probing correlations in quantum gases, and generating gauge potentials[16,23,30,31].

Realization of ultralong-range Rydberg molecules

The first realization of ultralong-range Rydberg molecules was made in an ultracold rubidium gas, in which a fraction of atoms were excited into Rb(ns) states[32]. The two-color excitation process photoassociated 3Σ Rydberg molecules. The electronic charge density for the Rb(5s − ns)3Σ Rydberg molecule should be spherical as was anticipated in ref. [9]. The BO potential in this electronic state contains several vibrational bound states, which were observed as peaks to the red of the atomic Rydberg Rb(35s-37s) lines. The molecular lines were measured by field ionizing the Rydberg molecules. Subsequent measurements observed ultralong-range Rydberg molecules in other species and in different initial atomic angular momentum states[33-40].

A more exotic class of homonuclear Rydberg molecules forms when excited electronic states with nearly degenerate high-angular momentum components mix. This hybridization of opposite parity electronic states results from the electric field of the Rydberg electron in scattering from the perturber atom[9,41]. The direct excitation of such molecular states with non-zero PEDMs, requires a large number of photons due to the dipole selection rules. The insight that led to the observation of the trilobite molecules with the usual two-color photoassociation schemes was the realization that some low-l states in alkali-metal atoms, such as Rb(ns) or Cs(ns) Rydberg states, have quantum defects, μ, with small non-integer parts[42]. For instance, in exciting Rb(ns) states with quantum defect μ = 3.13, the (n − 3)l ≥ 3 manifolds with vanishing μ lies close in energy to the ns state. An infinitesimal mixing of the degenerate (n − 3)l ≥ 3 manifold with the ns state introduces an opposite parity component into the ns wave function, rendering two-photon excitation of homonuclear Rydberg molecules with non-vanishing PEDM possible.

The mixing of opposite parity degenerate hydrogenic manifold into ns states in Rb was conclusively shown to lead to formation of Rb(5s − ns 3Σ) Rydberg molecules with Debye-size PEDM[43]. In rubidium, the mixing of the degenerate manifold into Rb(ns) states is at best 0.1%. The PEDMs scale as (n − μ)−2; this scaling was confirmed in the experiments[43].

Two unique features in Cs contribute to produce strong mixing and hence large PEDM: the Cs(ns) small non-integer quantum defect, μ = 4.05, causes the (n − 4)l ≥ 3 degenerate manifolds to lie energetically close to the ns states, and the p-wave resonances occur at much lower energies than the corresponding resonances in Rb, mixing in more strongly the Cs(ns) states with the close-by degenerate angular momentum manifolds, see the Section on higher-order partial-wave scattering. The resulting BO potential curves support bound vibrational levels to the blue of Cs(ns) thresholds with large PEDM[33]. Observed blue-shifted spectral features were confirmed in calculations of Rydberg molecular states above the dissociation thresholds. The first trilobite Rydberg molecules were formed in two-photon association of Cs(6s-37, 39, 40s 3Σ) Rydberg molecule states[34], see Fig. 2. Kilo-Debye PEDM for the trilobite Rydberg molecules were measured in field ionization spectra, indicating state mixing as high as 90%.

Fig. 2

Trilobite Rydberg molecules in an ultracold Cs gas. a Calculated electron density distributions for the 3Σ+ states of Cs correlating to 37s + 6s, 39s + 6s, and 40s + 6s. b The corresponding potential energy curves for the states shown in a with calculated vibrational levels. c The broadening of the vibrational levels indicated by an arrow in b as a function of applied electric field. The changes in the lineshape as a function of electric field are used to determine the dipole moments indicated in the figure. Figure from ref. [34]. Reprinted with permission from AAAS

Higher-order partial-wave scattering: The pseudopotential method in Eq. (3) has been particularly powerful in predicting the pressure broadening and shifts in high principal quantum number excitations, n > 20, when the ground state atom involved is a noble gas. However, chinks soon appeared in its potent armor. This situation was rectified by Omont[44], some 40 years after Fermi’s seminal contribution[5], who showed that the origin of deviations from Fermi’s prediction was in the electron-scattering phase shift. In cases with low n, when the electron kinetic energy, KE = E + e2/R, varies as a function of the perturber position R, a proper description of the electron–neutral interaction must include the energy dependence of the s-wave scattering length,  = , where is the electron momentum in atomic units.

In many but not all cases, the s-wave contact interaction is adequate to describe the observed Rydberg molecular features[32,43]. Omont showed how to incorporate higher partial-wave interactions, specifically, and most importantly, for alkali-metal atoms, where electron-scattering in the p-wave is resonant near threshold. In electron-alkali atom scattering, a metastable 3P negative ion exists at fairly low energies (e.g., ~ 0.03 eV in rubidium[45-47]).

The p-wave contribution to the electron-atom interaction potential takes the form[44],where  = − is the energy-dependent p-wave scattering volume, and δ(k) is the energy-dependent p-wave scattering phase shift. The arrows over the ∇ operators indicate the direction in which the operator acts. As with the s-wave interaction, the strength of the p-wave scattering length is dependent on the kinetic energy of the electron. Note that when the scattering phase shift becomes resonant, δ(k) = π/2, the p-wave contact interaction strength diverges. This divergence happens exactly at the energy of a metastable state of an electron bound to the neutral perturber. Further details on p-wave scattering effects on Rydberg molecular states are discussed in Box 2. An immediate consequence of the strong resonant p-wave interaction is the formation of a new class of Rydberg molecules, the so-called “butterfly” molecules[48,49]. The butterfly Rydberg molecules, like their larger and bigger brethren, the trilobite Rydberg molecules, possess PEDMs, albeit smaller ones. The butterfly Rydberg molecules were recently observed in ref. [39].

Spin coupling: In the preceding discussion, interactions between a ground state atom and a Rydberg electron were assumed to be spin-independent. While the scattering phase shifts did separately depend on the spin of the electrons, the scattering amplitudes added up incoherently. In alkali-metal atoms, the two valence electrons couple into singlet and triplet spinors: S = s + s = 0, and 1, with s and s, respectively, the Rydberg atom electron and ground atom electron spins. The spin–orbit (SO) interaction couples the Rydberg orbital angular momentum l to s, with a coefficient ASO. The hyperfine (HF) interaction couples the nuclear spin of the ground state atom, i, to s, with an isotropic coefficient AHF. The HF interaction in the Rydberg states and SO interaction in the Rydberg degenerate manifolds are ignored[11]. The scattering electron angular momentum, L, also couples to S and gives rise to partial-wave scattering amplitudes with the total electron spinor, J = L + S. For p-wave scattering, i.e., L = 1, three fine-structure phase shifts (J = 0, 1, 2) result. Typical BO interaction potential energies, which result from including the relativistic spin terms are shown in Fig. 3, where a new interaction curve whose spin character is mixed singlet/triplet, emerges.

Fig. 3

Singlet/triplet mixing in Cs Rydberg molecule BO potential energy curves. The interaction potentials including spin-dependent relativistic interactions (see Box 3 for details). The lowest potential is mainly a S = 1 (triplet) BO potential for a Cs(6s1/2) - Cs(32p3/2) Rydberg molecule in the ground hyperfine state, F = 3. The excited BO potential curve is predominantly a spin-mixed S = 0 and S = 1 (singlet and triplet) state. Unphysical cusps can be seen in both potentials at R ~ 1629 a.u. The cusps are the result of the semi-classical description of the energy-dependent scattering length at the transition from the classically allowed to the classically forbidden region. The cusps cannot be removed at the level of approximation discussed in this review, but they have little influence on the nature of molecular states of interest here. These potentials were recently employed in ref. [51] for interpretation of the observation of spin-mixed states in ref. [38]

The details of the spin-dependent Hamiltonian and the different interaction terms are given in Box 3. The first calculations of Rydberg molecules, including the hyperfine interaction, but neglecting singlet scattering, were performed by Anderson et al.[36] for Rb. Already, in the observed spectrum of Rb(35,37d5/2)F = 2 Rydberg molecules[50], there were features which matched the v = 0 vibrational levels in the mixed singlet/triplet potential energy curves. The observation and interpretation of singlet/triplet mixing due to the Cs ground-state HF interaction was made in ref. [38]. This interpretation neglected the important p-wave scattering contribution in Cs, which was remedied in ref. [51]. Mixed molecular states were observed and interpreted in Rb as well,[39,52-55]. Kleinbach et al.[53] observed blue-shifted states in Rb due to resonant spin–orbit effects; Tallant et al.[33] had observed the blue-shifted states in Cs due to strong non-adiabatic effects.

Box 2 P-wave resonances couple up potential energies

The effect of p-wave electron-atom scattering on the molecular Born-Oppenheimer potentials of an ultralong-range Rb2 Rydberg molecule is shown in the figure. A set of (on this scale) flat potentials exist at each Rydberg energy threshold. Plunging down through these potentials, a p-wave dominated potential is seen starting at the n = 30 hydrogenic threshold. This potential cuts through the 33s, 31d, and 32p thresholds before turning back up at shorter range. This new potential corresponds to an attractive 1/R Coulomb potential between the Rydberg Rb+ core and a metastable Rb− ion.

The fact that the p-wave scattering is not isotropic, leads to an interaction that depends on the 3-dimensional gradient of the Rydberg molecule wave function at the location of the perturber atom. The p-wave interaction maximizes this gradient of the electron wave function. The fact that the derivative acts in all three spatial dimensions gives rise to two types of potentials; those where the derivative is maximum along the internuclear axis resulting in Σ type molecules whose potentials are illustrated in panel a, and those where the derivative is maximum perpendicular to this axis resulting in Π type molecules whose potentials are illustrated in panel b. Σ molecules also experience the s-wave electron-perturber interaction meaning that the deep p-wave potentials also couple to s-wave dominated potentials seen near each Rydberg threshold. Π type molecules have a zero in the radial portion of the Rydberg electron wave function, and thus do not experience the s-wave electron-perturber interaction.

As the ionic potential plunges down, it goes through a broad crossing with the next nearly degenerate manifold of high-angular momentum (l > 2) hydrogenic Rydberg states (n = 29 in the case of the figure). The upper branch of this crossing is highlighted in the inset of panel a. In the Σ Rydberg molecules, the radial Rydberg oscillation forms wells that are deep enough to bind localized vibrational states creating a new kind of p-wave interaction dominated molecule, the so-called “butterfly” molecules[48,49].

Fabrikant et al.[106], included the different fine-structure terms of the 3P phase shifts in the ns3Σ states. The contributions from the J = 0 and J = 2 phase shifts in the ns3Σ curves produce oscillations in the potential curves, with mainly Σ molecular character, while the potential energy curve resulting from the 3P1 scattering retains mainly a Π molecular state, diabatically crossing other curves. The mostly Σ curves in the intermediate-R regions support the butterfly Rydberg molecular states

P-wave electron-atom interactions. The Rb Rydberg dimer potentials near the Rb(31d) threshold are shown for cases where the p-wave interaction is along (a), and perpendicular to (b) the internuclear axis. In a the normal Rydberg molecule potentials dominated by the s-wave electron-Rubidium interaction go through a set of narrow avoided crossings with a p-wave interaction dominated potential that falls sharply as the internuclear separation decreases. Similar behavior can be seen in b, however, because a node has been enforced in the angular portion of the Rydberg electron interaction, there is no s-wave electron-atom interaction and hence no oscillation in the potentials. The inset in a shows a new set of p-wave dominated potentials and the vibrational states created by them, which correspond to the butterfly type Rydberg molecule. Figure from ref. [48], © IOP Publishing. Reproduced with permission. All rights reserved

Box 3 Spins mix things up

The relevant spin-dependent relativistic terms in the Hamiltonian for the ultracold Rydberg-atom ground-atom system have been incorporated in successive approximations[50,51,55,106]. All of these contributions are encapsulated in the Hamiltonian,from ref.[55]. Here, M is the projection of the electron-scattering angular momentum on the internuclear axis; for Rb and Cs Rydberg excitations, scattering of the Rydberg electron from the ground state Rb or Cs atoms produces M = 0(Σ) and M = 1(Π) molecular symmetries, and is the M = 0, ±1 spherical component of the gradient operator acting in the indicated direction. The total spin angular momentum projection is given by M. The projection of spatial and spin electron angular momentum along the internuclear axis is given by M = M + M,

is the bare Rydberg Hamiltonian and the other terms in describe pairwise couplings between the various angular momentum quantities present in the system, shown schematically in the figure. The first coupling term is the spin–orbit interaction for the Rydberg electron and the next coupling term is the hyperfine interaction in the ground state atom. The remaining terms are sums over different internal angular momentum states, and incorporate the s- and p-wave electron-perturber interactions which are dependent on the total electronic spin S = s + s.

Because there is a non-zero relative orbital angular momentum between the Rydberg electron and the ground-state perturber, there is an additional spin–orbit effect. This additional spin–orbit effect results in a J = L + S dependent p-wave interaction that couples the total spin with the relative orbital angular momentum L[51,55,106]. It should be noted here that within the J dependent, p-wave interaction there is a subtle additional coupling between relative orbital angular momentum L and the total electron spin, S, that results from these two vectors precessing about J during the interaction. This coupling was first noted in ref. [106] and later reformulated in the form of contact interactions as discussed in ref. [55]. The precession leads to additional coupling terms between the Σ and Π type molecular potentials that becomes strong near the p-wave interaction resonances. One result of these pairwise angular-momentum coupling terms is that the only conserved angular momentum (in the body-fixed molecular frame) is the total projection along the internuclear axis, , where is the magnetic projection of the ground state nuclear spin

Pairwise angular momentum coupling. A diagram of relevant spin coupling terms in Rydberg atom - ground atom interaction. The red and blue balls indicate the positions of the Rydberg ionic core and the ground perturber atom, respectively. V signifies the spin-dependent scattering interactions. All other spin terms are indicated. Because the Rydberg molecules form at large internuclear distances, hyperfine interaction in the ground perturber atom can be thought of inducing long-range spin-mixing in the molecule. Here, the dashed curve encloses the spins coupled by the spin–orbit interaction, the dotted curve encloses the spins coupled by the hyperfine interaction, and the solid curve encloses the spins correlated by the spin-dependent Rydberg electron-ground state atom interaction

Realization of Rydberg macrodimers

Rydberg macrodimers, Fig. 1, are formed by pairing Rydberg atoms with attractive long-range and repulsive shorter range interactions such that potential wells can form. Due to the large polarizability of Rydberg atoms and high density of states, the bond lengths can be larger than 1 μm. The potential wells that bind the atoms together are broad on the molecular scale with well depths ranging from MHz to GHz. The multipolar interactions between pairs of Rydberg atoms that lead to the formation of macrodimers were treated in detail in refs. [14,15,19]. Broader description of Rydberg atom-Rydberg atom interactions can be found in refs. [16,21,22,56-58].

Rydberg macrodimers were first predicted by Boisseau et al.[23]. Subsequent calculations were carried out in[19,25,26,30,59]. Macrodimers were first observed by Overstreet et al.[24] and later by Saßmannshausen et al.[27]. Rydberg macrotrimers have also been predicted[60].

In the original Rydberg macrodimer description[23], interatomic potential energy wells supporting bound states were calculated to exist for a combination of attractive quadrupole (R−5) and repulsive vdW (R−6 and R−8) potentials, Fig. 4. The equilibrium distance for these molecules scales as 0.3n3a0; for n ~ 70 the equilibrium distances are >1 μm. The size of the Rydberg electronic wave function is smaller than the equilibrium distance, i.e., the potential wells lie outside the LeRoy radius[15,23]. As a consequence, exchange interactions between the Rydberg electrons can be neglected. In the case of Rydberg macrodimers the LeRoy radius is overly restrictive because the orbitals of the electrons are so diffuse.

Fig. 4

Cs Rydberg macrodimer BO potential energy curves. The Cs Rydberg macrodimer potentials in the region of the 89d + 89d asymptote in the presence of an electric field, ε. The total electronic angular momentum projection quantum number is m = 3. The electric field lies along the internuclear axis. Some potential energy curve crossings that support bound states are labeled, C1, C2, C1′, and C2′. The labels along the vertical axes give the angular momentum j for each atom of the pair. The potential energy curves are shown for 0 mV cm−1 (a), 30 mV cm−1 (b), 60 mV cm−1 (c) and 120 mV cm−1 (d). Figure from ref. [19], reprinted with permission from Taylor & Francis Ltd

A full description of Rydberg macrodimer formation requires that interactions between different potential energy curves be incorporated. The Rydberg macrodimers form as attractive and repulsive potentials cross and form avoided crossings[19]. The avoided crossings can occur naturally among potential curves with the same symmetry, or with an external field, Fig. 4. Complicated state crossings, arising from the interplay of different multipolar terms in the Hamiltonian, allow for tuning of the resonant Rydberg macrodimer potentials with small electric fields. Background fields of <1 V cm−1 are necessary for Rydberg states around n > 50[22,61]. Because the potential wells can be induced at larger R, the well depths are shallow, scaling as n−3, while the equilibrium distances scale as n8/3. As an illustrative example, this places the equilibrium distance >5 μm for n ≈ 90, Fig. 4[25,26].

The first observation of Rydberg macrodimers was by photoassociation and subsequent recoil ion momentum spectroscopy in ultracold Cs(n ~ 60)[24,62,63]. Cs(n ~ 41) Rydberg macrodimer states were later observed spectroscopically[27]. One of the difficulties in observing macrodimers is that the molecular potential wells are shallow and broad when compared to typical diatomic molecules. The vibrational frequencies are small, <10 MHz at n ~ 40 and <100 kHz at n ~ 90. As a consequence, the spectra of a Rydberg macrodimer appears as a continuum when spectral broadening effects such as rotation, radiative and blackbody decay, and laser spectral bandwidths are taken into account. Normally, rotation of a macrodimer can be neglected since the rotational constant of such a large molecule is small, ~ 3 Hz at n ~ 90, but because a large number of rotational states can be populated, there will be rotational broadening, which causes the vibrational states to overlap at larger n[19].

A continuous spectrum is not necessarily a difficulty, but in the case of macrodimers it becomes a challenge because the Rydberg atom interaction potentials are so flat that it becomes possible to excite Rydberg pairs on repulsive curves, inducing a photodissociative collision. The observed features in the spectra of ultracold Rydberg gases that asymptotically corresponded to the excitation of Rydberg atom pairs, or states shifted away from atomic Rydberg states, have been difficult to interpret, in particular to identify macrodimers[14,59,63-67]. Recoil ion momentum spectroscopy, Fig. 5 is able to distinguish the different cases, since photoinduced collisions convert potential energy of the complex into kinetic energy of the product atoms, while bound states remain at relatively well-defined R over the lifetime of the molecule. The signatures of these different behaviors can be distinguished by measuring the ion recoil after field ionization of the Rydberg atoms[16,24,63], Fig. 5. The spectral lineshapes of Cs Rydberg macrodimers have now been observed for n ~ 41[27]. This results in a lineshape where a peak exists at the potential minimum, decreases with vibrational excitation and increases as the Rydberg density of states grows at the dissociation limit[19,27].

Fig. 5

Kinetic energy spectra from the fragmentation of macrodimers. The kinetic energy as a function of the delay, τ, between when a particular molecular complex was excited and when it was probed via field ionization. In the leftmost panel a state correlating to 88d + 90d is first excited. The other figures correspond to states correlating with 63d + 65d, 64d + 66d, 65d + 67d, and 66d + 68d. The dashed line in the leftmost panel shows the kinetic energy expected for a pair of Cs atoms photoassociated at R = 2.8 μm and moving uniformly apart at a velocity of 17 cm s−1. The other four panels show the same quantity for macrodimer states for an electric field of 224 mV cm−1 for 63d + 65d, 205 mV cm−1 for 64d + 66d, 190 mV cm−1 for 65d + 67d and 158 mV cm−1 for 66d + 68d. The red points are the data while the dashed line is what is expected for a laser induced collision of the pair of Cs atoms with a recoil energy equal to the thermal energy of the atoms in the ultracold atom trap. Figure from ref. [24]. Reprinted with permission from Nature Publishing Group

Few-body and many-body Rydberg molecule spectroscopy

Even at the atomic densities used in the first observation of ultralong-range Rydberg molecules[32,68], there were hints that bound molecular states with more than one ground state atom could form. Calculations[69,70] showed that an additional ground state atom in the Rydberg electron orbit would bind to make a triatomic molecule in homonuclear systems such as with Rb or Cs. With increasing gas density, or increasing n, larger oligomeric molecules (tetramers, pentamers, …) may form[71,72]. These higher-order molecular states of Rb(5s − ns3Σ) were observed in, refs. [73,74], see Fig. 6a. Stretching and bending vibrations and electric field effects on triatomic Rydberg molecules were also studied[75,76]. Gaj et al.[73] found that with higher n Rydberg excitations, n ~ 80–110, so many perturber atoms could exist in a Rydberg orbit that a spectral overlap of molecular lines produced a density dependent broadening, Δ =  = , where ρ(R) is the local density from the position of the Rydberg ion and is the weighted average density, Fig. 6a.

Fig. 6

Polyatomic Rb molecules and many-body spectra. a The observation of polyatomic s-wave Rydberg molecules (dimers, trimers, and tetramers) in photoassociation of Rydberg molecules; figure reprinted from ref. [73], b density averaged absorption spectra for Rb(71s) excitations. In a, the gas is thermal at T = 0.5 μK and the MOT center density is ρ = 1.7 × 1012 cm−3. The laser profile illuminates the cloud over the full cylindrical radius (black) and in a small waist (dashed, blue); data is reported in Gaj et al.[73]. In b, the density averaged absorption spectrum at T/T = 0.47 for a partially condensed gas of (zero temperature) at peak density ρ = 2.3 × 1014 cm−3. Reprinted figure with permission from ref. [77] Copyright (2016) by the American Physical Society

A theory for the excitation of giant impurities in a quantum gas was developed by Schmidt et al.[77] using the functional determinant approach (FDA), Fig. 6b, c. This time-dependent method is capable of unmasking multiscale time dynamics in a quantum gas by accounting for quasiparticle excitations and molecular formation. The FDA simulations revealed the resolved molecular oligomeric excitations in a Rb BEC[78], as well as the unresolved broad features at much larger detuning, heralding the crossover from few-body to many-body polaronic states.

In Rydberg excitation in an ultracold gas, large electron-scattering phase shifts due to the presence of triplet p-wave shape resonances[47,48] introduce complicated non-adiabatic interactions between BO potential energy curves. A classical Monte Carlo simulation, sampling the energy shift from each Rydberg-perturber distance in the BO potential energies, can describe the observed spectral profiles[78]. However, this classical interpretation, which relies on p-wave resonances to produce large detuning shifts in the interaction potentials, cannot adequately describe the large n spectral tails emanating from the highest density condensed regions. The narrowing of the spectral tails is described numerically and analytically in terms of formation of a class of bosonic polarons, the Rydberg polarons. The classical interpretation misses the underlying many-particle Rydberg molecule formation.

In heavy alkaline-earth atoms (calcium, strontium, and barium), the electron-scattering in the p-wave channel binds (positive electron affinity) and there are no resonances. This atomic system provides a clean platform to study many-body Rydberg impurity dressing. Ultralong-range Sr Rydberg dimers have been observed[40]. The lifetimes of these Sr Rydberg molecules are longer than comparable states in Rb due to the aforementioned absence of p-wave shape resonances[79].

An FDA prediction for the formation of Rydberg polarons in a uniform density gas is the appearance of a Gaussian density profile, Fig. 6b, c[77]. It was shown that the Gaussian width narrows with increasing n, or the local density[80,81], as , with ρ0 the peak uniform density in the BEC. The emergence of the Gaussian profile is understood by decomposing the spectral line function in terms of excitations from a ground state BEC to interacting single-particle Rydberg molecule states. A binomial decomposition, including only two states, yields a Gaussian profile in the limit of large particle number. The observation of Rydberg excitations in Sr in ref. [81] and the detailed theory of FDA in ref. [80] reveal the intrinsic spectrum of Rydberg polarons which cannot be explained by mean field or classical interpretations[80].

The examples given in this section, show that the extreme properties of Rydberg atoms and their controllability can be used for experiments on many- and few-body physics. This line of investigation has only recently emerged.

Rydberg dimer lifetimes

Ultralong-range Rydberg molecular lifetimes for Rb, Cs, and Sr have been measured[32,34,81,82]. In Rb Rydberg molecules, two observed inelastic channels are molecular ions and l-mixing reactions[83,84]. Both reactions are largely due to the presence of strong p-wave shape resonances, producing an acceleration of the internuclear motion to the short-range “chemical” region (see Box 2 for details). The Rydberg core–atom interaction and the butterfly potential avoided crossings at short distances with degenerate manifold trilobite states produce the molecular ion and l-mixed states. At high principal quantum numbers, 90 < n < 110, a sudden increase in the lifetime of the Rydberg molecules was observed[83]. This threshold increase in lifetime could be due to many-body interaction of the Rydberg atom with the surrounding Bose gas and not influenced by the butterfly or trilobite interactions.

The absence of p-wave resonances in electron-scattering phase shifts from Sr(5s) atoms translates to generally longer lifetimes for Sr-photoassociated Rydberg molecules[81]. Even in higher density gases, where the average occupation number of perturber atoms in a Rydberg orbit is large, the molecular lifetimes are longer, on the order of the the Rydberg atom lifetime[85].

In Cs, for intermediate n quantum number excitations, the strong mixing of Cs(ns) states with the degenerate Cs((n − 4)l ≥ 3) manifolds, which produces trilobite molecules with kilo-Debye PEDM, dictates that the lifetimes are controlled by the highest nearly quantum-defect free angular momentum states. The Cs(5s − ns3Σ) molecular lifetimes were found to be influenced by the Cs((n − 4)f) state lifetimes[34].

The lifetimes of the Rydberg macrodimers are limited by the atomic Rydberg lifetimes, with spontaneous and blackbody decay[16,19,25]. Non-adiabatic processes and ionization can also limit the molecular lifetime[84], but in many cases the contribution of these decay mechanisms is much smaller than the atomic radiative lifetime[19]. Electronic autoionization is small because the equilibrium separation of the Rydberg atoms is typically large enough than there is little or no overlap of the electronic wave functions with the respective Rydberg ions. A quick look at Fig. 1 illustrates the relative size of the macrodimer equilibrium separation in comparison to the extent of the electronic wave function. Additionally, the spread of the Rydberg electron wave function and the short-range repulsive walls limit ionization. As long as the inner turning point in the bound state potential is outside the region where potentials correlating to high-angular momentum Rydberg states are, as shown in Fig. 4, the lifetimes of macrodimers will be on the order of the Rydberg atomic lifetimes. Calculations for Cs(89d) Rydberg macrodimers gave a liftetime of ~ 450 μs, determined entirely by Rydberg atom decay[19]. Measurements at densities of ~ 1010 cm−3 for photoinduced collisions and for macrodimers have affirmed this picture[24,63]. In cases where the the inner well can strongly interact with the curves with high-angular momentum character, Penning ionization can occur at a high rate, leading to molecular decay[86]. In ultracold gases where the density is larger, ~ 1012 cm−3, shorter macrodimer lifetimes have been observed[27]. The shorter lifetime was attributed to the presence of a third, nearby atom.

Giant polyatomic Rydberg molecules

The description in the Section on collisional interactions between a Rydberg atom and a ground-state atom involving the interaction of a Rydberg electron and a ground-state perturber atom invoked the concept of zero-range scattering, i.e., a(0). When the perturber is a dipolar molecule—a diatomic molecule with a PEDM, such as OH or KRb–the scattering length picture no longer holds. The first-order of interaction is charge-dipole[87],where d is the dipole moment of the molecular perturber, and N and B are the rotational angular momentum and constant for a rigid-rotor molecule, e.g., KRb. The electric field due to the Rydberg electron is , which polarizes the molecular dipole. For dipole moments above the critical dipole, known as the Fermi-Teller dipole moment, d > dcr ≈ 1.63 D, the electron may bind to the molecule, forming a negatively charged ion. When d < dcr, the anisotropic interaction of the Rydberg electron with the PEDM can form polyatomic Rydberg molecular states. These types of BO potential energy curves exhibit long-range wells and oscillatory behavior. Molecular states as described here have been calculated for a range of molecules with PEDM, e.g., Λ doublet and rigid-rotor molecules. These polyatomic potentials display two outer wells, whose electronic wave functions correspond to orientation of the molecular dipole toward and away from the Rydberg core[87].

Ultracold giant polyatomic Rydberg molecules have been studied in the context of electric field control[88], coherent control of orientation[42,87], hybridization of rotational states of dipolar molecules[89], and for conditional addressing and non-destructive readout of rotational quantum bits[90,91]. Polyatomic Rydberg molecules can potentially be realized using atomic Rb (or K) Rydberg excitation in an ultracold gas of KRb molecules.

Summary and outlook

Macrodimers and ultralong-range Rydberg molecules have exotic features: macroscopic bond length, kilo-Debye PEDM, tunability with laboratory sized fields, and fascinating electronic probability densities (trilobites and butterflies). Due to large equilibrium sizes, their rotational constants are in Hz to kHz regimes, and in experiments, localized rotational wave packets are excited. Efforts on both observational and theoretical fronts have evolved to utilize such distinct molecular properties for probes of quantum gases, for spin alignment and control, and for selective chemical reactivity in the ultracold regime. Cases where the Rydberg atom core is used to trap an electron inside a quantum gas are one example[92]. Spin-mixing of molecular potentials in ultralong-range Rydberg molecules allows for control of spin interactions and may open new vistas on the atomic analog of the celebrated central spin problem, in which a central spin interacts with a large number of strongly coupled spins in the environment. Rydberg molecules have already been investigated as probe of atomic occupation in optical lattices for quantum gas microscopy[93]. The many-body macroscopic occupation of Rydberg molecules in a BEC has been realized[77,81], but probing Pauli blocking in the many-body spectral Rydberg profile is worthy of investigation. Recently, it’s been proposed that trilobite Rydberg molecules still persist in random dense ultracold atomic gases[94].

Formation of Rydberg molecules places fundamental limits on applications in sensing and quantum information with Rydberg atoms. The interaction of Rydberg electrons with gas atoms can decohere electromagnetically induced transparency Rydberg polaritons and limit Rydberg blockade, light storage, and photon–photon interactions in ultracold Rydberg gases[95-97]. Rydberg macrodimer formation and Rydberg atom collisions can be a major contributor to dephasing and loss in experiments that require coherent light-matter interaction[28,29].

A major reason for investigating ultracold Rydberg molecules is their novel properties, such as their large dipole moments and relatively long lifetimes. Ultracold Rydberg molecules have binding energies, which are much less than their internal electronic energies. Because these two degrees of freedom are weakly coupled, ultracold Rydberg molecules can live for relatively long times allowing them to be used as spectroscopic probes of quantum gases. These properties might also be leveraged to use ultralong-range Rydberg molecules as an intermediate state for formation of optical Feshbach resonance molecules[98,99]. The extreme nature of ultralong-range molecules can give insights into quantum chemistry, as was studied in coherent photodissoaction and photoassociation with rotary echo techniques[100]. Butterfly molecules[39], by virtue or their favorable Franck-Condon overlap with low-lying Rydberg states, can be utilized to populate vibrationally bound heavy Rydberg systems, i.e., ion-pairs of ultracold atoms[101,102].

The Rydberg electrons are nearly perfect traps for the ions, so hybrid ion-atom trapping interactions could be envisioned with Rydberg molecules. It may be possible, with suitable density of Rydberg excitations, to realize ultralong-range Rydberg molecules and Rydberg macrodimers in the same trap, perhaps engineering the best properties of each type of molecule for tailored investigations. Another possible line of inquiry is to use Rydberg macrodimers to investigate gauge potentials[103-105]. Rydberg macrodimers have relatively easily controlled interaction potentials, so they offer interesting possibilities to engineer these interactions. Imaging particle dynamics can determine the exit channel trajectories if the initial state is aligned in the lab frame, thus measuring the effects of a pseudo-magnetic field (pseudo-Lorentz force on the colliding particles).

Ultracold Rydberg molecules are vehicles for applying low-energy quantum scattering techniques to the spectroscopy of quantum gases. Rydberg atom-Rydberg atom interactions are already receiving increased attention as a way to engineer collective states for many applications. As new kinds of exotic molecules, the exaggerated properties of Rydberg molecules are sure to be harnessed for new experiments.