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.

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

A tunnelling controlled reaction is studied with semiclassical transition state theory, rationalising the results of experiment.

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.

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 CH₃OH 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 CH₃OH and CH₂OH, 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'.

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.

Catalysis and tunnelling in the unimolecular decay of Criegee intermediates

Semi-classical Transition State theory can be applied to catalysed atmospheric reactions, but reaction mode anharmonicity must be treated carefully.

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.

Chemical reaction dynamics

Guest editors Xueming Yang, David Clary and Daniel Neumark introduce the chemical reaction dynamics themed issue of Chemical Society Reviews.

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.

100 Years of Atomic Theory

One hundred years ago, Niels Bohr's pioneering paper on the electronic structure of the hydrogen atom revolutionized atomic theory.

Quantum effects in the abstraction reaction by H atoms of primary and secondary hydrogens in n-C4H10: a test of a new potential energy surface construction method

Recently, von Horsten et al. suggested an efficient method to construct two dimensional potential energy surfaces (PESs) for use in quantum scattering simulations, in which they utilised the minimum energy path (MEP) and assumed harmonic behaviour of the PES near the MEP. In the same paper, the authors applied this method to various H-abstraction reactions from C1-C3 alkane molecules. In this work we demonstrate an alternative PES construction method, and apply it to the more challenging H-abstraction from n-C(4)H(10) reactions. The geometry optimizations and frequency calculations are done at the MP2/cc-pVTZ level of theory, while the energies are calculated with the CCSD(T) method with the same basis set. The calculations give adiabatic energy barrier heights of 45.9 kJ mol(-1) and 34.4 kJ mol(-1) for the primary and secondary hydrogens in n-C(4)H(10). When compared to purely classical transition state theory, quantum scattering calculations show that quantum tunnelling and zero-point effects have large contributions at low temperatures, typically below 500 K. The branching ratio study suggests that the abstraction of secondary hydrogen in n-C(4)H(10) dominates the overall reaction rate at low temperatures. The rate constants for the two abstraction channels become more comparable as the temperature increases.

Theoretical studies on bimolecular reaction dynamics

This perspective discusses progress in the theory of bimolecular reaction dynamics in the gas phase. The examples selected show that definitive quantum dynamical computations are providing insights into the detailed mechanisms of chemical reactions.

Quantum study on the branching ratio of the reaction NO2+OH

A reduced dimensionality (RD) approximation is developed for the title reaction which treats the angle of approach of the hydroxyl radical to the nitrogen dioxide molecule and the radial distance between the two species explicitly. All other degrees of freedom are treated adiabatically. Electronic structure calculations at the complete active space self-consistent field level are used to fit a potential energy surface (PES) in these two coordinates. Within this RD model the adiabatic capture centrifugal sudden approximation is used to calculate the high pressure limit rate constant. A correction for reflection from the PES due to rotationally nonadiabatic transitions is applied using the wave packet capture approximation. The branching ratio for the title reaction is calculated for the atmospherically significant temperature range of 200-400 K at 20 Torr without distinguishing between the conformers of HOONO. The result is k(HOONO)k(HNO(3) )=0.051 at 20 Torr and 300 K, which is in good agreement with the measured branching ratio between cis-cis-HOONO and nitric acid. This suggests that most of the different conformers of HOONO were converted to the most stable cis-cis conformer on the time scale of the measurements made.

Torsional anharmonicity in transition state theory calculations

We present a new application for the Torsional Path Integral Monte Carlo (TPIMC) method in which the TPI partition functions are introduced into the calculation of Transition State Theory (TST) rate constants. In this way, an explicit treatment of torsional anharmonicity is included in the TST calculations and the magnitude of these effects can be assessed. The new method is tested on the C(2)H(6) + H hydrogen abstraction reaction and concerted hydrogen transfer in the carbonic acid dimer, for which we have developed torsional potential energy surfaces. For the C(2)H(6) + H reaction the rate constants are halved at room temperature on including a treatment of torsional anharmonicity, while the effects are found to be much smaller for the hydrogen transfer reaction in the carbonic acid dimer.

Reduced dimensionality quantum dynamics of Cl + CH4? HCl + CH3 on an ab initio potential

We study the reaction Cl + CH(4)--> HCl + CH(3) using a 2-D potential energy surface obtained by fitting a double Morse analytical function to high level (CCSD(T)/cc-pVTZ//MP2/cc-pVTZ)ab initio data. Dynamics simulations are performed in hyperspherical coordinates with the close-coupled equations being solved using R-matrix propagation. Quantum contributions from spectator modes are included via a harmonic zero-point correction to the ab initio data prior to fitting the potential. This is the first time this method has been applied to a heavy-light-heavy reaction and the first time it has been used to study differential cross sections. We find thermal rate constants and state-to-state differential cross sections which are in good agreement with experimental data. We discuss the applicability of our method to the study of kinetic isotope effects (KIEs), which we derive for the CH(4)/CD(4) substitution. The calculated KIE compares favourably with experiment. Finally, we discuss the sensitivity of the results of dynamics simulations on the accuracy of the fitted potential.

Torsional anharmonicity in the conformational analysis of tryptamine

In this paper we calculate the relative conformer populations of the tryptamine molecule. Our approach combines high level electronic structure conformer energies with harmonic frequencies and an anharmonic treatment of the torsional motions using the torsional path integral Monte Carlo method. We have developed a 3-D potential energy surface as a function of the torsional coordinates at the B3LYP/6-31+G(d) level using 2535 grid points. Eight conformers of tryptamine were found to be significantly populated at 430 K as opposed to the experimental observation of seven. This, along with further comparisons with various experimental data, leads us to suppose that conformer interconversion occurs during the cooling phases of many of the experiments. The ordering of the calculated populations fits well with available experimental data. Torsional anharmonicity is found to affect conformer populations more significantly at 430 K than at 100 K (although overall the effects are small), while quantum mechanical effects are not important at either temperature.

Predicting Catalysis: Understanding Ammonia Synthesis from First-Principles Calculations

Here, we give a full account of a large collaborative effort toward an atomic-scale understanding of modern industrial ammonia production over ruthenium catalysts. We show that overall rates of ammonia production can be determined by applying various levels of theory (including transition state theory with or without tunneling corrections, and quantum dynamics) to a range of relevant elementary reaction steps, such as N(2) dissociation, H(2) dissociation, and hydrogenation of the intermediate reactants. A complete kinetic model based on the most relevant elementary steps can be established for any given point along an industrial reactor, and the kinetic results can be integrated over the catalyst bed to determine the industrial reactor yield. We find that, given the present uncertainties, the rate of ammonia production is well-determined directly from our atomic-scale calculations. Furthermore, our studies provide new insight into several related fields, for instance, gas-phase and electrochemical ammonia synthesis. The success of predicting the outcome of a catalytic reaction from first-principles calculations supports our point of view that, in the future, theory will be a fully integrated tool in the search for the next generation of catalysts.

Quantum-mechanical calculations on pressure and temperature dependence of three-body recombination reactions: application to ozone formation rates

A quantum-mechanical model is designed for the calculation of termolecular association reaction rate coefficients in the low-pressure fall-off regime. The dynamics is set up within the energy transfer mechanism and the kinetic scheme is the steady-state approximation. We applied this model to the formation of ozone O + O2 + M --> O3 + M for M = Ar, making use of semiquantitative potential energy surfaces. The stabilization process is treated by means of the vibrational close-coupling infinite order sudden scattering theory. Major approximations include the neglect of the O3 vibrational bending mode and rovibrational couplings. We calculated individual isotope-specific rate constants and rate constant ratios over the temperature range 10-1000 K and the pressure fall-off region 10(-7)-10(2) bar. The present results show a qualitative and semiquantitative agreement with available experiments, particularly in the temperature region of atmospheric interest.