Computational analyses of the vibrational spectra of fentanyl, carfentanil and remifentanil.

The infrared (IR) spectra of fentanyl, carfentanil and remifentanil, and protonated salts, are computed using quantum chemistry methods. New experimental FTIR spectra are also reported and compared to the calculations. The accuracy of two density functional theory methods, B3LYP and M06-2X, are tested against higher level theories (MP2) and the experimental data. Gas phase IR spectra are calculated for both the neutral and protonated molecules in order to compare with the experimental data measured for various salts of fentanyl and its analogues. Key vibrational modes are selected and studied in detail using a vibrational mode locality calculation. The main contributing atomic movements in these vibrational modes are identified.

Analytic Route to Tunneling Splittings Using Semiclassical Perturbation Theory.

We present an efficient, analytical, and simple route to approximating tunneling splittings in multidimensional chemical systems, directly from ab initio computations. The method is based on the Wentzel-Kramers-Brillouin (WKB) approximation combined with the vibrational perturbation theory. Anharmonicity and corner-cutting effects are implicitly accounted for without requiring a full potential energy surface. We test this method on the following three systems: a model one-dimensional double-well potential, the isomerization of malonaldehyde, and the isomerization of tropolone. The method is shown to be efficient and reliable.

Calculations on the unimolecular decomposition of the nerve agent VX.

It is very difficult to perform experiments on the physical parameters for the thermal decomposition of chemical nerve agents such as VX and computations, therefore, are useful. The reaction dynamics of the gas-phase pericyclic hydrogen transfer of the nerve agent VX is studied computationally. The geometries of the stationary structures are calculated at M06-2X/jul-cc-pVTZ level of theory. Single point energy calculations are carried out at the CBS/QB3 level to correct the energy barriers. Canonical reaction rate constants are calculated as a function of temperature. The one-dimensional semiclassical transition state theory is used to analyse the quantum tunneling effects. A reduced-dimensional hindered rotor model is proposed, tested, and applied to calculate the vibrational partition functions. It is found that the ester (O-side) and thioester (S-side) side chains of VX undergo pericyclic H-transfer reactions that result in decomposition of the molecule. The S-side reaction is favoured both kinetically and thermodynamically and dominates the pyrolysis over the temperature range from 600 K to 1000 K. It is predicted that VX completely decomposes in 2 s at temperatures above 750 K.

Hydrogen tunnelling in the rearrangements of carbenes: the role of dynamical calculations.

Tunnelling controlled chemical reactions are those which preferably proceed through pathways with high but narrow potential energy barriers, via quantum tunnelling, resulting in a product that would be disfavoured classically. These reactions are very sensitive to barrier width, height and temperature and so dynamical theoretical methods are required to describe these processes. Recent experimental work on charge-tagged phenyl pyruvic acid derivatives has found, in contrast to similar systems, no evidence of tunnelling control. Using semiclassical transition state theory, we rationalise these results and find tunnelling is significant in this system.

New Developments in Semiclassical Transition-State Theory.

This Feature Article describes some recent developments and applications of the Semiclassical Transition-State Theory (SCTST) for treating quantum tunneling in chemical reactions. A reduced dimensional form of the SCTST is discussed and is shown to be particularly efficient, as the required number of electronic structure calculations is reduced to an absolute minimum. We also describe how an alternative formulation of SCTST developed by Hernandez and Miller [ Chem. Phys. Lett. 1993 , 214 , 129 ], the SCTST-θ, has advantages in allowing for straightforward applications of the SCTST for any form of the potential expansion at the transition state. We also illustrate the power of SCTST in applications to more complex systems. We show how polyatomic modes such as internal rotations and torsions can be treated efficiently in SCTST calculations. We also describe some applications of the method to hydrogen atom tunneling in unimolecular reactions including the degradation of chemical nerve agents and the decay of the atmospherically important Criegee intermediates.

Theoretical Study of Gas-Phase Unimolecular Decomposition of Simulants of the Nerve Agent VX.

In order to further understand and support approaches for the degradation and destruction of toxic chemicals, the thermal decomposition of the nerve agent VX through possible pericyclic hydrogen transfer reactions is investigated using simulant molecules. A total of four simulant molecules are studied. Three of them have only one possible H-transfer site, while the other has two. They are chosen to bring physical insights into individual steps of the pericyclic reaction mechanism as well as the possible existence of competing mechanisms. The unimolecular reaction rate constants at the high-pressure limit are calculated. Geometries of stationary structures on the potential energy surfaces are calculated with the MP2 method as well as the B3LYP and M06-2X functionals and 6-311++G(d,p), jul-cc-pVTZ, and aug-cc-pVTZ basis sets. The barrier heights are corrected using energy values obtained at the CBS/QB3 level of theory. The contribution of the quantum tunneling effect to the reaction rate constants is included using one-dimensional semiclassical transition state theory. Adiabatic barrier heights, reaction rate constants, and branching ratio of the competing mechanisms are reported.

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. This includes quantum reactive scattering resonances, vibrational and rotational state effects, and quantum tunnelling in chemical reactions. Both experimental and theoretical advances are reviewed. Of particular emphasis is a description of the development of reduced dimensional theories which can highlight chemical reactions that are likely to be of interest for more accurate quantum reactive scattering studies and new experiments. Furthermore, the reduced dimensional models allow for the development and testing of computationally inexpensive procedures that enable calculations to be performed on quantum effects in reactions of larger polyatomic molecules.

Catalysis and tunnelling in the unimolecular decay of Criegee intermediates.

The unimolecular decay of Criegee intermediates is the major producer of OH radicals in the atmosphere. Here, Semi-Classical Transition State Theory (SCTST) in full and reduced dimensions is used to determine thermal rate constants for their unimolecular decay, as well as their decay catalysed by a single water molecule. These reactions shed light on the applicability of SCTST for catalysed hydrogen transfer reactions.

Application of one-dimensional semiclassical transition state theory to the CH3OH + H ⇌ CH2OH/CH3O + H2 reactions.

The rate constants of the two branches of H-abstractions from CH3OH by the H-atom and the corresponding reactions in the reverse direction are calculated using the one-dimensional semiclassical transition state theory (1D SCTST). In this method, only the reaction mode vibration of the transition state (TS) is treated anharmonically, while the remaining internal degrees of freedom are treated as they would have been in a standard TS theory calculation. A total of eight ab initio single-point energy calculations are performed in addition to the computational cost of a standard TS theory calculation. This allows a second-order Richardson extrapolation method to be employed to improve the numerical estimation of the third- and fourth-order derivatives, which in turn are used in the calculation of the anharmonic constant. Hindered-rotor (HR) vibrations are identified in the equilibrium states of CH3OH and CH2OH, and the TSs of the reactions. The partition function of the HRs are calculated using both a simple harmonic oscillator model and a more sophisticated one-dimensional torsional eigenvalue summation (1D TES) method. The 1D TES method can be easily adapted in 1D SCTST computation. The resulting 1D SCTST with 1D TES rate constants show good agreement to previous theoretical and experimental works. The effects of the HR on rate constants for different reactions are also investigated.This article is part of the theme issue 'Modern theoretical chemistry'.

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. These calculations are useful because they provide highly accurate information on the dynamics of chemical reactions which can be compared in detail with experimental results. They also serve as quantum mechanical benchmarks for testing approximate theories which can more readily be applied to more complicated reactions. This review includes theories for calculating quantities such as rate constants which have many important scientific applications.

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.

Rate constants of chemical reactions from semiclassical transition state theory in full and one dimension.

Semiclassical Transition State Theory (SCTST), a method for calculating rate constants of chemical reactions, offers gains in computational efficiency relative to more accurate quantum scattering methods. In full-dimensional (FD) SCTST, reaction probabilities are calculated from third and fourth potential derivatives along all vibrational degrees of freedom. However, the computational cost of FD SCTST scales unfavorably with system size, which prohibits its application to larger systems. In this study, the accuracy and efficiency of 1-D SCTST, in which only third and fourth derivatives along the reaction mode are used, are investigated in comparison to those of FD SCTST. Potential derivatives are obtained from numerical ab initio Hessian matrix calculations at the MP2/cc-pVTZ level of theory, and Richardson extrapolation is applied to improve the accuracy of these derivatives. Reaction barriers are calculated at the CCSD(T)/cc-pVTZ level. Results from FD SCTST agree with results from previous theoretical and experimental studies when Richardson extrapolation is applied. Results from our implementation of 1-D SCTST, which uses only 4 single-point MP2/cc-pVTZ energy calculations in addition to those for conventional TST, agree with FD results to within a factor of 5 at 250 K. This degree of agreement and the efficiency of the 1-D method suggest its potential as a means of approximating rate constants for systems too large for existing quantum scattering methods.

An investigation of one- versus two-dimensional semiclassical transition state theory for H atom abstraction and exchange reactions.

We investigate which terms in Reduced-Dimensionality Semiclassical Transition State Theory (RD SCTST) contribute most significantly in rate constant calculations of hydrogen extraction and exchange reactions of hydrocarbons. We also investigate the importance of deep tunneling corrections to the theory. In addition, we introduce a novel formulation of the theory in Jacobi coordinates. For the reactions of H atoms with methane, ethane, and cyclopropane, we find that a one-dimensional (1-D) version of the theory without deep tunneling corrections compares well with 2-D SCTST results and accurate quantum scattering results. For the "heavy-light-heavy" H atom exchange reaction between CH3 and CH4, deep tunneling corrections are needed to yield 1-D results that compare well with 2-D results. The finding that accurate rate constants can be obtained from derivatives of the potential along only one dimension further validates RD SCTST as a computationally efficient yet accurate rate constant theory.

Reduced-Dimensionality Semiclassical Transition State Theory: Application to Hydrogen Atom Abstraction and Exchange Reactions of Hydrocarbons.

Quantum mechanical methods for calculating rate constants are often intractable for reactions involving many atoms. Semiclassical transition state theory (SCTST) offers computational advantages over these methods but nonetheless scales exponentially with the number of degrees of freedom (DOFs) of the system. Here we present a method with more favorable scaling, reduced-dimensionality SCTST (RD SCTST), that treats only a subset of DOFs of the system explicitly. We apply it to three H abstraction and exchange reactions for which two-dimensional potential energy surfaces (PESs) have previously been constructed and evaluated using RD quantum scattering calculations. We differentiated these PESs to calculate harmonic frequencies and anharmonic constants, which were then used to calculate cumulative reaction probabilities and rate constants by RD SCTST. This method yielded rate constants in good agreement with quantum scattering results. Notably, it performed well for a heavy-light-heavy reaction, even though it does not explicitly account for corner-cutting effects. Recent extensions to SCTST that improve its treatment of deep tunneling were also evaluated within the reduced-dimensionality framework. The success of RD SCTST in this study suggests its potential applicability to larger systems.

Quantum dynamics of the abstraction reaction of H with cyclopropane.

The dynamics of the abstraction reaction of H atoms with the cyclopropane molecule is studied using quantum mechanical scattering theory. The quantum scattering calculations are performed in hyperspherical coordinates with a two-dimensional (2D) potential energy surface. The ab initio energy calculations are carried out with CCSD(T)-F12a/cc-pVTZ-F12 level of theory with the geometry and frequency calculations at the MP2/cc-pVTZ level. The contribution to the potential energy surface from the spectator modes is included as the projected zero-point energy correction to the ab initio energy. The 2D surface is fitted with a 29-parameter double Morse potential. An R-matrix propagation scheme is carried out to solve the close-coupled equations. The adiabatic energy barrier and reaction enthalpy are compared with high level computational calculations as well as experimental data. The calculated reaction rate constants shows very good agreement when compared with the experimental data, especially at lower temperature highlighting the importance of quantum tunnelling. The reaction probabilities are also presented and discussed. The special features of performing quantum dynamics calculation on the chemical reaction of a cyclic molecule are discussed.