The electronic structure of small lithium clusters

Ab initio molecular orbital calculations using the STO-6G and STO6-21G basis sets have been performed for the cluster series Li _n ⁺ , Li _n , and Li _n ^– (wheren=2–7). Thirty-two optimized structures are discussed and reported, many of which (especially for the anionic structures) have not yet been considered. The calculations suggest that for all three species the optimum geometries are planar. Of the two levels of theories that were investigated, STO-6G//STO-6G and STO6-21G//STO-6G, the latter hybrid theory was found to be less reliable. In particular, for the anionic structures these calculations should provide a platform from which more sophisticated, i.e., configuration interaction, geometry optimization can be performed.     

Comparison of the use of the MNDO and MINDO/3 methods with ab initio methods to study tautomerism in histamine, 2- and 4-methylhistamine

The semiempirical MNDO and MINDO/3 methods are used to study the various tautomeric forms of histamine, 2-methylhistamine, and 4-methylhistamine. Comparisons of the optimized structures and tautomerization energies are made with values obtained from ab initio Hartree-Fock calculations using the 3-21G and STO-3G basis sets. Based on these results and previous comparisons of STO-3G results with x-ray structures, the present results indicate that while there are some differences in the values of the structural parameters, the changes in structure upon tautomerization and/or protonation are very similar. Further analysis of the MNDO and MINDO/3 structures by means of their utilization in 3-21G and STO-3G calculations indicates that either of these semiempirical methods provides reliable values for the structural parameters. Both methods give good qualitative agreement with the ab initio calculations for the relative energies of the various tautomers in the three compounds. In these studies the MNDO method appears to give better quantitative agreement with the 3-21G and STO-3G results than the MINDO/3 method.     

An SCF and CI Study of the 1,3 Shift in the HX?CHY⇌XCH?YH Isoelectronic Series: X,YCH2, NH,and O

Ab initio molecular orbital calculations at SCF level with the 3-21G, 6-31G, and 6-31G** basis sets and CI level with the 6-31G basis set have been carried out for an isoelectronic series HX? CH?Y and X?CH? YH, where X, Y can be CH₂, NH, and O. Optimized structures (3-21G and 6-31G**) for both tautomers and the 1,3 hydrogen shift transition states are reported. The relative stabilities of the isomers and the barriers of the 1,3 shift are discussed in terms of proton affinities and bond orders. It is shown that both the relative stabilities of the tautomers and the relative barrier heights can be explained qualitatively using simple proton affinity arguments and that the barrier heights are quantitatively related to bond orders.     

Theoretical studies on pteridines. 2. Geometries,tautomer, ionization and reduction energies of substrates and inhibitors of dihydrofolate reductase

Ab initio SCF calculations with STO-3G and 3-21G basis sets are reported for approximately 85 pteridines of interest to the study of the reaction and inhibition mechanisms of dihydrofolate reductase. These include tautomeric, protonated, deprotonated, reduced and 6-substituted forms of the 2-amino-4-oxo- and 2,4-diamino-pteridines, many of which are not easily amenable to experimental investigation. Full geometry optimizations at the SCF/STO-3G level for 30 such pteridines have been performed. A step-wise computational protocol designed to identify the minimum level of theory necessary for reliable prediction of relative tautomer, reduction and ionization energies has been developed in an effort to minimize the cost of calculations for this reasonably large N-heterobicylic system. In general, SCF/STO-3G results were found to be inadequate while SCF/3-21G results obtained with STO-3G optimized geometries agreed with all available experimental evidence with the exception of the relative acidity of 6-methyl-7,8-dihydropterin. Correlation graphs relating experimental pK_a's to the calculated ionization energies are presented: these are of potential predictive value. An analysis is given of the importance of resonance substructures, such as the guanidinium and extended-guanidinium groups, in stabilizing some preferred tautomeric and ionized forms, and in explaining the observed geometry changes.     

Isolobal analogy between trivalent boron and divalent silicon

Singlet organosilylenes with a lone pair and an emptyp orbital are isolobal to trivalent borane if a B-H is equated to the lone pair on Si. Using this analogy, a particular isomer of CSi₂H₂ (24) is predicted to be a stable structure. MNDO calculations on24 and many of its possible isomers suggest that24 is at global minimum on the potential energy surface of CSi₂H₂.Ab initio calculations using a, minimal STO-3G basis set, on some selected structures also support these results.     

Geometry optimization and electronic structures of molecules H_3AXAH_3

The geometries of molecules H_3AXAH_3(X=O,S,Se and A=C,Si)have been optimizedusing STO-3G ab initio calculations and gradient method and the results are in good agreement withreported experimental values.From the STO-3G optimized geometries,we have also calculated theelectronic structures of these molecules using 4-31G and 6-31G basis sets to obtain the MO energies.atomic net charges and dipole moments.The ionization potentials calculated by 6-31G basis set are ingood agreement with experimental values.     

A theoretical approach to substituent effects. Interactions between directly bonded groups in the neutrals X?NH2, X?OH,and X?F and the anions X?NH− and X?O−

Ab initio molecular orbital calculations have been carried out for the neutrals X? NH₂, X? OH, and X? F and the anions X? NH^? and X? O^? with substituents X = Li, BeH, BH₂, CH₃, NH₂, OH, and F. All structures have been fully optimized with the 4-31G basis set which is found to perform considerably better than the minimal STO-3G basis in predicting the lengths of strongly polar bonds. A quantitative analysis of interactions between the directly bonded groups, utilizing energy changes in hydrogenation reactions, is presented and rationalized with the aid of perturbation molecular orbital theory. Favorable interactions occur when electron-donor groups bond to electron-acceptor groups. This applies to both σ and π interactions, the relative importance of which depends on the particular substituents.     

The Electronic Structure of Cubane (C8H8) as Revealed by Photoelectron Spectroscopy

High resolution He (Iα) and He (IIα) photoelectron spectra of cubane are reported. The assignments of the bands to different states of the cubane radical cation are made on the basis of ab initio STO-3G and MINDO/3 calculations, using geometries optimized within each treatment. The vibrational fine-structure observed supports the proposed assignment. An open shell MINDO/3 model for ground state cubane radical cation suggests that the Jahn-Teller distorted system fluctuates between twelve equivalent structures of C_2v-symmetry. Localized molecular orbitals derived from the STO-3G model of cubane indicate that the major feature which discriminates this molecule with respect to other hydrocarbons is the large interaction matrix element between the opposed CC-σ-orbitals of each face.     

A theoretical study on the structure and properties of phenothiazine derivatives and their radical cations

Semiempirical CNDO, AM1, PM3 and ab initio HF/STO-3G, HF/3-21G(d), and HF/6-31(d) methods were employed in the geometry optimization of the phenothiazine and the corresponding radical cation. The results obtained from the PM3 performances were as good as those from the ab initio calculations in the structure optimization of both phenothiazine and phenothiazine radical cation. The PM3 method was used to optimize the structures of a series of N-substituted phenothiazine derivatives and their radical cations. The PM3-optimized results were then analyzed with the ab initio calculation at the 6-311G(d,p) level, which yielded the total energy, frontier molecular orbitals, dipole moments, and charge and spin density distributions of the phenothiazine derivatives and their radical cations.     

Study of hydrogen bonding interactions relevant to biomolecular structure and function

Ab initio molecular orbital calculations were used to study hydrogen bonding interactions and interatomic distances of a number of hydrogen bonded complexes that are germane to biomolecular structure and function. The calculations were carried out at the STO-3G, 3-21G, 6-31G*, and MP2/6-31G* levels (geometries were fully optimized at each level). For anionic species, 6-31 + G* and MP2/6-31 + G* were also used. In some cases, more sophisticated calculations were also carried out. Whenever possible, the corresponding enthalpy, entropy, and free energy of complexation were calculated. The agreement with the limited quantity of experimental data is good. For comparison, we also carried out semiempirical molecular orbital calculations. In general, AM1 and PM3 give lower interaction enthalpies than the best ab initio results. With regard to structural results, AM1 tends to favor bifurcated structures for O? H-O and N? HO types of hydrogen bonds, but not for hydrogen bonds involving O-H? S and S-H? O, where the usual hydrogen bond patterns are observed. Overall, AM1 geometries are in general in poor agreement with ab initio structural results. On the other hand, PM3 gives geometries similar to the ab initio ones. Hence, from the structural point of view PM3 does show some improvement over AM1. Finally, insights into the formation of cyclic or open formate–water hydrogen bonded complexes are presented. © 1992 by John Wiley & Sons, Inc.     

The electronic structure of small sodium clusters

Ab initio molecular orbital calculations using the STO3-21G basis set has been carried out for the cluster series Na _n ⁺ , Na _n , and Na _n ^– (wheren=2–7). The basis set is shown to be reliable compared with more extensive basis sets at the Hartree-Fock level. Thirty-one optimized structures are reported and discussed, many of which (especially for the anions) have not been considered. The STO3-21G//STO3-21G calculations suggest that for most of the species the optimum geometries are planar. In particular, the optimized structures for the anionic species should provide a starting point for more sophisticated configuration interaction calculations.     

The ring-opening of an unsymmetrical tetrahedral intermediate - 2-hydroxy-1,3-oxathiolane

Ab initio calculations were performed on 2-hydroxy-1,3-oxathiolane: and the two products of its breakdown: Complete geometry optimizations were performed at minimal (STO-3G) and split-valence (3-21G) basis set levels. In addition, a single point calculation was performed at 6-31G* level withd orbitals added on sulfur only. The conformation of the oxathiolane intermediate and its stability relative to the breakdown products was investigated. The STO-3G basis set gave an envelope form while 3-21G gave the twist form of the five-membered ring as the most stable. For all three basis sets the ester product was more stable than thioester.     

A conformational analysis of the six isomers of thiobispyridine

The conformational behaviour of the six isomers of thiobispyridine has been investigated using ab initio STO-3G*//rigid-roto, STO-3G*//STO-3G* and 6–31G**//STO-3G* molecular orbital models. The analysis reveals both the importance of optimising critical structure parameters and the basis set dependence of calculated rotational barrier heights. The most reliable model (6–31G**//STO-3G*) clearly indicates that the minimum energy conformers are not planar and that energy barriers between 30–100 kJ mol^?1 restrict inter-conversion to planar structures, thereby preventing conjugation between the p-electrons of the sulfur atom and the π system of both pyridine rings. From the calculated barrier heights, two mechanisms can be employed to explain conformer interconversion about the C? S bond: a disrotatory one-ring flip or a conrotatory two-ring flip mechanism. Where comparisons can be made (eg. 2,2′-thiobispyridine), dipole moment calculations are shown to be in good agreement with experiment. Finally, of the six isomers, appropriately substituted 2,2′, 2,3′- and 2,4′-thiobispyridines are most prone to a Smiles rearrangment.     

A theoretical study of the structure and spectrum of 1,3,5-hexatriene,Part II

The results of ab initio STO-3G and STO-4G calculations are reported on tTt- and tCt-1,3,5-hexatriene, including a STO-3G partial geometry optimization of all C/C bond lengths in both isomers. Results are also reported from a partial INDO geometry optimization on.the tTt-isomer, and complete optimizations on both isomers using empirical consistent force Meld programs. Significant discrepancies are found to exist between these calculated structures and those based on the electron diffraction analysis of Traetteberg. In addition, further results of a CNDO/S spectral analysis are presented and compared with experiment.     

The electronic structure of peracids. Functional models for cytochrome p-450

An extensive set of ground state ab initio and semiempirical molecular orbital calculations has been performed on both peroxytrifluoroacetic and peroxyacetic acids. The equilibrium geometry of peroxyacetic acid was calculated with the STO-3G and MINDO/3 methods, and peroxide rotational barriers for both peracids were obtained with STO-3G and PCILO. Extended basis set calculations with the 6-31G^** basis were performed for both peracids to compare the electronic structures of these two compounds. Electrostatic potential maps in the region of the peroxide bonds of both peracids were also calculated using INDO wavefunetions. These results are discussed with respect to the enhanced reactivity of peroxytrifluoracetic acid relative to peroxyacetic acid and the nature of the oxygen electrophilicity in these compounds and by analogy in cytochrome P-450, for which peroxytrifluoroacetic acid is considered to be an effective chemical model.     

A charge-conserving approximation method for ab initio calculations

A new method is presented for approximate ab initio calculations in quantum chemistry. It is called CCAM (charge conserving approximation method). The calculation method does not include the use of empirical parameters. We use Slater type orbitals as basis set, replacing STO's by STO-2G functions to evaluate three- and four-center integrals and making the STO-2G two-orbital charge distributions have the same total charge as STO. The results are presented for test calculations on five molecules. In view of these results, CCAM is better than ab initio calculations over STO-6G in the results on total energies, kinetic energies and occupied orbital energies. In atomic populations, dipole moments and unoccupied orbital energies, CCAM is also satisfactory. We estimate that CCAM would be as fast as ab initio calculations over STO-2G in evaluating molecular integrals.     

MRD-CI studies of vertical excitation energies of unsaturated hydrocarbon molecules

Systematic MRD-CI calculations using the AM1 Hamiltonian have been carried out for two polyenes and eight aromatic hydrocarbons ranging from benzene to ovalene (C₃₂H₁₄). Twenty singlet–singlet excitation energies in these compounds were calculated and compared with experimental data and ab initio STO-3G results. On an absolute scale, the AM1/MRD-CI approach underestimates the excitation energies to states with dominant covalent character by an average of 1.1 eV, whereas the errors for ionic states are between ?1.0 and 1.0 eV. The STO-3G calculated data are much too high by ≈ 1 eV and ≈ 5 eV, respectively. The inclusion of σπ-correlation effects through second-order Epstein–Nesbet perturbation theory combined with the use of localized orbitals leads to a significant improvement of the ab initio calculated state energies. In an analogous AM1 treatment, negligible corrections for the σπ correlations are found, which is attributed to the implicit account in the parameters and approximation of the semiempirical Hamiltonian. The possible error sources of the calculational methods are discussed. © 1994 by John Wiley & Sons, Inc.     

Semiempirical NDDO calculations with STO-3G and 4-31G basis sets

Possible refinements of semiempirical methods include the use of larger basis sets and of correlated wave functions. These possibilities are investigated in semiempirical NDDO SCF calculations with the STO-3G and 4-31G basis sets, and in correlated calculations at the STO-3G level. The present approach is characterized by the analytical evaluation of all one-center terms and two-electron integrals, and the semiempirical adjustment of the remaining one-electron integrals and the nuclear repulsions. The NDDO SCF results tend to reproduce the correspondingab initio results more closely than experimental data, even if they are parametrized with respect to experiment. The explicit inclusion of electron correlation at the STO-3G level improves the calculated results only slightly.     

Relative energies of substituted benzene valence isomers

Ab initio calculations employing the STO-3G basis set are used to obtain the relative energies of the benzene valence isomers and some selected monosubstituted derivatives. We find that 3,3'-bicyclopropenyl, the least stable of the five (CH)₆ examined, is slightly more stable in the anti conformation than the gauche (Φ = 45°) conformation in agreement with experiment. Substituents are calculated to produce significant changes in the relative energies of the benzene valence isomers. The ground-state isomerization of 1-Dewar benzeneearbinyl cation to benzyl cation is more exothermic than the aromatization of Dewar benzene, but is, in contrast to the latter, symmetry-allowed.