Enthalpy of the gas-phase CO2 + Mg reaction from ab initio total energies

Various highly accurate ab initio composite methods of Gaussian-n (G1, G2, G3), their variations (G2(MP2), G3(MP2), G3//B3LYP, G3(MP2)//B3LYP), and complete basis set (CBS-Q, CBS-Q//B3LYP) series of models were applied to compute reaction enthalpies of the ground-state reaction of CO2 with Mg. All model chemistries predict highly endothermic reactions, with DeltaH(298) = 63.6-69.7 kcal x mol(-1). The difference between the calculated reaction enthalpies and the experimental value, evaluated with recommended experimental standard enthalpies of formation for products and reactants, is more than 20 kcal x mol(-1) for all methods. This difference originates in the incorrect experimental enthalpy of formation of gaseous MgO given in thermochemical databases. When the theoretical formation enthalpy for MgO calculated by a particular method is used, the deviation is reduced to 1.3 kcal x mol(-1). The performance of the methodologies used to calculate the heat of this particular reaction and the enthalpy of formation of MgO are discussed.     

Negative-ion photoelectron spectroscopy, gas-phase acidity, and thermochemistry of the peroxyl radicals CH(3)OO and CH(3)CH(2)OO

Methyl, methyl-d(3), and ethyl hydroperoxide anions (CH(3)OO(-), CD(3)OO(-), and CH(3)CH(2)OO(-)) have been prepared by deprotonation of their respective hydroperoxides in a stream of helium buffer gas. Photodetachment with 364 nm (3.408 eV) radiation was used to measure the adiabatic electron affinities: EA[CH(3)OO, X(2)A' '] = 1.161 +/- 0.005 eV, EA[CD(3)OO, X(2)A' '] = 1.154 +/- 0.004 eV, and EA[CH(3)CH(2)OO, X(2)A' '] = 1.186 +/- 0.004 eV. The photoelectron spectra yield values for the term energies: Delta E(X(2)A' '-A (2)A')[CH(3)OO] = 0.914 +/- 0.005 eV, Delta E(X(2)A' '-A (2)A')[CD(3)OO] = 0.913 +/- 0.004 eV, and Delta E(X(2)A' '-A (2)A')[CH(3)CH(2)OO] = 0.938 +/- 0.004 eV. A localized RO-O stretching mode was observed near 1100 cm(-1) for the ground state of all three radicals, and low-frequency R-O-O bending modes are also reported. Proton-transfer kinetics of the hydroperoxides have been measured in a tandem flowing afterglow-selected ion flow tube (FA-SIFT) to determine the gas-phase acidity of the parent hydroperoxides: Delta(acid)G(298)(CH(3)OOH) = 367.6 +/- 0.7 kcal mol(-1), Delta(acid)G(298)(CD(3)OOH) = 367.9 +/- 0.9 kcal mol(-1), and Delta(acid)G(298)(CH(3)CH(2)OOH) = 363.9 +/- 2.0 kcal mol(-1). From these acidities we have derived the enthalpies of deprotonation: Delta(acid)H(298)(CH(3)OOH) = 374.6 +/- 1.0 kcal mol(-1), Delta(acid)H(298)(CD(3)OOH) = 374.9 +/- 1.1 kcal mol(-1), and Delta(acid)H(298)(CH(3)CH(2)OOH) = 371.0 +/- 2.2 kcal mol(-1). Use of the negative-ion acidity/EA cycle provides the ROO-H bond enthalpies: DH(298)(CH(3)OO-H) = 87.8 +/- 1.0 kcal mol(-1), DH(298)(CD(3)OO-H) = 87.9 +/- 1.1 kcal mol(-1), and DH(298)(CH(3)CH(2)OO-H) = 84.8 +/- 2.2 kcal mol(-1). We review the thermochemistry of the peroxyl radicals, CH(3)OO and CH(3)CH(2)OO. Using experimental bond enthalpies, DH(298)(ROO-H), and CBS/APNO ab initio electronic structure calculations for the energies of the corresponding hydroperoxides, we derive the heats of formation of the peroxyl radicals. The "electron affinity/acidity/CBS" cycle yields Delta(f)H(298)[CH(3)OO] = 4.8 +/- 1.2 kcal mol(-1) and Delta(f)H(298)[CH(3)CH(2)OO] = -6.8 +/- 2.3 kcal mol(-1).     

Thermochemistry of the fluoroformyloxyl radical: a computational study based on coupled cluster theory

The standard enthalpy of formation of FCO(2) (X (2)B(2)) was determined by a computational approach based on coupled cluster theory [CCSD(T)] with energies extrapolated to the basis-set limit, with additional corrections accounting for core-valence correlation, scalar relativity, spin-orbit coupling, and zero-point vibrational motions. Utilizing a variety of independent reaction schemes, our best estimate is Delta(f)H(o)(0)(FCO(2)) = -86.0 +/- 0.6 kcal mol(-1) [Delta(f)H(o)(298) )(FCO(2)) = -86.7 +/- 0.6 kcal mol(-1)], which is shown to be more accurate than previous theoretical and experimental values. The chosen computational procedure was also applied to HCO (X (2)A'), where we find excellent agreement with experiment, and to FCO (X (2)A'), where we recommend an improved value of Delta(f)H(o)(0)(FCO) = -42.1 +/- 0.5 kcal mol(-1) [ Delta(f)H(o)(298)(FCO) = -42.0 +/- 0.5 kcal mol(-1)]. Further theoretical results concern the C-F bond dissociation energy, electron affinity, ionization energy, first and second excitation energies in FCO(2), fluoride ion affinity of CO(2), and equilibrium geometries of the molecules treated presently. For FCO (X (2)A') we propose an improved equilibrium structure: r(e)(CF) = 132.5(2) pm, r(e)(CO) = 116.7(2) pm, and theta(e)(FCO) = 127.8(2)(o).     

Thinking out of the black box: accurate barrier heights of 1,3-dipolar cycloadditions of ozone with acetylene and ethylene

Accurate barriers for the 1,3-dipolar cycloadditions of ozone with acetylene and ethylene have been determined via the systematic extrapolation of ab initio energies within the focal point approach of Allen and co-workers. Electron correlation has been accounted for primarily via coupled cluster theory, including single, double, and triple excitations, as well as a perturbative treatment of connected quadruple excitations [CCSD, CCSD(T), CCSDT, and CCSDT(Q)]. For the concerted [4 + 2] cycloadditions, the final recommended barriers are DeltaH(0K) = 9.4 +/- 0.2 and 5.3 +/- 0.2 kcal mol(-1) for ozone adding to acetylene and ethylene, respectively. These agree with recent results of Cremer et al. and Anglada et al., respectively. The reaction energy for O3 + C2H2 exhibits a protracted convergence with respect to inclusion of electron correlation, with the CCSDT/cc-pVDZ and CCSDT(Q)/cc-pVDZ values differing by 2.3 kcal mol-1. Recommended enthalpies of formation (298 K) for cycloadducts 1,2,3-trioxole and 1,2,3-trioxolane are +32.8 and -1.6 kcal mol(-1), respectively. Popular composite ab initio approaches [CBS-QB3, CBS-APNO, G3, G3B3, G3(MP2)B3, G4, G4(MP3), and G4(MP2)] predict a range of barrier heights for these systems. The CBS-QB3 computed barrier for ozone and acetylene, DeltaH(0K) = 4.4 kcal mol(-1), deviates by 5 kcal mol(-1) from the focal point value. CBS-QB3 similarly underestimates the barrier for the reaction of ozone and ethylene, yielding a prediction of only 0.7 kcal mol(-1). The errors in the CBS-QB3 results are significantly larger than mean errors observed in application to the G2 test set. The problem is traced to the nontransferability of MP2 basis set effects in the case of these reaction barriers. The recently published G4 and G4(MP2) approaches perform substantially better for O3 + C2H2, predicting enthalpy barriers of 9.0 and 8.4 kcal mol(-1), respectively. For the prediction of these reaction barriers, the additive corrections applied in the majority of the composite approaches considered lead to worse agreement with the reference focal point values than would be obtained relying only on single point energies evaluated at the highest level of theory utilized within each composite method.     

2-Pyridinethiol/2-pyridinethione tautomeric equilibrium. A comparative experimental and computational study

The gas phase and solvent dependent preference of the tautomerization between 2-pyridinethiol (2SH) and 2-pyridinethione (2S) has been assessed using variable temperature Fourier transform infrared (FTIR) experiments, as well as ab initio and density functional theory computations. No spectroscopic evidence (nu(S)(-)(H) stretch) for 2SH was observed in toluene, C(6)D(6), heptane, or methylene chloride solutions. Although, C(s)() 2SH is 2.61 kcal/mol more stable than C(s)() 2S (CCSD(T)/cc-pVTZ//B3LYP/6-311+G(3df,2p)+ZPE), cyclohexane solvent-field relative energies (IPCM-MP2/6-311+G(3df,2p)) favor 2S by 1.96 kcal/mol. This is in accord with the FTIR observations and in quantitative agreement with the -2.6 kcal/mol solution (toluene or C(6)D(6)) calorimetric enthalpy for the 2S/2SH tautomerization favoring the thione. As the intramolecular transition state for the 2S, 2SH tautomerization (2TS) lies 25 (CBS-Q) to 30 kcal/mol (CCSD/cc-pVTZ) higher in energy than either tautomer, tautomerization probably occurs in the hydrogen bonded dimer. The B3LYP/6-311+G(3df,2p) optimized C(2) 2SH dimer is 10.23 kcal/mol + ZPE higher in energy than the C(2)(h)() 2S dimer and is only 2.95 kcal/mol + ZPE lower in energy than the C(2) 2TS dimer transition state. Dimerization equilibrium measurements (FTIR, C(6)D(6)) over the temperature range 22-63 degrees C agree: K(eq)(298) = 165 +/- 40 M(-)(1), DeltaH = -7.0 +/- 0.7 kcal/mol, and DeltaS = -13.4 +/- 3.0 cal/(mol deg). The difference between experimental and B3LYP/6-311+G(3df,2p) [-34.62 cal/(mol deg)] entropy changes is due to solvent effects. The B3LYP/6-311+G(3df,2p) nucleus independent chemical shifts (NICS) are -8.8 and -3.5 ppm 1 A above the 2SH and 2S ring centers, respectively, and the thiol is aromatic. Although the thione is not aromatic, it is stabilized by the thioamide resonance. In solvent, the large 2S dipole, 2-3 times greater than 2SH, favors the thione tautomer and, in conclusion, 2S is thermodynamically more stable than 2SH in solution.     

Thermochemical studies of N-methylpyrazole and N-methylimidazole

The 351.1 nm photoelectron spectra of the N-methyl-5-pyrazolide anion and the N-methyl-5-imidazolide anion are reported. The photoelectron spectra of both isomers display extended vibrational progressions in the X2A' ground states of the corresponding radicals that are well reproduced by Franck-Condon simulations, based on the results of B3LYP/6-311++G(d,p) calculations. The electron affinities of the N-methyl-5-pyrazolyl radical and the N-methyl-5-imidazolyl radical are 2.054 +/- 0.006 eV and 1.987 +/- 0.008 eV, respectively. Broad vibronic features of the A(2)A' ' states are also observed in the spectra. The gas-phase acidities of N-methylpyrazole and N-methylimidazole are determined from measurements of proton-transfer rate constants using a flowing afterglow-selected ion flow tube instrument. The acidity of N-methylpyrazole is measured to be Delta(acid)G(298) = 376.9 +/- 0.7 kcal mol(-1) and Delta(acid)H(298) = 384.0 +/- 0.7 kcal mol(-1), whereas the acidity of N-methylimidazole is determined to be Delta(acid)G(298) = 380.2 +/- 1.0 kcal mol(-1) and Delta(acid)H(298)= 388.1 +/- 1.0 kcal mol(-1). The gas-phase acidities are combined with the electron affinities in a negative ion thermochemical cycle to determine the C5-H bond dissociation energies, D(0)(C5-H, N-methylpyrazole) = 116.4 +/- 0.7 kcal mol(-1) and D(0)(C5-H, N-methylimidazole) = 119.0 +/- 1.0 kcal mol(-1). The bond strengths reported here are consistent with previously reported bond strengths of pyrazole and imidazole; however, the error bars are significantly reduced.     

The influence of substituent groups on the resonance stabilization of benzene. An ab initio computational study

Accurate G3(MP2) calculations of the enthalpies of formation (Delta(f)H298) of organic molecules permit replication and extension of calculations that were formerly dependent on experimental thermochemical results. A case in point is Kistiakowski's classical calculation of the total stabilization enthalpy of benzene relative to that of cyclohexene, called for many years the "resonance energy". This paper investigates extension of the classical calculation to substituted benzenes. Slight modification of the usual procedure for Delta(f)H298 determination permits exclusion of all empirical information, leaving a purely ab initio result. Stabilization enthalpies relative to the corresponding 4-substituted cyclohexenes are presented for benzene, toluene, aniline, phenol, phenylacetylene, styrene, ethylbenzene, and phenylhydrazine. In the process of calculating these stabilization enthalpies, we have also obtained 42 values of Delta(f)H298 for monosubstituted benzenes, cyclohexenes, and cyclohexanes, 24 of which are not in the standard reference literature. For the remaining 18 G3(MP2) results, the unsigned mean difference between calculated Delta(f)H298 values and experimental results is +/-0.91 kcal x mol(-1).     

The ozonolysis of acetylene--a quantum chemical investigation.

The ozonolysis of acetylene was investigated using CCSD(T), CASPT2, and B3LYP-DFT in connection with a 6-311+G(2d,2p) basis set. The reaction is initiated by the formation of a van der Waals complex followed by a [4pi + 2pi] cycloaddition between ozone and acetylene (activation enthalpy DeltaH(a)(298) = 9.6 kcal/mol; experiment, 10.2 kcal/mol), yielding 1,2,3-trioxolene, which rapidly opens to alpha-ketocarbonyl oxide 5. Alternatively, an O atom can be transferred from ozone to acetylene (DeltaH(a)(298) = 15.6 kcal/mol), thus leading to formyl carbene, which can rearrange to oxirene or ketene. The key compound in the ozonolysis of acetylene is 5 because it is the starting point for the isomerization to the corresponding dioxirane 19 (DeltaH(a)(298) = 16.9 kcal/mol), for the cyclization to trioxabicyclo[2.1.0]pentane 10 (DeltaH(a)(298) = 19.5 kcal/mol), for the formation of hydroperoxy ketene 15 (DeltaH(a)(298) = 20.6 kcal/mol), and for the rearrangement to dioxetanone 9 (DeltaH(a)(298) = 23.6 kcal/mol). Compounds 19, 10, 15, and 9 rearrange or decompose with barriers between 13 and 16 kcal/mol to yield as major products formanhydride, glyoxal, formaldehyde, formic acid, and (to a minor extent) glyoxylic acid. Hence, the ozonolysis of acetylene possesses a very complicated reaction mechanism that deserves intensive experimental studies.     

Thermochemistry of methyl and ethyl nitro, RNO2, and nitrite, RONO, organic compounds

Computational quantum theory is employed to determine the thermochemical properties of n-alkyl nitro and nitrite compounds: methyl and ethyl nitrites, CH3ONO and C2H5ONO, plus nitromethane and nitroethane, CH3NO2 and C2H5NO2, at 298.15 K using multilevel G3, CBS-QB3, and CBS-APNO composite methods employing both atomization and isodesmic reaction analysis. Structures and enthalpies of the corresponding aci-tautomers are also determined. The enthalpies of formation for the most stable conformers of methyl and ethyl nitrites at 298 K are determined to be -15.64 +/- 0.10 kcal mol-1 (-65.44 +/- 0.42 kJ mol-1) and -23.58 +/- 0.12 kcal mol-1 (-98.32 +/- 0.58 kJ mol-1), respectively. DeltafHo(298 K) of nitroalkanes are correspondingly evaluated at -17.67 +/- 0.27 kcal mol-1 (-74.1 +/- 1.12 kJ mol-1) and -25.06 +/- 0.07 kcal mol-1 (-121.2 +/- 0.29 kJ mol-1) for CH3NO2 and C2H5NO2. Enthalpies of formation for the aci-tautomers are calculated as -3.45 +/- 0.44 kcal mol-1 (-14.43 +/- 0.11 kJ mol-1) for aci-nitromethane and -14.25 +/- 0.44 kcal mol-1 (-59.95 +/- 1.84 kJ mol-1) for the aci-nitroethane isomers, respectively. Data are evaluated against experimental and computational values in the literature with recommendations. A set of thermal correction parameters to atomic (H, C, N, O) enthalpies at 0 K is developed, to enable a direct calculation of species enthalpy of formation at 298.15 K, using atomization reaction and computation outputs.     

Heats of formation and singlet-triplet separations of hydroxymethylene and 1-hydroxyethylidene

Thermochemical parameters of hydroxymethylene (HC:OH) and 1-hydroxyethylidene (CH3C:OH) were evaluated by using coupled-cluster, CCSD(T), theory, in conjunction with the augmented correlation consistent, aug-cc-pVnZ, basis sets, with n = D, T, Q, and 5, extrapolated to the complete basis set limit. The predicted value at 298 K for Delta Hf(CH2O) is -26.0 +/- 1 kcal/mol, as compared to an experimental value of -25.98 +/- 0.01 kcal/mol, and for Delta Hf(CH:OH) it is 26.1 +/- 1 kcal/mol. The hydroxymethylene-formaldehyde energy gap is 52.1 +/- 0.5 kcal/mol, the singlet-triplet separation of hydroxymethylene is Delta E(ST)(HC:OH) = 25.3 +/- 0.5 kcal/mol, the proton affinity is PA(HC:OH) = 222.5 +/- 0.5 kcal/mol, and the ionization energy is IEa(HC:OH) = 8.91 +/- 0.04 eV. The predicted value at 298 K for Delta Hf(CH3CHO) is -39.1 +/- 1 kcal/mol as compared to an experimental value of -40.80 +/- 0.35 kcal/mol, and for Delta Hf(CH3C:OH) it is 11.2 +/- 1 kcal/mol. The hydroxyethylidene-acetaldehyde energy gap is 50.6 +/- 0.5 kcal/mol, the singlet-triplet separation of 1-hydroxyethylidene is Delta E(ST)(CH3C:OH) = 30.5 +/- 0.5 kcal/mol, the proton affinity is PA(CH3C:OH) = 234.7 +/- 0.5 kcal/mol, and the ionization energy is IEa(CH3C:OH) = 8.18 +/- 0.04 eV. The calculated energy differences between the carbene and aldehyde isomers, and, thus, the heats of formation of the carbenes, differ from the experimental values by 2.5 kcal/mol.     

Structures, internal rotor potentials, and thermochemical properties for a series of nitrocarbonyls, nitroolefins, corresponding nitrites, and their carbon centered radicals

Structures, enthalpy (Δ(f)H°(298)), entropy (S°(T)), and heat capacity (C(p)(T)) are determined for a series of nitrocarbonyls, nitroolefins, corresponding nitrites, and their carbon centered radicals using the density functional B3LYP and composite CBS-QB3 calculations. Enthalpies of formation (Δ(f)H°(298)) are determined at the B3LYP/6-31G(d,p), B3LYP/6-31+G(2d,2p), and composite CBS-QB3 levels using several work reactions for each species. Entropy (S) and heat capacity (C(p)(T)) values from vibration, translational, and external rotational contributions are calculated using the rigid-rotor-harmonic-oscillator approximation based on the vibration frequencies and structures obtained from the density functional studies. Contribution to Δ(f)H(T), S, and C(p)(T) from the analysis on the internal rotors is included. Recommended values for enthalpies of formation of the most stable conformers of nitroacetone cc(═o)cno2, acetonitrite cc(═o)ono, nitroacetate cc(═o)no2, and acetyl nitrite cc(═o)ono are -51.6 kcal mol(-1), -51.3 kcal mol(-1), -45.4 kcal mol(-1), and -58.2 kcal mol(-1), respectively. The calculated Δ(f)H°(298) for nitroethylene c═cno2 is 7.6 kcal mol(-1) and for vinyl nitrite c═cono is 7.2 kcal mol(-1). We also found an unusual phenomena: an intramolecular transfer reaction (isomerization) with a low barrier (3.6 kcal mol(-1)) in the acetyl nitrite. The NO of the nitrite (R-ONO) in CH(3)C(═O')ONO moves to the C═O' oxygen in a motion of a stretching frequency and then a shift to the carbonyl oxygen (marked as O' for illustration purposes).     

Accurate prediction of thermodynamic properties of alkyl peroxides by combining density functional theory calculation with least-square calibration

Owing to the significance in kinetic modeling of the oxidation and combustion mechanisms of hydrocarbons, a fast and relatively accurate method was developed for the prediction of Delta(f)H(298)(o) of alkyl peroxides. By this method, a raw Delta(f)H(298)(o) value was calculated from the optimized geometry and vibration frequencies at B3LYP/6-31G(d,p) level and then an accurate Delta(f)H(298)(o) value was obtained by a least-square procedure. The least-square procedure is a six-parameter linear equation and is validated by a leave-one out technique, giving a cross-validation squared correlation coefficient q(2) of 0.97 and a squared correlation coefficient of 0.98 for the final model. Calculated results demonstrated that the least-square calibration leads to a remarkable reduction of error and to the accurate Delta(f)H(298)(o) values within the chemical accuracy of 8 kJ mol(-1) except (CH(3))(2)CHCH(2)CH(2)CH(2)OOH which has an error of 8.69 kJ mol(-1). Comparison of the results by CBS-Q, CBS-QB3, G2, and G3 revealed that B3LYP/6-31G(d,p) in combination with a least-square calibration is reliable in the accurate prediction of the standard enthalpies of formation for alkyl peroxides. Standard entropies at 298 K and heat capacities in the temperature range of 300-1500 K for alkyl peroxides were also calculated using the rigid rotor-harmonic oscillator approximation.     

Thermochemistry, bond energies, and internal rotor potentials of dimethyl tetraoxide

Thermochemical properties of dimethyl tetraoxide (CH(3)OOOOCH(3)), the dimer of the methylperoxy radical, are studied using ab initio and density functional theory methods. Methylperoxy radicals are known to be important intermediates in the tropospheric ozone cycle, and the self-reaction of methylperoxy radicals, which is thought to proceed via dimethyl tetraoxide, leads to significant chain radical termination in this process. Dimethyl tetraoxide has five internal rotors, three of them unique; the potential energy profiles are calculated for these rotors, as well as for those in the CH(3)OO, CH(3)OOO, and CH(3)OOOO radicals. The dimethyl tetraoxide internal rotor profiles show barriers to rotation of 2-8 kcal mol(-1). Using B3LYP/6-31(d) geometries, frequencies, internal rotor potentials, and moments of inertia, we determine entropy and heat capacity values for dimethyl tetraoxide and its radicals. Isodesmic work reactions with the G3B3 and CBS-APNO methods are used; we calculate this enthalpy as -9.8 kcal mol(-1). Bond dissociation energies (BDEs) are calculated for all C-O and O-O bonds in dimethyl tetraoxide, again with the G3B3 and CBS-APNO theoretical methods, and we suggest the following BDEs: 46.0 kcal mol(-1) for CH(3)-OOOOCH(3), 20.0 kcal mol(-1) for CH(3)O-OOOCH(3), and 13.9 kcal mol(-1) for CH(3)OO-OOCH(3). From the BDE calculations and the isodesmic enthalpy of formation for dimethyl tetraoxide, we suggest enthalpies of 2.1, 5.8, and 1.4 kcal mol(-1) for the CH(3)OO, CH(3)OOO, and CH(3)OOOO radicals, respectively. We evaluate the suitability of 10 different density functional theory (DFT) methods for calculating thermochemical properties of dimethyl tetraoxide and its radicals with the 6-31G(d) and 6-311++G(3df,3pd) basis sets, using a variety of work reaction schemes. Overall, the best-performed DFT methods of those tested were TPSSh, BMK, and B1B95. Significant improvements in accuracy were made by moving from atomization to isodesmic work reactions, with most DFT methods yielding errors of less than 2 kcal mol(-1) with the 6-311++G(3df,3pd) basis set for isodesmic calculations on the dimethyl tetraoxide enthalpy. These isodesmic calculations were basis set consistent, with a considerable reduction in error found by using the 6-311++G(3df,3pd) basis set over the 6-31G(d) basis set. This was not the case, however, for atomization and bond dissociation work reactions, where the two basis sets returned similar results. Improved group additivity terms for the O-O-O moiety (O/O2 central atom group) are also determined.     

Transition state structure and energetics of the N(2)O + X (X = Cl,Br) reactions

The structural and vibrational properties of the transition state of the N(2)O + X (X = Cl,Br) reactions have been characterized by ab initio methods using density functional theory. We have employed Becke's hybrid functional (B3LYP), and transition state optimizations were performed with 6-31G(d), 6-311G(2d,2p), 6-311+G(3d,2p), and 6-311+G(3df,2p) basis sets. For the chlorine atom reaction the coupled-cluster method (CCSD(T)) with 6-31G(d) basis set was also used. All calculations resulted in transition state structures with a planar cis arrangement of atoms for both reactions. The geometrical parameters of transition states at B3LYP are very similar, and the reaction coordinates involve mainly the breaking of the N-O bond. At CCSD(T)/6-31G(d) level a contribution of the O-Cl forming bond is also observed in the reaction coordinate. In addition, several highly accurate ab initio composite methods of Gaussian-n (G1, G2, G3), their variations (G2(MP2), G3//B3LYP), and complete basis set (CBS-Q, CBS-Q//B3LYP) series of models were applied to compute reaction energetics. All model chemistries predict exothermic reactions. The G3 and G2 methods result in the smallest deviations from experiment, 1.8 and 0 kcal mol(-1), for the enthalpies of reaction for N(2)O reaction with chlorine and bromine, respectively. The G3//B3LYP and G1 methods perform best among the composite methods in predicting energies of the transition state, with a deviation of 1.9 and 3.0 kcal mol(-1), respectively, in the activation energies for the above processes. However, the B3LYP/6-311+G(3df,2p) method gives smaller deviations of 0.4 and -1.0 kcal mol(-1), respectively. The performance of the methodologies applied in predicting transition state energies was analyzed.     

Thermochemical properties, rotation barriers, and group additivity for unsaturated oxygenated hydrocarbons and radicals resulting from reaction of vinyl and phenyl radical systems with O2

Oxidation of unsaturated and aromatic hydrocarbons in atmospheric and combustion processes results in formation of linear and cyclic unsaturated, oxygenated-hydrocarbon intermediates. The thermochemical parameters delatafH degrees 298, S degrees 298, and C(p)(f298)(T) for these intermediates are needed to understand their stability and reaction paths in further oxidation. These properties are not available for a majority of these unsaturated oxy-hydrocarbons and their corresponding radicals, even via group additivity methods. Enthalpy, entropy, and heat capacity of a series of 40 oxygenated and non-oxygenated molecules, or radicals corresponding to hydrogen atom loss from the parent stable molecules are determined in this study. Enthalpy (delatafH degrees 298 in kcal mol(-1)) is derived from the density function calculations at the B3LYP/6-311g(d,p) calculated enthalpy of reaction (delatafH degrees rxn,298) and by use of isodesmic (work) reactions. Estimation of error in enthalpy delatafH degrees 298, from use of computational chemistry coupled with work reactions analysis, is presented using comparisons between the calculated and literature enthalpies of reaction. Entropies (S degrees 298) and heat capacities (C(p)(f298)(T)) were calculated using the B3LYP/6-311G(d,p) determined frequencies and geometries. Potential barriers for internal rotors in each molecule were determined and used (in place of torsion frequencies) to calculate contributions to S and C(p)(T) from the hindered rotors. Twenty-six groups for use in group additivity (GA) are also developed.     

Protonated acetylene revisited

The potential energy surface for protonated acetylene has been re-examined with large basis sets and highly correlated methods. The energy difference of 3.6–3.8 kcal/mol between the classical structure and non-classical (bridged) structure computed with CCSD (T)/cc-pVQZ, CCSD(T)/6-311+G(3df,2pd), BD(T)/cc- pVQZ, BD(T)/6-311+G(3df,2pd) and CBS-APNO methods is in very good agreement with the best previous calculations, 3.7–4.0 kcal/mol. In contrast, BLYP, B3LYP, PW91, PBE and TPSS density functional methods do rather poorly, yielding −0.52. 0.29, 1.81, 2.16 and 0.62 kcal/mol, respectively, with the 6-311+G(3df,2pd) basis. MP2 calculations predict the classical structure to be a transition state; however, frequency calculations at the CCSD/6-311+G(3df,2pd) level of theory show that the classical structure is a local minimum. CCSD(T), BD(T) and CBS-APNO energy calculations along the MP2/6-311+G(3df,2pd) reaction path indicate that the classical structure is a shallow local minimum separated from the non-classical structure by a very small barrier of 0.11–0.13 kcal/mol. Because the barrier for proton exchange between the non-classical isomers via the classical structure is broad and nearly flat at the top, the tunneling splitting should be reduced, possibly accounting for the 15% difference between the calculated and experimental barrier heights. Contribution to the Fernando Bernardi Memorial Issue.     

Cheletropic decomposition of cyclic nitrosoamines revisited: the nature of the transition states and a critical role of the ring strain

The cheletropic decompositions of 1-nitrosoaziridine (1), 1-nitroso-Delta(3)-pyrroline (2), 7-nitroso-7-azabicyclo[2.2. 1]hepta-2,5-diene (3), and 6-nitroso-6-azabicyclo[2.1.1]hexa-4-ene (4) have been studied theoretically using high level ab initio computations. Activation parameters of the decomposition of nitrosoaziridine 1 were obtained experimentally in heptane (DeltaH()(298) = 18.6 kcal mol(-)(1), DeltaS()(298) = -7.6 cal mol(-)(1) K(-)(1)) and methanol (20.3 kcal mol(-)(1), 0.3 cal mol(-)(1) K(-)(1)). Among employed theoretical methods (B3LYP, MP2, CCD, CCSD(T)//CCD), the B3LYP method in conjunction with 6-31+G, 6-311+G, and 6-311++G(3df,2pd) basis sets gives the best agreement with experimental data. It was found that typical N-nitrosoheterocycles 2-4 which have high N-N bond rotation barriers (>16 kcal mol(-)(1)) extrude nitrous oxide via a highly asynchronous transition state with a planar ring nitrogen atom. Nitrosoaziridine 1, with a low rotation barrier (<9 kcal mol(-)(1)) represents a special case. This compound can eliminate N(2)O via a low energy linear synperiplanar transition state (DeltaH()(298) = 20.6 kcal mol(-)(1), DeltaS()(298) = 2.5 cal mol(-)(1) K(-)(1)). Two higher energy transition states are also available. The B3LYP activation barriers of the cheletropic fragmentation of nitrosoheterocycles 2-4 decrease in the series: 2 (58 kcal mol(-)(1)) > 3 (18 kcal mol(-)(1)) > 4 (12) kcal mol(-)(1). The relative strain energies increase in the same order: 2 (0 kcal mol(-)(1)) < 3 (39 kcal mol(-)(1)) < 4 (52 kcal mol(-)(1)). Comparison of the relative energies of 2-4 and their transition states on a common scale where the energy of nitrosopyrroline 2 is assumed as reference indicates that the thermal stability of the cyclic nitrosoamines toward cheletropic decomposition is almost entirely determined by the ring strain.     

Substituent effects on singlet-triplet gaps and mechanisms of 1,2-rearrangements of 1,3-oxazol-2-ylidenes to 1,3-oxazoles

Electronic structures, partial atomic charges, singlet-triplet gaps (Delta E ST), substituent effects, and mechanisms of 1,2-rearrangements of 1,3-oxazol-2-ylidene ( 5) and 4,5-dimethyl- ( 6), 4,5-difluoro- ( 7), 4,5-dichloro- ( 8), 4,5-dibromo- ( 9), and 3-methyl-1,3-oxazol-2-ylidene ( 10) to the corresponding 1,3-oxazoles have been studied using complete-basis-set methods (CBS-QB3, CBS-Q, CBS-4M), second-order M?ller-Plesset perturbation method (MP2), hybrid density functionals (B3LYP, B3PW91), coupled-cluster theory with single and double excitations (CCSD) and CCSD plus perturbative triple excitations [CCSD(T)], and the quadratic configuration interaction method including single and double excitations (QCISD) and QCISD plus perturbative triple excitations [QCISD(T)]. The 6-311G(d,p), 6-31+G(d,p), 6-311+G(d,p), and correlation-consistent polarized valence double-xi (cc-pVDZ) basis sets were employed. The carbenes have singlet ground states, and the CBS-QB3 and CBS-Q methods predict Delta E ST values for 5- 8 and 10 of 79.9, 79.8, 74.7, 77.0, and 82.0 kcal/mol, respectively. CCSD(T), QCISD(T), B3LYP, and B3PW91 predict smaller Delta E ST values than CBS-QB3 and CBS-Q, with the hybrid density functionals predicting the smallest values. The concerted unimolecular exothermic out-of-plane 1,2-rearrangements of singlet 1,3-oxazol-2-ylidenes to their respective 1,3-oxazoles proceed via cyclic three-center transition states. The CBS-predicted barriers to the 1,2-rearrangements of singlet carbenes 5- 9 to their respective 1,3-oxazoles are 41.4, 40.4, 37.8, 40.4, and 40.5 kcal/mol, respectively. During the 1,2-rearrangements of singlet 1,3-oxazol-2-ylidenes 5- 9, there is a decrease in electron density at oxygen, N3 (the migration origin), and C5 and an increase in electron density at C2 (the migration terminus), C4, and the partially positive migrating hydrogen.     

Energetics of hydroxybenzoic acids and of the corresponding carboxyphenoxyl radicals. Intramolecular hydrogen bonding in 2-hydroxybenzoic acid

The energetics of the phenolic O-H bond in the three hydroxybenzoic acid isomers and of the intramolecular hydrogen O-H- - -O-C bond in 2-hydroxybenzoic acid, 2-OHBA, were investigated by using a combination of experimental and theoretical methods. The standard molar enthalpies of formation of monoclinic 3- and 4-hydroxybenzoic acids, at 298.15 K, were determined as Delta(f)(3-OHBA, cr) = -593.9 +/- 2.0 kJ x mol(-1) and Delta(f)(4-OHBA, cr) = -597.2 +/- 1.4 kJ x mol(-1), by combustion calorimetry. Calvet drop-sublimation calorimetric measurements on monoclinic samples of 2-, 3-, and 4-OHBA, led to the following enthalpy of sublimation values at 298.15 K: Delta(sub)(2-OHBA) = 94.4 +/- 0.4 kJ x mol(-1), Delta(sub)(3-OHBA) = 118.3 +/- 1.1 kJ x mol(-1), and Delta(sub)(4-OHBA) = 117.0 +/- 0.5 kJ x mol(-1). From the obtained Delta(f)(cr) and Delta(sub) values and the previously reported enthalpy of formation of monoclinic 2-OHBA (-591.7 +/- 1.3 kJ x mol(-1)), it was possible to derive Delta(f)(2-OHBA, g) = -497.3 +/- 1.4 kJ x mol(-1), Delta(f)(3-OHBA, g) = -475.6 +/- 2.3 kJ x mol(-1), and Delta(f)(4-OHBA, cr) = -480.2 +/- 1.5 kJ x mol(-1). These values, together with the enthalpies of isodesmic and isogyric gas-phase reactions predicted by density functional theory (B3PW91/aug-cc-pVDZ, MPW1PW91/aug-cc-pVDZ, and MPW1PW91/aug-cc-pVTZ) and the CBS-QMPW1 methods, were used to derive the enthalpies of formation of the gaseous 2-, 3-, and 4-carboxyphenoxyl radicals as (2-HOOCC(6)H(4)O(*), g) = -322.5 +/- 3.0 kJ.mol(-1) Delta(f)(3-HOOCC(6)H(4)O(*), g) = -310.0 +/- 3.0 kJ x mol(-1), and Delta(f)(4-HOOCC(6)H(4)O(*), g) = -318.2 +/- 3.0 kJ x mol(-1). The O-H bond dissociation enthalpies in 2-OHBA, 3-OHBA, and 4-OHBA were 392.8 +/- 3.3, 383.6 +/- 3.8, and 380.0 +/- 3.4 kJ x mol(-1), respectively. Finally, by using the ortho-para method, it was found that the H- - -O intramolecular hydrogen bond in the 2-carboxyphenoxyl radical is 25.7 kJ x mol(-1), which is ca. 6-9 kJ x mol(-1) above the one estimated in its parent (2-OHBA), viz. 20.2 kJ x mol(-1) (theoretical) or 17.1 +/- 2.1 kJ x mol(-1) (experimental).     

Peculiar structure of the HOOO(-) anion

The HOOO(-) anion (1) can adopt a triplet state (T-1) or a singlet state (S-1), where the former is 9.8 kcal/mol (DeltaH(298) = 10.3 kcal/mol) more stable than the latter. S-1 possesses a strong O-OOH bond with some double bond character and a weakly covalent OO-OH bond (1.80 A) according to CCSD(T)/6-311++G(3df,3pd) calculations (the longest O-O bond ever found for a peroxide). In aqueous solution, S-1 adopts a geometry closely related to that of HOOOH (OO(O), 1.388 A; (O)OO(H), 1.509 A; tau(OOOH), 78.3 degrees ), justifying that S-1 is considered the anion of HOOOH. Dissociation into HO anion and O(2)((1)Delta(g)) requires 15.4 (DeltaH(298) = 14.3; DeltaG(298) = 8.9) kcal/mol. Structure T-1 corresponds to a van der Waals complex between HO anion and O(2)((3)Sigma(g)(-)) having a binding energy of 2.7 (DeltaH(298) = 2.1) kcal/mol. Modes of generating S-1 in aqueous solution are discussed, and it is shown that S-1 represents an important intermediate in ozonation reactions.