# Probing interactions between Rydberg atoms with large electric dipole moments in amplitude modulated electric fields

## Abstract

Dipole-dipole interactions between helium atoms in Rydberg-Stark states with principal quantum number and approximately linear Stark energy shifts, resulting from induced electric dipole moments of approximately 7900 D, have been investigated experimentally. The experiments were performed in pulsed supersonic metastable helium beams, with particle number densities of up to cm. In the presence of amplitude-modulated, radio-frequency electric fields, changes in the spectral intensity distributions associated with the transitions to these states that are attributed to dipole-dipole interactions within the ensembles of excited atoms have been observed. The experimental results are in excellent agreement with calculations of the Rydberg energy level structure carried out using Floquet methods, and excitations shared by up to 4 atoms. The use of these Rydberg-Stark states as sensors for non-resonant broadband radio-frequency electrical noise is also discussed.

###### pacs:

37.10.De, 32.80.RmAtoms and molecules in Rydberg-Stark states with high principal quantum number, , can possess very large induced electric dipole moments, . For each value of , the maximum dipole moment is , where and are the electron charge and the Bohr radius, respectively gallagher94a (). These dipole moments exceed 1000 D for , and can be exploited to control the translational motion of gas-phase samples using inhomogeneous electric fields wing80a (); breeden81a (); townsend01a (); yamakita04a (); vliegen04a (); ko14a (). This has given rise to the realization of very efficient methods for acceleration, deceleration and electric trapping of atoms and molecules initially traveling in pulsed supersonic beams vliegen06a (); hogan08a (); hogan09a (); hogan12a (); allmendinger13a (); lancuba14a (). However, the large electric dipole moments of Rydberg-Stark states can also give rise to (i) strong electric dipole interactions within ensembles of atoms or molecules upon photoexcitation, or when confined in traps reinhard07a (); vogt07a (), and (ii) modifications to the energy level structure in the presence of radio-frequency (rf) electrical noise, or the effective time-dependent electric fields experienced by particles undergoing oscillatory motion in electrostatic traps hogan08a ().

We report here the results of experiments in which dipole-dipole interactions between helium (He) atoms in Rydberg-Stark states with electric dipole moments of D have been studied in amplitude-modulated electric fields. Detailed experimental studies of the interactions between atoms in these states have not previously been carried out but are of importance for the interpretation and refinement of experiments with trapped atoms and molecules. The results presented are also essential for experiments directed toward the deterministic preparation of single atoms in circular Rydberg states using the crossed-field method hyafil04a (); delande88a (); hare88a (); anderson13a (), and for planned experiments in which cold Rydberg atoms are to be transported and trapped close to superconducting microwave circuits for applications in hybrid quantum information processing hogan12b (); carter13a (); hermann14a (). In addition, the sensitivity of atoms in Rydberg-Stark states to weak rf electric fields, demonstrated here, highlights potential applications in broadband electrical noise sensing holloway14a ().

In experiments with atoms in Rydberg states that exhibit strong linear Stark energy shifts and large electric dipole moments, the interaction between pairs of atoms with dipole moments and , separated by a distance , corresponds to that of two interacting classical dipoles. The resulting interaction potential, , is therefore:

(1) |

where is the orientation of the position vector , and is the vacuum permittivity comparat10a (); reinhard07a (). This, effectively classical, dipole-dipole interaction differs from the resonant dipole-dipole interactions that arise from the large transition dipole moments associated with high Rydberg states anderson98a (); mourachko98a (); vogt06a (); saquet10 (), and the, second-order, van der Waals interactions often exploited in Rydberg excitation blockade experiments and studies of many-body effects in ensembles of cold atoms tong04a (). In the Rydberg-Stark states of He studied here, for which D, Eq. 1 leads to an average dipole-dipole interaction energy of 10 MHz for a mean atom-atom separation of m. This corresponds to a particle number density of cm. In contrast to electric-field induced dipole blockade experiments reported previously vogt07a (), the Rydberg-Stark states considered here are strongly -mixed ( is the orbital angular momentum quantum number of the Rydberg electron). It is this -mixing that gives rise to the approximately linear Stark energy shifts, and the particularly large electric dipole moments upon which the experiments rely.

A schematic diagram of the experimental apparatus is presented in Fig. 1. A pulsed supersonic beam of metastable He atoms is generated in a dc electric discharge at the exit of a pulsed valve halfmann00 (). The stagnation pressure on the high pressure side of the valve is 6 bar, and the pulse duration is 200 s. The resulting metastable He beam travels at 1950 m s. After the atoms exit the source vacuum chamber through a copper skimmer (diameter 2 mm) they pass through an ion deflection zone where ions produced in the discharge are completely filtered out. The atoms then enter a region between two parallel 70 mm 70 mm square metal plates, E1 and E2 in Fig. 1(a), separated in the dimension by 8.4 mm. In this region the atomic beam is crossed by two copropagating laser beams which drive a resonant two-photon transition. The first photoexcitation step is driven by an ultraviolet (UV) laser operated at 388.975 nm which is stabilized using a He discharge cell moron12 (). Atoms in the level are then excited to Rydberg states by an infrared (IR) laser operated at nm. Both laser beams are linearly polarized parallel to the electric field axis (axis) and have Gaussian intensity profiles with full-widths-at-half-maxima (FWHM) in the and dimensions of 95 and 125 m (204 and 80 m), respectively, in the case of the UV (IR) radiation. The UV laser power is set to 35 W (corresponding to an average intensity of 0.38 W cm) to minimize effects of saturation. The power of the IR beam is varied in the range mW (with corresponding average intensities of W cm).

For all experiments described here, photoexcitation was carried out to the outermost high-field-seeking state in the Stark manifold, which evolves adiabatically to the level in zero electric field. In an electric field V cm [red dot in Fig. 1(b)] this state has % s- and % d-character, an approximately linear Stark shift, and an electric dipole moment of D. After photoexcitation the Rydberg atoms fly into the detection region of the apparatus where a pulsed potential of +1.4 kV is applied to the metal plate E3 in Fig. 1(a), generating an electric field of V cm. This field ionizes the Rydberg atoms and accelerates the resulting He ions through an aperture in E4, which is grounded, and onto a microchannel plate (MCP) detector.

Because the electric dipole moment of this Rydberg state is so large, it is very sensitive to dc and rf electric fields. An example of this can be seen in Fig. 2 which contains spectra encompassing the transition to this state when V cm. In the absence of significant rf laboratory noise this transition exhibits a Lorentzian spectral profile [Fig. 2(a)]. However, if the electrical noise emanating from the He rf discharge cell used to stabilize the UV laser is not shielded, the spectral profile of the transition changes, with the appearance of a broadened pedestal at low spectral intensities beneath the main Lorentzian feature [Fig. 2(b)]. The origin of this lineshape can be understood by considering the effect of a weak rf amplitude modulation of the electric field on the Rydberg-Stark state. Such a field, , can be expressed as

(2) |

where is the amplitude of the modulation which occurs at an angular frequency . The oscillatory component of this field perturbs the time-independent Stark Hamiltonian and can be treated using Floquet methods shirley65 (); zhang94 (); spellmeyer97 (); spellmeyer98 (); vanditzhuijzen09 (); yoshida12 (); tretyakov14 (). This involves the construction of the Hamiltonian matrix in an basis, where is the azimuthal quantum number, and is an integer representing the number of rf photons associated with each Floquet sideband. The energies of the diagonal elements in this matrix were determined using the Rydberg formula with quantum defects , , and for the s, p, d, and f states, respectively martin87a (). States for which were included in the matrix, each with Floquet sidebands, ranging in energy from to ( is the maximum value of considered). For all calculations presented convergence was achieved for .

The dc electric field gives rise to off-diagonal matrix elements coupling states with , and , while the parallel rf field couples states with and spellmeyer97 (); spellmeyer98 (). The transition frequencies and spectral intensities were determined by calculating the eigenvalues and eigenvectors of the Hamiltonian matrix. The results of calculations carried out for V cm, MHz, and a range of values of can be seen in Fig. 2(c)-(e) (vertical bars), together with experimental data (blue dots) recorded for the same parameters. In these figures, the calculated spectral intensities have also been convoluted with Lorentzian functions with FWHM of MHz (red curves). The spectrum recorded for is that in Fig. 2(a). For low-amplitude modulation, e.g., mV cm [Fig. 2(c)], the unperturbed spectral feature is modified, with the intensity distributed over the low- rf sidebands. In this case, the sidebands dominate. As is increased to 18.0 mV cm [Fig. 2(d)], and 27.0 mV cm [Fig. 2(e)], the spectral intensity is further displaced toward higher values of .

The agreement between the results of the calculations and the experimental data in Fig. 2(c)-(e) indicates that the above approach can be extended to characterize ambient rf noise in the laboratory, such as that present when recording the spectrum in Fig. 2(b). By comparing this data to a set of calculated spectra we find that the experimentally recorded spectral profile is characteristic of the presence of an oscillating field with mV cm and MHz. We conclude that the source of this noise is the second harmonic of the rf driving field of the He discharge cell used to stabilize the UV laser. In this rf frequency regime our measurements are sensitive to changes in of V cm.

Inspired by electromagnetically induced transparency (EIT) studies of many-body effects in Rydberg ensembles mohapatra07 (); weatherill08 (), we have investigated how spectra recorded for selected values of and change depending upon the IR laser intensity and the particle number density (see Fig. 3). These experiments were performed with MHz and mV cm. Under these conditions and a low IR laser intensity ( W cm), the spectral intensity distribution exhibits a minimum at the unperturbed transition frequency arising from maximal suppression of the sideband, and intensity maxima in the sidebands [Fig. 3(c)]. As is increased to 718 W cm [Fig. 3(b)] and 1794 W cm [Fig. 3 (a)], the observed spectral intensity distributions change. Notably, in Fig. 3(a), the minimum at the unperturbed transition frequency is absent and the transitions to the most intense sidebands are saturated.

To quantify the effects contributing to the changes in the experimental spectra in Fig. 3, we calculate their dependence on . In these calculations the electric dipole transition moments from the intermediate level to each individual Floquet sideband were combined with the measured laser intensities and an unperturbed spectral linewidth of MHz. This linewidth was determined from experimental data in Fig. 3(a) where effects of saturation are negligible. Using this information, the saturation parameters (where is the single-atom Rabi frequency) were calculated for each Floquet sideband to determine the overall spectral intensity distribution function metcalf99a (). We note that this procedure does not involve any free or fitted parameters. However, to obtain quantitative agreement between the results of the calculation and the experimental data, it was necessary to introduce a constant factor of 2 by which the single-atom Rabi frequencies were scaled. This factor is independent of laser intensity and suggests that excitations in the atomic beam are shared by multiple atoms coupled via dipole-dipole interactions. By making the assumption that the resulting many-body Rabi frequency is stanojevic09 (); gaetan09 (); urban09 (), where is the effective number of atoms sharing each excitation, best agreement between the experimental data and the calculations in Fig. 3 was obtained for atoms. The deviation of the wings of the calculated spectral profiles from the experimental data in Fig. 3(a) and (b) is attributed to effects of the Gaussian spatial intensity distributions of the laser beams, and variations in the inter-particle spacing within the atomic beam. Only the spatially averaged IR laser intensities were considered in the calculations.

To confirm that the changes observed in the spectra in Fig. 3 do indeed arise as a result of many-body interactions, it is necessary to demonstrate that by reducing the atom number density the unsaturated Floquet sideband intensity distribution is recovered. This was achieved by adjusting the timings in the experiment so that the low density trailing component of the pulsed supersonic beam was probed. Spectra recorded for high and low atom number densities but a constant value of W cm) are displayed in Fig. 4 (blue dots). These measurements were performed for the same values of and as in Fig. 3, and approximately equal particle number density, in the case of Fig. 4(a). Comparison of the experimental data (blue dots) recorded at high density [Fig. 4(a)] with the calculated spectra indicates that, as in Fig. 3, each excitation is shared by atoms. In the low density regime, in which the amplitude of the signal is reduced by a factor of 34 [Fig. 4(b)], the unsaturated Floquet sideband intensity distribution is recovered with its characteristic minimum clearly visible at the unperturbed transition frequency. The calculations indicate that in the low density case excitations are associated with single atoms, i.e., . In this case the dip at the unperturbed transition frequency is approximately 80% of the intensity of the spectral maxima. We define this as the contrast which, in the experimental spectrum in Fig. 4(b), is % and in very good agreement with the results of the calculations. In the higher density case the calculated spectral intensity distributions for and 5 result in a contrast of 92, 96 and 98%, respectively. The contrast in the experimental data in Fig. 4(a) is 95% which is in closest agreement with the calculation for which .

Averaging over the angular dependence of the dipole-dipole interaction potential, Eq. 1, results in a blockade radius of 8.2 m for the unsaturated spectral half-width-at-half-maximum of MHz. If 4 atoms are considered to be located within a sphere of this radius a particle number density of cm results. This density agrees well with those previously achieved in similar pulsed supersonic atomic beams saquet10 (); halfmann00 ().

The results of the experiments and calculations reported here demonstrate the sensitivity of Rydberg-Stark states with large electric dipole moments to non-resonant rf electric fields, and to dipole-dipole interactions. The detection of changes in the spectral intensity distribution associated with Rydberg excitation in the presence of a low-frequency amplitude modulated electric fields, represents a direct spectroscopic observable which can be employed to identify many-body contributions to the Rydberg photoexcitation process, without the need to directly count the number of atoms contributing to each excitation. In addition, by acting as microscopic rf antennas, Rydberg atoms can be exploited as broadband probes of electrical noise. They are ideally suited to the characterization of rf noise at vacuum–solid-state interfaces of importance, for example, in hybrid approaches to quantum information processing hogan12b (); hermann14a (), and in chip-based ion traps turchette00 (). Combining the techniques presented here with Rydberg atom detection by EIT mohapatra07 () could, in the future, permit direct optical readout of such noise.

The observed effects of electric dipole interactions on the saturation of optical spectra encompassing transitions to Rydberg states with electric dipole moments of D demonstrates the possibility of investigating many-body physics in pulsed supersonic beams with particle number densities of cm. Of particular interest will be studies of many-body contributions in atom-molecule scattering at low collision energies. Elucidating the role, as reported here, played by these interactions within ensembles of Rydberg atoms and molecules is also crucial for the refinement of experiments directed toward investigations of the decay mechanisms of electrically trapped samples seiler11 ().

###### Acknowledgements.

This work was supported financially by the Department of Physics and Astronomy and the Faculty of Mathematical and Physical Sciences at University College London, and the Engineering and Physical Sciences Research Council under Grant No. EP/L019620/1.### References

- T. F. Gallagher, Rydberg Atoms (Cambridge University Press, Cambridge, 1994).
- W. H. Wing, Phys. Rev. Lett. 45, 631 (1980).
- T. Breeden and H. Metcalf, Phys. Rev. Lett. 47, 1726 (1981).
- D. Townsend, A. L. Goodgame, S. R. Procter, S. R. Mackenzie, and T. P. Softley, J. Phys. B: At. Mol. Opt. Phys. 34, 439 (2001).
- Y. Yamakita, S. R. Procter, A. L. Goodgame, T. P. Softley, and F. Merkt, J. Chem. Phys. 121, 1419 (2004).
- E. Vliegen, H. J. Wörner, T. P. Softley, and F. Merkt, Phys. Rev. Lett. 92, 033005 (2004).
- H. Ko and S. D. Hogan, Phys. Rev. A 89, 053410 (2014).
- E. Vliegen and F. Merkt, Phys. Rev. Lett. 97, 033002 (2006).
- S. D. Hogan and F. Merkt, Phys. Rev. Lett. 100, 043001 (2008).
- S. D. Hogan, Ch. Seiler, and F. Merkt, Phys. Rev. Lett. 103, 123001 (2009).
- S. D. Hogan, P. Allmendinger, H. Saßmannshausen, H. Schmutz, and F. Merkt, Phys. Rev. Lett. 108, 063008 (2012).
- P. Allmendinger, J. A. Agner, H. Schmutz, and F. Merkt, Phys. Rev. A 88, 043433 (2013).
- P. Lancuba and S. D. Hogan, Phys. Rev. A 90, 053420 (2014).
- A. Reinhard, T. C. Liebisch, B. Knuffman, and G. Raithel, Phys. Rev. A 75, 032712 (2007).
- T. Vogt, M. Viteau, A. Chotia, J. Zhao, D. Comparat, and P. Pillet, Phys. Rev. Lett. 99, 073002 (2007).
- P. Hyafil, J. Mozley, A. Perrin, J. Tailleur, G. Nogues, M. Brune, J. M. Raimond, and S. Haroche, Phys. Rev. Lett. 93, 103001 (2004).
- D. Delande and J. C. Gay, Europhys. Lett. 5, 303 (1988).
- J. Hare, M. Gross, and P. Goy, Phys. Rev. Lett. 61, 1938 (1988).
- D. A. Anderson, A. Schwarzkopf, R. E. Sapiro, and G. Raithel, Phys. Rev. A 88, 031401(R) (2013).
- S. D. Hogan, J. A. Agner, F. Merkt, T. Thiele, S. Filipp, and A. Wallraff, Phys. Rev. Lett. 108, 063004 (2012).
- J. D. Carter and J. D. D. Martin, Phys. Rev. A 88, 043429 (2013).
- C. Hermann-Avigliano, R. C. Teixeira, T. L. Nguyen, T. Cantat-Moltrecht, G. Nogues, I. Dotsenko, S. Gleyzes, J. M. Raimond, S. Haroche, and M. Brune, Phys. Rev. A 90, 040502(R) (2014).
- C. L. Holloway, J. A. Gordon, S. Jefferts, A. Schwarzkopf, D. A. Anderson, S. A. Miller, N. Thaicharoen, and G. Raithel, IEEE Transactions on Antennas and Propagation, 62, 6169 (2014).
- D. Comparat and P. Pillet, J. Opt. Soc. Am. B, 27, A208 (2010).
- W. R. Anderson, J. R. Veale, and T. F. Gallagher, Phys. Rev. Lett., 80, 249 (1998).
- I. Mourachko, D. Comparat, F. de Tomasi, A. Fioretti, P. Nosbaum, V. M. Akulin, and P. Pillet, Phys. Rev. Lett. 80, 253 (1998).
- T. Vogt, M. Viteau, J. Zhao, A. Chotia, D. Comparat, and P. Pillet, Phys. Rev. Lett., 97, 083003 (2006).
- N. Saquet, A. Cournol, J. Beugnon, J. Robert, P. Pillet, and N. Vanhaecke, Phys. Rev. Lett. 104, 133003 (2010).
- D. Tong, S. M. Farooqi, J. Stanojevic, S. Krishnan, Y. P. Zhang, R. Côté, E. E. Eyler, and P. L. Gould, Phys. Rev. Lett. 93, 063001 (2004).
- T. Halfmann, J. Koensgen, and K. Bergmann, Meas. Sci. Technol. 11, 1510-1514 (2000).
- F. Moron, A. L. Hoendervanger, M. Bonneau, Q. Bouton, A. Aspect, D. Boiron, D. Clèment, and C. I. Westbrook, Rev. Sci. Instrum. 83, 044705 (2012).
- J. H. Shirley, Phys. Rev. 138(4B), B979 (1965).
- Y. Zhang, M. Ciocca, L.-W. He, C. E. Burkhardt, and J. J. Leventhal, Phys. Rev. A 50, 1101 (1994).
- N. W. Spellmeyer, D. Kleppner, M. R. Haggerty, V. Kondratovich, J. B. Delos, and J. Gao, Phys. Rev. Lett. 79, 1650 (1997).
- N. W. Spellmeyer, PhD thesis, Massachusetts Institute of Technology, Cambridge, MA (1998).
- C. S. E. van Ditzhuijzen, A. Tauschinsky, and H. B. van Linden van den Heuvell, Phys. Rev. A 80, 063407 (2009).
- S. Yoshida, C. O. Reinhold, J. Burgdörfer, S. Ye, and F. B. Dunning, Phys. Rev. A 86, 043415 (2012).
- D. B. Tretyakov, V. M. Entin, E. A. Yakshina, I. I. Beterov, C. Andreeva, and I. I. Ryabtsev, Phys. Rev. A 90, 041403(R) (2014).
- W. C. Martin, Phys. Rev. A 36, 3575 (1987).
- A. K. Mohapatra, T. R. Jackson, and C. S. Adams, Phys. Rev. Lett. 98, 113003 (2007).
- K. J. Weatherill, J. D. Pritchard, R. P. Abel, M. G. Bason, A. K. Mohapatra, and C. S. Adams, J. Phys. B: At. Mol. Opt. Phys. 41, 201002 (2008).
- H. J. Metcalf and P. van der Straten, Laser Cooling and Trapping (Springer-Verlag New York, 1999).
- J. Stanojevic and R. Côté, Phys. Rev. A 80, 033418 (2009).
- A. Gaëtan, Y. Miroshnychenko, T. Wilk, A. Chotia, M. Viteau, D. Comparat, P. Pillet, A. Browaeys and P. Grangier, Nat. Phys. 5, 115 (2009).
- E. Urban, T. A. Johnson, T. Henage, L. Isenhower, D. D. Yavuz, T. G. Walker and M. Saffman, Nat. Phys. 5, 110 (2009).
- Q. A. Turchette, D. Kielpinski, B. E. King, D. Leibfried, D. M. Meekhof, C. J. Myatt, M. A. Rowe, C. A. Sackett, C. S. Wood, W. M. Itano, C. Monroe, and D. J. Wineland, Phys. Rev. A 61, 063418 (2000).
- Ch. Seiler, S. D. Hogan, H. Schmutz, J. A. Agner, and F. Merkt, Phys. Rev. Lett. 106, 073003 (2011).