Importance of long-range interactions in chemical reactions at cold and ultracold temperatures

Accurate diabatization based on combined-hyperbolic-inverse-power-representation: 1,2 2A' states of BeH2

A diabatic potential energy matrix (DPEM) for the two lowest states of BeH2+ has been constructed using the combined-hyperbolic-inverse-power-representation (CHIPR) method. By imposing symmetry constraints on the coefficients of polynomials, the complete nuclear permutation inversion symmetry is correctly preserved in the CHIPR functional form. The symmetrized CHIPR functional form is then used in the diabatization by ansatz procedure. The ab initio energies are reproduced with satisfactory accuracy. In addition, the CHIPR-based DPEM also reproduces the local topology of a conical intersection. Future work will focus on a complete four-state diabatic representation with emphasis on the long-range interactions and spin-orbit couplings, which will enable accurate quantum scattering calculations for the Be+(2P) + H2 → BeH+(X1Σ+) + H(2S) reaction.

Charge-Optimized Electrostatic Interaction Atom-Centered Neural Network Algorithm

Machine-learning algorithms have been proposed to capture electrostatic interactions by using effective partial charges. These algorithms often rely on a pretrained model for partial charge prediction using density functional theory-calculated partial charges as references, which introduces complexity to the force field model. The accuracy of the trained model also depends on the reliability of charge partition methods, which can be dependent on the specific system and methodology employed. In this study, we propose an atom-centered neural network (ANN) algorithm that eliminates the need for reference charges. Our algorithm requires only a single NN model for each element to obtain both atomic energy and charges. These atomic charges are then employed to compute electrostatic energies using the Ewald summation algorithm. Subsequently, the force field model is trained on total energy and forces, with the inclusion of electrostatic energy. To evaluate the performance of our algorithm, we conducted tests on three benchmark systems, including a Ge slab with an O adatom system, a TiO2 crystalline system, and a Pd-O nanoparticle system. Our results demonstrate reasonably accurate predictions of partial charges and electrostatic interactions. This algorithm provides a self-consistent charge prediction strategy and possibilities for robust and reliable modeling of electrostatic interactions in machine-learning potentials.

Signatures of Non-universal Quantum Dynamics of Ultracold Chemical Reactions of Polar Alkali Dimer Molecules with Alkali Metal Atoms: Li(2S) + NaLi(a3Σ+) → Na(2S) + Li2(a3Σu+)

Ultracold chemical reactions of weakly bound triplet-state alkali metal dimer molecules have recently attracted much experimental interest. We perform rigorous quantum scattering calculations with a new ab initio potential energy surface to explore the chemical reaction of spin-polarized NaLi(a³Σ⁺) and Li(²S) to form Li₂(a³Σ_u⁺) and Na(²S). The reaction is exothermic and proceeds readily at ultralow temperatures. Significantly, we observe strong sensitivity of the total reaction rate to small variations of the three-body part of the Li₂Na interaction at short range, which we attribute to a relatively small number of open Li₂(a³Σ_u⁺) product channels populated in the reaction. This provides the first signature of highly non-universal dynamics seen in rigorous quantum reactive scattering calculations of an ultracold exothermic insertion reaction involving a polar alkali dimer molecule, opening up the possibility of probing microscopic interactions in atom+molecule collision complexes via ultracold reactive scattering experiments.

Quantum-Controlled Collisions of H2 Molecules

The amount of information that can be obtained from a scattering experiment depends upon the precision with which the quantum states are defined in the incoming channel. By precisely defining the incoming states and measuring the outgoing states in a scattering experiment, we set up the boundary condition for experimentally solving the Schrödinger equation. In this Perspective we discuss cold inelastic scattering experiments using the most theoretically tractable H₂ and its isotopologues as the target. We prepare the target in a precisely defined rovibrational (v, j, m) quantum state using a special coherent optical technique called the Stark-induced adiabatic Raman passage (SARP). v and j represent the quantum numbers of the vibrational and rotational energy levels, and m refers to the projection of the rotational angular momentum vector j on a suitable quantization axis in the laboratory frame. Selection of the m quantum numbers defines the alignment of the molecular frame, which is necessary to probe the anisotropic interactions. For us to achieve the collision temperature in the range of a few degrees Kelvin, we co-expand the colliding partners in a mixed supersonic beam that is collimated to define a direction for the collision velocity. When the bond axis is aligned with respect to a well-defined collision velocity, SARP achieves stereodynamic control at the quantum scale. Through various examples of rotationally inelastic cold scattering experiments, we show how SARP coherently controls the dynamics of anisotropic interactions by preparing quantum superpositions of the orientational m states within a single rovibrational (v, j) energy state. A partial wave analysis, which has been developed for the cold scattering experiments, shows dominance of a resonant orbital that leaves its mark in the scattering angular distribution. These highly controlled cold collision experiments at the single partial wave limit allow the most direct comparison with the results of theoretical computations, necessary for accurate modeling of the molecular interaction potential.

Anisotropic dynamics of resonant scattering between a pair of cold aligned diatoms

The collision dynamics between a pair of aligned molecules in the presence of a partial-wave resonance provide the most sensitive probe of the long-range anisotropic forces important to chemical reactions. Here we control the collision temperature and geometry to probe the dynamics of cold (1-3 K) rotationally inelastic scattering of a pair of optically state-prepared D₂ molecules. The collision temperature is manipulated by combining the gating action of laser state preparation and detection with the velocity dispersion of the molecular beam. When the bond axes of both molecules are aligned parallel to the collision velocity, the scattering rate drops by a factor of 3.5 as collision energies >2.1 K are removed, suggesting a geometry-dependent resonance. Partial-wave analysis of the measured angular distribution supports a shape resonance within the centrifugal barrier of the l = 2 incoming orbital. Our experiment illustrates the strong anisotropy of the quadrupole-quadrupole interaction that controls the dynamics of resonant scattering.

Stereodynamical Control of Cold Collisions of Polyatomic Molecules with Atoms

Scattering between atomic and/or molecular species can be controlled by manipulating the orientation or alignment of the collision partners. Such stereodynamics is particularly pronounced at cold (∼1 K) collision temperatures because of the presence of resonances. Comparing to the extensively studied atomic and diatomic species, polyatomic molecules with strong steric anisotropy could provide a more sophisticated platform for studying such stereodynamics. Here, we provide the quantum mechanical framework for understanding state-to-state stereodynamics in rotationally inelastic scattering of polyatomic molecules with atoms and apply it to cold collision of oriented H₂O with He on a highly accurate potential energy surface. It is shown that strong stereodynamical control can be achieved near 1 K via shape resonances. Furthermore, quantum interference in scattering of a coherently prepared initial state of the H₂O species is explored, which is shown to be significant.

Bimolecular Chemistry in the Ultracold Regime

Advances in atomic, molecular, and optical physics techniques allowed the cooling of simple molecules down to the ultracold regime ([Formula: see text]1 mK) and opened opportunities to study chemical reactions with unprecedented levels of control. This review covers recent developments in studying bimolecular chemistry at ultralow temperatures. We begin with a brief overview of methods for producing, manipulating, and detecting ultracold molecules. We then survey experimental works that exploit the controllability of ultracold molecules to probe and modify their long-range interactions. Further combining the use of physical chemistry techniques such as mass spectrometry and ion imaging significantly improved the detection of ultracold reactions and enabled explorations of their dynamics in the short range. We discuss a series of studies on the reaction KRb + KRb → K₂ + Rb₂ initiated below 1 [Formula: see text]K, including the direct observation of a long-lived complex, the demonstration of product rotational state control via conserved nuclear spins, and a test of the statistical model using the complete quantum state distribution of the products. Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Cold and controlled chemical reaction dynamics

State-to-state chemical reaction dynamics, with complete control over the reaction parameters, offers unparalleled insight into fundamental reactivity.

Heuristic optimization of analytic laser pulses for vibrational stabilization of ultracold KRb

We proposed a methodology that allows to maximize the population transfer from a high vibrational state of the a ³Σ⁺ triplet state to the vibrational ground state of the X ¹Σ⁺ singlet state though the optimization of one pump and one dump laser pulses. The pump pulse is optimized using a fitness function, heuristically improved, that includes the effect of the spin-orbit coupling of the KRb [b-A]-scheme. The dump pulse is optimized to maximize the population transfer to the ground state. We performed a comparison with the case in which the pump and dump pulses are optimized to maximize the population transfer to the ground state employing a genetic algorithm with a single fitness function. The heuristic approach turned out to be 70% more efficient than a quantum optimal control optimization employing a single fitness function. The method proposed provides simple pulses that have an experimental realm.

Cold quantum-controlled rotationally inelastic scattering of HD with H2 and D2 reveals collisional partner reorientation

Molecular interactions are best probed by scattering experiments. Interpretation of these studies has been limited by lack of control over the quantum states of the incoming collision partners. We report here the rotationally inelastic collisions of quantum-state prepared deuterium hydride (HD) with H₂ and D₂ using a method that provides an improved control over the input states. HD was coexpanded with its partner in a single supersonic beam, which reduced the collision temperature to 0-5 K, and thereby restricted the involved incoming partial waves to s and p. By preparing HD with its bond axis preferentially aligned parallel and perpendicular to the relative velocity of the colliding partners, we observed that the rotational relaxation of HD depends strongly on the initial bond-axis orientation. We developed a partial-wave analysis that conclusively demonstrates that the scattering mechanism involves the exchange of internal angular momentum between the colliding partners. The striking differences between H₂/HD and D₂/HD scattering suggest the presence of anisotropically sensitive resonances.

Geometric Phase Appears in the Ultracold Hydrogen Exchange Reaction

Quantum reactive scattering calculations for the hydrogen exchange reaction H+H_{2}(v=4,j=0)→H+H_{2}(v^{'}, j^{'}) and its isotopic analogues are reported for ultracold collision energies. Because of the unique properties associated with ultracold collisions, it is shown that the geometric phase effectively controls the reactivity. The rotationally resolved rate coefficients computed with and without the geometric phase are shown to differ by up to 4 orders of magnitude. The effect is also significant in the vibrationally resolved and total rate coefficients. The dynamical origin of the effect is discussed and the large geometric phase effect reported here might be exploited to control the reactivity through the application of external fields or by the selection of a particular nuclear spin state.

Quantum dynamics study on the CHIPR potential energy surface for the hydroperoxyl radical: the reactions O + OH⇋O2 + H

Quantum scattering calculations of the O((3)P)+OH((2)Π)⇌O2((3)Σg (-))+H((2)S) reactions are presented using the combined-hyperbolic-inverse-power-representation potential energy surface [A. J. C. Varandas, J. Chem. Phys. 138, 134117 (2013)], which employs a realistic, ab initio-based, description of both the valence and long-range interactions. The calculations have been performed with the ABC time-independent quantum reactive scattering computer program based on hyperspherical coordinates. The reactivity of both arrangements has been investigated, with particular attention paid to the effects of vibrational excitation. By using the J-shifting approximation, rate constants are also reported for both the title reactions.

Tuning ultracold chemical reactions via Rydberg-dressed interactions

We show that ultracold chemical reactions with an activation barrier can be tuned using Rydberg-dressed interactions. Scattering in the ultracold regime is sensitive to long-range interactions, especially when weakly bound (or quasibound) states exist near the collision threshold. We investigate how, by Rydberg dressing a reactant, one enhances its polarizability and modifies the long-range van der Waals collision complex, which can alter chemical reaction rates by shifting the position of near-threshold bound states. We carry out a full quantum mechanical scattering calculation for the benchmark system H(2)+D, and show that resonances can be moved substantially and that rate coefficients at cold and ultracold temperatures can be increased by several orders of magnitude.

Dynamics of the D(+) + H2 → HD + H(+) reaction at the low energy regime by means of a statistical quantum method

The D(+) +H2(v = 0, j = 0, 1) → HD+H(+) reaction has been investigated at the low energy regime by means of a statistical quantum mechanical (SQM) method. Reaction probabilities and integral cross sections (ICSs) between a collisional energy of 10(-4) eV and 0.1 eV have been calculated and compared with previously reported results of a time independent quantum mechanical (TIQM) approach. The TIQM results exhibit a dense profile with numerous narrow resonances down to Ec ~ 10(-2) eV and for the case of H2(v = 0, j = 0) a prominent peak is found at ~2.5 × 10(-4) eV. The analysis at the state-to-state level reveals that this feature is originated in those processes which yield the formation of rotationally excited HD(v' = 0, j' > 0). The statistical predictions reproduce reasonably well the overall behaviour of the TIQM ICSs at the larger energy range (Ec ≥ 10(-3) eV). Thermal rate constants are in qualitative agreement for the whole range of temperatures investigated in this work, 10-100 K, although the SQM values remain above the TIQM results for both initial H2 rotational states, j = 0 and 1. The enlargement of the asymptotic region for the statistical approach is crucial for a proper description at low energies. In particular, we find that the SQM method leads to rate coefficients in terms of the energy in perfect agreement with previously reported measurements if the maximum distance at which the calculation is performed increases noticeably with respect to the value employed to reproduce the TIQM results.

Symmetry-adapted perturbation theory based on unrestricted Kohn-Sham orbitals for high-spin open-shell van der Waals complexes

Two open-shell formulations of the symmetry-adapted perturbation theory are presented. They are based on the spin-unrestricted Kohn-Sham (SAPT(UKS)) and unrestricted Hartree-Fock (SAPT(UHF)) descriptions of the monomers, respectively. The key reason behind development of SAPT(UKS) is that it is more compatible with density functional theory (DFT) compared to the previous formulation of open-shell SAPT based on spin-restricted Kohn-Sham method of Żuchowski et al. [J. Chem. Phys. 129, 084101 (2008)]. The performance of SAPT(UKS) and SAPT(UHF) is tested for the following open-shell van der Waals complexes: He···NH, H(2)O···HO(2), He···OH, Ar···OH, Ar···NO. The results show an excellent agreement between SAPT(UKS) and SAPT(ROKS). Furthermore, for the first time SAPT based on DFT is shown to be suitable for the treatment of interactions involving Π-state radicals (He···OH, Ar···OH, Ar···NO). In the interactions of transition metal dimers ((3)Σ(u)(+))Au(2) and ((13)Σ(g)(+))Cr(2) we show that SAPT is incompatible with the use of effective core potentials. The interaction energies of both systems expressed instead as supermolecular UHF interaction plus dispersion from SAPT(UKS) result in reasonably accurate potential curves.

Kinetics and dynamics of the S(1D2) + H2 → SH + H reaction at very low temperatures and collision energies

We report combined studies on the prototypical S(1D2) + H2 insertion reaction. Kinetics and crossed-beam experiments are performed in experimental conditions approaching the cold energy regime, yielding absolute rate coefficients down to 5.8 K and relative integral cross sections to collision energies as low as 0.68 meV. They are supported by quantum calculations on a potential energy surface treating long-range interactions accurately. All results are consistent and the excitation function behavior is explained in terms of the cumulative contribution of various partial waves.

Ion-molecule chemistry at very low temperatures: cold chemical reactions between Coulomb-crystallized ions and velocity-selected neutral molecules

The recent development of a range of techniques for producing cold atoms and molecules at very low translational temperatures T < or = 1 K has provided the opportunity to investigate collisional processes in a new physical regime. We have recently presented a new experimental method to study low-temperature reactive collisions between translationally cold ions and neutral molecules (S. Willitsch et al., Phys. Rev. Lett. 2008, 100, 043203). Our technique relies on the combination of a quadrupole-guide velocity selector for the generation of translationally cold neutral molecules with a facility to produce ordered structures of cold ions (Coulomb crystals) by laser cooling in a linear quadrupole ion trap. The strong localisation of the ions in the trap in combination with the high sensitivity of laser-induced-fluorescence detection enabled us to study chemical reactions on the single-particle level, down to temperatures of T approximately 1 K. In the current paper, we present a detailed characterisation of the scope and limitations of this method based on our study of the reaction between laser-cooled Ca+ ions and velocity-selected CH3F molecules. The properties of our cold-neutrals source and the dependence of the measured rate constant on the shape of the Coulomb crystals, trapping and laser-cooling parameters are discussed. An extension of our technique for the study of low-temperature reactions with sympathetically cooled molecular ions (translational temperature T > 10 mK) is presented and first results on the charge-transfer reaction between OCS+ and ND3 are discussed. Finally, perspectives for further developments of our method are explored.

Molecular collisions, from warm to ultracold

This introductory article contrasts molecular collisions, particularly reactive collisions, in the familiar "warm" domain with the ultracold regime where the relative deBroglie wavelengths become long compared with the range of interaction of the collision partners. Ultracold collisions have much greater sensitivity to entrance channel interactions, so offer the prospect of tuning by external fields to control onset of reaction. However, for ultracold collisions, kinematic constraints impose severe limitations on the observable dynamical properties. In the exit channel for appreciably exoergic reactions, the deBroglie wavelengths become short, so the exit dynamics are much like those for warm collisions. Reactions of alkali dimers, halides, and monoxide molecules are discussed that seem especially congenial for cold collision studies.

Chemical applications of laser- and sympathetically-cooled ions in ion traps

Ensembles of cold atomic and molecular ions in ion traps prepared at millikelvin temperatures by laser and sympathetic cooling have recently found considerable interest in both physics and chemistry. At very low temperatures the ions form ordered structures in the trap also known as "Coulomb crystals". Ion Coulomb crystals exhibit a range of intriguing properties which render them attractive systems for novel experiments in chemical dynamics, ultrahigh-resolution spectroscopy and quantum-information processing. In this article we review the methods used to prepare atomic and molecular ion Coulomb crystals and discuss some recent studies in mass spectrometry, low-temperature chemistry and precision spectroscopy to illustrate their scientific potential for chemical applications. Finally, we conclude with an outlook on outstanding challenges and prospective further developments in the field.

Vibrational cooling of spin-stretched dimer states by He buffer gas: quantum calculations for Li2(a 3Sigma(u)+) at ultralow energies

The interaction between the triplet state of the lithium dimer, (7)Li(2), with (4)He is obtained from accurate ab initio calculations where the vibrational dependence of the potential is newly computed. Vibrational quenching dynamics within a coupled-channel quantum treatment is carried out at ultralow energies, and large differences in efficiency as a function of the initial vibrational state of the targets are found as one compares the triplet results with those of the singlet state of the same target.

Cold Reactive Collisions between laser-cooled ions and velocity-selected neutral molecules

We report a new experimental method to study reactive ion-molecule collisions at very low temperatures. A source of laser-cooled ions in a linear Paul trap has been combined with a quadrupole-guide velocity selector to investigate the reaction of Ca+ with CH3F at collision energies E[over](coll)/k(B)> or =1 K with single-particle sensitivity. The technique represents a general approach to study reactive collisions between ions and polar molecules over a wide temperature range down to the cold regime.