Comparative study of surface cycloadditions of ethylene and 2-butene on the Si(100)-2 x 1 surface

Multireference wave functions were used to study the ethylene and 2-butene surface reactions on Si(100) in their lowest energy singlet states. In addition to the diradical pathway, a pi-complex pathway on the ethylene surface was found. The net barrier for the latter process is 4.5 kcal/mol higher than that for the former, making the pi-complex pathway kinetically less accessible. Therefore, although there is a competition between the two initial channels, the diradical path is slightly favored, and rotational isomerization is possible. However, since the initial potential energy surfaces of the two channels are different, depending on experimental conditions, the branching ratio between the two channels may change. Consequently, the combined effects that would favor one channel over the other may not derive directly from the initial reaction barrier. This provides an explanation of the experimental controversy. As a result, the final distributions of surface products may depend on the experimental kinetic environment, especially when the population change due to the rotational isomerization is expected to be very small. A significantly different reaction channel is found in the 2-butene surface reaction on Si(100), in which a methyl hydrogen easily transfers to the surface yielding a new type of surface product other than the expected [2 + 2] cycloaddition product, with a comparatively small activation barrier. Consequently, the overall surface reactions of ethylene and 2-butene may be quite different. Therefore, direct comparisons between ethylene and 2-butene experimental results would be very useful.     

Reactions of 1,3-cyclohexadiene with singlet oxygen. A theoretical study.

A thorough study of the reaction of singlet oxygen with 1,3-cyclohexadiene has been made at the B3LYP/6-31G(d) and CASPT2(12e,10o) levels. The initial addition reaction follows a stepwise diradical pathway to form cyclohexadiene endoperoxide with an activation barrier of 6.5 kcal/mol (standard level = CASPT2(12e,10o)/6-31G(d); geometries and zero-point corrections at B3LYP/6-31G(d)), which is consistent with an experimental value of 5.5 kcal/mol. However, as the enthalpy of the transition structure for the second step is lower than the diradical intermediate, the reaction might also be viewed as a nonsynchronous concerted reaction. In fact, the concertedness of the reaction is temperature dependent since entropy differences create a free energy barrier for the second step of 1.8 kcal/mol at 298 K. There are two ene reactions; one is a concerted mechanism (DeltaH(double dagger) = 8.8 kcal/mol) to 1-hydroperoxy-2,5-cyclohexadiene (5), while the other, which forms 1-hydroperoxy-2,4-cyclohexadiene (18), passes through the same diradical intermediate (9) as found on the pathway to endoperoxide. The major pathway from the endoperoxide is O-O bond cleavage (22.0 kcal/mol barrier) to form a 1,4-diradical (25), which is 13.9 kcal/mol less stable than the endoperoxide. From the diradical, two low-energy pathways exist, one to epoxyketone (29) and the other to the diepoxide (27), where both products are known to be formed experimentally with a product ratio sensitive to the nature of substitutents. A significantly higher activation barrier leads to C-C bond cleavage and direct formation of maleic aldehyde plus ethylene.     

New insight into the gas-phase bimolecular self-reaction of the HOO radical

The singlet and triplet potential energy surfaces (PESs) for the gas-phase bimolecular self-reaction of HOO*, a key reaction in atmospheric environments, have been investigated by means of quantum-mechanical electronic structure methods (CASSCF and CASPT2). All the reaction pathways on both PESs consist of a first step involving the barrierless formation of a prereactive doubly hydrogen-bonded complex, which is a diradical species lying about 8 kcal/mol below the energy of the reactants at 0 K. The lowest energy reaction pathway on both PESs is the degenerate double hydrogen exchange between the HOO* moieties of the prereactive complex via a double proton transfer mechanism involving an energy barrier of only 1.1 kcal/mol for the singlet and 3.3 kcal/mol for the triplet at 0 K. The single H-atom transfer between the two HOO* moieties of the prereactive complex (yielding HOOH + O2) through a pathway keeping a planar arrangement of the six atoms involves a conical intersection between either two singlet or two triplet states of A' and A" symmetries. Thus, the lowest energy reaction pathway occurs via a nonplanar cisoid transition structure with an energy barrier of 5.8 kcal/mol for the triplet and 17.5 kcal/mol for the singlet at 0 K. The simple addition between the terminal oxygen atoms of the two HOO* moieties of the prereactive complex, leading to the straight chain H2O4 intermediate on the singlet PES, involves an energy barrier of 7.3 kcal/mol at 0 K. Because the decomposition of such an intermediate into HOOH + O2 entails an energy barrier of 45.2 kcal/mol at 0 K, it is concluded that the single H-atom transfer on the triplet PES is the dominant pathway leading to HOOH + O2. Finally, the strong negative temperature dependence of the rate constant observed for this reaction is attributed to the reversible formation of the prereactive complex in the entrance channel rather than to a short-lived tetraoxide intermediate.     

Triplet ground state derivative of aza-m-xylylene diradical with large singlet-triplet energy gap

Organic molecules with a strong preference for triplet ground states, in which the triplet state is below the lowest singlet state by ≥10 kcal/mol, are typically short-lived and mostly detected as reactive intermediates. We now report a triplet ground state derivative of aza-m-xylylene diradical with a large singlet-triplet energy gap (ΔE(ST)) of ～10 kcal/mol, which is comparable to ΔE(ST) for the well-known reactive intermediate m-xylylene diradical. The aminyl diradical persists in solution at room temperature on the time scale of minutes.     

Ab initio calculations find 2,2-disilylcyclopentane-1,3-diyl is a singlet diradical with a high barrier to ring closure

Ab initio calculations on the lowest singlet and triplet states of 2,2-disilylcyclopentane-1,3-diyl find that the singlet lies well below the triplet. The C ₂ singlet diradical is calculated to be a minimum on the potential energy surface with an enthalpic barrier to ring closure of ΔH ^‡ ₂₉₈ = 13.5 kcal/mol at the CASPT2/6-31G* level of theory. The energy of the 1,3-divinyl-substituted singlet diradical is calculated to be only 0.8 kcal/mol higher than that of 5,5-disilyl-1,3-divinylbicyclo[2.1.0]pentane at this level of theory, but the transition state for their equilibration is computed to be 12.8 kcal/mol above the diradical in energy. Received: 2 July 1998 / Accepted: 4 August 1998 / Published online: 16 November 1998     

Theoretical investigation of the mechanisms of reaction of NCN with NO and NS

Quantum-chemical calculations were performed on the mechanisms of reaction of NCN with NO and NS. Possible mechanisms were classified according to four pathways yielding products in the following four possible groups: N2O/N2S + CN, N2 + NCO/NCS, N2 + CNO/CNS, and CNN + NO/NS, labeled in order from p1/p1s to p4/p4s. The local structures, transition structures, and potential-energy surfaces with respect to the reaction coordinates are calculated, and the barriers are compared. In the NCN + NO reaction, out of several adduct structures, only the nitroso adduct NCNNO lies lower in energy than the reactants, by 21.89 kcal/mol; that adduct undergoes rapid transformation into the products, in agreement with experimental observation. For the NS counterpart, both thionitroso NCNNS and thiazyl NCNSN adducts have energies much lower than those of the reactants, by 43 and 29 kcal/mol, respectively, and a five-membered-ring NCNNS (having an energy lower than those of the reactants by 36 kcal/mol) acts as a bridge in connecting these two adducts. The net energy barriers leading to product channels other than p4s are negative for the NS reaction, whereas those for the NO analogue are all positive. The channel leading to p1 (N2O + CN) has the lowest energy (3.81 kcal/mol), whereas the channels leading to p2 (N2 + NCO) and p2s (N2 + NCS) are the most exothermic (100.94 and 107.38 kcal/mol, respectively).     

Geometry optimization of the ring-opened oxirane diradical: mechanism of the addition reaction of the triplet oxygen atom to olefins

Geometry optimizations of several low-lying diradical states of the ring-opened oxirane (·CH₂CH₂O·) were performed by using the energy gradients of the UHF MINDO/3, STO-3G and 4-31G solutions. Both the STO-3G and 4-3 IG methods predict that the most stable form is the triplet state of the non-twisted σπ conformation in which the unpaired spins localized on the terminal carbon and oxygen atoms are oriented perpendicularly to each other. The singlet σσ diradical state in which both the radical-site p orbitals are coplanar with the molecular framework is only 2.3 (STO-3G) and 1.2 (4-31G) kcal/mol less stable than the triplet σπ diradical state. It is found that the geometry of the singlet σσ diradical is unique in that the C-C-O angle is noticeably small as compared with various other diradical states. Implications of these results to the mechanism of the oxirane-forming O(³P) + C₂H₄ reaction are discussed.     

The almost bottleable triplet carbene: 2,6-dibromo-4-tert-butyl-2',6'-bis(trifluoromethyl)-4'-isopropyldiphenylcarbene.

Computations on 2,6-dibromo-4-tert-butyl-2',6'-bis(trifluoromethyl)-4'-isopropyldiphenylcarbene (1) using ab initio and density functional theory methods underscore the unusual stability of the triplet over the singlet state. At the B3LYP/6-311G(d,p) level, the triplet state had a slightly bent central C-C-C bond angle of 167 degrees, whereas this angle in the singlet was 134 degrees. The B3LYP singlet-triplet splitting (12.2 kcal/mol) was larger than that of the parent molecule (5.8 kcal/mol), diphenylcarbene (2), which also has a triplet ground state. The energy of a suitable isodesmic reaction showed the triplet and singlet states of (1) to be destabilized, by 6.3 and 12.5 kcal/mol, respectively, due to the combined effects of the CF3, Br, and alkyl substituents. The linear-coplanar form of (3)(1), which might facilitate dimerization or electrophilic attack at the more exposed diradical center, was prohibitively (35.9 kcal/mol) higher in energy. Our results confirm Tomioka's conclusion that the triplet diarylcarbene, ortho-substituted with bulky CF3 and Br substituents, is persistent due to steric protection of the diradical center. Dimerization and other possible reaction pathways are inhibited, not only by the bulky ortho substituents but also by the para alkyl groups. The increase in stability of the triplet ((3)(1)) state relative to the singlet ((1)(1)) state does not influence the reactivity directly.     

2-Silyl group effect on the reactivity of cyclopentane-1,3-diyls. Intramolecular ring-closure versus silyl migration

Generation of singlet and triplet 2-silylcyclopentane-1,3-diyls and their reactivity have been investigated in the thermal and photochemical denitrogenation of 2,3-diaza-7-silylbicyclo[2.2.1]hept-2-ene. 5-Silylcyclopentene (silyl migration product) is quantitatively obtained, while 5-silylbicyclo[2.1.0]pentane (intramolecular ring-closure product) is not detected in the denitrogenation reactions. Deuterium labeling studies clarify that 5-silylcyclopentene is formed by a suprafacial [1,2] silyl migration in singlet 2-silylcyclopentane-1,3-diyl. UDFT calculations closely reproduce the observed reactivity of the singlet diradical: The enthalpic barriers of the intramolecular ring-closure are calculated to be DeltaH++exo468 = 5.8 kcal/mol and DeltaH++endo468 = 6.7 kcal/mol, which are much higher than the energy barrier for the [1,2] silyl migration, DeltaH++468 = 2.7 kcal/mol. The notable effect of the silyl group on raising the energy barrier of the intramolecular cyclization is rationalized by an electronic configuration of the lowest singlet state of 2-silylcyclopentane-1,3-diyls.     

Chemical behavior of the biradicaloid (HO...ONO) singlet states of peroxynitrous acid. The oxidation of hydrocarbons, sulfides, and selenides

Various high levels of theory have been applied to the characterization of two higher lying biradicaloid metastable singlet states of peroxynitrous acid. A singlet minimum (cis-2) was located that had an elongated O-O distance (2.17 A) and was only 12.2 kcal/mol [UB3LYP/6-311+G(3df,2p)+ZPVE] higher in energy than its ground-state precursor. A trans-metastable singlet (trans-2) was 10.9 kcal/mol higher in energy than ground-state HO-ONO. CASSCF(12,10)/6-311+G(d,p) calculations predict the optimized geometries of these cis- and trans-metastable singlets to be close to those obtained with DFT. Optimization of cis- and trans-2 within the COSMO solvent model suggests that both exist as energy minima in polar media. Both cis- and trans-2 exist as hydrogen bonded complexes with several water molecules. These collective data suggest that solvated forms of cis-2.3H(2)O and trans-2.3H(2)O represent the elusive higher lying biradicaloid minima that were recently (J. Am. Chem. Soc. 2003, 125, 16204) advocated as the metastable forms of peroxynitrous acid (HOONO). The involvement of metastable trans-2 in the gas phase oxidation of methane and isobutane is firmly established to take place on the unrestricted [UB3LYP/6-311+G(d,p)] potential energy surface (PES) with classical activations barriers for the hydrogen abstraction step that are 15.7 and 5.9 kcal/mol lower than the corresponding activation energies for producing products methanol and tert-butyl alcohol formed on the restricted PES. The oxidation of dimethyl sulfide and dimethyl selenide, two-electron oxidations, proceeds by an S(N)2-like attack of the heteroatom lone pair on the O-O bond of ground-state peroxynitrous acid. No involvement of metastable forms of HO-ONO was discernible.     

Characterization of singlet ground and low-lying electronic excited states of phosphaethyne and isophosphaethyne

The singlet ground ((approximate)X(1)Sigma1+) and excited (1Sigma-,1Delta) states of HCP and HPC have been systematically investigated using ab initio molecular electronic structure theory. For the ground state, geometries of the two linear stationary points have been optimized and physical properties have been predicted utilizing restricted self-consistent field theory, coupled cluster theory with single and double excitations (CCSD), CCSD with perturbative triple corrections [CCSD(T)], and CCSD with partial iterative triple excitations (CCSDT-3 and CC3). Physical properties computed for the global minimum ((approximate)X(1)Sigma+HCP) include harmonic vibrational frequencies with the cc-pV5Z CCSD(T) method of omega1=3344 cm(-1), omega2=689 cm(-1), and omega3=1298 cm(-1). Linear HPC, a stationary point of Hessian index 2, is predicted to lie 75.2 kcal mol(-1) above the global minimum HCP. The dissociation energy D0[HCP((approximate)X(1)Sigma+)-->H(2S)+CP(X2Sigma+)] of HCP is predicted to be 119.0 kcal mol(-1), which is very close to the experimental lower limit of 119.1 kcal mol(-1). Eight singlet excited states were examined and their physical properties were determined employing three equation-of-motion coupled cluster methods (EOM-CCSD, EOM-CCSDT-3, and EOM-CC3). Four stationary points were located on the lowest-lying excited state potential energy surface, 1Sigma- -->1A", with excitation energies Te of 101.4 kcal mol(-1) (1A"HCP), 104.6 kcal mol(-1)(1Sigma-HCP), 122.3 kcal mol(-1)(1A" HPC), and 171.6 kcal mol(-1)(1Sigma-HPC) at the cc-pVQZ EOM-CCSDT-3 level of theory. The physical properties of the 1A" state with a predicted bond angle of 129.5 degrees compare well with the experimentally reported first singlet state ((approximate)A1A"). The excitation energy predicted for this excitation is T0=99.4 kcal mol(-1) (34 800 cm(-1),4.31 eV), in essentially perfect agreement with the experimental value of T0=99.3 kcal mol(-1)(34 746 cm(-1),4.308 eV). For the second lowest-lying excited singlet surface, 1Delta-->1A', four stationary points were found with Te values of 111.2 kcal mol(-1) (2(1)A' HCP), 112.4 kcal mol(-1) (1Delta HPC), 125.6 kcal mol(-1)(2(1)A' HCP), and 177.8 kcal mol(-1)(1Delta HPC). The predicted CP bond length and frequencies of the 2(1)A' state with a bond angle of 89.8 degrees (1.707 A, 666 and 979 cm(-1)) compare reasonably well with those for the experimentally reported (approximate)C(1)A' state (1.69 A, 615 and 969 cm(-1)). However, the excitation energy and bond angle do not agree well: theoretical values of 108.7 kcal mol(-1) and 89.8 degrees versus experimental values of 115.1 kcal mol(-1) and 113 degrees. of 115.1 kcal mol(-1) and 113 degrees.     

Mechanisms of the reactions of W and W+ with NOx (x=1, 2): a computational study

The mechanisms of the reactions of W and W+ with NOx (x=1, 2) were studied at the CCSD(T)/[SDD+6-311G(d)]//B3LYP/[SDD+6-31G(d)] level of theory. It was shown that the insertion pathway of the reaction W(7S)+NO2(2A1) is a multistate process, which involves several lower lying electronic states of numerous intermediates and transition states, and leads to oxidation, WO(3Sigma)+NO(2Pi), and/or nitration, WN(4Sigma)+O2(3Sigmag-), of the W-center. Oxidation products WO(3Sigma)+NO(2Pi) lie 87.6 kcal/mol below the reactants, while the nitration channel is only 31.0 kcal/mol exothermic. Furthermore, it was shown that nitration of W with NO2 is kinetically less favorable than its oxidation. The addition-dissociation pathway of the reaction W(7S)+NO2(2A1) proceeds via the octet (ground) state potential energy surface of the reaction, requires 3.3 kcal/mol barrier, and leads exclusively to oxidation products. Calculations show that oxidation of the W+ cation by NO2 is a barrierless process in the gas phase, proceeds exclusively via the insertion pathway, and is exothermic by 82.9 kcal/mol. The nitration of W+ by NO2 is only 14.1 kcal/mol exothermic and could be accessible only under high-temperature conditions. Reactions of M=W/W+ with NO are also barrierless processes in the gas phase and lead to the N-O insertion product NMO, which are 105.4 and 77.4 kcal/mol lower than the reactants for W and W+, respectively.     

Is the protein surrounding the active site critical for hydrogen peroxide reduction by selenoprotein glutathione peroxidase? An ONIOM study

In this ONIOM(QM:MM) study, we evaluate the role of the protein surroundings in the mechanism of H2O2 reduction catalyzed by the glutathione peroxidase enzyme, using the whole monomer (3113 atoms in 196 amino acid residues) as a model. A new optimization scheme that allows the full optimization of transition states for large systems has been utilized. It was found that in the presence of the surrounding protein the optimized active site structure bears a closer resemblance to the one in the X-ray structure than that without the surrounding protein. H2O2 reduction occurs through a two-step mechanism. In the first step, the selenolate anion (E-Se(-)) formation occurs with a barrier of 16.4 kcal/mol and is endothermic by 12.0 kcal/mol. The Gln83 residue plays the key role of the proton abstractor, which is in line with the experimental suggestion. In the second step, the O-O bond is cleaved, and selenenic acid (R-Se-OH) and a water molecule are formed. The calculated barrier for this process is 6.0 kcal/mol, and it is exothermic by 80.9 kcal/mol. The overall barrier of 18.0 kcal/mol for H2O2 reduction is in reasonable agreement with the experimentally measured barrier of 14.9 kcal/mol. The protein surroundings has been calculated to exert a net effect of only 0.70 kcal/mol (in comparison to the "active site only" model including solvent effects) on the overall barrier, which is most likely due to the active site being located at the enzyme surface.     

A quantum mechanical study of the reactivity of (SiO)2-defective silica surfaces

The reactivity of the strained (SiO)(2)-four atom ring defect at the silica surfaces has been studied in a cluster approach adopting the ONIOM2[B3LYP6-31+G(d,p):MNDO] method to compute the ring opening reaction by interaction with H(2)O and NH(3). The vibrational "fingerprints" of the isolated defect are computed at 921, 930, and 934 cm(-1) in reasonable agreement with experimental evidence on amorphous silica outgassed at T>900 K. The opening of the (SiO)(2)-four-member ring by the considered molecules is exergonic and the actual value depends on the possible constraints enforced on the reaction products by the silica surrounding. The free kinetic energy barriers result from the interplay between the nucleophilic/electrophilic character of the adsorbed molecule and are 22 and 25 kcal mol(-1) for NH(3) and H(2)O, respectively. All free energy profiles envisage an activated complex in which the nucleophilic part of the molecule interacts on the coordinatively strained silicon atom of the (SiO)(2) defect followed by the proton transfer from the coordinated molecule towards the oxygen of the defective ring. Calculations show that this step can be speed up by the presence of more than one adsorbed molecule or even more (about seven orders of magnitude), by the copresence of water molecules acting as "proton transfer helpers." In these cases, the free energy barriers decrease to approximately 13 and 15 kcal mol(-1) for NH(3) and H(2)O, respectively. For the case of H(2)O adsorption, benchmark test calculations reveal that MP2, BLYP, and B3LYP energy profiles are in very good agreement with each other, whereas for PBE, both the reaction energy and the activation barrier are underestimated. Present data also show that the molecular model mimicking the (SiO)(2) defect is far less reactive than what appears to occur on the real defect at the surface of amorphous silica. So, only a combination of some further geometrical strains imparted by the solid on the (SiO)(2) defect, not accounted for by the cluster models, and higher adsorbate loadings are needed to reharmonize experiment and simulation. Notwithstanding, the vibrational features of the reaction products have been characterized and support the available experimental measurements.     

Singlet and triplet potential surfaces for the O2+C2H4 reaction

Electronic structure calculations at the CASSCF and UB3LYP levels of theory with the aug-cc-pVDZ basis set were used to characterize structures, vibrational frequencies, and energies for stationary points on the ground state triplet and singlet O(2)+C(2)H(4) potential energy surfaces (PESs). Spin-orbit couplings between the PESs were calculated using state averaged CASSCF wave functions. More accurate energies were obtained for the CASSCF structures with the MRMP2/aug-cc-pVDZ method. An important and necessary aspect of the calculations was the need to use different CASSCF active spaces for the different reaction paths on the investigated PESs. The CASSCF calculations focused on O(2)+C(2)H(4) addition to form the C(2)H(4)O(2) biradical on the triplet and singlet surfaces, and isomerization reaction paths ensuing from this biradical. The triplet and singlet C(2)H(4)O(2) biradicals are very similar in structure, primarily differing in their C-C-O-O dihedral angles. The MRMP2 values for the O(2)+C(2)H(4)→C(2)H(4)O(2) barrier to form the biradical are 33.8 and 6.1 kcal/mol, respectively, for the triplet and singlet surfaces. On the singlet surface, C(2)H(4)O(2) isomerizes to dioxetane and ethane-peroxide with MRMP2 barriers of 7.8 and 21.3 kcal/mol. A more exhaustive search of reaction paths was made for the singlet surface using the UB3LYP/aug-cc-pVDZ theory. The triplet and singlet surfaces cross between the structures for the O(2)+C(2)H(4) addition transition states and the biradical intermediates. Trapping in the triplet biradical intermediate, following (3)O(2)+C(2)H(4) addition, is expected to enhance triplet→singlet intersystem crossing.     

The 1deltag dioxygen ene reaction with propene: a density functional and multireference perturbation theory mechanistic study

This study aims to determine whether a balance between concerted and non-concerted pathways exists, and in particular to ascertain the possible role of diradical/zwitterion or peroxirane intermediates. Three non-concerted pathways, via 1) diradical or 2) peroxirane intermediates, and 3) by means of hydrogen-abstraction/radical recoupling, plus one concerted pathway (4), are explored. The intermediates and transition structures (TS) are optimized at the DFT(MPW1K), DFT(B3LYP) and CASSCF levels of theory. The latter optimizations are followed by multireference perturbative CASPT2 energy calculations. (1) The polar diradical forms from the separate reactants by surmounting a barrier (deltaE(++)(MPW1K)=12, deltaE++(B3LYP)=14, and deltaE(++)(CASPT2)=16 kcal mol(-1) and can back-dissociate through the same TS, with barriers of 11 (MPW1K) and 8 kcal mol(-1) (B3LYP and CASPT2). The diradical to hydroperoxide transformation is easy at all levels (deltaE(++)(MPW1K)<4, deltaE(++)(B3LYP)=1 and deltaE(++)(CASPT2)=1 kcal mol(-1)). (2) Peroxirane is attainable only by passing through the diradical intermediate, and not directly, due to the nature of the critical points involved. It is located higher in energy than the diradical by 12 kcal mol(-1), at all theory levels. The energy barrier for the diradical to cis-peroxirane transformation (deltaE(++)=14-16 kcal mol(-1)) is much higher than that for the diradical transformation to the hydroperoxide. In addition, peroxirane can very easily back-transform to the diradical (deltaE(++)<3 kcal mol(-1)). Not only the energetics, but also the qualitative features of the energy hypersurface, prevent a pathway connecting the peroxirane to the hydroperoxide at all levels of theory. (3) The last two-step pathway (hydrogen-abstraction by (1)O(2), followed by HOO-allyl radical coupling) is not competitive with the diradical mechanism. (4) A concerted pathway is carefully investigated, and deemed an artifact of restricted DFT calculations. Finally, the possible ene/[pi2+pi2] competition is discussed.     

Computational study on kinetics and mechanisms of unimolecular decomposition of succinic acid and its anhydride

The mechanisms and kinetics of unimolecular decomposition of succinic acid and its anhydride have been studied at the G2M(CC2) and microcanonical RRKM levels of theory. It was shown that the ZsgsZ conformer of succinic acid, with the Z-acid form and the gauche conformation around the central C-C bond, is its most stable conformer, whereas the lowest energy conformer with the E-acid form, ECGsZ, is only 3.1 kcal/mol higher in energy than the ZsgsZ. Three primary decomposition channels of succinic acid producing H2O + succinic anhydride with a barrier of 51.0 kcal/mol, H2O + OCC2H3COOH with a barrier of 75.7 kcal/mol and CO2 + C2H5COOH with a barrier of 71.9 kcal/mol were predicted. The dehydration process starting from the ECGCZ-conformer is found to be dominant, whereas the decarboxylation reaction starting from the ZsgsZ-conformer is only slightly less favorable. It was shown that the decomposition of succinic anhydride occurs via a concerted fragmentation mechanism (with a 69.6 kcal/mol barrier), leading to formation of CO + CO2 + C2H4 products. On the basis of the calculated potential energy surfaces of these reactions, the rate constants for unimolecular decomposition of succinic acid and its anhydride were predicted. In addition, the predicted rate constants for the unimolecular decomposition of C2H5COOH by decarboxylation (giving C2H6 + CO2) and dehydration (giving H3CCHCO + H2O) are in good agreement with available experimental data.     

Novel pathways for oxygen insertion into unactivated C-H bonds by dioxiranes. Transition structures for stepwise routes via radical pairs and comparison with the concerted pathway

The oxygen insertion into C-H bonds (of methane, isobutane, and acetone) by dioxiranes (parent dioxirane and dimethyldioxirane) to give alcohols was studied with the DFT theory, using both restricted and unrestricted B3LYP methods, and 6-31G(d) and 6-311+G(d,p) basis sets to evaluate the feasibility of stepwise mechanisms and their competition with the concerted counterpart. Confirming previous results by other authors, we have located, with the RB3LYP method, concerted TSs in which the oxygen bound to be inserted interacts very strongly with the hydrogen atom and very weakly with the carbon atom of the C-H bond. These TSs nicely explain all the experimental observations (e.g., configuration retention at the chiral centers), but all of them exhibit an RHF --> UHF wave function instability that preclude considering them as genuine transition structures. We also were able to characterize, with UB3LYP methods, two alternative two-step processes that can lead to final products (alcohol + carbonyl compound) via singlet radical pair intermediates. For the first step of both processes we located genuine diradicaloid TSs, namely, TSs rad,coll and TSs rad,perp, that have stable wave functions. In TSs rad,coll the alkane C-H bond tends to be collinear with the breaking O(1)- - -O(2) bond while in TSs rad,perp the alkane C-H bond is almost perpendicular to the O(1)- - -O(2) bond. The first step, of both processes, can represent an example of a "molecule induced homolysis" reaction: collision between alkane and dioxirane brings about the homolytic cleavage of the dioxirane O-O bond and the hydrogen abstraction follows afterward to produce the diradicaloid TS that then falls down to a singlet radical pair. This hypothesis was fully confirmed by IRC analysis in the case of TSs rad,coll. The possible pathways that lead from the intermediate radical pair to final products are discussed as well as the hypothesis that the radical collinear TSs may collapse directly to products in a "one-step nonconcerted" process. However, diradical mechanisms cannot explain the experimental data as satisfactorily as the concerted pathway does. As for computational predictions about competition of diradical vs concerted mechanisms, they strongly depend (i) on the alkane C-H type, (ii) on whether gas phase or solution is considered, and (iii) on the basis set used for calculations. In short, the concerted TS benefits, with respect to the corresponding diradicaloid TSs, of alkyl substitution at the C-H center, solvation effects, and basis set extension. Actually, in the case of DMD reactions with methane and acetone, the diradicaloid TSs are always (both in gas phase and in solution and with both the basis sets used) strongly favored over their concerted counterpart. In the case of DMD reaction with isobutane tertiary C-H bond the large favor for the diradicaloid TSs over the concerted TS, predicted in gas phase by the B3LYP/6-31G(d) method, progressively decreases as a result of basis set extension and introduction of solvent effects: the higher theory level [B3LYP/6-311+G(d,p)] suggests that in acetone solution TS conc has almost the same energy as TS rad,perp while TS rad,coll resides only 2 kcal/mol higher.