An Ab Initio Study of the Complexes Formed by the Positive Acetylene Ion with the Hydrogen and Nitrogen Molecules and the Argon Atom

The complexes formed by the positive acetylene ion with the hydrogen molecule, the nitrogen molecule, and the argon atom are investigated with ab initio calculations using the 6-311G** and the 6-31+G(2df,2pd) basis sets. MP2/6-311G** energies and optimum geometries are obtained, as well as single-point MP3, MP4, and QCISD(T) energies with the MP2/6-311G** optimized geometries. Single-point calculations are performed with the 6-31+G(2df,2pd) basis set at MP2/6-311G** optimized geometries.     

An ab initio study of the structure,dissociation energy,and heat of formation of Na2S

The equilibrium geometries and fundamental frequencies of Na₂S are calculated at HF, MP2(FC, FU), and MP3 with the 6–31G(d) basis set and at HF and MP2(FC, FU) with the 6–31G(d) basis set, respectively. The total energy at MP2(FU)/6–31G(d)-optimized geometry is computed at MP4 with 6–311G(d, p), 6–311 + G(d, p), and 6–311G(2df, p), at QCISD(T)/6–311G(d, p), and at MP2/6–311G(3df, 2p) levels, respectively. The dissociation energy, the atomization energy, and the heat of formation for Na₂S are evaluated using the G1 and G2 models. The calculated results indicated that Na₂S in its ground state was a bent structure (C_2v). Electron correlation corrections on the bending angle are very significant. The equilibrium geometrical parameters are R_e(Na-S) = 2.45 Å and ∠Na-S-Na = 111.13° at the MP2(FU)/6–31G(d) level. The theoretically estimated dissociation energy, total atomization energy, and heat of formation are 67.07, 117.55, and 0.35 kcal mol⁻¹, respectively, at 298.15 K. © 1997 John Wiley & Sons, Inc.     

The structure of 1-chlorosilatrane: An ab initio molecular orbital and a density functional theory study

The molecular geometries of the 1-chloro-, 1-fluoro-, 1-methyl-, and 1-hydrogenosilatranes were fully optimized by the restricted Hartree-Fock (HF) method supplemented with 3-21G, 3-21G(d), 6-31G(d), and CEP-31G(d) basis sets; by MP2 calculations using 6-31G(d) and CEP-31G(d) basis sets; and by GGA-DFT calculations using 6-31G(d5) basis set with the aim of locating the positions of the local minima on the energy hypersurface. The HF/6-31G(d) calculations predict long (>254 pm) and the MP2/CEP calculations predicted short (∼225 pm) equilibrium Si(SINGLE BOND)N distances. The present GGA-DFT calculations reproduce the available gas phase experimental Si(SINGLE BOND)N distances correctly. The solid phase experimental results predict that the Si(SINGLE BOND)N distance is shorter in 1-chlorosilatrane than in 1-fluorosilatrane. In this respect the HF results show a strong basis set dependence, the MP2/CEP results contradict the experiment, and the GGA-DFT results in electrolytic medium agree with the experiment. The latter calculations predict that 1-chlorosilatrane is more polarizable than 1-fluorosilatrane and also support a general Si(SINGLE BOND)N distance shortening trend for silatranes during the transition from gas phase to polar liquid or solid phase. The calculations predict that the ethoxy links of the silatrane skeleton are flexible. Consequently, it is difficult to measure experimentally the related bond lengths and bond and torsion angles. This is the probable origin of the surprisingly large differences for the experimental structural parameters. On the basis of experimental analogies, ab initio calculations, and density functional theory (DFT) calculations, a gas phase equilibrium (r_e) geometry is predicted for 1-chlorosilatrane. The semiempirical methods predict a so-called exo minimum (at above 310 pm Si(SINGLE BOND)N distance); however, the ab initio and GGA-DFT calculations suggest that this form is nonexistent. The GGA-DFT geometry optima were characterized by frequency analysis. © 1996 by John Wiley & Sons, Inc.     

Theoretical studies of organometallic compounds. II. All electron and pseudopotential calculations of M(CH3)nCl4 − n (M = C,Si, Ge,Sn, Pb; n = 0–4)

The performance of effective core potentials (ECP) for the main group elements of group IV has been studied by calculating the geometries and reaction energies of isodesmic reactions for the molecules M(CH₃)_nCl_{4 ? n} (M = C, Si, Ge, Sn, Pb; n = 0–4) at the Hartree–Fock level of theory. The results are compared with data from all electron calculations and experimental results as far as available. The all electron calculations were performed with a 3-21G(d) and a 6-31G(d) basis set for Si, a (43321/4321/41) basis set for Ge, and a (433321/43321/431) basis set for Sn. For the ECP calculations the potentials developed by Hay and Wadt with a configuration (n)s^a(n)p^b and the valence basis set (21/21), extended by a set of d functions, are employed. © 1992 by John Wiley & Sons, Inc.     

A comparative quantum mechanical study of bond separation energies as a measure of cyclic conjugation

Restricted Hartree-Fock (RHF), second-order Møller-Plesset (MP2), and density functional calculations [using the Becke/Lee-Yang-Parr (B-LYP) exchange/correlation gradient-corrected functionals] employing the 6-311G(d, p) and 6-311 + + G(d, p) basis sets have been carried out to calculate isodesmic bond separation energies for reactions involving a number of representative five- and six-membered ring organic compounds. The MP2 and density functional approaches yield reasonably good energies; the density functional method agrees particularly well with experiment, exhibiting a root-mean-square error of only 2.5 kcal/mol. Ring geometries are calculated satisfactorily in all approaches but are given particularly accurately by the MP2 approach. A comparison of the B-LYP bond separation energies with several other definitions of resonance energy shows that these different approaches correlate with each other in a reasonable fashion. © 1995 John Wiley & Sons, Inc.     

Theoretical studies of organometallic compounds. I. All electron and pseudopotential calculations of Ti(CH3)nCl4 − n (n = 0–4)

The performance of effective core potentials (ECP) and model potentials (MP) has been studied by calculating the geometries and reaction energies of isodesmic reactions for the molecules Ti(CH₃)_nCl_{4 ? n} (n = 0–4) at the Hartree–Fock level of theory. The results are compared with data from all electron calculations and experimental results as far as available. The all electron calculations were performed with a 3-21G basis set from Hehre and a (53321/521/41) basis set from Huzinaga. For the ECP calculations the potentials developed by Hay and Wadt, and for the MP calculations, the model potentials developed by Sakai and Huzinaga, are employed. © 1992 by John Wiley & Sons, Inc.     

An ab initio study of the conformational preferences of Hoechst 33258 in gas-phase and aqueous solution environments

A quantum mechanical study of the conformational preferences of Hoechst 33258, a synthetic minor groove-binding drug, has been performed in both gas-phase and aqueous solution. Gas-phase calculations were performed at the HF/6-31G(d) and MP2/6-31G(d) levels of theory, whereas calculations in the aqueous solution phase were performed using the PCM model with the 6-31G(d) basis set. The molecule was divided into three fragments, which were submitted to a systematic and detailed conformational study. The results clearly indicate that Hoechst 33258 does not adopt a planar conformation in either the gas-phase or aqueous solution. Thus, a folded conformation is not induced by binding of the molecule to DNA, but is an intrinsic property of the compound. Received: 3 March 1998 / Accepted: 29 May 1998 / Published online 19 August 1998     

Reaction energetics of tetrahedrane and other hydrocarbons: Ab initio and density functional treatments

Reaction energetics of the highly strained tetrahedrane molecule, as well as some smaller hydrocarbons, were computed at the ab initio MP4, QCISD (T), G1, and G2 levels and also by several density functional (DF) approaches: VWN(LDA), PW86, BLYP, and B3LYP. A variety of basis sets were used in the DF calculations. For atomization processes, nonlocal DF procedures are superior to MP4 and QCI. For the other reactions, B3LYP is the most reliable DF approach and is overall competitive with MP4 and QCISD (T), although the MP4/6-311G(2df, p) results are generally the closest to the experimental and/or the G2 values. There is no consistent correlation between the accuracy of the DF results and the size of the basis set. © 1996 John Wiley & Sons, Inc.     

Optical rotation and two chiral C atoms of aliphatic cyclic dipeptides

Accurate geometries structures and total energies have been determined for the conformers of cyclo(L-Pro-Gly), cyclo(L-Ala-L-Ala), and cyclo(L-Pro-Ala) in the gaseous phase, using HF and B3LYP correlation methods at 6−31++G(d), 6−311++G(d, p), 6−311++G(2d, 2p) and aug-cc-pvdz basis sets. High level computations MP2 with 6−311++G(2d, 2p) basis set indicate that the relative stabilities of the available conformers can be determined correctly at the B3LYP/6−311++g(2d, 2p) level of theory. We have also described the implementation of DFT and HF theory for calculations of the optical rotation at 589.3 nm. In L-Ala-L-Ala, and L-Pro-Ala molecules, they have two chiral C (C*), so we discuss the different effect of two chiral C to optical activity of cydo(L-Pro-Gly), cyclo(L-Ala-L-Ala), and cyclo(L-Pro-Ala).     

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     

Density functional calculations of difluorodiazete structures with Gaussian-orbital-type approach

The comparison of local nonlocal and hybrid DFT methods with RHF , MP 2, CCSD , and CPF ab initio methods in generating geometries and relative energies of cis- and trans-difluorodiazete, SVWN , BLYP , BP 86, BECKE 3LYP , and BECKE 3P 86 DFT methods with 6-311 + g(2d) and 6-311 + + g(3df) basis sets. The geometries generated with RHF ab initio models are quite different from experimental values and energy evaluation prefers the wrong isomer. The hybrid methods give superior geometries while energies evaluated with nonlocal DFT methods are better than the one obtained with MP 2 or CCDS ab initio methods. The results suggest DFT as the method of choice of studying similar systems. © 1996 John Wiley & Sons, Inc.     

The polarisability of Hg and the ground-state interaction potential of Hg2

Summary The Hg atom electric polarisability at MP3, MP4, QCISD or QCISD(T) level using a 20 valence-electron pseudopotential and a (9s8p7d2f)/[7s6p4d2f] basis set converges and agrees with experiment. The calculated Hg₂ X ¹ V (R) at the same level is roughly half of experiment atR _e but twice as large as the experimentalC ₆ at largeR.Dedicated to Inga Fischer-Hjalmars whose wit and charm made Quantum Chemistry a better place On leave of absence from: Department of Chemistry, Peking University, Beijing 100871, People's Republic of China     

Ab initio studies of tetrafluoroethylene and tetratrifluoromethylethylene molecules in the ground and “twisted” states

Ab initio calculations with full optimization of geometry have been carried out with the 6–31 G^* basis set on tetrafluoroethylene (with the unrestricted Hartree-Fock method—UHF and the second-order Moller-Plesset perturbation theory—MP2) and tetratrifluoromethylethylene (with UHF) molecules in the singlet ground and triplet biradical states. The symmetry of the tetrafluoroethylene molecule in the triplet biradical state was demonstrated to differ from that of ethylene (D _2d ) due to the deviation of fluorine atoms from CCFF plane. The MP2 optimized geometries of ethylene and tetrafluoroethylene were used for higher level calculations (MP3, MP4, CCSD). The energy of the ground state singlet-biradical triplet splitting decreases in the series: ethylene>tetrafluoroethylene> tetratrifluoromethylethylene. These data on energy splitting explain the increase in reactivity toward the [2+2]-cycloaddition on going from ethylene to tetrafluoroethylene. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 4, pp. 605–607, April, 1998.     

Density functional Gaussian-type orbital approach in theoretical study of S2F2 isomerization

The structures of two isomers, difluorodisulfane (FSSF) and thiothionylfluoride (SSF₂), and the corresponding transition structure were generated with density functional theory (DFT) methods. Three groups of DFT methods were used: local (Local Spin Density Approximation, LSDA), nonlocal (local with gradient corrections; BLYP and BP86), and hybrid methods that include a mixture of Hartree-Fock (HF) exchange with nonlocal correlation (Becke3BLYP, Becke3P86). An extended basis set [6-311 + + G(3df)] was used for all calculations, although satisfactory results can be obtained with the 6-311G(d) basis set. The geometries obtained were compared with both restricted Hartree-Fock (RHF) calculated and experimentally obtained values. The energy outcome and the activation barrier for the isomerization were evaluated. It was determined that excellent geometries can be obtained with the Becke3B86 hybrid method, whereas for reasonable energies MP2 single-point calculations on these geometries are necessary. © 1996 by John Wiley & Sons, Inc.     

Complexes of ammonia with propane and cyclopropane: electrostatic guidelines for ab initio treatment

A model based on the molecular electrostatic potential (MESP) is employed for the investigation of structures and energies of complexes of ammonia with propane and cyclopropane. The electrostatic model geometries are employed as starting points for an ab initio investigation at the self-consistent field and second-order M?ller-Plesset (MP2) levels. The most stable structures of C₃H₆..NH₃ and C₃H₈..NH₃ complexes have the interaction energies of 10.07 kJ/mol and 8.15 kJ/mol, respectively, at the MP2/6-31G(d,p) level. The energy rank order of the structures is not altered with the use of the 6-31++G(d,p) basis set, and the basis␣set superposition error has little effect. The interaction energy decomposition analysis shows that the electrostatic component is dominant over the other ones. MESP topography thus seems to offer valuable hints for predicting the structures of weakly bonded complexes. Received: 8 July 1998 / Accepted: 4 August 1998 / Published online: 2 November 1998     

On the additivity of basis set effects in some simple fluorine containing systems

The reactions F + H₂ → HF + H, HF → H + F, F → F⁺ + e^? and F + e^? → F^? were used as simple test cases to assess the additivity of basis set effects on reaction energetics computed at the MP4 level. The 6-31G and 6-311G basis sets were augmented with 1, 2, and 3 sets of polarization functions, higher angular momentum polarization functions, and diffuse functions (27 basis sets from 6-31Gd, p) to 6-31 ++ G(3df, 3pd) and likewise for the 6-311G series). For both series substantial nonadditivity was found between diffuse functions on the heavy atom and multiple polarization functions (e.g., 6-31 + G(3d, 3p) vs. 6-31 + G(d, p) and 6-31G(3d, 3p)). For the 6-311G series there is an extra nonadditivity between d functions on hydrogen and multiple polarization functions. Provided that these interactions are taken into account, the remaining basis set effects are additive to within ±0.5 kcal/mol for the reactions considered. Large basis set MP4 calculations can also be estimated to within ±0.5 kcal/mol using MP2 calculations, est. E_MP4(6-31 ++ G(3df, 3pd)) ≈ E_MP4(6-31G(d, p)) + E_MP2(6-31 ++ G(3df, 3pd)) – E_MP2(6-31G(d, p)) or E_MP4(6-31 + G(d, p) + E_MP2(6-31 ++ G(3df, 3pd)) – E_MP2(6-31 + G(d, p)) and likewise for the 6-311G series.     

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     

Analysis of fourth-order Møller–Plesset limit energies: the importance of three-electron correlation effects

Fourth-order M?ller–Plesset (MP4) correlation energies are computed for 28 atoms and simple molecules employing Dunning's correlation-consistent polarized-valence m-zeta basis sets for m=2, 3, 4, and 5. Extrapolation formulas are used to predict MP4 energies for infinitely large basis sets. It is shown that both total and partial MP4 correlation energies can be extrapolated to limit values and that the sum of extrapolated partial MP4 energies equals the extrapolated total MP4 correlation energy within calculational accuracy. Therefore, partial MP4 correlation energies can be presented in the form of an MP4 spectrum reflecting the relative importance of different correlation effects. Typical trends in calculated correlation effects for a given class of electron systems are independent of the basis set used. As first found by Cremer and He [(1996) J Phys Chem 100:6173], one can use MP4 spectra to distinguish between electron systems with well-separated electron pairs and systems for which electrons cluster in a confined region of atomic or molecular space. MP4 spectra for increasing size of the basis set reveal that smaller basis set calculations underestimate the importance of three-electron correlation effects for both classes by overestimating the importance of pair correlation effects. The minimum size of a basis set required for reliable MP4 calculations is given by a valence triple-zeta polarized basis, which even in the case of anions performs better than a valence double-zeta basis augmented by diffuse functions. Received: 14 June 2000 / Accepted: 16 June 2000 / Published online: 24 October 2000     

The effect of basis set superposition error on the convergence of interaction energies

 For the intermolecular interaction energies of ion-water clusters [OH⁻(H₂O) _n (n=1,2), F⁻(H₂O), Cl⁻(H₂O), H₃O⁺(H₂O) _n (n=1,2), and NH₄ ⁺(H₂O) _n (n=1,2)] calculated with correlation-consistent basis sets at MP2, MP4, QCISD(T), and CCSD(T) levels, the basis set superposition error is nearly zero in the complete basis set (CBS) limit. That is, the counterpoise-uncorrected intermolecular interaction energies are nearly equal to the counterpoise-corrected intermolecular interaction energies in the CBS limit. When the basis set is smaller, the counterpoise-uncorrected intermolecular interaction energies are more reliable than the counterpoise-corrected intermolecular interaction energies. The counterpoise-uncorrected intermolecular interaction energies evaluated using the MP2/aug-cc-pVDZ level is reliable. Received: 14 March 2001 / Accepted: 25 April 2001 / Published online: 9 August 2001     

Ab initio study of hydrogen bonding in the phenol–water system

Three hydrogen-bonded minima on the phenol-water, C₆H₅OH—H₂O, potential energy surface were located with 3–21G and 6–31G** basis sets at both Hartree–Fock and MP2 levels of theory. MP2 binding energies were computed using large “correlation consistent” basis sets that included extra diffuse functions on all atoms. An estimate of the effect of expanding the basis set to the triple-zeta level (multiple f functions on carbon and oxygen and multiple d functions on hydrogen) was derived from calculations on a related prototype system. The best estimates of the electronic binding energies for the three minima are –7.8, –5.0, and –2.0 kcal/mol. The consequences of uncertainties in the geometries and limitations in the level of correlation recovery are analyzed. It is suggested that our best estimates will likely underestimate the complete basis set, full CI values by 0.1–0.3 kcal/mol. Vibrational normal modes were determined for all three minima, including an MP2/6–31G** analysis for the most strongly bound complex. Computational strategies for larger phenol–water complexes are discussed. © John Wiley & Sons, Inc.