The effect of substitutents on the strain energies of small ring compounds

The effect of substitutents on the strain energy (SE) of cyclic molecules is examined at the CBS, G2, and G2(MP2) levels of theory. Alkyl substitutents have a meaningful effect upon the SE of small ring compounds. gem-Dimethyl substitution lowers the strain energy of cyclopropanes, cyclobutanes, epoxides, and dimethyldioxirane (DMDO) by 6-10 kcal/mol relative to an unbranched acyclic reference molecule. The choice of the reference compound is especially important for geminal electronegative substitutents. The SE of 1,1-difluorocyclopropane is estimated to be 20.5 kcal/mol relative to acyclic reference molecule 1,3-difluoropropane but is 40.7 kcal/mol with respect to the thermodynamically more stable (DeltaE = -20.2 kcal/mol) isomeric reference compound 2,2-difluoropropane. The SE of dioxirane (DO) is estimated to be approximately 18 kcal/mol while the SE of DMDO is predicted to be approximately equal to 11 kcal/mol by using homodesmotic reactions that maintain a balanced group equivalency. The total energy (CBS-APNO) of DMDO is 2.6 kcal/mol lower than that of isomeric 1,2-dioxacyclopentane that has an estimated SE of 5 kcal/mol. The thermodynamic stability of DMDO is a consequence of its relatively strong C-H (BDE = 102.7 kcal/mol) and C-CH(3) (BDE = 98.9 kcal/mol) bonds. By comparison, the calculated sec-C-H and -C-CH(3) G2 bond dissociation energies in propane are 100.3 and 90.5 kcal/mol.     

Relative reactivity of peracids versus dioxiranes (DMDO and TFDO) in the epoxidation of alkenes. A combined experimental and theoretical analysis

Comparative analysis of the calculated gas-phase activation barriers (DeltaE++) for the epoxidation of ethylene with dimethyldioxirane (DMDO) and peroxyformic acid (PFA) [15.2 and 16.4 kcal/mol at QCISD(T)// QCISD/6-31+G(d,p)] and E-2-butene [14.3 and 13.2 kcal/mol at QCISD(T)/6-31G(d)//B3LYP/6-311+G(3df,2p)] suggests similar oxygen atom donor capacities for both oxidants. Competition experiments in CH(2)Cl(2) solvent reveal that DMDO reacts with cyclohexene much faster than peracetic acid/acetic acid under scrupulously dried conditions. The rate of DMDO epoxidation is catalyzed by acetic acid with a reduction in the classical activation barrier of 8 kcal/mol. In many cases, the observed increase in the rate for DMDO epoxidation in solution may be attributed to well-established solvent and hydrogen-bonding effects. This predicted epoxidative reactivity for DMDO is not consistent with what has generally been presumed for a highly strained cyclic peroxide. The strain energy (SE) of DMDO has been reassessed and its moderated value (about 11 kcal/mol) is now more consistent with its inherent gas-phase reactivity toward alkenes in the epoxidation reaction. The unusual thermodynamic stability of DMDO is largely a consequence of the combined geminal dimethyl- and dioxa-substitution effects and unusually strong C-H and C-CH(3) bonds. Methyl(trifluoromethyl)dioxirane (TFDO) exhibits much lower calculated activation barriers than DMDO in the epoxidation reaction (the average DeltaDeltaE++ values are about 7.5 kcal/mol). The rate increase relative to DMDO of approximately 10(5), while consistent with the higher strain energy for TFDO (SE approximately 19 kcal/mol) is attributed largely to the inductive effect of the CF(3) group. We have also examined the effect of alkene strain on the rate of epoxidation with PFA. The epoxidation barriers are only slightly higher for the strained alkenes cyclopropene (DeltaE++ = 14.5 kcal/mol) and cyclobutene (DeltaE++ = 13.7 kcal/mol) than for cyclopentene (DeltaE++ = 12.1 kcal/mol), reflecting the fact there is little relief of strain in the transition state. Alkenes strained by twist or pi-bond torsion do exhibit much lower activation barriers.     

The effect of carbonyl substitution on the strain energy of small ring compounds and their six-member ring reference compounds

High level ab initio calculations have been applied to the estimation of ring strain energies (SE) of a series of three- and six-member ring compounds. The SE of cyclohexane has been estimated to be 2.2 kcal/mol at the CBS-APNO level of theory. The SE of cyclopropane has been increased to 28.6 kcal/mol after correction for the one-half of the SE of cyclohexane. The SEs of a series of carbonyl-containing three-member ring compounds have been estimated at the CBS-Q level by their combination with cyclopropane to produce a six-member ring reference compound. The SEs of cyclopropanone (5), the simplest alpha-lactone (6) [oxiranone], and alpha-lactam (7) [aziridinone] have been predicted to be 49, 47, and 55 kcal/mol, respectively, after correction for the SE of the corresponding six-member ring reference compound. The SEs of cyclohexanone, delta-valerolactone, and delta-valerolactam have been estimated to be 4.3, 11.3, and 5.1 kcal/mol, respectively. Marked increases in the SE of silacyclopropane and siladioxirane have been established, while significant decreases in the SEs of phosphorus, sulfur, dioxa- and diaza-containing three-member ring compounds were observed. The ring strain energies of the hydrocarbons (but not heterocycles) exhibit a strong correlation with their C-H bond dissociation energies.     

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.     

Strain energy of small ring hydrocarbons. Influence of C-h bond dissociation energies

Ab initio calculations at the G2, G3, and CBS-Q levels of theory have been applied to the question of the origin of ring strain in a series of unsaturated hydrocarbons. In addition to the angular ring strain germane to all three-membered ring hydrocarbons, a general trend is in evidence that suggests that the increased ring strain (SE) of unsaturated small ring alkenes may be attributed in part to their relatively weak allylic C-H bonds. The high strain energy of cyclopropene (54.1 kcal/ mol) is attributed largely to angular strain. The anomalously low SE of cyclobutene relative to cyclobutane (DeltaSE = 4 kcal/mol) is a consequence of normal C-H bond dissociation energies for cyclobutane (100.6 kcal/mol) and very strong vinyl C-H bonds (111.9 kcal/mol) and a relatively strong pi-bond energy (63.5 kcal/mol) for cyclobutene. The greater SE of methylenecyclopropane (39.5 kcal/ mol), relative to methylcyclopropane (29.8 kcal/mol), can be attributed to the strong ring C-H bonds of methylcyclopropane (110.5 kcal/mol) and relatively weak allylic C-H bonds (99.3 kcal/mol) of methylenecyclopropane. The increased SE of 1-methylcyclopropene relative to isomeric methylenecyclopropane is ascribed to its weak ring C-H bonds and to angular strain. The relative thermodynamic stability of a series of small ring alkenes is determined by a measure of their hydrogenation enthalpies. Independent confirmation of the SEs of a series of substituted cyclopropenes is provided by their dimerization/combination with cyclopropane to form a six-membered ring reference compound.     

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.     

Cage strain in nitro-substituted 1,3,5,7-tetraazacubanes

Homodesmotic reactions were designed for the computation of strain energies (SE) for four nitro-substituted 1,3,5,7-tetraazacubane derivatives. Total energies of the optimized geometric structures at the DFT-B3LYP/6-31G* and DFT-B3LYP/6-311G** levels were used to derive the SE. The variation of SE with respect to the number of substituents is similar with both basis sets. The SE value is 237.32 kcal/mol at the B3LYP/6-311G** level for 2,4,6,8-tetranitro-1,3,5,7-tetraazacubane, which is unexpectedly much larger than that of its cubane analogue. The SE increases remarkably with more nitro groups being attached to the cage skeleton of tetraazacubane. The ‘bending’ of the bonds within the cubic skeleton attributes to the increase of strains as the attached number of nitro groups increases.     

Formaldehyde oxime <--> nitrosomethane tautomerism

Formaldehyde oxime <--> nitrosomethane tautomerism, isomeric nitrone, and their common cations and anions are studied with Gaussian-2 theory using MP2(full)/6-31G geometries and with density functional theory using B3LYP/6-311+G**. Geometrical parameters, harmonic vibrational frequencies, relative stabilities, conformational stabilities, and ionization energies are compared with experimental gas-phase data when available. The formaldehyde oxime <--> nitrosomethane tautomerism is compared with the amide <--> imidol, imine <--> enamine, keto <--> enol, and nitro <--> aci-nitro tautomeric processes. Solvent effects are estimated by the self-consistent isodensity polarizable continuum model (SCIPCM). The influence of hydrogen bonding interactions with the solvent is addressed by including two water molecules. In the final evaluation, formaldehyde oxime is 15.8 kcal/mol more stable than nitrosomethane when the aqueous solvation correction of 3.8 kcal/mol is applied to the G2 energies. Unsolvated formaldehyde oxime is estimated to be 11.1 kcal/mol more stable than nitrone. The estimated gas-phase ionization energies (G2) are 362.5 kcal/mol for formaldehyde oxime, 350.6 kcal/mol for nitrosomethane, and 351.4 kcal/mol for nitrone.     

Improved modification for the density-functional theory calculation of thermodynamic properties for C-H-O composite compounds

A three-parametric modification equation and the least-squares approach are adopted to calibrating hybrid density-functional theory energies of C(1)-C(10) straight-chain aldehydes, alcohols, and alkoxides to accurate enthalpies of formation DeltaH(f) and Gibbs free energies of formation DeltaG(f), respectively. All calculated energies of the C-H-O composite compounds were obtained based on B3LYP6-311++G(3df,2pd) single-point energies and the related thermal corrections of B3LYP6-31G(d,p) optimized geometries. This investigation revealed that all compounds had 0.05% average absolute relative error (ARE) for the atomization energies, with mean value of absolute error (MAE) of just 2.1 kJ/mol (0.5 kcal/mol) for the DeltaH(f) and 2.4 kJ/mol (0.6 kcal/mol) for the DeltaG(f) of formation.     

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.     

Probing the bent bonds in cyclopropane systems for gas storage and separation process: A computational study

The hydrogen, carbon dioxide, and carbon monoxide gas adsorption and storage capacity of lithium-decorated cyclopropane ring systems were examined with quantum chemical calculations at density functional theory, DFT M06-2X functional using 6-31G(d) and cc-pVDZ basis sets. To examine the reliability of M06-2X DFT functional, a few representative systems are also examined with complete basis set CBS-QB3 method and CCSD-aug-cc-pVTZ level of theory. The cyclopropane systems can bind to one Li⁺ ion; however, the corresponding the methylated systems can bind with two Li⁺ ions. The cyclopropane systems can adsorb six hydrogen molecules with an average binding energy of 3.8 kcal/mol. The binding free energy (ΔG) values suggest that the hydrogen adsorption process is feasible at 273.15 K. The calculation of desorption energies indicates the recyclable property of gas adsorbed complexes. The same number of CO₂ and CO gas molecules can also be adsorbed with an average binding energy of −14.4 kcal/mol and −10.7 kcal/mol, respectively. The carbon dioxide showed ~3–4 kcal/mol better binding energy as compared to carbon monoxide and hence such designed systems can function as a potential candidate for the separation of these flue gas molecules. The nature of interactions in complexes was examined with atoms in molecules analysis revealed the electrostatic nature for the interaction of Li⁺ ion with cyclopropane rings. The chemical hardness and electrophilicity calculations showed that the gas adsorbed complexes are rigid and therefore robust as gas storage materials.     

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.     

Oxyfunctionalization of non-natural targets by dioxiranes. 5. Selective oxidation of hydrocarbons bearing cyclopropyl moieties

The powerful methyl(trifluoromethyl)dioxirane (1b) was employed to achieve the direct oxyfunctionalization of 2,4-didehydroadamantane (5), spiro[cyclopropane-1,2'-adamantane] (9), spiro[2.5]octane (17), and bicyclo[6.1.0]nonane (19). The results are compared with those attained in the analogous oxidation of two alkylcyclopropanes, i.e., n-butylcyclopropane (11) and (3-methyl-butyl)-cyclopropane (14). The product distributions observed for 11 and 14 show that cyclopropyl activation of alpha-C-H bonds largely prevails when no tertiary C-H are present in the open chain in the tether; however, in the oxyfunctionalixation of 14 cyclopropyl activation competes only mildly with hydroxylation at the tertiary C-H. The application of dioxirane 1b to polycyclic alkanes possessing a sufficiently rigid framework (such as 5 and 9) demonstrates the relevance of relative orientation of the cyclopropane moiety with respect to the proximal C-H undergoing oxidation. At one extreme, as observed in the oxidation of rigid spiro compound 9, even bridgehead tertiary C-H's become deactivated by the proximal cyclopropyl moiety laying in the unfavorable "eclipsed" (perpendicular) orientation; at the other end, a cyclopropane moiety constrained in a favorable "bisected" orientation (as for didehydroadamantane 5) can activate an "alpha" methylene CH2 to compete effectively with dioxirane O-insertion into tertiary C-H bonds. Comparison with literature reports describing similar oxidations by dimethyldioxirane (1a) demonstrate that methyl(trifluoromethyl)dioxirane (1b) presents similar selectivity and remarkably superior reactivity.     

Toward true DNA base-stacking energies: MP2, CCSD(T), and complete basis set calculations

Stacking energies in low-energy geometries of pyrimidine, uracil, cytosine, and guanine homodimers were determined by the MP2 and CCSD(T) calculations utilizing a wide range of split-valence, correlation-consistent, and bond-functions basis sets. Complete basis set MP2 (CBS MP2) stacking energies extrapolated using aug-cc-pVXZ (X = D, T, and for pyrimidine dimer Q) basis sets equal to -5.3, -12.3, and -11.2 kcal/mol for the first three dimers, respectively. Higher-order correlation corrections estimated as the difference between MP2 and CCSD(T) stacking energies amount to 2.0, 0.7, and 0.9 kcal/mol and lead to final estimates of the genuine stacking energies for the three dimers of -3.4, -11.6, and -10.4 kcal/mol. The CBS MP2 stacking-energy estimate for guanine dimer (-14.8 kcal/mol) was based on the 6-31G(0.25) and aug-cc-pVDZ calculations. This simplified extrapolation can be routinely used with a meaningful accuracy around 1 kcal/mol for large aromatic stacking clusters. The final estimate of the guanine stacking energy after the CCSD(T) correction amounts to -12.9 kcal/mol. The MP2/6-31G(0.25) method previously used as the standard level to calculate aromatic stacking in hundreds of geometries of nucleobase dimers systematically underestimates the base stacking by ca. 1.0-2.5 kcal/mol per stacked dimer, covering 75-90% of the intermolecular correlation stabilization. We suggest that this correction is to be considered in calibration of force fields and other cheaper computational methods. The quality of the MP2/6-31G(0.25) predictions is nevertheless considerably better than suggested on the basis of monomer polarizability calculations. Fast and very accurate estimates of the MP2 aromatic stacking energies can be achieved using the RI-MP2 method. The CBS MP2 calculations and the CCSD(T) correction, when taken together, bring only marginal changes to the relative stability of H-bonded and stacked base pairs, with a slight shift of ca. 1 kcal/mol in favor of H-bonding. We suggest that the present values are very close to ultimate predictions of the strength of aromatic base stacking of DNA and RNA bases.     

Computational study on the mechanisms and energetics of trimethylindium reactions with H2O and H2S

The reactions of trimethylindium (TMIn) with H2O and H2S are relevant to the chemical vapor deposition of indium oxide and indium sulfide thin films. The mechanisms and energetics of these reactions in the gas phase have been investigated by density functional theory and ab initio calculations using the CCSD(T)/[6-31G(d,p)+Lanl2dz]//B3LYP/[6-31G(d,p)+Lanl2dz] and CCSD(T)/[6-31G(d,p)+Lanl2dz] //MP2/[6-31G(d,p)+Lanl2dz] methods. The results of both methods are in good agreement for the optimized geometries and relative energies. When TMIn reacts with H2O and H2S, initial molecular complexes [(CH3)3In:OH2 (R1)] and [(CH3)3In:SH2 (R2)] are formed with 12.6 and 3.9 kcal/mol binding energies. Elimination of a CH4 molecule from each complex occurs with a similar energy barrier at TS1 (19.9 kcal/mol) and at TS3 (22.1 kcal/mol), respectively, giving stable intermediates (CH3)2InOH and (CH3)2InSH. The elimination of the second CH4 molecule from these intermediate products, however, has to overcome very high and much different barriers of 66.1 and 53.2 kcal/mol, respectively. In the case of DMIn with H2O and H2S reactions, formation of both InO and InS is exothermic by 3.1 and 30.8 kcal/mol respectively. On the basis of the predicted heats of formation of R1 and R2 at 0 K and -20.1 and 43.6 kcal/mol, the heats of formation of (CH3)2InOH, (CH3)2InSH, CH3InO, CH3InS, InO, and InS are estimated to be -20.6, 31.8, and 29.0 and 48.4, 35.5, and 58.5 kcal/mol, respectively. The values for InO and InS are in good agreement with available experimental data. A similar study on the reactions of (CH3)2In with H2O and H2S has been carried out; in these reactions CH3InOH and CH3InSH were found to be the key intermediate products.     

Accuracy of free energies of hydration for organic molecules from 6-31g*-derived partial charges

Absolute free energies of hydration have been computed for 13 diverse organic molecules using partial charges derived from ab initio 6-31G* wave functions. Both Mulliken charges and charges fit to the electrostatic potential surface (EPS) were considered in conjunction with OPLS Lennard–Jones parameters for the organic molecules and the TIP4P model of water. Monte Carlo simulations with statistical perturbation theory yielded relative free energies of hydration. These were converted to absolute quantities through perturbations to reference molecules for which absolute free energies of hydration had been obtained previously in TIP4P water. The average errors in the computed absolute free energies of hydration are 1.1 kcal/mol for the 6-31G* EPS charges and 4.0 kcal/mol for the Mulliken charges. For the EPS charges, the largest individual errors are under 2 kcal/mol except for acetamide, in which case the error is 3.7 kcal/mol. The hydrogen bonding between the organic solutes and water has also been characterized. © John Wiley & Sons, Inc.     

Accurate interaction energies of hydrogen-bonded nucleic acid base pairs

Hydrogen-bonded nucleic acids base pairs substantially contribute to the structure and stability of nucleic acids. The study presents reference ab initio structures and interaction energies of selected base pairs with binding energies ranging from -5 to -47 kcal/mol. The molecular structures are obtained using the RI-MP2 (resolution of identity MP2) method with extended cc-pVTZ basis set of atomic orbitals. The RI-MP2 method provides results essentially identical with the standard MP2 method. The interaction energies are calculated using the Complete Basis Set (CBS) extrapolation at the RI-MP2 level. For some base pairs, Coupled-Cluster corrections with inclusion of noniterative triple contributions (CCSD(T)) are given. The calculations are compared with selected medium quality methods. The PW91 DFT functional with the 6-31G basis set matches well the RI-MP2/CBS absolute interaction energies and reproduces the relative values of base pairing energies with a maximum relative error of 2.6 kcal/mol when applied with Becke3LYP-optimized geometries. The Becke3LYP DFT functional underestimates the interaction energies by few kcal/mol with relative error of 2.2 kcal/mol. Very good performance of nonpolarizable Cornell et al. force field is confirmed and this indirectly supports the view that H-bonded base pairs are primarily stabilized by electrostatic interactions.     

To what extent can methyl derivatives be regarded as stabilized tautomers of xanthine?

Taking xanthine (Xan) as an example, validity of an approach to experimental investigations of nucleotide bases' tautomeric equilibrium, based on the use of methyl derivatives corresponding to their prototropic tautomers, was studied by (1)H NMR in dimethylsulfoxide (DMSO) and by quantum chemical calculations at the B3LYP/6-311++G(d,p) level of theory. From (1)H NMR spectra of m(7)Xan, m(9)Xan, Xan and m(3)Xan conclusion was made that the N7H tautomeric forms of the last two compounds dominate in solution, which was supported by quantum chemical data. Calculated relative energies of the N9H tautomers of Xan and m(3)Xan (8.74 and 9.57 kcal/mol, accordingly) are rather close to the m(9)Xan one (9.11 kcal/mol). Nonspecific influence of DMSO modelled by the COSMO algorithm therewith reduces these values by approximately 2-3 kcal/mol. The data obtained imply that methyl derivatives are rather good models of high-energy tautomers of nucleotide bases, if their relative energies are not less than a few kcal/mol.     

Mechanism of OH formation from ozonolysis of isoprene: a quantum-chemical study

The formation and unimolecular reactions of primary ozonides and carbonyl oxides arising from the O(3)-initiated reactions of isoprene have been investigated using density functional theory and ab initio molecular orbital calculations. The activation energies of O(3) cycloaddition to the two double bonds of isoprene are found to be comparable (3.3-3.4 kcal mol(-1)), implying that the initial two O(3) addition pathways are nearly equally accessible. The reaction energies of O(3) addition to isoprene are between -47 and -48 kcal mol(-1). Cleavage of primary ozonides to form carbonyl oxides occurs with a barrier of 11-16 kcal mol(-1) above the ground state of the primary ozonide, and the decomposition energies range from -5 to -13 kcal mol(-1). OH formation is shown to occur primarily via decomposition of the carbonyl oxides with the syn-positioned methyl (alkyl) group, which is more favorable than isomerization to form dioxirane (by 1.1-3.3 kcal mol(-1)). Using the transition-state theory and master equation formalism, we determine an OH yield of 0.25 from prompt and thermal decomposition of the carbonyl oxides.