Ultracold, radiative charge transfer in hybrid Yb ion - Rb atom traps

Ultracold hybrid ion-atom traps offer the possibility of microscopic manipulation of quantum coherences in the gas using the ion as a probe. However, inelastic processes, particularly charge transfer can be a significant process of ion loss and has been measured experimentally for the Yb$^{+}$ ion immersed in a Rb vapour. We use first-principles quantum chemistry codes to obtain the potential energy curves and dipole moments for the lowest-lying energy states of this complex. Calculations for the radiative decay processes cross sections and rate coefficients are presented for the total decay processes. Comparing the semi-classical Langevin approximation with the quantum approach, we find it provides a very good estimate of the background at higher energies. The results demonstrate that radiative decay mechanisms are important over the energy and temperature region considered. In fact, the Langevin process of ion-atom collisions dominates cold ion-atom collisions. For spin dependent processes \cite{kohl13} the anisotropic magnetic dipole-dipole interaction and the second-order spin-orbit coupling can play important roles, inducing couplingbetween the spin and the orbital motion. They measured the spin-relaxing collision rate to be approximately 5 orders of magnitude higher than the charge-exchange collision rate \cite{kohl13}. Regarding the measured radiative charge transfer collision rate, we find that our calculation is in very good agreement with experiment and with previous calculations. Nonetheless, we find no broad resonances features that might underly a strong isotope effect. In conclusion, we find, in agreement with previous theory that the isotope anomaly observed in experiment remains an open question.

Yb ion immersed in a Rb vapour. We use first-principles quantum chemistry codes to obtain the potential energy curves and dipole moments for the lowest-lying energy states of this complex. Calculations for the radiative decay processes cross sections and rate coefficients are presented for the total decay processes; ν + → + +

Introduction
Charge transfer processes in ion−atom collisions are traditionally measured experimentally by their cross-sections and rate coefficient as a function of energy and temperature. For ambient temperatures, one can treat the relative motion of the ion and atom as a classical motion and focus on the quantum dynamics of the electrons. However, at ultracold temperatures the wave nature of the atomic motion is revealed. While the electronic and nuclear motion can still be adiabatically decoupled, the electronic and nuclear motions are strongly correlated. Under these conditions, the chemical pathways and scattering processes are highly sensitive to external fields that perturb the electronic structure which in turn transfer this effect coherently to the atomic motion. This means that the reaction processes are sensitive to external electric and magnetic fields and thus amenable to experimental control, for example by Feshbach resonances [2,3]. Such effects are extremely important in controlling coherence and correlation, with applications in molecular quantum information protocols and in hybrid quantum systems such as Coulomb crystals (ion arrays) embedded in a quantum degenerate gas [4][5][6].
In this paper, we are concerned with one aspect of ultracold ion-atom physics: the process of ion loss by charge transfer. This is critical in terms of the ultracold regime as to whether cooling, trapping and degeneracy can be achieved. It is extremely important in view of potential applications such as sympathetic ion cooling and micromotion minimization [7]. Furthermore, it allows the study of the fundamental process of ultracold charge transfer [8][9][10]. Quantum phenomena can dominate reaction dynamics at low temperatures. In such cold conditions the scattering process becomes sensitive to the isotopes [11] for example when the resonances are sharpened by tunnelling into long-lived metastable scattering states.
Interest has developed in expanding the range of quantum systems that can be trapped and manipulated on the quantum scale. Hybrid ion-atom systems are of great interest [12,13] since these are inherently strongly-interacting systems with a longer-range potential, and inelastic processes can be studied. Recently these systems have been explored considering two-body collisions, in which both collision partners are translationally cold [14], and on the many-body level [15], where the sympathetic cooling of the ion with ultracold atoms was observed. The study of these systems in the quantum regime can be applied to hybrid ion-atom devices [16] and, in addressing fundamental many-body effects of ionic impurities such as mesoscopic molecule formation [17] and density fluctuations [18]. These devices offer a unique opportunity to study reactive collisions (ultracold chemistry) [19] under controlled conditions, for example when external electric fields can be applied to modify the reaction rates/cross [15]. Unlike binary cold collisions between ground state neutral atoms, which are only elastic or inelastic in nature, reactive collisions (charge transfer) are a feature of Yb ions immersed in a gas of trapped alkali atoms. Consequently there has been increased interest in ultracold Yb-ion chemistry in the interactions with alkalis [20] and its resonant charge transfer process [21]. While non-adiabatic effects are strongly suppressed as the temperature falls towards zero, nonetheless the product Yb + + Rb is the thermodynamically favoured species [22] and thus the loss process can occur by spontaneous emission. It is this process that has received experimental and theoretical attention recently and we continue to study in detail in this paper.
Ultracold neutral atom interactions are characterized by pure s-wave scattering mediated at long-range by the dispersion forces [23,24]. Conversely, a bare ion creates a strong polarization force and hence the effective cross section is larger with significant contributions from higher-order partial waves [25]. Indeed the usual effective range expansion must be modified by logarithmic terms in the wavenumber expansion [26]. In the last few years theoretical studies of ultracold ion-atom collisions [27] included the investigation of the occurrence of magnetic Feshbach resonances with a view to examining the tunability of the ion-atom interaction focusing on the specific − + Ca Na 40 system [28,29], and calculations of the single-channel scattering properties of the + Ba ion with the Rb neutral atom [30] which suggest the possibility of sympathetic cooling of the barium ion by the buffer gas of ultracold rubidium atoms with a considerable efficiency.
In recent experiments [14,15,31], a single trapped ion of + Yb 174 in a Paul trap was immersed in a condensate of neutral Rb 87 atoms confined in a magneto-optical trap. A study of charge transfer cross sections showed that the simple classical Langevin model was inadequate to explain the reaction rates [15]. However, very little is known about the microscopic ultracold binary interactions between this ion and the rubidium atom [4,5].
The initial experimental study of the quantum coherence of charge transfer [15] was analysed using schematic energy curves as no accurate ab initio potentials existed. In particular, the potential energy curves (PECs) and couplings are not known with any accuracy. Thus the experimental study of the quantum coherence of charge transfer [15] was based on schematics of the energy curves. This prompted our in depth investigations [22] to map the lowest adiabatic states and the static properties of the molecular ion, in particular the turning points, potential minima, and crossing points of the lowest molecular energies. In addition to this, the dissociation energies and molecular constants provide useful spectroscopic data for dynamical investigations [32]. We have made a preliminary estimation of the pseudo-potential which approximates the ultracold interaction. This information is of great importance for modelling ultracold charge transfer, and in particular the quantum character of chemical reactivity and thus develop insights into ultracold quantum controlled chemistry [24], for example when external fields are applied to influence the reaction rates and reaction channels [15]. Of course, the presence of a bare charge in a dilute gas exposes many-body physics features such as exciton and polariton dynamics, which are also of great interest. It is also of great interest for laser manipulation of the collision to prevent losses through charge transfer or create translationally-cold trapped molecular ions via photoassociation. Since these processes are light-sensitive, then one can add an extra element of coherent control by using a laser to manipulate these processes [33,34].
In the present study we investigate radiative decay mechanisms, the charge-transfer process and the radiative association process using an optical-potential method.

Electronic structure calculation
Following our recent work on this molecular system [22] we extend those computations using a parallel version of the MOLPRO [35] suite of ab initio quantum chemistry codes (release MOLPRO 2010.1) to perform all the molecular structure calculations for this diatomic system + (Rb, Yb) . Low lying PECs as a function of internuclear distance out to R = 50 au are computed and in the present investigation we extend our earlier work [22] to calculate the transition dipole moments between the Σ + 1,3 states. As in our previous work we use an active two-electron model within a multi-reference configuration interaction (MRCI) and a full-configuration interaction (FCI) framework to calculate all the potentials. Briefly, within the MRCI model, the PECs are calculated using effective core potentials (ECP) to replace the nonvalence electrons (ECP68MDF for Yb, ECP36SDF for Rb), as a basis set for each atom, which allows for scalar-relativistic effects to be included explicitly. The scalar-relativistic effects are included by adding the corresponding terms of the Douglas-Kroll Hamiltonian to the one-electron integrals. To model the valence electrons, we use an augmented-correlation-consistent polarized valence basis set; aug-cc-pV6Z. We note that the basis set used yielded values consistent with those of Meyer and Bohn [36] for the neutral YbRb molecule. To take account of short-range interactions we employed the non-relativistic complete-active-space self consistent field (CASSCF)/MRCI method [37,38] available within the MOLPRO [35] ab initio quantum chemistry suite of codes. Figure 1 shows our adiabatic potential curves for this system as a function of the internuclear separation R. We note that the present quantum chemistry calculations and those of Sayfutyarova et al [32] use a similar approach. In summary, both calculations essentially use an effective core-potential and a multi-reference CI to cater for electron-correlation in the outer electrons. At short bond lengths all the results are obtained from the state-averaged CASSCF/MRCI approach. Our results are very similar to those obtained by Sayfutyarova et al [32] who used a total of 22 correlated electrons (14 electrons in doubly occupied orbitals). These authors also conducted CCSD(T) calculations on the Σ + X 1 , Σ + a 3 and Π b 3 states over 30 electrons to refine these potentials further. The effect of including extra electron correlation leads to a general reduction in the equilibrium bond lengths, most significant in the case of the Π b 3 state, where the combination of the CCSD (T) method and larger active electron calculation, predicts a much deeper potential well than found in our present work. However, this well nearly halved in size when spin-orbit effects were included. The Yb polarizability determined from the active two electron, MRCI and FCI potentials, is 128.5 au, (within 8% of the currently accepted value 139 ± 7 [39]) but smaller than the considerably more expensive CCSD(T) calculations of Sayfutyarova et al [32] who obtained α = 142.2. We note that both calculated values are within the experimental limits.
In figure 2 we illustrate the two singlet states involved in the radiative decay processes. The radiative charge transfer occurs along the Σ + A 1 state which has a shallow well. Comparing our results with the equivalent MRCI potential of Sayfutyarova et al [32] their results are in good agreement with ours. For this well, our earlier work [22] found a dissociation energy was . Although our excited state well is slightly deeper, there are significant differences in the Σ + X 1 ground-state. The dissociation energy of Sayfutyarova et al [32] is − 3496 cm 1 (0.4334 eV), almost twice the value of Lamb and co-workers [22], who obtained a value of 0.2202 eV. We note also that the equilibrium distance R e also  + YbRb potential energy curves (relative electronic energies) as a function of internuclear distance R (a ) 0 , (MRCI approximation) for the Σ + X 1 and Σ + A 1 states. The singlet Σ + A 1 is the entrance channel leading to the radiative charge-transfer channels: + + Rb Yb, the lowest energy ionic ytterbium states. The Σ + X 1 correlates for large R with the + Rb ion. The loss process from state Σ + A 1 can be through radiative association into the bound rovibrational manifold or above the dissociation threshold into the ion-atom charge exchange.
occurs at a shorter bond length of a 8.088 0 [32], compared to our finding which gave = R a 9.031 e 0 [22]. Figure 3 illustrates the dipole transition moment D (R) (au) as a function of internuclear separation R for the singlet and triplet Σ + states. Results for the dipole matrix elements between the Σ + X 1 -Σ + A 1 states and the Σ + a 3 -Σ + 2 3 states are illustrated. Comparing these results with the work of Sayfutyarova et al [32] we find very good agreement qualitatively for the − A X moment, although it is not possible to compare the triplet-triplet transition moment. Our results find a very smooth − A X singlet dipole which leads us to conclude that the sensitivity to the dynamics will be due to the wave function envelope. If there were some oscillatory behaviour in figure 3, then one could anticipate that this might be transferred to the radiative coupling. However, in the absence of structures in the moment, the resonance behaviour will be primarily potential scattering.
In figure 4 the transition rate calculated using equation (13) is presented as a function of the internuclear separation R. The decay rate Γ R ( ) decreases exponentially as R increases due to the exponential attenuation in the overlap of the atomic wave functions corresponding to charge transfer. Beyond R = 50 au, the potential of the Σ + A 1 state can be described by the long-range multipole expansion: where α d is the dipole polarizability of the neutral atom and where C 6 and C 8 are respectively the quadrupole and octupole polarizabilities, which have been evaluated in our previous study [22]. In our calculation of the phase shift we integrate into the asymptotic regime using the multipole series for the potential.

Theoretical method
In the simple classical model [40], the nuclear motion takes place on the incoming potential surface, V R ( ) A . Thus the motion is angular-momentum conserving, time-reversal invariant, and elastic-to a first approximation. Defining the collision energy, in the centre-of-mass frame, as E and the reduced mass of the nuclei as μ, then we can take the zero of potential energy at infinite separation in the incoming channel: Since angular momentum is conserved, then for an impact parameter b, the radial velocity can be written as: Thus the classical turning point will be the (largest) solution of the equation: The process of spontaneous emission has a rate Γ R ( ) which drives the charge transfer process. Consider a classical trajectory for a given collision energy, E, and impact parameter. Then for the decay process, letting t = 0 denotes the classical turning point where = ±∞ t are the end points. One can write for the probability of emission, for example as explained in [40]: R R c Then to a good approximation, in the case of weak coupling, states as a function of internuclear distance R (a ) 0 . The multi-reference-configurationinteraction (MRCI) approximation within the MOLPRO suite of codes is used to calculate the transition dipole moments. transitions.
we have: R R c Therefore, the semi-classical cross-section is simply,  (9)  isotopes (ignoring any resonant behaviour) then on dynamical grounds the cross-section is slightly higher. It is purely by coincidence that this energy dependence matches the classical Langevin model [42] for reactive collisions. For a long-range ion-atom potential we have the polarization potential, In the Langevin model, if the centrifugal barrier can be surmounted then the reaction proceeds with certainty and the cross-section is given by the simple formula: L d which displays the same energy dependence as (9) but based on completely different physics. Strictly speaking we have three quantum fields: the active electron, the nuclear motion, and the photon. One can construct the wavefunctions in the product (adiabatic) representation and then couple these through the Hamiltonian (including the vacuum photon states). However the process involves a weak-coupling, the irreversible spontaneous emission leading to charge transfer. Thus the collision of the + Yb ion with the Rb atom leading to loss of the + Yb ion can be considered as a second-order perturbation of the elastic lossless collision. The modified optical potential will have an imaginary (non-Hermitian) term proportional to the Einstein coefficient. This 'width' depends on the dipole moment and frequency of emission and is R-dependent. The optical potential method, in the context of radiative charge transfer, has been described in detail by Zygelman and Dalgarno [41]. We simply present the outline of the main equations and how it is modified for our application.
In the adiabatic approximation the dynamics occur on decoupled, centrally-symmetric PECs. The temperatures are so low that all non-adiabatic radial and rotational coupling are so weak that the vacuum coupling (by photoemission) is the only non-elastic process. Radiative charge transfer requires the optical dipole selection rules to be obeyed for transitions to the Σ + 1 state. Thus only the Σ + A 1 state has an allowed radiative charge transfer.
Using conventional notation, we use E to denote the collision energy in the centre-of-mass frame, and with m i and m a denoting the ion and atom masses, respectively, the reduced mass is defined: μ = + m m m m /( ) i a i a . Then the collision wavenumber is denoted by μ = k E 2 . Finally, all potential energies are with respect to the asymptotic incoming channel: . This means that the physics is essentially reduced to a single channel (effective complex radial potential) scattering problem: A A In the use of the simple optical potential we implicitly, and approximately, take into account both the process of radiative association and radiative transfer. That is the lower (exothermic charge exchange) state has an (infinite) number of bound (association) rovibrational levels and continuum states. This point is discussed in detail in previous applications [27,32] and its validity verified. In other terms, Γ, which is larger the higher the photon frequency, is taken as a vertical transition in analogy to the way that the 'reflection principle' is applied [43]. This approximation is better the larger the mass of the colliding atoms/ions. The problem can be summarized mathematically as [27] the solution of the Schrödinger equation, . Again using atomic units, we have that: and + Rb (4p S) 6 1 . Thus, Γ R ( ) is short range and exponentially damped with increasing R since it requires the electron to transfer from the atom to the ion. As the potential is central, even though it is complex, the usual separation in spherical coordinates applies, for example: , for the incoming channel A with angular momentum J, as follows: .
Then, without fear of ambiguity, we define the collision wavenumber: Then the corresponding radial functions, , will be the solutions of the equations: and δ J is the elastic phase shift. When the optical potential is used the radial equations for the functions in (14) are the same: apart from the modification for the complex wavenumber: We solve the problem directly integrating (19) using the Numerov method [45][46][47] and this is labelled the quantal approximation to distinguish it from the distorted-wave calculation and the semi-classical approximation discussed above.
The cross section for total collision-induced radiative decay from the entrance channel, the sum of the cross sections forr processes (1) and (2) can be obtained within the optical potential approximation. The cross section for collisioninduced radiative decay can then be written as, where k is given by (16), and g is the spin (statistical) weight.
Since the loss channel is via the Σ + X 1 state, and as the Yb ion and Rb atom combine to produce singlets, then only ion-atom collisions with singlet symmetry have a dipole-allowed spontaneous emission. So in this case, the statistical weight is, = g 1/4. At higher energy, a semi-classical approximation is invoked to calculate the cross sections for radiative decay.
The summation over the angular momentum in equation (22) can be replaced by an integral over the impact parameter, b, according to ≈ kb J. The JWKB approximation can then be used to obtain the wave function, This simplifies the calculation of the phase-shift, equation (21) [41,48,49] since the rapidly varying integrand gives us (in the classically allowed region): . Then using (21) we get the semiclassical approximation (9). The thermally averaged rate coefficient α , as a function of temperature T, is obtained by averaging over the Maxwell-Boltzmann distribution. That is, In the early work of Bates [49] an efficient and convenient procedure to evaluate the rate coefficient was outlined. In the present calculations for cross sections we start from − 10 12 μeV and extend these to higher energies, by invoking a semiclassical approximation above about − 10 2 eV up to 10 4 eV for the transition of interest.

Electronic states
Since the MRCI calculations do not explicitly include relativistic effects, although this is not important for the entrance collision channel or the lower Yb ( 1 S) + + Rb ( S) 1 asymptote as all the molecular states formed are of Σ + symmetry. This is borne out by the calculated energy of the asymptotic energies of the Σ + a 3 and Σ + A 1 states [22]. The asymptotes for the higher Π 3 and Σ + 3 states correlate to the Yb (6s6p 3 P o ) + + Rb (4p 6 S) 1 atomic products. The multiplet and its associated finestructure splitting in the triplet (Yb: 3 P o 0,1,2 ) is considerable: ∼0.3 eV. Only a fully relativistic treatment can accurately account for the spin-orbit interaction. In a magnetic trap of course the Zeeman splitting and hyperfine structure complicates matters further. Nonetheless, in our first analysis of this novel system, we can confidently say that a curve crossing will take place between the Σ + A 1 and Π b 3 states though at an energy above the + Yb ( S) 2 + Rb ( S) 2 asymptote. Such a crossing will facilitate a charge exchange reaction as observed in experiment at mK temperatures [14,15]. In our previous work on this complex [22] we have estimated the molecular constants for the four states that support bound rovibrational states.

Cross sections and collision rates
Cross sections were determined using the quantal optical potential approximation, for collision energies ranging from μ − 10 12 eV up to 10 eV. At higher collision energies a semiclassical approximation (equation (12)) was invoked for energies up to 10 keV in order to determine the cross sections.
In order to compare with experiment the statistical weight for the singlet, must be taken into account, this means taking = g 1/4 in (22). The cross section results presented in figures 5 and 6 for all the radiative decay process are the spinless results, g = 1. In figure 5 we show the optical potential results as they are mapped on to those obtained from the semi-classical approximation. Figure 6 illustrates the low partial-wave contributions to the total cross section in the optical potential approximation. The solid line and dashed line show the effect of the long-range interactions. The solid line is the full multipole expansion, while the dashed line only includes the quadrupole of polarization. At the higher energies, the corrections to leading-order polarization do not affect the positions of the resonances. This confirms that the resonance effect as short-range features, well described as potential scattering. At the very lowest energies however, as the wavelengths become extremely long and the centrifugal barrier much more significant, then, as is well-known, the long-range features of the potential take over. From these results one clearly sees that at collision energies below μ − 10 6 eV the cross section is totally dominated by s-wave scattering. However, owing to the sensitivity of the scattering to the potential, the estimation of the effective range parameters, including the complex scattering length, would be of great interest. So primarily, at intermediate energies the process is dominated by the short-range classical turning point. This may go some way to explaining why our dynamic results are in such good agreement with the more complex calculations of [32] although the agreement is surprising. Note, at higher energies non-adiabatic effects will naturally become more important, however, in this energy (temperature) range, the role of non-adiabatic coupling (radial and rotational) turns out to have very little importance indeed. Their influence is negligible for the radiative capture process, as has been shown recently in detailed studies by Sayfutyarova et al [32].
Finally, we consider whether thermal effects might be taken into account, to confirm the discrepancy between theory and experiment. In the experiments by Köhl et al [14,15] the kinetic energy of a single + Yb ion immersed in an ultracold Rb ensemble was varied by adding excess micromotion energy after displacement of an ion from the centre of a trap. The binary-collision ion-loss rate coefficient determined in this way does not correspond to a conventional thermally-averaged rate constant for the Maxwell collision energy distribution at a certain temperature which assumes thermal equilibrium. The relationship σ v is therefore used to designate an effective energy-dependent rate coefficient [32], where R E ( ) is given by, A X R We use this form to define a quasi-rate coefficient [24] rather than one averaged over a Maxwellian distribution defined in equation (24). In figure 7 we compare our calculations with experiment for this quasi-rate parameter R. The measured experimental value for the + Yb 174 isotope [14] shown in figure 7 indicate that the magnitude of R(E) is ± × − (4.0 0.3) 10 14 − cm s 3 1 , where as the calculations of Sayfutyarova et al [32] for this isotope give a value, after averaging the cross sections in the energy region 0.15-3.25 − cm 1 , through the savannah of resonances, a value of × − 2.9 10 14 − cm s 3 1 , which is just outside the experimental error. Carrying out a similar procedure with our cross sections results (solid black . The quantal optical potential calculations (22) are compared with the semi-classical approximation. In the figure we present the spinless cross sections, that is g = 1. The background of the quantal result (22) follows the semi-classical approximation (9) at collision energies above μeV and has the asymptotic − E 1/2 behaviour. . Again, in this case we present the spinless partial cross section, equation (22), with g = 1.
line, figure 7) yields a mean value slightly higher, in better agreement with experiment, but relatively close in magnitude to previous work [32]. For the . As found in previous studies [32] our theoretical predictions for this same ratio give a value of ≈1.03 indicating to the contrary.
Regarding the large isotope sensitivity, we do not observe as great a difference as reported in experimental measurements. Similar to the detailed calculations of Sayfutyarova et al [32] we find a dense forest of resonances, but no broad features that would lead to a strong isotope dependency as observed in the experiment. So, our conclusion is that this feature remains unexplained. It is possible to speculate that the experimental conditions, having the magnetic field present, create additional complications. For example the nuclear spin of the Yb ion. It is known from recent experiments that this has an important role in relaxation phenomena [1]. Further detailed theoretical studies and additional experiments would be essential in order to resolve this issue.
We note for + YbCa ultra-cold collisions [27] it is only at temperatures below a nano-Kelvin ( − 10 6 K) that a large isotope effect is seen. Above these temperatures there is a very small isotope effect. In the present work on + YbRb , the energy range is 0.15-3.25 − cm 1 (0.215-4.676 Kelvin), so a similar small isotope effect is seen as in + YbCa [27].

Conclusions
We have investigated the quantum nature of ultracold ion-atom collisions and calculated the relevant experimental processescross sections and rate coefficients for the different isotopes of + Yb ion. These calculations are important in the design and interpretation of the new generation of experiments involving ultracold ion-atom systems. PECs and transition dipole moments obtained from the MOLPRO suite of codes for lowlying molecular states of the diatomic molecular ionic system containing a ytterbium ion and a rubidium atom, with relevance to ultra-cold chemistry were used in our dynamical calculations. Cross sections as a function of energy, for the radiative decay, charge transfer and association processes involving Yb ions and Rb atoms are determined using an optical potential method. The MRCI approach is used to determine turning points, crossing points, potential minima and spectroscopic molecular constants for the lowest five molecular states. The long-range parameters, including the dispersion coefficients estimated from our ab initio data were used in our dynamical investigations. Quasienergy dependent rate coefficients are determined from our cross section for the radiative decay processes in ultracold collisions of a ytterbium ion and a rubidium atom based on our ab initio data and compared with the available experimental measurements and previous theoretical work [32]. The agreement is surprisingly good for the molecular electronic structure that gives rise to the complex optical potential. The smooth nature of the optical potential and validity of the semi-classical approximation indicates that one can accurately estimate the cross section with an elementary quadrature. The more complex quantal treatment, while exhibiting the expected potential resonances, does not give rise to a strong isotope effect, at least in the energy range we investigated. The estimates of the energy dependent collision rate are in suitable agreement with experiment [14,15] and with previous theoretical studies [32]. We find no broad resonances features that might underly a strong isotope effect. In conclusion, we find, in agreement with previous theoretical work [32] that the isotope anomaly observed in experimental studies remains unexplained.  isotope. Although theory shows suitable agreement with the experimental measurements, the strong isotope sensitivity observed in the experiment is not evident in the present calculations.
YbRb + MRCI potentials and dipole moments states. Calculations were performed at the multi-reference-configuration-interaction (MRCI) level as a function of internuclear separation R (a 0 ) with the MOLPRO suite of codes [35], see text for details. The potential energy curves (PECs) are calculated using an effective core potentials (ECP) to replace the non-valence electrons (ECP68MDF for Yb, ECP36SDF for Rb) and an AV6Z basis for the outer electrons.