## Engineering Rydberg-pair interactions in divalent atoms with hyperfine-split ionization thresholds

Frederic Hummel, 1, ∗ Sebastian Weber, 2, † Johannes M¨ ogerle, 2 Henri Menke, 3 Jonathan King, 1 Benjamin Bloom, 1 Sebastian Hofferberth, 4, ‡ and Ming Li 1, §

1 Atom Computing, Inc., Berkeley, California

2 Institute for Theoretical Physics III and Center for Integrated

Quantum Science and Technology, Universit¨t a Stuttgart, Germany

3 Department of Physics, Friedrich-Alexander-Universit¨t Erlangen-N¨rnberg, Germany a u

4 Institute of Applied Physics, University of Bonn, Germany

(Dated: August 2, 2024)

Quantum information processing with neutral atoms relies on Rydberg excitation for entanglement generation. While the use of heavy divalent or open-shell elements, such as strontium or ytterbium, has benefits due to their optically active core and a variety of possible qubit encodings, their Rydberg structure is generally complex. For some isotopes in particular, hyperfine interactions are relevant even for highly excited electronic states. We employ multi-channel quantum defect theory to infer the Rydberg structure of isotopes with non-zero nuclear spin and perform non-perturbative Rydberg-pair interaction calculations. We find that due to the high level density and sensitivities to external fields, experimental parameters must be precisely controlled. Specifically in 87 Sr, we study an intrinsic F¨rster resonance, unique to divalent atoms with hyperfine-split thresholds, o which simultaneously provides line stability with respect to external field fluctuations and enhanced long-range interactions. Additionally, we provide parameters for pair states that can be effectively described by single-channel Rydberg series. The explored pair states provide exciting opportunities for applications in the blockade regime as well as for more exotic long-range interactions such as largely flat, distance-independent potentials.

Keywords: Rydberg atoms, hyperfine interaction, pair-interaction potentials, F¨rster resonance o

## I. INTRODUCTION

The long-range interactions between Rydberg atoms cause exotic processes in cold quantum gases and enable applications in non-linear quantum optics [1-3]. Coherent Rydberg excitations are essential to quantum information processing using neutral atoms in order to introduce and mediate entanglement [4-6]. Ensembles of individually trapped atoms in optical tweezer arrays have enabled incredible progress for quantum simulation and computation with defect-free assembly [7, 8], reconfigurable lattice geometries [9-12], mid-circuit erasure conversion and readout [13, 14], continuous reloading [15, 16], and unseen connectivities by aid of rearrangement [17]. Recently, a variety of atomic species has been explored for quantum gas experiments and quantum technology applications exploiting different aspects of the atomic structure [18-22]. While alkali metals benefit from a long history of scientific research, the presence of a second valence electron in alkaline-earth and some lanthanide elements provides compelling control opportunities [23], e.g. the core electron can be utilized to match polarizabilities of ground- and Rydberg states for magic wavelength trapping [24] and isolated core excitations can be exploited to introduce coherent light shifts

∗ fhummel@atom-computing.com

† weber@itp3.uni-stuttgart.de

‡ hofferberth@iap.uni-bonn.de

§ mli@atom-computing.com

and probe Rydberg states [25, 26]. Some isotopes allow for the storage of quantum information in their nuclear spin, achieving coherence times on the order of tens of seconds [27].

Most schemes to entangle neutral atoms rely on strong dipolar interactions between Rydberg states, which are either coupled by near-resonant laser fields [28-31] or by dressing ground states with far off-resonant lasers [3235]. Yet, for the Rydberg structure, typically only one or at most a few energy levels are considered explicitly. Non-perturbative Rydberg-pair interaction calculations are required to confirm the absence of other pair states and leakage channels in the Rydberg spectrum [18, 36, 37].

In this article, we present calculations of nonperturbative Rydberg-pair interactions for divalent atoms with hyperfine-split ionization thresholds in order to study interaction tunabilty and the robustness of perturbative predictions. The complex spectra of these atoms are described within multi-channel quantum defect theory [38], which has been adapted to include nuclear spin using frame transformation [39]. Recently, 87 Sr was predicted to feature unusual properties not encountered in atoms with a simpler structure, namely, for some Rydberg states, simultaneous insensitivity to electric fields and resonant dipolar interaction in equal-state pairs [40]. We find this to be a general feature caused by the hyperfine interaction between different Rydberg series when the energy splitting of adjacent Rydberg manifolds becomes comparable to the threshold splitting in the ionic core. Specifically, the coupling of Rydberg series belong-

ing to different ionic-core hyperfine configurations gives rise to varying electron spin character of some Rydberg states as a function of the principal quantum number [41, 42] and to a corresponding energy dependence of their quantum defects. The resulting properties are appealing for quantum technology applications and Rydberg gas experiments: High degrees of level stability help preparing Rydberg states that minimize amplitude and phase errors. Intrinsic F¨rster resonances lead to strong o dipole-dipole interactions, allowing for the engineering of exotic many-body Hamiltonians and for the implementation of two- and multi-qubit gates [43]. These F¨ orster resonances do not require careful tuning of external fields, which is often the case for atoms with a single valence electron [44-48], but tuning of the principal quantum number is in general sufficient as we will demonstrate. Moreover, the F¨ orster resonances occur for two atoms being in the same state, allowing for global addressing of the atoms via a single laser beam.

This article is structured as follows: In Sec. II, we review the theory and methods, particularly the Hamiltonian for a pair of Rydberg atoms in Sec. II A. After this general discussion of Rydberg interactions, we provide the methods specific to divalent Rydberg atoms in Sec. II B, where we outline multi-channel quantum defect theory and frame transformation. Sec. II C details 87 Sr, with comments on exact diagonalization in Sec. II D and on the developed software in Sec. II E. Sec. III contains the results, discussing the single-atom spectral structure in Sec. III A and introducing field-free pair-interaction potentials for the F¨rster resonant Rydo berg series in Sec. III B. Sec. III C covers the response of single atoms and Rydberg pairs to magnetic fields, similarly for electric fields in Sec. III D. A comparison to a different, single-channel Rydberg series is provided in Sec. III E. Finally, Sec. IV contains the conclusion and an outlook.

## II. THEORY AND METHODS

## A. Rydberg-pair Hamiltonian

We consider a pair of atoms as illustrated in Fig. 1 and modelled by the Hamiltonian

$$H _ { 2 } = H _ { 1 } \otimes \mathbb { 1 } + \mathbb { 1 } \otimes H _ { 1 } + H _ { \text{int} } ( \mathbf R ), \quad \ ( 1 ) \quad \stackrel { \dots } { \mathbf T }$$

where H 1 is the single-atom Hamiltonian detailed further below and (using atomic units throughout)

$$H _ { \text{int} } ( R ) = \sum _ { \kappa _ { 1 }, \kappa _ { 2 } = 1 } ^ { \infty } V _ { \kappa _ { 1 } \kappa _ { 2 } } | R | ^ { - 1 - \kappa _ { 1 } - \kappa _ { 2 } } \quad \ ( 2 ) \quad \text{ while} \quad \text{quad } \text{quad } \text{$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$ \quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad\text{$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\ \text{$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$$$

is the multipole expansion of the electrostatic interaction between the two atoms suitable for internuclear separations exceeding the Le Roy radius and neglecting retardation effects. Here, V κ κ 1 2 is proportional to the spherical f = j r + f c

FIG. 1. The molecular system consisting of two Rydberg atoms with internuclear coordinate R and electronic coordinates r i for atom i , respectively. For each individual atom, the total hyperfine spin f is a good quantum number. In the fragmentation frame, f is obtained by coupling the ionic core's hyperfine angular momentum f c and the Rydberg electron's total angular momentum j r . The coupling scheme for short electron separations r i depends on the particular eigenstate.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00195v1-figure-c218e8a4-49f8-48c6-b532-88cf3d8cd759.jpg)

multipole operators on the Rydberg electron in atom 1 and 2, respectively. The exact form of V κ κ 1 2 depends on the choice of the coordinate system. We choose the internuclear axis R to align with the z axis and restrict Eq. 2 to dipole-dipole interactions:

$$V _ { 1 1 } = - \sum _ { q = - 1 } ^ { 1 } \binom { 2 } { 1 + q } r _ { 1 q } ^ { ( 1 ) } r _ { 1 - q } ^ { ( 2 ) }, \quad \ \ ( 3 )$$

where r ( ) i κq is the spherical multipole operator of order κ for atom i .

We find a matrix representation of H 2 by expanding in a set of appropriately symmetrized wave functions of pairs of atoms. In the absence of external fields, H 2 preserves reflection symmetry through a plane containing the internuclear axis. Since we consider homonuclear pairs, H 2 also preserves inversion symmetry. In the case of external field alignment with the internuclear axis, the total magnetic quantum number M = m 1 + m 2 is a conserved quantity, where we introduce an index, indicating m i corresponds to atom i . Finally, since only dipoledipole interaction are considered in the multipole expansion (Eq. 2), H 2 preserves permutation symmetry. We use these symmetries to decompose H 2 into blocks that can be treated independently of each other [49].

The Hamiltonian for a single Rydberg atom is given by

$$H _ { 1 } = H _ { 0 } + H _ { s } + H _ { z } + H _ { d },$$

where H 0 is the field-free Hamiltonian with the good quantum numbers f and m corresponding to the total hyperfine spin, as well as parity. The eigenvalues of H 0 ,

$$E _ { \nu f } = E _ { I } - \frac { R _ { \mu } } { \nu ^ { 2 } },$$

are conveniently expressed in terms of an effective principal quantum number ν , defined with respect to the

ionization threshold E I . Here, R µ is the mass scaled Rydberg constant. If an eigenstate is a superposition of multi-channel configurations as described in the following section, it is still useful to choose a specific E I in order to label eigenstates by a global ν and, hence, indicate their ordering. For hyperfine split-ionization thresholds, we choose E I to be the energetically higher one, which some single-channel Rydberg series converge to.

The effect of static and homogeneous external fields is considered by the following three terms: For the electric field F , the Stark Hamiltonian is given by

$$H _ { s } = - r \cdot F, \quad \ \ \ \ ( 6 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where r is the electronic coordinate (and, since the electron charge e = 1 in atomic units, corresponds to the electric dipole moment). For the magnetic field B , the Zeeman Hamiltonian is given by

$$H _ { z } = - \mu \cdot B, \quad \quad \quad \ ( 7 ) \quad \text{,}$$

where µ is the magnetic dipole moment taking into account the electronic and nuclear spins

$$\mu = - \mu _ { B } \left [ \ell _ { c } + \ell _ { r } + g _ { s } \left ( s _ { c } + s _ { r } \right ) \right ] + \mu _ { I } I, \quad ( 8 ) \quad \text{sinal}$$

where µ B and µ I are the Bohr magneton and the nuclear magnetic moment and g s is the electron spin g-factor. Here, the subscripts, c and r label core and Rydberg electron quantum numbers, while ℓ and s are orbital and electronic spin angular momentum, respectively, and I is the nuclear spin (cf. Fig. 1). The diamagnetic Hamiltonian is given by

$$H _ { \text{d} } = \frac { 1 } { 8 } \left ( r ^ { 2 } B ^ { 2 } - | r \cdot B | ^ { 2 } \right ). \quad \quad ( 9 ) \quad \text{scal} \\ \text{stat} \\ \text{to } 1.$$

This term is typically considered for ground-state atoms only when the magnetic field is on the order of 1 atomic unit ( ∼ 10 4 T). However, as it scales with the surface area of the atom, H d ∼ n 4 , for the Rydberg states considered here, diamagnetic interaction is comparable to Zeeman interaction for fields as low as ∼ 100 G.

In practise, before dealing with the interactions between two atom, the total single-atom Hamiltonian H 1 is constructed within a restricted basis (see Sec. II D) and diagonalized for a specific external field configuration; a process we call prediagonalization [49]. An appropriate two-atom basis is then constructed by pairwise, symmetrized combinations of the resulting eigenstates with the option to provide an energy cutoff. This is particularly relevant in order to only include well converged single-atom states in the two-atom basis. Finally, Rydberg pair-interaction potentials are obtained by constructing H 2 in the two-atom basis and subsequent diagonalization at the internuclear distances of interest.

## B. Multi-channel quantum defect theory

In order to find the eigensystem of H 0 , we follow the formalism presented by Robicheaux et al. combining multi-channel quantum defect theory (MQDT) and frame transformations to capture complex spin interactions in heavy atoms [39]. Generally, the perspective of quantum defect theory is to consider an atom as a collision of an electron and an ionic core [50]. In MQDT, the collision is described in terms of the real, symmetric K matrix that parameterizes the scattering between different spin configurations, so-called channels [38]. The frame transformation is the process of unitarily transforming K matrices according to different coupling schemes that depend on the relevant length scale of the Rydberg electronic coordinate [51]. Each frame has a different advantage: In the fragmentation frame, it is straight-forward to incorporate multiple ionization thresholds and account for atomic hyperfine structure (or in the case of molecules, for vibrational or rotational degrees of freedom). The close-coupling frame, on the other hand, provides a diagonal representation of the K matrix.

From a scattering perspective, the total hyperfine spin of the Rydberg atom f is formed by coupling the asymptotically separated spins of the Rydberg electron j r and the ionic core's hyperfine spin f c . However, at small separations, j r and f c are not meaningful and the fragmentation-frame coupling scheme has to be replaced by a more suitable close-coupling scheme, e.g. spin-orbit coupling. This provides the additional advantage that the the K matrix in the close-coupling frame, obtained from spectroscopic data, are applicable (to a good approximation) independently of the isotope and the corresponding nuclear spin configuration. An intuitive explanation for this is provided by comparing the energy scales of different spin couplings: Close to the groundstate, the hyperfine structure is only a small correction to the eigenstates dominated by spin-orbit coupling (or j j c r -coupling for some excited states, j c being the total angular momentum of the core electrons) and therefore has very limited impact on the coupling parameters.

The formalism presented in [39] and outlined below solves for the bound states of H 0 after transforming from the closed-coupling frame into the fragmentation frame. This is equivalent to, but leads to different equations than solving directly in the close-coupling frame, presented e.g. in [52]. Let us start in the close-coupling frame, in which we assume the K matrix to be diagonal, denoted as ¯ , with elements K

$$\bar { K } _ { \alpha } = \tan \pi \mu _ { \alpha }, \text{ \quad \ \ } ( 1 0 )$$

where µ α are (possibly energy dependent) quantum defects and α labels the close-coupling channels. We choose the transformation into the fragmentation frame to be expressed as a series of unitary matrices

$$\hat { U } _ { c c - f j } = \hat { U } _ { j j - f j } \hat { U } _ { L S - j j } \hat { U } _ { c c - L S }, \quad \ ( 1 1 )$$

where the subscripts label coupling schemes (detailed below), the first two matrices are solely determined by the corresponding angular momentum recouplings. ˆ U cc -LS accounts for corrections to pure spin-orbit ( LS ) coupling that is typically expressed as a series of simple rotations

between specific LS channels, the (possibly energy dependent) angles of which serve as an optimization parameter in order to fit experimental data 1 . We list here the representations of the coupling schemes that make for the frame transformation matrices analogously to references [39, 40]. When the Rydberg electron is close to the core, we write the wave function as

$$| L S \rangle = | ( ( ( ( s _ { c } s _ { r } ) S ( l _ { c } l _ { r } ) L ) J I ) f m \rangle, \quad ( 1 2 ) \quad \text{with}$$

where parentheses indicate which spins are coupled consecutively. This configuration is recoupled in an intermediate step to

$$| j j \rangle = | ( ( ( s _ { c } l _ { c } ) j _ { c } ( s _ { r } l _ { r } ) j _ { r } ) J I ) f m \rangle, \quad ( 1 3 )$$

before the final recoupling to the fragmentation frame

$$| f j \rangle = | ( ( ( s _ { c } l _ { c } ) j _ { c } I ) f _ { c } ( s _ { r } l _ { r } ) j _ { r } ) f m \rangle. \quad ( 1 4 ) \quad \stackrel { \dots } { \mathbb { T } }$$

This scheme can be summarized to be propagating the coupling of the valence electrons to the nuclear spin through different stages.

Finally, the K matrix in the fragmentation frame is obtained as

$$\hat { K } = \hat { U } _ { c c - f j } \bar { K } \hat { U } _ { c c - f j } ^ { \dagger }.$$

Generally, the K matrix is block diagonal in f and in parity. Thus, under the assumption ℓ c = 0, it further decomposes into blocks of equal orbital angular momentum ℓ r .

A bound state | ψ νf ⟩ at energy E νf is obtained by solving

$$\det \hat { M } = 0, \quad \quad ( 1 6 \tt a ) \quad \tt a s$$

$$\hat { M } = \tan \beta \, \hat { \mathbf i } + \hat { K }, \, \text{ \quad \ \ } ( 1 6 b )$$

where the elements of β , β i = π ν ( i -ℓ r,i ), are related to the relative principal quantum numbers ν i defined with respect to the associated ionization threshold of the i -th channel, i.e.

$$E = E _ { i } - \frac { R _ { \mu } } { \nu _ { i } ^ { 2 } }, \ \forall i, \quad \quad ( 1 7 ) \quad \sin ^ { \text{ca} } \text{size} \\ \text{ne}$$

for total energy E and thresholds E i . At a bound state energy E νf , the coefficients for each channel A i are defined via the one-dimensional null space

$$A = \ker \left ( ( \nu ^ { - 3 / 2 } \cos \beta ) ^ { \top }. \hat { M } \right ), \quad \ \ ( 1 8 \mathbf a ) \quad \$$

$$A ^ { 2 } = 1, \quad \quad ( 1 8 b ) \quad \text{in}$$

1 For light atoms, this is often not required and the K matrix is diagonal in the LS frame. As pointed out in [40], for 87 Sr, a coupling between the 5 snd series singlet and triplet states, 1 D 2 and 3 D 2 , can be accounted for by this additional rotation.

where the exponentiation of ν , containing elements ν i , as well as the matrix multiplication (.) are assumed to be element-wise operations. Then, the total MQDT wave function is given as

$$| \psi _ { \nu f } \rangle = \sum _ { i } A _ { i } ^ { ( \nu f ) } P _ { \nu _ { i } \ell _ { r, i } } ( r ) | \phi _ { i } \rangle, \quad \ \ ( 1 9 )$$

with the channel function | ϕ i ⟩ representing the ionic core as well as the spin and angular configuration of the Rydberg electron. The radial representation of the Rydberg electron is given by

$$P _ { \nu _ { i } \ell _ { r, i } } ( r ) = \left ( f _ { i } \cos \beta _ { i } + g _ { i } \sin \beta _ { i } \right ) \nu _ { i } ^ { - 3 / 2 }, \quad ( 2 0 )$$

with the regular and irregular Coulomb functions f i and g i , respectively. Eq. 16 guarantees the Rydberg wave function to go to zero asymptotically.

The total wave functions (Eq. 19) containing the radial representation only for the Rydberg electron (Eq. 20), but angular representation for both the Rydberg electron and the ionic core, are used for the calculation of all matrix elements, relevant for constructing the Hamiltonians discussed in the previous sections, particularly the electric multipole moments of order κ

$$\underset { \begin{array} { c } \text{in} \\ \text{le} \end{array} } { \left \langle \psi _ { \nu f } \right | r _ { \kappa q } \left | \psi _ { \nu ^ { \prime } f ^ { \prime } } \right \rangle } = \sum _ { i, j } A _ { i } ^ { ( \nu f ) } A _ { j } ^ { ( \nu ^ { \prime } f ^ { \prime } ) } R _ { i j } ^ { \kappa } \left \langle \phi _ { i } \right | Y _ { \kappa q } \left | \phi _ { j } \right \rangle, } \\ \text{in} \\ \text{le} \end{array} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad\text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in}, \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \text{in} \quad \end{array}$$

with the spherical harmonics Y κq and the radial integrals

$$R _ { i j } ^ { \kappa } = \int \mathrm d r \ r ^ { \kappa } P _ { \nu _ { i } \ell _ { r, i } } ( r ) P _ { \nu _ { j } ^ { \prime } \ell _ { r, j } ^ { \prime } } ( r ), \quad ( 2 2 )$$

as well as the magnetic dipole moment

$$\left \langle \psi _ { \nu f } | \, \mu \, | \psi _ { \nu ^ { \prime } f ^ { \prime } } \right \rangle = \sum _ { i, j } A _ { i } ^ { ( \nu f ) } A _ { j } ^ { ( \nu ^ { \prime } f ^ { \prime } ) } R _ { i j } ^ { 0 } \left \langle \phi _ { i } \right | \, \mu \, | \phi _ { j } \right \rangle. \quad ( 2 3 )$$

This assumes the approximation that contributions from the dipole operators acting on the core are negligible compared to contributions from the Rydberg electron because the corresponding wave function differ vastly in size. Here, the representation in the fragmentation channel is beneficial: Rydberg and core quantum numbers are well separated therein.

## C. Strontium 87

For the Rydberg structure of 87 Sr, we use the experimentally reported quantum defects for 88 Sr in the LS frame as listed in [40]. We are interested in the 5 snℓ Rydberg series close to the ionization threshold at relatively large principal quantum numbers n ⪆ 30. In Sr, some doubly excited states of the 5 pnp -, 4 dnd -, and 4 dnp -series lie energetically below the first ionization threshold and are particularly relevant in order to accurately describe the spectrum at low n ⪅ 15 [53]. However, for

our purposes, the influence of these series is sufficiently captured by considering the quantum defects of the associated 5 snℓ series as slowly varying functions of energy [40]. In the case of 87 Sr + with its fermionic core and neglecting electron excitations, i.e. ℓ c = 0, the finite nuclear spin I = 9 2 couples to the electronic angular momentum j c = s c = 1 2 giving rise to two core hyperfine configurations f c with energies separated by roughly 5 GHz. The energies of the ionic ground states with different hyperfine configurations serve as distinct ionization thresholds E f c for the Rydberg electron.

To visualize which pair-interaction potentials can possibly be excited via an intermediate level -c¿, we use a color code highlighting the potentials according to the dipole coupling of Rydberg-pair states | Ψ ⟩ to the pair state | c, t ⟩ , where | t ⟩ is the Rydberg target state to be specified and we have exemplary chosen | c ⟩ to be the 3 P 0 f = 9 2 clock state. Thus, the coloring corresponds to the matrix element D = | ⟨ Ψ | r ⊗ 1 1 | c, t ⟩ | . Note that the state | c, t ⟩ is not symmetrized, therefore the coloring is not sensitive to being in the dark- or bright-state sector of the Hamiltonian. In pure Born-Oppenheimer approximation, these sectors are well separated, however, considering motional degrees of freedom, there is some finite coupling. We consider this choice for coloration to be more instructive than coloring according to the overlap with a specified target state, because the latter neglects the impact of other dipole coupled states with for example differing hyperfine configuration or different orbital angular momentum. For e.g. blockade experiments, it is crucial to find an energy window as well as a region of internuclear distances, which is 'clean', i.e. free of dipole-coupled potentials. Unless otherwise specified, we choose for all pair-interaction potential figures in this article, D &gt; 10 -5 ea 0 for colored points and D &gt; 10 -9 ea 0 for gray, high-opacity points to inform about the encountered level densities.

## D. Remarks on basis sets

A major challenge for the exact diagonalization of Rydberg pair Hamiltonians is the choice of the basis set and the corresponding convergence of eigenstates. For van-der-Waals type interaction, the pair states relevant to converge a given pair potential can be energetically distant, defining a requirement for the considered energy window for diagonalization. In general external field configurations, when M is not a conserved quantity, the pair state level density becomes large and the number of required basis states is typically the limiting factor. Adding to this precarity, for multi-channel atoms, the single-atom level density is generally larger than e.g. in alkali atoms due to the relevance of the hyperfine interaction.

To illustrate this, let us consider the number of bound states within a single Rydberg manifold restricting the maximal orbital angular momentum to ℓ max = 3. In

87 Rb, where I = 3 2 , there are technically 32 levels with average degeneracy of 4 each. However, since the hyperfine interaction scales as n -3 , the nuclear spin can be disregarded for sufficiently large principal quantum numbers shrinking the number of levels to 7, or 32 in the case of resolved degeneracy. Contrary for 87 Sr, hyperfine interaction cannot be disregarded 2 and the number of levels is 64 with average degeneracy 10 each, due to the large nuclear spin I = 9 2 . A factor of 20 in the single atom level density has devastating consequences for the diagonalization of pair states, where the number of states scales quadratically and the number of matrix elements to construct the Hamiltonian scales with the power of 4.

In order to achieve convergence on the sub-MHz scale for the figures in this article, the typical basis sizes range from 10k to 25k states, after symmetrization. For external fields not aligned with the z axis (not discussed here), we were only able to converge potentials to the tens-of-MHz scale using around 45k basis states. Unless otherwise specified, we take into account 3 to 4 neighboring n manifolds, i.e. 7 to 9 manifolds in total, and orbital angular momenta ℓ from 0 to 3 in the single atom basis. States further separated from the target state typically do not have the coupling strength to significantly contribute, although even energetically distant states could, in principle, make a degenerate pair.

For calculations with finite electric field, we compare the results including higher ℓ up to 4 or 5, however this does not serve as a proper convergence check in the case of high electric fields, because contributions from degenerate highℓ states could, in principle, add up to modify the level structure, hence, the focus on small electric fields, where we expect those contributions to be small. Combinations of such single atom eigenstates are considered in the pair basis, if their total energy is within a specified energy window of the target state. We found 20 GHz to be sufficient for our convergence goal for the considered states with n ≈ 60, however, in arbitrary external field configurations, this may lead to intractable basis sizes and the convergence criterion may need to be loosened.

## E. Pairinteraction software

In order to account for the increased requirements associated with calculating pair-interaction potentials of multi-channel atoms discussed in the previous sections,

2 In alkali metals, the hyperfine interaction scales with the principal quantum number due to the overlap of the electron wave function with the nuclear core. While in alkaline-earths the Rydberg electron is similarly disconnected from the nucleus, the second valence electron is always close to the core probing its structure. The electron-spin coupling of the two valence electrons prevents decoupling of the hyperfine structure and the corresponding scaling.

we extended the pairinteraction software. A corresponding open-source release is currently pending. 3 Notably, the software supports sophisticated diagonalization routines such as the FEAST algorithm [54]. It is planned to offer cloud access to the pairinteraction software making use of CUDA for distributed computation on GPUs [55]. Furthermore, we propose a database standard for atomic state data. Such a database contains the bound state energies and appropriate state labeling according to the relevant quantum numbers as well as, importantly, the reduced matrix elements required for constructing Hamiltonians with high-performance [56].

Such a database separates the calculation of pairinteraction potentials from the more fundamental calculation of atomic spectra and, particularly, removes any requirement for the software to handle MQDT wave functions. Instead, we provide a Julia [57] implementation of MQDT using frame transformations following Robicheaux et al. [39] and outlined above as a utility tool for pairinteraction that will be available as open-source software. An MQDT spectrum is generally sensitive on the underlying quantum defect data and the considered channel couplings. Therefore any experimental advances in the spectroscopy of multi-channel atoms and improved efforts to fit quantum defects will, in principle, lead to new databases that can be used with the pairinteraction software.

## III. RESULTS

In the following we will show that the multi-channel couplings stemming from the hyperfine splitting of the ionic core and resulting in singlet-triplet mixed Rydberg series extends the complexity and versatility of Rydberg pair interactions. Specifically, we focus on an intrinsic equal-state F¨rster resonance in the 5 o sns f = 9 2 series of 87 Sr. We can control the corresponding Rydberg-pair interactions in four ways, namely by the choice of the principal quantum number, by the magnitude of homogeneous magnetic and electric fields, and by the internuclear distance between atoms. In contrast, the 5 sns f = 11 2 series of 87 Sr is effectively described in a single channel approximation and provides the features known from Rydberg states in alkali metals with properties regularly scaling with the principal quantum number.

## A. Overview of the F¨rster resonance o

Fig. 2 illustrates the level structure and coupling elements of 87 Sr relevant for the composition of Rydberg

3 We refer to the software release and the accompanying publication for a more complete description of the improvements and the corresponding capabilities. The software will be available at https://pairinteraction.github.io/ .

FIG. 2. Transition energies to the 5 snp manifold and reduced dipole matrix elements for the mixed singlet-triplet series 5 sns f = 9 2 . Each Rydberg manifold provides a total of ten different available states (colored points, some on top of each other), three of which are particularly strongly coupled. Rydberg scaling of energies and dipole elements is taken into account and the s state's energy is set to zero for reference. (a) Each p series has a characteristic energy dependence. (b-c) Zooms into (a), the three strongly coupled states are equidistant from the s state around ν = 60. Each combination of these series gives rise to a resonant pair with slightly different F¨ orster defect leading to complex Rydberg-pair interaction behavior.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00195v1-figure-1e33c91e-9783-4909-9354-3fdc11527d36.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00195v1-figure-e77c961e-df95-45a8-9129-3cb01d30b847.jpg)

pair-interaction potentials as functions of the effective principal quantum number ν . (For this series n ≈ ν +3 5, . where n counts the electronic orbital.) Accounting for the relevant scaling associated with these properties, we show the level separation of the 5 sns f = 9 2 series to its dipole-coupled states, i.e. states of 5 sn p f ′ ′ = { 7 2 , 9 2 , 11 2 } character. Those states determine both the polarizabilty and the character of the pair-interaction potentials. Each Rydberg manifold provides in total ten relevant states (not all of which are resolved in the figure), that can be grouped according to their total electronic character S , which can be mixed or pure. A subset of three states has very similar mixed singlet-triplet character to the 5 sns f = 9 2 series and their dipole coupling is correspondingly strong.

Importantly, the position of the states of neighboring manifolds within the spectrum relative to the s series slowly changes as a function of ν and around ν ≈ 61 the states of the neighboring manifolds are equidistant, cf. Fig. 2(b-c). In other words, the F¨rster defects, i.e. the o energy differences, of pairs ( n -1) p + np and ns + ns have zero crossings giving rise to a broad F¨rster resoo nance. Additionally, the contribution of these transitions to the Stark effect cancel giving rise to near zero polarizabilty. Both of these effects have been predicted by perturbative considerations [40]. Fig. 2 reveals the complex interplay of the contributing states. The more weakly coupled states as seen in Fig. 2(a), which have strong triplet character, contribute to the overall structure and illustrate the significant energy dependence of quantum defects, which originates in the hyperfine-split ionization

threshold of 87 Sr + .

In the following, we analyze the conditions for this F¨ orster resonance, indicating that it is intrinsic to an equal-state pair in the s series strongly coupling to the adjacent p levels at the zero crossing of the associated F¨ orster defect, which is unique to specific Rydberg series in divalent atoms. In Rydberg series with slowly varying energy dependence of the quantum defects, for example the 5 sns f = 11 2 series in 87 Sr, the functional behavior of the F¨rster defect features the regular o n -3 scaling. This scaling and the consequential absence of a zero crossing of the F¨rster defect prohibits the existence of a F¨rster reso o onance. Particularly, this scaling is generally present in alkali atoms, where F¨rster resonances are not typically o found in equal-state pairs and rarely without additional tuning of energy levels with external fields [44]. Instead, the existence of a F¨rster resonance requires the analyo sis of many spectral properties and the careful tuning of the quantum numbers to different values for the involved atoms in order to achieve degeneracy of a pair [47]. Examples are p series in which the fine structure sufficiently shifts the state to be resonant with adjacent s states or an nsn d ′ pair or inter-species resonances using different principal quantum numbers for each species [58].

In hyperfine-split ionic-core species, however, the existence of a mixed single-triplet s series featuring a F¨rster o resonance is guaranteed, because the hyperfine interaction leads to an exchange in electronic character as a function of n , when the intrinsic splitting of adjacent Rydberg states is comparable to the ionic hyperfine splitting. For example, we expect such a resonance for the 6 sns f = 1 2 series of 171 Yb, however, the uncertainty of quantum defects of the 5 snp series makes a quantitative prediction of its location in the spectrum currently dispensable . 4

The unique combination of the features makes the F¨ orster resonant state an interesting candidate for quantum technology applications: The bare Rydberg line is largely insusceptible to external fields and their fluctuations. Beyond the small polarizability, the mixed electron-spin character of the state also leads to a suppressed linear Zeeman splitting, because singlet states have a vanishing g factor, and the response to magnetic field is mostly quadratic driven by diamagnetic interaction requiring large fields. At the same time, the pair state is subject to dipole-dipole interactions that provide a larger blockade radius than typically encountered for pair states subject to induced-dipole (i.e. van-der-Waals) interactions. And finally, the equal-state pair can be easily addressed in a global beam setup requiring only a

4 On the other hand, the quantum defects for the s and d series in Yb are rather well known. It is, however, inherently difficult to probe the p series in heavy divalent atoms, in which Rydberg spectroscopy is typically performed using two-photon transitions, because the single photon transition is too deeply ultra violet. We hope to see new Yb spectroscopy data soon or that this article can serve as a motivation for additional experimental efforts.

single wavelength as opposed to F¨rster states involving o pairs of states differing in n or ℓ or even interspecies resonances as recently demonstrated [59].

## B. Field-free pair-interaction potentials

In Fig. 3 we show the Rydberg-pair interaction potentials for equal pairs in the 5 sns f = 9 2 series in the vicinity of the F¨rster resonance at o n = 64. Here, the asymptotic energy of the 5 sns +5 sns pair is set to zero and the coloring indicates dipolar coupling to the 3 P 0 f = 9 2 clock / state. Over a relatively short range of principal quantum numbers, n = { 59 62 64 66 69 , , , , } shown in subfigures (a) through (e), the potentials drastically change character. At low n far from resonance, the potential is repulsive and scales as R -6 with the internuclear distance, indicating induced-dipolar interactions, and the pair interaction is well described using perturbation theory. Closer to resonance, the scaling changes to R -3 , indicating permanentdipole interactions due to the mixing of pair states with s and p character. At n = 64, the s and p pairs are asymptotically degenerate and the molecular eigenstates are symmetric superpositions, which leads to symmetrically attractive and repulsive potentials as is characteristic for a F¨rster resonance. o However, the fact that several pairs of p character contribute, leads to some nongeneric level crossings and state mixings at large internuclear distances and at small detunings. Furthermore, the significant p state character of the potentials gives rise to yet more additional couplings to other pair states of either s or d character at larger energies (e.g. around E/h = 150 MHz and E/h = 250 MHz) that lead to avoided crossings. At larger n , the dominant potential remains attractive and the scaling transitions to dominantly R -6 again. At short internuclear distances R &lt; 4 µ m, the potential is modified due to significant dipolar interactions stemming from other pair states at smaller energies. The larger n , the more likely such modifications become, due to the increasing level density and the increasing dipolar interaction strength.

## C. Magnetic field response

The multi-channel character of the different Rydberg series in 87 Sr gives rise to diverse external field responses depending on the principal and orbital quantum numbers. High level densities and the requirement for magnetic field strengths sufficient to lift the degeneracy in the ground ( 1 S 0 ) or clock ( 3 P 0 ) manifold, which is often the case for quantum technology applications, can lead to non-generic level shifts and avoided crossings in the Rydberg spectra of single atoms and pairs.

Fig. 4 shows magnetic field maps for the F¨rster reso onant pair and surrounding states (a) asymptotically in the non-interacting limit and (b) at 4 µ m internuclear spacing. Here, the magnetic field is aligned with the in-

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00195v1-figure-d65126ba-c6bb-48cb-89d3-f38144308f09.jpg)

FIG. 3. Pair-interaction potentials for the 5 sns f = 9 2 equal-state pairs of 87 Sr for different n = { 59 62 64 66 69 , , , , } (a-e). Varying n , the potentials change character from repulsive to attractive, while on resonance (c), the potentials are symmetric and asymptotically degenerate. However, many other pair states contribute to the potential energy landscape introducing avoided crossings and possible channels for population transfer. The color indicates the dipole transition strength from the 3 P 0 clock state. Only the M = -9 symmetry sector is shown.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00195v1-figure-c1f7adbe-7080-4386-9c37-7371cb944b9e.jpg)

Multiple pair states with p character, which are close in energy at zero field, are shifted to smaller energies due to their strong linear response, i.e. Zeeman effect. Note that most of the states shown in Fig. 4 are field seeking because of M being negative. Only in the case of f &lt; I , such states are shifted to larger energies. However, some pairs of p states, while initially being shifted out of resonance for small fields, become resonant with the s pair again at B &gt; 150 G. This is striking and highlights the importance of diamagnetism: If only linear shifts were considered, all p pairs would be shifted out of resonance, which would drastically change the character of the pair interaction potentials.

FIG. 4. Magnetic field maps for Rydberg pairs in the vicinity of the F¨rster resonant pair. o (a) Asymptotically, in the noninteracting limit, pair states are subject to vastly differing field responses all subject to diamagnetic interaction. Black indicates pairs of s states, blue pairs of p states, and gray additional pairs involving s , p , and d states. (b) At 4 µ m internuclear distance and strong dipolar coupling, the varying F¨ orster defects cause avoided crossings of the potentials. The colorcode in (b) indicates coupling strength to the clock state.

ternuclear axis such that the total projection of the hyperfine state is still a conserved quantity and only the M = -9 symmetry sector is shown. For the 5 s 64 s equalstate pair, highlighted in black in Fig. 4(a), the linear magnetic field response is weak, due to a strong contribution from the singlet channel in that state, and the quadratic response due to diamagnetism dominates even at relatively small field strength. For such states that are field-seeking, i.e. they lower their energy with respect to the magnetic field strength due to Zeeman interaction, there exists a local minimum in the magnetic field map, where the pair state is insensitive to magnetic field fluctuations due to a balance between attractive Zeeman and repulsive diamagnetic interaction. For the state of interest this occurs at approximately B = 50 G.

At short internuclear distances in Fig. 4(b), dipoledipole interactions are additionally relevant, and the diverse magnetic field responses of the pair states leads to non-trivial functional behavior of interaction potentials, with a multitude of avoided crossings, and of their dipole addressability. In particular, at B = 200 G, one of the p pairs initially resonant at zero field, becomes degenerate with a different pair state. This pair, which is detuned by approximately -200 MHz at zero field (cf. subfigure (a)), is from the same 5 sns f = 9 2 Rydberg series, however, for different n , namely n ′ = n ± 2, and one of the two states has slightly different singlet-triplet configuration due to its multi-channel nature. In a magnetic field, this pair remains approximately equidistant to the F¨rster resoo nant state. However, at B = 200 G, the two s pairs become strongly coupled: Essentially, two F¨rster resoo nances occur simultaneously at approximately 200 MHz energy splitting and the involved p states mediate coupling between the s pairs.

The corresponding pair potentials are shown in Fig. 5 for B = 50 G (a) and B = 200 G (b). Dipolar interactions exceed the energy scale of the asymptotic splitting and the pair interaction potentials feature a broad avoided crossing (cf. Fig. 5(b)) giving rise to an essentially flat dipole-coupled potential over the full range of

FIG. 5. Competition of F¨rster resonances between pairs o 64 s + 64 , 63 s p + 64 , and 62 p s + 66 . s Their F¨rster defects o depend on the magnetic field strength. (a) At B = 50 G, the 64 s +64 s and several 63 p +64 p pairs are close to resonance, while 62 s +66 s is detuned by ∼ 200 MHz. (b) At B = 200 G, both s pairs are resonant with different p pairs while their respective splitting is approximately maintained. However, the p pairs mitigate a strong coupling between the s pairs, hence the large avoided crossing and dipole addressability.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00195v1-figure-b6746f3b-78a0-4541-90ab-52752daa0102.jpg)

R . While for blockade experiments, this situation needs to be avoided, such a potential is highly unusual and may lead to facilitation, where Rydberg pairs are addressed resonantly.

This example illustrates that the magnetic field can effectively act as a knob for competing F¨rster resonances. o The potentials in Fig. 5 can be qualitatively reproduced by a toy model that only includes the four properly symmetrized pair states that contribute dominantly:

$$\{ | n s n s \rangle \,, | ( n - 1 ) p n p \rangle \,, | n ^ { \prime } s n ^ { \prime \prime } s \rangle \,, | ( n - 1 ) p ^ { \prime } n p ^ { \prime } \rangle \}, \ \ ( 2 4 ) \quad t h e \overset { n } { f }$$

where only s and p pairs are coupled, respectively, by dipolar interaction (of similar strength) and two different F¨rster defects control the asymptotic detunings as o functions of the magnetic field.

## D. Non-perturbative polarizability calculations

After considering the effects of external magnetic fields in the previous section, we now turn our attention towards electric fields. Generally, electric fields pose a computational challenge because contributions from highorbital-angular-momentum states become relevant even for moderate field strengths. Additionally, the influence of electric fields steeply scales with the principal quantum number ( ∼ n 5 ). Therefore, here we restrict our analysis to electric field strengths of few hundreds of mV/cm. We find that even these small electric fields can cause responses not expected from perturbation theory.

This is illustrated in Fig. 6 for R = 6 µ m internuclear spacing and electric fields aligned with the internuclear axis. While pairs of the pure triplet state 5 s 64 s f = 11 2 show the expected linear response of the derivative of

FIG. 6. Electric field response of 5 s 64 s equal-state pairs in different Rydberg series at internuclear spacing R = 6 µ m. (a) The single channel f = 11 2 state follows the expectation from perturbation theory with linear behavior of the electric field derivative. (b) For the multi-channel f = 9 2 state, hyperpolarizability is relevant even at small electric fields, giving rise to configurations with zero polarizability (here: zero local slope). At resonance, the polarizability of the f = 9 2 pairs is larger than for f = 11 2 pairs because of significant p state admixture.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00195v1-figure-ad7cd682-42ae-4664-a36b-11cfc44c960c.jpg)

its eigenenergy with respect to the electric field strength d E /d F as a function of the electric field strength F , the behavior of the F¨rster resonant pairs 5 o s 64 s f = 9 2 indicates the presence of a hyperpolarizability that becomes relevant for relatively small field strengths F ∼ 150 mV/cm. In fact, this can be exploited to achieve extremely high degrees of line stability, when the electric field is tuned to a local extremum of d E /d F . Overall, the f = 9 2 pairs show a much larger polarizability than the f = 11 2 pairs. This can be attributed to the strong p -state admixture at resonance. Asymptotically, the pure 5 64 s s f = 9 2 has near zero polarizability, approximately 2 orders of magnitude smaller than the 5 64 s s f = 11 2 state.

## E. Single channel Rydberg series

As mentioned previously, the 5 sns f = 11 2 Rydberg series in 87 Sr provides, in a lot of aspects, less complexity than the f = 9 2 series. It is a pure triplet series, well described in a single channel approximation and, hence, shows a lot of the regular scaling behaviors with respect to the principal quantum number expected from alkali metals or other single-channel Rydberg series. In this series, an equal-state F¨rster resonance is not expected o and the polarizability scales as n 7 , while induced-dipolar interaction scales as n 11 .

Fig. 7 shows the pair interaction potentials (a) at zero magnetic field and (b) at B = 200 G aligned with the internuclear axis. Here, the M = -11 symmetry sector is shown and the coupling to the clock state is achieved by σ + driving. Essentially, the magnetic field does not have

FIG. 7. Pair potentials for the 5 s 64 s f = 11 2 equal-state pairs for finite magnetic fields B = 0 G (a) and B = 200 G (b). These potentials are spherically symmetric and the characterizing C6 coefficient is only about 5% smaller in (b) compared to (a).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00195v1-figure-12d544c8-363e-403b-a5df-0bcc95ca22cc.jpg)

a strong impact on the potential. At B = 200 G the C6 coefficient is reduced by approximately 5% and the spherically symmetric structure of the potential is maintained (not shown in the figure). However, to achieve, e.g., Rydberg blockade for pairs within that series, still care must be taken for the choice of the principal quantum number in order to avoid accidental resonances.

## IV. CONCLUSION AND OUTLOOK

We have analysed the structure of Rydberg pairinteraction potentials of 87 Sr, a multi-channel, divalent atomic species. The ionic core's hyperfine structure gives rise to intriguing features within the Rydberg series, particularly states with vanishing polarizability and resonant dipolar interaction. While these properties make those states very promising candidates for application in quantum technology, the arising complexity requires careful consideration of external field configurations. We have shown that magnetic fields can undermine Rydberg blockade because of competing F¨rster resonances that o

- [1] J. D. Pritchard, D. Maxwell, A. Gauguet, K. J. Weatherill, M. P. A. Jones, and C. S. Adams, Cooperative atomlight interaction in a blockaded Rydberg ensemble, Phys. Rev. Lett. 105 , 193603 (2010).
- [2] A. V. Gorshkov, J. Otterbach, M. Fleischhauer, T. Pohl, and M. D. Lukin, Photon-photon interactions via Rydberg blockade, Phys. Rev. Lett. 107 , 133602 (2011).
- [3] O. Firstenberg, C. S. Adams, and S. Hofferberth, Nonlinear quantum optics mediated by Rydberg interactions, J. Phys. B: At. Mol. Opt. Phys. 49 , 152003 (2016).
- [4] M. Saffman, T. G. Walker, and K. Mølmer, Quantum information with Rydberg atoms, Rev. Mod. Phys. 82 , 2313 (2010).
- [5] C. S. Adams, J. D. Pritchard, and J. P. Shaffer, Rydberg

give rise to avoided crossings in the potential energy landscape and to unanticipated potential shapes, such as an almost completely flat potential with respect to the internuclear distance. On the other hand, configurations can be found that offer insensitivity to electric fields with pair states that have zero polarizability.

The results from pair-interaction calculations are readily available as input for error models and fidelity estimations in the context of entangling gates in quantum computation [60]. Specifically, the results provide quantitative information about leakage into doubly excited states. Depending on the available channels, this undesired behavior can potentially be mitigated by appropriate pulse shaping.

We have integrated the tools to analyze the spectra of multi-channel divalent atoms based on quantum defect measurements into the pairinteraction software. Due to the increased level density of fermionic species, the calculation of pair potentials in arbitrary external field configurations remains challenging, particularly for 87 Sr due to its large nuclear spin.

Recently, 171 Yb has emerged as an interesting candidate for application in quantum computation. With a nuclear spin of I = 1 2 it is a natural choice for nuclear qubit encoding. Due to a similar multi-channel structure, we expect the Rydberg spectrum of 171 Yb to share a lot of features with 87 Sr. Quantitative predictions have been enabled by related work including spectroscopy of the 1 P 1 and 3 P J series in 171 Yb [61], which we became aware of during the compilation of this manuscript.

## ACKNOWLEDGMENTS

S.W. and J.M. acknowledge funding by the Federal Ministry of Education and Research (BMBF) under the Grants QRydDemo, MUNIQC-Atoms, and by the Horizon Europe program HORIZON-CL4-2021DIGITALEMERGING-01-30 via the project 101070144 (EuRyQa).

atom quantum technologies, J. Phys. B: At. Mol. Opt. Phys. 53 , 012002 (2019).

- [6] M. Morgado and S. Whitlock, Quantum simulation and computing with Rydberg-interacting qubits, AVS Quantum Science 3 , 023501 (2021).
- [7] M. Endres, H. Bernien, A. Keesling, H. Levine, E. R. Anschuetz, A. Krajenbrink, C. Senko, V. Vuletic, M. Greiner, and M. D. Lukin, Atom-by-atom assembly of defect-free one-dimensional cold atom arrays, Science 354 , 1024 (2016).
- [8] H. J. Manetsch, G. Nomura, E. Bataille, K. H. Leung, X. Lv, and M. Endres, A tweezer array with 6100 highly coherent atomic qubits (2024), arXiv:2403.12021 [quantph].

- [9] D. Barredo, S. de L´s´leuc, V. Lienhard, T. Lahaye, and e e A. Browaeys, An atom-by-atom assembler of defect-free arbitrary two-dimensional atomic arrays, Science 354 , 1021 (2016).
- [10] H. Kim, W. Lee, H.-g. Lee, H. Jo, Y. Song, and J. Ahn, In situ single-atom array synthesis using dynamic holographic optical tweezers, Nature Commun. 7 , 13317 (2016).
- [11] H. Bernien, S. Schwartz, A. Keesling, H. Levine, A. Omran, H. Pichler, S. Choi, A. S. Zibrov, M. Endres, M. Greiner, et al. , Probing many-body dynamics on a 51-atom quantum simulator, Nature 551 , 579 (2017).
- [12] S. Ebadi, T. T. Wang, H. Levine, A. Keesling, G. Semeghini, A. Omran, D. Bluvstein, R. Samajdar, H. Pichler, W. W. Ho, et al. , Quantum phases of matter on a 256-atom programmable quantum simulator, Nature 595 , 227 (2021).
- [13] S. Ma, G. Liu, P. Peng, B. Zhang, S. Jandura, J. Claes, A. P. Burgers, G. Pupillo, S. Puri, and J. D. Thompson, High-fidelity gates and mid-circuit erasure conversion in an atomic qubit, Nature 622 , 279 (2023).
- [14] M. A. Norcia, W. B. Cairncross, K. Barnes, P. Battaglino, A. Brown, M. O. Brown, K. Cassella, C.-A. Chen, R. Coxe, D. Crow, et al. , Midcircuit qubit measurement and rearrangement in a 171 Yb atomic array, Phys. Rev. X 13 , 041034 (2023).
- [15] R. Tao, M. Ammenwerth, F. Gyger, I. Bloch, and J. Zeiher, High-fidelity detection of large-scale atom arrays in an optical lattice, Phys. Rev. Lett. 133 , 013401 (2024).
- [16] M. A. Norcia, H. Kim, W. B. Cairncross, M. Stone, A. Ryou, M. Jaffe, M. O. Brown, K. Barnes, P. Battaglino, T. C. Bohdanowicz, et al. , Iterative assembly of 171 Yb atom arrays with cavity-enhanced optical lattices, PRX Quantum 5 , 030316 (2024).
- [17] D. Bluvstein, H. Levine, G. Semeghini, T. T. Wang, S. Ebadi, M. Kalinowski, A. Keesling, N. Maskara, H. Pichler, M. Greiner, et al. , A quantum processor based on coherent transport of entangled atom arrays, Nature 604 , 451 (2022).
- [18] A. Browaeys, D. Barredo, and T. Lahaye, Experimental investigations of dipole-dipole interactions between a few Rydberg atoms, J. Phys. B: At. Mol. Opt. Phys. 49 , 152001 (2016).
- [19] A. Jenkins, J. W. Lis, A. Senoo, W. F. McGrew, and A. M. Kaufman, Ytterbium nuclear-spin qubits in an optical tweezer array, Phys. Rev. X 12 , 021027 (2022).
- [20] N. Chen, L. Li, W. Huie, M. Zhao, I. Vetter, C. H. Greene, and J. P. Covey, Analyzing the Rydberg-based optical-metastable-ground architecture for 171 Yb nuclear spins, Phys. Rev. A 105 , 052438 (2022).
- [21] G. Unnikrishnan, P. Ilzh¨fer, o A. Scholz, C. H¨ olzl, A. G¨ otzelmann, R. K. Gupta, J. Zhao, J. Krauter, S. Weber, N. Makki, H. P. B¨chler, T. Pfau, and F. Meinert, u Coherent control of the fine-structure qubit in a single alkaline-earth atom, Phys. Rev. Lett. 132 , 150606 (2024).
- [22] A. Cao, W. J. Eckner, T. L. Yelin, A. W. Young, S. Jandura, L. Yan, K. Kim, G. Pupillo, J. Ye, N. D. Oppong, and A. M. Kaufman, Multi-qubit gates and 'Schr¨dinger o cat' states in an optical clock (2024), arXiv:2402.16289 [quant-ph].
- [23] F. B. Dunning, T. C. Killian, S. Yoshida, and J. Burgd¨ orfer, Recent advances in Rydberg physics using alkaline-earth atoms, J. Phys. B: At. Mol. Opt. Phys. 49 , 112003 (2016).
- [24] J. T. Wilson, S. Saskin, Y. Meng, S. Ma, R. Dilip, A. P. Burgers, and J. D. Thompson, Trapping alkaline earth Rydberg atoms optical tweezer arrays, Phys. Rev. Lett. 128 , 033201 (2022).
- [25] A. P. Burgers, S. Ma, S. Saskin, J. Wilson, M. A. Alarc´n, o C. H. Greene, and J. D. Thompson, Controlling Rydberg excitations using ion-core transitions in alkaline-earth atom-tweezer arrays, PRX Quantum 3 , 020326 (2022).
- [26] K.-L. Pham, T. F. Gallagher, P. Pillet, S. Lepoutre, and P. Cheinet, Coherent light shift on alkaline-earth Rydberg atoms from Isolated Core Excitation without autoionization, PRX Quantum 3 , 020327 (2022).
- [27] K. Barnes, P. Battaglino, B. J. Bloom, K. Cassella, R. Coxe, N. Crisosto, J. P. King, S. S. Kondov, K. Kotru, S. C. Larsen, et al. , Assembly and coherent control of a register of nuclear spin qubits, Nature Commun. 13 , 2779 (2022).
- [28] D. Jaksch, J. I. Cirac, P. Zoller, S. L. Rolston, R. Cˆt´, o e and M. D. Lukin, Fast quantum gates for neutral atoms, Phys. Rev. Lett. 85 , 2208 (2000).
- [29] M. D. Lukin, M. Fleischhauer, R. Cote, L. M. Duan, D. Jaksch, J. I. Cirac, and P. Zoller, Dipole blockade and quantum information processing in mesoscopic atomic ensembles, Phys. Rev. Lett. 87 , 037901 (2001).
- [30] H. Levine, A. Keesling, G. Semeghini, A. Omran, T. T. Wang, S. Ebadi, H. Bernien, M. Greiner, V. Vuleti´, c H. Pichler, and M. D. Lukin, Parallel implementation of high-fidelity multiqubit gates with neutral atoms, Phys. Rev. Lett. 123 , 170503 (2019).
- [31] S. Jandura and G. Pupillo, Time-Optimal Twoand Three-Qubit Gates for Rydberg Atoms, Quantum 6 , 712 (2022).
- [32] I. Bouchoule and K. Mølmer, Spin squeezing of atoms by the dipole interaction in virtually excited Rydberg states, Phys. Rev. A 65 , 041803 (2002).
- [33] J. E. Johnson and S. L. Rolston, Interactions between rydberg-dressed atoms, Phys. Rev. A 82 , 033412 (2010).
- [34] G. Pupillo, A. Micheli, M. Boninsegni, I. Lesanovsky, and P. Zoller, Strongly correlated gases of rydberg-dressed atoms: Quantum and classical dynamics, Phys. Rev. Lett. 104 , 223002 (2010).
- [35] A. Mitra, M. J. Martin, G. W. Biedermann, A. M. Marino, P. M. Poggi, and I. H. Deutsch, Robust MølmerSørensen gate for neutral atoms using rapid adiabatic Rydberg dressing, Phys. Rev. A 101 , 030301 (2020).
- [36] S. de L´s´leuc, S. Weber, V. Lienhard, D. Barredo, H. P. e e B¨ uchler, T. Lahaye, and A. Browaeys, Accurate mapping of multilevel Rydberg atoms on interacting spin-1 / 2 particles for the quantum simulation of Ising models, Phys. Rev. Lett. 120 , 113602 (2018).
- [37] S. Weber, S. de L´s´leuc, V. Lienhard, D. Barredo, T. Lae e haye, A. Browaeys, and H. P. B¨chler, Topologically prou tected edge states in small Rydberg systems, Quantum Sci. Technol. 3 , 044001 (2018).
- [38] M. Aymar, C. H. Greene, and E. Luc-Koenig, Multichannel Rydberg spectroscopy of complex atoms, Rev. Mod. Phys. 68 , 1015 (1996).
- [39] F. Robicheaux, D. W. Booth, and M. Saffman, Theory of long-range interactions for Rydberg states attached to hyperfine-split cores, Phys. Rev. A 97 , 022508 (2018).
- [40] F. Robicheaux, Calculations of long range interactions for 87Sr Rydberg states, J. Phys. B: At. Mol. Opt. Phys. 52 , 244001 (2019).
- [41] P. Liao, R. Freeman, R. Panock, and L. Humphrey,

Hyperfine-induced singlet-triplet mixing in 3He, Optics Commun. 34 , 195 (1980).

- [42] R. Beigang, E. Matthias, and A. Timmermann, Influence of singlet-triplet mixing on the hyperfine structure of 5 snd Rydberg states in 87 Sr, Phys. Rev. Lett. 47 , 326 (1981).
- [43] D. Petrosyan, F. Motzoi, M. Saffman, and K. Mølmer, High-fidelity Rydberg quantum gate via a two-atom dark state, Phys. Rev. A 96 , 042306 (2017).
- [44] T. G. Walker and M. Saffman, Consequences of Zeeman degeneracy for the van der Waals blockade between Rydberg atoms, Phys. Rev. A 77 , 032723 (2008).
- [45] I. I. Ryabtsev, D. B. Tretyakov, I. I. Beterov, and V. M. Entin, Observation of the Stark-tuned F¨rster resonance o between two Rydberg atoms, Phys. Rev. Lett. 104 , 073003 (2010).
- [46] D. B. Tretyakov, V. M. Entin, E. A. Yakshina, I. I. Beterov, C. Andreeva, and I. I. Ryabtsev, Controlling the interactions of a few cold Rb Rydberg atoms by radiofrequency-assisted F¨rster resonances, Phys. Rev. A o 90 , 041403 (2014).
- [47] I. I. Beterov and M. Saffman, Rydberg blockade, F¨rster o resonances, and quantum state measurements with different atomic species, Phys. Rev. A 92 , 042710 (2015).
- [48] A. Paris-Mandoki, H. Gorniaczyk, C. Tresp, I. Mirgorodskiy, and S. Hofferberth, Tailoring Rydberg interactions via F¨rster resonances: state combinations, hopping and o angular dependence, J. Phys. B: At. Mol. Opt. Phys. 49 , 164001 (2016).
- [49] S. Weber, C. Tresp, H. Menke, A. Urvoy, O. Firstenberg, H. P. B¨chler, and S. Hofferberth, Tutorial: u Calculation of Rydberg interaction potentials, J. Phys. B: At. Mol. Opt. Phys. 50 , 133001 (2017).
- [50] M. J. Seaton, Quantum defect theory I. General formu-
- lation, Proc. Phys. Soc. 88 , 801 (1966).
- [51] U. Fano, Quantum defect theory of l uncoupling in H 2 as an example of channel-interaction treatment, Phys. Rev. A 2 , 353 (1970).
- [52] O. Robaux and M. Aymar, A program for analysing the Rydberg series of highly excited discrete spectra by M.Q.D.T., Comp. Phys. Comm. 25 , 223 (1982).
- [53] C. L. Vaillant, M. P. A. Jones, and R. M. Potvliege, Multichannel quantum defect theory of strontium bound Rydberg states, J. Phys. B: At. Mol. Opt. Phys. 47 , 155001 (2014).
- [54] E. Polizzi, Density-matrix-based algorithm for solving eigenvalue problems, Phys. Rev. B 79 , 115112 (2009).
- [55] Nvidia, CUDA toolkit (2024).
- [56] D. Foundation, DuckDB (2024).
- [57] J. Bezanson, A. Edelman, S. Karpinski, and V. B. Shah, Julia: A fresh approach to numerical computing, SIAM Review 59 , 65 (2017).
- [58] P. M. Ireland, D. M. Walker, and J. D. Pritchard, Interspecies F¨ orster resonances for Rb-Cs Rydberg d -states for enhanced multi-qubit gate fidelities, Phys. Rev. Res. 6 , 013293 (2024).
- [59] S. Anand, C. E. Bradley, R. White, V. Ramesh, K. Singh, and H. Bernien, A dual-species Rydberg array (2024), arXiv:2401.10325 [quant-ph].
- [60] A. Pagano, S. Weber, D. Jaschke, T. Pfau, F. Meinert, S. Montangero, and H. P. B¨chler, Error budgeting for a u controlled-phase gate with strontium-88 Rydberg atoms, Phys. Rev. Res. 4 , 033019 (2022).
- [61] M. Peper, Y. Li, D. Y. Knapp, M. Bileska, S. Ma, G. Liu, P. Peng, B. Zhang, S. P. Horvath, A. P. Burgers, and J. D. Thompson, Spectroscopy and modeling of 171 Yb Rydberg states for high-fidelity two-qubit gates (2024), arXiv:2406.01482 [physics.atom-ph].
Frederic Hummel Sebastian Weber Johannes Moegerle Henri Menke Jonathan King Benjamin Bloom Sebastian Hofferberth Ming Li
Engineering Rydberg-pair interactions in divalent atoms with hyperfine-split ionization thresholds
Quantum information processing with neutral atoms relies on Rydberg
excitation for entanglement generation. While the use of heavy divalent or
open-shell elements, such as strontium or ytterbium, has benefits due to their
optically active core and a variety of possible qubit encodings, their Rydberg
structure is generally complex. For some isotopes in particular, hyperfine
interactions are relevant even for highly excited electronic states. We employ
multi-channel quantum defect theory to infer the Rydberg structure of isotopes
with non-zero nuclear spin and perform non-perturbative Rydberg-pair
interaction calculations. We find that due to the high level density and
sensitivities to external fields, experimental parameters must be precisely
controlled. Specifically in ${}^{87}$Sr, we study an intrinsic F\"orster
resonance, unique to divalent atoms with hyperfine-split thresholds, which
simultaneously provides line stability with respect to external field
fluctuations and enhanced long-range interactions. Additionally, we provide
parameters for pair states that can be effectively described by single-channel
Rydberg series. The explored pair states provide exciting opportunities for
applications in the blockade regime as well as for more exotic long-range
interactions such as largely flat, distance-independent potentials.