Quadricyclane radical cation rearrangements: a computational study of the transformations to 1,3,5-cycloheptatriene and norbornadiene

An alternative skeletal rearrangement of the quadricyclane radical cation (Q*+) explains the side products formed in the one-electron oxidation to norbornadiene. First, the bicyclo[2.2.1]hepta-2-ene-5-yl-7-ylium radical cation, with an activation energy of 14.9 kcal mol(-1), is formed. Second, this species can further rearrange to 1,3,5-cycloheptatriene through two plausible paths, that is, a multistep mechanism with two shallow intermediates and a stepwise path in which the bicyclo[3.2.0]hepta-2,6-diene radical cation is an intermediate. The multistep rearrangement has a rate-limiting step with an estimated activation energy of 16.5 kcal mol(-1), which is 2.8 kcal mol(-1) lower in energy than the stepwise mechanism. However, the lowest activation energy is found for the Q*+ cycloreversion to norbornadiene that has a transition structure, in close correspondence with earlier studies, and an activation energy of 10.1 kcal mol(-1), which agrees well with the experimental estimate of 9.3 kcal mol(-1). The computational estimates of activation energies were done using the CCSD(T)/6-311+G(d,p) method with geometries optimized on the B3LYP/6-311+G(d,p) level, combined with B3LYP/6-311+G(d,p) frequencies.     

The rearrangement of the cubane radical cation in solution

The rearrangement of the cubane radical cation (1*+) was examined both experimentally (anodic as well as (photo)chemical oxidation of cubane 1 in acetonitrile) and computationally at coupled cluster, DFT, and MP2 [BCCD(T)/cc-pVDZ//B3LYP/6-31G* ZPVE as well as BCCD(T)/cc-pVDZ//MP2/6-31G* + ZPVE] levels of theory. The interconversion of the twelve C2v degenerate structures of 1*+ is associated with a sizable activation energy of 1.6 kcalmol(-1). The barriers for the isomerization of 1*- to the cuneane radical cation (2*+) and for the C-C bond fragmentation to the secocubane-4,7-diyl radical cation (10*+) are virtually identical (deltaH0++ = 7.8 and 7.9 kcalmol(-1), respectively). The low-barrier rearrangement of 10*+ to the more stable syn-tricyclooctadiene radical cation 3*+ favors the fragmentation pathway that terminates with the cyclooctatetraene radical cation 6*+. Experimental single-electron transfer (SET) oxidation of cubane in acetonitrile with photoexcited 1,2,4,5-tetracyanobenzene, in combination with back electron transfer to the transient radical cation, also shows that 1*+ preferentially follows a multistep rearrangement to 6*+ through 10*+ and 3*+ rather than through 2*+. This was confirmed by the oxidation of syn-tricyclooctadiene (3), which, like 1, also forms 6 in the SET oxidation/rearrangement/electron-recapture process. In contrast, cuneane (2) is oxidized exclusively to semibullvalene (9) under analogous conditions. The rearrangement of 1*+ to 6*+ via 3*+, which was recently observed spectroscopically upon ionization in a hydrocarbon glass matrix, is also favored in solution.     

A theoretical study of the effect of a tetraalkylammonium counterion on the hydrogen bond strength in Z-hydrogen maleate.

High-level ab initio calculations (B3LYP/6-31+G and QCISD(T)/6-311+G**) were carried out to resolve the disagreement between recent experimental and computational estimates of the relative strength of the intramolecular hydrogen bond in Z-hydrogen maleate anion with respect to the normal hydrogen bond in maleic acid. The computational estimates for the strength of the intramolecular hydrogen bond in the gas-phase maleate anion are in a range of 14-28 kcal/mol depending on the choice of the reference structure. Computational data suggest that the electrostatic influence of a counterion such as a tetraalkylammonium cation can considerably weaken the hydrogen bonding interaction (by 1.5-2 times) in the complexed hydrogen maleate anion relative to that in the naked anion. The estimated internal H-bonding energies for a series of Z-maleate/R4N+ salts (R = CH3, C2H5, CH3CH2CH2CH2) range from 8 to 13 kcal/mol. The calculated energy differences between the E- and Z-hydrogen maleates complexed to Me4N+, Et4N+, and Bu4N+ cation are 4.9 (B3LYP/6-31+G(d,p)) and 5.7 and 5.8 kcal/mol (B3LYP/6-31G(d)). It is also demonstrated that the sodium cation exerts a similar electrostatic influence on the hydrogen bond strength in bifluoride anion (FHF-). The present study shows that while low-barrier short hydrogen bonds can exist in the gas phase (the barrier for the hydrogen transfer in maleate anion is only 0.2 kcal/mol at the QCISD(T)/6-311+G//QCISD/6-31+G level), whether they can also be strong in condensed media or not depends on how their interactions with their immediate environment affect their strength.     

Gaseous reaction mechanism of C2F radical with water

The kinetic properties of the carbon-fluorine radicals are little understood except those of CFn (n =1-3). In this article, a detailed mechanistic study was reported on the gas-phase reaction between the simplest pi-bonded C2F radical and water as the first attempt to understand the chemical reactivity of the C2F radical. Various reaction channels are considered. The most kinetically competitive channel is the quasi-direct hydrogen-abstraction route forming P5 HCCF + OH. At the CCSD(T)/6-311+G(2d,2p)//B3LYP/6-311G(d,p)+ZPVE, CCSD(T)/6-311+G(3df,2p)//QCISD/6-311G(d,p)+ZPVE and Gaussian-3//B3LYP/6-31G(d) levels, the overall H-abstraction barriers (4.5, 4.7, and 4.2 kcal/mol) for the C2F + H2O reaction are comparable to the corresponding values (5.5, 3.7, and 5.7 kcal/mol) for the analogous C2H + H2O reaction. This suggests that C2F is a reactive radical like the extensively studied C2H, in contrast to the situation of the CF and CF2 radicals that have much lower reactivity than the corresponding hydrocarbon species. Thus, the C2F radical is expected to play an important role in the combustion processes of the carbon-fluorine chemistry. Furthermore, addition of a second H2O can catalyze the reaction with the H-abstraction barrier significantly reduced to a marginally zero value (0.5 kcal/mol). This is also indicative of the potential relevance of the title reactions in the low-temperature atmospheric chemistry.     

EPR studies of amine radical cations. Part 2. Thermal and photo-induced rearrangements of propargylamine and allylamine radical cations in low-temperature freon matrices

Matrix EPR studies and quantum chemical calculations have been used to characterize the consecutive H-atom shifts undergone by the nitrogen-centered parent radical cations of propargylamine (1b*+) and allylamine (5*+) on thermal or photoinduced activation. The radical cation rearrangements of these unsaturated parent amines occur initially by a 1,2 H-atom shift from C1 to C2 with pi-bond formation at the positively charged nitrogen; this is followed by a consecutive reaction involving a second H-atom shift from C2 to C3. Thus, exposure to red light (lambda > 650 nm) converts 1b*+ to the vinyl-type distonic radical cation 2*+ which in turn is transformed on further photolysis with blue-green light (lambda approximately 400-600 nm) to the allene-type heteroallylic radical cation 3*+. Calculations show that the energy ordering is 1b*+ > 2*+ > 3*+, so that the consecutive H-atom shifts are driven by the formation of more stable isomers. Similarly, the parent radical cation of allylamine 5*+ undergoes a spontaneous 1,2-hydrogen atom shift from C1 to C2 at 77 K with a t1/2 of approximately 1 h to yield the distonic alkyl-type iminopropyl radical cation 6*+; this thermal reaction is attributed largely to quantum tunneling, and the rate is enhanced on concomitant photobleaching with visible light. Subsequent exposure to UV light (lambda approximately 350-400 nm) converts 6*+ by a 2,3 H-shift to the 1-aminopropene radical cation 7*+, which is confirmed to be the lowest-energy isomer derived from the ionization of either allylamine or cyclopropylamine. Although the parent radical cations of N, N-dimethylallylamine (9*+) and N-methylallylamine (11*+) are both stabilized by the electron-donating character of the methyl group(s), the photobleaching of 9*+ leads to the remarkable formation of the cyclic 1-methylpyrrolidine radical cation 10*+. The first step of this transformation now involves the migration of a hydrogen atom to C2 of the allyl group from one of the methyl groups (rather than from C1); the reaction is then completed by the cyclization of the generated MeN + (=CH2) CH2CH2CH2* distonic radical cation, possibly in a concerted overall process. In contrast to the ubiquitous H-atom transfer from carbon to nitrogen that occurs in the parent radical cations of saturated amines, the alternate rearrangements of either 1b*+ or 5*+ to an ammonium-type radical cation by a hypothetical H-atom shift from C1 to the ionized NH2 group are not observed. This is in line with calculations showing that the thermal barrier for this transformation is much higher (approximately 120 kJ mol-1) than those for the conversion of 1b*+ --> 2*+ and 5*+--> 6*+ (approximately 40-60 kJ mol-1).     

Corannulene as a Lewis base: computational modeling of protonation and lithium cation binding.

A computational modeling of the protonation of corannulene at B3LYP/6-311G(d,p)//B3LYP/6-311G(d,p) and of the binding of lithium cations to corannulene at B3LYP/6-311G(d,p)//B3LYP/6-31G(d,p) has been performed. A proton attaches preferentially to one carbon atom, forming a sigma-complex. The isomer protonated at the innermost (hub) carbon has the best total energy. Protonation at the outermost (rim) carbon and at the intermediate (bridgehead rim) carbon is less favorable by ca. 2 and 14 kcal mol(-)(1), respectively. Hydrogen-bridged isomers are transition states between the sigma-complexes; the corresponding activation energies vary from 10 to 26 kcal mol(-)(1). With an empirical correction obtained from calculations on benzene, naphthalene, and azulene, the best estimate for the proton affinity of corannulene is 203 kcal mol(-)(1). The lithium cation positions itself preferentially over a ring. There is a small energetic preference for the 6-ring over the 5-ring binding (up to 2 kcal mol(-)(1)) and of the convex face over the concave face (3-5 kcal mol(-)(1)). The Li-bridged complexes are transition states between the pi-face complexes. Movement of the Li(+) cation over either face is facile, and the activation energy does not exceed 6 kcal mol(-)(1) on the convex face and 2.2 kcal mol(-)(1) on the concave face. In contrast, the transition of Li(+) around the corannulene edge involves a high activation barrier (24 kcal mol(-)(1) with respect to the lowest energy pi-face complex). An easier concave/convex transformation and vice versa is the bowl-to-bowl inversion with an activation energy of 7-12 kcal mol(-)(1). The computed binding energy of Li(+) to corannulene is 44 kcal mol(-)(1). Calculations of the (7)Li NMR chemical shifts and nuclear independent chemical shifts (NICS) have been performed to analyze the aromaticity of the corannulene rings and its changes upon protonation.     

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.     

The Curtius rearrangement of cyclopropyl and cyclopropenoyl azides. A combined theoretical and experimental mechanistic study

A combined experimental and theoretical study addresses the concertedness of the thermal Curtius rearrangement. The kinetics of the Curtius rearrangements of methyl 1-azidocarbonyl cycloprop-2-ene-1-carboxylate and methyl 1-azidocarbonyl cyclopropane-1-carboxylate were studied by (1)H NMR spectroscopy, and there is close agreement between calculated and experimental enthalpies and entropies of activation. Density functional theory (DFT) calculations (B3LYP/6-311+G(d,p)) on these same acyl azides suggest gas phase barriers of 27.8 and 25.1 kcal/mol. By comparison, gas phase activation barriers for the rearrangement of acetyl, pivaloyl, and phenyl azides are 27.6, 27.4, and 30.0 kcal/mol, respectively. The barrier for the concerted Curtius reaction of acetyl azide at the CCSD(T)/6-311+G(d,p) level exhibited a comparable activation energy of 26.3 kcal/mol. Intrinsic reaction coordinate (IRC) analyses suggest that all of the rearrangements occur by a concerted pathway with the concomitant loss of N2. The lower activation energy for the rearrangement of methyl 1-azidocarbonyl cycloprop-2-ene-1-carboxylate relative to methyl 1-azidocarbonyl cyclopropane-1-carboxylate was attributed to a weaker bond between the carbonyl carbon and the three-membered ring in the former compound. Calculations on the rearrangement of cycloprop-2-ene-1-oyl azides do not support pi-stabilization of the transition state by the cyclopropene double bond. A comparison of reaction pathways at the CBS-QB3 level for the Curtius rearrangement versus the loss of N2 to form a nitrene intermediate provides strong evidence that the concerted Curtius rearrangement is the dominant process.     

Isomerization and fragmentation reactions of gaseous phenylarsane radical cations and phenylarsanyl cations. A study by tandem mass spectrometry and theoretical calculations

The unimolecular reactions of radical cations and cations derived from phenylarsane, C6H5AsH2 (1) and dideutero phenylarsane, C6H5AsD2 (1-d2), were investigated by methods of tandem mass spectrometry and theoretical calculations. The mass spectrometric experiments reveal that the molecular ion of phenylarsane, 1*+, exhibits different reactivity at low and high internal excess energy. Only at low internal energy the observed fragmentations are as expected, that is the molecular ion 1*+ decomposes almost exclusively by loss of an H atom. The deuterated derivative 1-d2 with an AsD2 group eliminates selectively a D atom under these conditions. The resulting phenylarsenium ion [C6H5AsH]+, 2+, decomposes rather easily by loss of the As atom to give the benzene radical cation [C6H6]*+ and is therefore of low abundance in the 70 eV EI mass spectrum. At high internal excess energy, the ion 1*+ decomposes very differently either by elimination of an H2 molecule, or by release of the As atom, or by loss of an AsH fragment. Final products of these reactions are either the benzoarsenium ion 4*+, or the benzonium ion [C6H7]+, or the benzene radical cation, [C6H6]*+. As key-steps, these fragmentations contain reductive eliminations from the central As atom under H-H or C-H bond formation. Labeling experiments show that H/D exchange reactions precede these fragmentations and, specifically, that complete positional exchange of the H atoms in 1*+ occurs. Computations at the UMP2/6-311+G(d)//UHF/6-311+G(d) level agree best with the experimental results and suggest: (i) 1*+ rearranges (activation enthalpy of 93 kJ mol(-1)) to a distinctly more stable (DeltaH(r)(298) = -64 kJ mol(-1)) isomer 1 sigma*+ with a structure best represented as a distonic radical cation sigma complex between AsH and benzene. (ii) The six H atoms of the benzene moiety of 1 sigma*+ become equivalent by a fast ring walk of the AsH group. (iii) A reversible isomerization 1+<==>1 sigma*+ scrambles eventually all H atoms over all positions in 1*+. The distonic radical cation 1*+ is predisposed for the elimination of an As atom or an AsH fragment. The calculations are in accordance with the experimentally preferred reactions when the As atom and the AsH fragment are generated in the quartet and triplet state, respectively. Alternatively, 1*(+) undergoes a reductive elimination of H2 from the AsH2 group via a remarkably stable complex of the phenylarsandiyl radical cation, [C6H5As]*+ and an H2 molecule.     

o-Quinone methide as alkylating agent of nitrogen, oxygen, and sulfur nucleophiles. The role of H-bonding and solvent effects on the reactivity through a DFT computational study

The reactivity of the alkylating agent o-quinone methide (o-QM) toward NH(3), H(2)O, and H(2)S, prototypes of nitrogen-, oxygen-, and sulfur-centered nucleophiles, has been studied by quantum chemical methods in the frame of DF theory (B3LYP) in reactions modeling its reactivity in water with biological nucleophiles. The computational analysis explores the reaction of NH(3), H(2)O, and H(2)S with o-QM, both free and H-bonded to a discrete water molecule, with the aim to rationalize the specific and general effect of the solvent on o-QM reactivity. Optimizations of stationary points were done at the B3LYP level using several basis sets [6-31G(d), 6-311+G(d,p), adding d and f functions to the S atom, 6-311+G(d,p),S(2df), and AUG-cc-pVTZ]. The activation energies calculated for the addition reactions were found to be reduced by the assistance of a water molecule, which makes easier the proton-transfer process in these alkylation reactions by at least 12.9, 10.5, and 6.0 kcal mol(-1) [at the B3LYP/AUG-cc-pVTZ//B3LYP/6-311+G(d,p) level], for ammonia, water, and hydrogen sulfide, respectively. A proper comparison of an uncatalyzed with a water-catalyzed reaction mechanism has been made on the basis of activation Gibbs free energies. In gas-phase alkylation of ammonia and water by o-QM, reactions assisted by an additional water molecule H-bonded to o-QM (water-catalyzed mechanism) are favored over their uncatalyzed counterparts by 5.6 and 4.0 kcal mol(-1) [at the B3LYP/6-311+G(d,p) level], respectively. In contrast, the hydrogen sulfide alkylation reaction in the gas phase shows a slight preference for a direct alkylation without water assistance, even though the free energy difference (DeltaDeltaG(#)) between the two reaction mechanisms is very small (by 1.0 kcal mol(-1) at the B3LYP/6-311+G(d,p),S(2df) level of theory). The bulk solvent effect, evaluated by the C-PCM model, significantly modifies the relative importance of the uncatalyzed and water-assisted alkylation mechanism by o-QM in comparison to the case in the gas phase. Unexpectedly, the uncatalyzed mechanism becomes highly favored over the catalyzed one in the alkylation reaction of ammonia (by 7.0 kcal mol(-1)) and hydrogen sulfide (by 4.0 kcal mol(-1)). In contrast, activation induced by water complexation still plays an important role in the o-QM hydration reaction in water as solvent.     

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.     

Nucleofugality of phenyl and methyl carbonates

A series of X,Y-substituted benzhydryl phenyl carbonates 1 and X,Y-substituted benzhydryl methyl carbonates 2 were subjected to solvolysis in different methanol/water, ethanol/water, and acetone/water mixtures at 25 degrees C. The LFER equation, log k = sf(Ef + Nf), was used to derive the nucleofuge-specific parameters (Nf and sf) for phenyl carbonate (1LG) and methyl carbonate (2LG) leaving groups in a given solvent in SN1 type reaction. Kinetic measurements showed that phenyl carbonates solvolyze one order of magnitude faster than methyl carbonates. Optimized geometries of 1LG and 2LG at B3LYP/6-311G(d,p), B3LYP/6-311++G(d,p), and MP2(full)/6-311++G(d,p) levels revealed that negative charge delocalization in carbonate anions to all three oxygen atoms occurs due to negative hyperconjugation. Phenyl carbonate (1LG) is a better leaving group (Nf = -0.84 +/- 0.07 in 80% v/v aq EtOH) than methyl carbonate 2LG (Nf = -1.84 +/- 0.07 in 80% v/v aq EtOH) because of more pronounced negative hyperconjugation, which is characterized with a more elongated RO-C bond and more increased RO-C-CO angle in 1LG than in 2LG. Calculated affinities of benzhydryl cation toward methyl and phenyl carbonate anions (DeltaDeltaEaff = 11.7 kcal/mol at the B3LYP/6-311++G(d,p) level and DeltaDeltaEaff = 2.7 kcal/mol at the PCM-B3LYP/6-311++G(d,p) level in methanol, respectively) showed that 1LG is more stabilized than 2LG, which is in accordance with greater solvolytic reactivity of 1 than 2.     

The guanine cation radical: investigation of deprotonation states by ESR and DFT

This work reports ESR studies that identify the favored site of deprotonation of the guanine cation radical (G*+) in an aqueous medium at 77 K. Using ESR and UV-visible spectroscopy, one-electron oxidized guanine is investigated in frozen aqueous D2O solutions of 2'-deoxyguanosine (dGuo) at low temperatures at various pHs at which the guanine cation radical, G*+ (pH 3-5), singly deprotonated species, G(-H)* (pH 7-9), and doubly deprotonated species, G(-2H)*- (pH > 11), are found. C-8-deuteration of dGuo to give 8-D-dGuo removes the major proton hyperfine coupling at C-8. This isolates the anisotropic nitrogen couplings for each of the three species and aids our analyses. These anisotropic nitrogen couplings were assigned to specific nitrogen sites by use of 15N-substituted derivatives at N1, N2, and N3 atoms in dGuo. Both ESR and UV-visible spectra are reported for each of the species: G*+, G(-H)*, and G(-2H)*-. The experimental anisotropic ESR hyperfine couplings are compared to those obtained from DFT calculations for the various tautomers of G(-H)*. Using the B3LYP/6-31G(d) method, the geometries and energies of G*+ and its singly deprotonated state in its two tautomeric forms, G(N1-H)* and G(N2-H)*, were investigated. In a nonhydrated state, G(N2-H)* is found to be more stable than G(N1-H)*, but on hydration with seven water molecules G(N1-H)* is found to be more stable than G(N2-H)*. The theoretically calculated hyperfine coupling constants (HFCCs) of G*+, G(N1-H)*, and G(-2H)*- match the experimentally observed HFCCs best on hydration with seven or more waters. For G(-2H)*-, the hyperfine coupling constant (HFCC) at the exocyclic nitrogen atom (N2) is especially sensitive to the number of hydrating water molecules; good agreement with experiment is not obtained until nine or 10 waters of hydration are included.     

Structure of the boronic acid dimer and the relative stabilities of its conformers

Despite the widespread use of boronic acids in materials science and as pharmaceutical agents, many aspects of their structure and reactivity are not well understood. In this research the boronic acid dimer, [HB(OH)(2)](2), was studied by second-order M?ller-Plesset (MP2) perturbation theory and coupled cluster methodology with single and double excitations (CCSD). Pople split-valence 6-31+G*, 6-311G**, and 6-311++G** and Dunning-Woon correlation-consistent cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ basis sets were employed for the calculations. A doubly hydrogen-bonded conformer (1) of the dimer was consistently found to be lowest in energy; the structure of 1 was planar (C(2h)) at most computational levels employed but was significantly nonplanar (C(2)) at the MP2/6-311++G** and CCSD/6-311++G** levels, the result of an intrinsic problem with Pople-type sp-diffuse basis functions on heavy atoms. The dimerization energy, enthalpy, and free energy for the formation of (1) from the exo-endo conformer of the monomer were -10.8, -9.2, and +1.2 kcal/mol, respectively, at the MP2/aug-cc-pVTZ level. Several other hydrogen-bonded conformers of the dimer were local minima on the potential energy surface (PES) and ranged from 2 to 5 kcal/mol higher in energy than 1. Nine doubly OH-bridged conformers, in which the boron atoms were tetracoordinated, were also local minima on the PES, but they were all greater than 13 kcal/mol higher in energy than 1; doubly H-bridged structures proved to be transition states. MP2 and CCSD results were compared to those from the BLYP, B3LYP, OLYP, O3LYP, PBE1PBE, and TPSS functionals with the 6-311++G** and aug-cc-pVTZ basis sets; the PBE1PBE functional performed best relative to the MP2 and CCSD results. Self-consistent reaction field (SCRF) calculations predict that boronic acid dimerization is less favorable in solution than in vacuo.     

Kinetics of the multichannel reaction of methanethiyl radical (CH3S*) with 3O2

The CH3S* + O2 reaction system is considered an important process in atmospheric chemistry and in combustion as a pathway for the exothermic conversion of methane-thiyl radical, CH3S*. Several density functional and ab initio computational methods are used in this study to determine thermochemical parameters, reaction paths, and kinetic barriers in the CH3S* + O2 reaction system. The data are also used to evaluate feasibility of the DFT methods for higher molecular weight oxy-sulfur hydrocarbons, where sulfur presents added complexity from its many valence states. The methods include: B3LYP/6-311++G(d,p), B3LYP/6-311++G(3df,2p), CCSD(T)/6-311G(d,p)//MP2/6-31G(d,p), B3P86/6-311G(2d,2p)//B3P86/6-31G(d), B3PW91/6-311++G(3df,2p), G3MP2, and CBS-QB3. The well depth for the CH3S* + 3O2 reaction to the syn-CH3SOO* adduct is found to be 9.7 kcal/mol. Low barrier exit channels from the syn-CH3SOO* adduct include: CH2S + HO2, (TS6, E(a) is 12.5 kcal/mol), CH3 + SO2 via CH3SO2 (TS2', E(a) is 17.8) and CH3SO + O (TS17, E(a) is 24.7) where the activation energy is relative to the syn-CH3SOO* stabilized adduct. The transition state (TS5) for formation of the CH3SOO adduct from CH3S* + O2 and the reverse dissociation of CH3SOO to CH3S* + O2 is relatively tight compared to typical association and simple bond dissociation reactions; this is a result of the very weak interaction. Reverse reaction is the dominant dissociation path due to enthalpy and entropy considerations. The rate constants from the chemical activation reaction and from the stabilized adduct to these products are estimated as functions of temperature and pressure. Our forward rate constant and CH3S loss profile are in agreement with the experiments under similar conditions. Of the methods above, the G3MP2 and CBS-QB3 composite methods are recommended for thermochemical determinations on these carbon-sulfur-oxygen systems, when they are feasible.     

Quantum mechanical study of the ring-closing reaction of the hexatriene radical cation

The ring-closing reaction of hexatriene radical cation 1(*)(+) to 1,3-cyclohexadiene radical cation 2(*)(+) was studied computationally at the B3LYP/6-31G* and QCISD(T)/6-311G*//QCISD/6-31G* levels of theory. Both, concerted and stepwise mechanisms were initially considered for this reaction. Upon evaluation at the B3LYP level of theory, three of the possible pathways-a concerted C(2)-symmetric via transition structure 3(*)(+) and stepwise C(1)-symmetric pathways involving three-membered ring intermediate 5(*)(+) and four-membered ring intermediate 6(*)(+)-were rejected due to high-energy stationary points along the reaction pathway. The two remaining pathways were found to be of competing energy. The first proceeds through the asymmetric, concerted transition structure 4(*)(+) with an activation barrier E(a) = 16.2 kcal/mol and an overall exothermicity of -23.8 kcal/mol. The second pathway, beginning from the cis,cis,trans rotamer of 1(*)(+), proceeds by a stepwise pathway to the cyclohexadiene product with an overall exothermicity of -18.6 kcal/mol. The activation energy for the rate-determining step in this process, the formation of the intermediate bicyclo[3.1.0]hex-2-ene via transition structure 9(*)(+), was found to be 20.4 kcal/mol. More rigorous calculations of a smaller subsection of the potential energy hypersurface at the QCISD(T)//QCISD level confirmed these findings and emphasized the importance of conformational control of the reactant.     

A new scheme for determining the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies of glycine and alanine peptides

In this paper a new scheme was proposed to calculate the intramolecular hydrogen-bonding energies in peptides and was applied to calculate the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies of the glycine and alanine peptides. The density-functional theory B3LYP6-31G(d) and B3LYP6-311G(d,p) methods and the second-order Moller-Plesset perturbation theory MP26-31G(d) method were used to calculate the optimal geometries and frequencies of glycine and alanine peptides and related structures. MP26-311++G(d,p), MP26-311++G(3df,2p), and MP2/aug-cc-pVTZ methods were then used to evaluate the single-point energies. It was found that the B3LYP6-31G(d), MP26-31G(d), and B3LYP6-311G(d,p) methods yield almost similar structural parameters for the conformers of the glycine and alanine dipeptides. MP2/aug-cc-pVTZ predicts that the intramolecular seven-membered ring N-H...O=C hydrogen-bonding strength has a value of 5.54 kcal/mol in glycine dipeptide and 5.73 and 5.19 kcal/mol in alanine dipeptides, while the steric repulsive interactions of the seven-membered ring conformers are 4.13 kcal/mol in glycine dipeptide and 6.62 and 3.71 kcal/mol in alanine dipeptides. It was also found that MP26-311++G(3df,2p) gives as accurate intramolecular N-H...O=C hydrogen-bonding energies and steric repulsive interactions as the much more costly MP2/aug-cc-pVTZ does.     

Theoretical study on structures and stability of Si2CP isomers

The structures, energetics, spectroscopies, and isomerization of various doublet Si2CP species are explored theoretically. In contrast to the previously studied SiC2N and SiC2P radicals that have linear SiCCN and SiCCP ground states, the title Si2CP radical has a four-membered-ring form cSiSiPC 1 (0.0 kcal/mol) with Si-C cross-bonding as the ground-state isomer at the CCSD(T)/6-311G(2df)//B3LYP/6-311G(d)+ZPVE level, similar to the Si2CN radical. The second low-lying isomer 2 at 11.6 kcal/mol has a SiCSiP four-membered ring with C-P cross-bonding, yet it is kinetically quite unstable toward conversion to 1 with a barrier of 3.5 kcal/mol. In addition, three cyclic species with divalent carbene character, i.e., cSiSiCP 7, 7' with C-P cross-bonding and cSiCSiP 8 with Si-Si cross-bonding, are found to possess considerable kinetic stability, although they are energetically high lying at 44.4, 46.5, and 41.4 kcal/mol, respectively. Moreover, a linear isomer SiCSiP 5 at 44.3 kcal/mol also has considerable kinetic stability and predominantly features the interesting cumulenic /Si=C=Si=P/* form with a slight contribution from the silicon-phosphorus triply bonded form /Si=C*-Si[triple bond]P/. The silicon-carbon triply bonded form *Si[triple bond]C-Si[triple bond]P/ has negligible contribution. All five isomers are expected to be observable in low-temperature environments. Their bonding nature and possible formation strategies are discussed. For relevant species, the QCISD/6-311G(d) and CCSD(T)/6-311+G(2df) (single-point) calculations are performed to provide more reliable results. The calculated results are compared to those of the analogous C3N, C3P, SiC2N, and Si2CN radicals with 17 valence electrons. Implications in interstellar space and P-doped SiC vaporization processes are also discussed.     

Singlet hydrocarbon carbenes with high barriers toward isomerization: a computational investigation

A prerequisite for a stable singlet hydrocarbon carbene is the existence of high barriers toward isomerization. Four derivatives of cyclopentylidene (1-4) with rigid and varying carbon cages are examined computationally at the B3LYP/6-311+G(d,p) level of theory. Singlet ground states are predicted for carbenes 1-4, with DeltaE(ST)'s = 7-22 kcal/mol. The rearrangement paths considered are 1,3-hydrogen shift, 1,2-carbon shift and beta-CC bond-cleavage. Carbenes 3 and 4 lie in relatively shallow potential-energy wells (around 4 and 6 kcal/mol, respectively) and are expected to rearrange via 1,3-hydrogen shifts to cyclopropane derivatives. For 1 and 2, the lowest energy rearrangement path is beta-CC bond-cleavage requiring about 12 and 20 kcal/mol, respectively, placing 2 in the deepest potential-energy well among the four carbenes.