Ab initio investigation on stability and properties of XYCO... HZ complexes. II: Post hartree-fock studies on H2CO... HF

Ab initio MP2 level of theory in conjunction with three basis sets of a triple-zeta quality was applied to study the molecular geometry and stability of the H₂CO... HF complex. An interaction energy predicted for this system at the highest, MP4(SDTQ)/6-311 + +G(2df, 2pd)//MP2/6-311 + +G(2df, 2pd), level corrected for the BSSE and ZPE contributions amounts to -4.85 kcal/ mol. BSSE contributes significantly to the interaction energies at all applied levels. Reliable MP2/ 6-311 + +G(2df, 2pd) level harmonic vibrational frequencies, IR intensities, and the predicted isotopic shifts upon deuteration and¹⁸O substitution are presented in order to facilitate experimental studies on the IR spectrum of the title complex.     

Hetero‐π‐systems from 2 + 2 cycloreversion,part 2.1Ab initio thermochemical study of heterocyclobutanes 2 + 2 cycloreversion to form heteroethenes H2C=X (X=NH,O, SiH2, PH,S)

Ab initio and DFT thermochemical study of diradical mechanism of 2 + 2 cycloreversion of parent heterocyclobutanes and 1,3‐diheterocyclobutanes, cyclo‐(CH₂CH₂CH₂X), and cyclo‐(CH₂XCH₂X), where X = NH, O, SiH₂, PH, S, was undertaken by calculating closed‐shell singlet molecules at three levels of theory: MP4/6‐311G(d)//MP2/6‐31G(d)+ZPE, MP4/6‐311G(d,p)//MP2/6‐31G (d,p)+ZPE, and B3LYP/6‐311+G(d,p)+ZPE. The enthalpies of 2 + 2 cycloreversion decrease on going from group 14 to group 16 elements, being substantially higher for the second row elements. Normally endothermic 2 + 2 cycloreversion is predicted to be exothermic for 1,3‐diazetidine and 1,3‐dioxtane. Strain energies of the four‐membered rings were calculated via the appropriate homodesmic reactions. The enthalpies of ring opening via the every possible one‐bond homolysis that results in the formation of the corresponding 1,4‐diradical were found by subtracting the strain energies from the central bond dissociation energies of the heterobutanes CH₃CH₂—CH₂XH, CH₃CH₂—XCH₃, and HXCH₂—XCH₃. The latter energies were determined via the enthalpies of the appropriate dehydrocondensation reactions, using C—H and X—H bond energies in CH₃XH calculated at G2 level of theory. Except 1,3‐disiletane, in which ring‐opening enthalpy attains 69.7 kcal/mol, the enthalpies of the most economical ring openings do not exceed 60.7 kcal/mol. The 1,4‐diradical decomposition enthalpies found as differences between 2 + 2 cycloreversion and ring‐opening enthalpies were negative, the least exothermicity was calculated for ⋅ CH₂SiH₂CH₂CH₂. The only exception was 1,3‐disiletane, which being diradical, CH₂SiH₂CH₂SiH₂, decomposed endothermically. Since decomposition of the diradical containing two silicon atoms required extra energy, raising the enthalpy of the overall reaction to 78.9 kcal/mol, 1,3‐disiletane was predicted to be highly resisting to 2 + 2 cycloreversion. © 2007 Wiley Periodicals, Inc. Heteroatom Chem 18:704–720, 2007; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/hc.20377     

Ab initio/GIAO-CCSD(T) study of structures, energies, and 13C NMR chemical shifts of C4H7(+) and C5H9(+) ions: relative stability and dynamic aspects of the cyclopropylcarbinyl vs bicyclobutonium ions

The structures and energies of the carbocations C 4H 7 (+) and C 5H 9 (+) were calculated using the ab initio method. The (13)C NMR chemical shifts of the carbocations were calculated using the GIAO-CCSD(T) method. The pisigma-delocalized bisected cyclopropylcarbinyl cation, 1 and nonclassical bicyclobutonium ion, 2 were found to be the minima for C 4H 7 (+) at the MP2/cc-pVTZ level. At the MP4(SDTQ)/cc-pVTZ//MP2/cc-pVTZ + ZPE level the structure 2 is 0.4 kcal/mol more stable than the structure 1. The (13)C NMR chemical shifts of 1 and 2 were calculated by the GIAO-CCSD(T) method. Based on relative energies and (13)C NMR chemical shift calculations, an equilibrium involving the 1 and 2 in superacid solutions is most likely responsible for the experimentally observed (13)C NMR chemical shifts, with the latter as the predominant equilibrating species. The alpha-methylcyclopropylcarbinyl cation, 4, and nonclassical bicyclobutonium ion, 5, were found to be the minima for C 5H 9 (+) at the MP2/cc-pVTZ level. At the MP4(SDTQ)/cc-pVTZ//MP2/cc-pVTZ + ZPE level ion 5 is 5.9 kcal/mol more stable than the structure 4. The calculated (13)C NMR chemical shifts of 5 agree rather well with the experimental values of C 5H 9 (+).     

Theoretical studies on electron delocalisation in selenourea

Ab initio and density functional calculations have been performed on the different possible structures of selenourea(su), urea(u) and thiourea(tu) to understand the extent of delocalisation in selenourea in comparison to urea and thiourea. Selenourea(su-1) withC ₂ symmetry has the minima on the potential energy surface at MP2(fu)/6-31+G* level. The C-N rotational barrier in selenourea is 8.69 kcal/mol, which is 0.29 and 0.11 kcal/mol more than that of urea and thiourea respectively at MP2(fu)/6-31+G* level. N-inversion barrier is 0.55 kcal/mol at MP2(fu)6-31+G* level. NBO analysis has been carried out to understand the nature of different interactions responsible for the electron delocalisation.     

Theoretical investigation of the equilibrium structure of the FSO<Subscript>3</Subscript> radical

Six different theoretical approaches (MP2, DFT-B3LYP/MPW3LYP, QCISD, UHF/ROHF-CCSD, UHF/ROHF-CCSD[T], and EOMIP-CCSD) were employed for a detailed theoretical investigation of a minimum on the potential surface of the fluorosulfate radical FSO₃ ground state. Optimized geometries, spin distributions, vibrational frequencies, and centrifugal distortion constants were calculated within C _3v and C ₁ initial symmetry constraints.     

Degenerate lithium-hydrogen exchange reactions: Ab initio models for metallation mechanisms involving H2, CH4, NH3, H2O,and HF

Hydrogen exchange reactions between lithium and sodium compounds, MX (M=Li: X=H, CH₃, NH₂, OH, F; M=Na: X=CH₃), and the corresponding hydrides, HX, have been modelled by means of ab initio calculations including electron correlation and zero point energy (ZPE) corrections. Small or no activation barriers (from the initial complexes) are encountered in systems involving lone pairs (10.8, 2.4, 0.0 kcal/mol for X=NH₂, OH, F, respectively). Since the association energies of the initial complexes are much larger (21.0, 20.4, 23.5 kcal/mol, respectively; MP2/6–31+G*/6–31+G* + ZPE), such exchange reactions should occur spontaneously in the gas phase. The methyl systems (X=CH₃) have the largest barriers: 26.7 (M=Li) and 31.7 (M=Na) kcal/mol (MP2/6–31+G*/6–31G* + ZPE), and the initial complexes are only weakly bound. The significance of these systems as models for hydrogen exchange reactions in complexes of electropositive transition metals is discussed. However, the gegenion-free exchange of hydrogen between CH₃ and CH₄ has a much lower, 11.8 kcal/mol barrier (MP2/6–31+G*/6–31+G* + ZPE). All the transition structures are highly ionic (charges on the metals > +0.8). The effect of aggregation has been considered by examining the hydrogen exchange between (LiX)₂ and HX(X=H, CH₃, NH₂, OH). Although these dimer reactions formally involve six, instead of four electrons, no “aromatic” preference is observed.     

Anomeric effect and theoretical vibrational spectra compared to the experiment in 2-chloro-1,3,2-dioxaphosphorinane-2-oxide, -sulfide,and -selenide. I. Anomeric effect and structures

2-Chloro-1,3,2-dioxaphosphorinane-2-oxide, -sulfide, and -selenide are studied at the DFT/B3LYP level and several ab initio methods using a 6-311G** basis set. Our energy optimizations by all these methods show that for oxide DFT and ab initio methods are not much different, while for the sulfide and the selenide the DFT relative energies are higher by about a kcal/mol as compared to those of MP2, MP3, MP4(SDTQ)//MP2, and CCSD(T)//MP2 (//MP2 indicates that a single-point calculation based on the MP2 optimized geometry is performed). However, regardless of rather large relative energies, that does not change the fact that in all three cases the conformational equilibrium mixture contains more than 95% of the lowest, chair-equatorial conformer (this indicates that the P=X bond is in equatorial position). This one and the next higher conformer (chair-axial) are confirmed to be real conformers (energy minima) in all cases. The energetically much higher twist and boat forms are probably just stationary states and local maxima because in many cases, geometry optimizations do not converge to them. Only for MP2 and the selenide do all optimizations converge to the desired stationary state. The relative energies could all be explained in terms of anomeric effects and ring strains. The decreasing covalent character of the P=X bond, with X changing from O to S and to Se, shows itself in the increasing bond lengths and the decreasing strength of anomeric effects.     

Ab Initio Calculations of Structure and Energetics of the Hypermetalated Hydrogen Oxide Molecules: Li2NaOH and LiNa2OH

We made ab initio electronic calculations of the structure and energetics of mixed hypermetalated hydrogen oxides, Li₂NaOH and LiNa₂OH. There exist five equilibrium geometries for each complex. In all levels of calculation the global minimum structure for Li₂NaOH has C_2v symmetry and a large distance between sodium and oxygen, 4.24 Å (MP2/6-31G*). The dissociation energies to all possible products were also calculated. Li₂NaOH → Na + Li₂OH δH = +25.33 kcal/mol (at MP4/6-311++G**//6-31G* + ZPE scaled by 0.9). All other dissociation processes are highly endothermic. Similar procedures were applied to LiNa₂OH. The global minimum structure for LiNa₂OH belongs to point group C_s. It is also endothermic to all possible dissociation paths. LiNa₂OH →Na + LiNaOH δH = +12.72 kcal/mol (at MP4/6-311++G*//6-31G* + ZPE scaled by 0.9). The nuclear repulsion energy is crucial in energetics of the structures. The distribution of electron density and bonding properties for these equilibrium structures were analyzed.     

The gas phase hydrogen‐bonded dimers of HOCl: A high‐level quantum chemical study

The O···H? O and Cl···H? O hydrogen bonding interactions were analyzed for HOCl dimers by using B3LYP, MP2, CCSD, and MP4(SDTQ) methods in conjunction with the various basis sets. Five isomers were found for the HOCl dimer. The ZPE and BSSE corrected binding energies were computed at the different levels of theory. At the optimized geometries obtained at CCSD/AUG‐cc‐pVDZ level, energies were re‐evaluated at MP4(SDTQ)/AUG‐cc‐pVTZ and CCSD(T)/cc‐pVTZ levels of theory. We found an average of ?20.9 and ?9.6 kJ/mol for the strength of the O···H and Cl···H hydrogen bonding interactions, respectively. Excitation and vertical ionization energies as well as rotational constants were computed at different levels of theory. The quantum theory of atoms in molecules (QTAIM) and natural bond orbital (NBO) analysis were used to elucidate the nature of the interactions of HOCl dimers. The interaction energies were decomposed by Morokuma methodology. We have computed Δ_fH°(HOCl) and Δ_fH°(HOCl⁺) using the atomization reactions. The Δ_fH°₂₉₈(HOCl) values are ?17.85 and ?18.05 kcal/mol by using CBS‐Q and CBS‐QB3 extrapolation models, respectively, in good agreement with the results given in JANAF tables. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Theoretical study on structures and stability of SiCP2 isomers

The structures, energetics, spectral parameters and stability of the singlet SiCP₂ isomers are explored at the density functional theory and ab initio levels. Eight isomers connected by ten interconversion transition states are located at the CCSD(T)/6-311G(2d)//B3LYP/6-311G(d)level. The kinetically stable isomers and their relevant interconversion transition states are further refined at CCSD(T)/6-311+G(2df)//QCISD/6-311G(d) level. At QCISD/6-311G(d) level, one four-membered ring isomer cSiPCP and two linear structures PSiCP, SiCPP possess considerable kinetic stability (more than 15 kcal/mol). The valence bond structures of three kinetically stable SiCP₂ isomers are analyzed. The similarities and discrepancies in structure, energy and stability between SiCP₂ and its analogous C₂P₂, Si₂P₂, SiCN₂ and CSiNP molecules are also discussed. The predicted structures and spectroscopic properties are expected to be informative for the identification of the SiCP₂ in the laboratory and space.     

Does CH prefer a C2v rather than a Cs structure?

The closely related C_s ( 1 ) and C_2v ( 3 ) structures of CH have been reinvestigated at many ab initio levels using MP2/6-31G** and MP2/6-311 + + G(2df, 2pd) geometries. The largest basis sets employed were 6-311G(3df, 2p), 6-311 + + G(3df, 3pd), and the Dunning “correlation consistent” polarized triple-split valence basis set (cc-pVTZ). Electron correlation was probed at the MP4 level, but the QCISD method was also used with the largest basis sets. While electron correlation favors 3 over 1 by about 2 kcal/mol, the correlated relative energies with all basis sets employed range from 0.36–1.03 kcal/mol in favor of 1 . The best estimate of this difference, 0.86 kcal/mol, is essentially identical with the (scaled) zero-point energy difference, 0.84 kcal/mol, favoring 3 over 1 . These results indicate that 1 and 3 have almost exactly the same energy at 0 K. Our best value for the dissociation energy of CH is 42.0 kcal/mol [QCISD(T)/6-311 + + G(3df, 3pd)//MP2(fu)/6-311 + + G(2df, 2pd), corrected to 298 K], which agrees very well with the experimental value. © 1992 by John Wiley & Sons, Inc.     

Theoretical study of thermal decomposition of azomethane

Summary.  Thermal one- and two-bond dissociation processes of cis- and trans-azomethane were studied by ab initio computation with DZP and TZ2P basis sets, using the d(N–C) bond lengths as the reaction coordinates. The geometries were optimized at the MP2 level, and the dissociation energies obtained exploiting a single-point, fourth-order M?ller–Plesset calculations [MP4SDTQ/TZ2P]. At this level of theory including zero-point energies, the trans-isomer is by 9.3 kcal/mol more stable than the cis-isomer. The results show that the energetically more favourable one-bond cleavage proceeds without transition state with the predicted bond dissociation energy D ₀ of 47.8 kcal/mol for trans-azomethane and 38.5 kcal/mol for cis-azomethane. With calculated barrier heights the unimolecular dissociation rate constants have been determined by means of the RRKM theory. The second-order saddle points localized for synchronous decomposition pathways lie 13 (trans)-23(cis) kcal/mol above the one-bond dissociation energies [MP2/DZP]. Received May 28, 1996/Final version received November 1, 1996 / Accepted November 1, 1996     

Computational study on the mechanism and thermodynamic of atmospheric oxidation of HCN with ozone

Oxidation of hydrogen cyanide with ozone on the singlet potential energy surface is carried out using the MP2 and CCSD(T)//MP2 theoretical approaches in connection with the 6-311++G(d,p) basis set. In this reaction, energy barrier of transition states are low, so the reaction of HCN with ozone can occur easily at atmospheric condition. With variety of pre-reactive complex, five types of products are obtained at the MP2 method which four types of them have enough thermodynamic stability at the standard condition. CO₂ + HNO are final adducts in a process that is computed to be exothermic by ?132.605 kcal/mol in standard enthalpy and spontaneous by ?144.166 kcal/mol in Gibbs free energy of reaction. In kinetic viewpoint, the production of OCN + HO₂ adducts path with one low level transition state is the most favored path.     

Geometries and stabilities of various configurations of benzene dimer: details of novel V-shaped structure revealed

Benzene dimer configurations namely T-shaped, parallel-displaced, sandwich, and V-shaped that were proposed by experimental studies are investigated using second- and fourth-order Møller–Plesset perturbation theory. The MP2 method with aug-cc-pVDZ and aug-cc-pVTZ basis sets unequivocally shows that the parallel-displaced configuration is considerably more stable than T-shaped structure. On the other hand, the MP4(SDTQ)/aug-cc-pVDZ level predicts that the T-shaped and parallel-displaced configurations are nearly isoenergetic, which is parallel to the previous results of estimated CCSD(T)/CBS level reported recently. The lowest energy T-shaped configuration is stabilized by 0.17 kcal/mol over the parallel-displaced configuration at the MP4(SDTQ)/aug-cc-pVDZ level. Although the structures of all the four different types of configurations are found to be stable at both MP2 and full MP4 methods, the V-shaped configuration is the least stable among them. The calculated interaction energy of ?2.3 kcal/mol for the lowest energy T-shaped structure at the MP4(SDTQ)/aug-cc-pVDZ level is in good agreement with the experimental value of ?2.4 ± 0.4 kcal/mol. We conclude that the MP4(SDTQ) with a reasonably good basis set can be used for systems involving π–π interactions to obtain qualitative and quantitative results.     

Post-Hartree-Fock and DFT level studies on the Cl2CO …︁ Cl2 complex: Accurate molecular parameters,harmonic vibrational frequencies,and interaction energies

The Cl₂CO …︁ Cl₂ complex was studied using ab initio post-Hartree-Fock theory at the MP2 and MP4 levels and, for comparison, the DFT method with 6-311G(2d), 6-311 + G(2d), and Sadlej's medium-size polarized (MSPBS) basis sets. A potential energy search recovered a planar minimum-energy structure characterized by a bent conformation. For this weakly bound complex, the interaction energy corrected for the basis set superposition error amounted to − 0.88, − 1.09, − 1.43, and − 0.38 kcal/mol at the MP4(SDTQ)/6-311G(2d), MP4(SDTQ)/6-311 + G(2d), MP4(SDTQ)/MSPBS, and DFT(Becke3LYP)/6-311 + G(2d) levels of theory, respectively. Two highly symmetrical forms, linear and T-shaped, correspond to transition-state conformers. The analysis of harmonic vibrational frequencies and potential energy distribution was performed at the MP2 and DFT levels with the 6-311 + G(2d) basis set. © 1996 John Wiley & Sons, Inc.     

Theoretical investigation of the cyclic GaO2 and GaS2 molecules at DFT and correlated wave function levels

The geometries and the bonding properties have been predicted for cyclic GaO₂ and GaS₂ species at density functional theory (DFT), MPn (n=2,3,4 with different substitutions), QCISD(T), and CCSD(T) all‐electron correlation levels with 6‐311+G* basis set. The geometrical optimizations and the harmonic vibrational frequency analysis are performed using DFT and second‐order Møller–Plesset (MP2) methods. The relevant energy quantities are also calibrated at the high‐order electron correlation levels [MP3, MP4, quadratic configuration interaction (QCI), and coupled cluster (CC)]. Each species possesses a ²A₂ ground state with a higher energy level ²A₁ state. The corresponding state–state separations are about 32 kcal/mol for GaO₂ species and about 20 kcal/mol for GaS₂ species at the QCISD(T)/6‐311+G* level. The QCISD(T) and CCSD(T) calculations yield dissociation energies of 42.0 and 59.0 kcal/mol for two species, respectively, and other methods yield dissociation energies within ∼5 kcal/mol. Result analysis has indicated that the cyclic GaO₂ should be classified as superoxide and the GaS₂ species should be classified as supersulfide in their ground state, and those in the excited state (²A₁) should not be. However, the cyclic GaS₂ (²A₂) is less ionic than the GaO₂ (²A₂) and they are far less ionic than NaO₂. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 222–231, 2001     

Estimates of the structure and dimerization energy of polyacetylene from ab initio calculations on finite polyenes

Ab initio calculations at the Hartree-Fock (HF) and the second-order Møller-Plesset (MP2) levels are performed for finite polyenes C_2nH_2n+2 to estimate the structure and dimerization energy (E_dim) of polyacetylene. The effect of electron correlation on the structure of finite polyenes is analyzed in detail. The MP3/6–31G* C(DOUBLE BOND)C and C(SINGLE BOND)C bond lengths in polyacetylene are estimated by a simple extrapolation method using empirical corrections for the MP2 deficiencies, yielding values [C(DOUBLE BOND)C(MP3) ∼ 1.36 Å and C(SINGLE BOND)C(MP3) ∼ 1.44 Å] that are in a good agreement with experiment (C(DOUBLE BOND)C (DOUBLE BOND) 1.36 Å and C(SINGLE BOND)C (DOUBLE BOND) 1.44–1.45 Å). Comparison is also made with other theoretical estimates of polyacetylene structure. E_dim is approximated by the energy difference between the equilibrium and hypothetical polyenic structures. It is estimated that E_dim is ∼ 1.4–1.5 kcal/mol (0.06–0.07 eV) per carbon-carbon bond at the HF level with 4–21G and 6–31G* basis sets and ∼ 0.3–0.5 kcal/mol (0.013–0.022 eV) at the MP2 level with the 6–31G* basis set. It is concluded that E_dim is very sensitive to the level of approximation employed so that a proper treatment of electron correlation is essential to obtain a reliable estimate of the dimerization energy. © 1997 John Wiley & Sons, Inc.     

Theoretical refinements of theendo and theexo structures of tetraborane(8) carbonyl

The molecular structures of theendo (1a) andexo (1b) isomers of B₄H₈CO have been optimized at the ab initio MP2(Full)/6-31G^* level of theory. The agreement of the computed geometrical parameters with the recently published electron-diffraction (GED) data is very good, even though a number of geometrical constraints were applied in the experimental determination. The IGLO (individual gauge for localized orbitals)¹¹B NMR chemical shifts, calculated at the II//MP2/6-31G^* level, are also in accord with experiment. The formation of1a and1b by association of B₄H₈ and CO is computed to be exothermic by 22.8 and 22.2 kcal/mol, respectively, at the MP2(Full)/6-31G^*//MP2(Full)/6-31G^* + ZPE(6-31G^*) level of theory. The Lewis acid strength of B₄H₈ toward CO is comparable to that of BH₃.     

Ab initio study of the inversion barrier in NF3+

Ab initio Hartree–Fock, Møller–Plesset perturbation theory (MP 2), and quadratic configuration interaction, using single and double substitutions (QCISD ), calculations were carried out for the NF₃⁺ ion. Optimized structures were examined at the various levels of theory. Calculation of the inversion barrier height shows the importance of optimizing the geometry at the post-Hartree–Fock level and the inclusion of polarization functions. The best calculated inversion barrier was 13.3 kcal/mol, compared to an experimental value of 17.3 kcal/mol. The dissociation transition state was computed to determine the well depth of the NF₃⁺ ion and its stability toward dissociation. The computed well depth was 28 and 48 kcal/mol at the SCF and MP 2 levels, respectively. © 1994 John Wiley & Sons, Inc.     

Structures and Stability of NHS2 Isomers

Potential energy surface of HNS2 is investigated by means of second-order Moller-Plesset Perturbation theory(MP2) and QCISD(T) (single-point)methods,At final QCISD (T)/6-311 G(3df,2p)//MP2/6-311 G(d,p),level with zero-point vibrational energies included,cis-NHSS is found to be global minimum on the potential energy sufrace,followed by low-lying trans-HNSS ,HN(S)S(C2v).cis-HSNS,cis-NSSN,trans-HSNS,trans-NSSN,and HN(S)S(Cs) by 13.46,66.92,78.25,80.38,81.22,81.38 and 86.40 kJ/mol,respectively.A new high-lying HS(N)S isomer with Cs symmetry is located on the potential energy surface ,The kinetic stabilities of all isomers are predicted.Comprisons are made for HNS2 with its analogues,NHO2,HPS2 and HPO2,The causes that lead to the differences between HNS2 and its analogues are hypervalent capacity of phosphorus and distinct electronegativities of hydrogen,nitrogen and phosphorus.