Quantum wavepacket method for state-to-state reactive cross sections in hyperspherical coordinates

Quantum Dynamics of O(3P)+HBr→OH+Br Reaction: Integral Cross Sections and Rate Constants and their Initial State Dependence

The integral cross sections (ICSs) and rate constants (RCs) for the title reaction are calculated by means of accurate quantum wave packet method employing the recent ab initio potential energy surface developed by Peterson and co-workers. Hundreds of partial wave contributions (up to J=150) are calculated explicitly taking Coriolis coupling into account. The ICSs are found to increase monotonically with energy and their energy dependences are smooth except in the energy range below the potential energy barrier, where several resonance peaks are observed. The temperature dependences of the RCs obey in general the Arrhenius law. These behaviors manifest the dominance of abstract mechanism. Initial rotational state (ji) excitations are found to greatly enhance the reactivity starting from ji=4, while exhibiting a complicated dependence for ji<4. From Boltzmann average of the initial state specified RCs (for ji=0-15) we obtained thermal RCs which are in reasonable agreement with available experimental results.

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.

New Diabatic Potential Energy Surfaces for the Li + H2 Reaction and Time-Dependent Quantum Wave Packet Studies

New global diabatic potential energy surfaces (DPESs) for the ground (12A') and first excited (22A') states for the Li + H2 system were developed, with more than 30,000 energy points at the IC-MRCI+Q level of theory, utilizing the aug-cc-pV5Z basis set for the H atoms and the cc-pCV5Z basis set for the Li atom, fitted by a single neural network (NN) with symmetry. Product state-resolved quantum dynamics calculations of the nonadiabatic reaction Li (2P) + H2 (X 1 ∑g+, v0 = 0, j0 = 0) → LiH (X 1∑+) + H(2S) were carried out using these new DPESs and also the previous HYLC-DPESs. The numerical results suggested that our newly constructed DPESs provided an accurate description of the LiH2 system.

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.

Dynamical calculations of O(3P) + OH(2Π) reaction on the CHIPR potential energy surface using the fully coupled time-dependent wave-packet approach in hyperspherical coordinates

We have carried out quantum dynamics calculations for the O + OH → H + O₂ reaction on the CHIPR [A. J. C. Varandas, J. Chem. Phys., 2013, 138, 134117] potential energy surface (PES) for ground state HO₂ using the fully coupled 3D time-dependent wavepacket formalism [S. Adhikari and A. J. C. Varandas, Comput. Phys. Commun., 2013, 184, 270] in hyperspherical coordinates. Reaction probabilities for J > 0 are calculated for different initial rotational states of the OH radical (v = 0; j = 0, 1). State-to-state as well as total integral cross sections and rate-coefficients are evaluated and compared with previous theoretical calculations and available experimental studies. Using the rate constant for the forward (hereinafter considered to be H + O₂ → O + OH) and backward (O + OH → H + O₂) reactions of this reactive system, the equilibrium constant of the reversible process [H + O₂ ⇌ O + OH] is calculated as a function of temperature and compared with previous experimental measurements.

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.

Glories, hidden rainbows and nearside-farside interference effects in the angular scattering of the state-to-state H + HD → H2 + D reaction

Yuan et al. [Nat. Chem., 2018, 10, 653] have reported state-of-the-art measurements of differential cross sections (DCSs) for the H + HD → H2 + D reaction, measuring for the first time fast oscillations in the small-angle forward region of the DCSs. We theoretically analyse the angular scattering dynamics in order to quantitatively understand the physical content of structure in the DCSs. We study the H + HD(vi = 0, ji = 0, mi = 0) → H2(vf = 0, jf = 0,1,2,3, mf = 0) + D reaction for the whole range of scattering angles from θR = 0° to θR = 180°, where v, j, m are the vibrational, rotational and helicity quantum numbers respectively for the initial and final states. The restriction to mf = 0 arises because states with mf ≠ 0 have DCSs that are identically zero in the forward (θR = 0°) and backward (θR = 180°) directions. We use accurate quantum scattering matrix elements computed by Yuan et al. at a translational energy of 1.35 eV for the BKMP2 potential energy surface. The following theoretical techniques are employed to analyse the DCSs: (a) full and nearside-farside (NF) partial wave series (PWS) and local angular momentum theory, including resummations of the full PWS up to third order. We also use window representations of the scattering matrix, which give rise to truncated PWS, (b) six asymptotic (semiclassical) small-angle glory theories and four N rainbow theories, (c) we introduce "CoroGlo" tests, which let us distinguish between glory and corona scattering at small angles for Legendre PWS, (d) the semiclassical optical model (SOM) of Herschbach is employed to understand structure in the DCSs at intermediate and large angles. Our conclusions are: (a) the small-angle peaks in the DCSs arise mainly from glory scattering. For the 000 → 020 transition, there is also a contribution from a broad, or hidden, N rainbow, (b) at larger angles, the fast oscillations in the DCSs arise from NF interference, (c) the N scattering in the fast oscillation region contains a hidden rainbow for the 000, 020, 030 cases. For the 000 → 020 transition, the rainbow extends up to θR ≈ 60°; for the 000 and 030 cases, the angular ranges containing a N rainbow are smaller, (d) at intermediate and backward angles, the slowly varying DCSs, which merge into slow oscillations, are explained by the SOM. Physically it shows this structure in a DCS arises from direct scattering and is a distorted mirror image of the corresponding probability versus total angular momentum quantum number plot.

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.

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.

Non-adiabatic dynamics studies of the K(4p2P) + H2(X1Σ) reaction based on new diabatic potential energy surfaces

Global diabatic potential energy surfaces (PESs) of the KH2 system corresponding to the ground (12A') and first excited (22A') states were constructed for the first time. In ab initio calculations, the MRCI-F12 method with AVTZ and def2-QZVP basis sets was adopted and 17 865 ab initio energy points were calculated. The mixing angle, which is used to obtain the diabatic energies, was calculated by the molecular properties of the transition dipole moment. The diabatic PESs were fitted individually by the permutation invariant polynomial neural network method and the topographical features of the diabatic PESs are discussed in detail. The non-adiabatic dynamics studies of the K(4p2P) + H2(v0 = 0, 1, j0 = 0) reaction were carried out using the APH method based on the new diabatic PESs. The collision reaction processes K(4p2P) + H2(v0 = 0, 1, j0 = 0) → H + KH and the quenching processes K(4p2P) + H2(v0 = 0, 1, j0 = 0) → K(4s2S) + H2 were studied at the state-to-state level of theory. For the reaction process, the dynamics results indicated that the vibrational excitation of H2 was significantly more effective at promoting the reaction than the translational energy. In addition, the differential cross-sections were forward-biased scattering, which indicated that the direct abstraction mechanism plays a dominant role in the reaction. For the quenching process, the vibrational excitation of H2 molecules could improve the quenching efficiency obviously.

Observation of the geometric phase effect in the H+HD→H2+D reaction below the conical intersection

It has long been known that there is a conical intersection (CI) between the ground and first excited electronic state in the H₃ system. Its associated geometric phase (GP) effect has been theoretically predicted to exist below the CI since a long time. However, the experimental evidence has not been established yet and its dynamical origin is waiting to be elucidated. Here we report a combined crossed molecular beam and quantum reactive scattering dynamics study of the H+HD → H₂+D reaction at 2.28 eV, which is well below the CI. The GP effect is clearly identified by the observation of distinct oscillations in the differential cross section around the forward direction. Quantum dynamics theory reveals that the GP effect arises from the phase alteration of a small part of the wave function, which corresponds to an unusual roaming-like abstraction pathway, as revealed by quasi-classical trajectory calculations.

The geometric phase effect associated with a conical intersection between the ground and first excited electronic state has been predicted in the H₃ system below the conical intersection energy. The authors, by a crossed molecular beam technique and quantum dynamic calculations, provide experimental evidence and insight into its origin.

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.

Fully coupled (J > 0) time-dependent wave-packet calculations using hyperspherical coordinates for the H + O2 reaction on the CHIPR potential energy surface

ICS calculation by time dependent wavepacket approach for H + O₂ reaction using non-zero J values.