A reactant-coordinate-based wave packet method for full-dimensional state-to-state quantum dynamics of tetra-atomic reactions: Application to both the abstraction and exchange channels in the H + H2O reaction

Differential Cross-Sections for the Vibrationally Excited H + HOD(vOH = 1-4) → H2 + OD Reactions

Using the time-dependent wave-packet approach, we calculate the first fully converged state-to-state differential cross-sections for the H + HOD(vOH = 1-4) → H2 + OD reactions on a highly accurate neural network PES. It is found that, unlike the loss of memory effect observed in the product distributions for low vibrational excitation reactions, high initial OH vibrational excitation significantly influences not only the product vibrational distribution but also the angular distribution. Furthermore, for the H + HOD(vOH = 3,4) reactions, the total integral cross-sections maintain the pronounced oscillatory structures in the J = 0 probabilities at low collision energies, which originate from the prereactive van der Waals resonances. Notably, the product angular distributions exhibit forward-backward peaked behavior at these energies, akin to those observed in the complex-forming system with deep potential wells. This is attributed to the much longer lifetimes of these resonances compared to those of conventional transition state resonances.

Study on Quantum Dynamics of the Na + Na2 Exchanged Reaction and Lifetime Prediction of Na3 Complex Based on the Neural Network Potential Energy Surface

A high-precision global potential energy surface (PES) is constructed for the Na3 system based on high-level ab initio calculations and the fundamental invariant neural network (FI-NN) method. The root-mean-square error (RMSE) of the PES is 2.88 cm-1. The state-resolved quantum dynamics of the ground-state Na + Na2 (v = 0, j = 0) → Na2 (v', j') + Na reaction is studied using the time-dependent wave packet (TDWP) method on the new PES. Analysis of the relevant integral cross sections revealed a complicated energy-transfer mechanism during collisions. Similarly, the characteristics of the differential cross sections indicate that the complex-forming mechanism plays a dominant role in the reaction, providing conditions for a comprehensive exploration of the lifetimes of the complexes. Based on the Rice-Ramsperger-Kassel-Marcus (RRKM) theory, the calculated lifetime of the Na3 complex is approximately 3.9 ns.

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.

Nonadiabatic dynamics studies of the H(2S) + RbH(X1Σ+) reaction: based on new diabatic potential energy surfaces

The global diabatic potential energy surfaces (PESs) that correspond to the ground (1²A') and first excited states (2²A') of the RbH₂ system PES are constructed based on 17 786 ab initio points. The neural network method is used to fit the PESs and the topographic features of the new diabatic PESs are discussed in detail. Based on the newly constructed diabatic PESs, the dynamics calculations of the H(²S) + RbH(X¹Σ⁺) → Rb(5²S) + H₂(X¹Σ_g ⁺)/Rb(5²P) + H₂(X¹Σ_g ⁺) reactions are performed using the time-dependent wave packet method. The dynamics properties of these two channels such as the reaction probabilities, integral cross sections, and differential cross sections (DCSs) are calculated at state-to-state level of theory. The nonadiabatic effects are discussed in detail, and the results indicate that the adiabatic results are overestimated from the dynamics values. The DCSs of these two channels are forward biased, which indicates that the abstraction mechanism plays a dominant role in the reaction.

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.

Permutationally Invariant Polynomial Expansions with Unrestricted Complexity

A general strategy is presented for constructing and validating permutationally invariant polynomial (PIP) expansions for chemical systems of any stoichiometry. Demonstrations are made for three categories of gas-phase dynamics and kinetics: collisional energy-transfer trajectories for predicting pressure-dependent kinetics, three-body collisions for describing transient van der Waals adducts relevant to atmospheric chemistry, and nonthermal reactivity via quasiclassical trajectories. In total, 30 systems are considered with up to 15 atoms and 39 degrees of freedom. Permutational invariance is enforced in PIP expansions with as many as 13 million terms and 13 permutationally distinct atom types by taking advantage of petascale computational resources. The quality of the PIP expansions is demonstrated through the systematic convergence of in-sample and out-of-sample errors with respect to both the number of training data and the order of the expansion, and these errors are shown to predict errors in the dynamics for both reactive and nonreactive applications. The parallelized code distributed as part of this work enables the automation of PIP generation for complex systems with multiple channels and flexible user-defined symmetry constraints and for automatically removing unphysical unconnected terms from the basis set expansions, all of which are required for simulating complex reactive systems.

A new global potential energy surface of the SH2+(X4A'') system and quantum calculations for the S+ + H2(v = 0-3, j = 0) reaction

A new global potential energy surface (PES) for the ground state of the SH2+(X4A'') system is constructed using a permutation invariant polynomial neural network method. In ab initio calculations, the MRCI-F12 method with the AVTZ basis set is used. Furthermore, the dynamics calculations of the S+ + H2(v = 0-3, j = 0) → SH+ + H reaction are carried out based on the new PES. The reaction probabilities and integral cross sections are compared with available theoretical calculations. Present values are in general good agreement with the previous theoretical studies. However, some discrepancies can still be found due to different PESs used in the calculation. Furthermore, the vibrational energy of the reactant molecule can significantly enhance the reactivity compared to the translational energy. The differential cross sections indicated that the reaction mechanism is changed from the "head-on" rebound mechanism to the tripping mechanism with the increasing number of initial vibrational excitation state.

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.

Charge Transfer Processes for H + H2+ Reaction Employing Coupled 3D Wavepacket Approach on Beyond Born-Oppenheimer Based Ab Initio Constructed Diabatic Potential Energy Surfaces

The dynamics of the H + H₂⁺ reaction has been analyzed from the electronically first excited state of diabatic potential energy surfaces constructed by employing the Beyond Born-Oppenheimer theory [J. Chem. Phys. 2014, 141, 204306]. We have employed the coupled 3D time-dependent wavepacket formalism in hyperspherical coordinates for multisurface reactive scattering problems. To be specific, the charge transfer processes have been investigated extensively by calculating state-to-state as well as total reaction probabilities and integral cross sections, when the reaction process is initiated from the first excited electronic state (2¹A'). We have depicted the convergence profiles of reaction probabilities for the competing charge transfer processes, namely, reactive charge transfer (RCT) and nonreactive charge transfer (NRCT) processes for different total energies with respect to total angular momentum, J. Total and state-to-state integral cross sections are calculated as a function of total energy for the initial rovibrational state, namely, v = 0, j = 0 level of H₂⁺ (²Σ_g⁺) molecule and are compared with previous theoretical calculations. Finally, we have calculated temperature-dependent rate constants using our presently evaluated cross sections and compared their average with the experimentally measured one.

Dynamics studies of the H + HBr reaction: Based on a new potential energy surface

The initial state specific quantum wave packet dynamics studies of the H + HBr (v₀ = 0, j₀ = 0-2) reaction were performed using a new global potential energy surface (PES) of the ground state of the BrH₂ system for the collision energy ranging from 0.01 to 2.0 eV. The PES was constructed using the permutation invariant polynomial neural network method based on approximately 63 000 ab initio points, which were calculated by the multireference configuration interaction method with AVTZ and AVQZ basis sets. To improve the accuracy of the PES, Davidson's correction and spin-orbit coupling effects were considered in the ab initio calculation and the basis set was extrapolated to complete basis set limit. The new PES was compared with the previous ones and also the available experimental data, which suggests that the new PES is more accurate. The state-to-state quantum wave packet dynamics was carried out using the reactant-coordinate based approach. The reaction probabilities, integral and differential cross sections, rovibrational state distributions of product and rate constants, etc., were compared with the available theoretical and experimental studies. In general, the present work is in better agreement with the available experimental data. The quantum dynamics studies suggest that the rotational excitation of HBr has little effect on the reaction.

Fermi resonance controlled product branching in the H + HOD reaction

Excitation of the first overtone of bending mode results in a significant enhancement in the HD + OH channel due to the 1 : 2 Fermi resonance between the fundamental OD stretch and the first overtone of the bend.

A new potential energy surface of the OH2+ system and state-to-state quantum dynamics studies of the O+ + H2 reaction

A new global potential energy surface of the O⁺ + H₂ system was constructed with neural network method, using about 63000 ab initio points, which were calculated by employing the multi-reference configuration interaction method with aug-cc-pVTZ and aug-cc-pVQZ basis sets.

State-to-state mode specificity in H + DOH(νOH = 1) → HD + OH(ν2 = 0) reaction: vibrational non-adiabaticity or local-mode excitation?

State-of-the-art full dimensional state-to-state quantum dynamics reveal a startling observation in which the DOH(ν_OH = 1) molecule reacts with a H atom to produce a vibrationless OH product. This interesting observation is attributed to a small OD excited local-mode component in the reactant wavefunction.

A global potential energy surface and dynamics study of the Au+ + H2 → H + Au+H reaction

A global potential energy surface (PES) of the ground state of the Au⁺H₂ system was constructed using a neural network method with permutation invariant polynomials.