A new approach to theab initio energy of the homodesmic reaction for the resonance energy of benzene

Summary A scheme of the basis set superposition error (BSSE) correction is first proposed and introduced to determine theab initio energy of the homodesmic reaction for the resonance energy of benzene. Calculations with 6-31G*(5D) and 6-31G*(6D) basis sets at the complete fourth-order Møller-Plesset perturbation level furnish the energy value of 21.35 kcal/mol after the correction, which is in complete agreement with the experimental value of 21.3±0.2 kcal/mol. The energy values at the lower theoretical levels are generally underestimated but they are superior to the uncorrected values. The inclusion of triple excitations displays the dominant contribution of the correlation energy. Detailed analysis of the results reveals some of the similarities between the homodesmic reaction of benzene and the interaction of van der Waals molecule, which provides further justification of the BSSE correction scheme presented in this study.     

A comparative ab initio study of ethylene,acetylene and benzene

Ground state properties have been calculated by use of a medium-sized Gaussian basis set and comparison with other bases has been made. Contraction to double-zeta of a comparatively small basis is found to be superior to a large set of primitive Gaussians contracted to minimal basis. Molecular optimization is not important for double-zeta bases. Inclusion of a balanced set of polarization functions is essential in all cases studied. Population analysis gives a certain insight in molecular properties but contour maps are found to be significantly superior. This is demonstrated on bonding properties of corresponding orbitals within the series. In case of benzene Slater's energyband plot is shown to be useful for classifying bonding properties.     

Systematic regularities of molecular SCF energies

New formulae for the approximate computation of total molecular energies are developed based on ab initio calculations of n-alkanes. Their application to various kinds of molecules reveals that good expectation values for total molecular energies can be obtained by considering only the one-electron terms h _i and the nuclear repulsion energy. It is further shown that very good agreement with SCF total energies is obtained by a relationship which connects the total energy with the sum of inner-shell (core) orbital energies. The results turn out to be better than those obtained using Ruedenberg's approximation, which takes both inner-shell and valence-shell orbital energies into account.     

The matching of structural elements in reactions for evaluating stabilization energies for benzene and monosilabenzene

Using the 6-31G^* basis set the stabilization energy for monosilabenzene has been calculated as E _T for the homodesmotic reactions in which (i) ethylene is the other reactant with trans 1,3-butadiene and trans 2-sila-1,3-butadiene as the products, and (ii) trans 1,3-butadiene is the other reactant with trans 1,3,5-hexatriene and trans 3-sila-1,3,5-hexatriene as the products. The values are 17.2 and 17.6 kcal mol^–1 respectively, compared to 24.8 and 23.5 kcal mol^–1 respectively for the corresponding reactions of benzene. The difference is E _T for the interchange of =SiH- and =CH- groups between the aromatic and aliphatic structural environments, and is thus more informative than the quotient. The extent of the matching of structural elements in these reactions, and in the benzene reactions using longer chain polyenes, is assessed in terms of the structural environment about each heavy atom in its entirety — namely the nearest neighbor atoms (NN), the next nearest neighbor atoms (NNN), and the more distant atoms (NNNN and NNNNN). Matching by number alone is designated isoplesiotic, and if by both number and kind, homoplesiotic.     

Ab initio program for treatment of related molecules. I. Integral part

Anab initio integral program is described. It utilizes the local symmetries to avoid the redundant computation of integrals over spatially equivalent subsets of the basis. The integrals are grouped in a particular way to facilitate their transfer. The program is very suitable for the treatment of related systems with model geometries. The computing times of different programs are compared and the efficiency of the presented one is demonstrated.     

Ab initio calculation of the band structure of some boron polymers

Ab initio calculations using a STO-3G basis set have been performed on the polymer systems (HBX) _n where X = Be, BH, CH₂, NH, and O. Energy band diagrams and accompanying density-of-states plots have been obtained. The highest filled orbital of (HBNH) _n and (HBO) _n occurs at the X-point and possesses character while the - framework orbital at the X- point is the highest filled level for (HBBe) _n , (HBBH) _n , and (HBCH₂) _n . The conduction band for all five species has symmetry and the band gap of the (HBX) _n species increases in the order X = Be < BH < NH < O < CH₂. An estimate of the energy of polymerisation of the (HBX) _n systems suggests that HBNH is particularly stabilised by polymerisation. The electron distribution in (HBBe) _n shows a - electron drift towards boron, while in the other four systems the net electron transfer is directed away from boron. There is significant electron back-donation to boron in (HBO) _n and (HBNH) _n .     

Ab initio molecular orbital calculations of the infrared spectra of hydrogen-bonded complexes of water,ammonia and hydroxylamine

Summary The geometries of three hydrogen-bonded dimers of hydroxylamine have been optimized, at the MP2 level of theory, using the 6-31G** basis set. These calculations yielded three separate local minima on the dimer potential energy surface. The interaction energies of these three species have been calculated, and corrected for basis set superposition error. The infrared band wavenumbers and intensities have been computed, and the monomer-dimer wavenumber shifts and intensity enhancements rationalized in terms of the types and strengths of hydrogen bonds present. The predicted wavenumbers have been correlated with those measured in a recent matrix isolation spectroscopic study, and an argument for the structure of the preferred dimer has been presented.     

Dissociation Energies and Rotational Barriers About CC Single,Double, and Triple Bonds: A Hybrid HF-DFT Approach (Becke3LYP/6-311++G**)

Structural Chemistry -     

Ab initio study of oxaphosphaalkyne P=C-OH and its isomers

Anab initio study on HP=C=O PC-OH isomerization has been performed at the 3-21G(d(P)), 4–31G(d(P)) and 6–31G** levels. In contrast to analogous transformation of 1-aza-3-phosphaallene, this process is endothermic in the gas phase. It proceedsvia a [1,3]H shift without an intermediate. The route, including two subsequent [1,2]H shifts, leads to phosphinidene.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 1, pp. 172–173, January, 1994.     

DFT and ab initio quantum chemical studies on p-cyanobenzoic acid

The Fourier transform infrared (FTIR) and FT-Raman spectra of p-cyanobenzoic acid (CBA) have been recorded in the range 4000-400 and 4000-100 cm(-1), respectively. The complete vibrational assignment and analysis of the fundamental modes of the compound were carried out using the observed FTIR and FT-Raman data. The vibrational frequencies determined experimentally were compared with theoretical wavenumbers obtained from ab initio HF and DFT-B3LYP gradient calculations employing 6-31G**, 6-311++G** and cc-pVTZ basis sets for the optimised geometry of the compound. The geometry and normal modes of vibration obtained from the HF and DFT methods are in good agreement with the experimental data. The normal coordinate analysis was also carried out with ab initio force fields utilising Wilson's FG matrix method. The interactions of cyano and carboxylic acid groups with the skeletal vibrational modes were investigated.     

Problems and prospects in the ab initio treatment of pure and defective crystals

Accurate Hartree-Fock LCAO calculations for moderately complex crystalline systems are now feasible; a number of important applications may be envisaged in the areas of material science and technology. Some critical aspects of the corresponding computer schemes are discussed which are of fundamental importance in determining the cost of the computation. Data are provided concerning actual computations which are indicative of the kind of periodic systems that can (or cannot) be treated at present. The result of a perfect-crystal ab initio HF study can be used as an input for treating with the same approximation local-defect problems, by use of suitable embedding techniques. A scheme of this kind is presented, and its computational implications are discussed: due to the intrinsic complexity of this problem, it may be foreseen that the study of defects in crystals will be a typical application of supercomputers in the area of quantum chemistry.This paper was presented at the International Conference on The Impact of Supercomputers on Chemistry, held at the University of London, London, UK, 13–16 April 1987     

Hydration of the carboxylate group: Anab initio molecular orbital study of acetate-water complexes

The hydration of the carboxylate group in the acetate anion has been investigated by performingab initio molecular orbital calculations on selected conformers of complexes with the form CH₃CO₂ ^– ·nH₂O·mH₂O, wheren andm denote the number of water molecules in the first and second hydration spheres around the carboxylate group, andn + m 7. The results of RHF/6–31G^* optimizations for all the complexes and MP2/6–31+G^** optimizations for several one-water complexes are reported. The primary consequence of hydration on the structure of the acetate anion is a decrease in the length of the C-C bond. Enthalpy and free energy changes calculated at the MP2/6–31+G^** and MP2/6–311+ +G^** levels are reported for the reactions CH₃CO₂ ^– + [H₂O] _P CH₃CO₂ ^– ·nH₂O ·mH₂O where [H₂O] _P is a water cluster containingp water molecules andp=n+m 7. The calculations show that conformers with the lowest enthalpy change on complex formation are often not those with the lowest free energy change, due to a greater entropic loss in complexes with tighter and more favorable enthalpic interactions. Hydrogen bonding of six water molecules directly to the carboxylate group in CH₃CO₂ ^– is found to account for approximately 40% of the enthalpy change and 37% of the free energy change associated with bulk solvation.     

Ab initio calculations of electrostatic potentials and deformation densities for a series of choline ester model systems

Ab initio LCAO-MO-SCF calculations using a double zeta basis set have been performed for the methyl esters of acetic acid, carbamic acid, methylcarbonic acid, and trifluoracetic acid, in order to model the corresponding choline esters. The systems have been compared by means of population analyses, electron density differences, electrostatic potentials and potential differences. The significance of the electrostatic potential in connection with crystal structure and packing has been studied. The differences in the proton affinity of the compounds have been correlated to differences in the potentials.     

Ab initio study of the conformational equilibrium of ethylene glycol

Summary The conformational equilibrium of ethylene glycol (CH₂OHCH₂OH) has been examined by performing geometry optimizations at the 6-31G*, MP2/6-31G* and 6-31G** levels. Final energies have been calculated at the MP3 level with the optimized geometries. The two most stable conformers are atGg andgGg but it is verified that the inclusion of electronic correlations reduces their energy difference of 0.6 kcal/mol at the HF level to less than 0.2 kcal/mol. The possible coexistence of the two most stable conformers is in agreement with some previous studies of Frei et al. For thetXg conformer a detailed analysis of the intramolecular potential as a function of rotation around the C-C bond is also reported.     

Ab initio MO studies of nuclear spin-spin coupling constants in CH4, SiH4, AlH 4 − and GeH4 systems

Summary Ab initio molecular orbital calculations of electron coupled nuclear spin-spin coupling constants are performed for CH₄, SiH₄, AlH ₄ ^– and GeH₄ systems using the SCF perturbation theory. Basis set dependence of the major contributing terms such as orbital diamagnetic, orbital paramagnetic, spin dipolar and Fermi contact terms are studied. The study also illustrates the relative importance of bond centred functions and nuclear centred polarization functions in predicting the directly bonded and geminal couplings in the systems selected. Basis sets having uncontracted cores functions and augmented with bond functions seem to predict most of these couplings fairly satisfactorily when compared to the experimental values.     

Gradient optimization of polarization exponents inab initio MO calculations on H2SO → HSOH and CH3SH → CH2SH2

Summary Analytical gradients were used to optimize the polarization function exponents in the 6-31G(d) and 6-31G(d, p) basis sets for the reactants, transition structures and products in the reactions H₂SO HSOH and CH₃SH CH₂SH₂. The optimizedd exponents on the heavy atoms change by ±10% in the course of the reactions and depend on the bonding of the heavy atoms. Thep exponents on the hydrogens change by as much as a factor of 5 and depend on the element to which the hydrogen is bonded and its valency. The effect of exponent optimization on the relative energies is small (±3 kcal/mol). With the 6-31G(d, p) basis set, optimization of the polarization exponents can make some of the bonds significantly more polar, as judged by the Mulliken charges.     

Ab initio calculations on large molecules using molecular fragments. Characterization of the zwitterion of glycine

The ab initio molecular fragment approach is applied to a characterization study of the ground state of the zwitterion of glycine. Included among the properties studied are the — conformational energy surface, the electronic structure, and the magnitude and direction of the dipole moment. The results of the present study are compared to the results of other theoretical and experimental studies.This work was supported in part by the National Science Foundation, the University of Kansas, and the Upjohn Company, Kalamazoo, Michigan 49001.NSF Trainee (1969-1972).Alfred P. Sloan Research Fellow (1971–1973).     

New molecular species of potential interest to interstellar chemistry: A theoretical study of MgSiN, MgNSi and related species

An ab initio study has been performed to characterize the probable magnesium containing interstellar species MgSiN, MgNSi and their ionized, hydrogenated and protonated forms. We are able to locate four protonated and four hydrogenated magnesium species with planar geometry at MP2(Full)/cc-pVTZ level of theory. MgNSi is found to be more stable than MgSiN and their connecting transition state is also located. The ionization potential for both MgSiN and MgNSi are small, 7.91 eV and 7.01 eV, respectively. All possible protonation sites are considered for these two species but the preferred protonation sites are found to be silicon for MgNSi and nitrogen for MgSiN. Enthalpies of formation at 0 K (Δ_fH°) and bond dissociation energies (D^o(X–Y)) are computed for all the species at G3 and G3MP2 level of theory. Finally, the reaction enthalpies for ion–molecule processes are calculated and most of the processes are found to be exothermic and hence thermodynamically favorable in interstellar region.     

Desolvation effects on the dissociation energy of diatomic molecules:Ab initio study of the dissociation of Li-F in polar media

Summary The potential curve of the ground state dissociation of Li-F in water has been studied by a combination of a standardab initio Hartree-Fock procedure and a perturbative reaction field approach. The electrostatic solute-solvent interaction is accounted for by the generalized Born formalism introduced through a perturbation approach. The calculations were carried out at a 6–311+G* basis set level. Diffuse functions ofs symmetry were included to model a desolvation potential. A double well potential curve was obtained for the dissociation of this molecule in the presence of a highly polarizable medium. The first minimum, corresponding to an ion pair, electrostatically bound, is found at aR(Li-F)<6.0 Å distance. As the two ions come together, a desolvation barrier of about 30 kcal/mol is to be overcome before the formation of the neutral Li-F at 1.56 Å. The barrier to ionization towards the ion pair is calculated to be about 14 kcal/mol. The dissociation of the ion pair towards the free ions is discussed in terms of the electrostatic solvation entropy changes.Contribution No 6 from Centro de Mecánica Cuantica Aplicada (CMCA)     

Theoretical studies of molecular conformation. Derivation of an additive procedure for the computation of intramolecular interaction energies. Comparison withab initio SCF computations

An additive procedure (SIBFA) is developed for the rapid computation of conformational energy variations in very large molecules. The macromolecule is built out of constitutive molecular fragments and the intramolecular energy is computed as a sum of interaction energies between the fragments. The electrostatic and the polarization components are calculated using multicenter multipole expansions of theab initio SCF electron density of the fragments. The repulsion component is obtained as a sum of bond and lone pair interactions.Tests of the procedure on a series of model compounds containing ether oxygens and pyridine-like nitrogens are reported and compared with the results of correspondingab initio SCF calculations. The resulting methodology is compatible with the simultaneous computation of intermolecular interactions.