Molecular orbital calculations based on linear combinations of fragment orbitals

A semiempirical MO method based on localized fragment orbitals has been developed, which is particularly suited for the construction of orbital correlation diagrams for the discussion of the electronic structure of complex molecules in terms of fragments and their interactions. The method allows for the inclusion of experimental ionization potentials and electron affinities of the fragments within the calculation of the Fock matrix elements and may thus form the basis of an interpretation of photoelectron spectra, comparable to the interpretation of UV spectra by means of the MIM method of Longuet-Higgins and Murrell. Several levels of approximation are discussed using the acrolein molecule as an example.     

Molecular orbital studies on the conformation of GABA (γ-aminobutyric acid)

In distinction to Extended Hückel Theory which predicts as the most stable conformation of free zwitterionic GABA a totally extended form, PCILO and SCF ab initio studies show that the intrinsically preferred conformation of the isolated molecule is a highly folded one, resulting from strong interactions between the two charged ends. Computations are also carried out for hydrated GABA in the supermolecule approach allowing moreover for the flexibility of binding of some of the water molecules of the first hydration shell. They predict the coexistence in solution of a large number of conformations showing different degrees of folding (or extension), a result confirmed by recent NMR studies. This and a number of similar results show that we have to adapt our thinking on the role of conformations in pharmacological activity to this situation, which was frequently obscured by the more abundant results of X-ray crystallography yielding a single conformation.     

A new approach to variational CI calculations on systems containing non-interacting molecules

A new method to perform variational CI calculations on systems containing non-interacting molecules is presented.     

On the use of pseudopotentials in molecular calculations

A series of tests was performed of the Kahn-Goddard-Melius-Topiol pseudopotentials in view of their utilization with small contracted basis sets in molecular computations. The effects of inner-shell separability and of basis set contraction are underlined. The utilizability of Topiol's valence least-squares fitted Gaussian basis sets is studied.     

A non-empirical molecular orbital study on the relative stabilities of adenine and guanine tautomers

The relative stabilities of a series of adenine and guanine tautomers have been calculated using anab initio Hartree-Fock-Roothaan SCF MO method. The calculated relative stabilities agree in general with the results of earlier semiempirical studies. According to the present study, tautomeric forms with regular Kekulé structure for the six-membered purine ring are the most stable. The amine-imine tautomerization of purine bases is not likely to be responsible for spontaneous mutations in DNA.     

Topological properties of benzenoid systems

The sextet polynomial of non-branched cata-condensed benzenoid molecules is proved to be related (Eq. (3)) to the characteristic polynomial of a tree.     

Comments on the electronic properties of cyclopropane,cyclobutane, and cyclohexane

Spectral quantities of cyclopropane, cyclobutane, cyclohexane, and of several derivatives, have been calculated by a semiempirical all-valence electron SCF-CI MO method. In cyclopropane, HOMO is practically localized in the carbon-frame, and LVMO is purely so. In cyclobutane, these two MO's are based on C-H bonds, while cyclohexane holds an intermediate position. Despite the overall similarity-experimental and computed-of the spectra of these molecules, assignments are non-parallel. Like cyclopropane, cyclobutane can extend conjugation, but to a diminished degree; cyclohexane behaves in this respect like an acyclic alkane. An interpretation of this gradation, in terms of the nature of high-lying MO's, is proposed.     

Non-empirical molecular orbital theory of the electronic structure of molecular crystals

A non-empirical molecular orbital treatment of molecular crystals, based on SCF perturbation theory and matrix partitioning methods is presented.     

Non-empirical calculations and basis set effects on the inversion geometry of the aniline molecule

The geometry of the amine group and the barrier to internal conversion in aniline have been studied by single-determinantab initio SCF calculations using several basis sets from minimal to double-zeta quality. The results obtained from different types and sizes of basis sets are discussed. Calculations performed with the two most flexible basis sets yield inversion barriers of 0.9–1.1 kcal/mole and angles of pyramidalization at the nitrogen atom of 38–39 ° which are in good agreement with the experimental results (1.3 kcal/ mole and 38 °). Orbital and overlap population analyses are performed and compared with the expected mesomeric and inductive effect. The calculated dipole moment 1.48–1.49 D also agrees with the experimental values (1.48–1.50 D).Dedicated to Professor O. E. Polansky, Mülheim/Ruhr, on the occasion of his 60th birthday.     

Structure and properties of non-classical polymers II. Band structure and spin densities

Theoretical investigation of the band structure of three types of nonclassical polymers, namely alternant (one- and two-dimensional), nonalternant and heteroatomic, are carried out. Although polyradicals, these polymers have a considerable delocalization energy which may determine their relative stability.The spin-density distribution of the alternant type of non-classical polymers corresponds to a ferrimagnetic ground state at 0 K.The non-classical polymers represent a new class of organic systems as their band structure and magnetic properties essentially differ from those of common polymers.     

Ab initio crystal orbital studies on linear chains of H atoms

Anab initio crystal orbital method is used to calculate the energies of an infinite chain of H atoms and of linear arrangements of H₂ molecules with different interatomic distances. The H₂ arrangements are not stable in respect to isolated molecules. The cohesive energy of an optimized arrangement of H atoms chain is 0.0354 a.u.     

Transferability of some properties of localized molecular orbitals

The transferability of the first and second order electric moments of localized orbitals have been shown. Some moment characteristics have also been investigated. The standard deviation of first order localized moments is less than 2.0%. The standard deviation of second order electric moment components are larger: although they do not exceed 4.0% for bond orbitals, for lone pair ones they vary about 1.5–10.1%.     

Extended Hückel calculations on the chemisorption of acetylene on tungsten

The chemisorption of acetylene on various faces of bcc tungsten, modelled by clusters of various sizes, has been studied by the Extended Hückel method. Relative binding energies on the various sites have been obtained, and the bonding modes of acetylene are discussed, also in comparison with experiment. It is concluded that, whenever possible, acetylene is di- bonded to the surface, and sp ²hybridization is considered likely in all cases.This paper is dedicated to Professor Hermann Hartmann on the occasion of his 65th birthday.     

Complementary molecular orbital investigations on the conformation of choline derivatives

Quantum-mechanical computations by the PCILO method, applied previously to the study of the conformational properties of acetylcholine and its derivatives modified in the central part of this molecule, are extended to modifications involving its cationic head and its ester terminal. The replacement of the methyl groups of the cationic head by hydrogens or ethyl groups leads to a steep decline in parasympathomimetic activity. It is shown that the triethyl derivative conserves the gauche form as the most stable one. The redistribution of the electronic charges at the onium group implies, however, a transition from an ionic to a hydrophobic binding. The replacement of the methyls by two or three hydrogens leads to a different preferred gauche-gauche conformation. The replacement of the methyl group at the ester terminal by a phenyl ring enables a comparison with the conformational properties of local anesthetics. The study brings about evidence, substantiated by NMR spectroscopy, that acetylcholine analogs and protonated local anesthetics are conformationally similar. Choline ethers also show a general preference for a gauche conformation. Nevertheless, biological studies do not indicate a constant correlation between conformation and biological potency. Conformational analogies or discrepancies alone cannot thus account for the fine details of the biological activity which must depend also on the electronic structure.This work was supported by the A.T.P. N A 655-2303 of the C.N.R.S.     

Theoretical studies on the conformation of saccharides

Conformations of 2-methoxytetrahydropyran as a model for the six-membered ring in aldopyranosides have been calculated by the PCILO method using the algorithm of the conjugated gradient to optimize the geometry. The calculated geometry of the fourteen basic forms of 2-methoxytetrahydropyran was found to be in agreement with the available data obtained by X-ray diffraction of pyranosides. The results indicate differences in the geometry of 2-methoxytetrahydropyran resulting from the change of the axial vs. equatorial position of the methoxyl group. These changes are particularly meaningful in the values of bond angles and they are in agreement with the anomeric and exoanomeric effects. The experimentally found differences in the energies of an axial (⁴ C ₁) and equatorial (¹ C ₄) conformer, G = 2.9–3.7 kJ/mol, and the dipole moment, = 1.20 ± 0.05 D (1D = 3.33 10^–30mAs) agree well with the calculated values E = 3.18 kJ/mol and <> = 1.18 D which, in turn, suggest that the axial conformer is preferred over the equatorial one by a ratio a:e = 78:22.     

A theoretical study of electrophilic reactions

The electronic structures of protonated formyl and acetylium cations and their deprotonation paths leading to HCO⁺, COH⁺ and CH₃CO⁺have been studied by means of ab initio calculations. The results support Olah's theory that dipositive species can be the de facto reagents in electrophilic reactions.     

Quantum mechanical studies of environmental effects on biomolecules

LCAO SCF calculations with a (7s, 3p/3s) atomic basis are performed on urea and hydrated urea in order to determine the preferred sites of water fixation using the supermolecule approach. It appears that up to five water molecules can be directly bound to urea in the first hydration shell.     

Screened potential in π-electron systems and its application to the estimation of excitation energies

A diagrammatic technique was developed for the estimation of the screened potential of -electron systems. The screened potential was expanded in terms of the polarization propagators which were constructed from either the singlet, , or triplet vertex part, . These vertex parts correspond to the singlet or triplet excitations, respectively, in the Random Phase Approximation (RPA) containing exchange diagrams. The excitation energies were calculated by using the screened potential in the framework of RPA with exchange. The excitation energies of several conjugated molecules with or without a hetero atom are in agreement with the experimental data.     

On the configuration interaction method for singlet states

The direct CI method, which avoids explicit calculation of the Hamiltonian matrix, is presented in a new form. The method is linked with Davidson's algorithm for iterative evaluation of the ground state eigenvector. The viability of the method is indicated by the test calculations on water which are described.