Quantum Dynamics and Kinetics of the F + H2 and F + D2 Reactions at Low and Ultra-Low Temperatures

Stereodynamical control of cold HD + D2 collisions

We report full-dimensional quantum calculations of stereodynamic control of HD(v = 1, j = 2) + D2 collisions that has been probed experimentally by Perreault et al. using the Stark-induced adiabatic Raman passage (SARP) technique. Computations were performed on two highly accurate full-dimensional H4 potential energy surfaces. It is found that for both potential surfaces, rotational quenching of HD from with concurrent rotational excitation of D2 from is the dominant transition with cross sections four times larger than that of elastically scattered D2 for the same quenching transition in HD. This process was not considered in the original analysis of the SARP experiments that probed ΔjHD = -2 transitions in HD(vHD = 1, jHD = 2) + D2 collisions. Cross sections are characterized by an l = 3 resonance for ortho-D2(jD2 = 0) collisions, while both l = 1 and l = 3 resonances are observed for the para-D2(jD2 = 1) partner. While our results are in excellent agreement with prior measurements of elastic and inelastic differential cross sections, the agreement is less satisfactory with the SARP experiments, in particular for the transition for which the theoretical calculations indicate that D2 rotational excitation channel is the dominant inelastic process.

Quantum Dynamics Studies of the Li + Na2 (V = 0, j = 0) → Na + NaLi Reaction on a New Neural Network Potential Energy Surface

The ultracold reaction offers a unique opportunity to elucidate the intricate microscopic mechanism of chemical reactions, and the Na2Li system serves as a pivotal reaction system in the investigation of ultracold reactions. In this work, a high-precision potential energy surface (PES) of the Na2Li system is constructed based on high-level ab initio energy points and the neural network (NN) method, and a proper asymptotic functional form is adopted for the long-range interaction, which is suitable for the study of cold or ultracold collisions. Based on the new NN PES, the dynamics of the Li + Na2 (v = 0, j = 0) → Na + NaLi reaction are studied in the collision energy range of 10-7 to 80 cm-1. In the high collision energy range of 8 to 80 cm-1, the dynamics of the reaction is studied using the time-dependent wave packet method and the statistical quantum mechanical (SQM) method. Comparing the results of the two methods, it is found that the SQM method provides a rough description of the product ro-vibrational state distribution but overestimates the integral cross-section values. With the decrease of collision energy, the reaction differential cross section gradually changes from forward-backward symmetric scattering to predominant forward scattering. In the low collision energy range from 10-7 to 8 cm-1, the SQM method is used to study the reaction dynamics, and the rate constant in the Wigner threshold region is estimated to be 2.87 × 10-10 cm3/s.

Coupled 3D (J ≥ 0) Time-Dependent Wave Packet Calculation for the F + H2 Reaction on Accurate Ab Initio Multi-State Diabatic Potential Energy Surfaces

We had calculated adiabatic potential energy surfaces (PESs), nonadiabatic, and spin-orbit (SO) coupling terms among the lowest three electronic states (12A', 22A', and 12A″) of the F + H2 system using the multireference configuration interaction (MRCI) level of theory, and the adiabatic-to-diabatic transformation equations were solved to formulate the diabatic Hamiltonian matrix [J. Chem. Phys. 2020, 153, 174301] for the entire region of the nuclear configuration space. The accuracy of such diabatic PESs is explored by performing scattering calculations to evaluate integral cross sections (ICSs) and rate constants. The nonadiabatic and SO effects are studied by utilizing coupled 3D time-dependent wave packet formalism with zero and nonzero total angular momentum on multiple adiabatic/diabatic surfaces calculation. We depict the convergence profiles of reaction probabilities for the reactive as well as nonreactive processes on various electronic states at different collision energies with respect to total angular momentum including all helicity quantum numbers. Finally, total ICSs are calculated as functions of collision energies for the initial rovibrational state (v = 0, j = 0) of the H2 molecule along with the temperature-dependent rate coefficient, where those quantities are compared with previous theoretical and experimental results.

A single resonance Regge pole dominates the forward-angle scattering of the state-to-state F + H2 → FH + H reaction at Etrans = 62.09 meV

The aim of the present paper is to bring clarity, through simplicity, to the important and long-standing problem: does a resonance contribute to the forward-angle scattering of the F + H2 reaction? We reduce the problem to its essentials and present a well-defined, yet rigorous and unambiguous, investigation of structure in the differential cross sections (DCSs) of the following three state-to-state reactions at a translational energy of 62.09 meV: F + H2(vi = 0, ji = 0, mi = 0) → FH(vf = 3, jf = 0, 1, 2, mf = 0) + H, where vi, ji, mi and vf, jf, mf are the initial and final vibrational, rotational and helicity quantum numbers respectively. Firstly, we carry out quantum-scattering calculations for the Fu-Xu-Zhang potential energy surface, obtaining accurate numerical scattering matrix elements for indistinguishable H2. The calculations use a time-independent method, with hyperspherical coordinates and an enhanced Numerov method. Secondly, the following theoretical techniques are employed to analyse structures in the DCSs: (a) full and Nearside-Farside (NF) partial wave series (PWS) and local angular momentum theory, including resummations of the full PWS up to second order. (b) The recently introduced "CoroGlo" test, which lets us distinguish between glory and corona scattering at forward angles for a Legendre PWS. (c) Six asymptotic (semiclassical) forward-angle glory theories and three asymptotic farside rainbow theories, valid for rainbows at sideward-scattering angles. (d) Complex angular momentum (CAM) theories of forward and backward scattering, with the Regge pole positions and residues computed by Thiele rational interpolation. Thirdly, our conclusions for the three PWS DCSs are: (a) the forward-angle peaks arise from glory scattering. (b) A broad (hidden) farside rainbow is present at sideward angles. (c) A single Regge pole contributes to the DCS across the whole angular range, being most prominent at forward angles. This proves that a resonance contributes to the DCSs for the three transitions. (d) The diffraction oscillations in the DCSs arise from NF interference, in particular, interference between the Regge pole and direct subamplitudes.

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.

Time-dependent wave packet dynamics study of the resonances in the H + LiH+(v = 0, j = 0) → Li+ + H2 reaction at low collision energies

The depletion process of LiH⁺ by H collision plays an important role in the evolution of the early universe and astrophysical processes, including the eventual charge-states, abundances of atomic and molecular species and ensuing astrochemistry. Here, a quantum dynamics study on the H + LiH⁺(v = 0, j = 0) → Li⁺ + H₂ reaction is performed at the low collision energy range from 0.1 meV to 10 meV using the time-dependent wave packet method. A Feshbach resonance peak is observed near 0.8 meV collision energy on the total reaction probability curves. This resonance originates from the coupling with the v = 0, j = 1 energy level of the reactant LiH⁺, and it is dominated by the contributions of J = 0-4 partial waves. Another partial wave resonance is also found on the total integral cross section at 1.2 meV, which is closely connected to the opening of the J = 7 partial wave. The opening of the J = 7 partial wave generates a notable forward scattering peak, and the Feshbach resonance can promote both the forward and backward scatterings. Moreover, the total and product vibrational state-resolved rate coefficients for the temperature range of 1-100 K are also reported.

Accurate Calculation of Rate Constant and Isotope Effect for the F + H2 Reaction by the Coupled 3D Time-Dependent Wave Packet Method on the Newly Constructed Ab Initio Ground Potential Energy Surface

We employ coupled three-dimensional (3D) time dependent wave packet formalism in hyperspherical coordinates for reactive scattering problem on the newly constructed ab initio calculated ground adiabatic potential energy surface for the F + H₂/D₂ reaction. The convergence profiles for various reactive channels are depicted at low collision energy regimes with respect to the total angular momentum (J) quantum numbers. For two different reactant diatomic molecules (H₂ and D₂) initially at their respective ground roto-vibrational state (v = 0, j = 0), calculated state-to-state as well as total integral cross sections as a function of collision energy, temperature dependent rate constants, and the kinetic isotope effect for various reactivity profiles of F + H₂ and F + D₂ reactions are presented along with previous theoretical and experimental results.

Wave Packet Approach to Adiabatic and Nonadiabatic Dynamics of Cold Inelastic Scatterings

Due to the extremely large de Broglie wavelength of cold molecules, cold inelastic scattering is always characterized by the time-independent close-coupling (TICC) method. However, the TICC method is difficult to apply to collisions of large molecular systems. Here, we present a new strategy for characterizing cold inelastic scattering using wave packet (WP) method. In order to deal with the long de Broglie wavelength of cold molecules, the total wave function is divided into interaction, asymptotic and long-range regions (IALR). The three regions use different numbers of ro-vibrational basis functions, especially the long-range region, which uses only one function corresponding to the initial ro-vibrational state. Thus, a very large grid range can be used to characterize long de Broglie wavelengths in scattering coordinates. Due to its better numerical scaling law, the IALR-WP method has great potential in studying the inelastic scatterings of larger collision systems at cold and ultracold regimes.

Feshbach resonances in D + HD(v = 1, j = 0) reaction at low collision energies

Feshbach resonances in D + HD(v = 1, j = 0) reaction are studied by using the time-independent quantum method. The integral cross section (ICS) results present three Feshbach resonance peaks, which are different from H + HD(v = 1, j = 0) reaction dominated by only one peak. These resonances are attributed to coupling with adiabatic effective potentials of D + HD(v = 1, j = 1) reaction, and the most obvious peak is contributed by J = 1 at 83.16 cm^-1 collision energy. For J = 0 and 2, the resonances are related with the same L partial wave and present a double-peak structure in total ICS. The characteristics of product angular distribution show that the resonance of J = 1 is long-lived, while the lifetimes are relatively shorter for the resonance of J = 0 and 2.

Quantum Wave Packet Treatment of Cold Nonadiabatic Reactive Scattering at the State-To-State Level

Cold and ultracold collisions are dominated by quantum effects, such as resonances, tunneling, and nonadiabatic transitions between different electronic states. Due to the extremely long de Broglie wavelength in such processes, quantum reactive scattering is most conveniently characterized using the time-independent close-coupling (TICC) methods. However, the TICC approach is difficult for systems with a large number of channels because of its steep numerical scaling laws. Here, a recently proposed quantum wave packet (WP) approach for solving adiabatic reactive scattering problems at low collision energies is extended to include nonadiabatic transitions. To impose the outgoing boundary conditions, the total scattering wavefunction is split into three parts, the interaction, the asymptotic, and the long-range regions. Each region is associated with a different set of basis functions, which could be optimized separately. In this way, an extremely long grid can be used to accommodate the characteristic long de Broglie wavelengths in the scattering coordinate. The better numerical scaling laws of the WP approach have the potential for handling larger nonadiabatic reactive systems at low temperatures in the future.

The F + HD(v = 0, 1; j = 0, 1) reactions: stereodynamical properties of orbiting resonances

The excitation functions (reaction cross-section as a function of collision energy) of the F + HD(v = 0, 1; j = 0, 1) benchmark system have been calculated in the 0.01-6 meV collision energy interval using a time-independent hyperspherical quantum dynamics methodology. Special attention has been paid to orbiting resonances, which bring about detailed information on the three-atom interaction during the reactive encounter. The location of the resonances depends on the rovibrational state of the reactants HD(v,j), but is the same for the two product channels HF + D and DF + H, as expected for these resonances that are linked to the van der Waals well at the entrance. The resonance intensities depend both on the entrance and on the exit channels. The peak intensities for the HF + D channel are systematically larger than those for DF + H. Vibrational excitation leads to an increase of the peak intensity by more than an order of magnitude, but rotational excitation has a less drastic effect. It deceases the resonance intensity of the F + HD(v = 1) reaction, but increases somewhat that of F + HD(v = 0). Polarization of the rotational angular momentum with respect to the initial velocity reveals intrinsic directional preferences in the F + HD(v = 0, 1; j = 1) reactions that are manifested in the resonance patterns. The helicities (Ω = 0, Ω = ±1) possible for j = 1 contribute to the resonances, but that from Ω± 1 is, in general, dominant and in some cases exclusive. It corresponds to a preferential alignment of the HD internuclear axis perpendicular to the initial direction of approach and, thus, to side-on collisions. This work also shows that external preparation of the reactants, following the intrinsic preferences, would allow the enhancement or reduction of specific resonance features, and would be of great help for their eventual experimental detection.

Insights into the kinetics and molecular mechanism of the Newman–Kwart rearrangement

Kinetics and molecular mechanism of the Newman–Kwart rearrangement (NKR) of N,N-dimethyl O-arylthiocarbamate into N,N-dimethyl S-arylcarbamate have been investigated theoretically.

Advances and New Challenges to Bimolecular Reaction Dynamics Theory

Dynamics of bimolecular reactions in the gas phase are of foundational importance in combustion, atmospheric chemistry, interstellar chemistry, and plasma chemistry. These collision-induced chemical transformations are a sensitive probe of the underlying potential energy surface(s). Despite tremendous progress in past decades, our understanding is still not complete. In this Perspective, we survey the recent advances in theoretical characterization of bimolecular reaction dynamics, stimulated by new experimental observations, and identify key new challenges.

From the Kinetic Theory of Gases to the Kinetics of Rate Processes: On the Verge of the Thermodynamic and Kinetic Limits

A variety of current experiments and molecular dynamics computations are expanding our understanding of rate processes occurring in extreme environments, especially at low temperatures, where deviations from linearity of Arrhenius plots are revealed. The thermodynamic behavior of molecular systems is determined at a specific temperature within conditions on large volume and number of particles at a given density (the thermodynamic limit): on the other side, kinetic features are intuitively perceived as defined in a range between the extreme temperatures, which limit the existence of each specific phase. In this paper, extending the statistical mechanics approach due to Fowler and collaborators, ensembles and partition functions are defined to evaluate initial state averages and activation energies involved in the kinetics of rate processes. A key step is delayed access to the thermodynamic limit when conditions on a large volume and number of particles are not fulfilled: the involved mathematical analysis requires consideration of the role of the succession for the exponential function due to Euler, precursor to the Poisson and Boltzmann classical distributions, recently discussed. Arguments are presented to demonstrate that a universal feature emerges: Convex Arrhenius plots (super-Arrhenius behavior) as temperature decreases are amply documented in progressively wider contexts, such as viscosity and glass transitions, biological processes, enzymatic catalysis, plasma catalysis, geochemical fluidity, and chemical reactions involving collective phenomena. The treatment expands the classical Tolman’s theorem formulated quantally by Fowler and Guggenheim: the activation energy of processes is related to the averages of microscopic energies. We previously introduced the concept of “transitivity”, a function that compactly accounts for the development of heuristic formulas and suggests the search for universal behavior. The velocity distribution function far from the thermodynamic limit is illustrated; the fraction of molecules with energy in excess of a certain threshold for the description of the kinetics of low-temperature transitions and of non-equilibrium reaction rates is derived. Uniform extension beyond the classical case to include quantum tunneling (leading to the concavity of plots, sub-Arrhenius behavior) and to Fermi and Bose statistics has been considered elsewhere. A companion paper presents a computational code permitting applications to a variety of phenomena and provides further examples.

Benchmark Quantum Kinetics at Low Temperatures toward Absolute Zero and Role of Entrance Channel Wells on Tunneling, Virtual States, and Resonances: The F + HD Reaction

This paper reports a study of the quantum reaction dynamics and kinetics of the F + HD reaction at low and ultralow temperatures, focusing on the range from the Wigner limit up to 50 K. Close coupling time-independent quantum reactive scattering calculations for the production of HF and DF molecules have been carried out on two potential energy surfaces differing in the description of the reaction entrance channel. This case is computationally more demanding than the cases of F with H₂ and D₂ ( De Fazio et al. Frontiers in Chemistry 2019 , 7 , 328 ) but offers a wider phenomenology regarding the roles of quantum mechanical effects of tunneling, of virtual states, and of resonances. The results show that at the temperatures in the cold and ultracold regimes small changes in the entrance channel long-range interaction induce surprising near threshold features. The presence of a virtual state close to the reactive threshold gives rise to a marked anti-Arrhenius behavior of the rate constants below 100 mK. This effect enhances reaction rates by about 2 orders of magnitude, making them of the same order as those at room temperature and confining the onset of the Wigner regime in the microkelvin region.

Orbiting resonances in the F + HD (v = 0, 1) reaction at very low collision energies. A quantum dynamical study

Time-independent, fully converged, quantum dynamical calculations have been performed for the F + HD (v = 0, j = 0) and F + HD (v = 1, j = 0) reactions on an accurate potential energy surface down to collision energies of 0.01 meV. The two isotopic exit channels, HF + D and DF + H, have been investigated. The calculations reproduce satisfactorily the Feshbach resonance structures for collision energies between 10 and 40 meV, previously reported in the literature for the HF + D channel. Contrary to the results of a former literature work, vibrational excitation of HD is found to enhance reactivity in all cases down to the lowest collision energy investigated. Shape-type orbiting resonances are found for collision energies lower than 2 meV. The resonances appear as peaks in the reaction cross sections that are associated to specific values of the total angular momentum, J. In contrast with the Feshbach resonances at higher energies, the orbiting resonance structure, which is caused by the van der Waals well of the entrance channel, is identical for the HF + D and DF + H exit channels. The orbiting resonance peaks for F + HD (v = 0) are very small, but those for F + HD (v = 1) could be observed, in principle, with a combination of Raman pumping and merged beams methods.