Quantum scattering calculations on chemical reactions

Higher-Order Split Operator Schemes for Solving Tetratomic Reactions Using the Time-Dependent Wave Packet Method

In this work, using the time-dependent quantum wave packet method, quite a few typical higher-order split operators (HOSOs) were for the first time applied to calculate the tetratomic reactive scattering processes in the hyperspherical coordinate. It was found that the HOSOs were hardly efficient for a tetratomic reaction calculation, unlike those for a triatomic reactive scattering calculation. We proposed an efficient HOSO with a force gradient (denoted as 2G1 in the main text) for efficiently and accurately calculating a tetratomic reaction using the quantum wave packet method. Several typical tetratomic reactions, such as H2 + OH, HF + OH, and H2 + OH+, are calculated for demonstrating the effectiveness of the proposed 2G1 in terms of (product state-resolved) reaction probability and inelastic probability, by comparing with the performance of the previously reported various HOSOs. We suggest that the 2G1 propagator could be applied to efficiently calculate a general tetratomic reaction.

Simulation of Chemical Reactions on a Quantum Computer

Studying chemical reactions, particularly in the gas phase, relies heavily on computing scattering matrix elements. These elements are essential for characterizing molecular reactions and accurately determining reaction probabilities. However, the intricate nature of quantum interactions poses challenges, necessitating the use of advanced mathematical models and computational approaches to tackle the inherent complexities. In this study, we develop and apply a quantum computing algorithm for the calculation of scattering matrix elements. In our approach, we employ the time-dependent method based on the Møller operator formulation where the S-matrix element between the respective reactant and product channels is determined through the time correlation function of the reactant and product Møller wavepackets. We successfully apply our quantum algorithm to calculate scattering matrix elements for 1D semi-infinite square well potential and on the colinear hydrogen exchange reaction. As we navigate the complexities of quantum interactions, this quantum algorithm is general and emerges as a promising avenue, shedding light on new possibilities for simulating chemical reactions on quantum computers.

Reactant-Product Decoupling Technique Using the Intermediate Coordinate Method

Although the reactant-product decoupling (RPD) technique was proposed over two decades ago, it remains an efficient approach for calculating product state-resolved information on some simple direct reactions using the quantum wave packet method. In the past, usually the RPD technique employed the collocation method to transform the wave function between reactant and product arrangements, which requires quite large computational efforts. In this work, the intermediate coordinate (IC) method is employed to realize the RPD technique. Numerical examples demonstrate that this new IC RPD (IRPD) technique has superior computational efficiency compared with the original method employing the collocation method. Especially, the new IRPD technique significantly saves disk space and computer memory. To illustrate the features of our new method, the total reaction probabilities of the H + H2, H + Br2, and F + H2 reactions with J = 0 and the differential cross sections of the H + H2 and F + H2 reactions at a series of collision energy are calculated and presented. With this efficient and effective new RPD technique, the Li + HF reaction, which involves sharp resonances with long-range wave functions in the van der Waals wells in both the reactant and product arrangements, is also calculated with several J at the product state-resolved level to reveal the ability of the RPD technique for describing resonance wave functions. With these numerical examples, it is found that, for the reaction with resonances, the RPD approach should be applied carefully. Otherwise, it is very possible that the resonances could disappear with the application of the RPD technique.

Simulating chemical reaction dynamics on quantum computer

The electronic energies of molecules have been successfully evaluated on quantum computers. However, more attention is paid to the dynamics simulation of molecules in practical applications. Based on the variational quantum eigensolver (VQE) algorithm, Fedorov et al. proposed a correlated sampling (CS) method and demonstrated the vibrational dynamics of H2 molecules [J. Chem. Phys. 154, 164103 (2021)]. In this study, we have developed a quantum approach by extending the CS method based on the VQE algorithm (labeled eCS-VQE) for simulating chemical reaction dynamics. First, the CS method is extended to the three-dimensional cases for calculation of first-order energy gradients, and then, it is further generalized to calculate the second-order gradients of energies. By calculating atomic forces and vibrational frequencies for H2, LiH, H+ + H2, and Cl- + CH3Cl systems, we have seen that the approach has achieved the CCSD level of accuracy. Thus, we have simulated dynamics processes for two typical chemical reactions, hydrogen exchange and chlorine substitution, and obtained high-precision reaction dynamics trajectories consistent with the classical methods. Our eCS-VQE approach, as measurement expectations and ground-state wave functions can be reused, is less demanding in quantum computing resources and is, therefore, a feasible means for the dynamics simulation of chemical reactions on the current noisy intermediate-scale quantum-era quantum devices.

Theoretical Development of the Interaction-Asymptotic Region Decomposition Method for Tetratomic Reactive Scattering

An accurate and efficient time-dependent wave packet method is proposed for solving the product state-resolved reaction probabilities of the tetratomic reactive system. In this method, the entire scattering process is divided into the interaction region and multiple asymptotic regions, sharing the same spirit as the interaction-asymptotic region decomposition (IARD) approach in a triatomic reactive scattering process. The hyperspherical coordinate is adopted in the interaction region, while the corresponding Jacobi coordinate is employed in each asymptotic region. Therefore, in this IARD method, the "coordinate problem", the difficulty of expressing the wave function in the entire region using a single coordinate system, can be effectively avoided, and only a very small number of the grid points (or the basis functions) are required. For the numerical illustration, the typical tetratomic reaction H2 + OH with zero total angular momentum is calculated, and compared with other quantum wave packet methods. Our proposed IARD method for the tetratomic reactive system is much more efficient and accurate.

Accurate fundamental invariant-neural network representation of ab initio potential energy surfaces

Highly accurate potential energy surfaces are critically important for chemical reaction dynamics. The large number of degrees of freedom and the intricate symmetry adaption pose a big challenge to accurately representing potential energy surfaces (PESs) for polyatomic reactions. Recently, our group has made substantial progress in this direction by developing the fundamental invariant-neural network (FI-NN) approach. Here, we review these advances, demonstrating that the FI-NN approach can represent highly accurate, global, full-dimensional PESs for reactive systems with even more than 10 atoms. These multi-channel reactions typically involve many intermediates, transition states, and products. The complexity and ruggedness of this potential energy landscape present even greater challenges for full-dimensional PES representation. These PESs exhibit a high level of complexity, molecular size, and accuracy of fit. Dynamics simulations based on these PESs have unveiled intriguing and novel reaction mechanisms, providing deep insights into the intricate dynamics involved in combustion, atmospheric, and organic chemistry.

ChemPotPy: A Python Library for Analytic Representations of Potential Energy Surfaces and Diabatic Potential Energy Matrices

Constructing analytic representations of global and semiglobal potential energy surfaces is difficult and can be laborious, and it is even harder when one needs coupled potential energy surfaces and their electronically nonadiabatic couplings. When accomplished, however, the resulting potential functions are a valuable resource. To facilitate the convenient use of potentials that have been developed, we provide a collection of existing surfaces in a library with consistent units and formats. A potential energy surface library of this type, namely PotLib, was built more than 20 years ago. However, that library only provided pristine Fortran subroutines for each potential energy surface, and therefore, it is not as user-friendly as would be desirable. Here, we report the creation of ChemPotPy, a CHEMical library of POTential energy surfaces in PYthon. ChemPotPy is a user-friendly library for analytic representation of single-state and multistate potential energy surfaces and couplings. A given entry in the library contains an analytic potential energy function or analytic functions for a set of coupled potential energy surfaces, and depending on the case, it may also include analytic or numerical gradients, nonadiabatic coupling vectors, and/or diabatic potential energy matrices and their gradients. Only three inputs, namely, the chemical formula of the system, the name of the potential energy surface or surface set, and the Cartesian geometry, are required. ChemPotPy uses the same units for input and output quantities of all surfaces and surface sets to facilitate general interfaces with the dynamics programs. The initial version of the library contains 338 entries, and we anticipate that more will be added in the future.

Theoretical Insights into the Dynamics of Gas-Phase Bimolecular Reactions with Submerged Barriers

Much attention has been paid to the dynamics of both activated gas-phase bimolecular reactions, which feature monotonically increasing integral cross sections and Arrhenius kinetics, and their barrierless capture counterparts, which manifest monotonically decreasing integral cross sections and negative temperature dependence of the rate coefficients. In this Perspective, we focus on the dynamics of gas-phase bimolecular reactions with submerged barriers, which often involve radicals or ions and are prevalent in combustion, atmospheric chemistry, astrochemistry, and plasma chemistry. The temperature dependence of the rate coefficients for such reactions is often non-Arrhenius and complex, and the corresponding dynamics may also be quite different from those with significant barriers or those completely dominated by capture. Recent experimental and theoretical studies of such reactions, particularly at relatively low temperatures or collision energies, have revealed interesting dynamical behaviors, which are discussed here. The new knowledge enriches our understanding of the dynamics of these unusual reactions.

Synthesis of Hidden Subgroup Quantum Algorithms and Quantum Chemical Dynamics

We describe a general formalism for quantum dynamics and show how this formalism subsumes several quantum algorithms, including the Deutsch, Deutsch-Jozsa, Bernstein-Vazirani, Simon, and Shor algorithms as well as the conventional approach to quantum dynamics based on tensor networks. The common framework exposes similarities among quantum algorithms and natural quantum phenomena: we illustrate this connection by showing how the correlated behavior of protons in water wire systems that are common in many biological and materials systems parallels the structure of Shor's algorithm.

ABC+D: A time-independent coupled-channel quantum dynamics program for elastic and ro-vibrational inelastic scattering between atoms and triatomic molecules in full dimensionality

We discuss the details of a time-independent quantum mechanical method and its implementation for full-dimensional non-reactive scattering between a closed-shell triatomic molecule and a closed-shell atom. By solving the time-independent Schrödinger equation within the coupled-channel framework using a log-derivative method, the state-to-state scattering matrix (S-matrix) can be determined for inelastic scattering involving both the rotational and vibrational modes of the molecule. Various approximations are also implemented. The ABC+D code provides an important platform for understanding an array of physical phenomena involving collisions between atoms and molecules.

A Machine Learning Approach for Rate Constants. III. Application to the Cl(2P) + CH4 → CH3 + HCl Reaction

The temperature dependence of the thermal rate constant for the reaction Cl(³P) + CH₄ → HCl + CH₃ is calculated using a Gaussian Process machine learning (ML) approach to train on and predict thermal rate constants over a large temperature range. Following procedures developed in two previous reports, we use a training data set of approximately 40 reaction/potential surface combinations, each of which is used to calculate the corresponding database of rate constant at approximately eight temperatures. For the current application, we train on the entire data set and then predict the temperature dependence of the title reaction employing a "split" data set for correction at low and high temperatures to capture both tunneling and recrossing. The results are an improvement on recent RPMD calculations compared to accurate quantum ones, using the same high-level ab initio potential energy surface. Both tunneling at low temperatures and significant recrossing at high temperatures are observed to influence the rate constants. The recrossing effects, which are not described by TST and even sophisticated tunneling corrections, do appear in experiment at temperatures above around 600 K. The ML results describe these effects and in fact merge at 600 K with RPMD results (which can describe recrossing), and both are close to experiment at the highest experimental temperatures. These results are in accord with a recent high-level experiment-theory study of this reaction.

Stereodynamics-Controlled Product Branching in the Nonadiabatic H + NaD → Na(3s, 3p) + HD Reaction at Low Temperatures

Nonadiabatic processes play an important role at energies near or higher than conical intersection of adiabatic potential energy surfaces in chemical reactions. In this work, dynamics of the nonadiabatic H + NaD reaction at low temperatures are studied by using the quantum wave packet method based on an improved L-shaped grid. The nonadiabatic H + NaD reaction has two exothermic reaction channels: Na(3s) + HD and Na(3p) + HD; the latter can only occur via nonadiabatic transition. The dynamics results show that the product branching of the H + NaD reaction at collision energies ranging from 20 to 80 cm^-1 is controlled by stereodynamics. The Na(3s) and Na(3p) reaction channels occur through collinear collision and side-on collision, respectively. When the collision energy is lower than 20 cm^-1, the resonance-mediated reaction mechanism is dominant in both the Na(3s) and Na(3p) reaction channels.

Toward a Comprehensive Understanding of Mode-Specific Dynamics of Polyatomic Reactions: A Full-Dimensional Quantum Dynamics Study of the H + NH3 Reaction

Mode specificity not only sheds light on reaction dynamics but also opens the door for chemical reaction control. This work reports a state-of-the-art full-dimensional quantum dynamics study on the prototypical hydrogen abstraction reaction of hydrogen with ammonia, which serves as a benchmark for advancing our fundamental understanding of polyatomic reaction dynamics. By taking advantage of the (3 + 1) Radau-Jacobi coordinates, the bond-specific probabilities are resolved with the reactant NH₃ initiated from either a non-degenerate or degenerate stretching vibrational state. The observed different atom-specific abstraction probabilities from individual states of the degenerate pair are rationalized in the local mode representation according to the different vibrational energy deposited in each N-H bond. It is verified that the three H atoms in NH₃ have the same abstraction probability as that from the degenerate pair and the linear combination of the degenerate pair gives the same reaction probability. In addition, the symmetric and asymmetric stretching modes of the reactant NH₃ have comparable efficacies on driving the reaction.

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.

High-fidelity first principles nonadiabaticity: diabatization, analytic representation of global diabatic potential energy matrices, and quantum dynamics

Nonadiabatic dynamics, which goes beyond the Born-Oppenheimer approximation, has increasingly been shown to play an important role in chemical processes, particularly those involving electronically excited states. Understanding multistate dynamics requires rigorous quantum characterization of both electronic and nuclear motion. However, such first principles treatments of multi-dimensional systems have so far been rather limited due to the lack of accurate coupled potential energy surfaces and difficulties associated with quantum dynamics. In this Perspective, we review recent advances in developing high-fidelity analytical diabatic potential energy matrices for quantum dynamical investigations of polyatomic uni- and bi-molecular nonadiabatic processes, by machine learning of high-level ab initio data. Special attention is paid to methods of diabatization, high fidelity construction of multi-state coupled potential energy surfaces and property surfaces, as well as quantum mechanical characterization of nonadiabatic nuclear dynamics. To illustrate the tremendous progress made by these new developments, several examples are discussed, in which direct comparison with quantum state resolved measurements led to either confirmation of the observation or sometimes reinterpretation of the experimental data. The insights gained in these prototypical systems greatly advance our understanding of nonadiabatic dynamics in chemical systems.

Benchmarking an improved statistical adiabatic channel model for competing inelastic and reactive processes

Inelastic collisions and elementary chemical reactions proceeding through the formation and subsequent decay of an intermediate collision complex, with an associated deep well on the potential energy surface, pose a challenge for accurate fully quantum mechanical approaches, such as the close-coupling method. In this study, we report on the theoretical prediction of temperature-dependent state-to-state rate coefficients for these complex-mode processes, using a statistical quantum method. This statistical adiabatic channel model is benchmarked by a direct comparison using accurate rate coefficients from the literature for a number of systems (H₂ + H⁺, HD + H⁺, SH⁺ + H, and CH⁺ + H) of interest in astrochemistry and astrophysics. For all of the systems considered, an error of less than factor 2 was found, at least for the dominant transitions and at low temperatures, which is sufficiently accurate for applications in the above mentioned disciplines.

Quantum Spin State Selectivity and Magnetic Tuning of Ultracold Chemical Reactions of Triplet Alkali-Metal Dimers with Alkali-Metal Atoms

We demonstrate that it is possible to efficiently control ultracold chemical reactions of alkali-metal atoms colliding with open-shell alkali-metal dimers in their metastable triplet states by choosing the internal hyperfine and rovibrational states of the reactants as well as by inducing magnetic Feshbach resonances with an external magnetic field. We base these conclusions on coupled-channel statistical calculations that include the effects of hyperfine contact and magnetic-field-induced Zeeman interactions on ultracold chemical reactions of hyperfine-resolved ground-state Na and the triplet NaLi(a^{3}Σ^{+}) producing singlet Na_{2}(^{1}Σ_{g}^{+}) and a Li atom. We find that the reaction rates are sensitive to the initial hyperfine states of the reactants. The chemical reaction of fully spin-polarized, high-spin states of rotationless NaLi(a^{3}Σ^{+},v=0,N=0) molecules with fully spin-polarized Na is suppressed by a factor of 10-100 compared to that of unpolarized reactants. We interpret these findings within the adiabatic state model, which treats the reaction as a sequence of nonadiabatic transitions between the initial nonreactive high-spin state and the final low-spin states of the reaction complex. In addition, we show that magnetic Feshbach resonances can similarly change reaction rate coefficients by several orders of magnitude. Some of these resonances are due to resonant trimer bound states dissociating to the N=2 rotational state of NaLi(a^{3}Σ^{+},v=0) and would thus exist in systems without hyperfine interactions.

A Time-Independent Quantum Approach to Ro-vibrationally Inelastic Scattering between Atoms and Triatomic Molecules

A full-dimensional time-independent quantum mechanical theory for ro-vibrationally inelastic scattering of triatomic molecules with atoms is formulated. The Jacobi-Radau coordinate system used in the calculation allows not only a near perfect description of the vibrational problem but also the adaptation of the exchange symmetry for A₂B type triatoms. The S-matrix elements are obtained by solving the close-coupling equations with contracted basis using the log-derivative method. This method is applied to the inelastic scattering of the water molecule by a chlorine atom, which sheds light on the energy gap law in energy transfer in atom-triatom collisions.

Modeling Spin-Crossover Dynamics

In this article, we review nonadiabatic molecular dynamics (NAMD) methods for modeling spin-crossover transitions. First, we discuss different representations of electronic states employed in the grid-based and direct NAMD simulations. The nature of interstate couplings in different representations is highlighted, with the main focus on nonadiabatic and spin-orbit couplings. Second, we describe three NAMD methods that have been used to simulate spin-crossover dynamics, including trajectory surface hopping, ab initio multiple spawning, and multiconfiguration time-dependent Hartree. Some aspects of employing different electronic structure methods to obtain information about potential energy surfaces and interstate couplings for NAMD simulations are also discussed. Third, representative applications of NAMD to spin crossovers in molecular systems of different sizes and complexities are highlighted. Finally, we pose several fundamental questions related to spin-dependent processes. These questions should be possible to address with future methodological developments in NAMD.

Dynamics of the isotope exchange reaction of D with H3 +, H2D+, and D2H. J Chem Phys 2021;154:084307. [PMID: 33639774 DOI: 10.1063/5.0038434]

We have measured the merged-beams rate coefficient for the titular isotope exchange reactions as a function of the relative collision energy in the range of ∼3 meV-10 eV. The results appear to scale with the number of available sites for deuteration. We have performed extensive theoretical calculations to characterize the zero-point energy corrected reaction path. Vibrationally adiabatic minimum energy paths were obtained using a combination of unrestricted quadratic configuration interaction of single and double excitations and internally contracted multireference configuration interaction calculations. The resulting barrier height, ranging from 68 meV to 89 meV, together with the various asymptotes that may be reached in the collision, was used in a classical over-the-barrier model. All competing endoergic reaction channels were taken into account using a flux reduction factor. This model reproduces all three experimental sets quite satisfactorily. In order to generate thermal rate coefficients down to 10 K, the internal excitation energy distribution of each H₃ ⁺ isotopologue is evaluated level by level using available line lists and accurate spectroscopic parameters. Tunneling is accounted for by a direct inclusion of the exact quantum tunneling probability in the evaluation of the cross section. We derive a thermal rate coefficient of <1×10^-12 cm³ s^-1 for temperatures below 44 K, 86 K, and 139 K for the reaction of D with H₃ ⁺, H₂D⁺, and D₂H⁺, respectively, with tunneling effects included. The derived thermal rate coefficients exceed the ring polymer molecular dynamics prediction of Bulut et al. [J. Phys. Chem. A 123, 8766 (2019)] at all temperatures.

Interaction-Asymptotic Region Decomposition Method for an Insertion Reaction: Application to the S(1D) + H2 Reaction

With adjusting principal axes hyperspherical (APH) coordinate in the interaction region, and the Jacobi coordinates in the asymptotic regions, an efficient multidomain interaction-asymptotic region decomposition (IARD) method has been developed to solve the "coordinate problem" in a product-state-resolved reactive scattering calculation using the quantum wave packet method. Although the APH coordinate treats with all three channels equally, and is efficient for describing the interaction region for some direct reactions, it is inefficient for describing the insertion-type reaction due to the singularity problem, such as the S(¹D) + H₂ reaction. To deal with this issue, in this work, the channel-dependent Delves hyperspherical (DH) coordinate is proposed to describe the interaction region using the IARD method. The proposed DH-IARD method was applied to calculate the product-state-resolved reaction probabilities of the H + HD reaction, and the differential and integral cross sections of the typical insertion reaction S(¹D) + H₂. It is found that the new DH-IARD method is much more efficient than the previous APH-IARD method for dealing with insertion reactions. The partial wave resonance structures were observed in the integral cross section. It is found that at a low collision energy, the position of the initial wave packet has to be put far away. Otherwise, the partial wave resonance structures could not be correctly reproduced due to the reef well arising with a large total angular momentum J.

Time-dependent quantum mechanical wave packet dynamics

Starting from a model study of the collinear (H, H₂) exchange reaction in 1959, the time-dependent quantum mechanical wave packet (TDQMWP) method has come a long way in dealing with systems as large as Cl + CH₄. The fast Fourier transform method for evaluating the second order spatial derivative of the wave function and split-operator method or Chebyshev polynomial expansion for determining the time evolution of the wave function for the system have made the approach highly accurate from a practical point of view. The TDQMWP methodology has been able to predict state-to-state differential and integral reaction cross sections accurately, in agreement with available experimental results for three dimensional (H, H₂) collisions, and identify reactive scattering resonances too. It has become a practical computational tool in predicting the observables for many A + BC exchange reactions in three dimensions and a number of larger systems. It is equally amenable to determining the bound and quasi-bound states for a variety of molecular systems. Just as it is able to deal with dissociative processes (without involving basis set expansion), it is able to deal with multi-mode nonadiabatic dynamics in multiple electronic states with equal ease. We present an overview of the method and its strength and limitations, citing examples largely from our own research groups.

Experimental study of the proton-transfer reaction C + H2+ → CH+ + H and its isotopic variant (D2+)

We report absolute integral cross section (ICS) measurements using a dual-source merged-fast-beams apparatus to study the titular reactions over the relative translational energy range of Er ∼ 0.01-10 eV. We used photodetachment of C- to produce a pure beam of atomic C in the ground electronic 3P term, with statistically populated fine-structure levels. The H2+ and D2+ were formed in an electron impact ionization source, with well known vibrational and rotational distributions. The experimental work is complemented by a theoretical study of the CH2+ electronic system in the reactant and product channels, which helps to clarify the possible reaction mechanisms underlying the ICS measurements. Our measurements provide evidence that the reactions are barrierless and exoergic. They also indicate the apparent absence of an intermolecular isotope effect, to within the total experimental uncertainties. Capture models, taking into account either the charge-induced dipole interaction potential or the combined charge-quadrupole and charge-induced dipole interaction potentials, produce reaction cross sections that lie a factor of ∼4 above the experimental results. Based on our theoretical study, we hypothesize that the reaction is most likely to proceed adiabatically through the 14A' and 14A'' states of CH2+via the reaction C(3P) + H2+(2Σ+g) → CH+(3Π) + H(2S). We also hypothesize that at low collision energies only H2+(v ≤ 2) and D2+(v ≤ 3) contribute to the titular reactions, due to the onset of dissociative charge transfer for higher vibrational v levels. Incorporating these assumptions into the capture models brings them into better agreement with the experimental results. Still, for energies ⪅0.1 eV where capture models are most relevant, the modified charge-induced dipole model yields reaction cross sections with an incorrect energy dependence and lying ∼10% below the experimental results. The capture cross section obtained from the combined charge-quadrupole and charge-induced dipole model better matches the measured energy dependence but lies ∼30-50% above the experimental results. These findings provide important guidance for future quasiclassical trajectory and quantum mechanical treatments of this reaction.

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.

Computation and analysis of bound vibrational spectra of the neon tetramer using row orthonormal hyperspherical coordinates

This paper presents the first implementation of the row-orthonormal hyperspherical coordinate formalism for the computation of the vibrational spectrum of a tetratomic system. The wavefunction of Ne₄ is expanded on a large basis set of hyperspherical harmonics generated numerically. This method not only provides spectra with reasonable accuracy, but also gives physical insight into the vibrational dynamics of the system. The characteristics of the spectra are related to the symmetry and localization of the wavefunction in configuration space.

Calculations of quantum tunnelling rates for muonium reactions with methane, ethane and propane

Thermal rate constants for Mu + CH₄, Mu + C₂H₆ and Mu + C₃H₈ and their equivalent reactions with H were evaluated with ab initio instanton rate theory. The potential-energy surfaces are fitted using Gaussian process regression to high-level electronic-structure calculations evaluated around the tunnelling pathway. This method was able to successfully reproduce various experimental measurements for the rate constant of these reactions. However, it was not able to reproduce the faster-than-expected rate of Mu + C₃H₈ at 300 K reported by Fleming et al. [Phys. Chem. Chem. Phys., 2015, 17, 19901 and Phys. Chem. Chem. Phys., 2020, 22, 6326]. Analysis of our results indicates that the kinetic isotope effect at this temperature is not significantly influenced by quantum tunnelling. We consider many possible factors for the discrepancy between theory and experiment but conclude that in each case, the instanton approximation is unlikely to be the cause of the error. This is in part based on the good agreement we find between the instanton predictions and new multiconfigurational time-dependent Hartree (MCTDH) calculations for Mu + CH₄ using the same potential-energy surface. Further experiments will therefore be needed to resolve this issue.

A General User-Friendly Tool for Kinetic Calculations of Multi-Step Reactions within the Virtual Multifrequency Spectrometer Project

We discuss the implementation of a computer program for accurate calculation of the kinetics of chemical reactions integrated in the user-friendly, multi-purpose Virtual Multifrequency Spectrometer tool. The program is based on the ab initio modeling of the involved molecular species, the adoption of transition-state theory for each elementary step of the reaction, and the use of a master-equation approach accounting for the complete reaction scheme. Some features of the software are illustrated through examples including the interconversion reaction of hydroxyacetone and 2-hydroxypropanal and the production of HCN and HNC from vinyl cyanide.

A Machine Learning Approach for Prediction of Rate Constants

We report a machine learning approach to train and predict bimolecular thermal rate constants over a large temperature range. The approach uses Gaussian process (GP) regression to evaluate the difference between accurate quantum results and Eckart-corrected conventional transition state theory, mostly for collinear reactions. Training is done on a database of rate constants for 13 reaction/potential surface combinations, and testing is performed on a set of 39 reaction/potential surface combinations. Averaged over all test reactions, the GP method is within 80% of the accurate answer, whereas transition state theory (TST) is only within 330% and Eckart-corrected TST (ECK) is within 110%. In the tunneling region, GP is generally (with a few exceptions) more accurate and sometimes much more accurate. In the high-temperature recrossing region, GP is significantly more accurate than either TST or ECK, neither of which addresses the possibility of recrossing. The GP predictions for the 3D reactions O(³P) + H₂, OH + H₂, O(³P) + CH₄, and H + CH₄, for which accurate quantum results are available, provide further encouragement to the machine learning approach.

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

Integral cross sections and rate constants for the prototypical chemical reactions of the fluorine atom with molecular hydrogen and deuterium have been calculated over a wide interval of collision energy and temperature ranging from the sub-thermal (50 K) down to the ultra-cold regimes (0.5 mK). Rigorous close coupling time-independent quantum reactive scattering calculations have been carried out on two potential energy surfaces, differing only at long-range in the reactants' channel. The results show that tunnel, resonance and virtual state effects enhance under-barrier reactivity giving rise to pronounced deviations from the Arrhenius law as temperature is lowered. Within the ultra-cold domain (below 1 mK), the reactivity is governed by virtual state effects and by tunneling through the reaction barrier; in the cold regime (1 mK–1 K), the shape resonances in the entrance channel of the potential energy surface make the quantum tunneling contribution larger so enhancing cross sections and rate constants by about one order of magnitude; at higher temperatures (above 10 K), the tunneling pathway enhanced by the constructive interference between two Feshbach resonances trapped in the reaction exit channel competes with the thermally activated mechanism, as the energy gets closer to the reaction barrier height. The results show that at low temperatures cross sections and rate constants are extremely sensitive to small changes in the long-range intermolecular interaction in the entrance channel of the potential energy surface, as well as to isotopic substitution.

Conical intersection-regulated intermediates in bimolecular reactions: Insights from C(1D) + HD dynamics

The importance of conical intersections (CIs) in electronically nonadiabatic processes is well known, but their influence on adiabatic dynamics has been underestimated. Here, through combined experimental and theoretical studies, we show that CIs induce a barrier and regulate conversion from a precursor metastable intermediate (CI-R) to a deep well one. This results in bond-selective activation, influencing the adiabatic dynamics markedly in the C(1D) + HD reaction. Theory is validated by experiment; quantum dynamics calculations on highly accurate ab initio potential energy surfaces yield rate coefficients and product branching ratios in excellent agreement with the experiment. Quasi-classical trajectory calculations reveal that the CI-R intermediate leads to unusual reaction mechanisms (designated as C─H activation complex conversion and cyclic complex), which are responsible for large branching ratios. We also reveal that CI-R intermediates exist in other reactive systems, and the dynamical effects uncovered here may have general significance.

Kinetic isotope effects in the water forming reaction H2/D2 + OH from rigorous close-coupling quantum dynamics simulations

Rigorous quantum dynamics simulations of thermal rate constants and kinetic isotope effects for the water-forming H₂/D₂ + OH reaction are presented, which show increased tunneling below 300 K and can serve as benchmarks for approximate methods.

Modified Gaussian Wave Packet Method for Calculating Initial State Wave Functions in Photodissociation

A modified Gaussian wave packet relaxation method is proposed to calculate the ground state wave function using an expansion of frozen Gaussian wave packets. This new procedure consists of two steps. In the first step, a multidimensional Gaussian product placed at the ground state equilibrium geometry is propagated in imaginary time. The relaxation optimizes the widths of the one-dimensional Gaussians. In the second step, additional Gaussian wave packets with the same widths are placed near the equilibrium geometry, and the corresponding expansion coefficients are optimized using the same relaxation method. This new algorithm is tested in photodissociation of NOCl and NH₃, and the results show good agreement with the exact results in the energy, wave function, and absorption spectrum. In particular, the highly structured absorption spectrum of NH₃ is reproduced, underscoring the accuracy of both the initial wave packet and the excited state propagation.

Low-temperature chemistry using the R-matrix method

Techniques for producing cold and ultracold molecules are enabling the study of chemical reactions and scattering at the quantum scattering limit, with only a few partial waves contributing to the incident channel, leading to the observation and even full control of state-to-state collisions in this regime. A new R-matrix formalism is presented for tackling problems involving low- and ultra-low energy collisions. This general formalism is particularly appropriate for slow collisions occurring on potential energy surfaces with deep wells. The many resonance states make such systems hard to treat theoretically but offer the best prospects for novel physics: resonances are already being widely used to control diatomic systems and should provide the route to steering ultracold reactions. Our R-matrix-based formalism builds on the progress made in variational calculations of molecular spectra by using these methods to provide wavefunctions for the whole system at short internuclear distances, (a regime known as the inner region). These wavefunctions are used to construct collision energy-dependent R-matrices which can then be propagated to give cross sections at each collision energy. The method is formulated for ultracold collision systems with differing numbers of atoms.

Active vs. spectator modes in nonadiabatic photodissociation dynamics of the hydroxymethyl radical via the 22A(3s) Rydberg state

The choice of the active degrees of freedom (DOFs) is a pivotal issue in a reduced-dimensional model of quantum dynamics when a full-dimensional one is not feasible. Here, several five-dimensional (5D) models are used to investigate the nonadiabatic photodissociation dynamics of the hydroxymethyl (CH₂OH) radical, which possesses nine internal DOFs, in its lowest absorption band. A normal-mode based scheme is used to identify the active and spectator modes, and its predictions are confirmed by 5D quantum dynamical calculations. Our results underscore the important role of the CO stretching mode in the photodissociation dynamics of CH₂OH, originating from the photo-induced promotion of an electron from the half-occupied π*_CO antibonding orbital to a carbon Rydberg orbital.

Angular momentum-scattering angle quantum correlation: a generalized deflection function

A natural generalization of the classical deflection function, the functional dependence of the deflection angle on the angular momentum (or the impact parameter), is the joint probability density function of these two quantities, revealing the correlation between them. It provides, at a glance, detailed information about the reaction mechanisms and how changes in the impact parameter affect the product angular distribution. It is also useful to predict the presence of quantum phenomena such as interference. However, the classical angular momentum-scattering angle correlation function has a limited use whenever quantum effects become important. Rigorously speaking, there is not a quantum equivalent of the classical joint distribution, as the differential cross section depends on the coherences between the different values of J caused by the cross terms in the expansion of partial waves. In this article, we present a simple method to calculate a quantum analog of this correlation, a generalized deflection function that can shed light onto the reaction mechanism using just quantum mechanical results. Our results show that there is a very good agreement between the quantum and classical correlation functions as long as quantum effects are not all relevant. When this is not the case, it will also be shown that the quantum correlation function is most useful to observe the extent of quantum effects such as interference among different reaction mechanisms.