Complexes formed in collisions between ultracold alkali-metal diatomic molecules and atoms

We explore the properties of 3-atom complexes of alkali-metal diatomic molecules with alkali-metal atoms, which may be formed in ultracold collisions. We estimate the densities of vibrational states at the energy of atom-diatom collisions, and find values ranging from 2.2 to 350~K$^{-1}$. However, this density does not account for electronic near-degeneracy or electron and nuclear spins. We consider the fine and hyperfine structure expected for such complexes. The Fermi contact interaction between electron and nuclear spins can cause spin exchange between atomic and molecular spins. It can drive inelastic collisions, with resonances of three distinct types, each with a characteristic width and peak height in the inelastic rate coefficient. Some of these resonances are broad enough to overlap and produce a background loss rate that is approximately proportional to the number of outgoing inelastic channels. Spin exchange can increase the density of states from which laser-induced loss may occur.

Characterizing quasibound states and scattering resonances

Characterizing quasibound states from coupled-channel scattering calculations can be a laborious task, involving extensive manual iteration and fitting. We present an automated procedure, based on the phase shift or S-matrix eigenphase sum, that reliably converges on a quasibound state (or scattering resonance) from some distance away. It may be used for both single-channel and multichannel scattering. It produces the energy and width of the state and the phase shift of the background scattering, and hence the lifetime of the state. It also allows extraction of partial widths for decay to individual open channels. We demonstrate the method on a very narrow state in the Van der Waals complex Ar--H$_2$, which decays only by vibrational predissociation, and on near-threshold states of $^{85}$Rb$_2$, whose lifetime varies over 4 orders of magnitude as a function of magnetic field.

Prospects of forming high-spin polar molecules from ultracold atoms

We have investigated Feshbach resonances in collisions of high-spin atoms such as Er and Dy with closed-shell atoms such as Sr and Yb, using coupled-channel scattering and bound-state calculations. We consider both low-anisotropy and high-anisotropy limits. In both regimes we find many resonances with a wide variety of widths. The wider resonances are suitable for tuning interatomic interactions, while some of the narrower resonances are highly suitable for magnetoassociation to form high-spin molecules. These molecules might be transferred to short-range states, where they would have large magnetic moments and electric dipole moments that can be induced with very low electric fields. The results offer the opportunity to study mixed quantum gases where one species is dipolar and the other is not, and open up important prospects for a new field of ultracold high-spin polar molecules.

Microwave shielding of ultracold polar molecules with imperfectly circular polarization

We investigate the use of microwave radiation to produce a repulsive shield between pairs of ultracold polar molecules and prevent collisional losses that occur when molecular pairs reach short range. We carry out coupled-channels calculations on RbCs+RbCs and CaF+CaF collisions in microwave fields. We show that effective shielding requires predominantly circular polarization, but can still be achieved with elliptical polarization that is around 90% circular.

Time delays in ultracold atomic and molecular collisions

We study the behavior of the Eisenbud-Wigner collisional time delay around Feshbach resonances in cold and ultracold atomic and molecular collisions. We carry out coupled-channels scattering calculations on ultracold Rb and Cs collisions. In the low-energy limit, the time delay is proportional to the scattering length, so exhibits a pole as a function of applied field. At high energy, it exhibits a Lorentzian peak as a function of either energy or field. For narrow resonances, the crossover between these two regimes occurs at an energy proportional to the square of the resonance strength parameter $s_textrm{res}$. For wider resonances, the behavior is more complicated and we present an analysis in terms of multichannel quantum defect theory.

User Manual for MOLSCAT, BOUND and FIELD, Version 2020.0: programs for quantum scattering properties and bound states of interacting pairs of atoms and molecules

MOLSCAT is a general-purpose package for performing non-reactive quantum scattering calculations for atomic and molecular collisions using coupled-channel methods. Simple atom-molecule and molecule-molecule collision types are coded internally and additional ones may be handled with plug-in routines. Plug-in routines may include external magnetic, electric or photon fields (and combinations of them). Simple interaction potentials are coded internally and more complicated ones may be handled with plug-in routines. BOUND is a general-purpose package for performing calculations of bound-state energies in weakly bound atomic and molecular systems using coupled-channel methods. It solves the same sets of coupled equations as MOLSCAT, and can use the same plug-in routines if desired, but with different boundary conditions. FIELD is a development of BOUND that locates external fields at which a bound state exists with a specified energy. One important use is to locate the positions of magnetically tunable Feshbach resonance positions in ultracold collisions. Versions of these programs before version 2019.0 were released separately. However, there is a significant degree of overlap between their internal structures and usage specifications. This manual therefore describes all three, with careful identification of parts that are specific to one or two of the programs.

MOLSCAT: a program for non-reactive quantum scattering calculations on atomic and molecular collisions

MOLSCAT is a general-purpose program for quantum-mechanical calculations on nonreactive atom-atom, atom-molecule and molecule-molecule collisions. It constructs the coupled-channel equations of atomic and molecular scattering theory, and solves them by propagating the wavefunction or log-derivative matrix outwards from short range to the asymptotic region. It then applies scattering boundary conditions to extract the scattering matrix (S matrix). Built-in coupling cases include atom + rigid linear molecule, atom + vibrating diatom, atom + rigid symmetric top, atom + asymmetric or spherical top, rigid diatom + rigid diatom, rigid diatom + asymmetric top, and diffractive scattering of an atom from a crystal surface. Interaction potentials may be specified either in program input (for simple cases) or with user-supplied routines. For the built-in coupling cases, MOLSCAT can loop over partial wave (or total angular momentum) to calculate elastic and inelastic cross integral sections and spectroscopic line-shape cross sections. Post-processors are available to calculate differential cross sections, transport, relaxation and Senftleben-Beenakker cross sections, and to fit the parameters of scattering resonances. MOLSCAT also provides an interface for plug-in routines to specify coupling cases (Hamiltonians and basis sets) that are not built in; plug-in routines are supplied to handle collisions of a pair of alkali-metal atoms with hyperfine structure in an applied magnetic field. For low-energy scattering, MOLSCAT can calculate scattering lengths and effective ranges and can locate and characterize scattering resonances as a function of an external variable such as the magnetic field.

BOUND and FIELD: programs for calculating bound states of interacting pairs of atoms and molecules

The BOUND program calculates the bound states of a complex formed from two interacting particles using coupled-channel methods. It is particularly suitable for the bound states of atom-molecule and molecule-molecule Van der Waals complexes and for the near-threshold bound states that are important in ultracold physics. It uses a basis set for all degrees of freedom except $R$, the separation of the centres of mass of the two particles. The Schrodinger equation is expressed as a set of coupled equations in $R$. Solutions of the coupled equations are propagated outwards from the classically forbidden region at short range and inwards from the classically forbidden region at long range, and matched at a point in the central region. Built-in coupling cases include atom + rigid linear molecule, atom + vibrating diatom, atom + rigid symmetric top, atom + asymmetric or spherical top, rigid diatom + rigid diatom, and rigid diatom + asymmetric top. Both programs provide an interface for plug-in routines to specify coupling cases (Hamiltonians and basis sets) that are not built in. With appropriate plug-in routines, BOUND can take account of the effects of external electric, magnetic and electromagnetic fields, locating bound-state energies at fixed values of the fields. The related program FIELD uses the same plug-in routines and locates values of the fields where bound states exist at a specified energy. As a special case, it can locate values of the external field where bound states cross scattering thresholds and produce zero-energy Feshbach resonances. Plug-in routines are supplied to handle the bound states of a pair of alkali-metal atoms with hyperfine structure in an applied magnetic field.

Microwave shielding of ultracold polar molecules

We use microwaves to engineer repulsive long-range interactions between ultracold polar molecules. The resulting shielding suppresses various loss mechanisms and provides large elastic cross sections. Hyperfine interactions limit the shielding under realistic conditions, but a magnetic field allows suppression of the losses to below 10-14 cm3 s-1. The mechanism and optimum conditions for shielding differ substantially from those proposed by Gorshkov et al. [Phys. Rev. Lett. 101, 073201 (2008)], and do not require cancelation of the long-range dipole-dipole interaction that is vital to many applications.

Hyperfine structure of 2Sigma molecules containing alkaline-earth atoms

Ultracold molecules with both electron spin and an electric dipole moment offer new possibilities in quantum science. We use density-functional theory to calculate hyperfine coupling constants for a selection of molecules important in this area, including RbSr, LiYb, RbYb, CaF and SrF. We find substantial hyperfine coupling constants for the fermionic isotopes of the alkaline-earth and Yb atoms. We discuss the hyperfine level patterns and Zeeman splittings expected for these molecules. The results will be important both to experiments aimed at forming ultracold open-shell molecules and to their applications.