Quantum scattering study of electronic Coriolis and nonadiabatic coupling effects in O(1D)+H2→OH+H

An effective approximation of Coriolis coupling in reactive scattering: application to the time-dependent wave packet calculations

Coriolis coupling plays a crucial role in reactive scattering, but dynamics calculations including the complete Coriolis coupling significantly increase the difficulty of numerical evolution due to the corresponding expensive matrix processing. The coupled state approximation that completely ignores the off-diagonal Coriolis coupling saves computational cost significantly but its error is usually unacceptable. In this paper, an improved coupled state approximation inspired by recently published results [D. Yang, X. Hu, D. H. Zhang and D. Xie, J. Chem. Phys., 2018, 148, 084101.] of the inelastic scattering problem is extended to deal with the reactive scattering. The calculations using the time-dependent wave packet method reveal that the new method can accurately reproduce the rigorous results of the H + HD (j0 < 3) → D + H2 reaction and immensely improve the computational efficiency. Additionally, we extend the new method to the non-adiabatic Li(2p) + H2 (v0 = 0, j0 = 0, 1) → H + LiH reaction, showcasing its advantages of low computational cost and high accuracy.

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.

A combined theoretical and experimental investigation of the kinetics and dynamics of the O(1D) + D2 reaction at low temperature

The O(¹D) + H₂ reaction is a prototype for simple atom-diatom insertion type mechanisms considered to involve deep potential wells. While exact quantum mechanical methods can be applied to describe the dynamics, such calculations are challenging given the numerous bound quantum states involved. Consequently, efforts have been made to develop alternative theoretical strategies to portray accurately the reactive process. Here we report an experimental and theoretical investigation of the O(¹D) + D₂ reaction over the 50-296 K range. The calculations employ three conceptually different approaches - mean potential phase space theory, the statistical quantum mechanical method and ring polymer molecular dynamics. The calculated rate constants are in excellent agreement over the entire temperature range, exhibiting only weak temperature dependence. The agreement between experiment and theory is also very good, with discrepancies smaller than 26%. Taken together, the present and previous theoretical results validate the hypothesis that long-lived complex formation dominates the reaction dynamics at low temperature.

The effect of intersystem crossings in N((2)D) + H2 collisions

The transitions between quartet and doublet states of the NH2 molecule are studied for the first time, allowing the evaluation of the N((4)S) + H2 reactive channel. High level ab initio calculations of the spin-orbit coupling are performed over the whole configurational space of the NH2 molecule and fitted to a proposed analytic form. Quasiclassical trajectories coupled with the surface hopping method are employed to calculate reaction cross section and rate constants. The reaction is largely affected by the initial rovibrational states of H2, while the formation of long-lived complexes enhances the reaction probability.

A finite-element visualization of quantum reactive scattering. II. Nonadiabaticity on coupled potential energy surfaces

This is the second in a series of papers detailing a MATLAB based implementation of the finite element method applied to collinear triatomic reactions. Here, we extend our previous work to reactions on coupled potential energy surfaces. The divergence of the probability current density field associated with the two electronically adiabatic states allows us to visualize in a novel way where and how nonadiabaticity occurs. A two-dimensional investigation gives additional insight into nonadiabaticity beyond standard one-dimensional models. We study the F((2)P) + HCl and F((2)P) + H2 reactions as model applications. Our publicly available code (http://www2.chem.umd.edu/groups/alexander/FEM) is general and easy to use.

Quantum wave packet and quasiclassical trajectory studies of the reaction H(2S) + CH(X2 Π; v = 0, j = 1) → C(1D) + H2 (X1 Σg+): Coriolis coupling effects and stereodynamics

The quantum mechanics (QM) and quasiclassical trajectory (QCT) calculations have been carried out for the title reaction with the ground minimal allowed rotational state of CH (j = 1) on the 1 (1)A' potential energy surface. For the reaction probability at total angular momentum J = 0, a similar trend of the QM and QCT calculations is observed, and the QM results are larger than the latter almost in the whole considered energy range (0.1-1.5 eV). The QCT integral cross sections are larger than the QM results with centrifugal sudden approximation, while smaller than those from QM method including Coriolis coupling for collision energies bigger than 0.25 eV. The quantum wave-packet computations show that the Coriolis coupling effects get more and more pronounced with increasing of J. In addition to the scalar properties, the stereodynamical properties, such as the average rotational alignment factor <P2 (j'•k)>, the angular distributions P(θr ), P(ϕr ), P(θr ,ϕr ), and the polarization-dependent generalized differential cross sections have been explored in detail by QCT approach.

Coriolis coupling and nonadiabaticity in chemical reaction dynamics

The nonadiabatic quantum dynamics and Coriolis coupling effect in chemical reaction have been reviewed, with emphasis on recent progress in using the time-dependent wave packet approach to study the Coriolis coupling and nonadiabatic effects, which was done by K. L. Han and his group. Several typical chemical reactions, for example, H+D(2), F+H(2)/D(2)/HD, D(+)+H(2), O+H(2), and He+H(2)(+), have been discussed. One can find that there is a significant role of Coriolis coupling in reaction dynamics for the ion-molecule collisions of D(+)+H(2), Ne+H(2)(+), and He+H(2)(+) in both adiabatic and nonadiabatic context.

Nonadiabatic state-to-state reactive collisions among open shell reactants with conical intersections: the OH((2)Pi) + F((2)P) example

Accurate wave packet calculations on the OH((2)Pi) + F((2)P) → O((3)P) + HF((1)Sigma(+)) reactive collisions are performed using a recently proposed coupled diabatic states. Adiabatic and nonadiabatic dynamics are compared in detail, analyzing the final state distribution of products. It is found that with the new surfaces a significant increase of the rate constant is obtained, with noticeable nonadiabatic effects. The inclusion of the spin-orbit splittings for the calculation of the electronic partition function produces an important increase of the reaction rate constants, yielding a rather good agreement with the experimental results. It is also concluded that spin-orbit couplings are also necessary in the entrance channel to describe this reaction.

Effect of Coriolis coupling in chemical reaction dynamics

It is essential to evaluate the role of Coriolis coupling effect in molecular reaction dynamics. Here we consider Coriolis coupling effect in quantum reactive scattering calculations in the context of both adiabaticity and nonadiabaticity, with particular emphasis on examining the role of Coriolis coupling effect in reaction dynamics of triatomic molecular systems. We present the results of our own calculations by the time-dependent quantum wave packet approach for H + D2 and F(2P3/2,2P1/2) + H2 as well as for the ion-molecule collisions of He + H2 +, D(-) + H2, H(-) + D2, and D+ + H2, after reviewing in detail other related research efforts on this issue.

Dynamics of insertion reactions of H2 molecules with excited atoms

Recent progress in the study of insertion reactions of hydrogen molecules with excited atoms is reviewed in this article. In particular, the dynamics of the reaction of O(1D), N(2D), C(1D), and S(1D) with H2 and its isotopomers, which have received a great deal of attention over the past decade, are examined in detail. All of these systems have in common the existence of several potential energy surfaces (PES) correlating with the reagents' states, and consequently, they can give rise to reaction following different adiabatic and nonadiabatic pathways. The main contribution, however, arises from their ground singlet PESs which feature the existence of deep wells with small or null barriers for insertion. Accordingly, these reactions proceed mainly via formation of relatively long-lived collision complexes and display an overall nearly statistical behavior. In spite of their similarities, the various reactions have peculiar characteristics caused by important differences of their respective PESs. The contribution of excited PES to the global reactivity, which has also become an important issue and a challenge both for theory and experiment, is also examined. The different theoretical approaches are discussed in the text, along with the experimental results obtained by a variety of techniques. The recent exact quantum treatments of these reactive systems together with the development of a rigorous statistical model have contributed to a very accurate description which in many cases matches very well the detailed measurements. The quasi-classical trajectory (QCT) method has also provided a fairly accurate description of the reaction dynamics for these systems. In particular, the analysis in terms of collision times has yielded interesting clues about the reaction mechanisms.

Experimental and quantum dynamical study on an asymmetric insertion reaction: state-to-state dynamics of O(1D) + HD(1Sigmag+, v' = 0, j' = 0)-->OH(Pi, v", N") + D(2S)

Quantum-state-resolved differential cross sections of the O((1)D) + HD --> OH + D reaction at the collision energy of 7.11 kJ/mol has been determined experimentally and theoretically. The results of the time-dependent wave-packet calculations are overall in good agreement with the crossed beam scattering data, providing a benchmark example of an asymmetric insertion reaction at the state-to-state scattering level. The good agreement between experiment and theory suggests that the underlying ground potential energy surface is generally correct and that the nonadiabatic effect involving the electronic excited pathway is apparently small in this system.

Resonances of CH2(ãA11) and their roles in unimolecular and bimolecular reactions

Low-lying resonances of the CH2(a 1A1) system (J=0) in an accurate ab initio potential energy surface are studied using a filter-diagonalization method. The width of these resonances fluctuates by more than two orders of magnitude and on average increases with the energy. Analysis of the resonance states concludes that the unimolecular decay of the excited molecular system near the dissociation threshold is neither mode specific nor statistical state specific. This is apparently due to remnant regularity embedded in the largely chaotic classical phase space, as evidenced by periodic orbit analysis. As a result, the Rice-Ramsperger-Kassel-Marcus and statistical adiabatic channel models overestimate the average unimolecular decay rate. The implications of the resonances for the bimolecular C(1D)+H2 reaction are also discussed.

Product multiplet branching in the O(1D)+H2→OH(2Π)+H reaction

The statistical model of atom-diatom insertion reactions is combined with coupled-states capture theory and used to calculate product multiplet-resolved integral cross sections for the title reaction. This involves an ab initio determination of the four electronic potential energy surfaces that correlate with the products ((1,3)A(') and (1,3)A(")), and an accurate description of the electronic and spin-orbit couplings between them. The dependence of the resulting cross sections on the final-state rotational quantum number shows a statistical behavior similar to that observed in earlier studies of the reaction in which only the lowest ((1)A(')) potential was retained. In addition, however, the present calculations provide information on the branching between the OH((2)Pi) multiplet levels. Although the two spin-orbit manifolds are predicted to be equally populated, we find a strong propensity for the formation of the Pi(A(')) Lambda-doublet states. These two predictions confirm the experimental results of Butler, Wiesenfeld, Gericke, Brouard, and their co-workers. The nonstatistical population of the OH Lambda-doublet levels is a consequence of the bond breaking in the intermediate H(2)O complex and is preserved through the multiple curve crossings as the products separate. This exit-channel coupling is correctly described by the present theory.

A field theoretical approach to calculate electronic Born-Oppenheimer coupling terms

In this paper we suggest to consider the spatial distribution of the Born-Oppenheimer nonadiabatic coupling terms as fields which are created by sources, located at degeneracy points, and which can be derived using the ordinary mathematical tools of field theory. It is shown that the curl-divergence equations as formed within a given Hilbert space [M. Baer, Chem. Phys. Lett. 35, 112 (1975)] can be converted into a set of inhomogeneous coupled Poisson equations which are solved for a given set of boundary conditions. The method is applied to the three-state Hilbert subspace of the H(3) system. The numerical results are compared with ab initio calculations for which a very encouraging fit is found.

The effect of spin–orbit coupling on fast neutral chemical reaction O(3P)+CH3→CH3O

The effect of nonadiabatic transitions through the spin-orbit couplings has been investigated on the fast neutral reaction, O((3)P)+CH(3)-->CH(3)O. Adiabatic potential energies and the spin-orbit coupling terms have been evaluated for the four electronic states of CH(3)O ((2)E, (2)A(2), (4)E, and (4)A(2)) that correlate with the O((3)P)+CH(3) asymptote, as a function of CO distance and OCH angle under the C(3v) symmetry, by ab initio electronic structure calculations using multireference internally contracted single and double excitation configuration interaction method with the 6-311G(2df,2pd) basis sets. Multistate quantum reactive scattering calculations have been carried out with the use of thus obtained potential energies and spin-orbit coupling matrices, based on the generalized R-matrix propagation method. The calculated thermal rate constants show a slight positive dependence on temperature in a range between 50 and 2000 K, supporting the previous experimental results. It is shown that the spin-orbit coupled excited states give rise to reflections over the centrifugal barrier due to the quantum interference. Classical capture calculations yield larger rate constants due to the neglect of quantum reflections. It is concluded that the effect of nonadiabatic transitions is of minor importance on the overall reactivity in this reaction.

Quantum scattering calculations on chemical reactions

This review discusses recent quantum scattering calculations on bimolecular chemical reactions in the gas phase. This theory provides detailed and accurate predictions on the dynamics and kinetics of reactions containing three atoms. In addition, the method can now be applied to reactions involving polyatomic molecules. Results obtained with both time-independent and time-dependent quantum dynamical methods are described. The review emphasises the recent development in time-dependent wave packet theories and the applications of reduced dimensionality approaches for treating polyatomic reactions. Calculations on over 40 different reactions are described.