Temperature and pressure dependent rate constants of the reactions of OH• with cyclopentene from variational TST and SS-QRRK methods.

In this work, the pressure- and temperature-dependent reaction rate constants for the hydrogen abstraction and addition of hydroxyl radicals to the unsaturated cyclopentene were studied. Geometries and vibrational frequencies of reactants, products, and transition states were calculated using density functional theory, with single-point energy corrections determined at the domain-based local pair natural orbital-coupled-cluster single double triple/cc-pVTZ-F12 level. The high-pressure limit rate constants were calculated using the canonical variational transition state theory with the small-curvature tunneling approximation. The vibrational partition functions were corrected by the effects of torsional and ring-puckering anharmonicities of the transition states and cyclopentene, respectively. Variational effects are shown to be relevant for all the hydrogen abstraction reactions. The increasing of the rate constants by tunneling is significant at temperatures below 500 K. The pressure dependence on the rate constants of the addition of OH• to cyclopentene was calculated using the system-specific quantum Rice-Ramsperger-Kassel model. The high-pressure limit rate constants decrease with increasing temperature in the range 250-1000 K. The falloff behavior was studied at several temperatures with pressures varying between 10-3 and 103 bar. At temperatures below 500 K, the effect of the pressure on the addition rate constant is very modest. However, at temperatures around and above 1000 K, taking pressure into account is mandatory for an accurate rate constant calculation. Branching ratio analyses reveal that the addition reaction dominates at temperatures below 500 K, decreasing rapidly at higher temperatures. Arrhenius parameters are provided for all reactions and pressure dependent Arrhenius parameters are given for the addition of OH• to cyclopentene.

Assessing the Molecular Basis of the Fuel Octane Scale: A Detailed Investigation on the Rate Controlling Steps of the Autoignition of Heptane and Isooctane.

N-Heptane and 2,2,4-trimethylpentane (isooctane) are the key species in the modeling of ignition of hydrocarbon-based fuel formulations. Isooctane is knock-resistant whereas n-heptane is a very knock-prone hydrocarbon. It has been suggested that interconversion of their associated alkylperoxy and hydroperoxyalkyl species via hydrogen-transfer isomerization reaction is the key step to understand their different knocking behavior. In this work, the kinetics of unimolecular hydrogen-transfer reactions of n-heptylperoxy and isooctylperoxy are determined using canonical variational transition-state theory and multidimensional small curvature tunneling. Internal rotation of involved molecules is taken explicitly into account in the molecular partition function. The rate coefficients are calculated in the temperature range 300-900 K, relevant to low-temperature autoignition. The concerted HO2 elimination is an important reaction that competes with some H-transfer and is associated with chain termination. Thus, the branching ratio between these reaction channels is analyzed. We show that variational and multidimensional tunneling effects cannot be neglected for the H-transfer reaction. In particular, the pre-exponential Arrhenius fitting parameter derived from our rate constants shows a strong dependence on the temperature, because tunneling increases quickly at temperatures below 500 K. On the basis of our results, the existing qualitative model for the reasons for different knock behavior observed for n-heptane and isooctane is quantitatively validated at the molecular level.

Chemical bonding and the equilibrium composition of Grignard reagents in ethereal solutions.

A thorough analysis of the electronic structure and thermodynamic aspects of Grignard reagents and its associated equilibrium composition in ethereal solutions is performed. Considering methylmagnesium halides containing fluorine, chlorine, and bromine, we studied the neutral, charged, and radical species associated with their chemical equilibrium in solution. The ethereal solvents considered, tetrahydrofuran (THF) and ethyl ether (Et(2)O), were modeled using the polarizable continuum model (PCM) and also by explicit coordination to the Mg atoms in a cluster. The chemical bonding of the species that constitute the Grignard reagent is analyzed in detail with generalized valence bond (GVB) wave functions. Equilibrium constants were calculated with the DFT/M06 functional and GVB wave functions, yielding similar results. According to our calculations and existing kinetic and electrochemical evidence, the species R(•), R(-), (•)MgX, and RMgX(2)(-) must be present in low concentration in the equilibrium. We conclude that depending on the halogen, a different route must be followed to produce the relevant equilibrium species in each case. Chloride and bromide must preferably follow a "radical-based" pathway, and fluoride must follow a "carbanionic-based" pathway. These different mechanisms are contrasted against the available experimental results and are proven to be consistent with the existing thermodynamic data on the Grignard reagent equilibria.

On the generalized valence bond description of the anomeric and exo-anomeric effects: an ab initio conformational study of 2-methoxytetrahydropyran.

An ab initio conformational study of the alpha- and beta-glycosidic C1-O1 bonds has been carried out on the axial and equatorial forms of 2-methoxytetrahydropyran (2-MTHP) at the HF/6-31G(d,p) and GVB-PP/6-31G(d,p) levels of calculation. Six conformers of 2-MTHP were fully optimized at both levels. The calculations have shown that the conformer containing the (+sc) orientation around the axial C1-O1 bond is the most stable one and is favored over that bearing the (-sc) arrangement about the equatorial C1-O1 bond by 1.39 (HF) and 1.52 (GVB-PP)kcal/mol. The potential energy surfaces for rotating about the axial and equatorial C1-O1 bonds were constructed at the HF and GVB-PP levels. For each form of 2-MTHP the HF and GVB-PP potential curves exhibit similar profiles. This shows that both methods provide similar descriptions for the position of the conformational minima and for the values and location of the rotational barriers. In addition to the conformational study, a discussion concerning the nature of the chemical bond in acetal fragments and the origin of the anomeric and exo-anomeric effects is presented in terms of optimized non-orthogonal GVB orbitals of 2-MTHP. The intramolecular factors that respond for the order of stability and conformational changes in bond lengths of the conformers of 2-MTHP are examined in light of the GVB description. The problems associated with the use of the NBOs (natural bond orbitals) to analyze chemical bonding in the acetal fragments are discussed, and the choice for the GVB-PP description is justified.

Nature of the chemical bond in protonated methane.

Protonated methane, CH(5)(+), is a key reactive intermediate in hydrocarbon chemistry and a borderline case for chemical structure theory, being the simplest example of hypercoordinated carbon. Early quantum mechanical calculations predicted that the properties of this species could not be associated with only one structure, because it presents serious limitations of the Born-Oppenheimer approximation. However, ab initio molecular dynamics and diffusion Monte Carlo calculations showed that the most populated structure could be pictured as a CH(3) tripod linked to a H(2) moiety. Despite this controversy, a model for the chemical bonds involved in this ion still lacks. Here we present a modern valence bond model for the electronic structure of CH(5)(+). The chemical bond scheme derived directly from our calculations pictures this ion as H(3)C...H(2)(+). The fluxionality can be seen as the result of a proton transfer between C-H bonds. A new insight on the vibrational bands at approximately 2400 and approximately 2700 cm(-1) is suggested. Our results show that the chemical bond model can be profitably applied to such intriguing systems.

Unified mechanistic concept of electrophilic aromatic nitration: convergence of computational results and experimental data.

The mechanism of electrophilic aromatic nitration was revisited. Based on the available experimental data and new high-level quantum chemical calculations, a modification of the previous reaction mechanism is proposed involving three separate intermediates on the potential energy diagram of the reaction. The first, originally considered an unoriented pi-complex or electron donor acceptor complex (EDA), involves high electrostatic and charge-transfer interactions between the nitronium ion and the pi-aromatics. It explains the observed low substrate selectivity in nitration with nitronium salts while maintaining high positional selectivity, as well as observed oxygen transfer reactions in the gas phase. The subsequent second intermediate originally considered an oriented "pi-complex" is now best represented by an intimate radical cation-molecule pair, C(6)H(6)(+)(*)()/NO(2), that is, a SET complex, indicative of single-electron transfer from the aromatic pi-system to NO(2)(+). Subsequently, it collapses to afford the final sigma-complex intermediate, that is, an arenium ion. The proposed three discrete intermediates in electrophilic aromatic nitration unify previous mechanistic proposals and also contribute to a better understanding of this fundamentally important reaction. The previously obtained ICR data of oxygen transfer from NO(2)(+) to the aromatic ring are also accommodated by the proposed mechanism. The most stable intermediate of this reaction on its potential energy surface is a complex between phenol and NO(+). The phenol.NO(+) complex decomposes affording C(6)H(6)O(+)(*)/PhOH(+) and NO, in agreement with the ICR results.