Quantum mechanical angular distributions for the F+H2 reaction

A Transition State Resonance Radically Reshapes Angular Distributions of the F + H 2 → FH(v f = 3) + H Reaction in the 62-102 meV Energy Range

Reactive angular distributions of the benchmark F + H 2(v i = 0) → FH(v f = 3) + H reaction show unusual propensity toward small scattering angles, a subject of a long debate in the literature. We use Regge trajectories to quantify the resonance contributions to state-to-state differential cross sections. Conversion to complex energy poles allows us to attribute the effect almost exclusively to a transition state resonance, long known to exist in the F + H 2 system and its isotopic variant F + HD. For our detailed analysis of angular scattering we employ the package DCS_Regge, recently developed for the purpose [Akhmatskaya E.; Sokolovski D.Comput. Phys. Commun.2022, 277, 108370].

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.

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.

New Potential Energy Surface for the H + Cl2 Reaction and Quantum Dynamics Studies

The reaction of H + Cl2 → HCl + Cl plays a crucial role in various fields. However, no previous study has investigated this reaction using accurate quantum mechanical methods. In this paper, we construct a global potential energy surface (PES) using the neural network method with more than 20,000 ab initio energies obtained by the MRCI-F12+Q method with the aug-cc-pV5Z basis and extrapolated to the complete basis set limit. The spin-orbit coupling of the Cl atom has been considered in the PES. With this new PES, product state-resolved quantum dynamics calculations for the H + Cl2 (v0 = 0, j0 = 0-2) → HCl + Cl reaction was carried out. Numerical results show that the initial rotational excitation of the Cl2 has negligible effects on the reactivity. Product state-resolved integral cross sections (ICS) and rate constants reveal that the HCl is most favorably produced in its v' = 2 vibrational state. The calculated product vibrational state-resolved and total reaction rate constants suggest that the new global PES is accurate enough, as compared with the available experimental measurements.

Vibration to Vibration: Product Energy Distribution of F + HD Crossed Molecular Beam Experiments

We investigated the F + HD(v = 1, j = 0) → HF + D reaction using the crossed molecular beam technique combined with the D atom Rydberg tagging time-of-flight spectroscopy. By detecting the products at various scattering angles for different collision energies in the range of 0.8-1.2 kcal/mol, we observed the forward-scattering products of HF(v' = 4) and determined the threshold energy for the opening of this reaction channel. Similar experiments were conducted for the F + HD(v = 0, j = 0) → HF + D reaction within the range of 1.1-1.6 kcal/mol, where forward-scattering products of HF(v' = 3) were observed, and the threshold energy for this reaction channel was determined as well. Furthermore, we measured the differential cross-sections for the F + HD → HF + D reaction in both the vibrational ground state and the excited state of HD and analyzed the vibrational quantum-state distribution of the HF products. It was found that the population of vibrational quantum states of the HF products increases synchronously with the excitation of the reactant HD vibrationally.

Coupled 3D (J ≥ 0) Time-Dependent Wave Packet Calculation for the F + H2 Reaction on Accurate Ab Initio Multi-State Diabatic Potential Energy Surfaces

We had calculated adiabatic potential energy surfaces (PESs), nonadiabatic, and spin-orbit (SO) coupling terms among the lowest three electronic states (12A', 22A', and 12A″) of the F + H2 system using the multireference configuration interaction (MRCI) level of theory, and the adiabatic-to-diabatic transformation equations were solved to formulate the diabatic Hamiltonian matrix [J. Chem. Phys. 2020, 153, 174301] for the entire region of the nuclear configuration space. The accuracy of such diabatic PESs is explored by performing scattering calculations to evaluate integral cross sections (ICSs) and rate constants. The nonadiabatic and SO effects are studied by utilizing coupled 3D time-dependent wave packet formalism with zero and nonzero total angular momentum on multiple adiabatic/diabatic surfaces calculation. We depict the convergence profiles of reaction probabilities for the reactive as well as nonreactive processes on various electronic states at different collision energies with respect to total angular momentum including all helicity quantum numbers. Finally, total ICSs are calculated as functions of collision energies for the initial rovibrational state (v = 0, j = 0) of the H2 molecule along with the temperature-dependent rate coefficient, where those quantities are compared with previous theoretical and experimental results.

A single resonance Regge pole dominates the forward-angle scattering of the state-to-state F + H2 → FH + H reaction at Etrans = 62.09 meV

The aim of the present paper is to bring clarity, through simplicity, to the important and long-standing problem: does a resonance contribute to the forward-angle scattering of the F + H2 reaction? We reduce the problem to its essentials and present a well-defined, yet rigorous and unambiguous, investigation of structure in the differential cross sections (DCSs) of the following three state-to-state reactions at a translational energy of 62.09 meV: F + H2(vi = 0, ji = 0, mi = 0) → FH(vf = 3, jf = 0, 1, 2, mf = 0) + H, where vi, ji, mi and vf, jf, mf are the initial and final vibrational, rotational and helicity quantum numbers respectively. Firstly, we carry out quantum-scattering calculations for the Fu-Xu-Zhang potential energy surface, obtaining accurate numerical scattering matrix elements for indistinguishable H2. The calculations use a time-independent method, with hyperspherical coordinates and an enhanced Numerov method. Secondly, the following theoretical techniques are employed to analyse structures in 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 second order. (b) The recently introduced "CoroGlo" test, which lets us distinguish between glory and corona scattering at forward angles for a Legendre PWS. (c) Six asymptotic (semiclassical) forward-angle glory theories and three asymptotic farside rainbow theories, valid for rainbows at sideward-scattering angles. (d) Complex angular momentum (CAM) theories of forward and backward scattering, with the Regge pole positions and residues computed by Thiele rational interpolation. Thirdly, our conclusions for the three PWS DCSs are: (a) the forward-angle peaks arise from glory scattering. (b) A broad (hidden) farside rainbow is present at sideward angles. (c) A single Regge pole contributes to the DCS across the whole angular range, being most prominent at forward angles. This proves that a resonance contributes to the DCSs for the three transitions. (d) The diffraction oscillations in the DCSs arise from NF interference, in particular, interference between the Regge pole and direct subamplitudes.

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.

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.

Beyond Born-Oppenheimer constructed diabatic potential energy surfaces for F + H2 reaction

First principles based beyond Born-Oppenheimer theory has been implemented on the F + H₂ system for constructing multistate global diabatic Potential Energy Surfaces (PESs) through the incorporation of Nonadiabatic Coupling Terms (NACTs) explicitly. The spin-orbit (SO) coupling effect on the collision process of the F + H₂ reaction has been included as a perturbation to the non-relativistic electronic Hamiltonian. Adiabatic PESs and NACTs for the lowest three electronic states (1²A', 2²A', and 1²A″) are determined in hyperspherical coordinates as functions of hyperangles for a grid of fixed values of the hyperradius. Jahn-Teller (JT) type conical intersections between the two A' states translate along C_2v and linear geometries in F + H₂. In addition, A' and A″ states undergo Renner-Teller (RT) interaction at collinear configurations of this system. Both JT and RT couplings are validated by integrating NACTs along properly chosen contours. Subsequently, we have solved adiabatic-to-diabatic transformation (ADT) equations to evaluate the ADT angles for constructing the diabatic potential matrix of F + H₂, including the SO coupling terms. The newly calculated diabatic PESs are found to be smooth, single-valued, continuous, and symmetric and can be invoked for performing accurate scattering calculations on the F + H₂ system.

Competing Dynamical Mechanisms in Inelastic Collisions of H + HF

The dynamics of inelastic collisions between HF and H has been investigated in detail by means of time-independent quantum mechanical calculations on the LWA-78 potential energy surface ( Li , G. ; et al. J. Chem. Phys. 2007 , 127 , 174302 ). Reaction probabilities, differential cross sections, and three-vector correlations have been calculated and analyzed. Our results show that there are two competing collision mechanisms that correlate with low and high impact parameters and show very different stereodynamical preferences. The mechanism promoted by high impact parameters is the only one present at low collision energies. We also observe the presence of an apparent threshold in the inelastic cross section for relatively high initial HF rotational quantum numbers, which is associated with the larger energy difference between adjacent rotational quantum states with increasing rotation.

Quantum Dynamics and Kinetics of the F + H2 and F + D2 Reactions at Low and Ultra-Low Temperatures

Integral cross sections and rate constants for the prototypical chemical reactions of the fluorine atom with molecular hydrogen and deuterium have been calculated over a wide interval of collision energy and temperature ranging from the sub-thermal (50 K) down to the ultra-cold regimes (0.5 mK). Rigorous close coupling time-independent quantum reactive scattering calculations have been carried out on two potential energy surfaces, differing only at long-range in the reactants' channel. The results show that tunnel, resonance and virtual state effects enhance under-barrier reactivity giving rise to pronounced deviations from the Arrhenius law as temperature is lowered. Within the ultra-cold domain (below 1 mK), the reactivity is governed by virtual state effects and by tunneling through the reaction barrier; in the cold regime (1 mK–1 K), the shape resonances in the entrance channel of the potential energy surface make the quantum tunneling contribution larger so enhancing cross sections and rate constants by about one order of magnitude; at higher temperatures (above 10 K), the tunneling pathway enhanced by the constructive interference between two Feshbach resonances trapped in the reaction exit channel competes with the thermally activated mechanism, as the energy gets closer to the reaction barrier height. The results show that at low temperatures cross sections and rate constants are extremely sensitive to small changes in the long-range intermolecular interaction in the entrance channel of the potential energy surface, as well as to isotopic substitution.

Dynamical resonances in chemical reactions

The transition state is a key concept in the field of chemistry and is important in the study of chemical kinetics and reaction dynamics.

Steepest-entropy-ascent nonequilibrium quantum thermodynamic framework to model chemical reaction rates at an atomistic level

The steepest entropy ascent (SEA) dynamical principle provides a general framework for modeling the dynamics of nonequilibrium (NE) phenomena at any level of description, including the atomistic one. It has recently been shown to provide a precise implementation and meaning to the maximum entropy production principle and to encompass many well-established theories of nonequilibrium thermodynamics into a single unifying geometrical framework. Its original formulation in the framework of quantum thermodynamics (QT) assumes the simplest and most natural Fisher-Rao metric to geometrize from a dynamical standpoint the manifold of density operators, which represent the thermodynamic NE states of the system. This simplest SEAQT formulation is used here to develop a general mathematical framework for modeling the NE time evolution of the quantum state of a chemically reactive mixture at an atomistic level. The method is illustrated for a simple two-reaction kinetic scheme of the overall reaction F+H_{2}⇔HF+F in an isolated tank of fixed volume. However, the general formalism is developed for a reactive system subject to multiple reaction mechanisms. To explicitly implement the SEAQT nonlinear law of evolution for the density operator, both the energy and the particle number eigenvalue problems are set up and solved analytically under the dilute gas approximation. The system-level energy and particle number eigenvalues and eigenstates are used in the SEAQT equation of motion to determine the time evolution of the density operator, thus effectively describing the overall kinetics of the reacting system as it relaxes toward stable chemical equilibrium. The predicted time evolution in the near-equilibrium limit is compared to the reaction rates given by a standard detailed kinetic model so as to extract the single time constant needed by the present SEA model.

A review of dynamical resonances in A + BC chemical reactions

The concept of the transition state has played an important role in the field of chemical kinetics and reaction dynamics. Reactive resonances in the transition-state region can dramatically enhance the reaction probability; thus investigation of the reactive resonances has attracted great attention from chemical physicists for many decades. In this review, we mainly focus on the recent progress made in probing the elusive resonance phenomenon in the simple A  +  BC reaction and understanding its nature, especially in the benchmark F/Cl  +  H₂ and their isotopic variants. The signatures of reactive resonances in the integral cross section, differential cross section (DCS), forward- and backward-scattered DCS, and anion photodetachment spectroscopy are comprehensively presented in individual prototype reactions. The dynamical origins of reactive resonances are also discussed in this review, based on information on the wave function in the transition-state region obtained by time-dependent quantum wave-packet calculations.

Complex angular momentum theory of state-to-state integral cross sections: resonance effects in the F + HD → HF(v′ = 3) + D reaction

State-to-state reactive integral cross sections (ICSs) are often affected by quantum mechanical resonances, especially in the neighborhood of a reactive threshold.

Isotope-Dependent Rotational States Distributions Enhanced by Dynamic Resonance States: A Comparison Study of the F + HD → HF(vHF = 2) + D and F + H2 → HF(vHF = 2) + H Reaction

An interesting trimodal structure in the HF (v' = 2) rotational distribution produced by the F + HD (v = 0, j = 0) reaction, but monomodal structure in the HF (v' = 2) rotational distribution produced by the F + H2 (v = 0, j = 0) reaction, were observed using a high-resolution crossed molecular beam apparatus. The rotational states of product HF (v' = 2) are much hotter in the F + HD reaction. It is uncovered that the observations are due to the dominant role of the dynamical resonance states in these two isotopic reactions. The angular potential well in the region of the resonance state of the F + HD reaction is much deeper and supports wave function with high angular kinetic energy, which in turn comes from different H tunneling processes in the F + HD and F + H2 reaction.

Resonances in rotationally inelastic scattering of OH(X2Π) with helium and neon

We present detailed calculations on resonances in rotationally and spin-orbit inelastic scattering of OH (X(2)Π, j = 3/2, F(1), f) radicals with He and Ne atoms. We calculate new ab initio potential energy surfaces for OH-He, and the cross sections derived from these surfaces compare well with the recent crossed beam scattering experiment of Kirste et al. [Phys. Rev. A 82, 042717 (2010)]. We identify both shape and Feshbach resonances in the integral and differential state-to-state scattering cross sections, and we discuss the prospects for experimentally observing scattering resonances using Stark decelerated beams of OH radicals.

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.

CHEMICAL REACTIONS IN THE O(1D) + HCl SYSTEM III

Using the accurate global potential energy surfaces for the 11A′′ and 21A′ states reported in the previous sister Paper I, detailed quantum dynamics calculations are performed for these adiabatic surfaces separately for J = 0 (J: total angular momentum quantum number). In addition to the significant overall contributions of these states to the title reactions reported in the second Paper II of this series, quantum dynamics on these excited potential energy surfaces (PES) are clarified in terms of the PES topographies, which are quite different from that of the ground PES. The reaction mechanisms are found to be strongly selective and nicely explained as vibrationally nonadiabatic transitions in the vicinity of potential ridge.

A STATE-TO-STATE QUANTUM DYNAMICAL STUDY OF THE H + HBr REACTION

The time-dependent wave packet method was used to calculate the state-to-state differential cross sections for abstraction and exchange processes in the title reaction on the Kurosaki–Takayanagi potential energy surface [Kurosaki Y, Takayanagi T, J Chem Phys119:7838, 2003], with the reactant HBr initially in the ground rovibrational state. It is found that the trend in the product distributions is similar for abstraction and exchange processes, but the differential cross sections are very different. For the exchange reaction, the product is mainly scattered in the backward hemisphere for collision energy up to 2.0 eV, although forward scattering gradually shows up in high collision energies. While for abstraction reaction, the differential cross section changes substantially with the collision energy, moving from predominantly backward peaked at low collision energy to predominantly forward peaked at high collision energy. The rovibrational state resolved differential cross section at collision energy of 2.0 eV exhibits two peaks for the abstraction reaction, one is around the angle of 50°, and the other at 0°. It is found that the peaks around 50°, are below the corresponding maximum j' lines provided by the kinematic constraint model, while the forward-scattered peaks straddle both sides of the kinematic limit, and are likely contributed from both the direct and the migratory reaction mechanisms as proposed by Zare and coworkers.