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.     

SCF MO CI perturbation calculations of inner shell and valence shell vertical ionization potentials

A CI method for calculating inner and valence shell vertical ionization potentials is presented. It is based on ab initio SCF MO calculations for the neutral closedshell ground state followed by CI perturbation calculations for the ground and ion states including all spin and symmetry adapted singly and doubly excited configurations with respect to the main configurations of the state of interest. The state energy is computed by performing a CI calculation for a set of selected configurations, and then adding the contributions of the remaining configurations as estimated by second order Brillouin-Wigner perturbation theory. The use of the same set of MO's for all states together with the CI perturbation method makes the method rather rapid. The numerical results are, in spite of the limited Gaussian basis sets used, in good agreement with experiment.     

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.     

On the evaluation of CI matrix elements for a canonically ordered basis

Using the unitary group approach it is shown that the amount of storage needed for the construction of symbolic CI matrix element lists for N-electron basis functions with large numbers of open shells and arbitrary multiplicities may substantially be reduced compared to methods currently available in the literature.     

Gauss' theorem approach to delta function terms in solid spherical harmonic expansions

The nature of the Dirac delta-function singularity in the expansion theorem for an irregular solid spherical harmonic about another centre is discussed for the casel=2. An alternative derivation, motivated by Hobson's derivative expression for solid spherical harmonics and utilizing Gauss' Divergence Theorem, is presented. The orientation dependence is then simply derived from the rotational properties of spherical harmonics.     

Iterative natural orbitals for configuration interaction using perturbation theory

Approximate natural orbitals are determined iteratively from CI expansions constructed using first-order perturbation theory in order to investigate the possibility of eliminating the complete transformation of MO integrals on each iteration. Results on LiH and H₂O are compared with fully variationally determined NO's to assess questions of convergence.     

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.     

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 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.     

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.     

Acidity scale for strong Bronsted acids referred to their complexes with water and in the Gaseous phase

CNDO/2 calculations on AH...B complexes in the gaseous phase (AH = strong acids, B = water) gives the following order H₂O < HF < HCl < HSO₃CF₃ < H₂SO₄ < HSO₃F < HClO₄ < H₃O⁺ if we take as reference the Wiberg residual indexp(H... O). Other properties and in particular the charge transfer behave in the same fashion. Charge variations for each atom give an explanation for the Pauling rule for the oxygenated acids.

     

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.     

Theoretical study of the potential energy curves of the series of diatomic radicals MeIIX

The potential energy curves of some low lying electronic states of the diatomic radicals BeCl and MgF have been calculated. The calculation has been performed according to a stepwise procedure, outlined in a previous work. The potential energy curves are very similar to those of the mercury halide radicals, the electronic transitions of which can be employed for efficient laser apparatus.     

A quantum-mechanical study of the interaction of methylglyoxal with guanine

Ab initio self-consistent field molecular orbital computations on the relative stabilities of the different possible intermediate adducts for the reactions between methylglyoxal and guanine, as well as the evaluation of the relative stabilities of the two different possible final cyclic products (IIIb and IIIc) point all to the conclusion that it is the addition product in which the methyl group is close to the amino nitrogen which is the most stable one.     

Conformational analysis of dimethylphosphate with quantum-mechanical and classical methods

In order to shed light on the conformational behavior of polynucleotide chains, and in particular to clarify the origins of the barriers to internal rotation in the phosphodiester linkage, we computed, with a quantum-mechanical ab initio procedure, the energies associated to 86 combinations of the two torsion angles in the dimethylphosphate anion (CH₃O)₂PO₂ ^–, and then we sought for an analytical expression apt to reproduce these energies with the highest possible accuracy. An excellent agreement (standard deviation of the fitted energies from the ab initio energies 0.28 kcal/mole) with the quantum-mechanical calculations was reached with a potential consisting of four terms: 1) a 6–12 Lennard-Jones contribution, in which different parameters are used to describe the interactions of methyls with the ester oxygens and with the anionic oxygens; 2) a contribution with twofold periodicity, accounting for the anomeric effects connected to the interactions between the lone pair electrons and the polar bonds of phosphorus with the anionic oxygens; 3) a contribution with threefold periodicity, representing the usual bond-staggering term; and 4) a Coulombic contribution, arising from electrostatic interactions between partially charged atoms.     

Theoretical interpretation of the valence X-ray photoelectron spectrum of TiC

The valence X-ray photoelectron spectrum of TiC is interpreted in terms of the density of states and in terms of a simple model for transition matrix elements. The conceptual difficulties arising with both the initial and the final state are discussed.     

The electronic structure and one-electron properties of BH computed from SCF and CI wave functions

The SCF and CI wave functions for BH, obtained in calculations described in detail elsewhere [2], are compared through their electron distributions and electron moments.Taken in part from a Ph.D. thesis submitted to the University of Toronto in 1971.     

The Virtual Charge Model of a polarizable medium as a basis for Hückel calculations with the ω-technique

It is shown that the consideration of a strongly polarizable medium using the Virtual Charge Model (VCM) allows for the possibility of treating independently each electronic subsystem of a molecule by a Hückel-type -technique. The case of a finite polarizable environment can be treated by a first-order perturbation calculation.