Rubén Meana-Pañeda

New Pathways for Formation of Acids and Carbonyl Products in Low-Temperature Oxidation: The Korcek Decomposition of γ-Ketohydroperoxides

We present new reaction pathways relevant to low-temperature oxidation in gaseous and condensed phases. The new pathways originate from γ-ketohydroperoxides (KHP), which are well-known products in low-temperature oxidation and are assumed to react only via homolytic O-O dissociation in existing kinetic models. Our ab initio calculations identify new exothermic reactions of KHP forming a cyclic peroxide isomer, which decomposes via novel concerted reactions into carbonyl and carboxylic acid products. Geometries and frequencies of all stationary points are obtained using the M06-2X/MG3S DFT model chemistry, and energies are refined using RCCSD(T)-F12a/cc-pVTZ-F12 single-point calculations. Thermal rate coefficients are computed using variational transition-state theory (VTST) calculations with multidimensional tunneling contributions based on small-curvature tunneling (SCT). These are combined with multistructural partition functions (Q(MS-T)) to obtain direct dynamics multipath (MP-VTST/SCT) gas-phase rate coefficients. For comparison with liquid-phase measurements, solvent effects are included using continuum dielectric solvation models. The predicted rate coefficients are found to be in excellent agreement with experiment when due consideration is made for acid-catalyzed isomerization. This work provides theoretical confirmation of the 30-year-old hypothesis of Korcek and co-workers that KHPs are precursors to carboxylic acid formation, resolving an open problem in the kinetics of liquid-phase autoxidation. The significance of the new pathways in atmospheric chemistry, low-temperature combustion, and oxidation of biological lipids are discussed.

Least-Action Tunneling Transmission Coefficient for Polyatomic Reactions.

We present a new least-action variational approximation for tunneling in polyatomic reactions based on the procedure developed by Garrett and Truhlar for atom-diatom reactions. (63) The method calculates the semiclassical ground-state tunneling probability at every tunneling energy by minimizing the value of imaginary action integral along a family of paths ranging from the minimum energy path to the straight path. The method is illustrated by applications to two hydrogen-atom abstraction reactions from methane using analytical potential energy surfaces.

High-level direct-dynamics variational transition state theory calculations including multidimensional tunneling of the thermal rate constants, branching ratios, and kinetic isotope effects of the hydrogen abstraction reactions from methanol by atomic hydrogen

We report a detailed theoretical study of the hydrogen abstraction reaction from methanol by atomic hydrogen. The study includes the analysis of thermal rate constants, branching ratios, and kinetic isotope effects. Specifically, we have performed high-level computations at the MC3BB level together with direct dynamics calculations by canonical variational transition state theory (CVT) with the microcanonically optimized multidimensional tunneling (μOMT) transmission coefficient (CVT/μOMT) to study both the CH(3)OH+H→CH(2)OH+H(2) (R1) reaction and the CH(3)OH+H→CH(3)O+H(2) (R2) reaction. The CVT/μOMT calculations show that reaction R1 dominates in the whole range 298≤T (K)≤2500 and that anharmonic effects on the torsional mode about the C-O bond are important, mainly at high temperatures. The activation energy for the total reaction sum of R1 and R2 reactions changes substantially with temperature and, therefore, the use of straight-line Arrhenius plots is not valid. We recommend the use of new expressions for the total R1 + R2 reaction and for the R1 and R2 individual reactions.

MSTor version 2013: A new version of the computer code for the multi-structural torsional anharmonicity, now with a coupled torsional potential

We present an improved version of the MSTor program package, which calculates partition functions and thermodynamic functions of complex molecules involving multiple torsions; the method is based on either a coupled torsional potential or an uncoupled torsional potential. The program can also carry out calculations in the multiple-structure local harmonic approximation. The program package also includes seven utility codes that can be used as stand-alone programs to calculate reduced moment of inertia matrices by the method of Kilpatrick and Pitzer, to generate conformational structures, to calculate, either analytically or by Monte Carlo sampling, volumes for torsional subdomains defined by Voronoi tessellation of the conformational subspace, to generate template input files for the MSTor calculation and Voronoi calculation, and to calculate one-dimensional torsional partition functions using the torsional eigenvalue summation method.

Potential energy surface of triplet N2O2

We present a global ground-state triplet potential energy surface for the N2O2 system that is suitable for treating high-energy vibrational-rotational energy transfer and collision-induced dissociation. The surface is based on multi-state complete-active-space second-order perturbation theory/minimally augmented correlation-consistent polarized valence triple-zeta electronic structure calculations plus dynamically scaled external correlation. In the multireference calculations, the active space has 14 electrons in 12 orbitals. The calculations cover nine arrangements corresponding to dissociative diatom-diatom collisions of N2, O2, and nitric oxide (NO), the interaction of a triatomic molecule (N2O and NO2) with the fourth atom, and the interaction of a diatomic molecule with a single atom (i.e., the triatomic subsystems). The global ground-state potential energy surface was obtained by fitting the many-body interaction to 54 889 electronic structure data points with a fitting function that is a permutationally invariant polynomial in terms of bond-order functions of the six interatomic distances.

An Ancient Fingerprint Indicates the Common Ancestry of Rossmann-Fold Enzymes Utilizing Different Ribose-Based Cofactors

Nucleoside-based cofactors are presumed to have preceded proteins. The Rossmann fold is one of the most ancient and functionally diverse protein folds, and most Rossmann enzymes utilize nucleoside-based cofactors. We analyzed an omnipresent Rossmann ribose-binding interaction: a carboxylate side chain at the tip of the second β-strand (β2-Asp/Glu). We identified a canonical motif, defined by the β2-topology and unique geometry. The latter relates to the interaction being bidentate (both ribose hydroxyls interacting with the carboxylate oxygens), to the angle between the carboxylate and the ribose, and to the ribose’s ring configuration. We found that this canonical motif exhibits hallmarks of divergence rather than convergence. It is uniquely found in Rossmann enzymes that use different cofactors, primarily SAM (S-adenosyl methionine), NAD (nicotinamide adenine dinucleotide), and FAD (flavin adenine dinucleotide). Ribose-carboxylate bidentate interactions in other folds are not only rare but also have a different topology and geometry. We further show that the canonical geometry is not dictated by a physical constraint—geometries found in noncanonical interactions have similar calculated bond energies. Overall, these data indicate the divergence of several major Rossmann-fold enzyme classes, with different cofactors and catalytic chemistries, from a common pre-LUCA (last universal common ancestor) ancestor that possessed the β2-Asp/Glu motif.

Chloroform as a Hydrogen Atom Donor in Barton Reductive Decarboxylation Reactions

The utility of chloroform as both a solvent and a hydrogen atom donor in Barton reductive decarboxylation of a range of carboxylic acids was recently demonstrated (Ko, E. J. et al. Org. Lett. 2011, 13, 1944). In the present work, a combination of electronic structure calculations, direct dynamics calculations, and experimental studies was carried out to investigate how chloroform acts as a hydrogen atom donor in Barton reductive decarboxylations and to determine the scope of this process. The results from this study show that hydrogen atom transfer from chloroform occurs directly under kinetic control and is aided by a combination of polar effects and quantum mechanical tunneling. Chloroform acts as an effective hydrogen atom donor for primary, secondary, and tertiary alkyl radicals, although significant chlorination was also observed with unstrained tertiary carboxylic acids.

Prediction of experimentally unavailable product branching ratios for biofuel combustion: The role of anharmonicity in the reaction of isobutanol with OH

Isobutanol is a prototype biofuel, and sorting out the mechanism of its combustion is an important objective where theoretical modeling can provide information that is unavailable and not easily obtained by experiment. In the present work the rate constants and branching ratios for the hydrogen abstraction reactions from isobutanol by hydroxyl radical have been calculated using multi-path variational transition-state theory with small-curvature tunneling. We use hybrid degeneracy-corrected vibrational perturbation theory to show that it is critical to consider the anharmonicity difference of high-frequency modes between reactants and transition states. To obtain accurate rate constants, we must apply different scaling factors to the calculated harmonic vibrational frequencies at the reactants and at the transition states. The factors determining the reaction rate constants have been analyzed in detail, including variational effects, tunneling contributions, the effect of multiple reaction paths on transmission coefficients, and anharmonicities of low- and high-frequency vibrational modes. The analysis quantifies the uncertainties in the rate calculations. A key result of the paper is a prediction for the site dependence of hydrogen abstraction from isobutanol by hydroxyl radical. This is very hard to measure experimentally, although it is critical for combustion mechanism modeling. The present prediction differs considerably from previous theoretical work.

Tunneling and Conformational Flexibility Play Critical Roles in the Isomerization Mechanism of Vitamin D

The thermal isomerization reaction converting previtamin D to vitamin D is an intramolecular [1,7]-sigmatropic hydrogen shift with antarafacial stereochemistry. We have studied the dynamics of this reaction by means of the variational transition-state theory with multidimensional corrections for tunneling in both gas-phase and n-hexane environments. Two issues that may have important effects on the dynamics were analyzed in depth, i.e., the conformations of previtamin D and the quantum effects associated with the hydrogen-transfer reaction. Of the large number of conformers of previtamin D that were located, there are 16 that have the right disposition to react. The transition-state structures associated with these reaction paths are very close in energy, so all of them should be taken into account for an accurate calculation of both the thermal rate constants and the kinetic isotope effects. This issue is particularly important because the contribution of each of the reaction paths to the total thermal rate constant is quite sensitive to the environment. The dynamics results confirm that tunneling plays an important role and that model systems that were considered previously to study the hydrogen shift reaction cannot mimic the complexity introduced by the flexibility of the rings of previtamin D. Finally, the characterization of the conformers of both previtamin D and vitamin D allowed the calculation of the thermal equilibrium constants of the isomerization process.

Mechanisms of Double Proton Transfer. Theory and Applications

Abstract An analytical two-dimensional (2D) potential-energy surface based on two equal hydrogen bonds coupled by a correlation term, recently introduced [J. Chem. Phys. 127 (2007) 174513] to describe the dynamics of double proton transfer, is reviewed and generalized. It is then applied to the evaluation of proton transfer dynamics in a number of realistic systems, namely several molecules and dimers that exhibit various degrees of correlation between the motions of the two protons. The three parameters required to generate this 2D potential are derived from electronic structure and force field calculations, such that they include implicitly the effect of coupled skeletal modes. It follows that explicit introduction of such coupled modes is not required to obtain the basic relations that define the stationary points of the 2D surface, and thereby the reaction mechanism. Based on these relations, a detailed analysis is reported of a variety of systems exhibiting double proton transfer, including, apart from previously investigated porphine and porphycene, representing weak correlation, and the formic and benzoic acid dimers, representing strong correlation, two newly investigated systems which shed light on the hitherto not represented intermediate correlation category, namely naphthazarin, and the 4-bromopyrazole dimer.