New Method To Extract Final-State Information of Polyatomic Reactions Based on Normal Mode Analysis

Theoretical studies on the kinetics and dynamics of the BeH+ + H2O reaction: comparison with the experiment

The reaction of BeH+ with background gaseous H2O may play a role in qubit loss for quantum information processing with Be+ as trapped ions, and yet its reaction mechanism has not been well understood until now. In this work, a globally accurate, full-dimensional ground-state potential energy surface (PES) for the BeH+ + H2O reaction was constructed by fitting a total of 170 438 ab initio energy points at the level of RCCSD(T)-F12/aug-cc-pVTZ using the fundamental invariant-neural network method. The total root-mean-square error of the final PES was 0.178 kcal mol-1. For comparison, quasi-classical trajectory calculations were carried out on the PES at an experimental temperature of 150 K. The obtained thermal rate constant and product branching ratio of the BeD+ + H2O reaction agreed quite well with experimental results. In addition, the vibrational state distributions and energy disposals of the products were calculated and rationalized using the sudden vector projection model.

Dissection of the Multichannel Reaction O(3P) + C2H2: Differential Cross-Sections and Product Energy Distributions

The O(³P) + C₂H₂ reaction plays an important role in hydrocarbon combustion. It has two primary competing channels: H + HCCO (ketenyl) and CO + CH₂ (triplet methylene). To further understand the microscopic dynamic mechanism of this reaction, we report here a detailed quasi-classical trajectory study of the O(³P) + C₂H₂ reaction on the recently developed full-dimensional potential energy surface (PES). The entrance barrier TS1 is the rate-limiting barrier in the reaction. The translation of reactants can greatly promote reactivity, due to strong coupling with the reaction coordinate at TS1. The O(³P) + C₂H₂ reaction progress through a complex-forming mechanism, in which the intermediate HCCHO lives at least through the duration of a rotational period. The energy redistribution takes place during the creation of the long-lived high vibrationally (and rotationally) excited HCCHO in the reaction. The product energy partitioning of the two channels and CO vibrational distributions agree with experimental data, and the vibrational state distributions of all modes of products present a Boltzmann-like distribution.

Kinetic and Dynamic Studies of the F(2P) + ND3 → DF + ND2 Reaction

The fast F reaction with NH₃ poses a big challenge to experimental studies because of secondary chemical and collisional reactions. The quasi-classical trajectory method is utilized to investigate the mode specificity, product energy disposal, and temperature dependence of the thermal rate coefficient of F + ND₃ → DF + ND₂ on a recently developed potential energy surface. The effect of isotopic substitution is explored by comparing the F + ND₃ reaction with the F + NH₃ reaction. The computed results permit a better understanding of the F + ammonia reaction. The DF vibrational state has a Λ-type distribution, in accordance with the experimental measurement by the fast flow reactor technique. The product ND₂ is dominantly populated in the ground state, and a considerable amount of ND₂ is produced in the fundamental states of the bending mode. The similar vibrational state distributions of HF and NH₂ in the F + NH₃ reaction indicate a weak isotopic substitution effect on the product energy disposal. Exciting the umbrella mode of ND₃ suppresses the reaction at low energies below 5 kcal mol^-1, in sharp contrast to the observation in the F + NH₃ reaction. These dynamical behaviors can be partially explained by the sudden vector projection model. In addition, the thermal rate coefficient of F + ND₃ shows no temperature dependence in the range between 150 and 2000 K. There exists an inverse kinetic isotope effect at temperatures from 150 to 1500 K.

Effects of bending excitation on the reaction dynamics of fluorine atoms with ammonia

Vibrational excitation has been established as an efficient way to control the chemical reaction outcome. Stretching vibration of polyatomic molecules is believed to be efficient to promote abstraction reactions since energy is placed directly into the breaking bond. In this work, we report on a counterexample showing that exciting the low-frequency umbrella bending mode of ammonia enhances its reaction with fluorine atoms much more than exciting the high-frequency symmetric or asymmetric stretching mode over a wide range of collision energy, validated using both quasiclassical trajectory simulations and full-dimensional quantum dynamics calculations under the centrifugal-sudden approximation. This interesting mode-specific reaction dynamic originates from the increased chance of capturing the fluorine atom by ammonia due to the enlarged attractive interaction between them and the enhancement of the direct stripping reaction mediated by two submerged barriers.

Final-State-Resolved Dynamics of the H3+ + CO → H2 +HCO+/HOC+ Reaction: A Quasi-Classical Trajectory Study

The ion-molecule reaction H₃⁺ + CO → H₂ + HCO⁺/HOC⁺, which initiates the formation of crucial organic molecules, plays a key role in interstellar and circumstellar environments. In this work, the quasi-classical trajectory method is employed to study the reaction dynamics on a recently developed full-dimensional global potential energy surface (PES). The calculated product internal energy distributions and relative internal excited fractions agree reasonably well with the experimental measurements. For the two reaction channels, most of the available energy flows into the vibrational modes of HCO⁺ or HOC⁺ at low collision energies, followed by the translational mode and the rotational modes of HCO⁺ or HOC⁺. As the collision energy increases, the proportion of the product translational energy increases while the proportion of the product vibrational energy decreases. Furthermore, the CH and CO stretching modes and their combination bands are effectively excited for the product HCO⁺ while the bending mode is remarkably excited for the product HOC⁺.

When classical trajectories get to quantum accuracy: II. The scattering of rotationally excited H2 on Pd(111)

The classical trajectory method in a quantum spirit assigns statistical weights to classical paths on the basis of two semiclassical corrections: Gaussian binning and the adiabaticity correction.

When Classical Trajectories Get to Quantum Accuracy: The Scattering of H2 on Pd(111)

When elementary reactive processes occur at such low energies that only a few states of reactants and/or products are available, quantum effects strongly manifest and the standard description of the dynamics within the classical framework fails. We show here, for H₂ scattering on Pd(111), that by pseudoquantizing in the spirit of Bohr the relevant final actions of the system, along with adequately treating the diffraction-mediated trapping of the incoming wave, classical simulations achieve an unprecedented agreement with state-of-the-art quantum dynamics calculations.

Theoretical Study of Barrierless Chemical Reactions Involving Nearly Elastic Rebound: The Case of S(1D) + X2, X = H, D

For some values of the total angular momentum consistent with reaction, the title processes involve nonreactive trajectories proceeding through a single rebound mechanism during which the internal motion of the reagent diatom is nearly unperturbed. When such paths are in a significant amount, the classical reaction probability is found to be markedly lower than the quantum mechanical one. This finding was recently attributed to an unusual quantum effect called diffraction-mediated trapping, and a semiclassical correction was proposed in order to take into account this effect in the classical trajectory method. In the present work, we apply the resulting approach to the calculation of opacity functions as well as total and state-resolved integral cross sections (ICSs) and compare the values obtained with exact quantum ones, most of which are new. As the title reactions proceed through a deep insertion well, mean potential statistical calculations are also presented. Seven values of the collision energy, ranging from 30 to 1127 K, are considered. Two remarkable facts stand out: (i) The corrected classical treatment strongly improves the accuracy of the opacity function as compared to the usual classical treatment. When the entrance transition state is tight, however, those trajectories crossing it with a bending vibrational energy below the zero point energy must be discarded. (ii) The quantum opacity function, particularly its cutoff, is finely reproduced by the statistical approach. Consequently, the total ICS is also very well described by the two previous approximate methods. These, however, do not predict state-resolved ICSs with the same accuracy, proving thereby that (i) one or several genuine quantum effects involved in the dynamics are missed by the corrected classical treatment and (ii) the dynamics are not fully statistical.

Theoretical study of the F(2P) + NH3→ HF + NH2 reaction on an accurate potential energy surface: dynamics and kinetics

The highly exothermic hydrogen abstraction reaction of the F atom with NH₃ is investigated using the quasi-classical trajectory method on a newly developed potential energy surface (PES) for the ground electronic state. The full-dimensional PES is constructed by fitting 41 282 ab initio energy points at the level of UCCSD(T)-F12/aug-cc-pVTZ. The flexible fundamental invariant-neural network method is applied in the fitting, resulting in a total root mean square error of 0.13 kcal mol^-1. On one hand, the calculated differential cross sections agree reasonably well with the experimental results and indicate that the reaction is dominated by the direct abstraction and stripping mechanisms while a considerable amount of reaction takes place by the indirect "yo-yo" mechanism. The product energy partition also reproduces well the experimental result, which can be understood according to the geometry change along the minimum energy path. On the other hand, the obtained vibrational state distribution of the product HF follows P_νHF=2≈P_νHF=1 > P_νHF=0 > P_νHF=3, less consistent with the scattered experimental results. In addition, the calculated thermal rate coefficients have practically no temperature dependence within the statistical errors.

Bridging the Gap between Direct Dynamics and Globally Accurate Reactive Potential Energy Surfaces Using Neural Networks

Direct dynamics simulations become increasingly popular in studying reaction dynamics for complex systems where analytical potential energy surfaces (PESs) are unavailable. Yet, the number and/or the propagation time of trajectories are often limited by high computational costs, and numerous energies and forces generated on-the-fly become wasted after simulations. We demonstrate here an example of reusing only a very small portion of existing direct dynamics data to reconstruct a 90-dimensional globally accurate reactive PES describing the interaction of CO₂ with a movable Ni(100) surface based on a machine learning approach. In addition to reproducing previous results with much better statistics, we predict scattering probabilities of CO₂ at the state-to-state level, which is extremely demanding for direct dynamics. We propose this unified way to investigate gaseous and gas-surface reactions of medium size, initiating with hundreds of preliminary direct dynamics trajectories, followed by low-cost and high-quality simulations on full-dimensional analytical PESs.

Uncovering the role of the stationary points in the dynamics of the F- + CH3I reaction

We describe an analysis method which assigns geometries to stationary points along (quasi)classical trajectories. The method is applied to the F^- + CH₃I reaction, thereby uncovering the role of the minima and transition states in the dynamics of the S_N2 inversion, S_N2 retention via front-side attack and double inversion, induced inversion, and proton-transfer channels. Stationary-point probability distributions, stationary-point-specific trajectory orthogonal projections, root-mean-square distance distributions, transition probability matrices, and time evolutions of the stationary points reveal long-lived front-side (F^-ICH₃) and hydrogen-bonded (F^-HCH₂I) complexes in the entrance channel and significant post-reaction ion-dipole complex (FCH₃I^-) formation in the S_N2 exit channel. Most of the proton-transfer stationary points (FHCH₂I^-) participate in all the reaction channels with larger distance deviations than the double-inversion transition state. Significant forward-backward transitions are observed between the minima and transition states indicating complex, indirect dynamics. The utility of distance and energy constraints is also investigated, thereby restricting the assignment into uniform configuration or energy ranges around the stationary points.

Dynamics and kinetics of the reaction OH + H2S → H2O + SH on an accurate potential energy surface

Mode specificity and product energy disposal are unveiled in the reaction OH + H₂S → H₂O + SH.