A theoretical study on C2HXSi silylenes (X = H,CN, NH2, and OMe)

Relative stabilities and structural characters of 30 silylenic C₂HXSi species (X = H, NH₂, CN, and OMe), with singlet (s) and/or triplet (t) states, are calculated at six levels of theory: HF/6‐311++G**, MP3/6‐31G*, B1LYP/6‐311++G**, B3LYP/6‐311++G**, MP2/6‐311++G**, and MP4(SDTQ)/6‐311++G**. The four possible isomers considered for C₂SiHX are (i) 3‐X‐1‐silacyclopropenylidene ( 1 _s‐X and 1 _t‐X ), (ii) X‐vinilydensilylene ( 2 _s‐X and 2 _t‐X ), (iii) ethynyl‐X‐silylene ( 3 _s‐X and 3 _t‐X ), and (iv) (X‐ethynyl)silylene ( 4 _s‐X and 4 _t‐X ). The GIAO–NICS calculations show that singlet cyclic structures, 1 _s‐X , are considerably more aromatic than benzene. Conversely, triplet cyclic C₂HCNSi breaks down through optimization, and transforms into a novel high‐spin acyclic carbenosilylene minimum ( 1 _t‐CN ). Singlet 3 and triplet 3 cross at a divalent angle (|XSiC) of 152°. This angle narrows to 137° for crossing of singlet 3 _s‐CN and triplet, 3 _t‐CN . The smallest |XSiC occurs at 132° for crossing of 3 _s‐H and 3 _t‐H . © 2007 Wiley Periodicals, Inc. Heteroatom Chem 18:283–293, 2007; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/hc.20297     

Ab initio and DFT energetics of silylenic X–CNSi (X=H, F, Cl, and Br)

Singlet–triplet energy splitting for 24 silylenic reactive intermediates, X–CNSi (where X=H, F, Cl and Br), are compared and contrasted at 11 levels of theory: B1LYP/6-31++G**, B3LYP/6-31++G**, B1LYP/6-311++G**, B3LYP/6-311++G**, MP3/6-31G*, MP3/6-311++G**, MP2/6-31+G**, MP2/6-311++G**, MP4 (SDTQ)/6-311++G**, QCISD(T)/6-311++G** and CCSD(T)/6-311++G**. Each X-substituted silylenic species may either be singlet (s) or triplet (t), with one of the following three structures: 3-X-2-aza-1-silacyclopropenylidene (1_s-X, 1_t-X); [(X-imino)methylene]silylene (2_s-X, 2_t-X); and X-cyanosilylene (3_s-X, 3_t-X). For all X–CNSi species studied, orders of singlet–triplet energy separations (ΔE_s-t,X), appear as a function of electro-negativity (F>Cl>Br>H). For the six H–CNSi isomers (X=H), stability order is: 3_s-H>1_s-H>2_t-H>3_t-H>2_s-H>1_t-H. Likewise, stability order for the six isomers with X=F, is: 3_s-F>3_t-F>1_s-F>1_t-F>2_s-F>2_t-F. For X=Cl, the order of stability is: 3_s-Cl>1_s-Cl>3_t-Cl>2_t-Cl>1_t-Cl>2_t-Cl. Finally, the order of stability for six isomers of Br–CNSi is: 3_s-Br>3_t-Br>1_s-Br>2_s-Br>2_t-Br>1_t-Br. The lowest energy minimum, among all 24 species scrutinized, appears to be the singlet acyclic 3_s-X. Triplet silylene 2_t-H is suggested to be more stable than its corresponding 2_s-H at MP3, MP2 and DFT levels of theory. Comparisons between relative stabilities; multiplicities and geometrical parameters of 1–3 are discussed.     

Theoretical study of light-metal tetrahydroborates

Ab initio calculations of the structural, energetic, vibrational, and magnetic characteristics of the lowest-lying structures for isolated molecules and ions of light-metal tetrahydroborates (Li, Na, Be, Mg, and Al) have been performed by the perturbation theory (MP2), quadratic configuration interaction (QCISD(T)), coupled cluster (CCSD(T)), and density functional theory (B3LYP) methods using the 6-31G*, 6-31G**, 6-311+G**, and 6-311++G** basis sets. The trends in the behavior of the structural parameters, the energies of different decomposition pathways, barriers to internal rotation of BH₄ groups, normal mode frequencies, magnetic shielding constants, and spin density distribution (in radicals) have been analyzed in various related series of these compounds. The results obtained by ab initio methods and at the DFT level are compared. The economical approximation B3LYP/6-311++G**//B3LYP/6-311+G** adequately reproduces the results obtained at the higher level of theory CCSD(T)/6-311++G*s*//MP2/6-31G* even though it requires considerably shorter CPU times and smaller amounts of memory.     

A quantum mechanical approach to the kinetics of the hydrogen abstraction reaction H2O2 + •OH → HO2 + H2O

The kinetics of the hydrogen abstraction from H₂O₂ by ^?OH has been modeled with MP2/6‐31G*//MP2/6‐31G*, MP2‐SAC//MP2/6‐31G*, MP2/6‐31+G**//MP2/6‐31+G**, MP2‐SAC// MP2/6‐31+G**, MP4(SDTQ)/6‐311G**//MP2/6‐31G*, CCSD(T)/6‐31G*//CCSD(T)/6‐31G*, CCSD(T)/6‐31G**//CCSD(T)/6‐31G**, CCSD(T)/6‐311++G**//MP2/6‐31G* in the gas phase. MD simulations have been used to generate initial geometries for the stationary points along the potential energy surface for hydrogen abstraction from H₂O₂. The effective fragment potential (EFP) has been used to optimize the relevant structures in solution. Furthermore, the IEFPCM model has been used for the supermolecules generated via MD calculations. IEFPCM/MP2/6‐31G* and IEFPCM/CCSD(T)/6‐31G* calculations have also been performed for structures without explicit water molecules. Experimentally, the rate constant for hydrogen abstraction by ^?OH drops from 1.75 × 10^?12 cm³ molecule^?1 s^?1 in the gas phase to 4.48 × 10^?14 cm³ molecule^?1 s^?1 in solution. The same trend has been reproduced best with MP4 (SDTQ)/6‐311G**//MP2/6‐31G* in the gas phase (0.415 × 10^?12 cm³ molecule^?1 s^?1) and with EFP (UHF/6‐31G*) in solution (3.23 × 10^?14 cm³ molecule^?1 s^?1). © 2005 Wiley Periodicals, Inc. Int J Chem Kinet 37: 502–514, 2005     

Theoretical study of light-metal tetrahydroaluminates

Ab initio calculations of the structural, energetic, vibrational, and magnetic characteristics of the lowest-lying structures for isolated molecules and ions of light-metal tetrahydroaluminates (Li, Na, Be, Mg, and Al) have been performed by the perturbation theory (MP2), quadratic configuration interaction (QCISD(T)), coupled cluster (CCSD(T)), and density functional theory (B3LYP) methods using the 6-31G*, 6-31G**, 6-311+G**, and 6-311++G** basis sets. The trends in the behavior of the molecular characteristics have been analyzed in various related series of these compounds. The results are compared with the data on analogous light-metal tetrahydroborates calculated at the same levels of theory. The differences in structure and stability between analogous hydroborate and alanate complexes are examined. The economical approximation B3LYP/6-311++G**//B3LYP/6-311+G** has been shown to adequately reproduce the results obtained at the higher level of theory CCSD(T)/6-311++G**//MP2/6-31G* even though it requires considerably shorter CPU times and smaller amounts of memory.     

A quest for triplet silylenes XHSi3 at ab initio and DFT levels (X = H, F, Cl and Br)

Four ground state triplet silylenes are found among 30 possible silylenic XHSi₃ structures (X = H, F, Cl and Br), at seven ab initio and DFT levels including: B3LYP/6-311++G^∗∗, HF/6-311++G^∗∗, MP3/6-311G^∗, MP2/6-311+G^∗∗, MP4(SDTQ)/6-311++G^∗∗, QCISD(T)/6-311++G^∗∗ and CCSD(T)/6-311++G^∗∗. The latter six methods indicate that the triplet states of 3-flouro-1,2,3-trisilapropadienylidene, 1-chloro-1,2,3-trisilapropargylene and 3-chloro-1,2,3-trisilapropargylene are energy minima. These triplets appear more stable than their corresponding singlet states which cannot even exist for showing negative force constants. Also, triplet state of 1-flouro-1,2,3-trisilapropargylene is possibly accessible for being an energy minimum, since its corresponding singlet state is not a real isomer. Some discrepancies are observed between energetic and/or structural results of DFT vs. ab initio data.     

Theoretical Studies on Intermolecular Hydrogen-bond Interactions between Hexamethylenetetramine and Nitric Acid

The structures of the complexes generated by hexamethylenetetramine and nitric acid have been fully optimized by B3LYP method at the 6-311++G** and aug-cc-pVTZ levels. The intermolecular hydrogen-bonding interactions have been calculated by the B3LYP/6-311++G**, B3LYP/aug-cc-pVTZ, MP2(full)/6-311++G** and CCSD(T)/6-311++G** methods, respectively. The NBO (nature bond orbital), AIM (atom in molecule), temperature effect and solvation effect have been analyzed to reveal the origin of the interactions. The results indicate that the stable hydrogen-bonded complexes could be generated by hexamethylenetetramine and nitric acid. The interactions follow the order of (a)>>(e)>(b)>(c)>(d)>(f)>(g). The C-N bonds which are adjacent to the methylene involving the hydrogen bonds tend to break in the chemical reaction. Due to the exothermic process, low temperature is conducive to the formation of the composition, which tallies with the experimental result.     

Tetrahedrane--dossier of an unknown

To probe whether tetrahedrane should be isolable the thermodynamics and kinetics of C4H4 singlet and triplet structures were studied extensively at the CCSD(T)/cc-pVTZ//CCSD(T)/cc-pVDZ, CCSD(T)/cc-pVDZ, CCSD(T)/cc-pVDZ//B3 LYP/6-311G**, and B3 LYP/6-311G** levels of theory. The reaction of cyclopropene with atomic carbon, which was previously suggested to involve tetrahedrane as a reactive intermediate, was re-examined experimentally with low-temperature matrix-isolation techniques. While experimental and theoretical results exclude the intermediacy of tetrahedrane in the above reaction, it is predicted to be an isolable molecule. Among the many C4H4 species, we pay special attention to the electronic effects on the ground state multiplicity of the respective carbenes.     

Theoretical study of [Si,O,C,O] species: Prediction of new species on triplet potential energy surface

The intermediates [Si,O,C,O] of the Si + CO₂ reaction have been studied in detail using high level ab iniitio methods. Both singlet and triplet [Si,O,C,O] species are characterized structurally and energetically. On the singlet potential energy surface (PES), the vdw‐OSi–CO isomer and in the triplet PES, the bent‐SiOCO isomer is found to be thermodynamically as well as kinetically most stable species. All possible isomerization transition states (TS) are located on both singlet and triplet potential surfaces. On the triplet surface, the stability of the trans‐OSiCO isomer is comparable with that of the bent‐SiOCO isomer. A non‐planar cis‐SiOCO isomer is located on the triplet PES, which is predicted for the first time. Heats of formation at 0 K (Δ_fH°, 0 K) for all singlet and triplet species are computed using G3B3, G3MP2, and CBS‐Q theories. The discrepancy between G3B3 and the other two methods for the heat of formation value for triplet trans‐OSiCO is discussed. The PESs for singlet as well as triplet species with their dissociation asymptotes are explored at the CCSD(T)/6‐311G(d,p)//MP2/6‐311G(d,p) level of theory. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Structural studies on ethyl isovalerate by microwave spectroscopy and quantum chemical calculations

We observed the microwave spectrum of ethyl isovalerate by molecular beam Fourier transform microwave spectroscopy. The rotational and centrifugal distortion constants of the most abundant conformer were determined. Its structure was investigated by comparison of the experimental rotational constants with those obtained by ab initio methods. In a first step, the rotational constants of various conformers were calculated at the MP2/6-311++G** level of theory. Surprisingly, no agreement with the experimental results was found. Therefore, we concluded that in the case of ethyl isovalerate more advanced quantum chemical methods are required to obtain a reliable molecular geometry. Ab initio calculations carried out at MP3/6-311++G**, MP4/6-311++G**, and CCSD/6-311++G** levels and also density functional theory calculations using the B3LYP/6-311++G** method gave similar results for the rotational constants, but they were clearly distinct from those obtained at the MP2/6-311++G** level. With use of these more advanced methods, the rotational constants of the lowest energy conformer were in good agreement with those obtained from the microwave spectrum.     

Direct dynamics study on mechanism and kinetics of the biradical self‐reaction of HOO

A theoretical study of the mechanism and the kinetics for the hydrogen abstraction reaction of the biradical hydroperoxy radical has been presented at the CCSD(T)/6‐311++G(3d,2p)//CCSD/6‐31+G(d,p) level of theory. Our theoretical calculations suppose a stepwise mechanism involving the formation of a postreactant complex in the triplet and singlet entrance channels. Four transition states of the six‐membered chain complexes (³TS1 and ¹TS1) and six‐membered ring complexes (³TS2 and ¹TS2) are located at the high dual level CCSD(T)/6‐311++G(3d,2p)//CCSD/6‐31+G(d,p) method. The rate constants of Path 1 ～ Path 4 at the CCSD(T)/6‐311++G(3d,2p)//CCSD/6‐31+G (d,p) level are calculated by means of the conventional transition state theory (TST) and canonical variational TST without and with small‐curvature tunneling (SCT) correction within the temperature range of 200–2,500 K. The calculated results show that the triplet channel is the dominating reaction channel and Path 2 is found to be the most favorable pathway. The rate constants of Path 2 are in good agreement with the experimental values at the experimentally measured temperatures. Moreover, the variational effect is not obvious in the low temperature range but is not neglectable in the high temperature range. The SCT plays an important role particularly in the low temperature range. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Computational mechanistic study of methanol and molecular oxygen reaction on the triplet and singlet potential energy surfaces

The reaction of CH₃OH with the O₂ on the triplet and singlet potential energy surfaces (PES) was carried out using the B3LYP, MP2, and CCSD(T)//B3LYP theoretical approaches in connection with the 6-311++G(3df–3pd) basis set. Three pre-reactive complexes, ¹C1, ¹C2, and ³C1, on the singlet and triplet PES were formed between methanol and molecular oxygen. From a variety of the complexes, seven types of products are obtained, of which four types are found to be thermodynamically stable. Results reveal that there exists one intersystem crossing between triplet and singlet PES. For P4 adduct that is the main and kinetically the most favorable product, the rate constants are calculated in the temperature range of 200–1,000 K in the reliable pathway.     

Theoretical study on the potential energy surface of NC<Subscript>3</Subscript>P isomers

Density functional theory (DFT) calculations have been used to study the isomerization process in the NC₃P system. At the DFT/B3LYP/6-311G(d) level, 28 triplet and 28 singlet minima were obtained on their respective potential energy surfaces. The linear triplet ³NCCCP is the lowest-energy structure among the isomers. On the triplet PES, only linear isomers ³NCCCP, ³CNCCP, ³CCCNP, and ³CCNCP possess great kinetic and thermodynamic stabilities to exist under low-temperature conditions (such as in the dense interstellar clouds). At the same time, one chain-like and four three-membered-ring isomers on the singlet PES have been located with high kinetic and thermodynamic stabilities. Further CCSD(T)/6-311G(2df)//QCISD/6-311G(d), CCSD(T)/cc-pVTZ//DFT/B3LYP/cc-pVTZ, and CASPT2(14,12)/cc-pVQZ//CASSCF(14,12)/cc-p VQZ calculations are performed on the structures, frequencies, and energies of the relevant species. The bonding natures were analyzed and the results were compared with the analogous NC₃N and NC₂P molecules so as to aid their future experimental or astrophysical detection.     

Rate coefficients of the CF3CHFCF3 + H → CF3CFCF3 + H2 reaction at different temperatures calculated by transition state theory with ab initio and DFT reaction paths

The minimum energy path (MEP) of the reaction, CF₃CHFCF₃ + H → transition state (TS) → CF₃CFCF₃ + H₂, has been computed at different ab initio levels and with density functional theory (DFT) using different functionals. The computed B3LYP/6‐31++G**, BH&HLYP/cc‐pVDZ, BMK/6‐31++G**, M05/6‐31+G**, M05‐2X/6‐31+G**, UMP2/6‐31++G**, PUMP2/6‐31++G**//UMP2/6‐31++G**, RCCSD(T)/aug‐cc‐pVDZ//UMP2/6‐31++G**, RCCSD(T)/aug‐cc‐pVTZ(spd,sp)//UMP2//6‐31++G**, RCCSD(T)/CBS//M05/6‐31+G**, and RCCSD(T)/CBS//UMP2/6‐31++G** MEPs, and associated gradients and Hessians, were used in reaction rate coefficient calculations based on the transition state theory (TST). Reaction rate coefficients were computed between 300 and 1500 K at various levels of TST, which include conventional TST, canonical variational TST (CVT) and improved CVT (ICVT), and with different tunneling corrections, namely, Wigner, zero‐curvature, and small‐curvature (SCT). The computed rate coefficients obtained at different ab initio, DFT and TST levels are compared with experimental values available in the 1000–1200 K temperature range. Based on the rate coefficients computed at the ICVT/SCT level, the highest TST level used in this study, the BH&HLYP functional performs best among all the functionals used, while the RCCSD(T)/CBS//MP2/6‐31++G** level is the best among all the ab initio levels used. Comparing computed reaction rate coefficients obtained at different levels of theory shows that, the computed barrier height has the strongest effect on the computed reaction rate coefficients as expected. Variational effects on the computed rate coefficients are found to be negligibly small. Although tunneling effects are relatively small at high temperatures (～1500 K), SCT corrections are significant at low temperatures (～300 K), and both barrier heights and the magnitudes of the imaginary frequencies affect SCT corrections. © 2012 Wiley Periodicals, Inc.     

Ab Initio Study on the Potential Energy Surfaces of B4 Cluster

The singlet and triplet potential energy surfaces (PES) for the isomerization and dissociation reactions of B₄ isomers have been investigated using ab initio methods. Ten B4 isomers have been identified and of these 10 species, 4 have not been reported previously. The singlet rhombic structure ¹1 is found to be the most stable on the B₄ surface, in agreement with the results of previous reports. Several isomerization and dissociation pathways have been found. On the singlet PES, the linear ¹3b can rearrange to rhombus ¹1 directly, while ¹3c rearranges to ¹1 through two‐step reactions involving a cyclic intermediate. On the triplet PES, the capped triangle structure ³2 undergoes ring opening to the linear isomer ³3b with a barrier of 34.8 kcal/mol and 44.9 kcal/mol, and the latter undergoes ring closure to the square structure ³1 with a barrier of 30.4 kcal/mol and 33.0 kcal/mol at the MP4/6–311+G(3df)//MP2/6–311G(d) and CCSD/aug‐cc‐pVTZ//MP2/6–311G(d) levels of theory, respectively. The direct decomposition of singlet B₄ yielding to B₃+B is shown to have a large endothermicity of 87.3 kcal/mol (CCSD), and that producing 2B₂ to have activation energy of 133.4 kcal/mol (CCSD).     

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

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

The C7-C10 cycloalkanes revisited

The conformers of cycloheptane through cyclodecane have been examined at the B3LYP/6-311+G* and MP2/6-311+G* theoretical levels, with some additional calculations at the CCD/6-311+G* and CCSD(T)/6-311++G** levels. With cyclooctane, B3LYP predicts that the boat-chair and crown conformers have similar energies, whereas MP2 and CCSD(T) predict that the crown conformer is 2 kcal/mol higher in energy. The latter is in agreement with the electron diffraction data. With cyclononane, B3LYP predicts that two of the higher-energy conformers found in molecular mechanics calculations should convert to one of the lower-energy conformers. However, MP2/6-311+G* optimizations find them to be true minima on the potential energy surface. B3LYP systematically predicts larger C-C-C bond angles for these compounds than either MP2 or CCD. The results of molecular mechanics MM4 calculations are generally in good agreement with those obtained using MP2.     

Tautomerism and the maximum hardness principle

Through the application of the atom–bond electronegativity equalization method (ABEEM) to the calculation of the hardnesses of more than 300 tautomers, it can be seen that the maximum hardness principle is nearly useless to account for their relative stabilities. Moreover, by calculating the energies of these tautomers with the HF, B3LYP, B3PW91, and MP2 methods at the 6‐31G, 6‐31G*, 6‐31G**, 6‐31+G**, 6‐311G**, or 6‐311++G** level, it is found that all these methods may not be always reliable in predicting their relative stabilities. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Singlet–triplet energy differences in divalent five membered cyclic conjugated Arduengo-type carbenes XC<Subscript>2</Subscript>HN<Subscript>2</Subscript>M (M = C,Si, Ge,Sn, and Pb; X = F,Cl, Br,and I)

Singlet-triplet energy differences in Arduengo-type carbenes XC₂HN₂C compared and contrasted with their sila, germa, stana and plumba analogues; at B3LYP/6-311++G** level of theory. Free Gibbs energy differences between triplet (t) and singlet (s) states (ΔG(t–s)) change in the following order: plumbylenes > stannylenes > germylenes > silylenes > carbenes. The singlet states in XC₂HN₂C are generally more stable when the electron withdrawing groups such as–F was used at β-position. However, the singlet states in XC₂N₂HM (M = Si, Ge, Sn, and Pb) are generally more stable when the withdrawing groups such as–F was placed. The puckering energy is investigated for each the singlet and triplet states. The DFT calculations found the linear correlation to size of the group 14 divalent element (M), the ∠N–M–N angle, and the Δ(LUMO–HOMO) of XC₂HN₂M.     

Thermochemical properties of molecules and cations of MgNnHm (n=1,2 and m=1–6): enthalpies of formation, proton affinities and ionization energies

Ab initio molecular orbital calculations are reported for small neutral molecules and cations containing magnesium, nitrogen and hydrogen. Structures have been optimized using gradient techniques at B3LYP/6-31+G(d) and at MP2(full)/6-311++G(d,p). Single-point calculations are reported at QCISD(T)(full)/6-311++G(2df,p) and at CCSD(T)(full)/6-311++G(2df,p) levels using geometries optimized at MP2(full)/6-311++G(d,p). Standard enthalpies of formation at 298 K have been calculated at these two higher levels of theory. Other thermochemical properties calculated include ionization energies and proton affinities. The binding enthalpies of ammonia to Mg⁺, MgNH₂⁺ and MgNH₃⁺ are also reported.