We use quantum and classical adiabatic capture theories to study the chemical reaction Li + CaH LiH + Ca. Using a recently developed ab initio potential energy surface, which provides an accurate representation of long-range interactions in the entrance reaction channel, we calculate the adiabatic channel potentials by diagonalizing the Li-CaH Hamiltonian as a function of the atom-molecule separation. The resulting adiabatic channel potentials are used to calculate both the classical and quantum capture probabilities as a function of collision energy, as well as the temperature dependencies of the partial and total reaction rates. The calculated reaction rate agrees well with the measured value at 1 K (V Singh et al 2012 Phys. Rev. Lett. 108 203201), suggesting that the title reaction proceeds without an activation barrier. The calculated classical adiabatic capture rate agrees well with the quantum result in the multiple-partial-wave regime of relevance to the experiment. Significant differences are found only in the ultracold limit ( mK), demonstrating that adiabatic capture theories can predict the reaction rates with nearly quantitative accuracy in the multiple-partial-wave regime.
In this thesis we investigate theoretically the cold collision properties of metastable neon (Ne*). The collisions of metastable rare gases (Rg*) exhibit a plethora of interesting phenomena: Due to the high internal excitation energies of the metastable atoms ionization processes are very likely. The produced ions can be detected with single-ion efficiency in the experiment. In experiments with atomic Rg* gases at cold temperatures T~1 mK, this allows for in-situ monitoring of the real-time dynamics of the atomic gas and for example for the study of two-particle correlations. Furthermore, in cold Rg* gases one can examine the prospects of creating a Bose-Einstein condensate or a degenerate Fermi gas of Rg*. Ne* is a promising candidate for Bose-Einstein condensation and in the group of G. Birkl the properties of cold gases of Ne* are investigated experimentally. Two-body loss rate coefficients due to the ionization processes and elastic cross sections have been measured. It is useful to complement the experimental measurements with a theoretical study in order to obtain a better understanding of the cold Ne* gases and the collision physics of Ne*. At the low temperatures of the atomic gases, these gases are dilute with particle number densities of n~10^9 cm^{-3}. For these low densities, the behavior of the gas is determined mainly by two-body collision physics which is solved by quantum scattering theory in terms of the S matrix. The collision physics is given by short-range interactions and long-range interactions of the atoms. The short-range and long-range interaction potentials of Ne* have been calculated in the molecular basis of Ne2. In this thesis we first calculate short-range and long-range interactions of atoms and molecular basis states of diatomic molecules with a simpler electronic structure than Ne* and Ne2 in order to introduce the notation and to discuss the characteristics of the calculated molecular interaction potentials and in order to employ them in the scattering calculations for Ne*. Instead of solving the full scattering equations by taking into account all the molecular interaction potentials, we demonstrate that the collision physics of Ne* can be described in terms of a coupled two-channel model with a single interaction potential only which describes the elastic scattering of Ne*. In this model, ionization in Ne* collisions is described by the transition of the upper elastic interaction channel to the lower channel, representing the loss or ionization channel. We introduce two versions of the two-channel model in this work. In the first version, the two channels are given by square-well potentials. We can solve the scattering equations analytically and for complex wave numbers k and study the solutions in the complex k plane. With an expansion of the S matrix in terms of its poles in the complex k plane we find a parametrization of the two-body loss rate coefficients. We show that these coefficients describe ionizing collisions of Ne* by comparing them to the experimental measurements and to the numerical results obtained by the two-channel model in the second version. In this second version of the two-channel model, the elastic scattering channel is given by a realistic interaction potential of Ne*, consisting of a calculated short-range and long-range molecular potential, and the loss channel by a model ionization potential. For this model we calculate the two-body loss rate coefficients and the elastic cross sections for the isotope mixtures of Ne* which have been measured experimentally and demonstrate that the free potential parameters can be optimized to the experimental data to achieve very good agreement of the numerical results with the experimental measurements. We discuss the validity of the two-channel model by comparing the results to existing models for cold Rg* collisions and find that the two-channel model of this work is a useful extension to these models.
The last decade has seen substantial advances in the exploration of inelastic and reactive collisions at collision energies in the cold ( <
< 1 K) and ultracold ( < 1
<1 mK) regime.1-13 These advances ushered in an understanding of cold and ultracold collisions that is impacting research areas ranging from molecular spectroscopy to quantum simulation and computation to quantum control of chemical reaction dynamics. Cold collisions are dominated by resonances, as amply manifest in most of the systems studied.1, 14 Resonances are also present in the thermal ( >
> 100 K) and hyperthermal ( >
> 1000 K) collision energy regimes,15 but obscured by background (non-resonant) scattering and the presence of many partial waves pertaining to higher angular momenta that blur the partial waves' individual contributions. Moreover, the outcome of thermal and hyperthermal collisions is somewhat insensitive to the details of the interaction potential – the potential energy surface (PES) – especially to its attractive long-range region. Only in a handful of cases, such as the F+HD exchange reaction, resonances in the thermal regime have been experimentally detected and theoretically accounted for by quantum mechanical scattering calculations.16, 17 In contrast, in the cold regime, collisions are dominated by only a few partial waves whose effects on the integral and differential scattering cross-sections can be singled out and linked to particular resonances. The cold collision energy regime has been attained by making use of crossed molecular beams that were either Stark- or Zeeman-decelerated or whose intersection angle could be reduced all the way down to the merging-beams limit. A beam emerging from a Stark or Zeeman decelerator has a variable, well-defined laboratory velocity. When crossed with a secondary beam at small angles (5 ∘
∘ –10 ∘
∘ ), collision energies in the 0.3–15 K could be readily achieved, as demonstrated by the group of van de Meerakker.9, 11 The technique of merged beams makes use of the ability to bend one of the beams by a magnetic or electric field such that the crossing angle becomes zero, thereby affording collision energies in the range 0.1 K–10 K. Merging beams have been used in studies of Penning ionization by the groups of Narevicius5, 12 and Osterwalder18 and made it possible to unravel shape and Feshbach resonances arising in this process. The techniques mentioned so far require the reactants to be polar or paramagnetic. For reactants that are neither, collision energies as low as 5 K could be achieved by taking advantage of a favorable combination of low intersection angles and low reduced mass of the reactants. This approach has been demonstrated by Bergeat, Costes, Naulin, and coworkers in a series of reactive3, 19 and inelastic collision experiments8, 20 of H 2
2 . Yet, another approach to achieving low relative collision velocities relies on the coexpansion of the collision partners in a single molecular beam.21 If their masses are similar, both partners will end up with nearly the same velocities, within a few ms − 1
−1 . Using this technique, Zare, Mukherjee, and coworkers carried out a series of experiments in which rotational quenching of HD by H 2
2 , D 2
2 , and He was investigated. Moreover, in these experiments, the HD molecule was prepared in v
v =1 and j
j =2 using Stark-field-induced adiabatic Raman passage that allowed, in addition, to align the rotational angular momentum of the HD reactant molecule with respect to the initial relative velocity of the HD + He reactants.22-24 Collisional energy transfer between rotation and translation, in which a molecule is excited or relaxed via scattering by an atom or another molecule, has been of particular significance for developing our understanding of inelastic collisions. Most notably, the scattering of NO molecules by rare gas (Rg) atoms has served as a benchmark system for both experimentalists and theorists for decades. In its ground electronic state, NO is a 2 Π
2Π molecule with a single unpaired electron in an antibonding π ∗
π∗ orbital. In this state, the rotational levels are split into two spin–orbit manifolds, defined by the projection | Ω |
∣Ω∣ ( = 1 2
=21 or 3 2
23 ) of the electronic angular momentum (orbital and spin) on the internuclear axis. The rotational levels with the same j
j in the two manifolds are separated by approximately 123 cm − 1
−1 for v
v =0 (117 cm − 1
−1 for v
v =10). In addition, each rotational level within each manifold is further split into two Λ
Λ -doublet states ( ε = ± 1
ϵ=±1 , e/f), which are nearly degenerate, but differ in parity, p = ( − 1 ) j − ε / 2
p=(−1)j−ϵ/2 . As the Rg atom approaches the NO molecule, the degeneracy of the Λ
Λ doublets is lifted, giving rise to two PESs, A ′
A′ and A ′ ′
A′′ , in which the orbital of the unpaired electron lies, respectively, in the plane of the three atoms or perpendicular to it. In computational practice, it is more convenient to use linear combinations (sums and differences) of the A ′ ′
A′′ and A ′
A′ PESs. Using hexapole focusing or Stark deceleration combined with detection via velocity ion imaging, fully state-to-state-resolved integral and differential cross sections for rotationally elastic and inelastic collisions could be measured with high resolution.9, 11, 25, 26 More than that, the orientation and alignment of the scattered (outgoing) NO molecules could be measured as a function of the scattering angle27-30 and, quite recently, the effects of orientation of the reactant (incoming) NO molecules on the differential cross section in collisions with Rg, mostly Ar, could be established.31-35 On the theory side, quantum mechanical close coupling (QMCC) scattering calculations on accurate two-dimensional ab initio PESs using the benchmark coupled-cluster, CCSD(T), method with the NO internuclear distance fixed at its equilibrium value36 were performed and found to account for all the details revealed by the experiments. Although state of the art, most of these experiments except for those described in Refs. 9, 11 were carried out at thermal energies for NO molecules in their lowest level, | v = 0 , j = 1 2 , Ω = 1 2 , f ⟩
∣v=0,j=21,Ω=21,f⟩ . Scattering of highly vibrationally excited molecules in the cold energy regime remained an unexplored frontier until the recent work by the groups of Arthur Suits, Ad van der Avoird, Gerrit Groenenboom, and Hua Guo who added spin–orbit and rotational state selection and looked for spin- and rotational-state changing collisions to boot. Published in the current issue of Natural Sciences,37 this work is a culmination of a 3-year collaboration that included both implementation of new experimental techniques and boosting of theoretical methods. Working with the NO+Ar system, the authors made use of stimulated emission pumping via the electronically excited A 2 Σ +
A2Σ+ state to promote NO( 2 Π 1 2
2Π21 ) molecules to the vibrationally excited state with v = 10
v=10 and a selected spin–orbit state | Ω |
∣Ω∣ and rotational state j
j , | v = 10 , j , Ω , f ⟩
∣v=10,j,Ω,f⟩ .38, 39 Low (cold) collision energies were achieved in a novel experimental setup of two beams (dual-beam arrangement) crossing at 4 ∘
∘ and intersected by the pump and dump laser beams of the stimulated-emission pumping scheme. Downstream, the collision products were probed using 1+1 ′
′ REMPI with the ionization laser tuned just above the ionization threshold. Fully state-resolved angular distributions of the scattered NO molecules were measured by velocity ion slice imaging and detected by a charge-coupled device (CCD). In order to be able to detect the product (outgoing) molecules resulting from the cold NO + Ar collisions, some excess energy had to be made available to the products. This came in the form of the initial spin–orbit and/or rotational state energy. In the study highlighted herein,37 angular distributions leading to the spin–orbit changing transition | v = 10 , Ω = 3 2 , j = 3 2 f ⟩ → | v ′ = 10 , Ω ′ = 1 2 , j ′ = 9 2 , e ⟩
∣v=10,Ω=23,j=23f⟩→∣v′=10,Ω′=21,j′=29,e⟩ were measured at energies in the 2 K–5 K range. Although intrabeam scattering could have been used as an alternative as it had been in previous work,21 the authors demonstrated that a dual-beam setup provides a superior angular and energy resolution. The highlighted experiment challenged the existing theoretical treatments of the NO+ Rg systems: first, the approximation of fixing NO at its equilibrium distance was bound to fail for NO in v = 10
v=10 due to the concomitant elongation of the NO bond. Second, for such a highly excited vibrational state, both the long- and short-range interactions are affected, triggering nonadiabatic effects involving rotational and spin–orbit degrees of freedom. The above rendered the single-reference coupled-cluster method used in previous treatments of NO( 2 Π 1 2 , v = 0
2Π21,v=0 ) inadequate. Therefore, the authors of the highlighted paper37 constructed a new coupled-cluster PES as well as a multireference configuration-interaction (MRCI) PES and carried out full-fledged (3D) QMCC calculations on both PESs. These tour-de-force calculations reproduced quite well the experimental differential cross sections and defined a new state of the art in scattering calculations. However, some discrepancies remain as at low collision energies, and the accuracy of the long-range potential needed is a challenge for any ab initio method. Moreover, treating spin-changing collisions that are dominated by the essentially attractive difference potential, V A ′ ′ − V A ′
VA′′−VA′ , puts additional demands on the accuracy of the A ′ ′
A′′ and A ′
A′ PESs. The calculated integral cross section in the 3–4 cm − 1
−1 energy interval exhibits a bump that may indicate the presence of resonances that show up as backscattering in the measured differential cross section. A partial-wave analysis points to overlapping resonances associated with contributions from several of the highest partial waves involved. In any case, tracing the resonances may require additional improvements of the existing PESs. All in all, the highlighted combined experimental and theoretical study sets a new standard for our understanding of cold collisions. There is no conflict of interest. The peer review history for this article is available at https://publons.com/publon/10.1002/ntls.20210088. This Highlight has been internally reviewed and edited by Bretislav Friedrich and Gerard Meijer. No data available.
Multichannel Quantum Defect Theory (MQDT) is shown to be capable of producing quantitatively accurate results for low-energy atom-molecule scattering calculations. With a suitable choice of reference potential and short-range matching distance, it is possible to define a matrix that encapsulates the short-range collision dynamics. Multichannel quantum defect theory can provide an efficient alternative to full coupled-channel calculations for low-energy molecular collisions. However, the efficiency relies on interpolation of the Y matrix that encapsulates the short-range dynamics. It is shown how the phases of the MQDT reference functions may be chosen so as to remove such poles from the vicinity of a reference energy and dramatically increase the range of interpolation. For the test cases of Mg+NH and Li+NH, the resulting optimized Y matrix may be interpolated smoothly over an energy range of several Kelvin and a magnetic field range of over 1000G. Calculations at additional energies and fields can then be performed at a computational cost that is proportional to the number of channels N and not to N cubed. MQDT thus provides a promising method for carrying out low-energy molecular scattering calculations on systems where full exploration of the energy and the field dependence is currently impractical.
We investigate the ultracold collisions of rotationally excited dipolar molecules in free-space, taking the hetero-nuclear bi-alkali molecule of KRb as an example. We show that we can sharply tune the elastic, inelastic and reactive rate coefficients of lossy molecular collisions when a second rotationally excited colliding channel crosses the threshold of the initial colliding channel, with the help of an applied electric field, as found by Avdeenkov et al for non-lossy molecules ( Phys. Rev. A 73 022707 ). We can increase or decrease the inelastic and reactive processes whether the second channel is above or below the initial channel. This is seen for both bosonic and fermionic molecules. Additionally, we include the electric quadrupole and octopole moment to the dipole moment in the expression of the long-range multipole–multipole interaction. We found that for processes mediated by the incident channel, such as elastic and reactive collisions, the inclusion of quadrupole and octopole moments is not important at ultralow energies. The moments are important for processes mediated by state-to-state transitions like inelastic collisions.