The Interplay between the Temperature and Pressure on the Reaction Pathways of the Prenol Oxidation by Hydroxyl Radicals.

Despite prenol emerging as a next-generation biofuel, some questions about its mechanism still need to be adequately proposed to rationalize its consumption and evaluate its efficiency in spark-ignition (SI) engines. Here, we present new insights into the reaction mechanism of prenol (3-methyl-2-buten-1-ol) with OH radicals as a function of temperature and pressure. We have determined that the different temperature and pressure conditions control the preferred products. At combustion temperatures and low pressures, OH-addition adducts are suppressed, increasing the formation of α and δ allylic radicals responsible for the auto-ignition.

Predicting pressure-dependent rate constants for the furan + OH reactions and their impact under tropospheric conditions.

In this work, we have evaluated the influence of temperature and pressure on the mechanism of furan oxidation by the OH radical. The stationary points on the potential energy surface were described at the M06-2X/aug-cc-pVTZ level of theory. In the kinetic treatment at the high-pressure limit (HPL), we have combined the multistructural canonical variational theory with multidimensional small-curvature tunneling corrections and long-range transition state theory. The system-specific quantum Rice-Ramsperger-Kassel theory was employed to estimate the pressure-dependent rate. In the HPL, the OH addition on the α carbon is the dominant pathway in the mechanism, producing a product via the ring-opening process, also confirmed by the product branching ratio calculations. The overall rate constant, obtained by a kinetic Monte Carlo simulation, reads the form koverall=5.22×10-13T/3001.10⁡exp1247(K/T) and indicates that the furan oxidation by OH radicals is a pressure-independent reaction under tropospheric conditions.

Understanding the reactivity and selectivity of Diels-Alder reactions involving furans.

The reactivity and endo/exo selectivity of the Diels-Alder cycloaddition reactions involving furan and substituted furans as dienes have been computationally explored. In comparison to cyclopentadiene, it is found that furan is comparatively less reactive and also less endo-selective in the reaction with maleic anhydride as the dienophile. Despite that, both the reactivity and the selectivity can be successfully modified by the presence of substituents at either 2- or 3-positions of the heterocycle. In this sense, it is found that the presence of strong electron-donor groups significantly increases the reactivity of the system while the opposite is found in the presence of electron-withdrawing groups. The observed trends in both the reactivity and selectivity are analyzed quantitatively in detail by means of the activation strain model of reactivity in combination with the energy decomposition analysis methods.

Multistructural partition function truncation and its effect on the thermal rate constants.

Thermal rate constants for the hydrogen abstraction reaction of methyl pentanoate were calculated using the multistructural canonical variational theory with small-curvature tunneling (MS-CVT/SCT). The conformational search for the stationary points generated by these reactions was performed with an algorithm that combines systematic and stochastic searches at dual-level. At the high-level (MPWB1K/6-31+G(d,p)), 244 geometries for methyl pentanoate and transition states were found. A systematic estimative of the efficiency of the truncation of multistructural rovibrational partition functions and its effects on thermal rate constants was carried out. In this analysis, we observed that rotation about -OCH3 group (ϕ4) of the ester generates unstable structures with little contribution to the magnitude of the partition function. We also estimate that the truncation of the partition functions in a small set of conformations can produce deviations up to 75% compared to converged theoretical results.

Prenol as a Next-Generation Biofuel or Additive: A Comprehension of the Hydrogen Abstraction Reactions by a H Atom.

Thermal rate coefficients for the hydrogen abstraction reactions of prenol (3-methyl-2-butenol) by a hydrogen atom were calculated with the multipath canonical variational theory with small-curvature tunneling (MP-CVT/SCT). The conformational search was performed with a dual-level approach, and the multistructural torsional anharmonicity effects were corrected through the rovibrational partition function calculated with the multistructural method based on a coupled torsional potential (MS-T(C)). This methodology allows us to estimate the thermal rate constants in the temperature range of 200-2500 K and fit them into two analytical expressions. Differences between the number of conformations on the torsional potential energy surfaces for prenol and the transition state decrease the thermal rate constants for the H-abstraction at the α carbon. An opposite behavior was detected for the abstractions on the δ site. The product branching ratios were calculated using single-structure and multipath approaches. The product distributions from the former are shown to be inadequate for studying the mechanism under combustion conditions. The values estimated from MP-CVT/SCT rate coefficients indicated that the radicals from (Rα) and (Rδ)/(Rδ') are formed in considerable amounts. These species are fundamental in comprehending the inhibition and promotion of the autoignition phenomena.

Differences in the torsional anharmonicity between reactant and transition state: the case of 3-butenal + H abstraction reactions.

Thermal rate coefficients for the hydrogen-abstraction reactions of 3-butenal by a hydrogen atom were obtained applying multipath canonical variational theory with small-curvature tunneling (MP-CVT/SCT). Torsional anharmonicity due to the hindered rotors was taken into account by calculating the rovibrational partition function using the extended two-dimensional torsional (E2DT) method. For comparison, rovibrational partition functions were also estimated using the multistructural method with torsional anharmonicity based on a coupled torsional potential (MS-T(C)). By contrast, with (E)-2-butenal reactions, the abstraction reactions of 3-butenal proceed via five reaction channels (R1)-(R5). In a conformational search, 45 distinguishable structures of transition states were found, including enantiomers, which were separated into six conformational reaction channels (CRCs). The individual reactive paths were constructed, the recrossing and semiclassical transmission coefficients estimated, and the multipath rate constants were obtained. High torsional barriers between the wells of CRC2/CRC6 indicate a harmonic behavior. Consequently, a difference between the torsional anharmonicity of 3-butenal and the transition states is responsible for the increase in the thermal rate constants for channel (R2). Analysis of the contributions of each conformer of the transition state shows an important contribution of the high-energy rotamers in the total flux of (R1)-(R5). After fitting the rate constants in a four-parameter equation, the activation energy estimation showed a strong temperature dependence.

Rate coefficients and product branching ratios for (E)-2-butenal + H reactions.

Thermal rate constants for the hydrogen abstraction reactions of (E)-2-butenal by hydrogen atoms were calculated, for the first time, using the multipath canonical variational theory with small-curvature tunneling (MP-CVT/SCT). After a torsional potential energy surface exploration, ten conformations of the transition states (including the mirror images) were found and separated into four conformational reaction channels (CRCs). Individual energy paths of each CRC were built, recrossing and quantum tunneling effects estimated, and the thermal rate constants obtained. Due to the hindered rotors, the torsional anharmonicity was incorporated in the rate coefficient through the calculations of the rovibrational partition functions using the extended two-dimensional torsional method (E2DT). For comparison, the one-well (1W-CVT/SCT) and harmonic multipath (MP-CVT/SCT) thermal rate constants were also estimated. In addition, kinetic Monte Carlo (KMC) simulations were performed to predict the product branching ratios. For all kinetic approaches, the formation of products of (R1) is predominant. Compared to the harmonic multipath estimation, the percentage of reaction (R4) increases by approximately 9% when the torsional anharmonicity is taken into account. For the reactions (R2) and (R3), the product branching ratio is slightly decreased when compared with the harmonic simulation.

Variational transition state theory rate constants and H/D kinetic isotope effects for CH3 + CH3 OCOH reactions.

The rate constants and H/D kinetic isotope effect for hydrogen abstraction reactions involving isotopomers of methyl formate by methyl radical are computed employing methods of the variational transition state theory (VTST) with multidimensional tunneling corrections. The energy paths were built with a dual-level method using the moller plesset second-order perturbation theory (MP2) method as the low-level and complete basis set (CBS) extrapolation as the high-level energy method. Benchmark calculations with the CBSD-T approach give an enthalpy of reaction at 0 K for R1 (-4.5 kcal/mol) and R2 (-4.2 kcal/mol) which are in good agreement with the experiment, that is, -4.0 and - 4.8 kcal/mol. For the reactional paths involving the isotopomers CH3 + CH3 OCOH → CH4 + CH3 OCO and CH3 + CH3 OCOD → CH3 D + CH3 OCO, the value of kH /kD (T = 455 K) using the canonical VTST/small-curvature tunneling approximation method is 6.7 in close agreement with experimental value (6.2). © 2019 Wiley Periodicals, Inc.

Anharmonicity of Coupled Torsions: The Extended Two-Dimensional Torsion Method and Its Use To Assess More Approximate Methods.

In this work we present the extended two-dimensional torsion (E2DT) method and use it to analyze the performance of several methods that incorporate torsional anharmonicity more approximately for calculating rotational-vibrational partition functions. Twenty molecules having two hindered rotors were studied for temperatures between 100 and 2500 K. These molecules present several kinds of situations; they include molecules with nearly separable rotors, molecules in which the reduced moments of inertia change substantially with the internal rotation, and molecules presenting compound rotation. Partition functions obtained by the rigid-rotor harmonic oscillator approximation, a method involving global separability of torsions and the multistructural methods without explicit potential coupling [MS-T(U)] and with explicit potential coupling [MS-T(C)] of torsions, are compared to those obtained with a quantized version - called the extended two-dimensional torsion (E2DT) method - of the extended hindered rotor approximation of Vansteenkiste et al. ( Vansteenkiste et al. J. Chem. Phys. 2006 , 124 , 044314 ). In the E2DT method, quantum effects due to the torsional modes were incorporated by the two-dimensional nonseparable method, which is a method that is based on the solution of the torsional Schrödinger equation and that includes full coupling in both the kinetic and potential energy. By comparing other methods to the E2DT method and to experimental thermochemical data, this study concludes that the harmonic approximation yields very poor results at high temperatures; the global separation of torsions from the rest of the degrees of freedom is not justified even when an accurate method to treat the torsions is employed; it is confirmed that methods based on less complete potential energy coupling of torsions, such as MS-T(U), are not accurate when dealing with rotors with different barrier heights, and more complete inclusion of torsional coupling to the method in MS-T(C) improves substantially the results in such a way that it could be used in cases where the E2DT method is unaffordable.

Hindered rotor tunneling splittings: an application of the two-dimensional non-separable method to benzyl alcohol and two of its fluorine derivatives.

In this work we present a novel application of the two-dimensional non-separable (2D-NS) method to the calculation of torsional tunneling splittings in systems with two hindered internal rotors. This method could be considered an extension of one-dimensional methods for the case of compounds with two tops. The 2D-NS method includes coupling between torsions in the kinetic and potential energy. Specifically, it has been applied to benzyl alcohol (BA) and two of its fluorine derivatives: 3-fluorobenzyl alcohol (3FBA) and 4-fluorobenzyl alcohol (4FBA). These molecules present two torsions, i.e., about the -CH2OH (ϕ1) and -OH (ϕ2) groups. The electronic structure calculations to build the two-dimensional torsional potential energy surface were performed at the DF-LMP2-F12//DF-LMP2/cc-pVQZ level of theory. For BA and 4FBA the calculated ground-state vibrational level splittings are 429 and 453 MHz, respectively, in good agreement with the experimental values of 337.10 and 492.82 MHz, respectively. In these two cases there are four equivalent wells and the tunneling splitting is the result of transitions between the two closer minima along ϕ1. The analysis of the wavefunctions, as well as the previous experimental work on the system, supports this conclusion. For 3FBA the observed ground-state splitting is 0.82 MHz, whereas in this case the calculated value amounts only to 0.02 MHz. The 2D-NS method, through the analysis of the wavefunctions, shows that this tiny tunneling splitting occurs between the two most stable minima of the potential energy surface. Additionally, we predict that the first vibrationally excited tunneling splitting will also be small and exclusively due to the interconversion between the second lowest minima.

Kinetic Isotope Effects in Multipath VTST: Application to a Hydrogen Abstraction Reaction.

In this work we apply multipath canonical variational transition state theory with small-tunneling corrections (MP-CVT/SCT) to the hydrogen abstraction reaction from ethanol by atomic hydrogen in aqueous solution at room temperature. This reaction presents two transition states which can interconvert by internal rotations about single bonds and another two transition states that are non-interconvertible enantiomers to the former structures. The study also includes another three reactions with isotopically substituted species for which there are experimental values of thermal rate constants and kinetic isotope effects (KIEs). The agreement between the MP-CVT/SCT thermal rate constants and the experimental data is good. The KIEs obtained by the MP-CVT/SCT methodology are factorized in terms of individual transition state contributions to facilitate the analysis. It was found that the percentage contribution of each transition state to the total KIE is independent of the isotopic substitution.

Exploring the electronic states of iodocarbyne: a theoretical contribution.

A manifold of electronic states correlating with the two lowest-lying dissociation channels of the iodocarbyne (CI) species is theoretically characterized for the first time in the literature. A contrast between the Λ + S and the relativistic (Ω) descriptions clearly shows the effect of perturbations on electronic states above 20 000 cm(-1) and the potential difficulties to detect them experimentally. For the bound states, spectroscopic parameters were evaluated, as well as the dipole moment functions. Similarly to CO, the polarity predicted for this iodocarbyne is C(δ-)I(δ+); as illustrated in the text, this is also the case for the other halocarbynes. As a potential mechanism for the experimental spectroscopic characterization of CI, we suggest the radiative association between C and I atoms, with light emitted in the red region of the visible spectra. Transition probabilities were also evaluated predicting very weak intensities. For the states 1/2(II) and 3/2(II), we have estimated radiative lifetimes of 7.1 and 714 ms, respectively.

Calcium-containing diatomic dications in the gas phase.

Sputtering (ion surface bombardment) of various calcium-containing powder samples with an energetic (17 keV), high-current (16)O(-) beam has produced the diatomic dications of CaSi(2+), CaP(2+), CaF(2+), CaH(2+), CaCl(2+), CaBr(2+) and CaI(2+). These molecular gas-phase species have been identified in positive ion mass spectra at half-integer m/z values; their ion flight times through a magnetic-sector mass spectrometer were roughly 10(-5) s. Most of them appear to be novel molecular ions; the stability of the latter four (CaH(2+), CaCl(2+), CaBr(2+) and CaI(2+)) had been demonstrated in previous theoretical studies, whereas only CaF(2+) and CaBr(2+) had been observed before. Here we combine the results of our experimental search with a detailed theoretical study of the remaining three systems CaSi(2+), CaP(2+) and CaF(2+). All electronic states correlating with the first dissociation channel are characterized using high level ab initio electronic structure calculations. In their ground states, we find CaSi(2+) to be a long-lived metastable molecule, whereas CaF(2+) and CaP(2+) are thermodynamically stable, with respective equilibrium internuclear distances of 6.253, 4.740, and 5.731 a(0). CaSi(2+) has a well depth of 7116 (0.88) cm(-1) (eV) and a dissociation asymptote 7956 (0.99) cm(-1) (eV) below the ground state minimum. The dissociation energy of CaF(2+) is estimated to be 3404 (0.42) cm(-1) (eV), whereas for CaP(2+) we found 2547 (0.32) cm(-1) (eV), and a barrier height of 8118 (1.01) cm(-1) (eV). Their adiabatic double ionisation energies are 22.87, 16.91, and 17.32 eV, respectively, for the F, Si, and P containing dications.