Activation energies of pericyclic reactions: performance of DFT, MP2, and CBS-QB3 methods for the prediction of activation barriers and reaction energetics of 1,3-dipolar cycloadditions, and revised activation enthalpies for a standard set of hydrocarbon pericyclic reactions

Activation barriers and reaction energetics for the three main classes of 1,3-dipolar cycloadditions, including nine different reactions, were evaluated with the MPW1K and B3LYP density functional methods, MP2, and the multicomponent CBS-QB3 method. The CBS-QB3 values were used as standards for 1,3-dipolar cycloaddition activation barriers and reaction energetics, and the density functional theory (DFT) and MP2 methods were benchmarked against these values. The MPW1K/6-31G* method and basis set performs best for activation barriers, with a mean absolute deviation (MAD) value of 1.1 kcal/mol. The B3LYP/6-31G* method and basis set performs best for reaction enthalpies, with a MAD value of 2.4 kcal/mol, while the MPW1K method shows large errors for reaction energetics. The MP2 method gives the expected systematic underestimation of barriers. Concerted and nearly synchronous transition structures are predicted by all DFT and MP2 methods. Also reported are revised estimated 0 K experimental activation enthalpies for a standard set of hydrocarbon pericyclic reactions and updated comparisons to experiment for DFT, ab initio, and multicomponent methods. B3LYP and MPW1K methods with MAD values of 1.5 and 2.1 kcal/mol, respectively, fortuitously outperform the multicomponent CBS-QB3 method, which has a MAD value of 2.3. The MAD value of the O3LYP functional improves to 2.4 kcal/mol from the previously reported 3.0 kcal/mol.     

Theory of 1,3-dipolar cycloadditions: distortion/interaction and frontier molecular orbital models

Quantum chemical calculations of activation barriers and reaction energies for 1,3-dipolar cycloadditions by the high-accuracy CBS-QB3 method reveal previously unrecognized quantitative trends in activation barriers. The distortion/interaction model of reactivity explains why (1) there is a monotonic decrease of approximately 6 kcal/mol in the activation energy along the series oxides, imine, and ylide for the diazonium, nitrilium, and azomethine betaine classes of 1,3-dipoles; (2) nitrilium and azomethine betaines with the same trio of atoms have almost identical cycloaddition barrier heights; (3) barrier heights for the cycloadditions of a given 1,3-dipole with ethylene and acetylene have the same activation energies (mean absolute deviation of 0.6 kcal/mol) in spite of very different reaction thermodynamics (Delta DeltaH(rxn) range = 14-43 kcal/mol) and frontier molecular orbital (FMO) energy gaps. The energy to distort the 1,3-dipole and dipolarophile to the transition state geometry, rather than FMO interactions or reaction thermodynamics, controls reactivity for cycloadditions of 1,3-dipoles with alkenes or alkynes. A distortion/interaction energy analysis was also carried out on the transition states for the cycloadditions of diazonium dipoles with a set of substituted alkenes (CH2CHX, X = OMe, Me, CO 2Me, Cl, CN) and reveals that FMO interaction energies between the 1,3-dipole and the dipolarophile differentiate reactivity when transition state distortion energies are nearly constant.     

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.     

Thermolysis of alkyl sulfoxides and derivatives: a comparison of experiment and theory.

Gas-phase activation data were obtained for model sulfoxide elimination reactions. The activation enthalpy for methyl 3-phenylpropyl sulfoxide is 32.9 +/- 0.9 kcal/mol. Elimination by methyl vinyl sulfoxide to form acetylene has an enthalpic barrier of 41.6 +/- 0.8 kcal/mol and that of 3-phenylpropyl methanesulfinate to form hydrocinnamaldehyde is 34.6 +/- 0.6 kcal/mol. Calculations at the MP2/6-311+G(3df,2p)//MP2/6-31G(d,p) level for simplified models of these reactions provide barriers of 32.3, 40.3, and 32.7 kcal/mol, respectively. A series of other compounds are examined computationally, and it is shown that the substituent effects on the sulfoxide elimination reaction are much more straightforward to interpret if DeltaH data are available in addition to the usually determined DeltaH++. The activation enthalpy of the reverse addition reaction is also subject to structural variation and can usually be rationalized on the basis of nucleophilicity of the sulfur or polarity matching between the sulfenic acid and olefin derivative.     

Computational methods to calculate accurate activation and reaction energies of 1,3-dipolar cycloadditions of 24 1,3-dipoles

Theoretical calculations were performed on the 1,3-dipolar cycloaddition reactions of 24 1,3-dipoles with ethylene and acetylene. The 24 1,3-dipoles are of the formula X≡Y(+)-Z(-) (where X is HC or N, Y is N, and Z is CH(2), NH, or O) or X═Y(+)-Z(-) (where X and Z are CH(2), NH, or O and Y is NH, O, or S). The high-accuracy G3B3 method was employed as the reference. CBS-QB3, CCSD(T)//B3LYP, SCS-MP2//B3LYP, B3LYP, M06-2X, and B97-D methods were benchmarked to assess their accuracies and to determine an accurate method that is practical for large systems. Several basis sets were also evaluated. Compared to the G3B3 method, CBS-QB3 and CCSD(T)/maug-cc-pV(T+d)Z//B3LYP methods give similar results for both activation and reaction enthalpies (mean average deviation, MAD, < 1.5 kcal/mol). SCS-MP2//B3LYP and M06-2X give small errors for the activation enthalpies (MAD < 1.5 kcal/mol), while B3LYP has MAD = 2.3 kcal/mol. SCS-MP2//B3LYP and B3LYP give the reasonable reaction enthalpies (MAD < 5.0 kcal/mol). The B3LYP functional also gives good results for most 1,3-dipoles (MAD = 1.9 kcal/mol for 17 common 1,3-dipoles), but the activation and reaction enthalpies for ozone and sulfur dioxide are difficult to calculate by any of the density functional methods.     

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.     

CCSDT and CCSD(T) calculations on model H-bonded and stacked complexes

The CCSD(T) and CCSDT interaction energies were determined for model planar H-bonded complexes (formamide…formamide, formamidine…formamidine) and stacked complexes (ethylene…ethylene, formaldehyde…formaldehyde). Various basis sets from the 6-31G*(0.25) to aug-cc-pVDZ were used. Difference between CCSD(T) and CCSDT interaction energies were small and become negligible (bellow 0.1 kcal/mol) if the aug-cc-pVDZ (or aug-cc-pVDZ/cc-pVDZ) basis set was applied. This result strongly supports the use of the CCSD(T) method for determination of true stabilization energies of extended complexes.     

Transition states for the dimerization of 1,3-cyclohexadiene: a DFT, CASPT2, and CBS-QB3 quantum mechanical investigation

Quantum mechanical calculations using restricted and unrestricted B3LYP density functional theory, CASPT2, and CBS-QB3 methods for the dimerization of 1,3-cyclohexadiene (1) reveal several highly competitive concerted and stepwise reaction pathways leading to [4 + 2] and [2 + 2] cycloadducts, as well as a novel [6 + 4] ene product. The transition state for endo-[4 + 2] cycloaddition (endo-2TS, DeltaH(double dagger)(B3LYP(0K)) = 28.7 kcal/mol and DeltaH(double dagger)(CBS-QB3(0K)) = 19.0 kcal/mol) is not bis-pericyclic, leading to nondegenerate primary and secondary orbital interactions. However, the C(s) symmetric second-order saddle point on the B3LYP energy surface is only 0.3 kcal/mol above endo-2TS. The activation enthalpy for the concerted exo-[4 + 2] cycloaddition (exo-2TS, DeltaH(double dagger)(B3LYP(0K)) = 30.1 kcal/mol and DeltaH(double dagger)(CBS-QB3(0K)) = 21.1 kcal/mol) is 1.4 kcal/mol higher than that of the endo transition state. Stepwise pathways involving diallyl radicals are formed via two different C-C forming transition states (rac-5TS and meso-5TS) and are predicted to be competitive with the concerted cycloaddition. Transition states were located for cyclization from intermediate rac-5 leading to the endo-[4 + 2] (endo-2) and exo-[2 + 2] (anti-3) cycloadducts. Only the endo-[2 + 2] (syn-3) transition state was located for cyclization of intermediate meso-5. The novel [6 + 4] "concerted" ene transition state (threo-4TS, DeltaH(double dagger)(UB3LYP(0K)) = 28.3 kcal/mol) is found to be unstable with respect to an unrestricted calculation. This diradicaloid transition state closely resembles the cyclohexadiallyl radical rather than the linked cyclohexadienyl radical. Several [3,3] sigmatropic rearrangement transition states were also located and have activation enthalpies between 27 and 31 kcal/mol.     

Predictions of substituent effects in thermal azide 1,3-dipolar cycloadditions: implications for dynamic combinatorial (reversible) and click (irreversible) chemistry

Substituent effects in 1,3-dipolar cycloadditions of azides with alkenes and alkynes were investigated with the high-accuracy CBS-QB3 method. The possibilities for noncatalytic activation and the reversibility or irreversibility of these reactions was explored; the possibilities for uses in dynamic combinatorial chemistry (DCC) or click chemistry were explored. The activation enthalpies for reactions of ethylene and acetylene with hydrazoic acid, formyl, phenyl-, methyl-, and methanesulfonylazides exhibit modest variation, with Delta H++ ranging from 17 to 20 kcal/mol. A detailed study of formylazide cycloadditions with various alkenes and alkynes reveals a narrow range of activation enthalpies (17-21 kcal/mol). The activation enthalpies for the reactions of azides with alkenes and alkynes are similar. FMO theory and distortion/interaction energy control have been used to rationalize the rates and regiochemistries of cycloadditions involving alkene dipolarophiles. Significantly, triazoles, formed from alkynes, are 30-40 kcal/mol more stable than tetrazolines formed from alkenes. On the basis of initial reactant concentrations, kinetic and thermodynamic values are suggested for the identification of reversible reactions that approach equilibrium over 24 h, as well as for fast irreversible reactions. Although azide cycloadditions are suitable for irreversible chemistry and are typically unsuitable for reversible applications, theoretical procedures established by these studies have provided guidelines for the prediction of useful reversible libraries.     

Thermochemical properties, DeltafH degrees (298), S degrees (298), and Cp degrees (T), for n-butyl and n-pentyl hydroperoxides and the alkyl and peroxy radicals, transition states, and kinetics for intramolecular hydrogen shift reactions of the peroxy radicals

Alkyl radicals in atmospheric and combustion environments undergo a rapid association with molecular oxygen (3O2) to form an alkyl peroxy radical (ROO*). One important reaction of these peroxy radicals is the intramolecular H-shift (intramolecular abstraction) to form a hydroperoxide alkyl radical (R'*COOH), where the hydroperoxide alkyl radical may undergo chemical activation reaction with O2 and result in chain branching at moderate to low temperatures. The thermochemistry and trends in kinetic parameters for the hydrogen shift reactions from each carbon (4-8-member-ring TST's) in n-butyl and n-pentyl peroxy radicals (CCCCOO* and CCCCCOO*) are analyzed using density functional and ab initio calculation methods. Thermochemical properties, DeltafH degrees (298 K), C-H bond energies, S degrees (298 K), and Cp degrees (T) of saturated linear C4 and C5 aliphatic peroxides (ROOH), as well as the corresponding hydroperoxide alkyl radicals (R'*COOH), are determined. DeltafH degrees (298 K) are obtained from isodesmic reactions and the total energies of the CBS-QB3 and B3LYP computational methods. Contributions to the entropy and the heat capacity from translation, vibration, and external rotation are calculated using the rigid-rotor-harmonic-oscillator approximation based on the CBS-QB3 frequencies and structures. The results indicate that pre-exponential factors, A(T), decrease with the increase of the ring size (4-8-member-ring TS, H-atom included). The DeltaH for 4-, 5-, 6-, and 7-member rings in n-butyl (and n-pentyl) peroxy are 40.8 (40.8), 31.4 (31.5), 20.5 (20.0), 22.6-p (19.4) kcal mol(-1), respectively. The DeltaH for the 8-member ring in n-pentylperoxy is 23.8-p kcal mol(-1), All abstractions are from secondary (-CH2-) groups except those marked (-p), which are from primary sites. Enthalpy and barrier values from the B3LYP/6-311++G(2d,p) and BHandHLYP/6-311G(d,p) methods are compared with CBS-QB3 results. The B3LYP results show good agreement with the higher level CBS-QB3 calculation method; the BHandH barriers for the intramolecular peroxy H-shifts are not acceptable.     

Reactivity of alkyl versus silyl peroxides. The consequences of 1, 2-silicon bridging on the epoxidation of alkenes with silyl hydroperoxides and bis(trialkylsilyl)peroxides

The bond dissociation energies for a series of silyl peroxides have been calculated at the G2 and CBS-Q levels of theory. A comparison is made with the O-O BDE of the corresponding dialkyl peroxides, and the effect of the O-O bond strength on the activation barrier for oxygen atom transfer is discussed. The O-O bond dissociation enthalpies (DeltaH(298)) for bis (trimethylsilyl) peroxide (1) and trimethylsilyl hydroperoxide (2) are 54.8 and 53.1 kcal/mol, respectively at the G2 (MP2) and CBS-Q levels of theory. The O-O bond dissociation energies computed at G2 and G2(MP2) levels for bis(tert-butyl) peroxide and tert-butyl hydroperoxide are 45.2 and 48.3 kcal/mol, respectively. The barrier height for 1,2-methyl migration from silicon to oxygen in trimethylsilyl hydroperoxide is 47.9 kcal/mol (MP4//MP2/6-31G). The activation energy for the oxidation of trimethylamine to its N-oxide by bis(trimethylsilyl) peroxide is 28.2 kcal/mol (B3LYP/6-311+G(3df,2p)// B3LYP/6-31G(d)). 1,2-Silicon bridging in the transition state for oxygen atom transfer to a nucleophilic amine results in a significant reduction in the barrier height. The barrier for the epoxidation of E-2-butene with bis(dimethyl(trifluoromethyl))silyl peroxide is 25.8 kcal/mol; a reduction of 7.5 kcal/mol relative to epoxidation with 1. The activation energy calculated for the epoxidation of E-2-butene with F(3)SiOOSiF(3) is reduced to only 2.2 kcal/mol reflecting the inductive effect of the electronegative fluorine atoms.     

Toward subchemical accuracy in computational thermochemistry: focal point analysis of the heat of formation of NCO and [H,N,C,O] isomers

In continuing pursuit of thermochemical accuracy to the level of 0.1 kcal mol(-1), the heats of formation of NCO, HNCO, HOCN, HCNO, and HONC have been rigorously determined using state-of-the-art ab initio electronic structure theory, including conventional coupled cluster methods [coupled cluster singles and doubles (CCSD), CCSD with perturbative triples (CCSD(T)), and full coupled cluster through triple excitations (CCSDT)] with large basis sets, conjoined in cases with explicitly correlated MP2-R12/A computations. Limits of valence and all-electron correlation energies were extrapolated via focal point analysis using correlation consistent basis sets of the form cc-pVXZ (X=2-6) and cc-pCVXZ (X=2-5), respectively. In order to reach subchemical accuracy targets, core correlation, spin-orbit coupling, special relativity, the diagonal Born-Oppenheimer correction, and anharmonicity in zero-point vibrational energies were accounted for. Various coupled cluster schemes for partially including connected quadruple excitations were also explored, although none of these approaches gave reliable improvements over CCSDT theory. Based on numerous, independent thermochemical paths, each designed to balance residual ab initio errors, our final proposals are DeltaH(f,0) ( composite function )(NCO)=+30.5, DeltaH(f,0) ( composite function )(HNCO)=-27.6, DeltaH(f,0) ( composite function )(HOCN)=-3.1, DeltaH(f,0) ( composite function )(HCNO)=+40.9, and DeltaH(f,0) ( composite function )(HONC)=+56.3 kcal mol(-1). The internal consistency and convergence behavior of the data suggests accuracies of +/-0.2 kcal mol(-1) in these predictions, except perhaps in the HCNO case. However, the possibility of somewhat larger systematic errors cannot be excluded, and the need for CCSDTQ [full coupled cluster through quadruple excitations] computations to eliminate remaining uncertainties is apparent.     

Experimental determination of the heat of hydrogenation of phenylcyclobutadiene

The heat of hydrogenation of phenylcyclobutadiene (DeltaH degrees (hyd) = 57.4 +/- 4.9 kcal mol(-1)) was determined via a thermodynamic cycle by carrying out gas-phase measurements on 1-phenylcyclobuten-3-yl cation. This leads to an antiaromatic destabilization energy of 27 +/- 5 kcal mol(-1), a difference of 9.6 +/- 4.9 kcal mol(-1) for the first and second C-H bond dissociation energies of 1-phenylcyclobutene, and an estimate of 96 +/- 5 kcal mol(-1) for the heat of formation of cyclobutadiene. These results are compared to G3, G3(MP2), and B3LYP computations and represent the first experimental measurements of the energy of a monocyclic cyclobutadiene.     

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.     

10]Annulene: bond shifting and conformational mechanisms for automerization

We report density-functional and coupled-cluster calculations on conformation change and degenerate bond shifting in [10]annulene isomers 1-5. At the CCSD(T)/cc-pVDZ//CCSD/6-31G level, conversion of the twist (1) to the heart (2) has a barrier of 10.1 kcal/mol, compared to Ea = 16.2 kcal/mol for degenerate "two-twist" bond shifting in 1. Pseudorotation in the all-cis boat isomer (3) proceeds with a negligible barrier. The naphthalene-like isomer 4 has a 3.9 kcal/mol barrier to degenerate bond shifting. The azulene-like isomer 5 is the only species for which the nature of the bond-equalized form (5-eq) depends on the method. At the CCSD(T)/cc-pVDZ//CCSD/6-31G level, 5-eq is 1.2 kcal/mol more stable than the bond-alternating form 5-alt. Conversion of 5-eq to 4 has a barrier of 12.6 kcal/mol. Despite being significantly nonplanar, both 5-eq and the transition state for bond shifting in 4 are highly aromatic based on magnetic susceptibility exaltations. On the basis of a detailed consideration of these mechanisms and barriers, we can now, with greater confidence, rule out 4 and 5 as candidates to explain the NMR spectra observed by Masamune. Our results support Masamune's original assignments for both isolated isomers.     

Quantum mechanical study of the potential energy surface of the ClO + NO2 reaction

In the study of the reaction pathways of the ClO + NO2 reaction including reliable structures of the reactants, products, intermediates, and transition states as well as energies the MP2/6-311G(d), B3LYP/6-311G(d), and G2(MP2) methods have been employed. Chlorine nitrate, ClONO2, is formed by N-O association without an entrance barrier and is stabilized by 29.8 kcal mol(-1). It can undergo either a direct 1,3 migration of Cl or OCl rotation to yield an indistinguishable isomer. The corresponding barriers are 45.8 and 7.1 kcal mol(-1), respectively. ClONO2 can further decompose into NO3 + Cl with an endothermicity of 46.4 kcal mol(-1). The overall endothermicity of the NO2 + ClO --> NO3 + Cl reaction is calculated to be 16.6 kcal mol(-1). The formation of cis-perp and trans-perp conformer of chlorine preoxynitrite, ClOONO(cp) and ClOONO(tp), are exothermic by 5.4 and 3.8 kcal mol(-1), respectively. Calculations on the possible reaction pathways for the isomerization of ClOONO to ClONO2 showed that the activation barriers are too high to account for appreciable nitrate formation from peroxynitrite isomerization. All quoted relative energies are related to G2(MP2) calculations.     

Theoretical determinations of the ambient conformational distribution and unimolecular decomposition of n-propylperoxy radical

The conformational distribution and unimolecular decomposition pathways for the n-propylperoxy radical have been generated at the CBS-QB3, B3LYP/6-31+G and mPW1K/6-31+G levels of theory. At each of the theoretical levels, the 298 K Boltzmann distributions and rotational profiles indicate that all five unique rotamers of the n-propylperoxy radical can be expected to be present in significant concentrations at thermal equilibrium. At the CBS-QB3 level, the 298 K distribution of rotamers is predicted to be 28.1, 26.4, 19.6, 14.0, and 11.9% for the gG, tG, gT, gG', and tT conformations, respectively. The CBS-QB3 C-OO bond dissociation energy (DeltaH298 K) for the n-propylperoxy radical has been calculated to be 36.1 kcal/mol. The detailed CBS-QB3 potential energy surface for the unimolecular decomposition of the n-propylperoxy radical indicates that important bimolecular products could be derived from two 1,4-H transfer mechanisms available at T < 500 K, primarily via an activated n-propylperoxy adduct.     

Conformational energies for 2-substituted butanes

The conformational free energies for some 2-substituted butanes where X = F, Cl, CN, and CCH were calculated using G3-B3, CBS-QB3, and CCSD(T)/6-311++G(2d,p) as well as other theoretical levels. The above methods gave consistent results with free energies relative to the trans conformers as follows: X = CCH, g+ = 0.77 +/- 0.05 kcal/mol. g- = 0.88 +/- 0.05 kcal/mol; X = CN, g+ = 0.85 +/- 0.05 kcal/mol, g- = 0.75 +/- 0.05 kcal/mol; X = Cl, g+ = 0.70 +/- 0.05 kcal/ml, g- = 0.80 +/- 0.05 kcal/mol; and X = F, g+ = 0.53 +/- 0.05 kcal/mol, g- = 0.83 +/- 0.05 kcal/mol. The conformational free energies also were estimated using the observed liquid phase IR spectra and intensities calculated using B3LYP/6-311++G** and MP2/6-311++G**. The rotational free energy profiles for all of the compounds were estimated at the G3-B3 level.     

An ab-initio study of some homolytic substitution reactions of acyl radicals at silicon, germanium and tint

Ab initio calculations using the 6-311G**, cc-pVDZ, and (valence) double-zeta pseudopotential (DZP) basis sets, with (MP2, QCISD, CCSD(T)) and without (UHF) the inclusion of electron correlation, and density functional (BHandHLYP, B3LYP) calculations predict that homolytic substitution reactions of acetyl radicals at the silicon atoms in disilane can proceed via both backside and frontside attack mechanisms. At the highest level of theory (CCSD(T)/cc-pVDZ//MP2/cc-pVDZ), energy barriers (deltaE double dagger) of 77.2 and 81.9 kJ mol(-1) are calculated for the backside and frontside reactions respectively. Similar results are obtained for reactions involving germanium and tin with energy barriers (deltaE double dagger) of 53.7-84.2, and 55.0-89.7 kJ mol(-1) for the backside and frontside mechanisms, respectively. These data suggest that both homolytic substitution mechanisms are feasible for homolytic substitution reactions of acetyl radicals at silicon, germanium, and tin. BHandHLYP calculations provide geometries and energy barriers for backside and frontside transition states in good agreement with those obtained by traditional ab initio techniques.