Highly Strained Pn(CH)3 (Pn = N, P, As, Sb, Bi) Tetrahedranes: Theoretical Characterization.

Recent experimental research by Cummins and co-workers has established the existence of a tetrahedrane molecule with one CH moiety replaced by phosphorus. We present here the first theoretical studies of the entire Pn(CH)3 (Pn = N, P, As, Sb, Bi) class of molecules. Geometries are obtained at the highly reliable CCSD(T)/aug-cc-pwCVTZ(-PP) level of theory. Harmonic vibrational frequencies are determined and analyzed to confirm the nature of each stationary point and provide helpful findings that may aid in the detection of each species. Most notable is the result that the geometric parameters associated with the (CH)3 moiety in the tetrahedranes exhibit little change under pnictogen substitution, while the Pn-C bonds and C-Pn-C bond angles greatly increase and decrease, respectively. Strain energies are predicted and range from 122.3 kcal mol-1 (N(CH)3) to 99.4 kcal mol-1 (Bi(CH)3) at the DF-CCSD(T)//B3LYP-D3/aug-cc-pV(T+d)Z(-PP) level of theory. The obtained geometries are further analyzed with Natural Bond Orbital (NBO) methods to understand the bonding and electronic structure of each species. We also provide insight into how different substituents can help make the tetrahedrane structure more energetically favorable due to electron delocalization into substituent antibonding orbitals. The effect of additional delocalization also weakens the Pn-C bonds, especially for the heavier pnictogens. This work concludes with a list of considerations that summarize our key findings and motivate future work aimed at producing novel pnictogen-substituted tetrahedrane molecules.

The Recently Synthesized Dimagnesiabutadiene and the Analogous Dimetalla-Beryllium, -Calcium, -Strontium, and -Barium Compounds.

The 2014 synthesis of the remarkable dimagnesium compound Mg2 [C4 (CH3 )2 (Si(CH3 )3 )2 ](C3 H7 )2 (C4 H8 O)2 may point the way to a new chapter in alkaline earth organometallic chemistry. Accordingly, we have studied the known Mg compound and the analogous Be, Ca, Sr, and Ba structures. Although most of our theoretical predictions come from density functional methods, the latter have been benchmarked using coupled cluster theory including single, double, and perturbative triplet excitations, CCSD(T) using cc-pVTZ basis sets. Among our most important predictions are the energies for dissociation to the butadiene plus the RM-MR [R=(C3 H7 )2 (C4 H8 O)2 ; M=Be, Mg, Ca, Si, and Ba] entities. The most reliable predictions for the dissociation energies are 99-104 (Be), 85-93 (Mg), 90-99 (Ca), 83-92 (Sr), and 83-94 (Ba) kcal mol-1 . Thus, there is reason to anticipate that the four unknown compounds should be achievable synthetically. The predicted metal-metal distances (not single bonds) are 2.89 Š(Mg⋅⋅⋅Mg), 3.46 Š(Ca⋅⋅⋅Ca), 3.75 Š(Sr⋅⋅⋅Sr), and 4.04 Š(Ba⋅⋅⋅Ba). The separated RM-MR compounds have longer M-M distances but genuine metal-metal single bonds. This perhaps counter intuitive result is due to the presence of the bridging carbons in the alkaline earth butadiene compounds. All five compounds incorporate metal-carbon ionic interactions.

Characterizing a nonclassical carbene with coupled cluster methods: cyclobutylidene.

Carbenes represent a special class of reactive compounds that possess a lone pair of electrons on a carbon atom. Among the myriad examples of carbenes in the literature, cyclobutylidene stands out as a unique nonclassical compound that includes transannular interaction between opposing C1 and C3 carbon atoms within a four-membered ring. On its lowest potential energy surface (X[combining tilde](1)A'), cyclobutylidene quickly rearranges, following three reaction paths: (i) 1,2-H migration; (ii) 1,2-C migration; and, (iii) 1,3-H migration. Herein, this reactivity is examined with high-level coupled-cluster methods [up to CCSDT(Q)]. At this level of theory, combined with extrapolation techniques to obtain energies at the complete basis set (CBS) limit, the long-standing disparity between theoretical and experimental results is resolved. Specifically, cyclobutylidene is predicted to prefer 1,2-C migration rather than 1,2-H migration. Rate constants for the three reaction paths are obtained from canonical variational transition state theory (CVT) and yield reasonable agreement with existing experimental results. Further characterization of cyclobutylidene is also reported: the singlet-triplet gap (ΔES-T) is found to be -9.3 kcal mol(-1) at the CCSDT(Q)/CBS level of theory, and anharmonic vibrational frequencies are determined with second-order vibrational perturbation theory (VPT2).

The methylsulfinyl radical CH3SO examined.

Methylsulfinyl radical, a key intermediate in marine atmospheric chemistry, plays a central role in the oxidation of dimethyl sulfide. CH3SO has been extensively studied here with ab initio quantum mechanical methods, with methods as complete as CCSDT(Q) in conjunction with basis sets as large as cc-pV(5+d)Z. In this research, we report high-level computations for the ground and first excited electronic states of the methylsulfinyl radical. The structures of the X[combining tilde] (2)A'' and à (2)A' states are quite different with S-O distances of 1.499 and 1.652 Å, respectively. The X[combining tilde] to à adiabatic energy difference is predicted to be 45.1 kcal mol(-1), compared to 21.1 kcal mol(-1) for the analogous well-characterized methylperoxy radical CH3OO. The CH3SO barrier to internal rotation is 0.92 kcal mol(-1). The unknown X[combining tilde] (2)A'' torsional vibrational frequency τ is predicted to be 142 cm(-1) (harmonic) and 128 cm(-1) (anharmonic). Our predictions of the à (2)A' excited state vibrational frequencies are the first to be reported.

Exploring the effect of axial ligand substitution (X = Br, NCS, CN) on the photodecomposition and electrochemical activity of [MnX(N-C)(CO)3] complexes.

The synthesis, electrochemical activity, and relative photodecomposition rate is reported for four new Mn(i) N-heterocyclic carbene complexes: [MnX(N-ethyl-N'-2-pyridylimidazol-2-ylidine)(CO)3] (X = Br, NCS, CN) and [MnCN(N-ethyl-N'-2-pyridylbenzimidazol-2-ylidine)(CO)3]. All compounds display an electrocatalytic current enhancement under CO2 at the potential of the first reduction, which ranges from -1.53 V to -1.96 V versus the saturated calomel electrode. Catalytic CO production is observed for all species during four-hour preparative-scale electrolysis, but substantial H2 is detected in compounds where X is not Br. All species eventually decompose under both 350 nm and 420 nm light, but cyanide substituted complexes (X = CN) last significantly longer (up to 5×) under 420 nm light as a result of a blue-shifted MLCT band.

Butterfly and rhombus structures for binuclear cobalt carbonyl sulfur and phosphinidene complexes of the type Co2(CO)6E2 (E = S, PX).

Theoretical studies on Co(2)(CO)(6)(PX)(2) derivatives (X = H, Cl, OH, OMe, NH(2), NMe(2)) predict the lowest energy structures to be butterfly structures containing five two-electron two-center bonds in the central Co(2)P(2) unit. Among these butterfly structures the energy increases as the unique bond forming the "body" of the butterfly changes from Co-Co to Co-P and then P-P. Higher energy rhombus structures are also found for Co(2)(CO)(6)(PX)(2) with only Co-P bonds in the Co(2)P(2) framework without any Co-Co or P-P bonds. In addition, for Co(2)(CO)(6)(POR)(2) (R = H, Me) still higher energy "diphosphine" structures are also found containing only three rather than four Co-P bonds, one P-P bond, and no Co...Co bond. For the isoelectronic Co(2)(CO)(6)S(2) a rhombus structure is competitive in energy with the butterfly structures with five structures lying within approximately 4 kcal/mol thereby predicting a fluxional system. A tetrahedrane structure was not found for Co(2)(CO)(6)S(2) in contrast to the tetrahedrane structure known experimentally for the related Fe(2)(CO)(6)S(2) with one less electron per metal atom.

Binuclear manganese and rhenium carbonyls M2(CO)n (n = 10, 9, 8, 7): comparison of first row and third row transition metal carbonyl structures.

The equilibrium geometries, thermochemistry, and vibrational frequencies of the homoleptic binuclear rhenium carbonyls Re2(CO)n (n = 10, 9, 8, 7) were determined using the MPW1PW91 and BP86 methods from density functional theory (DFT) with the effective core potential basis sets LANL2DZ and SDD. In all cases triplet structures for Re2(CO)n were found to be unfavorable energetically relative to singlet structures, in contrast to corresponding Mn2(CO)n derivatives, apparently owing to the larger ligand field splitting of rhenium. For M2(CO)10 (M = Mn, Re) the unbridged structures (OC)5M-M(CO)5 are preferred energetically over structures with bridging CO groups. For M2(CO)9 (M = Mn, Re) the two low energy structures are (OC)4M(micro-CO)M(CO)4 with an M-M single bond and a four-electron donor bridging CO group and (OC)4M[double bond, length as m-dash]M(CO)5 with no bridging CO groups and an M[double bond, length as m-dash]M distance suggesting a double bond. The lowest energy structures for Re2(CO)8 have Re[triple bond, length as m-dash]Re distances in the range 2.6-2.7 A suggesting the triple bonds required to give the Re atoms the favored 18-electron configuration. Low energy structures for Re2(CO)7 are either of the type (OC)(4)M[triple bond, length as m-dash]M(CO)3 with short metal-metal distances suggesting triple bonds or have a single four-electron donor bridging CO group and longer M-M distances consistent with single or double bonds. The 18-electron rule thus appears to be violated in these highly unsaturated Re2(CO)7 structures.

Understanding electron attachment to the DNA double helix: the thymidine monophosphate-adenine pair in the gas phase and aqueous solution.

Electron attachment to the 2'-deoxythymidine-5'-monophosphate-adenine pairs (5'-dTMPH-A and 5'-dTMP(-)-A) has been investigated at a carefully calibrated level of theory (B3LYP/DZP++) to investigate the electron-accepting properties of thymine (T) in the DNA double helix under physiological conditions. All molecular structures have been fully optimized in vacuo and in solution. The adiabatic electron affinity of 5'-dTMPH-A in the gas phase has been predicted to be 0.67 eV. Solvent effects greatly increase the electron capture ability of 5'-dTMPH-A. In fact, the adiabatic electron affinity increases to 2.04 eV with solvation. The influence of the solvent environment on the electron-attracting properties of 5'-dTMPH-A arises not only from the stabilization of the corresponding radical anion through charge-dipole interactions, but also by changing the distribution of the unpaired electron in the molecular system. The unpaired electron is covalently bound even during vertical attachment, due to the solvent effects. Solvent effects also weaken the pairing interaction in the thymidine monophosphate-adenine complexes. The phosphate deprotonation is found to have a relatively minor influence on the capture of electrons by the 5'-dTMPH-A species in aqueous solution. The electron distributions, natural population analysis, and geometrical features of the models examined illustrate that the influence of the phosphate deprotonation is limited to the phosphate moiety in aqueous solution. Therefore, it is reasonable to expect that electron attachment to nucleotides will be independent of monovalent counterions in the vicinity of the phosphate group in aqueous solution.

Effects of fluorine on the structures and energetics of the propynyl and propargyl radicals and their anions.

[reaction: see text] The adiabatic electron affinity (EA(ad)) of the CH(3)-C[triple bond]C(*) radical [experiment = 2.718 +/- 0.008 eV] and the gas-phase basicity of the CH(3)-C[triple bond]C:(-) anion [experiment = 373.4 +/- 2 kcal/mol] have been compared with those of their fluorine derivatives. The latter are studied using theoretical methods. It is found that there are large effects on the electron affinities and gas-phase basicities as the H atoms of the alpha-CH(3) group in the propynyl system are substituted by F atoms. The predicted electron affinities are 3.31 eV (FCH(2)-C[triple bond]C(*)), 3.86 eV (F(2)CH-C[triple bond]C(*)), and 4.24 eV (F(3)C-C[triple bond]C(*)), and the predicted gas-phase basicities of the fluorocarbanion derivatives are 366.4 kcal/mol (FCH(2)-C[triple bond]C:(-)), 356.6 kcal/mol (F(2)CH-C[triple bond]C:(-)), and 349.8 kcal/mol (F(3)C-C[triple bond]C:(-)). It is concluded that the electron affinities of fluoropropynyl radicals increase and the gas-phase basicities decrease as F atoms sequentially replace H atoms of the alpha-CH(3) in the propynyl system. The propargyl radicals, lower in energy than the isomeric propynyl radicals, are also examined and their electron affinities are predicted to be 0.98 eV ((*)CH(2)-C[triple bond]CH), 1.18 eV ((*)CFH-C[triple bond]CH), 1.32 eV ((*)CF(2)-C[triple bond] CH), 1.71 eV ((*)CH(2)-C[triple bond]CF), 2.05 eV ((*)CFH-C[triple bond]CF), and 2.23 eV ((*)CF(2)-C[triple bond]CF).

The dichotomy of dimetallocenes: coaxial versus perpendicular dimetal units in sandwich compounds.

Theoretical studies of the dimetallocene (eta5-C5H5)2Zn2 lead to optimized D5h or D5d structures in which the Zn-Zn bond is coaxial with the C5 axes of the two Cp rings, with a Zn-Zn distance of 2.33 A, corresponding to a Zn-Zn single bond. (eta5-C5H5)2Ni2 (singlet state) and (eta5-C5H5)2Cu2 (triplet) have similar structures with a NiNi triple bond (2.06 A) and a Cu=Cu double bond (2.22 A). However, DFT computations on (C5H5)2Ni2 and (C5H5)2Cu2 (both singlet states) lead to a totally different type of optimized structure (Ci symmetry) lying at significantly lower energies, with the metal-metal bonds perpendicular to the C5 axes of the Cp rings.

The 2'-deoxyadenosine-5'-phosphate anion, the analogous radical, and the different hydrogen-abstracted radical anions: molecular structures and effects on DNA damage.

The 2'-deoxyadenosine-5'-phosphate (5'-dAMP) anion and its related radicals have been studied by reliably calibrated theoretical approaches. This study reveals important physical characteristics of 5'-dAMP radical related processes. One-electron oxidation of the 5'-dAMP anion is found on both the phosphoryl group and the adenine base with electron detachment energies close to that of phosphate. Partial removal of electron density from the adenine fragment leads to an extended pi system which includes the amine group of the adenine. Although the radical-centered carbon increases the extent of bonding with its adjacent atoms, it usually weakens the chemical bonds between the atoms at the alpha- and beta-positions. This tendency should be important in predicting the reactivity of the sugar-based radicals. The overall stability sequence of the H-abstracted 5'-dAMP anionic radicals is consistent with the analogous results for the H-abstracted neutral radicals of the adenosine nucleoside: aliphatic radicals > aromatic radicals. The negatively charged phosphoryl group attached to atom C(5)' of the ribose does not change this energetic sequence. All the H-abstraction produced 5'-dAMP radical anions are distonic radical anions. Studies have shown that the charge-radical-separating feature of the distonic radical anions is biologically relevant. This result should be important in understanding the reactive properties of these H-abstraction-produced anion radicals.

Assessing alkyl-, silyl-, and halo-substituent effects on the electron affinities of silyl radicals.

Neutral anion energy differences for a large class of alpha-substituted silyl radicals have been computed to determine the effect of alkyl, silyl, and halo substituents on their electron affinities. In particular, we report theoretical predictions of the adiabatic electron affinities (AEAs), vertical electron affinities (VEAs), and vertical detachment energies (VDEs) for a series of methyl-, silyl-, and halo-substituted silyl radical compounds. This work utilizes the carefully calibrated DZP++ basis set, in conjunction with the pure BLYP and OLYP functionals, as well as with the hybrid B3LYP, BHLYP, PBE1PBE, MPW1K, and O3LYP functionals. Bromine has the largest effect in stabilizing the anions, and the BLYP/DZP++ AEA for SiBr(3) is 3.29 eV. The other predicted electron affinities are for SiH(3) (1.37 eV), SiH(2)CH(3) (1.09 eV), SiH(2)F (1.54 eV), SiH(2)Cl (1.94 eV), SiH(2)Br (2.05 eV), SiH(2)(SiH(3)) (1.77 eV), SiH(CH(3))(2) (0.92 eV), SiHF(2) (1.86 eV), SiHCl(2) (2.53 eV), SiHBr(2) (2.67 eV), Si(CH(3))(3) (0.86 eV), SiF(3) (2.66 eV), SiCl(3) (3.21 eV), Si(SiH(3))(3) (2.25 eV), and SiFClBr (3.13 eV). For the five silyl radicals where experimental data are available, the BLYP functional gives the most accurate determination of AEAs; the average absolute error is 0.04(1) eV, whereas the corresponding errors for the O3LYP, MPW1K, PBE1PBE, B3LYP, OLYP, and BHLYP functionals are 0.05(8), 0.06(0), 0.06(3), 0.08(5), 0.11(5), and 0.15(3) eV, respectively.

DNA nucleosides and their radical anions: molecular structures and electron affinities.

The deoxyribonucleosides have been studied to determine the properties of combinations of 2-deoxyribose with each of the isolated DNA bases for both neutral and anionic species. We have used a carefully calibrated theoretical method [Chem. Rev. 2002, 102, 231], employing the B3LYP hybrid Hartree-Fock/DFT functional with the DZP++ basis set. Predictions are made of the geometric parameters, adiabatic electron affinities, charge distributions based on natural population analysis, and decomposition enthalpy for the neutral and anionic forms of the four 2'-deoxyribonucleosides in DNA: 2'-deoxyriboadenosine (dA), 2'-deoxyribocytidine (dC), 2'-deoxyriboguanosine (dG), and 2'-deoxyribothymidine (dT). Geometric changes in the anions show that the glycosidic bond exhibits little change with excess charge for the guanosine but significant shortening for the adenosine and for the pyrimidines. The zero-point corrected adiabatic electron affinities in eV for each of the 2'-deoxyribonucleosides are as follows: 0.06, dA; 0.09, dG; 0.33, dC; and 0.44, dT. These values are uniformly greater than those of the corresponding isolated bases (-0.28, A; -0.07, G; 0.03, C; 0.20, T) and the isolated 2-deoxyribose (-0.38) at the same level of theory. The vertical detachment energies of dT and dC are substantial, 0.72 and 0.94 eV, and these anions should be observable. A high VDE, 0.91 eV, is also found for dA but its anion is unlikely to be stable due to the small AEA of 0.06 eV. The high VDE reflects the fact that the molecular structures of the anions and the corresponding neutral species are quite different. Valence character is displayed for the SOMOs of dA, dC, and dT, while some component of diffuse character is visible in the SOMO of dG. Further analysis of electronic changes upon electron attachment include an examination of the NPA charges, which show that in the neutral 2'-deoxyribonucleosides the sum of NPA charges for every base is the same, -0.28 with the sum of 2-deoxyribose charges being positive, +0.28. In the anions, the trend in charge division varies based on the nature of the excess electron in the anions. Thermodynamically, the overall enthalpy change for the reaction of water with the neutral nucleosides to give bases and ribose is approximately zero. The analogous decomposition is exothermic by 8 to 11 kcal mol-1 for the anions, indicating possible challenges for anionic gas-phase nucleoside exploration in the presence of water.

Definitive ab initio studies of model SN2 reactions CH(3)X+F- (X=F, Cl, CN, OH, SH, NH(2), PH(2)).

The energetics of the stationary points of the gas-phase reactions CH(3)X+F(-)-->CH(3)F+X(-) (X=F, Cl, CN, OH, SH, NH(2) and PH(2)) have been definitively computed using focal point analyses. These analyses entailed extrapolation to the one-particle limit for the Hartree-Fock and MP2 energies using basis sets of up to aug-cc-pV5Z quality, inclusion of higher-order electron correlation [CCSD and CCSD(T)] with basis sets of aug-cc-pVTZ quality, and addition of auxiliary terms for core correlation and scalar relativistic effects. The final net activation barriers for the forward reactions are: E (b/F,F)=-0.8, E (b/F, Cl)=-12.2, E (b/F,OH)=+13.6, E b/F,OH=+16.1, E b/F,SH=+2.8, Eb/F, NH=+32.8, and E b/F,PH =+19.7 kcal x mol(-1). For the reverse reactions E b/F,F= -0.8, Eb/Cl,F =+18.3, E b/CN,F=+12.2, E b/OH,F =-1.8, E b/SH,F =+13.2, E b/NH(2),=-1.5, and E b/PH(2) =+9.6 kcal x mol(-1). The change in energetics between the CCSD(T)/aug-cc-pVTZ reference prediction and the final extrapolated focal point value is generally 0.5-1.0 kcal mol(-1). The inclusion of a tight d function in the basis sets for second-row atoms, that is, utilizing the aug-cc-pV(X+d)Z series, appears to change the relative energies by only 0.2 kcal x mol(-1). Additionally, several decomposition schemes have been utilized to partition the ion-molecule complexation energies, namely the Morokuma-Kitaura (MK), reduced variational space (RVS), and symmetry adapted perturbation theory (SAPT) techniques. The reactant complexes fall into two groups, mostly electrostatic complexes (FCH(3).F(-) and ClCH(3).F(-)), and those with substantial covalent character (NCCH(3).F(-), CH(3)OH.F(-), CH(3)SH.F(-), CH(3)NH(2).F(-) and CH(3)PH(2).F(-)). All of the product complexes are of the form FCH(3).X(-) and are primarily electrostatic.

Odd carbon long linear chains HC2n+1H (n = 4-11): properties of the neutrals and radical anions.

The optimized geometries, adiabatic electron affinities, vertical electron affinities, vertical electron detachment energies (for the anions), and IR-active vibrational frequencies have been predicted for the long linear carbon chains HC(2)(n)()(+1)H (n = 4-11). The B3LYP density functional combined with DZP and TZ2P basis sets was used in this theoretical study. These methods have been extensively calibrated versus experiment for the prediction of electron affinities (Chem. Rev. 2002, 102, 231). The computed physical properties are discussed and compared with the even carbon chains HC(2)(n)()H. The predicted electron affinities form a remarkably regular sequence: 2.12 eV (HC(9)H), 2.42 eV (HC(11)H), 2.66 eV (HC(13)H), 2.85 eV (HC(15)H), 3.01 eV (HC(17)H), 3.14 eV (HC(19)H), 3.25 eV (HC(21)H), and 3.35 eV (HC(23)H). These electron affinities are as much as 0.4 eV higher than those for analogous even carbon chains. The predicted structures display an intermediate cumulene-polyacetylene type of bonding, with the inner carbons appearing cumulenic and the outer carbons polyacetylenic.

Resurrection of neutral tris-homoaromaticity.

Neutral in-plane tris-homoaromaticity is evaluated in tris(bismethano)benzene (15) and modifications of this parent structure in which the pi-orbitals might interact in the plane established by the unsaturated carbon atoms (in-plane conjugation). On the basis of magnetic susceptibility exaltation, nucleus-independent shift (NICS), and aromatic stabilization energy (ASE, evaluated via homodesmotic and isodesmic equations using B3LYP/6-311+G + ZPVE energies, as well as by MM3 and MM4 force field computations), we identified triene 17, a triply bridged analogue of 15, as the system where homoaromaticity is most effective. The NICS(total) in the center of 17 is -30.1 ppm and the diatropic pi-contribution is -18.0 ppm. This structure possesses more than one-third of the aromatic stabilization of benzene and is the best candidate for neutral tris-homoaromaticity ever proposed. The previously described tris-(bismethano)-benzene (15) also shows homoaromaticity but to a smaller extent compared to 17. Structure 18, which is closely related to 17, also is significantly homoaromatic but, as evaluated by MM3, strain partially counteracts the stabilizing effects from homoconjugation. Such a counteracting increase in strain largely cancels or even overwhelms the stabilization from homoconjugation in all other species considered in this study.

Examination of the Stabilities of Group 14 (C, Si, Ge, Sn, Pb) Congeners of Dihydroxycarbene and Dioxirane. Comparison to Formic Acid and Hydroperoxycarbene Congeners.

The relative energetics of four XH(2)O(2) (X = C, Si, Ge, Sn, Pb) isomers, dihydroxycarbene, formic acid, dioxirane, and hydroperoxycarbene, were determined using the BLYP and B3LYP density functionals with DZP and TZ2P basis sets, as well as CCSD and CCSD(T) single-point energies at the BLYP/TZ2P optimized geometries. Relative to dihydroxycarbene, formic acid was 41.8 kcal/mol lower in energy while dioxirane and hydroperoxycarbene were 51.3 and 63.6 kcal/mol higher, respectively, with CCSD(T). Furthermore, using an effective core potential (ECP) the dihydroxy congener was shown to be the most stable isomer for X = Si-Pb. The formic acid and dioxirane congeners become increasingly less stable as one descends group 14. Our results show that divalency is preferred for Si-Pb (dihydroxy congeners are the most stable) but the tetravalent formic acid congeners remain more stable than the hydroperoxy congeners, showing that divalency is not universally preferred among these isomers.

Aza[10]annulene: Next Higher Aromatic Analogue of Pyridine.

The aza[10]annulene isomers (CH)(9)N derived from "twist", "heart", "naphthalene-like", and "azulene-like" [10]annulenes by replacement of a CH group by nitrogen were investigated at the CCSD(T)/DZd//B3LYP/6-311+G level of theory. Except for the twist structure, all aza[10]annulenes studied were found to be aromatic. As the transannular repulsion is reduced in the aza[10]annulenes, the aromatic isomers become competitive energetically with the twist form.

Structure, Spectra, and Reaction Energies of the Aluminum-Nitrogen (HAl-NH)(2) and (H(2)Al-NH(2))(2) Rings and the (HAl-NH)(4) Cluster.

Rings of four-coordinate aluminum and nitrogen are easily synthesized and well studied, as are clusters of four-coordinate Al and N. Only recently, however, have rings that are derivatives of the model compounds (HAl-NH)(n)() (n = 2, 3) with three-coordinate Al and N been synthesized. Ab initio investigations of the structure, bonding, vibrational spectra, and reaction energies for the three-coordinate ring (HAl-NH)(2), the four-coordinate ring (H(2)Al-NH(2))(2), and the four-coordinate cluster (HAl-NH)(4) are presented. Even in the absence of differences in steric factors, the four-membered ring has longer Al-N bonds than either the six-membered ring or the unassociated aluminum amide, H(2)Al-NH(2). This is due to both rehybridization and pi interactions. Theoretical reaction energies for formation of the (HAl-NH)(4) cluster from the (H(2)Al-NH(2))(2) ring are consistent with intermolecular loss of hydrogen, or the necessity of surface catalysis.

Carbene Rearrangements Unsurpassed: Details of the C(7)H(6) Potential Energy Surface Revealed.

The rearrangement of phenylcarbene (1) to 1,2,4,6-cycloheptatetraene (3) has been studied theoretically, using SCF, CASSCF, CASPT2N, DFT (B3LYP), CISD, CCSD, and CCSD(T) methods in conjunction with the 6-31G, 6-311+G, 6-311G(2d,p), cc-pVDZ, and DZd basis sets. Stationary points were characterized by vibrational frequency analyses at CASSCF(8,8)/6-31G and B3LYP/6-31G. Phenylcarbene (1) has a triplet ground state ((3)A") with a singlet-triplet separation (DeltaE(ST)) of 3-5 kcal mol(-)(1). In agreement with experiment, chiral 3 is the lowest lying structure on this part of the C(7)H(6) potential energy surface. Bicyclo[4.1.0]hepta-2,4,6-triene (2) is an intermediate in the rearrangement of 1 into 3, but it is unlikely to be observable experimentally due to a barrier height of only 1-2 kcal mol(-)(1). The enantiomers of 3 interconvert via the (1)A(2) state of cycloheptatrienylidene (4) with an activation energy of 20 kcal mol(-)(1). The "aromatic" (1)A(1) state, previously believed to be the lowest singlet state of 4, is roughly 10 kcal mol(-)(1) higher in energy than the (1)A(2) state, and, in violation of Hund's rule, (3)A(2) is also calculated to lie above (1)A(2) in energy. Thus, even if (3)A(2) were populated, it is likely to undergo rapid intersystem crossing to (1)A(2). We suggest (3)B(1)-4 is the metastable triplet observed by EPR.