Time-dependent wavepacket investigation of state-to-state reactive scattering of Cl with para-H2 including the open-shell character of the Cl atom

Diabatic Potential Energy Surfaces of SrH2+ and Dynamics Studies of the Sr+(5s2S) + H2 Reaction

Based on the ab initio energy points of both ground and excited states, a neural network fitting method combined with a specific function was successfully used to construct the diabatic potential energy surfaces (PESs) of the SrH2+ system. The topographical features of the diabatic PESs were examined in detail. The results indicate that the nonadiabatic transition characteristics between ground and excited states are accurately described by the newly constructed diabatic PESs. To verify the validity and applicability of the diabatic PESs, as well as the nonadiabatic effects during the reaction process, the quantum dynamics studies of the Sr+(5s2S) + H2 reaction were performed based on both adiabatic and diabatic PESs. The dynamics results indicate that adiabatic dynamics results are dozens of times larger than those of nonadiabatic. This illustrates the significant effect of nonadiabaticity, indicating that adiabatic dynamics results often overestimate the actual values. The integral cross sections (ICSs) were calculated and compared with the experimental data. The diabatic ICSs are in good agreement with the experimental results. The reasonable dynamics results indicate that the newly constructed diabatic PESs are suitable for the relevant dynamics studies.

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.

Development of the Time-Independent Methods for the Cl + H2/F + HD Reaction Using Hyper-Spherical Coordinates Including (Full) Spin-Orbit Characteristics

Recently, a combined study of high-resolution molecular crossed beam experiment and accurate full-dimensional time-dependent theory, including full spin-orbit characteristics on the effect of electronic spin and orbital angular momenta in the F + HD reaction, was reported by some of us, focusing on the partial wave resonance phenomenon (Science 2021, 371, 936-940). It revealed that the time-dependent theory could explain all of the details observed in the high-resolution experiment. Here, we develop two time-independent close-coupling methods using hyperspherical coordinates, including the two-state model, where only a part of the spin-orbit characteristics is considered, and the six-state model, where the full spin-orbit characteristics is considered. With these two newly developed theoretical models and the adiabatic theoretical model, the detailed reaction dynamics of the F + HD (v = 0, j = 0) reaction and the Cl + H2 (v = 0, j = 0) reaction are investigated and compared. Some of the results are compared with the time-dependent quantum wave packet theory and the experimental observations, and good agreements have been obtained, which suggests the validity of the pure-procession approximation in the six-state model using different theoretical methods. This work demonstrates the ability of the reactive scattering theory including full spin-orbit characteristics for describing the reactions of a halogen atom plus hydrogen molecule and its isotopologues.

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.

Significant Isotope Effects from the Nonadiabatic Couplings in the Cl(2P) + HD(v = 0, j = 0) Reaction

The impact of non-Born-Oppenheimer couplings on the isotopic effects in the reaction of the Cl(2P) atom with the HD (v = 0, j = 0) molecule is investigated with our recently developed nonadiabatic time-independent quantum scattering methods, where the full open-shell characteristics are included in the six-state model, and also with the recently developed two-state model solving by time-independent methods, where part of the open-shell characteristic is included. The same reaction is also calculated with the simple adiabatic model using the lowest adiabatic potential energy surface. Compared with the results from different models, it is found that the reactivity of the Cl + HD → HCl + D channel is significantly overestimated in the adiabatic model. In contrast, the reactivity of the other channel agrees well with the nonadiabatic models. This is due to the van der Waals well in the reactant channel being changed a lot by including the nonadiabatic couplings. These quantum dynamics calculations suggest that sometimes the adiabatic model should be used with caution; otherwise, it may result in significant deviations for some reactions.

Quantum Dynamics of Nonadiabatic Renner-Teller Effects in Atom + Diatom Collisions

We review the quantum nonadiabatic dynamics of atom + diatom collisions due to the Renner-Teller (RT) effect, i.e., to the Hamiltonian operators that contain the total spinless electronic angular momentum L̂. As is well-known, this rovibronic effect is large near collinear geometries when at least one of the interacting states is doubly degenerate. In general, this occurs in insertion reactions and at short-range, where the potential wells exhibit deep minima and support metastable complexes. Initial-state-resolved reaction probabilities, integral cross sections, and thermal rate constants are calculated via the real wavepacket method, solving the equation of motion with an approximated or with an exact spinless RT Hamiltonian. We present the dynamics of 10 single-channel or multichannel reactions showing how RT effects depend on the product channels and comparing with the Born-Oppenheimer (BO) approximation or coexisting conical-intersection (CI) interactions. RT effects not only can significantly modify the adiabatic dynamics or correct purely CI results, but also they can be very important in opening collision channels which are closed at the BO or CI level, as in electronic-quenching reactions. In the OH(A²Σ⁺) + Kr electronic quenching, where both nonadiabatic effects (CI and RT) coexist, they are in competition because CI dominates the reactivity but RT couplings reduce the large CI cross section and open a CI-forbidden evolution toward products, so that CI + RT quantum results are in good agreement with experimental or semiclassical findings. The different roles of these couplings are due to the unlike nuclear geometries where they are large: rather far from or near to linearity for CI or RT, respectively. The OH(A²Σ⁺) + Kr electronic quenching was investigated with the exact RT Hamiltonian, validating the approximated one, which was employed for all other collisions.

Quantum interference between spin-orbit split partial waves in the F + HD → HF + D reaction

The effect of electron spin-orbit interactions on chemical reaction dynamics has been a topic of much research interest. Here we report a combined experimental and theoretical study on the effect of electron spin and orbital angular momentum in the F + HD → HF + D reaction. Using a high-resolution imaging technique, we observed a peculiar horseshoe-shaped pattern in the product rotational-state-resolved differential cross sections around the forward-scattering direction. The unusual dynamics pattern could only be explained properly by highly accurate quantum dynamics theory when full spin-orbit characteristics were considered. Theoretical analysis revealed that the horseshoe pattern was largely the result of quantum interference between spin-orbit split-partial-wave resonances with positive and negative parities, providing a distinctive example of how spin-orbit interaction can effectively influence reaction dynamics.

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.

Quantum interference in H + HD → H2 + D between direct abstraction and roaming insertion pathways

Understanding quantum interferences is essential to the study of chemical reaction dynamics. Here, we provide an interesting case of quantum interference between two topologically distinct pathways in the H + HD → H2 + D reaction in the collision energy range between 1.94 and 2.21 eV, manifested as oscillations in the energy dependence of the differential cross section for the H2 (v′ = 2, j′ = 3) product (where v′ is the vibrational quantum number and j′ is the rotational quantum number) in the backward scattering direction. The notable oscillation patterns observed are attributed to the strong quantum interference between the direct abstraction pathway and an unusual roaming insertion pathway. More interestingly, the observed interference pattern also provides a sensitive probe of the geometric phase effect at an energy far below the conical intersection in this reaction, which resembles the Aharonov–Bohm effect in physics, clearly demonstrating the quantum nature of chemical reactivity.

Imaging the State-to-State Dynamics of the H + D2 → HD + D Reaction at 1.42 eV

High-resolution state-resolved differential cross sections (DCSs) are of great importance in understanding quantum reaction dynamics, and they are the most detailed observables that can be experimentally measured. Here we report a synergic crossed molecular beam and quantum reaction dynamics study on the H + D₂ reaction. With the time-sliced velocity map ion imaging (VMI) technique and the near-threshold ionization scheme, we acquired the product rovibrational state-resolved DCSs of the H + D₂ (v = 0, j = 0) → HD (v', j') + D reaction at a collision energy of 1.42 eV. For HD products with small j' quantum numbers, significant forward scattering with clear angular oscillations was observed. The forward scattering disappears for the rotational states with large j' quantum numbers. Interestingly, as the j' number increases, the peak of the DCS shifts from backward to sideways systematically. The experimental observation agrees very well with theoretical quantum mechanical dynamics results, which reveals that the systematic shift of the peak in the DCS from backward scattering to sideways scattering can be understood very well with the strong correlation between the product rotational quantum number j' and the specific partial waves (J = 3-12), whereas the forward angular oscillations are from the coherent summation of larger partial waves.

Conservation of Angular Momentum in Direct Nonadiabatic Dynamics

Direct nonadiabatic dynamics is used to study processes involving multiple electronic states from small molecules to materials. Compared with dynamics with fitted analytical potential energy surfaces, direct dynamics is more user-friendly in that it obtains all needed energies, gradients, and nonadiabatic couplings (NACs) by electronic structure calculations. However, the NAC that is usually used does not conserve angular momentum or the center of mass in widely used mixed quantum-classical nonadiabatic dynamics algorithms, in particular, trajectory surface hopping, semiclassical Ehrenfest, and coherent switching with decay of mixing. We show that by using a projection operator to remove the translational and rotational components of the originally computed NAC, one can restore the conservation.

Nonadiabatic Ultracold Quantum Reactive Scattering of Hydrogen with Vibrationally Excited HD(v = 5-9)

The results from electronically non-adiabatic and adiabatic quantum reactive scattering calculations are presented for the H + HD(v = 5-9) → H + HD(v', j') reaction at ultracold collision energies from 10 nK to 60 K. Several experimentally verifiable signatures of the geometric phase are reported in the total and vibrationally and rotationally resolved rate coefficients. Most notable is the predicted 2 orders of magnitude enhancement of the rotationally resolved ultracold rates of odd symmetry relative to those of even symmetry. Prominent shape resonances appear at higher collision energies (100 mK to 20 K), which could be measured experimentally. Significant geometric phase effects are also reported on the resonance energies and lifetimes. In particular, an enhancement (suppression) of the l = 1 (l = 2) shape resonances for HD(v = 5, 6) is predicted for even symmetry relative to those of odd symmetry.

First-principles dynamics of collisional intersystem crossing: resonance enhanced quenching of C(1D) by N2

Intersystem crossing is a common and important nonadiabatic process in molecular systems. Revealed with a first-principles investigation, the quenching of C(¹D) by N₂ is efficient due to multiple passages via long-lived collisional resonances, despite relatively small spin–orbit couplings.

Influence of rovibrational excitation on the non-diabatic state-to-state dynamics for the Li(2p) + H2 → LiH + H reaction

The non-adiabatic state-to-state dynamics of the Li(2p) + H2 → LiH + H reaction has been studied using the time-dependent wave packet method, based on a set of diabatic potential energy surfaces recently developed by our group. Integral cross sections (ICSs) can be increase more than an order of magnitude by the vibrational excitation of H2, whereas the ICSs are barely affected by the rotational excitation of H2. Moreover, ICSs of the title reaction with vibrationally excited H2 decrease rapidly with increasing collision energy, which is a typical feature of non-threshold reaction. This phenomenon implies that the title reaction can transformed from an endothermic to an exothermic reaction by vibrational excitation of H2. With the increase of the collision energy, the sideways and backward scattered tendencies of LiH for the Li(2p) + H2(v = 0, j = 0, 1) → LiH + H reactions are enhanced slightly, while the backward scattering tendency of LiH for the Li(2p) + H2(v = 1, j = 0) → LiH + H reaction becomes remarkably weakened. For the reaction with vibrationally excited H2 molecule, both direct and indirect reaction mechanism exist simultaneously.

A new potential energy surface of LiHCl system and dynamic studies for the Li(2S) + HCl(X1Σ+) → LiCl(X1Σ+) + H(2S) reaction

A new global potential energy surface (PES) is constructed for the ground state of LiHCl system based on high-quality ab initio energy points calculated using multi-reference configuration interaction calculations with the Davidson correction. The AVQZ and WCVQZ basis sets are employed for H and Li atoms, respectively. To compensate the relativistic effects of heavy element, the AWCVQZ-DK basis set is employed for Cl atom. The neural network method is used for fitting the PES, and the root mean square error is small (1.36 × 10^-2 eV). The spectroscopic constants of the diatoms obtained from the new PES agree well with experimental data. The geometric characteristics of the transition state and the complex are examined and compared with the previous theoretical values. To study the reaction dynamics of the Li(²S) + HCl(X¹Σ⁺) → LiCl(X¹Σ⁺) + H(²S) reaction, quantum reactive scattering dynamics calculations using collection reactant-coordinate-based wave packet method are conducted based on the new PES. The results of the reaction probabilities indicate that a small barrier exists along the reaction path as observed from the PES. The integral cross section curves reveal that the product molecule LiCl is easily excited. In addition, the reaction is dominated by forward scattering, and similar pattern is observed from Becker's experiment.

Global diabatic potential energy surfaces and quantum dynamical studies for the Li(2p) + H2(X(1)Σ(+)g) → LiH(X(1)Σ(+)) + H reaction

The global diabatic potential energy surfaces which are correlated with the ground state 1A' and the excited state 2A' of the Li(2p) + H2 reaction are presented in this study. The multi-reference configuration interaction method and large basis sets (aug-cc-pVQZ for H atom and cc-pwCVQZ for Li atom) were employed in the ab initio single-point energy calculations. The diabatic potential energies were generated by the diabatization scheme based on transition dipole moment operators. The neural network method was utilized to fit the matrix elements of the diabatic energy surfaces, and the root mean square errors were extremely small (3.69 meV for , 5.34 meV for and 5.06 meV for ). The topographical features of the diabatic potential energy surfaces were characterized and the surfaces were found to be sufficiently smooth for the dynamical calculation. The crossing seam of the conical intersections between the and surfaces were pinpointed. Based on this new analytical diabatic potential energy surfaces, time-dependent wave packet calculation were conducted to investigate the mechanism of the title reaction. At low collision energies, the product LiH molecule tends to forward scattering, while at high collision energies, the forward and backward scatterings exist simultaneously.

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 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 reactive scattering of O(3P)+H2 at collision energies up to 4.4 eV

We report the results of quantum scattering calculations for the O((3)P)+H2 reaction for a range of collision energies from 0.4 to 4.4 eV, important for astrophysical and atmospheric processes. The total and state-to-state reactive cross sections are calculated using a fully quantum time-independent coupled-channel approach on recent potential energy surfaces of (3)A' and (3)A″ symmetry. A larger basis set than in the previous studies was used to ensure single-surface convergence at higher energies. Our results agree well with the published data at lower energies and indicate the breakdown of reduced dimensionality approach at collision energies higher than 1.5 eV. Differential cross sections and momentum transfer cross sections are also reported.

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.

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.

INVESTIGATION OF THE CONTRIBUTION FOR THE NONCOLLINEAR CHANNEL OF THE F(2P3/2,2P1/2) + H2/D2 REACTIONS ON FOUR DIABATIC POTENTIAL ENERGY SURFACES

We performed the nonadiabatic time-dependent wave packet calculation on the four diabatic potential energy surfaces, which have the different barrier height, to investigate the contribution of the noncollinear channel for the F (2P) + H2/D2 (v = j = 0) reactions. The reaction probabilities, integral cross-sections, and rate constants are presented. The results indicate that the probabilities as the function of the collision energy have an obvious translation. The reactive activity of the reactions comes from the noncollinear reactive channel. The bent barrier height would decrease the reactive activity. The integral cross-sections are in the order of AWS < LWA-5 < LWA-78 ≈ MASW, which is opposite to that of the bent barrier height. At the lower temperature, the difference of the rate constants is unambiguous. As the temperature increases, the difference reduces. At the higher temperature, the rate constants computed on the four potential energy surfaces are close.