Quantum wavepacket method for state-to-state reactive cross sections

Reaction Path-Resolved Quantum Transition State Framework Using Hyperspherical (APH) Coordinates: The Geometric Phase Effects in the H + H2 Reaction

The quantum transition state framework was developed to calculate the reaction path-resolved scattering matrix for atom-diatom reactions in hyperspherical (APH) coordinates. This approach allows for simply and directly calculating the reaction path-resolved scattering matrix, especially when the encircling reaction path is negligible. It could be used to determine the reactivities of specific pathways in a chemical reaction, providing insights into phenomena such as geometric phase effects. To validate our method, detailed calculations for the H + H2 reaction at J = 0, 1, 2 were carried out, and the results were compared with those from our previous theoretical models including the geometric phase effects. Analysis of the path-resolved reaction probabilities reveals that encircling paths play a very minor role in this reaction, even at energies exceeding the conical intersection minimum. Furthermore, the analysis suggests that the geometric phase effects arising from the interferences between different reactive pathways could become more significant for highly rotationally excited reactants and products.

The Interaction-Asymptotic Region Decomposition Method in Jacobi Coordinates: Triatomic Reactive Scatterings

The interaction-asymptotic region decomposition (IARD) technique has been proven to be a good solution to the long-standing coordinate problem in reactive scattering calculations. In this work, the IARD technique was further developed using Jacobi coordinates for the interaction region, instead of the previously used hyperspherical coordinates. Although the Jacobi coordinate may not be as optimal as the hyperspherical coordinates for describing the interaction region in reactive scattering processes, it has simpler kinetic operators and provides a more physically intuitive picture. By developing an intermediate interpolation method, which could efficiently transform the overlapped wave functions from the asymptotic regions to the interaction region, the new implementation of the IARD technique for triatomic reactive scatterings is similarly efficient and accurate. The differential cross sections of the H+H2, and product state-resolved reaction probabilities of the F+HD and16O+36O2 reactions, which involve products of extremely low translational energy and are challenging for a single coordinate-based method, were calculated as numerical examples to show the ability of the new method.

An improved method for reactive scatterings in ultra-cold conditions using the time-dependent approach

A new efficient method for considering the long-range effect of reactive scattering processes in ultra-cold conditions has been developed using the time-dependent quantum wave packet theory, where the initial wave packet could be placed at a position near the interaction region. This is in contrast to previous methods, where the initial wave packet has to be placed far from the interaction region. The new method reduces the numerical effort significantly. Typical reactions, such as S(1D) + H2, D+ + H2, and 7Li + 7Li2 (v0 = 1, j0 = 0), under cold or ultra-cold conditions, are used to demonstrate the numerical efficiency of the new method.

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.

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.

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.

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.

Effect of Reagent Vibration and Rotation on the State-to-State Dynamics of the Hydrogen Exchange Reaction, H + H2 → H2 + H

State-to-state dynamics of the benchmark hydrogen exchange reaction H + H₂ (v = 0-4, j = 0-3) → H₂ (v', j') + H is investigated with the aid of the real wave packet approach of Gray and Balint-Kurti (J. Chem. Phys. 1998, 108, 950-962) and electronic ground BKMP2 potential energy surface of Boothroyd et al. (J. Chem. Phys. 1996, 104, 7139-7152). Initial state-selected and product state-resolved reaction probabilities, integral cross section, and product diatom vibrational and rotational level populations at a few collision energies are reported to elucidate the energy disposal mechanism. State-specific thermal rate constants are also calculated and compared with the available literature results. Coriolis coupling terms of the nuclear Hamiltonian are included, and calculations are parallelized over the helicity quantum number, Ω'. Attempts are made, in particular, to study the effect of reagent vibrational and rotational excitations on the dynamical attributes. It is found that the calculations become computationally expensive with reagent vibrational and rotational excitation. Reagent vibrational excitation is found to enhance the reactivity and has significant impact on the energy disposal to the vibrational and rotational degrees of freedom of the product. The interplay of reagent translational and vibrational energy on the product vibrational distribution unfolds an important aspect of the energy disposal mechanism. The effect of reagent rotation on the state-to-state dynamics is found not to be very significant, and the weak effect turns out to be specific to v'.

High Order Split Operators for the Time-Dependent Wavepacket method of Triatomic Reactive Scattering in Hyperspherical Coordinates

Since the introduction of a series of methods for solving the time-dependent Schrödinger equation (TDSE) in the 80s of the last centry, such as the Fourier transform, the split operator (SO), the Chebyshev polynomial propagator, and complex absorbing potential, investigation of the molecular dynamics within quantum mechanics principle have become popular. In this paper, the application of the time-dependent wave packet (TDWP) method using high-order SO propagators in hyperspherical coordinates for solving triatomic reactive scattering was investigated. The fast sine transform was applied to calculate the derivatives of the wave function of the radial degree of freedom. These high-order SO propagators are examined in different forms, i.e., TVT (Kinetic–Potential–Kinetic) and VTV (Potential–Kinetic–Potential) forms with three typical triatomic reactions, H + H 2, O + O 2 and F + HD. A little difference has been observed among the performances of high-order SO propagators in the TVT and VTV representations in the hyperspherical coordinate. For obtaining total reaction probabilities with 1% error, some of the S class high-order SO propagators, which have symmetric forms, are more efficient than second order SO for reactions involving long lived intermediate states. High order SO propagators are very efficient for obtaining total reaction probabilities.

Recent advances in quantum scattering calculations on polyatomic bimolecular reactions

This review surveys quantum scattering calculations on chemical reactions of polyatomic molecules in the gas phase published in the last ten years.

Recent Advances in Quantum Dynamics of Bimolecular Reactions

In this review, we survey the latest advances in theoretical understanding of bimolecular reaction dynamics in the past decade. The remarkable recent progress in this field has been driven by more accurate and efficient ab initio electronic structure theory, effective potential-energy surface fitting techniques, and novel quantum scattering algorithms. Quantum mechanical characterization of bimolecular reactions continues to uncover interesting dynamical phenomena in atom-diatom reactions and beyond, reaching an unprecedented level of sophistication. In tandem with experimental explorations, these theoretical developments have greatly advanced our understanding of key issues in reaction dynamics, such as microscopic reaction mechanisms, mode specificity, product energy disposal, influence of reactive resonances, and nonadiabatic effects.

State-to-State Quantum Dynamics of Reactions O((3)P) + HD (v = 0-1, j = 0) → OH+D and OD+H: Reaction Mechanism and Vibrational Excitation

Time-dependent quantum wave packet dynamics calculations have been performed in order to characterize the dynamics and mechanism of O((3)P) + HD (v = 0-1, j = 0) → OH+D and OD+H reactive collisions using the adiabatic potential energy surface by Rogers et al. [J. Phys. Chem. A 2000, 104, 2308] Special attention has been paid to the calculations and discussion of the state resolved integral and differential cross sections and the product state distributions. In addition, the intramolecular isotopic branching ratio has been determined. The results revealed that the OD + H is the favored product channel and the product OH has the same quantum number v as the reactant HD. For low collision energy, the product angular distributions concentrate in the backward region being consistent with a rebounding mechanism. In the case of higher collision energy, the stripping collisions with larger impact parameters tend to produce sideways and forward scatterings, especially for the HD vibrationally excited state. The cross section and intramolecular isotopic branching ratio are in agreement with the previous theoretical results. A cartoon depiction collision model is built and works well for our calculation results.

Quantum, Statistical, and Quasiclassical Trajectory Studies For the Ne + HeH(+) → NeH(+) + He Reaction on the Ground Electronic State

Real wave packet, statistical quantum, and quasiclassical trajectory methods were employed to study the dynamics of Ne + HeH(+)(v0,j0) → He + NeH(+) reaction on an ab initio potential energy surface [J. Phys. Chem. A 2013, 117, 13070-13078]. Quantum and statistical quantum calculations were performed within the centrifugal sudden (CS) approximation as well as including the Coriolis coupling (CC). Dense oscillatory structures of the quantum reaction probabilities and fair agreement between quantum and statistical cross sections suggest a complex forming mechanism for the reaction. No significant differences between cross sections obtained within the CS and CC approaches are observed. Quasiclassical trajectory results give an excellent average description of the quantum CC results. At low collision energies, there is a substantial decrease in reactivity for the reaction upon rovibrational excitation. Initial state selected rate constants for the title reaction are calculated between 20 and 1000 K, and the calculated value at 300 K agrees quite well with the available experimental result. Reaction cross sections and rate constants are also compared with those calculated via the Langevin capture model for exothermic reactions.

Collision time of a triatomic chemical reaction A + BC

The collision time is an important quantity of an elementary chemical reaction and describes the speed of the collision process in a collision reaction. In this study, we present a generalized method to calculate the collision time of a triatomic reaction in which the collision time is defined by the sum of the incoming time, the intermediate complex time, and the outgoing time. Two variables including the total distance Rtotal and Ravg, the average value of Rtotal over time, are used to compute the three components of the collision time. We compute three triatomic reactions including Ca + HCl → CaCl + H, O + HCl → OH + Cl/OCl + H, and O + HF → OH + F at different collision energies and initial diatomic vibrational levels using the quasi-classical trajectory method to confirm that the method could be reliable and reasonable. The time evolutions of Rtotal could efficiently classify the direct and indirect reactive mechanisms and reveal a distinct discrepancy of the two mechanisms. As the collision energy and initial diatomic vibrational level increase, the percentage of direct reaction trajectories increases. At the same time, the average and maximal values of collision time decrease. Comparing the maximal collision time and the reactive probability distributions of the products, it could be found that most reactive trajectories’ collision time is less than 2 ps. Moreover, the present calculations indicate that the method could be applicable to estimate the lifetime of the intermediate complex for the reaction systems with deep potential wells and the collision time of the reactions with a direct abstraction mechanism.

Wave packet and statistical quantum calculations for the He + NeH⁺ → HeH⁺ + Ne reaction on the ground electronic state

A real wave packet based time-dependent method and a statistical quantum method have been used to study the He + NeH(+) (v, j) reaction with the reactant in various ro-vibrational states, on a recently calculated ab initio ground state potential energy surface. Both the wave packet and statistical quantum calculations were carried out within the centrifugal sudden approximation as well as using the exact Hamiltonian. Quantum reaction probabilities exhibit dense oscillatory pattern for smaller total angular momentum values, which is a signature of resonances in a complex forming mechanism for the title reaction. Significant differences, found between exact and approximate quantum reaction cross sections, highlight the importance of inclusion of Coriolis coupling in the calculations. Statistical results are in fairly good agreement with the exact quantum results, for ground ro-vibrational states of the reactant. Vibrational excitation greatly enhances the reaction cross sections, whereas rotational excitation has relatively small effect on the reaction. The nature of the reaction cross section curves is dependent on the initial vibrational state of the reactant and is typical of a late barrier type potential energy profile.