Aromaticity in metallabenzenes

The electronic structure and bonding situation in 21 metallabenzenes (metal=Os, Ru, Ir, Rh, Pt, and Pd) were investigated at the DFT level (BP86/TZ2P) by using an energy decomposition analysis (EDA) of the interaction energy between various fragments. The aim of the work is to estimate the strength of the pi bonding and the aromatic character of the metallacyclic compounds. Analysis of the electronic structure shows that the metallacyclic moiety has five occupied pi orbitals, two with b1 symmetry and three with a2 symmetry, which describe the pi-bonding interactions. The metallabenzenes are thus 10 pi-electron systems. This holds for 16-electron and for 18-electron complexes. The pi bonding in the metallabenzenes results mainly from the b1 contribution, but the a2 contribution is not negligible. Comparison of the pi-bonding strength in the metallacyclic compounds with acylic reference molecules indicates that metallabenzenes should be considered as aromatic compounds whose extra stabilization due to aromatic conjugation is weaker than in benzene. The calculated aromatic stabilization energies (ASEs) are between 8.7 kcal mol(-1) for 13 and 37.6 kcal mol(-1) for 16 which is nearly as aromatic as benzene (ASE=42.5 kcal mol(-1)). The classical metallabenzene model compounds 1 and 4 exhibit intermediate aromaticity with ASE values of 33.4 and 17.6 kcal mol(-1). The greater stability of the 5d complexes compared with the 4d species appears not to be related to the strength of pi conjugation. From the data reported here there is no apparent trend or pattern which indicates a correlation between aromatic stabilization and particular ligands, metals, coordination numbers or charge. The lower metal-C5H5 binding energy of the 4d complexes correlates rather with weaker sigma-orbital interactions.     

McLafferty rearrangement of the radical cations of butanal and 3-fluorobutanal: a theoretical investigation of the concerted and stepwise mechanisms

The stepwise and concerted pathways for the McLafferty rearrangement of the radical cations of butanal (Bu(+)) and 3-fluorobutanal (3F-Bu(+)) are investigated with density functional theory (DFT) and ab initio methods in conjunction with the 6-311+G(d,p) basis set. A concerted transition structure (TS) for Bu(+), (H), is located with a Gibbs barrier height of 37.7 kcal/mol as computed with CCSD(T)//BHandHLYP. Three pathways for the stepwise rearrangement of Bu(+) have been located, which are all found to involve different complexes. The barrier height for the H(gamma) transfer is found to be 2.2 kcal/mol, while the two most favorable TSs for the C(alpha)-C(beta) cleavage are located 8.9 and 9.2 kcal/mol higher. The energies of the 3F-Bu(+) system have been calculated with the promising hybrid meta-GGA MPWKCIS1K functional of DFT. Interestingly, the fluorine substitution yields a barrier height of only 20.5 kcal/mol for the concerted TS, (3F-H). A smaller computed dipole moment, 12.1 D, for (3F-H) compared with 103.2 D for (H) might explain the stabilization of the substituted TS. The H(gamma) transfer, with a barrier height of 4.9 kcal/mol, is found to be rate-determining for the stepwise McLafferty rearrangement of 3F-Bu(+), in contrast to the unsubstituted case. By inspection of the spin and charge distributions of the stationary points, it is noted that the bond cleavages in the concerted rearrangements are mainly of heterolytic nature, while those in the stepwise channels are found to be homolytic.     

The ozonolysis of acetylene--a quantum chemical investigation.

The ozonolysis of acetylene was investigated using CCSD(T), CASPT2, and B3LYP-DFT in connection with a 6-311+G(2d,2p) basis set. The reaction is initiated by the formation of a van der Waals complex followed by a [4pi + 2pi] cycloaddition between ozone and acetylene (activation enthalpy DeltaH(a)(298) = 9.6 kcal/mol; experiment, 10.2 kcal/mol), yielding 1,2,3-trioxolene, which rapidly opens to alpha-ketocarbonyl oxide 5. Alternatively, an O atom can be transferred from ozone to acetylene (DeltaH(a)(298) = 15.6 kcal/mol), thus leading to formyl carbene, which can rearrange to oxirene or ketene. The key compound in the ozonolysis of acetylene is 5 because it is the starting point for the isomerization to the corresponding dioxirane 19 (DeltaH(a)(298) = 16.9 kcal/mol), for the cyclization to trioxabicyclo[2.1.0]pentane 10 (DeltaH(a)(298) = 19.5 kcal/mol), for the formation of hydroperoxy ketene 15 (DeltaH(a)(298) = 20.6 kcal/mol), and for the rearrangement to dioxetanone 9 (DeltaH(a)(298) = 23.6 kcal/mol). Compounds 19, 10, 15, and 9 rearrange or decompose with barriers between 13 and 16 kcal/mol to yield as major products formanhydride, glyoxal, formaldehyde, formic acid, and (to a minor extent) glyoxylic acid. Hence, the ozonolysis of acetylene possesses a very complicated reaction mechanism that deserves intensive experimental studies.     

Conformation-dependent state selectivity in O-O cleavage of ONOONO: an "inorganic cope rearrangement" helps explain the observed negative activation energy in the oxidation of nitric oxide by dioxygen

ONOONO has been proposed as an intermediate in the oxidation of nitric oxide by dioxygen to yield nitrogen dioxide. The O-O bond breaking reactions of this unusual peroxide, and subsequent rearrangements, were evaluated using CBS-QB3 and B3LYP/6-311G hybrid density functional theory. The three stable conformers (cis,cis-, cis,trans-, and trans,trans-ONOONO, based on the O-N-O-O dihedral angles of either approximately 0 degrees or approximately 180 degrees ) are predicted to have very different O-O cleavage barriers: 2.4, 13.0, and 29.8 kcal/mol, respectively. These large differences arise because bond breaking leads to correlation of the nascent NO(2) fragments with either the ground (2)A(1) state or the excited (2)B(2) state of NO(2), depending on the starting ONOONO conformation. A cis-oriented NO(2) fragment correlates with the (2)A(1) state, whereas a trans-oriented NO(2) fragment correlates with the (2)B(2) state. Each NO(2) fragment that correlates with (2)A(1) lowers the O-O homolysis energy by approximately 15 kcal/mol, similar to the approximately 17-25 kcal/mol (2)A(1) --> (2)B(2) energy difference in NO(2). Hence, this provides an unusual example of conformation-dependent electronic state selectivity. The O-O bond homolysis of cis,cis-ONOONO is particularly interesting because it has a very low barrier and arises from the most stable ONOONO conformer, and also due to obvious similarities to the well-known [3,3]-sigmatropic shift of 1,5-hexadiene, i.e., the Cope rearrangement. As an additional proof of our state selectivity postulate, a comparison is also made to breakage of the O-O bond of cis,cis-formyl peroxide, where no significant stabilization of the transition state is available because the (2)A(1) and (2)B(2) states of formyloxy radical are near-degenerate in energy. In the case of trans,trans-ONOONO, the O-O bond breaking transition state is a concerted rearrangement yielding O(2)NNO(2), whereas for cis,cis- and cis,trans-ONOONO, the initially formed NO(2) radical pairs can undergo further rearrangement to yield ONONO(2). It is proposed that previous spectroscopic observations of certain N=O stretching frequencies in argon-matrix-isolated products from the reaction of NO with O(2) (or (18)O(2)) are likely from ONONO(2), not the OONO radical as reported.     

Competition between alpha-cleavage and energy transfer in alpha-azidoacetophenones

Molecular modeling demonstrates that the first excited state of the triplet ketone (T1K) in azide 1b has a (pi,pi*) configuration with an energy that is 66 kcal/mol above its ground state and its second excited state (T2K) is 10 kcal/mol higher in energy and has a (n,pi*) configuration. In comparison, T1K and T2K of azide 1a are almost degenerate at 74 and 77 kcal/mol above the ground state with a (n,pi*) and (pi,pi*) configuration, respectively. Laser flash photolysis (308 nm) of azide 1b in methanol yields a transient absorption (lambdamax=450 nm) due to formation of T1K, which decays with a rate of 2.1 x 105 s-1 to form triplet alkylnitrene 2b (lambdamax=320 nm). The lifetime of nitrene 2b was measured to be 16 ms. In contrast, laser flash photolysis (308 nm) of azide 1a produced transient absorption spectra due to formation of nitrene 2a (lambdamax=320 nm) and benzoyl radical 3a (lambdamax=370 nm). The decay of 3a is 2 x 105 s-1 in methanol, whereas nitrene 2a decays with a rate of approximately 91 s-1. Thus, T1K (pi,pi*) in azide 1b leads to energy transfer to form nitrene 2b; however, alpha-cleavage is not observed since the energy of T2K (n,pi*) is 10 kcal/mol higher in energy than T1K, and therefore, T2K is not populated. In azide 1a both alpha-cleavage and energy transfer are observed from T1K (n,pi*) and T2K (pi,pi*), respectively, since these triplet states are almost degenerate. Photolysis of azide 1a yields mainly product 4, which must arise from recombination of benzoyl radicals 3a with nitrenes 2a. However, products studies for azide 1b also yield 4b as the major product, even though laser flash photolysis of azide 1b does not indicate formation of benzoyl radical 3b. Thus, we hypothesize that benzoyl radicals 3 can also be formed from nitrenes 2. More specifically, nitrene 2 does undergo alpha-photocleavage to form benzoyl radicals and iminyl radicals. The secondary photolysis of nitrenes 2 is further supported with molecular modeling and product studies.     

Reactivity of 13,13-dibromo-2,4,9,11-tetraoxadispiro[5.0.5.1]tridecane toward organolithiums: remarkable resistance to the DMS rearrangement

Reactions of dibromocyclopropane 2a, containing two spiro-fused 1,3-dioxane rings, with MeLi gave only the methylation products 8 and 9 even at elevated temperatures. In contrast, the cyclohexane analogue 2b treated with MeLi underwent a smooth rearrangement to bicyclo[1.1.0]butane 11b at -78, -10, or +35 degrees C. Treatment of 2a with PhLi gave the alpha-Ph anion 13 as the only product, which underwent smooth methylation with MeI to give 14. Under the same conditions, 2b with PhLi gave bicyclo[1.1.0]butane 11b accompanied by bromophenyl derivative 8b. Treatment of either dibromide with t-BuLi gave a mixture of products including debrominated cyclopropanes 12. Experimental results were augmented with DFT calculations for salts 23 and MP2//DFT-level calculations for carbenes 22. They demonstrated a higher stability of the dioxane alpha-bromo anion with respect to alpha-elimination by 4.8 kcal/mol and also a lower tendency of the carbene 22a to undergo rearrangement by 4.0 kcal/mol than the cyclohexane analogues. These differences have been attributed to the inductive effect of the four oxygen atoms, which results in lower LUMO energy, the higher positive charge at the carbenic center, and the overall more electrophilic character of carbene 22a as compared to the cyclohexane derivative 22b. The rearrangement of carbenes 22 to the corresponding allenes 1, the thermodynamic products, requires a higher activation energy DeltaG(double dagger)(298) by 4.2 kcal/mol for dioxane and 6.4 kcal/mol for cyclohexane derivatives than for the formation of the bicyclo[1.1.0]butanes 11. The DeltaG(double dagger)(298) for intramolecular insertions to the CH bond is low and calculated as 6.0 kcal/mol for dioxane 22a and 2.0 kcal/mol for the formation of cyclohexane 22b.     

A comparison of acetyl- and methoxycarbonylnitrenes by computational methods and a laser flash photolysis study of benzoylnitrene

Density functional theory (DFT), CCSD(T), and CBS-QB3 calculations were performed to understand the chemical and reactivity differences between acetylnitrene (CH(3)C(=O)N) and methoxycarbonylnitrene (CH(3)OC(=O)N) and related compounds. CBS-QB3 theory alone correctly predicts that acetylnitrene has a singlet ground state. We agree with previous studies that there is a substantial N-O interaction in singlet acetylnitrene and find a corresponding but weaker interaction in methoxycarbonylnitrene. Methoxycarbonylnitrene has a triplet ground state because the oxygen atom stabilizes the triplet state of the carbonyl nitrene more than the corresponding singlet state. The oxygen atom also stabilizes the transition state of the Curtius rearrangement and accelerates the isomerization of methoxycarbonylnitrene relative to acetylnitrene. Acetyl azide is calculated to decompose by concerted migration of the methyl group along with nitrogen extrusion; the free energy of activation for this concerted process is only 27 kcal/mol, and a free nitrene is not produced upon pyrolysis of acetyl azide. Methoxycarbonyl azide, on the other hand, does have a preference for stepwise Curtius rearrangement via the free nitrene. The bimolecular reactions of acetylnitrene and methoxycarbonylnitrene with propane, ethylene, and methanol were calculated and found to have enthalpic barriers that are near zero and free energy barriers that are controlled by entropy. These predictions were tested by laser flash photolysis studies of benzoyl azide. The absolute bimolecular reaction rate constants of benzoylnitrene were measured with the following substrates: acetonitrile (k = 3.4 x 10(5) M(-1) (s-1)), methanol (6.5 x 10(6) M(-1) s(-1)), water (4.0 x 10(6) M(-1) s(-1)), cyclohexane (1.8 x 10(5) M(-1) s(-1)), and several representative alkenes. The activation energy for the reaction of benzoylnitrene with 1-hexene is -0.06 +/- 0.001 kcal/mol. The activation energy for the decay of benzoylnitrene in pentane is -3.20 +/- 0.02 kcal/mol. The latter results indicate that the rates of reactions of benzoylnitrene are controlled by entropic factors in a manner reminiscent of singlet carbene processes.     

Electroreduction of a series of alkylcobalamins: mechanism of stepwise reductive cleavage of the Co-C bond

The electrochemical (EC) reduction mechanism of methylcobalamin (Me-Cbl) in a mixed DMF/MeOH solvent in 0.2 M tetrabutylammonium fluoroborate electrolyte was studied as a function of temperature and solvent ratio vs a nonaqueous Ag/AgCl/Cl(-) reference electrode. Double-potential-step chronoamperometry allowed the rate constant of the subsequent homogeneous reaction to be measured over the temperature range from 0 to -80 degrees C in 40:60 and 50:50 DMF:MeOH ratios. Activation enthalpies are 5.8 +/- 0.5 and 7.6 +/- 0.3 kcal/mol in the 40:60 and 50:50 mixtures of DMF/MeOH, respectively. Digital simulation and curve-fitting for an EC mechanism using a predetermined homogeneous rate constant of 5.5 x 10(3) s(-1) give E degrees' = -1.466 V, k degrees = 0.016 cm/s, and alpha = 0.77 at 20 degrees C for a quasi-reversible electrode process. Digital simulation of the results of Lexa and Savéant (J. Am. Chem. Soc. 1978, 100, 3220-3222) shows that the mechanism is a series of stepwise homogeneous equilibrium processes with an irreversible step following the initial electron transfer (ET) and allows estimation of the equilibrium and rate constants of these reactions. An electron coupling matrix element of H(kA) = (4.7 +/- 1.1) x 10(-4) eV ( approximately 46 J/mol) is calculated for the nonadiabatic ET step for reduction to the radical anion. A reversible bond dissociation enthalpy for homolytic cleavage of Me-Cbl is calculated as 31 +/- 2 kcal/mol. The voltammetry of the ethyl-, n-propyl-, n-butyl-, isobutyl-, and adenosyl-substituted cobalamin was studied, and estimated reversible redox potentials were correlated with Co-C bond distances as determined by DFT (B3LYP/ LANL2DZ) calculations.     

Free energy profiles for monomer capture in Grubbs- and SHOP-type olefin polymerization catalysts: a constraint ab initio molecular dynamics study

Density functional theory together with Car-Parrinello ab initio molecular dynamics simulation has been used to investigate the free energy profiles (FEP) of monomer capture in Grubbs- and SHOP-type olefin polymerization catalysts. The FEPs along the reaction coordinates at 300 K were determined directly by a point wise thermodynamic integration technique. Comparison between potential energy profile (PEP) and the FEP has been made. The results show that, for both catalysts, the PEP for the monomer ethylene uptake by the metal center is a typical Morse curve without energy barrier. However, a small barrier (1.8 kcal/mol for Grubbs catalyst and 2.4 kcal/mol for SHOP catalyst) exists on the FEP. The pi complexation energy on the FES at 300 K is higher by 10-12 kcal/mol over that on the PES. The differences between FES and PES are due to entropy contribution. Slow growth simulations on the ethylene capture process show that the ethylene attacks the metal center by an asynchronous mode. This indicates that the forming of the pi-bonding between the metal and ethylene is initiated by electrophilic attack of the metal to one of the ethylene carbons.     

DFT study of Co-C bond cleavage in the neutral and one-electron-reduced alkyl-cobalt(III) phthalocyanines

Density functional theory (DFT) has been applied to the analysis of the structural and electronic properties of the alkyl-cobalt(III) phthalocyanine complexes, [CoIIIPc]-R (Pc = phthalocyanine, R = Me or Et), and their pyridine adducts. The BP86/6-31G(d) level of theory shows good reliability for the optimized axial bond lengths and bond dissociation energies (BDEs). The mechanism of the reductive cleavage was probed for the [CoIIIPc]-Me complex which is known as a highly effective methyl group donor. In the present analysis, which follows a recent study on the reductive Co-C bond cleavage in methylcobalamin (J. Phys. Chem. B 2007, 111, 7638-7645), it is demonstrated that addition of an electron and formation of the pi-anion radical [CoIII(Pc*)]-Me- significantly lowers the energetic barrier required for homolytic Co-C bond dissociation. Such BDE lowering in [CoIII(Pc*)]-Me- arises from the involvement of two electronic states: upon electron addition, a quasi-degenerate pi*Pc state is initially formed, but when the cobalt-carbon bond is stretched, the unpaired electron moves to a sigma*Co-C state and the final cleavage involves the three-electron (sigma)2(sigma*)1 bond. As in corrin complexes, the pi*Pc-sigma*Co-C states crossing does not take place at the equilibrium geometry of [CoIII(Pc*)]-Me- but only when the Co-C bond is stretched to approximately 2.3 A. The DFT computed Co-C BDE of 23.3 kcal/mol in the one-electron-reduced phthalocyanine species, [CoIII(Pc*)]-Me-, is lowered by approximately 37% compared to the neutral Py-[CoIIIPc]-Me complex where BDE = 36.8 kcal/mol. A similar comparison for the corrin-containing complexes shows that a DFT computed BDE of 20.4 kcal/mol for [CoIII(corrin*)]-Me leads to approximately 45% bond strength reduction, in comparison to 37.0 kcal/mol for Im-[CoIII(corrin)]-Me+. These results suggest some preference by the alkylcorrinoids for the reductive cleavage mechanism.     

Experimental and theoretical studies on the thermal decomposition of heterocyclic nitrosimines

A series of substituted 2-nitrosiminobenzothiazolines (2) were synthesized by the nitrosation of the corresponding 2-iminobenzothiazolines (6). Thermal decomposition of 2a--f and of the seleno analogue 7 in methanol and of 3-methyl-2-nitrosobenzothiazoline (2a) in acetonitrile, 1,4-dioxane, and cyclohexane followed first-order kinetics. The activation parameters for thermal deazetization of 2a were measured in cyclohexane (Delta H(++) = 25.3 +/- 0.5 kcal/mol, Delta S(++) = 1.3 +/- 1.5 eu) and in methanol (Delta H(++) = 22.5 +/- 0.7 kcal/mol, Delta S(++) = -12.9 +/- 2.1 eu). These results indicate a unimolecular decomposition and are consistent with a proposed stepwise mechanism involving cyclization of the nitrosimine followed by loss of N(2). The ground-state conformations of the parent nitrosiminothiazoline (9a) and transition states for rotation around the exocyclic C==N bond, electrocyclic ring closure, and loss of N(2) were calculated using ab initio molecular orbital theory at the MP2/6-31G* level. The calculated gas-phase barrier height for the loss of N(2) from 9a (25.2 kcal/mol, MP4(SDQ, FC)/6-31G*//MP2/6-31G* + ZPE) compares favorably with the experimental barrier for 2a of 25.3 kcal/mol in cyclohexane. The potential energy surface is unusual; the rotational transition state 9a-rot-ts connects directly to the orthogonal transition state for ring-closure 9aTS. The decoupling of rotational and pseudopericyclic bond-forming transition states is contrasted with the single pericyclic transition state (15TS) for the electrocyclic ring-opening of oxetene (15) to acrolein (16). For comparison, the calculated homolytic strength of the N--NO bond is 40.0 kcal/mol (MP4(SDQ, FC)/6-31G*//MP2/6-31G* + ZPE).     

Photoinduced C-N bond cleavage in 2-azido-1,3-diphenyl-propan-1-one derivatives: photorelease of hydrazoic acid

Photolysis of 3-azido-1,3-diphenyl-propan-1-one (1a) in toluene yields 1,3-diphenyl-propen-1-one (2), whereas irradiation of 3-azido-2,2-dimethyl-1,3-diphenyl-propan-1-one (1b) results in the formation of mainly 2,2-dimethyl-1,3-diphenyl-propan-1-one. Laser flash photolysis (308 nm) of 1a,b in acetonitrile reveals a transient absorption (lambda max = approximately 310 nm) due to the formation of radicals 4a and 4b, respectively, which have lifetimes of approximately 14 micros at ambient temperature. TD-DFT calculations (B3LYP/6-31+G(d)) reveal that the first and second excited states of the triplet ketone (T1K (n,pi*) and T2K (pi,pi*)) in azide 1a are almost degenerate, at approximately 74 and 76 kcal/mol above the ground state (S0), respectively. Similarly, azide 1b has T1K and T2K 75 and 82 kcal/mol above S0, respectively. The calculated transition state for cleaving the C-N bond is located 71 and 74 kcal/mol above S0 in azides 1a and 1b, respectively. The calculated bond dissociation energies for breaking the C-N bond are 55 and 58 kcal/mol for azides 1a and 1b, respectively, making C-N bond breakage accessible from T1K in azides 1 at ambient temperature. In comparison, the irradiation of azides 1 in argon matrices at 14 K lead to the formation of the corresponding triplet alkyl nitrenes (1-n), via intramolecular energy transfer from T2K. The characterization of 1-n was supported by isotope labeling, IR spectroscopy, and molecular modeling.     

The Exceptional Conformational Bias of 1,5-Diene-3,4-diols Is Controlled by Electrostatic Interactions: A Theoretical Elucidation of the Origin in Rotational Isomer Stability

A conformational study of 1,5-hexadiene-3,4-diol (1) and 3-buten-1-ol (2) with the ab initio molecular orbital methods reveals that electrostatic interactions, rather than steric or tosional effects, control the conformational preferences of these compounds. It is found that the 1,2-dioxygen function prefers the anti arrangement due to the lone pair electron repulsion. A pseudo 1,3-diaxial oxygen/pi bond repulsion and a 1,3-diaxial attraction between an oxygen lone pair and a vinyl H are found to be responsible for the profound preference for the CO-eclipsed form in diol 1. Quantitatively, the repulsion between two gauche oxygen atoms is calculated to be approximately 1.5 kcal/mol at the MP2/6-31G//HF/6-31G level of theory. The repulsion between an oxygen atom and a pi bond at the pseudo 1,3-diaxial position is approximately 1 kcal/mol, and the attractive interaction between an oxygen atom and a vinyl H at the pseudo 1,3-diaxial position is approximately 0.5 kcal/mol, respectively, at the MP2/6-31G//MP2/6-31G level.     

Conformational selectivity in the diels-alder cycloaddition: predictions for reactions of s-trans-1,3-butadiene

Diels-Alder cycloaddition of s-trans-1,3-butadiene (1) should yield trans-cyclohexene (7), just as reaction of the s-cis conformer gives cis-cyclohexene (9). Investigation of this long-overlooked process with Hartree-Fock, Moller-Plesset, CASSCF, and DFT methods yielded in every case a C(2)-symmetric concerted transition state. At the B3LYP/6-31G (+ZPVE) level, this structure is predicted to be 42.6 kcal/mol above reactants, while the overall reaction is endothermic by 16.7 kcal/mol. A stepwise diradical process has been studied by UBLYP/6-31G theory and found to have barriers of 35.5 and 17.7 kcal/mol for the two steps. Spin correction lowers these values to 30.1 and 13.0 kcal/mol. The barrier to pi-bond rotation in cis-cyclohexene (9) is predicted (B3LYP theory) to be 62.4 kcal/mol, with trans-cyclohexene (7) lying 53.3 kcal/mol above cis isomer 9. Results suggest that pi-bond isomerization and concerted reaction may provide competitive routes for Diels-Alder cycloreversion. It is concluded that full understanding of the Diels-Alder reaction requires consideration of both conformers of 1,3-butadiene.     

Novel substituent effects on the mechanism of the thermal denitrogenation of 2,3-diazabicyclo[2.2.1]hept-2-ene derivatives, stepwise versus concerted

The substituent effect on the thermal denitrogenation mechanism of 7,7-disubstituted 2,3-diazabicyclo[2.2.1]hept-2-enes, concerted versus stepwise, has been investigated in detail. Unrestricted DFT calculations at the B3LYP/6-31G(d) level of theory suggest that azoalkanes that possess electron-withdrawing substituents at C(7) prefer to expel the nitrogen molecule in a stepwise manner. The activation energy is calculated to be ca. 36 kcal/mol for the dihydroxy-substituted azoalkane. In contrast, the preferred mechanism of the concerted denitrogenation is predicted for azoalkanes that possess electron-donating substituents at C(7). The activation energy is computed to be ca. 28 kcal/mol for the silyl-substituted azoalkane. The theoretical prediction of the substituent effects on the mechanistic change is supported by analyzing the activation parameters of the azoalkane decompositions. The activation enthalpy for the decomposition of the 7,7-diethoxy-substituted azoalkane is determined to be 39.1 kcal/mol, which is 13.1 kcal/mol higher in energy for the denitrogenation of the 7-silyl-substituted azoalkane. These dramatic substituent effects can be reasonably explained by the preferred electronic configuration of the lowest singlet state of the cyclopentane-1,3-diyls produced during the denitrogenation of the azoalkanes.     

Besides N(2), What Is the Most Stable Molecule Composed Only of Nitrogen Atoms?

Polynitrogen molecules have been studied systematically at high levels of ab initio and density functional theory (DFT). Besides N(2), the thermodynamically most stable N(n)() molecules, located with the help of a newly developed energy increment system, are all based on pentazole units. The geometric, energetic, and magnetic criteria establish pentazole (2) and its anion (3) to be as aromatic as their isoelectronic analogues, e.g., furan, pyrrole, and the cyclopentadienyl anion. The bond lengths in 2 and 3 are equalized; both have large aromatic stabilization energies (ASE) and also substantial magnetic susceptibility exaltations (Lambda). The C(s)() symmetric azidopentazole (14), a candidate for experimental investigation, is the lowest energy N(8) isomer but is still 196.7 kcal/mol higher in energy than four N(2) molecules. Octaazapentalene (12) with 10 pi electrons also is aromatic. The D(2)(d)() symmetric bispentazole (21) is the lowest energy N(10) minimum but is 260 kcal/mol higher in energy than five N(2) molecules. For strain-free molecules, the average deviation is +/-2.6 kcal/mol between the DFT energies and those based on the increment scheme. The increment scheme also provides estimates of the strain energies of polynitrogen compounds, e.g., tetraazatetrahedrane (8, 48.2 kcal/mol), octaazacubane (11, 192.6 kcal/mol), and N(20) (27, 294.6 kcal/mol), and is useful in searching for new high-energy-high-density materials.     

Chemisorption of (CHx and C2Hy) hydrocarbons on Pt(111) clusters and surfaces from DFT studies

We used the B3LYP flavor of density functional theory (DFT) to study the chemisorption of all CH(x) and C(2)H(y) intermediates on the Pt(111) surface. The surface was modeled with the 35 atom Pt(14.13.8) cluster, which was found to be reliable for describing all adsorption sites. We find that these hydrocarbons all bind covalently (sigma-bonds) to the surface, in agreement with the studies by Kua and Goddard on small Pt clusters. In nearly every case the structure of the adsorbed hydrocarbon achieves a saturated configuration in which each C is almost tetrahedral with the missing H atoms replaced by covalent bonds to the surface Pt atoms. Thus, (Pt(3))CH prefers a mu(3) hollow site (fcc), (Pt(2))CH(2) prefers a mu(2) bridge site, and PtCH(3) prefers mu(1) on-top sites. Vinyl leads to (Pt(2))CH-CH(2)(Pt), which prefers a mu(3) hollow site (fcc). The only exceptions to this model are ethynyl (CCH), which binds as (Pt(2))C=CH(Pt), retaining a CC pi-bond while binding at a mu(3) hollow site (fcc), and HCCH, which binds as (Pt)HC=CH(Pt), retaining a pi bond that coordinates to a third atom of a mu(3) hollow site (fcc) to form an off center structure. These structures are in good agreement with available experimental data. For all species we calculated heats of formation (DeltaH(f)) to be used for considering various reaction pathways on Pt(111). For conditions of low coverage, the most strongly bound CH(x) species is methylidyne (CH, BE = 146.61 kcal/mol), and ethylidyne (CCH(3), BE = 134.83 kcal/mol) among the C(2)H(y) molecules. We find that the net bond energy is nearly proportional to the number of C-Pt bonds (48.80 kcal/mol per bond) with the average bond energy decreasing slightly with the number of C ligands.     

Nature of the oxomolybdenum-thiolate pi-bond: implications for Mo-S bonding in sulfite oxidase and xanthine oxidase

The electronic structure of cis,trans-(L-N(2)S(2))MoO(X) (where L-N(2)S(2) = N,N'-dimethyl-N,N'-bis(2-mercaptophenyl)ethylenediamine and X = Cl, SCH(2)C(6)H(5), SC(6)H(4)-OCH(3), or SC(6)H(4)CF(3)) has been probed by electronic absorption, magnetic circular dichroism, and resonance Raman spectroscopies to determine the nature of oxomolybdenum-thiolate bonding in complexes possessing three equatorial sulfur ligands. One of the phenyl mercaptide sulfur donors of the tetradentate L-N(2)S(2) chelating ligand, denoted S(180), coordinates to molybdenum in the equatorial plane such that the OMo-S(180)-C(phenyl) dihedral angle is approximately 180 degrees, resulting in a highly covalent pi-bonding interaction between an S(180) p orbital and the molybdenum d(xy) orbital. This highly covalent bonding scheme is the origin of an intense low-energy S --> Mo d(xy) bonding-to-antibonding LMCT transition (E(max) approximately 16000 cm(-)(1), epsilon approximately 4000 M(-)(1) cm(-)(1)). Spectroscopically calibrated bonding calculations performed at the DFT level of theory reveal that S(180) contributes approximately 22% to the HOMO, which is predominantly a pi antibonding molecular orbital between Mo d(xy) and the S(180) p orbital oriented in the same plane. The second sulfur donor of the L-N(2)S(2) ligand is essentially nonbonding with Mo d(xy) due to an OMo-S-C(phenyl) dihedral angle of approximately 90 degrees. Because the formal Mo d(xy) orbital is the electroactive or redox orbital, these Mo d(xy)-S 3p interactions are important with respect to defining key covalency contributions to the reduction potential in monooxomolybdenum thiolates, including the one- and two-electron reduced forms of sulfite oxidase. Interestingly, the highly covalent Mo-S(180) pi bonding interaction observed in these complexes is analogous to the well-known Cu-S(Cys) pi bond in type 1 blue copper proteins, which display electronic absorption and resonance Raman spectra that are remarkably similar to these monooxomolybdenum thiolate complexes. Finally, the presence of a covalent Mo-S pi interaction oriented orthogonal to the MOO bond is discussed with respect to electron-transfer regeneration in sulfite oxidase and Mo=S(sulfido) bonding in xanthine oxidase.     

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.     

Revision of the dissociation energies of mercury chalcogenides--unusual types of mercury bonding.

Mercury chalcogenides HgE (E=O, S, Se, etc.) are described in the literature to possess rather stable bonds with bond dissociation energies between 53 and 30 kcal mol(-1), which is actually difficult to understand in view of the closed-shell electron configuration of the Hg atom in its ground state (...4f(14)5d(10)6s(2)). Based on relativistically corrected many body perturbation theory and coupled-cluster theory [IORAmm/MP4, Feenberg-scaled IORAmm/MP4, IORAmm/CCSD(T)] in connection with IORAmm/B3LYP theory and a [17s14p9d5f]/aug-cc-pVTZ basis set, it is shown that the covalent HgE bond is rather weak (2-7 kcal mol(-1)), the ground state of HgE is a triplet rather than a singlet state, and that the experimental bond dissociation energies have been obtained for dimers (or mixtures of monomers, dimers, and even trimers) Hg2E2 rather than true monomers. The dimers possess association energies of more than 100 kcal mol(-1) due to electrostatic forces between the monomer units. The covalent bond between Hg and E is in so far peculiar as it requires a charge transfer from Hg to E (depending on the electronegativity of E) for the creation of a single bond, which is supported by electrostatic forces. However, a bonding between Hg and E is reduced by strong lone pair-lone pair repulsion to a couple of kcal mol(-1). Since a triplet configuration possesses somewhat lower destabilizing lone pair energies, the triplet state is more stable. In the dimer, there is a Hg-Hg pi bond of bond order 0.66 without any a support. Weak covalent Hg-O interactions are supported by electrostatic bonding. The results for the mercury chalcogenides suggests that all experimental dissociation energies for group-12 chalcogenides have to be revised because of erroneous measurements.