Quantum Dynamics Study of the F + CH4 → HF + CH3 Reaction on an Ab Initio Potential Energy Surface

Role of the Vibrational and Translational Energies in the CN(v)+C2H6(ν1, ν2, ν5 and ν9) Reactions. A Theoretical QCT Study

Quasi-classical trajectory (QCT) calculations were conducted on the newly developed full-dimensional potential energy surface, PES-2023, to analyse two critical aspects: the influence of vibrational versus translational energy in promoting reactivity, and the impact of vibrational excitation within similar vibrational modes. The former relates to Polanyi's rules, while the latter concerns mode selectivity. Initially, the investigation revealed that independent vibrational excitation by a single quantum of ethane's symmetric and asymmetric stretching modes (differing by only 15 cm-1) yielded comparable dynamics, reaction cross-sections, HCN(v) vibrational product distributions, and scattering distributions. This observation dismisses any significant mode selectivity. Moreover, an equivalent amount of energy provided as translational energy (at total energies of 9.6 and 20.0 kcal mol-1) gave rise to slightly lower reactivity compared to the same amount of energy provided as vibrational energy. This effect is more evident at low energies, presenting a counterintuitive scenario in an 'early transition state' reaction. These findings challenge the straightforward application of Polanyi's rules in polyatomic systems. Regarding CN(v) vibrational excitation, our calculations reveal that the reaction cross-section remains practically unaffected by this vibrational excitation, suggesting that the CN stretching mode is a spectator mode. The results were rationalized by considering several factors: the strong coupling between different vibrational modes, and between vibrational modes and the reaction coordinate; and a significant vibrational energy redistribution within the ethane reactant before collision. This redistribution creates an unphysical energy flow, resulting in loss of adiabaticity and vibrational memory before the reactants' collision. These theoretical findings require future confirmation through experimental or theoretical quantum mechanical studies, which are currently unavailable.

Feshbach resonances in the F + CHD3 → HF + CD3 reaction

The signature of dynamics resonances was observed in the benchmark polyatomic F + CH₄/CHD₃ reactions more than a decade ago; however, the dynamical origin of the resonances is still not clear due to the lack of reliable quantum dynamics studies on accurate potential energy surfaces. Here, we report a six-dimensional state-to-state quantum dynamics study on the F + CHD₃ → HF + CD₃ reaction on a highly accurate potential energy surface. Pronounced oscillatory structures are observed in the total and product rovibrational-state-resolved reaction probabilities. Detailed analysis reveals that these oscillating features originate from the Feshbach resonance states trapped in the peculiar well on the HF(v' = 3)-CD₃ vibrationally adiabatic potential caused by HF chemical bond softening. Most of the resonance structures on the reaction probabilities are washed out in the well converged integral cross sections (ICS), leaving only one distinct peak at low collision energy. The calculated HF vibrational state-resolved ICS for CD₃(v = 0) agrees quantitatively with the experimental results, especially the branching ratio, but the theoretical CD₃ umbrella vibration state distribution is found to be much hotter than the experiment.

Theory-Based Mechanism for Fluoromethane Combustion I: Thermochemistry and Abstraction Reactions

A new detailed chemical kinetic mechanism is presented for small fluorinated hydrocarbons. Ab initio electronic structure theory is used to provide heats of formation with subchemical accuracy. The ANL0 method is extended to include fluorine. The resulting heats of formation at 0 K are in excellent agreement with 36 benchmark species in the Active Thermochemical Tables, with a mean error of μ = -0.02 kJ/mol and a standard deviation of σ = 0.91 kJ/mol. The thermophysical properties for 92 small-molecule H/C/O/F species are computed. The rate coefficients for 40+ H-abstraction reactions involving H, O, F, OH, OF, HO₂, and various methyl radicals with CH₄, CH₃F, CH₂F₂, CHF₃, CH₂O, and CHFO are discussed. The computed rate constants are in excellent agreement with the available literature. Additionally, 30+ rate constants are provided for F abstraction, which are several orders of magnitude smaller than H abstraction. The thermophysical properties and rate constants are provided in a mechanism. This mechanism is the first in a series of theory-based investigations into the thermal destruction of per- and polyfluorinated species.

Time-dependent quantum dynamics study of the F + C2H6 → HF + C2H5 reaction

The reaction probabilities, integral cross sections, energy efficiency and rate constants are investigated for the F + C₂H₆ reaction using the quantum reaction dynamics, wave packet method. The ground-state integral cross section calculated using a six-degree-of-freedom approach is in very good agreement with the quasi-classical trajectory results. We find that the H-CH₂CH₃ stretching motion has the largest enhancement to reactivity, followed by the H-CH₂-CH₃ bending motion. However, the stretching motion between CH₂ and CH₃ slightly hinders the reactivity. The energy-form efficacy based on an equal amount of total energy shows that translational energy is more effective in enhancing the reactivity than vibrational energy of the H-CH₂CH₃ stretching motion at a relatively lower translational energy, while the reverse is true at a relatively high translational energy. An energy-shifting method is employed to calculate the full-dimensional rate constants. The quantum rate constants agree well with one of the two main experimental measurements, and the activation energy has an excellent agreement with the one calculated using canonical variational transition-state theory.

Using quantum dynamics to study the effect of energy efficiency on the reactivity of the OH + DBr reaction

We report a time-dependent, full dimensional, wave-packet calculation for the reaction of OH + DBr to examine the effect of the energy efficiency on the reactivity. This study shows that the vibrational excitations of the OH and DBr enhance the reaction. However, the rotational excitations of OH and DBr both hinder the reaction. As a result, the vibrational energies of both the OH and DBr reactants are more efficient at promoting the reactivity than the translational energy, while the rotational energies of OH and DBr are less effective than the translational energy. By analyzing the state population of the vibrational and rotational states along the reaction pathway, we also developed an approach in order to explain the enhancement of the vibrational excitation and the hindrance of the rotational excitation of the reaction. We found that the initial-state selected vibrational excited states of OH and DBr are the dominant components, respectively, for surmounting the barrier. However, the initial-state selected rotational excited states of OH and DBr are no longer the dominant states for surmounting the transition state owing to their population changes in the van der Waals well. This quantitative analysis demonstrates the potential well in the entrance valley plays an important role in the energy efficiency with regards to the reactivity.

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.

From Reactive Rainbow to Dynamic Resonance Well

Scattering resonance is a fascinating phenomenon which manifests as a peak or a dip in an observable as a function of collisional energy (E_c). Its occurrence requires a potential well to support the resonance states. In this regard, reactive resonance is unusual in that it can exist in a reaction with unbound Born-Oppenheimer potential energy surface, on which the quasi-bound states are dynamically trapped-meaning that some energy is temporarily tied to the other degrees of freedom than the reaction coordinate. The concept of vibrational adiabaticity has been the cornerstone in understanding this phenomenon, for which the vibrationally adiabatic well depth is of primary concern. Recent studies on the F + CH₃D reaction have accumulated compelling evidence for a dominant resonance-mediated pathway at low E_c as well as for a rainbow feature in pair-correlated angular distribution at higher E_c. Here, we report an in-depth study to not only substantiate both claims but also, more importantly, make the first attempt to quantify the vibrationally adiabatic well depth directly from the observed rainbow structure and then compare with the theoretical prediction.

Imaging pair-correlated reaction cross sections in F + CH3D(νb = 0, 1) → CH2D(ν4 = 1) + HF(ν)

The title reactions were studied in a crossed-beam experiment at collisional energies (Ec) from 0.5 to 4.7 kcal mol-1. The νb (ν4) vibrational mode denotes the bending (umbrella) motion of the CH3D reactant (CH2D product). Using a time-sliced, velocity-map imaging technique, we extracted the state-specific, pair-correlated integral and differential cross sections. As with other isotopically analogous ground-state reactions, an inverted vibrational population of the HF coproduct was observed. Both the step-like excitation function near the threshold and the oscillatory forward-backward peakings in the Ec-evolution of the two dominant pair-correlated angular distributions at lower Ec suggest a resonance-mediated, time-delay mechanism. As Ec increases, the angular distribution of the HF(ν = 2) product evolves into a smooth and broad swath in the backward hemisphere, indicative of a direct rebound mechanism. One quantum excitation of the bending modes of CH3D(νb = 1) promotes the reaction rate by two- to three-fold up to Ec = 2.1 kcal mol-1. Broadly speaking, all major findings are qualitatively in line with previous results in the reactions of the F atom with other isotopologues. However, the rainbow feature recently observed in the CH2D(00) + HF(ν = 3) product channel is entirely absent. A possible rationale is put forward, which reinforces the previous reactive rainbow conjecture and calls for future theoretical investigations.

Understanding mode-specific dynamics in the local mode representation

Mode specificity is a main characteristic of transition state control of reaction dynamics. The normal mode representation has been widely employed to describe the mode specificity in elementary chemical reactions. However, spectroscopists have demonstrated that the local mode representation has advantages in analyzing the overtone and combination band spectra. In this work, the mode-specific reaction dynamics between the hydrogen atom and the molecules H2S and H2O is studied using a full-dimensional quantum scattering model in the (2 + 1) Radau-Jacobi coordinates. The mode specificities in the reactions that violates our physical intuition in the normal mode representation are well rationalized in the local mode representation. The energy flow between different XH bonds resulting from the intramolecular interaction and/or intermolecular interaction is unveiled, together with its impacts on dynamics of the abstraction and exchange reactions.

A neural network potential energy surface for the F + CH4reaction including multiple channels based on coupled cluster theory

We report here a new global and full dimensional potential energy surface (PES) for the F + CH₄reaction.

Recent advances in quantum scattering calculations on polyatomic bimolecular reactions

This review surveys quantum scattering calculations on chemical reactions of polyatomic molecules in the gas phase published in the last ten years.

Crossed beam polyatomic reaction dynamics: recent advances and new insights

This review summarizes the developments in polyatomic reaction dynamics, focusing on reactions of unsaturated hydrocarbons with O-atoms and methane with atoms/radicals.

Energy efficiency in surmounting the central energy barrier: a quantum dynamics study of the OH + CH3 → O + CH4 reaction

A quantum dynamics study of the OH + CH₃ with a slightly early barrier shows that vibrational energy is more effective in promoting the reactivity than translational energy, which is just opposite to the Polanyi rules.