Delocalization corrections to the strictly localized molecular orbitals: A linearized SCF approximation

An explicit formula is derived for calculating the delocalization corrections (tails) to be added to the strictly localized bond orbitals. It was obtained by solving analytically the SCF problem for the interbond interactions in a linearized approximation. The model calculations at the CNDO/2 level show that this simple approach is sufficient to account for the molecular conformations.     

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.     

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.     

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.     

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.     

An efficient ab initio gradient program

An ab initio Hartree-Fock gradient program is described. It is characterized by (1) efficiency of the gradient evaluation, and (2) capability of handling higher angular momentum (d andf) basis functions. The latter are constructed from shifted Cartesian Gaussian p-type primitives. A satisfactory solution is presented for the problems connected with the neglect of small integrals in a gradient program. Methods for increasing the efficiency of the SCF procedure are discussed.     

Chemical applications of group theory and topology

The following procedure is described for investigating the qualitative dynamics of simple chemical systems: 1) A so-called influence diagram is generated representing the relationships between the reference reactants (phase-determining intermediates); 2) This influence diagram is used to generate a truth table indicating possible transitions between state vectors representing the signs of the time derivatives of of the reference reactant concentrations; 3) The truth table is used to determine a state transition diagram representing the flow topology around unstable equilibrium points; 4) The characteristic equation of the adjacency matrix of the influence diagram is solved in order to determine the presence of such unstable equilibrium points. The two types of qualitative dynamics possible for chemical systems containing two reference reactants and one feedback circuit are bifurcation between two attracting regions (bistability) and limit cycle oscillation. However, in two reference reactant systems oscillation requires an additional self-activating loop to generate the unstable equilibrium point required for its realization. Bistability and limit cycle oscillation are also two of the possible types of qualitative dynamics for chemical systems containing three reference reactants. However, chemical systems with three reference reactants and two or more feedback circuits can also contain interlocking limit cycles, which can lead to toroidal oscillations or chaos. The influence diagrams are given for the systems exhibiting these various types of dynamic behavior along with a summary of the important properties of all 729 possible influences for simple chemical systems containing three reference reactants.     

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.     

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.     

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.

     

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.     

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.     

Ab initio SCF and CI study of the electronic spectrum of pyridine N-oxide

Configuration interaction (CI) studies of ground, n *, * * electronically excited states are reported for pyridine N-oxide. The transition energy to the lowest * excited ¹ B ₂ state is calculated at 4.35 eV, compared to the experimental spectrum range of 3.67–4.0 eV. This state lies below the lowest n * excited ¹ A ₂ state calculated at 4.81 eV above the ground state. The only experimentally reported triplet state at 2.92 eV above the ground state is predicted to be the ³ A ₁ (*) state. The calculated energy lies at 3.27 eV. Numerous other high-lying singlet states as well as the triplet states have also been calculated. The intramolecular charge transfer character of the ground and the excited states have been studied in terms of the calculated dipole moment and other physical properties.     

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.     

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.     

Ab initio calculation of some low-lying electronic excited states of methane

The energies of some low-lying electronic excited states of methane are calculated by using wave functions built up in terms of plane waves modulated by multicenter Gaussian factors. The wave functions of the various states are evaluated by a two steps iterative process. In the first step, each excited orbital is determined while keeping all other rigid; in the second, rearrangement effects are introduced. Final results are in good agreement with experimental data and allow to enhance an assignement hypothesis for the first electronic transitions.     

Dissociation of diimide

Parts of the potential energy surface of N₂H₂ have been studied using CASSCF- and contracted CI-methods. Of particular interest was the concerted dissociation of cis- and trans-diimide into N₂ and H₂, since the trans-dissociation is symmetry allowed and the cis-dissociation forbidden. Three different saddle points were located, of which only one, of C ₂- symmetry, is a true transition state. Elaborate numerical gradient methods using exact Hessians and update procedures had to be used to find these saddle points on the unexpectedly complex N₂H₂-surface. The barrier height with respect to trans-diimide is 61 kcal/mol after vibration correction. Since this energy is higher than the barrier for interconversion, cis- and trans-diimide have the same transition state. It is further found that diimide preferably dissociates stepwise, by losing one hydrogen at a time, rather than in a concerted way. This conclusion is drawn basically because the geometry of the transition state for the concerted dissociation has a very long H-H distance of 5.6 a.u. The N-H bond energy in trans-diimide is 56 kcal/mol after vibration correction.     

On the mathematical formulation and interpretation of LACO MO theory

The mathematical basis of LCAO MO theory is studied, both within the Hartree-Fock approximation and in more exact formulations. The basic LCAO expansion for molecular orbitals ¦> in terms of atomic orbitals ¦x> is conveniently written ¦> = ¦x> S ^–1 B where S is the overlap matrix for atomic orbitals and B is the matrix of atomic orbital-molecular orbital overlaps. It is suggested that matrices P and Q, defined by P=B B and Q=BnB where n is the matrix of molecular orbital occupation numbers, are appropriate to the interpretation of molecular calculations in terms of atomic orbital components, electronic populations and the degree of bonding. Implications for Hartree-Fock calculations are investigated.