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.

Spiers Memorial Lecture: New directions in molecular scattering

The field of molecular scattering is reviewed as it pertains to gas-gas as well as gas-surface chemical reaction dynamics. We emphasize the importance of collaboration of experiment and theory, from which new directions of research are being pursued on increasingly complex problems. We review both experimental and theoretical advances that provide the modern toolbox available to molecular-scattering studies. We distinguish between two classes of work. The first involves simple systems and uses experiment to validate theory so that from the validated theory, one may learn far more than could ever be measured in the laboratory. The second class involves problems of great complexity that would be difficult or impossible to understand without a partnership of experiment and theory. Key topics covered in this review include crossed-beams reactive scattering and scattering at extremely low energies, where quantum effects dominate. They also include scattering from surfaces, reactive scattering and kinetics at surfaces, and scattering work done at liquid surfaces. The review closes with thoughts on future promising directions of research.

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.

Globally Accurate Gaussian Process Potential Energy Surface and Quantum Dynamics Studies on the Li(2S) + Na2 → LiNa + Na Reaction at Low Collision Energies

The LiNa2 reactive system has recently received great attention in the experimental study of ultracold chemical reactions, but the corresponding theoretical calculations have not been carried out. Here, we report the first globally accurate ground-state LiNa2 potential energy surface (PES) using a Gaussian process model based on only 1776 actively selected high-level ab initio training points. The constructed PES had high precision and strong generalization capability. On the new PES, the quantum dynamics calculations on the Li(2S) + Na2(v = 0, j = 0) → LiNa + Na reaction were carried out in the 0.001–0.01 eV collision energy range using an improved time-dependent wave packet method. The calculated results indicate that this reaction is dominated by a complex-forming mechanism at low collision energies. The presented dynamics data provide guidance for experimental research, and the newly constructed PES could be further used for ultracold reaction dynamics calculations on this reactive system.

An accurate many-body expansion potential energy surface for AlH2 (22A') and quantum dynamics in Al(3P) + H2 (v0 = 0-3, j0 = 0, 2, 4, 6) collisions

An accurate potential energy surface is constructed for the excited state of AlH₂ by fitting extensive ab initio points calculated at the multi-reference configuration interaction level based on aug-cc-pV(Q+d)Z and aug-cc-pV(5+d)Z basis sets. All the calculated energies are corrected via the many-body expansion method and extrapolated to the complete basis set limit. The various topographic features of the new potential energy surface are investigated to demonstrate the correct behavior of Al(³P) + H₂(X¹Σ_g⁺) and AlH(a³Π) + H(²S) dissociation limits. By employing the time-dependent wave packet approach, the integral scattering cross-sections obtained from the Coriolis coupling calculation and the centrifugal sudden approximation, respectively, are compared in detail and show that the former has a higher effect on the reaction. Moreover, the thermal rate constants for Al(³P) + H₂ (v₀ = 0-3, j₀ = 0, 2, 4, 6) in the temperature range of 0-5000 K are calculated, thereby providing insights into the influence of ro-vibrational quantum numbers on the thermal rate constants.

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.

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.

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.

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.

A new potential energy surface for the ground electronic state of the LiH2 system, and dynamics studies on the H((2)S) + LiH(X(1)Σ(+)) → Li((2)S) + H2(X(1)Σg(+)) reaction

A new global potential energy surface (PES) is obtained for the ground electronic state of the LiH2 system based on high-level energies. The energy points are calculated at the multireference configuration interaction level with aug-cc-pVXZ (X = Q, 5) basis sets, and these energies are extrapolated to the complete basis set limit. The neural network method and hierarchical construction scheme are applied in the fitting process and the root mean square error of the fitting result is very small (0.004 eV). The dissociation energies and equilibrium distances for LiH(X(1)Σ(+)) and H2(X(1)Σg(+)) obtained from the new PES are in good agreement with the experimental data. On the new PES, time-dependent wave packet studies for the H((2)S) + LiH(X(1)Σ(+)) → Li((2)S) + H2(X(1)Σg(+)) reaction have been carried out. In this reaction, no threshold is found due to the absence of an energy barrier on the minimum energy path. The calculated integral cross sections are high at low collision energy and will decrease with the increase of the collision energy. The product molecule H2 tends to be forward scattering due to direct reactive collisions, which becomes more evident at higher collision energies.

State-resolved time-dependent wave packet and quasiclassical trajectory studies of the adiabatic reaction S(3P) + HD on the (1(3)A″) state

Time-dependent wave packet (TDWP) and quasiclassical trajectory (QCT) calculations have been carried out for the reaction S(3P) + HD(X1Σg+) at the lowest 13A″ state with both rotational and vibrational excitations of reactant HD. The calculated integral cross sections from QCT agree fairly well with the TDWP calculations. The reaction probability results from TDWP show that the reaction displays a strong tendency to the SD channel. When the reactant HD is vibrationally excited, both channels are promoted apparently. The vibration of the HD bond tends to reduce the difference of reactivity between the two channels. The detailed state-to-state differential cross sections (DCSs) are calculated. These distributions show some significant characters of the barrier-type reactions. At the same time, the scattering width of product SD has a certain relationship with its rotation excitation. For the vector properties, P(θr), P(r), and P(θr,r) distributions are calculated by QCT, and the increased collision energy weakens the rotational polarization of the SD molecule.

EXPONENTIAL PROPAGATORS (INTEGRATORS) FOR THE TIME-DEPENDENT SCHRÖDINGER EQUATION

The time-dependent Schrödinger Equation (TDSE) is a parabolic partial differential equation (PDE) comparable to a diffusion equation but with imaginary time. Due to its first order time derivative, exponential integrators or propagators are natural methods to describe evolution in time of the TDSE, both for time-independent and time-dependent potentials. Two splitting methods based on Fer and/or Magnus expansions allow for developing unitary factorizations of exponentials with different accuracies in the time step △t. The unitary factorization of exponentials to high order accuracy depends on commutators of kinetic energy operators with potentials. Fourth-order accuracy propagators can involve negative or complex time steps, or real time steps only but with gradients of potentials, i.e. forces. Extending the propagators of TDSE's to imaginary time allows to also apply these methods to classical many-body dynamics, and quantum statistical mechanics of molecular systems.

INVESTIGATION OF THE EXCHANGE REACTION H + H′S → HS + H′ ON THE 1A′ STATE POTENTIAL ENERGY SURFACE

The quantum time dependent wave packet (TDWP) and quasiclassical trajectory (QCT) calculations were carried out to study the exchange reaction H(2S) + H′S(2Π) → HS(2Π) + H′(2S) on the 1A′ potential energy surface (PES). The integral cross sections of the H + H′S (v = j = 0) → HS + H′ reaction calculated by the two methods were presented. The results reveal that the integral cross sections (ICS) decrease with the collision energy increasing. The result of the QCT calculations is reasonably consistent with the time-dependent wave packet. Moreover, the differential cross sections (DCS) were calculated by the QCT method at the four different collision energies, which display a forward–backward symmetry. A long-lifetime H2S intermediate complex of the exchange reaction was found according to the trajectories. In the stereodynamics investigation, the polar and dihedral angle distribution functions were calculated, which have the distinct oscillations. The oscillations could be attributed to the deep well on the 1A′ PES. However, based on the polar-angle and dihedral angle distribution functions, it could be predicted that the main product rotational angular momentum preferentially point to the positive or negative direction of y-axes.