Comparison of the performance of dispersion-corrected density functional theory for weak hydrogen bonds

Potential energy curves for five complexes with weak to medium strong hydrogen bonds have been computed with dispersion corrected DFT methods. The electronic density based vdW-DF2 and VV10 van der Waals density functionals have been tested, as well as an atom pair-wise correction method (DFT-D3). The short-range exchange-correlation components BLYP and rPW86-PBE together with the extended aug-cc-pVQZ basis sets have been employed. Reference data have been computed at the estimated CCSD(T)/CBS(aQ-a5) level of theory. The investigated systems are CH(4)·NH(3), Cl(3)CH·NH(3), NH(3)·NH(3), CH(3)F·C(2)H(2) and CH(3)F·H(2)O with binding energies ranging from -0.7 kcal mol(-1) to -5.5 kcal mol(-1). We find that all dispersion corrected methods perform reasonably well for these hydrogen bonds, but also observe distinct differences. The BLYP-D3 method provides the best results for three out of five systems. For the fluorinated complexes, the VV10 method gives remarkably good results. The vdW-DF2 method yields good interaction energies similar to the other methods (mean average deviation of 0.2-0.3 kcal mol(-1)), but fails to provide accurate equilibrium separations. Based on these results and previous experience with the computation of non-covalent interactions, for large-scale applications we can recommend DFT-D3 based structure optimizations with subsequent checking of interaction energies by single-point VV10 computations. Comparison of the DFT-D3 and VV10 results leads to the conclusion that the short-range exchange-correlation functional and not the dispersion correction mainly determines the achievable accuracy.     

Complex formation of trimethylaluminum and trimethylgallium with ammonia: evidence for a hydrogen-bonded adduct

We have investigated the formation of gas-phase adducts of trimethylaluminum and trimethylgallium with ammonia using room-temperature Fourier transform infrared experiments and density functional theory calculations. Our results indicate for the first time that, at higher partial pressures, a product distinct from the well-known (CH3)3M:NH3 adduct grows in for both M = Al and M = Ga. Comparison of the experimental and calculated IR spectra, along with calculations of the energetics, indicates that this second product is the result of hydrogen bonding of a second NH3 molecule to the (CH3)3M:NH3 adduct and can be written as (CH3)3M:NH3...NH3. The binding energy of this hydrogen-bonded adduct is calculated to be 26.8 kcal/mol for M = Al and 18.4 kcal/mol for M = Ga and is lower in energy (more stable) relative to the 1:1 (CH3)3M:NH3 adduct by 7.2 kcal/mol for M = Al and 6.6 kcal/mol for M = Ga. In contrast, an alternative complex involving the formation of two separate M-N donor-acceptor bonds, which is written as H3N:(CH3)3M:NH3, is calculated to be lower in energy relative to (CH3)3M:NH3 by only 0.1 kcal/mol for M = Al and 0.2 kcal/mol for M = Ga and is not observed experimentally. These results show that hydrogen bonding plays an important role in the interaction of ammonia with metal organic precursors involving Al, Ga, and In, under typical metal organic chemical vapor deposition AlGaInN growth conditions.     

Accurate thermochemical properties for energetic materials applications. I. Heats of formation of nitrogen-containing heterocycles and energetic precursor molecules from electronic structure theory

The heats of formation of 1H-imidazole, 1H-1,2,4-trizazole, 1H-tetrazole, CH3NO2, CH3N3, CH3NH2, CH2CHNO2, HClO4, and phenol, as well as cations and anions derived from some of the molecules have been calculated using ab initio molecular orbital theory. These molecules are important as models for compounds used for energetic materials synthesis. The predicted heats of formation of the heterocycle-based compounds are in excellent agreement with available experimental values and those derived from proton affinities and deprotonation enthalpies to <1 kcal/mol. The predicted value for the tetrazolium cation differs substantially from the experimental value, likely due to uncertainty in the measurement. The heats of formation of the nitro and amino molecules, as well as phenol/phenolate, also are in good agreement with the experimental values (<1.5 kcal/mol). The heat of formation of CH3N3 is predicted to be 72.8 kcal/mol at 298 K with an estimated error bar of +/-1 kcal/mol on the basis of the agreement between the calculated and experimental values for DeltaH(f)(HN3). The heat of formation at 298 K of HClO4 is -0.4 kcal/mol, in very good agreement with the experimental value, as well as a W2 literature study. An extrapolation of the CCSD(T)/aug-cc-pV(Q,5) energies was required to obtain this agreement. This result suggests that very large basis sets (> or =aug-cc-pV5Z) may be needed to fully recover the valence correlation energy contribution in compounds containing elements with high formal oxidation states at the central atom. In addition tight d functions are needed for the geometry predictions. Douglas-Kroll-Hess (DKH) scalar relativistic corrections for HClO4 and ClO4- at the MP2 level with correlation-consistent DKH basis sets were predicted to be large, likely due to the high formal oxidation state at the Cl.     

H-Atom abstraction from CH(3)NHNH(2) by NO(2): CCSD(T)/6-311++G(3df,2p)//MPWB1K/6-31+G(d,p) and CCSD(T)/6-311+G(2df,p)//CCSD/6-31+G(d,p) calculations

Stationary points of paths for H atom abstraction from CH(3)NHNH(2) (monomethylhydrazine) by NO(2) were characterized via CCSD(T)/6-311++G(3df,2p)//MPWB1K/6-31+G(d,p) and CCSD(T)/6-311+G(2df,p)//CCSD/6-31+G(d,p) calculations. Five transition states connecting CH(3)NHNH(2)-NO(2) complexes to a manifold that includes CH(3)NHNH-HONO, CH(3)NNH(2)-HONO, CH(3)NNH(2)-HNO(2), and CH(3)NHNH-HNO(2) complexes were identified. Transition states that connect CH(3)NHNH-HONO, CH(3)NNH(2)-HONO, CH(3)NNH(2)-HNO(2), and CH(3)NHNH-HNO(2) complexes to each other via H atom exchange and/or hindered internal rotation were also identified. The high point in the minimum energy path from the CH(3)NHNH(2) + NO(2) reactant asymptote to the manifold of HONO-containing product states is a transition state 8.6 kcal/mol above the reactant asymptote. From a kinetics standpoint, this value is considerably higher than the 5.9 kcal/mol value that was estimated for it based on theoretical results for H atom abstraction from NH(3) by NO(2).     

Favorable pendant-amino metal chelation in VX nerve agent model systems

We have performed DFT computational studies [B3LYP, 6-31+G] to obtain metal ion coordination isomers of VX-Me [MeP(O)(OMe)(SCH2CH2NMe2)], a model of two of the most lethal nerve agents: VX [MeP(O)(OEt)(SCH2CH2N(iPr)2)] and Russian-VX [MeP(O)(OCH2CHMe2)(SCH2CH2N(Et)2)]. Our calculations involved geometry optimizations of the neutral VX-Me model as well as complexes with H+, Li+, Na+, K+, Be2+, Mg2+, and Ca2+ that yielded 2-8 different stable chelation modes for each ion that involved mainly mono- and bidentate binding. Importantly, our studies revealed that the [O(P),N] bidentate binding mode, long thought to be the active mode in differentiating the hydrolytic path of VX from other nerve agents, was the most stable for all ions studied here. Binding energy depended mainly on ionic size as well as charge, with binding energies ranging from 364 kcal mol(-1) for Be2+ to 33 kcal mol(-1) for K+. Furthermore, calculated NMR shifts for VX-Me correlate to experimental values of VX.     

The cysteine radical cation: structures and fragmentation pathways

A theoretical study on the structures, relative energies, isomerization reactions and fragmentation pathways of the cysteine radical cation, [NH(2)CH(CH(2)SH)COOH].+, is reported. Hybrid density functional theory (B3LYP) has been used in conjunction with the 6-311++G(d,p) basis set. The isomer at the global minimum, Captodative-1, has the structure NH(2)C.(CH(2)SH)C(OH)(2)+; the stability of this ion is attributed to the captodative effect in which the NH(2) functions as a powerful pi-electron donor and C(OH)(2)+ as a powerful pi-electron acceptor. Ion Distonic-S-1, H(3)N(+)CH(CH(2)S.)COOH, in which the radical is formally situated on the S atom, is higher in enthalpy (DeltaH degrees (0)) than Captodative-1 by 6.1 kcal mol(-1), but is lower in enthalpy than another isomer Distonic-C-1, H(3)N(+)C.(CH(2)SH)COOH, by 8.2 kcal mol(-1). Isomerization of the canonical radical cation of cysteine, [H(2)NCH(CH(2)SH)COOH].+, (Canonical-1), to Captodative-1 has an enthalpy of activation of 25.8 kcal mol(-1), while the barrier against isomerization of Canonical-1 to Distonic-S-1 is only 9.6 kcal mol(-1). Two additional transient tautomers, one with the radical located at C(alpha) and the charge on SH(2), and the other a carboxy radical with the charge on NH(3), are reported. Plausible fragmentation pathways (losses of small molecules, CO(2), CH(2)S, H(2)S and NH(3), and neutral radicals COOH. , HSCH(2). and NH(2).) from Canonical-1 are examined.     

EPR and electronic spectral studies of bis(dihalcogenocarbamato)copper(II)-solvent interactions

The interaction of bis(diethyldithiocarbamato)copper(II), Cu(Et2dtc)2, and bis(diethyldiselenocarbamato)copper(II), Cu(Et2dsc)2, complexes with solvents is studied by EPR and electronic spectroscopy. The solvents used are CCl4, CHCl3, CH2Cl2, C6H5 x CH3, DMFA and DMSO. It is found that Cu(Et2dsc)2 is destroyed in a first order reaction in CCl4 with an activation energy of 5.2 kcal/mol. The other complex, Cu(Et2dtc)2, is only destroyed in DMSO. The observed effects and reaction pathways are discussed in terms of solute-solvent donor-acceptor interactions taking into account the differences in the electronic structures of both complexes.     

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.     

A 3-fold "butterfly valve" in command of the encapsulation's kinetic stability. Molecular baskets at work

Molecular basket 1, composed of a semirigid tris-norbornadiene framework and three revolving pyridine-based gates at the rim, has been built to "dynamically" enclose space and as such regulate molecular encapsulation. The gates were shown to fold via intramolecular hydrogen bonding and thereby form a C3nu symmetrical receptor: the 1H NMR resonance for the amide N-H protons of the pyridine gates appeared downfield (delta= 10.98 ppm), and the N-H vibrational stretch (IR) was observed at 3176 cm(-1). Accordingly, density functional theory (DFT, B3LYP) investigations revealed for the closed conformers of 1 to be energetically the most stable and dominant. The gearing of the pyridine "gates", about their axis, led to the interconversion of two dynamic enantiomers 1A and 1B comprising the clockwise and counterclockwise seam of intramolecular hydrogen bonds. Dynamic 1H NMR spectroscopic measurements and line-shape simulations suggested that the energy barrier of 10.0 kcal/mol (DeltaG++(A/B), 298 K) is required for the 1A/B interconversion, when CCl4 occupies the cavity of 1. Likewise, the activation free energy for CCl4 departing the basket was found to be 13.1 kcal/mol (DeltaG++, 298 K), whereas the thermodynamic stability of 1:CCl4 complex was -2.7 kcal/mol (DeltaGdegrees, 298 K). In view of that, CCl4 (but also (CH3)3CBr) was proposed to escape from, and a molecule of solvent to enter, the basket when the gates rotate about their axis: the exit of CCl4 requires the activation energy of 12.7 kcal/mol (DeltaG++(A/B) + DeltaGdegrees), similar to the experimentally found 13.1 kcal/mol (DeltaG++).     

13C NMR, infrared, solvation and theoretical investigation of the conformational isomerism in 1-haloacetones (X = Cl, Br and I)

The solvent dependence of the 13C NMR spectra of chloroacetone (CA), bromoacetone (BA) and iodoacetone (IA) are reported and the 3J(CH) couplings analysed using ab initio calculations and solvation theory. In CA the energy difference (E(cis) - E(gauche)) between the cis (Cl-C-C=O 0 degrees) and gauche (Cl-C-C=O 155 degrees) conformers is 1.7 kcal mol(-1) in the vapour, decreasing to 0.8 kcal mol(-1) in CCl4 solution and to -1.0 kcal mol(-1) in the pure liquid. The conformational equilibrium, in BA, is between the more polar cis (Br-C-C=O 0 degrees) and gauche (Br-C-C=O 132 degrees) conformations. The energy difference (E(cis) - E(gauche)) is 1.8 kcal mol(-1) in the vapour, decreasing to 0.9 kcal mol(-1) in CCl4 solution and to -0.4 kcal mol(-1) in the pure liquid. The energy difference (E(cis) - E(gauche)), in IA, between the cis (I-C-C=O 0 degrees) and gauche (I-C-C=O 104 degrees) conformers is 1.1 kcal mol(-1) in the vapour phase, decreasing to 0.5 kcal mol(-1) in CCl4 solution and to -0.5 kcal mol(-1) in the pure liquid. The vapour state energy difference for BA [1.4 kcal mol(-1) at B3LYP/6-311++G(d,p)] and for IA [1.6 kcal mol(-1) at B3LYP/6-311++G(d,p)/LANL2DZ)] are in very good agreement with the above values. For CA the agreement is also satisfactory [1.4 kcal mol(-1) at B3LYP/6-311++G(d,p)].     

Free energies for degradation reactions of 1,2,3-trichloropropane from ab initio electronic structure theory

Electronic structure methods were used to calculate the gas and aqueous phase reaction energies for reductive dechlorination (i.e., hydrogenolysis), reductive β-elimination, dehydrochlorination, and nucleophilic substitution by OH? of 1,2,3-trichloropropane. The thermochemical properties ΔH(f)°(298.15 K), S°(298.15 K, 1 bar), and ΔG(S)(298.15 K, 1 bar) were calculated by using ab initio electronic structure calculations, isodesmic reactions schemes, gas-phase entropy estimates, and continuum solvation models for 1,2,3-trichloropropane and several likely degradation products: CH3?CHCl?CH2Cl, CH2Cl?CH2?CH2Cl, C?H2?CHCl?CH2Cl, CH2Cl?C?H?CH2Cl, CH2═CCl?CH2Cl, cis-CHCl═CH?CH2Cl, trans-CHCl═CH?CH2Cl, CH2═CH?CH2Cl, CH2Cl?CHCl?CH2OH, CH2Cl?CHOH?CH2Cl, CH2═CCl?CH2OH, CH2═COH?CH2Cl, cis-CHOH═CH?CH2Cl, trans-CHOH═CH?CH2Cl, CH(═O)?CH2?CH2Cl, and CH3?C(═O)?CH2Cl. On the basis of these thermochemical estimates, together with a Fe(II)/Fe(III) chemical equilibrium model for natural reducing environments, all of the reactions studied were predicted to be very favorable in the standard state and under a wide range of pH conditions. The most favorable reaction was reductive β-elimination (ΔG(rxn)° ≈ ?32 kcal/mol), followed closely by reductive dechlorination (ΔG(rxn)° ≈ ?27 kcal/mol), dehydrochlorination (ΔG(rxn)° ≈ ?27 kcal/mol), and nucleophilic substitution by OH? (ΔG(rxn)° ≈ ?25 kcal/mol). For both reduction reactions studied, it was found that the first electron-transfer step, yielding the intermediate C?H2?CHCl?CH2Cl and the CH2Cl?C?H?CH2Cl species, was not favorable in the standard state (ΔG(rxn)° ≈ +15 kcal/mol) and was predicted to occur only at relatively high pH values. This result suggests that reduction by natural attenuation is unlikely.     

A computational investigation of the reactions of methylene, chlorocarbene, and dichlorocarbene with cyclopropane

The reactions of CH(2), CHCl, and CCl(2) with cyclopropane, 1, have been examined computationally. In all cases the lowest energy reaction between the carbene and 1 is predicted to be C-H insertion. In the reaction of CH(2) with 1, the transition state for C-C insertion leading to cyclobutane is 1.7 kcal/mol higher in enthalpy than the transition state for C-H insertion at the G3B3 level. A pathway higher in energy than C-H insertion in the reactions of CHCl and CCl(2) with 1 involves two-bond cleavages generating ethylene along with chloro and dichloroethylene, respectively.     

Dimesitylketone O-oxide: spectroscopic characterization, conformation, and reaction modes: OH formation and OH capture.

Dimesitylketone O-oxide 1b was synthesized by photolysis of dimesityldiazomethane dissolved in an oxygen saturated CCl3F solution at 140 K. Conformation and geometry of 1b were determined by comparing measured NMR chemical shifts with the corresponding chemical shifts calculated at the DFT-IGLO level of theory where it had to be considered that the molecule exists in two enantiomeric forms. Measured and calculated 1H chemical shifts agree within 0.1 ppm while the calculated 13C shift of the COO carbon (210.6 ppm) differs by only 0.4 ppm from the measured shift of 211.0 ppm. The two mesityl rings are perpendicular to each other and enclose angles of 40 and 57 degrees with the COO plane. The preferred rearrangement process of 1b is an H migration from one of the ortho-methyl groups to the terminal O atom of the COO unit. The calculated activation enthalpy of this process is 12.7 kcal/mol (B3LYP/cc-pVTZ). In contrast, the activation enthalpy for isomerization to dioxirane is 5 kcal/mol higher. In CCl3F, the activation barrier for the thermal decay was determined to be 13.8 +/- 0.2 kcal/mol and in acetonitrile 13.1 +/- 0.4 kcal/mol. H migration initiates cleavage of the OO bond and the production of an OH and a benzyl radical. Recombination of the latter in the solvent cage leads to the formation of 2-methylhydroxy-pentamethylbenzophenone, while escape of the OH radical from the solvent cage yields a ketone. These results confirm the possibility of OH production from carbonyl oxides in the solution phase.     

Dimerization and double proton transfer-induced tautomerism of 4(3H)-pyrimidinone in solution studied by IR spectroscopy and quantum chemical calculations

The tautomerism and dimerization of 4(3H)-pyrimidinone (4(3H)Pyr) in carbon tetrachloride (CCl4) and chloroform (CHCl3) solutions were investigated using IR spectroscopy and quantum chemical calculations. The observed IR spectra in the NH and OH stretching regions clearly revealed the predominance of the keto tautomer in both solvent systems. The enol form only exists in a very small proportion in the CCl4 solution. The tautomeric constant for the two monomers KT[OH/NH] = 0.012 and DeltaE = 2.62 kcal/mol were estimated at 25 degrees C. This result was supported by the self-consistent reaction field/polarizable continuum (SCRF/PCM) calculation at the MP4(full, SDQ)/aug-cc-pVDZ level, which predicted DeltaE = 3.06 kcal/mol in CCl4. In the C=O stretching region, two bands were observed, suggesting the coexistence of two keto structures at equilibrium. The calculated IR spectra indicated that the bands at 1711 and 1675 cm(-1) arise from the keto monomer and keto-keto (KK) ring dimer, respectively. At elevated temperature, the populations of both the keto and enol monomers increased for the CCl4 solution. The present study revealed that the keto <--> enol tautomerization does not occur in the isolated monomer molecule. The double proton transfer (DPT) reaction in the KK ring dimer presumably plays a substantial role in the population increase of the enol monomer. To our knowledge, this may be the first observation of the tautomerization in a model base pair via the temperature-induced ground-state DPT reaction under a nonpolar liquid environment reported so far. This tautomerism can serve as a mimic circumstance for the spontaneous mutations induced by proton transfer in the DNA base pairs.     

Theoretical study of (CH...C)- hydrogen bonds in CH(4-n)X(n) (X = F, Cl; n = 0, 1, 2) systems complexed with their homoconjugate and heteroconjugate carbanions

This work deals with a theoretical study of the (CH...C)- hydrogen bonds in CH4, CH3X, and CH2X2 (X = F, Cl) complexed with their homoconjugate and heteroconjugate carbanions. The properties of the complexes are calculated with the B3LYP method using the 6-311++G(d,p) or 6-311++G(2df,2p) basis sets. The deprotonation enthalpies (DPE) of the CH bond or the proton affinities of the carbanions (PA(C-) are calculated as well. All the systems with the exception of the CH4...CHCl2(-) one are characterized by a double minimum potential. In some of the complexes, the (CH(b)...C)- hydrogen bond is linear. In other systems, such as CH3F...CH2F- and CH3F...CHF2(-), there is a large departure from linearity, the systems being stabilized by electrostatic interactions between the nonbonded H of the neutral molecule and the F atom of the carbanion. In the transition state, the (CH(b)...C)- bond is linear, and there is a large contraction of the intermolecular C...C distance. The binding energies vary within a large range, from -1.4 to -11.1 kcal mol(-1) for the stable complexes and -8.6 to -44.1 kcal mol(-1) for the metastable complexes. The energy barriers to proton transfer are between 5 and 20 kcal mol(-1) for the heteroconjugate systems and between 3.8 and 8.3 kcal mol(-1) for the homoconjugate systems. The binding energies of the linear complexes depend exponentially on 1.5DPE - PA(C-), showing that the proton donor is more important than the proton acceptor in determining hydrogen bond strength. The NBO analysis indicates an important electronic reorganization in the two partners. The elongations of the CH bond resulting from the interaction with the carbanion depend on the occupation of the sigma*(CH(b)) antibonding orbitals and on the hybridization of the C bonded to H(b). The frequency shifts of the nu(CH)(A1) stretching vibration range between 15 and 1150 cm(-1). They are linearly correlated to the elongation of the CH(b) bond.     

Computational study on mechanistic details of the aminoethanol rearrangement catalyzed by the vitamin B12-dependent ethanolamine ammonia lyase: His and Asp/Glu acting simultaneously as catalytic auxiliaries

The rearrangement of aminoethanol catalyzed by ethanolamine ammonia lyase is investigated by computational means employing DFT (B3LYP/6-31G) and ab initio molecular orbital theory (QCISD/cc-pVDZ). The study aims at providing a detailed account on various crucial aspects, in particular a distinction between a direct intramolecular migration of the partially protonated NH(2) group vs elimination of NH(4)(+). Three mechanistic scenarios were explored: (i) According to the calculations, irrespective of the nature of the protonating species, intramolecular migration of the NH(3) group is energetically less demanding than elimination of NH(4)(+). However, all computed activation enthalpies exceed the experimentally derived activation enthalpy (15 kcal/mol) associated with the rate-determining step, i.e., the hydrogen abstraction from the 5'-deoxyadenosine by the product radical. For example, when imidazole is used as a model system for His interacting with the NH(3) group of the substrate, the activation enthalpy for the migration process amounts to 27.4 kcal/mol. If acetic acid is employed to mimic Asp or Glu, the activation enthalpy is somewhat lower, being equal to 24.2 kcal/mol. (ii) For a partial deprotonation of the substrate 2 at the OH group, the rearrangement mechanism consists of the dissociation of an NH(2) radical from C(2) and its association at C(1) atom. For all investigated proton acceptors (i.e., OH(-), HCOO(-), CH(3)COO(-), CH(2)NH, imidazole), the activation enthalpy for the dissociation step also exceeds 15 kcal/mol. Typical data are 20.2 kcal/mol for Ac(-) and 23.8 kcal/mol for imidazole. (iii) However, in a synergistic action of partial protonation of the NH(2) group and partial deprotonation of the OH group by the two conceivable catalytic auxiliaries Asp/Glu and His, the activation enthalpy computed is compatible with the experimental data. For imidazole and acetate as model systems, the activation enthalpy is equal to 13.7 kcal/mol. This synergistic action of the two catalytic groups is expected to take place in a physiologically realistic pH range of 6-9.5, and the present computational findings may help to further characterize the yet unknown structural details of the ethanolamine ammonia lyase's active site.     

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

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

Ab Initio Kinetics for Thermal Decomposition of CH(3)N(•)NH(2), cis-CH(3)NHN(•)H, trans-CH(3)NHN(•)H, and C(•)H(2)NNH(2) Radicals

The thermal decomposition of the CH(3)N(?)NH(2), cis-CH(3)NHN(?)H, trans-CH(3)NHN(?)H, and C(?)H(2)NNH(2) radicals, which are the four radical products from the H-abstraction reactions of monomethylhydrazine, were theoretically studied by using ab initio Rice-Ramsperger-Kassel-Marcus (RRKM) transition-state theory and master equation analysis. Various decomposition pathways were identified by using either the QCISD(T)/cc-pV∞Z//CASPT2/aug-cc-pVTZ or the QCISD(T)/cc-pV∞Z//B3LYP/6-311++G(d,p) quantum chemistry methods. The results reveal that the β-scission of NH(2) to form methyleneimine is the predominant channel for the decomposition of the C(?)H(2)NNH(2) radical due to its small energy barrier of 13.8 kcal mol(-1). The high pressure limit rate coefficient for the reaction is fitted by 3.88 × 10(19)T(-1.672) exp(-9665.13/T) s(-1). In addition, the pressure dependent rate coefficients exhibit slight temperature dependence at temperatures of 1000-2500 K. The cis-CH(3)NHN(?)H and trans-CH(3)NHN(?)H radicals are the two distinct spatial isomers with an energy barrier of 26 kcal mol(-1) for their isomerization. The β-scission of CH(3) from the cis-CH(3)NHN(?)H radical to form trans-diazene has an energy barrier of 35.2 kcal mol(-1), and the β-scission of CH(3) from the trans-CH(3)NHN(?)H radical to form cis-diazene has an energy barrier of 39.8 kcal mol(-1). The CH(3)N(?)NH(2) radical undergoes the β-scission of methyl hydrogen and amine hydrogen to form CH(2)═NNH(2), trans-CH(3)N═NH, and cis-CH(3)N═NH products, with the energy barriers of 42.8, 46.0, and 50.2 kcal mol(-1), respectively. The dissociation and isomerization rate coefficients for the reactions were calculated via the E/J resolved RRKM theory and multiple-well master equation analysis at temperatures of 300-2500 K and pressures of 0.01-100 atm. The calculated rate coefficients associated with updated thermochemical property data are essential components in the development of kinetic mechanisms for the pyrolysis and oxidation of MMH and its derivatives.     

Effect of geminal substitution at silicon on 1-sila- and 1,3-disilacyclobutanes' strain energies, their 2+2 cycloreversion enthalpies, and Si=C pi-bond energies in silenes.

To evaluate the effect of geminal substitution at silicon on 1-sila- and 1,3-disilacyclobutanes' strain energies, their 2+2 cycloreversion enthalpies, and Si=C pi-bond energies in silenes, an ab initio MO study of silenes, R2Si=CH2 (1), 1-silacyclobutanes, cyclo-R2Si(CH2)3 (2), and 1,3-disilacyclobutanes, cyclo-(R2SiCH2)2 (3), was performed using the level of theory denoted MP4/TZ(d)//MP2/6-31G(d) (TZ means the 6-311G(d) basis set for elements of the second period and hydrogen, and the McLean-Chandler (12s,9p)/[6s,5p](d) basis set for the third period elements). In the series R = H, CH3, SiH3, CH3O, NH2, Cl, F, the growth of the reaction enthalpies and strain energies is proportional to the substituents' electronegativities. 2+2 cycloreversion of 2 is endothermic by 40.6-63.1 kcal/mol, whereas that of 3 is endothermic by 72.7-114.2 kcal/mol. On going from a silicon to a fluorine substituent at the sp2-hybridized silicon atom, the pi-bond energy in 1 weakens by 11.3 kcal/mol, and the Si=C bond length shortens by 0.053 A. The effect of substituents' electronegativities at the double-bonded silicon atom in silenes is formulated as follows: the higher electronegativity, the shorter and the weaker the Si=C pi-bond. The latter is rationalized in terms of more strained geometry resulting from the energetic cost for planarizing the R2SiC moiety. The enthalpies of the ring-opening reaction are 68.0-80.1 kcal/mol (a cleavage of the Si-C bond in 3), 65.0-76.4 kcal/mol (a cleavage of the Si-C bond in 2), and 58.0-64.9 kcal/mol (a cleavage of the C-C bond in 2). The pronounced difference in the enthalpies of 2+2 cycloreversion of 1-sila- and 1,3-disilacyclobutanes is mainly due to the difference in the enthalpies of diradicals' decomposition. The decomposition of diradicals resulting from a cleavage of C-C and Si-C bonds in 2 is exothermic by 24.3-3.3 kcal/mol (apart from the difluoro derivative which is endothermic by 5.1 kcal/mol) and 27.0-13.3 kcal/mol, respectively. The decomposition of a 1,4-diradical resulting from ring opening of 3, apart from the disilyl derivative, is the endothermic process for which the enthalpy varies from 10.6 to 40.4 kcal/mol.     

Theoretical studies on alkyne addition to molybdenum alkylidenes: the preference of alpha- and beta-additions.

The addition of alkynes HC=CR to Mo(NH)(CH(2))(OR')(2) (R = H, Me, Ph; R' = CH(3), CF(3)) has been studied with both ab initio molecular orbital and density functional calculations. Geometry optimizations were carried out with the HF/3-21G, HF/HW3, and B3LYP/HW3 methods. The transition structures for these addition reactions are in distorted trigonal bipyramidal geometries, similar to those of alkene additions. The calculated activation enthalpy for HC=CH addition to Mo(NH)(CH(2))(OR')(2) is about 10.3 kcal/mol for R' = CH(3) and about 2.3 kcal/mol for R' = CF(3), indicating a significant preference for acetylene addition to Mo(NH)(CH(2))(OCF(3))(2) over Mo(NH)(CH(2))(OCH(3))(2). These barriers are higher than those of the corresponding ethylene addition by about 2-4 kcal/mol, even though the reaction of acetylene is much more exothermic. The alpha-addition of HC=CR (R = Me, Ph) is found to be considerably more favorable than the beta-addition to Mo(NH)(CH(2))(OR')(2). Interestingly, the alpha-addition has a lower activation energy, while the beta-addition has a higher activation energy, compared to that of the parent acetylene addition. Thus, alpha-addition is intrinsically favored over beta-addition by over 4 kcal/mol. This preference is reduced by solvent effect. All these can be explained by a destabilizing interaction between the nonreacting pi-orbital of alkyne and one of the lone pairs on the imido nitrogen. The steric effect of the bulky ligands in the real catalysts is also investigated qualitatively by the PM3 method. These studies give results in good accord with the experimentally observed regioselectivity.