Ab Initio Reaction Rate Constants Computed Using Semiclassical Transition-State Theory: HO + H2 → H2O + H and Isotopologues

Spiers Memorial Lecture : Introductory lecture: quantum dynamics of chemical reactions

This Spiers Memorial Lecture discusses quantum effects that can be calculated and observed in the chemical reactions of small molecules.

Theoretical Studies on the Kinetics of Multi-Channel Gas-Phase Unimolecular Decomposition of Acetaldehyde

Theoretical kinetic studies are performed on the multichannel thermal decomposition of acetaldehyde. The geometries of the stationary points on the potential energy surface of the reaction are optimized at the MP2(full)/6-311++G(2d,2p) level of theory. More accurate energies are obtained by single point energy calculations at the CCSD(T,full)/augh-cc-pVTZ+2df, CBS-Q and G4 levels of theory. Here, by application of steady-state approximation to the thermally activated species CH₃CHO* and CH₂CHOH* and performance of statistical mechanical manipulations, expressions for the rate constants for different product channels are derived. Special attempts are made to compute accurate energy-specific rate coefficients for different channels by using semiclassical transition state theory. It is found that the isomerization of CH₃CHO to the enol-form CH₂CHOH plays a significant role in the unimolecular reaction of CH₃CHO. The possible products of the reaction are formed via unimolecular decomposition of CH₃CHO and CH₂CHOH. The computed rate coefficients reveal that the dominant channel at low temperatures and high pressures is the formation of CH₂CHOH due to the low barrier height for CH₃CHO → CH₂CHOH isomerization process. However, at high temperatures, the product channel CH₃ + CHO becomes dominant.

A Combined Theoretical and Experimental Study of Sarin (GB) Decomposition at High Temperatures

Theoretical and experimental results are presented for the pyrolytic decomposition of the nerve agent sarin (GB) in the gas phase. High-level quantum chemistry calculations are performed together with a semiclassical transition-state theory for describing quantum mechanical tunneling. The experimental and theoretical results for the temperature dependence of the survival times show very good agreement, as does the calculated and measured activation energy for thermal decomposition. The combined results suggest that the thermal decomposition of GB, for temperature ranging from 350 to 500 °C, goes through a pericyclic reaction mechanism with a transition state consisting of a six-membered ring structure.

How does substitution affect the unimolecular reaction rates of Criegee intermediates?

Unimolecular reaction rates of Criegee intermediates show substitution effect.

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.

A ring polymer molecular dynamics study of the OH + H2(D2) reaction

Using ring polymer molecular dynamics we have calculated the rate coefficients for the OH + H₂ reaction.

Semiclassical Transition-State Theory Based on Fourth-Order Vibrational Perturbation Theory: The Symmetrical Eckart Barrier

Semiclassical transition-state theory based on fourth-order vibrational perturbation theory (VPT4-SCTST) is applied to compute the barrier transmission coefficient for the symmetric Eckart potential. For a barrier parametrized to mimic the H2 + H exchange reaction, the results obtained are in excellent agreement with exact quantum calculations over a range of energy that extends down to roughly 1% of the barrier height, V0, where tunneling is negligible. The VPT2-SCTST treatment, which is commonly used in chemical kinetics studies, also performs quite well but already shows an error of a few percent at ca. 0.8 V0 where tunneling is still important. This suggests that VPT4-SCTST could offer an improvement over VPT2-SCTST in applications studies. However, the computational effort for VPT4-SCTST treatments of molecules is excessive, and any improvement gained is unlikely to warrant the increased effort. Nevertheless, the treatment of the symmetric Eckart barrier problem here suggests a simple modification of the usual VPT2-SCTST protocol that warrants further investigation.

Theoretical studies on the kinetics of hydrogen abstraction reactions of H and CH3 radicals from CH3OCH3 and some of their H/D isotopologues

The hydrogen abstraction reactions by H and CH3 radicals from CH3OCH3 and some of their H/D isotopologues are studied by semiclassical transition state theory. Many high-level density functional, ab initio, and combinatory quantum chemical methods, including B3LYP, BB1K, MP2, MP4, CCSD(T), CBS-Q, and G4 methods, are employed to compute the energies and rovibrational properties of the stationary points for the title reactions. Xij vibrational anharmonicity coefficients, used in semiclassical transition state theory, are computed at the B3LYP, BB1K, and MP2 levels of theory. Thermal rate coefficients and kinetic isotope effects are computed over the temperature range from 200 to 2500 K and compared with available experimental data. The computed rate constants for the title reactions are represented as the equation k(T) = ATn exp[−E(T + T0)/(T2 + T02)].

Instanton and noninstanton tunneling in periodically perturbed barriers: semiclassical and quantum interpretations

In multidimensional barrier tunneling, there exist two different types of tunneling mechanisms, instanton-type tunneling and noninstanton tunneling. In this paper we investigate transitions between the two tunneling mechanisms from the semiclassical and quantum viewpoints taking two simple models: a periodically perturbed Eckart barrier for the semiclassical analysis and a periodically perturbed rectangular barrier for the quantum analysis. As a result, similar transitions are observed with change of the perturbation frequency ω for both systems, and we obtain a comprehensive scenario from both semiclassical and quantum viewpoints for them. In the middle range of ω, in which the plateau spectrum is observed, noninstanton tunneling dominates the tunneling process, and the tunneling amplitude takes the maximum value. Noninstanton tunneling explained by stable-unstable manifold guided tunneling (SUMGT) from the semiclassical viewpoint is interpreted as multiphoton-assisted tunneling from the quantum viewpoint. However, in the limit ω→0, instanton-type tunneling takes the place of noninstanton tunneling, and the tunneling amplitude converges on a constant value depending on the perturbation strength. The spectrum localized around the input energy is observed, and there is a scaling law with respect to the width of the spectrum envelope, i.e., the width ∝ℏω. In the limit ω→∞, the tunneling amplitude converges on that of the unperturbed system, i.e., the instanton of the unperturbed system.

Reaction of HO with CO: Tunneling Is Indeed Important

The potential energy surface and chemical kinetics for the reaction of HO with CO, which is an important process in both combustion and atmospheric chemistry, were computed using high-level ab initio quantum chemistry in conjunction with semiclassical transition state theory under the limiting cases of high and zero pressure. The reaction rate constants calculated from first principles agree extremely well with all available experimental data, which range in temperature over a domain that covers both combustion and terrestrial atmospheric chemistry. The role of quantum tunneling is confirmed to be extremely important, which supports recent work by Continetti and collaborators regarding the loss of hydrogen atoms from vibrationally excited states of HOCO. A sensitivity analysis has been carried out and serves as the basis for a plausible estimate of uncertainty in the calculations.