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.     

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.     

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.     

Multiplicity of the ground state of large alternant organic molecules with conjugated bonds

The multiplicity and the full spin of the ground state of large alternate molecules with conjugated bonds are considered. It is strictly shown that if the numbers of starred and unstarred atoms (say, carbon) differ from each other the full spin of the molecule is more than zero. Some possible planar and linear molecules having the full spin to be proportional to their sizes are presented. Particularly, they would be ferromagnets at infinite sizes.     

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.     

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.     

Band structure of even and ions of odd polyenes

The analysis of experimental data for singlet transitions (E _n) of even polyenes (I), cations (II) and anions (III) of odd polyenes show that for infinite chains E (I)/E (II)=E (I)/E (III) = 2:1. It is shown that the energy gap is equal for the three systems. In cases (II) and (III) there is a level (NBMO) in the gap which is vacant in (II) and occupied in (III). That is why the first optical transition in (II) and (III) depends on the semiwidth of the gap.     

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.     

Band structure of nonclassical polymers

It is shown that fully conjugated alternant nonclassical polymers are comparatively stable systems as a result of the considerable delocalisation energy.The energy characteristics of the polymers are strongly determined by the electron spin distribution of the degenerate nonbonding molecular orbitals; the full spin configuration,S>0, is favoured. The spin densities depend on the electron correlation and alternate; this corresponds to a ferrimagnetic state of the polymer at 0°K.     

Configuration interaction methods for improving unrestricted Hartree-Fock spin densities

Limited Configuration Interaction wave functions based on Unrestricted Hartree-Fock natural orbitals are found to be easy to compute and to give much more satisfactory spin densities than are provided by techniques currently in use.     

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.     

Clusters of hydrogen molecules: Ab initio SCF calculations corrected semiempirically for correlation energies

Stabilization energy of the (H₂) _n clusters (n = 2–8) was calculated as a sum of the SCF interaction energy and the semiempirical interaction correlation energy estimated according to Sinanolu and Pamuk. Optimum successive attachment of hydrogen molecules leads to the formation of a gas-phase solvation shell consisting of seven hydrogen molecules. Basis set effect has been found to be important with all clusters under study. The non-additivity effect was investigated with the (H₂)₄ cluster. Vertical ionization potentials of the clusters considered are predicted to be 0.4–0.6 eV lower than the ionization potential of the parent H₂ molecule.     

Etude théorique de la conformation des complexes moléculaires de type donneur accepteur d'électrons

The stabilisation energy of donor acceptor molecular complexes given by tetracyanoethylene with thiazole and phenylthiazoles have been calculated by several semiempirical approaches. The results obtained with the Extended Hückel theory, the CNDO/2 method and the PCILO method, the latter a perturbative one, are compared and this comparison shows that the PCILO method is well adapted to the study of molecular complexes. The different contributions to the stabilisation energy and the origin of the intermolecular forces have been discussed in relation to the geometry of the complex. For all the studied complexes, the calculated results obtained by the PCILO method agree very well with experiments.

     

Numerical integration of overlap and ligand contributions to the electric field gradient

The total electric field gradient (EFG) tensor V _pq is calculated by numerical integration of threedimensional integrals. Each of them is solved a) by integrating over one dimension analytically and b) by integrating over the remaining two dimensions on the basis of a Gauss-type integration rule. The use of 100 abscissas in the twodimensional numerical integration scheme yields satisfactory accuracy which was checked by evaluating overlap integrals; an increase to 400 abscissas does not increase the result drastically. Calculating quadrupole splittings E _Q from numerically integrated electric field gradient tensors V _pq we observe that depending a) on the amount of covalency and b) on the amount of deviation from octahedral or tetrahedral symmetry, involved in a molecular system, overlap and ligand contributions to V _pq play an important role. Especially for the sandwich compound ferrocene, Fe(C₅H₅)₂, we find a significant difference between E _Q ^{num. int.} which follows from the numerical integration method, and E _Q ^conventional which is derived from effective charges.Supported by Deutsche Forschungsgemeinschaft     

On the affinity of cytosine towards electrophiles

Ab initio SCF computations on the intrinsic preferences of the H⁺, CH ₃ ⁺ and C₂H ₅ ⁺ cations towards the two principal sites of protonation or alkylation on cytosine, N₃ or O₂, show that this preference undergoes a continuous modification with the increase in size and complexity of the cation. N₃ is the preferred site of fixation of H⁺, O₂ the preferred site of C₂H ₅ ⁺ , while CH ₃ ⁺ has no marked preference. The exchange repulsion term of the binding energy appears responsible for the preference of C₂H ₅ ⁺ for O₂.This work was supported by the Ligue Francaise contre le Cancer and the National Foundation for Cancer Research (USA)     

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.     

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.     

Calculations of the electronic structures of cage molecules using free-electron orbitals as a basis

A non-empirical molecular orbital method, particularly suitable for calculations on cage-like molecules, is described. The method uses as basis functions the set of free-electron functions which are the solutions of Schrödinger's equation for an electron confined between two concentric, spherical potential energy barriers. Application of the theory to the SCF calculation of the energies of the delocalized electrons in benzene and tetrasulphur tetranitride shows that the model is capable of interpreting the properties of such systems. However, it does highlight a difficulty in the calculation of excited state energies with one-centre models which appears to be largely unrecognized.Extension of the method to a consideration of all the valence electrons, using P₄ as an example, reveals problems the origin of which is an inadequate treatment of the core electrons. It is suggested that these problems may best be dealt with by use of a suitable pseudo potential.     

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.