The F+HD→DF(HF)+H(D) reaction revisited: Quasiclassical trajectory study on anabinitiopotential energy surface and comparison with molecular beam experiments

Classification of the HCN isomerization reaction dynamics in Ar buffer gas via machine learning

The effect of the presence of Ar on the isomerization reaction HCN ⇄ CNH is investigated via machine learning. After the potential energy surface function is developed based on the CCSD(T)/aug-cc-pVQZ level ab initio calculations, classical trajectory simulations are performed. Subsequently, with the aim of extracting insights into the reaction dynamics, the obtained reactivity, that is, whether the reaction occurs or not under a given initial condition, is learned as a function of the initial positions and momenta of all the atoms in the system. The prediction accuracy of the trained model is greater than 95%, indicating that machine learning captures the features of the phase space that affect reactivity. Machine learning models are shown to successfully reproduce reactivity boundaries without any prior knowledge of classical reaction dynamics theory. Subsequent analyses reveal that the Ar atom affects the reaction by displacing the effective saddle point. When the Ar atom is positioned close to the N atom (resp. the C atom), the saddle point shifts to the CNH (HCN) region, which disfavors the forward (backward) reaction. The results imply that analyses aided by machine learning are promising tools for enhancing the understanding of reaction dynamics.

Li + HF and Li + HCl Reactions Revisited I: QCT Calculations and Simulation of Experimental Results

The Li + HF and Li + HCl reactions share some common features. They have the same kinematics, relatively small barrier heights, bent transition states, and are both exothermic when the zero point energy is considered. Nevertheless, the pioneering crossed beam experiments by Lee and co-workers in the 80s (Becker et al., J. Chem. Phys. 1980, 73, 2833) revealed that the dynamics of the two reactions differ significantly, especially at low collision energies. In this work, we present theoretical simulations of their results in the laboratory frame (LAB), based on quasiclassical trajectories and obtained using accurate potential energy surfaces. The calculated LAB angular distributions and time-of-flight spectra agree well with the raw experimental data, although our simulations do not reproduce the experimentally derived center-of-mass (CM) differential cross section and velocity distributions. The latter were derived by forward convolution fitting under the questionable assumption that the CM recoil velocity and scattering angle distribution were uncoupled, while our results show that the coupling between them is relevant. Some important insights into the reaction mechanism discussed in the article by Becker et al. had not been contrasted with those that can be extracted from the theoretical results. Among them, the correlation between the angular momenta involved in the reactions has also been examined. Given the kinematics of both systems, the reagent orbital angular momentum, l , is almost completely transformed into the rotation of the product diatom, j'. However, contrary to the coplanar mechanism proposed in the original paper, we find that the initial and final relative orbital angular momenta are not necessarily parallel. Both reactions are found to be essentially direct, although about 15% of the LiFH complexes live longer than 200 fs.

Quantum Dynamics Studies of the Significant Intramolecular Isotope Effects on the Nonadiabatic Be+(2P) + HD → BeH+/BeD+ + D/H Reaction

Quantum time-dependent wave packet dynamics studies on the nonadiabatic Be⁺(²P) + HD → BeH⁺/BeD⁺ + D/H reaction are performed for the first time employing recently constructed diabatic potential energy surfaces. Strong intramolecular isotope effects and unusual results are presented, which are attributed to the dynamic effects of shallow wells induced by avoided crossing on the diagonal V₂₂^d surface. The BeH⁺ + D and BeD⁺ + H channels are dominated by high-J and low-J partial waves, respectively. The BeD⁺/BeH⁺ branching ratio is larger than 10 at low energy and gradually decreases with increasing collision energy. The BeH⁺ product is primarily distributed at low vibrational states, whereas there exists an obvious population inversion of vibrational states on the BeD⁺ product. The results of differential cross sections suggest that the formation of the BeH⁺ + D channel favors a direct reaction process, while the BeD⁺ + H channel is mainly generated by the complex-forming mechanism.

Orbiting resonances in the F + HD (v = 0, 1) reaction at very low collision energies. A quantum dynamical study

Time-independent, fully converged, quantum dynamical calculations have been performed for the F + HD (v = 0, j = 0) and F + HD (v = 1, j = 0) reactions on an accurate potential energy surface down to collision energies of 0.01 meV. The two isotopic exit channels, HF + D and DF + H, have been investigated. The calculations reproduce satisfactorily the Feshbach resonance structures for collision energies between 10 and 40 meV, previously reported in the literature for the HF + D channel. Contrary to the results of a former literature work, vibrational excitation of HD is found to enhance reactivity in all cases down to the lowest collision energy investigated. Shape-type orbiting resonances are found for collision energies lower than 2 meV. The resonances appear as peaks in the reaction cross sections that are associated to specific values of the total angular momentum, J. In contrast with the Feshbach resonances at higher energies, the orbiting resonance structure, which is caused by the van der Waals well of the entrance channel, is identical for the HF + D and DF + H exit channels. The orbiting resonance peaks for F + HD (v = 0) are very small, but those for F + HD (v = 1) could be observed, in principle, with a combination of Raman pumping and merged beams methods.

Competition between the H- and D-atom transfer channels in the H2O+ + HD reaction: reduced-dimensional quantum and quasi-classical studies

The ion-molecule reaction between a water cation and a hydrogen molecule has recently attracted considerable interest due to its importance in astrochemistry. In this work, the intramolecular isotope effect of the H₂O⁺ + HD reaction is investigated using a seven-dimensional initial state-selected time-dependent wave packet approach as well as a full-dimensional quasi-classical trajectory method on a full-dimensional ab initio global potential energy surface. The calculated branching ratios for the formation of H₃O⁺ and H₂DO⁺via H- and D-transfer agree reasonably well with the experimental values. The preference to the formation of the H₃O⁺ product observed using the experiment at low collision energies is reproduced by theoretical calculations and explained by a one-dimensional effective potential model.

A comparative study of the Au + H₂, Au⁺ + H₂, and Au⁻ + H₂ systems: Potential energy surfaces and dynamics of reactive collisions

In order to study the Au(-) + H2 collision, a new global potential energy surface (PES) describing the ground electronic state of AuH2(-) system is developed and compared with the PESs of the neutral [Zanchet et al., J. Chem. Phys. 132, 034301 (2010)] and cationic systems [Anaís et al., J. Chem. Phys. 135, 091102 (2011)]. We found that Au(-) - H2 presents a H-Au-H insertion minimum attributed to the stabilization of the LUMO 3b2 orbital, which can be considered as the preamble of the chemisorption well appearing in larger gold clusters. While the LUMO orbital is stabilized, the HOMO 6a1 is destabilized, creating a barrier at the geometry where the energy orbitals' curves are crossing. In the anion, this HOMO is doubly occupied, while in the neutral system is half-filled and completely empty in the cation, explaining the gradual disappearance of the well and the barrier as the number of electrons decreases. The cation presents a well in the entrance channel partially explained by electrostatic interactions. The three systems' reactions are highly endothermic, by 1.66, 2.79, and 3.23 eV for AuH, AuH(+), and AuH(-) products, respectively. The reaction dynamics is studied using quasi-classical trajectory method for the three systems. The one corresponding to the anionic system is new in this work. Collision energies between 1.00 and 8.00 eV, measured for the cation, are in good agreement with the simulated cross section for the AuH(+). It was also found that the total fragmentation, in three atoms, competes becoming dominant at sufficiently high energy. Here, we study the competition between the two different reaction pathways for the anionic, cationic, and neutral species, explaining the differences using a simple model based on the topology of the potential energy surfaces.

The Cl + O3 reaction: a detailed QCT simulation of molecular beam experiments

QCT calculations have been carried out to determine angle–velocity differential cross-sections to simulate the results of molecular beam experiments.

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.

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.

USING THE QUASI-CLASSICAL TRAJECTORY METHOD TO STUDY THE F + H2 AND ITS ISOTOPIC VARIANTS: VECTOR CORRELATIONS

The vector correlations between products and reagents for the reactions F + HH → HF + H , F + HD → HF + D and F + HT → HF + T , have been investigated using the quasi-classical trajectory (QCT) method on the Stark-Werner (SW) ab initio potential energy surface. The distribution P(θr) of angle between k′ and j′, the distribution P(ϕr) of dihedral angle denoting k-k′-j′ correlation are calculated. The polarization dependent generalized differential cross sections have also been studied. The evident influence of isotope substitution on the product polarization is revealed. This effect may be derived from the different mass factor of the three reactions and the barrier height of the F + H2 PES.

Quasiclassical trajectory study of the SiH(4)+H-->SiH(3)+H(2) reaction on a global ab initio potential energy surface

The SiH(4)+H-->SiH(3)+H(2) reaction has been investigated by the quasiclassical trajectory (QCT) method on a recent global ab initio potential energy surface [M. Wang et al., J. Chem. Phys. 124, 234311 (2006)]. The integral cross section as a function of collision energy and thermal rate coefficient for the temperature range of 300-1600 K have been obtained. At the collision energy of 9.41 kcalmol, product energy distributions and rovibrational populations are explored in detail, and H(2) rotational state distributions show a clear evidence of two reaction mechanisms. One is the conventional rebound mechanism and the other is the stripping mechanism similar to what has recently been found in the reaction of CD(4)+H [J. P. Camden et al., J. Am. Chem. Soc. 127, 11898 (2005)]. The computed rate coefficients with the zero-point energy correction are in good agreement with the available experimental data.

Interacting resonances in the F+H2 reaction revisited: Complex terms, Riemann surfaces, and angular distributions

We study the effect of overlapping resonances on the angular distributions of the reaction F+H2(v=0,j=0)-->HF(v=2,j=0)+H in the collision energy range from 5 to 65 meV, i.e., under the reaction barrier. Reactive scattering calculations were performed using the hyperquantization algorithm on the potential energy surface of Stark and Werner [J. Chem. Phys. 104, 6515 (1996)]. The positions of the Regge and complex energy poles are obtained by Pade reconstruction of the scattering matrix element. The Sturmian theory is invoked to relate the Regge and complex energy terms. For two interacting resonances, a two-sheet Riemann surface is contracted and inverted. The semiclassical complex angular momentum analysis is used to decompose the scattering amplitude into the direct and resonance contributions.

Cumulative reaction probabilities: A comparison between quasiclassical and quantum mechanical results

This article presents a quasiclassical trajectory (QCT) method for determining the cumulative reaction probability (CRP) as a function of the total energy. The method proposed is based on a discrete sampling using integer values of the total and orbital angular momentum quantum numbers for each trajectory and on the development of equations that have a clear counterpart in the quantum mechanical (QM) case. The calculations comprise cumulative reaction probabilities at a given total angular momentum J, as well as those summed over J. The latter are used to compute QCT rate constants. The method is illustrated by comparing QCT and exact QM results for the H+H2, H+D2, D+H2, and H+HD reactions. The agreement between QCT and QM results is very good, with small discrepancies between the two data sets indicating some genuine quantum effects. The most important of these involves the value of the CRP at low energies which, due to the absence of tunneling, is lower in the QCT calculations, causing the corresponding rate constants to be smaller. The second is the steplike structure that is clearly displayed in the QM CRP for J = 0, which is much smoother in the corresponding QCT results. However, when the QCT density of reactive states, i.e., the derivatives of the QCT CRP with respect to the energy, is calculated, a succession of maxima and minima is obtained which roughly resembles those found in the QM calculations, although the latter are considerably sharper. The analysis of the broad peaks in the QCT density of reactive states indicates that the distributions of collision times associated with the maxima are somewhat broader, with a tail extending to larger collision times, than those associated with the minima. In addition, the QM and QCT dynamics of the isotopic variants mentioned above are compared in the light of their CRPs. Issues such as the compliance of the QCT CRP with the law of microscopic reversibility, as well as the similarity between the CRPs for ortho and para species in the QM and QCT cases, are also addressed.

A crossed-beam study of the F+HD→DF+H reaction: The direct scattering channel

State-to-state differential cross sections of the title reaction are presented at four collision energies, ranging from 1.18 to 4.0 kcal /mol. Product angular distributions are predominantly backscattered at low energies and shift toward sideways (peaking near 150 degrees ) at higher energies. Experimental evidence for contributions from migratory trajectories was found in the more detailed angle-specific internal state distributions. The dynamics of this reaction is mostly governed by classical mechanics, and several major findings can qualitatively be rationalized. These "classical" behaviors serve as "references" and are to be contrasted to the attributes observed for the other isotopic product channel, HF+D, in a forthcoming paper.

Influence of rotation and isotope effects on the dynamics of the N(D2)+H2 reactive system and of its deuterated variants

Integral cross sections and thermal rate constants have been calculated for the N((2)D)+H(2) reaction and its isotopic variants N((2)D)+D(2) and the two-channel N((2)D)+HD by means of quasiclassical trajectory and statistical quantum-mechanical model methods on the latest ab initio potential-energy surface [T.-S. Ho et al., J. Chem. Phys. 119, 3063 (2003)]. The effect of rotational excitation of the diatom on the dynamics of these reactions has been investigated and interesting discrepancies between the classical and statistical model calculations have been found. Whereas a net effect of reagent rotation on reactivity is always observed in the classical calculations, only a very slight effect is observed in the case of the asymmetric N((2)D)+HD reaction for the statistical quantum-mechanical method. The thermal rate constants calculated on this Potential-Energy Surface using quasiclassical trajectory and statistical model methods are in good agreement with the experimental determinations, although the latter are somewhat larger. A reevaluation of the collinear barrier of the potential surface used in the present study seems timely. Further theoretical and experimental studies are needed for a full understanding of the dynamics of the title reaction.

Effect of rotational energy on the reaction Li+HF(υ=0,j)→LiF+H: An experimental and computational study

In a crossed molecular-beam study we have measured angular and time-of-flight distributions of the product LiF from the reaction Li + HF(upsilon = 0)-->LiF + H at various collision energies ranging from 97 to 363 meV for three markedly different rotational state distributions of HF obtained at nozzle temperatures close to 315, 510, and 850 K. Particularly, for the low and intermediate collision energies we observe significant effects of the varying j-state populations on the shape of the product angular distributions. At 315 K an additional feature appears in the angular distributions which is interpreted as being due to scattering from HF dimers. The experimental data are compared with simulations of the monomer reaction based on extensive quasiclassical trajectory calculations on a new state-of-the-art ab initio potential energy surface. We find an overall good agreement between the theoretical simulations and the experimental data for the title reaction, especially at the highest HF nozzle temperature.

Comparison of experimental time-of-flight spectra of the HF products from the F+H2 reaction with exact quantum mechanical calculations

High resolution HF product time-of-flight spectra measured for the reactive scattering of F atoms from n-H2(p-H2) molecules at collision energies between 69 and 81 meV are compared with exact coupled-channel quantum mechanical calculations based on the Stark-Werner ab initio ground state potential energy surface. Excellent agreement between the experimental and computed rotational distributions is found for the HF product vibrational states v'=1 and v'=2. For the v'=3 vibrational state the agreement, however, is less satisfactory, especially for the reaction with p-H2. The results for v'=1 and v'=2 confirm that the reaction dynamics for these product states is accurately described by the ground electronic state 1 (2)A' potential energy surface. The deviations for HF(v'=3, j' > or =2) are attributed to an enhancement of the reaction resulting from the 25% fraction of excited ((2)P(12)) fluorine atoms in the reactant beam.

Angular distributions for the F+H2→HF+H reaction: The role of the F spin-orbit excited state and comparison with molecular beam experiments

We report quantum mechanical calculations of center-of-mass differential cross sections (DCS) for the F+H(2)-->HF+H reaction performed on the multistate [Alexander-Stark-Werner (ASW)] potential energy surfaces (PES) that describe the open-shell character of this reaction. For comparison, we repeat single-state calculations with the Stark-Werner (SW) and Hartke-Stark-Werner (HSW) PESs. The ASW DCSs differ from those predicted for the SW and HSW PES in the backward direction. These differences arise from nonadiabatic coupling between several electronic states. The DCSs are then used in forward simulations of the laboratory-frame angular distributions (ADs) measured by Lee, Neumark, and co-workers [J. Chem. Phys. 82, 3045 (1985)]. The simulations are scaled to match experiment over the range 12 degrees <Theta(lab)<80 degrees. As a natural consequence of the reduced backward scattering, the ASW ADs are more forward and sideways scattered than predicted by the HSW PES. At the two higher collision energies (2.74 and 3.42 kcal/mol) the enhanced sideways scattering of HF v(')=2 products bring the ASW ADs in very good agreement with the experiment. At the lowest collision energy (1.84 kcal/mol), the simulations, for all three sets of PESs consistently underestimate the sideways scattering. The residual disagreements, particularly at the lowest collision energy, may be due to the known deficiencies in the PESs.

The investigation of spin–orbit effect for the F(2P)+HD reaction

In this paper, we employ the time-dependent quantum wave packet method to study the reaction of F((2)P(3/2), (2)P(1/2)) with HD on the Alexander-Stark-Werner potential energy surface. The reaction probabilities and total integral cross sections of the spin-orbit ground and excited states for the two possible products of the system are calculated. Because the reaction channel of the excited spin-orbit state is closed at the resonance energy, the resonance feature does not appear in the reaction probabilities and cross section for the F((2)P(1/2))+HD(v=j=0)-->HF+D reaction, in contrast with that found for the ground spin--orbit state. We also compare the average cross sections of the two possible products with the experimental measurement. The resonance peak in the present average cross section for the HF+D product is slightly larger than the experimental result, but much smaller than that of the single-state calculations on the potential energy surface of Stark and Werner. It seems that the spin--orbit coupling would play a relatively important role in this reaction. Moreover, the isotope effects of the ground and excited spin--orbit states and the reactivity of the two product channels from the excited spin--orbit state are presented.