Quantum gas microscopy of Rydberg macrodimers

See allHide authors and affiliations

Science  17 May 2019:
Vol. 364, Issue 6441, pp. 664-667
DOI: 10.1126/science.aaw4150

Bonding's outer limit

In a Rydberg state, an atom has been very nearly, but not quite, ionized. This puts the electron relatively far from the nucleus, and two atoms in such a state can thus form a rather long-range bond. Hollerith et al. observed this phenomenon in fine detail by exciting pairs of ultracold rubidium atoms along the diagonal of an optical lattice. The authors resolved the vibrational state structure spectroscopically and showed that the Rydberg dimers manifested bond lengths exceeding 700 nanometers.

Science, this issue p. 664


The subnanoscale size of typical diatomic molecules hinders direct optical access to their constituents. Rydberg macrodimers—bound states of two highly excited Rydberg atoms—feature interatomic distances easily exceeding optical wavelengths. We report the direct microscopic observation and detailed characterization of such molecules in a gas of ultracold rubidium atoms in an optical lattice. The bond length of about 0.7 micrometers, comparable to the size of small bacteria, matches the diagonal distance of the lattice. By exciting pairs in the initial two-dimensional atom array, we resolved more than 50 vibrational resonances. Using our spatially resolved detection, we observed the macrodimers by correlated atom loss and demonstrated control of the molecular alignment by the choice of the vibrational state. Our results allow for rigorous testing of Rydberg interaction potentials and highlight the potential of quantum gas microscopy for molecular physics.

A quantitative determination of the structure of molecules is an essential goal of physical chemistry and is crucial to revealing and understanding their properties. The high level of quantum control and the ultracold temperatures achieved in atomic physics provide tools to study molecules and their structure (1). Prominent examples include the observations of weakly bound Feshbach molecules (2, 3), the controlled photoassociation of individual molecules in a microtrap (4), molecules comprising ground-state atoms bound to a highly excited Rydberg atom (57), or pure long-range molecules (8, 9) that are bound purely electrostatically. The binding mechanism of the latter, where the electron orbitals of the constituents do not overlap, has also been predicted to pertain between two Rydberg atoms. These “Rydberg macrodimers” (1014) are truly remarkable in their macroscopic bond lengths, which are 10,000 times those of usual diatomic molecules and thus reach typical interparticle distances in magneto-optical traps (15), optical lattices (16), or optical tweezers (17, 18). Their enormous size not only enables direct optical access to individual constituents, but also allows for the controlled binding of two atoms optically pinned at the correct distance. First signatures of Rydberg macrodimers have been observed in systems of laser-cooled atoms by spectroscopy (15) and pulsed-field ionization (19), but a vibrationally and spatially resolved detection has been lacking so far.

Here, we present a precise study of Rydberg macrodimers starting from ground-state 87Rb atoms deterministically arranged in an optical lattice. We probed the vibrational levels by two-photon spectroscopy, resolving more than 50 excited vibrational resonances. The observed spectrum agrees quantitatively with ab initio calculations of Rydberg interaction potentials (2022), providing a stringent test for their accuracy. Using the site-resolved detection and single-atom sensitivity of our quantum gas microscope (23, 24), we identified the macrodimer signal microscopically as a loss of pairs of ground-state atoms at a distance of a bond length. Furthermore, we controlled the spatial orientation of the photoassociated molecules by the parity of the vibrational wave function and the polarization of the excitation laser.

Rydberg interaction potentials are of van der Waals type for asymptotically large interatomic separations R (25), but the situation is more complicated at smaller distances. In a generic situation, repulsive interactions may increase the energy of a lower-lying pair state, |α˜, while attractive interactions decrease the energy of a higher-lying one, |β˜. At some specific distance, the two potentials become degenerate, and any finite coupling between them opens a gap (Fig. 1A), which is large when |α˜ and |β˜ contain large amplitudes of mutually dipole-dipole–coupled states. The resulting potential minimum hosts a series of bound macrodimer states Φν(R) (10, 11, 15), where ν denotes the vibrational quantum number of the nuclear motion. In our experiment, we chose such an avoided crossing of two 0g+ potentials with gerade symmetry and zero angular momentum projection on the interatomic axis (1, 20). For large distances, the selected pair states transform into the states |α˜|e,e|35P1/2,35P1/2 and |β˜|e,e|35P1/2,35P3/2 (26), which can be optically coupled from the ground state |g, g〉 by a two-photon transition. Here, |g〉 = |5S1/2F = 2, mF = 0〉 and the bond length of this macrodimer state is predicted to be 720 nm, close to the diagonal lattice spacing of Rin = 752 nm in our optical lattice. With the atoms initially prepared in the motional ground state of the optical lattice, this coincidence of length scales results in a strong optical coupling due to the large wave function overlap.

Fig. 1 Schematic of the experiment.

(A) The two-photon and two-atom excitation from the ground state |g, g〉 occurs off-resonantly with detuning Δ via intermediate states where one atom is excited to the Rydberg state |e〉 ≡ |35P1/2〉. The avoided crossing of two coupled Rydberg interaction potentials |α˜ and |β˜ (blue and red solid lines) leads to the formation of a binding potential, whose molecular states can be laser-excited if 2Δ matches the interaction shift from the asymptotic pair state |e, e〉. (B) The potential well hosts bound states (gray horizontal lines) with an energy spacing of ~6 MHz. The vibrational wave functions, indicated in orange, are much narrower in the internuclear distance R than the RMS width of 2σlat of the initial relative wave packet Φg(R) in the lattice. The lower panel shows the amplitudes cα˜,β˜(R) resulting from a spatially dependent decomposition of the binding potential. (C) The atoms are initially arranged in a regular grid of lattice sites with a lattice constant alat = 532 nm and illuminated by the UV laser with wave vector k. (D) An exemplary fluorescence image of an atomic cloud illustrates the correlated losses due to molecule formation (left), indicated as orange symbols in the reconstructed image (right).

Our experiments started with a two-dimensional atomic Mott insulator of 87Rb with a lattice filling of 94 ± 1% in the atomic ground state. The atoms were pinned in a deep optical lattice with a root mean square (RMS) width σlat = 68 nm of the motional ground state in the atomic plane (Fig. 1B) with a temperature below the on-site trapping frequency. The molecules were photoassociated by an ultraviolet (UV) Rydberg excitation laser at a wavelength of 298 nm propagating along the diagonal direction of the optical lattice, with linear polarization aligned either in or orthogonal to the lattice plane. Our typical optical Rabi coupling from the state |g〉 to |e〉 was Ω/2π = 1.2 ± 0.1 MHz. We detected the excited macrodimers as missing pairs of ground-state atoms as these were ejected from the optical lattice very efficiently as a result of kinetic energy release in the macrodimer decay. The remaining ground-state atoms were then imaged with near-unity fidelity using a quantum gas microscope (Fig. 1D).

We first aimed to identify the presence of the bound macrodimers by their spectroscopic fingerprint. To this end, we illuminated the atomic ensemble for tUV = 100 ms with the UV light polarized in the atomic plane for varying detuning Δ from the bare Rydberg state |e〉. In each step, we swept the optical frequency of the excitation laser by 480 kHz during the illumination time to ensure coverage of the full spectral region between neighboring data points. We observed the coupling to the bare Rydberg state |e〉 as a very broad saturated loss resonance (Fig. 2). At negative detunings, the resonance features an asymmetric broadening due to coupling to attractively interacting pair-state potentials (27). At positive detunings around Δ/2π = 180 MHz, the first macrodimer bound state becomes two-photon resonant. A nonsaturated high-resolution spectrum of this line (Fig. 2C) features a full width at half maximum of 139 ± 5 kHz, which is of the same order as the measured linewidth of the bare Rydberg resonance (28). For higher vibrational resonances, we observed a reduction of the line strength, which we attribute to a combination of increased intermediate-state detuning and reduced overlap of the spatial wave functions. In addition, we observed a suppression of the excitation to odd vibrational states ν due to the approximate odd parity symmetry of the molecular wave functions with respect to the equilibrium distance. This is consistent with the Franck-Condon principle, which predicts the coupling to be proportional to the overlap integrals of the broad initial and the tightly confined final spatial wave functions Φg(R) and Φν(R) (Fig. 1B). However, a closer inspection of the experimental data shows that this simple picture needs to be refined. Repeating the spectroscopy with orthogonal polarization (so that the optical electric field oscillates out of the atomic plane) results in a suppression of the line strength of the even lines.

Fig. 2 High-resolution spectroscopy detuned from the Rydberg resonance.

(A) Wide-range spectroscopy signal. We find no signatures of bound states on the interaction-broadened red-detuned side of the Rydberg resonance [the two isolated features are due to lattice-induced Raman resonances (26)]. On the blue-detuned side, we observe clear dips with regular spacing and alternating line strength due to the coupling to macrodimers. For very large detuning, a second series of molecular bound states belonging to a 1u potential becomes resonant. (B) Zooming into the frequency region between 180 and 320 MHz reveals a spacing of the vibrational resonances of ~3 MHz, which slightly decreases for higher vibrational states. We find very good agreement with the theoretical predictions (orange lines), but the line strengths are saturated in the experiment (see fig. S8 for a nonsaturated measurement). (C) High-resolution spectroscopy of the lowest vibrational level. Error bars denote SEM.

Our microscopic access provides direct in situ information about the spatial alignment of the associated molecules and valuable insights into the underlying coupling mechanism (26, 29, 30). We compared different molecular lines by illuminating the cloud with UV light, resonant with a given vibrational state ν, until the filling of the lattice decreased to roughly 87 ± 1%. For a quantitative analysis, we evaluated spatially averaged density-density correlations g(2)(i,j)=n^k+i,l+jn^k,ln^k+i,l+jn^k,lk,l for the measured spatial atom distributions on a region of interest of 9 × 9 lattice sites (Fig. 3). Correlations show a clear peak at a distance of a lattice diagonal, revealing the bond length of the molecule. Moreover, we controlled the orientation of the photoassociated molecules by choosing a combination of vibrational quantum number and polarization of the light field. For even oscillator states ν, the correlations are stronger along the lattice diagonal parallel to the polarization of the excitation light. The molecular orientation, however, flips when considering odd oscillator states for which the dimers form predominantly along the direction perpendicular to the polarization.

Fig. 3 Macrodimer imaging and molecular orientation on the lattice.

(A) Evaluating g(2)(i, j) after illuminating the atoms with UV light polarized in the atomic plane for various macrodimer lines ν reveals a directionality in the excitation rate. For even vibrational states, the excitation rate is stronger along the polarization of the excitation laser, whereas odd states favor molecule formation perpendicular to the polarization. For both cases, we show an exemplary image from our microscope. The origin of the alternating molecular orientation is the polarization of the light with respect to the quantization axis of the dimer. The light is π-polarized with both axes aligned and σ±-polarized in the perpendicular case. (B) For out-of-plane polarization, there is no longer a preferred direction and we observe equal correlations for both diagonal directions.

The key to understanding this striking alternation in the orientation of the molecules is the interplay of electronic and motional degrees of freedom. The Born-Oppenheimer wave function of the macrodimer can be expressed as |ΨMolν(R)=Φν(R)|χel(R). Whereas at large distances, the electronic part |χel(R)=cα˜(R)|α˜+cβ˜(R)|β˜ (11, 12) is dominated by the state |β˜, for short distances the |α˜ contribution dominates (Fig. 1B), and this parametric dependence of the electronic potential has to be taken into account in the optical excitation. The two-photon Rabi coupling Ω˜ν from the ground state to the macrodimer states thus splits into two terms. Neglecting the weak spatial dependence of the two-photon Rabi couplings Ω˜α˜ (Ω˜β˜) to the states |α˜ (|β˜), we obtain Ω˜νΩ˜α˜fα˜ν+Ω˜β˜fβ˜ν with the generalized Franck-Condon integrals fα˜ν=Φν*(R)cα˜*(R)Φg(R)dR and fβ˜ν (26, 29). Because of the change of the pair state amplitudes around the potential minimum (Fig. 1B), fα˜ν and fβ˜ν have the same sign for even ν but an opposite sign for odd ν. The electronic contributions Ω˜α˜ (Ω˜β˜) depend on the alignment of the polarization relative to the molecular axis because the states |α˜ (|β˜) obey molecular symmetry constraints. In the case where the axes are parallel, Ω˜α˜ and Ω˜β˜ have the same sign and therefore both terms in Ω˜ν add constructively for even vibrational states, leading to a dominating signal in g(2)(1, 1). For the perpendicular case, the sign of Ω˜β˜ flips and constructive interference occurs for odd vibrational states, resulting in a stronger value for g(2)(1, –1). Moreover, for the measurement with light polarized out of plane (Fig. 1A), the spatial signal exhibits isotropic correlations because in that case the polarization is perpendicular to both lattice diagonals.

Although most of the line positions of the measured spectrum shown in Fig. 2 agree with the theoretical model, we find deviations for low-lying oscillator states. A finer scan of that region is shown in Fig. 4A. These deviations originate from a third pair state |γ˜ asymptotically corresponding to the optically uncoupled state |32D3/2, 37P1/2〉 intersecting the binding potential around the potential minimum. At the degeneracy point, the weak dipole-quadrupole coupling between the intersecting potentials opens another gap energetically comparable to the vibrational energy. As a consequence, a separation of the vibrational motion and interatomic interaction is no longer possible. For the theory used to describe the coarse vibrational structure shown in Fig. 2, we only accounted for the crossing formed by the coupling between |α˜ and |β˜. We extended our theory by allowing for the vibronic coupling between the vibrational modes and the electronic states |α˜, |β˜, and |γ˜. The modified eigenenergies are indicated as orange lines in Fig. 4B. Using the refined theory, we indeed can identify almost all observed lines. To confirm the effect of the intersection, we repeated our spectroscopic measurements for the lowest states in the analogous potential for n = 36, where such an additional crossing is absent. In this case, we observed a pure harmonic oscillator–like spectrum in excellent agreement with the calculations. The remarkably high sensitivity of the measured line structure to even weak modifications of the interaction potentials underlines the promise of Rydberg macrodimer spectroscopy for benchmarking Rydberg interaction potentials. This also holds for Rydberg interactions in the presence of applied magnetic fields, where accurate calculations are more difficult (26).

Fig. 4 Breakdown of the Born-Oppenheimer approximation.

(A) A closer look at the potential for n = 35 reveals a 12-MHz-wide gap, which is absent for n = 36. (B) We find a modification of the coarse spectrum discussed in Fig. 2B in the vicinity of the gap (gray shaded region). A theoretical treatment including the gap (orange lines) allows for an assignment of most of the observed lines. (C) For n = 36, the measured spectrum agrees well with the theoretical expectation. The excitation light was polarized in the atomic plane; error bars denote SEM.

In the future, the coupling to macrodimers could be used to realize quantum gates at well-defined qubit distance or to enhance Rydberg dressing schemes (28, 31), where the pair state admixture of the strongly interacting doubly excited state is enhanced relative to the singly excited intermediate states. The strongly spatially dependent loss revealed in the correlation measurements and also by modulating the initial atom distribution (26) could be used to engineer dissipatively stabilized few- or many-body states (32, 33). Furthermore, the approach demonstrated here can readily be extended to study multi-atom bound states (34, 35). Finally, bringing the coupling rate to the macrodimers closer to the decay rate of the individual Rydberg atoms may allow for the observation of novel many-body physics arising from spatial constraints and coherent interactions.

Supplementary Materials

Materials and Methods

Figs. S1 to S9

Tables S1 and S2

References (3741)

References and Notes

  1. See supplementary materials.
  2. The data displayed in the figures of this manuscript are available from the Edmond repository of the Max Planck Society, doi 10.17617/3.29.
Acknowledgments: We thank all contributors to the open-source programs “pair interaction” and “ARC” as well as R. Côté, W. D. Phillips, N. Šibalić, and J. Deiglmayr for valuable discussions. Funding: We acknowledge support by the DNRF through a Niels Bohr Professorship for T.P. and funding by MPG. This project has received funding from the European Union’s Horizon 2020 research and innovation program under grant agreement 817482 (PASQuanS) and the European Research Council (ERC) 678580 (RyD-QMB) and also from project 319278 (UQUAM) of the Seventh Framework Programme. We also acknowledge funding from Deutsche Forschungsgemeinschaft (project BL 574/15-1) within SPP 1929 (GiRyd). Author contributions: S.H. and J.R. acquired the data underlying this study and, together with J.Z. and A.R.-A., they maintained and improved the experimental setup. V.W. and T.P. contributed the theoretical predictions. I.B. and C.G. supervised the study. All authors worked on the interpretation of the data and contributed to the final manuscript. Competing interests: The authors declare no competing interests. Data and materials availability: The data that support the plots presented in this paper are publicly available from the Open Access Data Repository of the Max Planck Society (36).

Stay Connected to Science

Navigate This Article