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.     

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.     

Coupling constants in the direct configuration interaction method

A general algorithm of evaluation of the coefficients of molecular integrals (coupling constants) appearing in the direct configuration interaction method is derived. The configurations are assumed to be spin-adapted antisymmetrized products of orthonormal orbitals. No limitation is imposed either upon the reference state (the number of the singly occupied orbitals may be arbitrary) or upon the excitation multiplicity.     

Calculation of polarizabilities by the CNDO method

The polarizabilities of a number of small and medium size molecules are calculated using CNDO/2 wave functions. The first-order and the second-order perturbation equations are solved by variational procedures. The results are encouraging for the component along the internuclear axes and in the molecular plane. The effect of deorthogonalization of the CNDO/2 orbitals is discussed.     

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

Conformations of the thiathiophthenes and related molecules

A conformational ab initio MO study has been carried out for the thiathiophthene molecule (TTP) and two related model compounds, thiomalonaldehyde (TMA) and its conjugate base (TMA(-)). The conformational energy surfaces for TMA, TMA(-) and TTP were generated using a least squares fit to the calculated data and plotted on a CALCOMP plotter. The results of the calculations showed that the cis-cis planar conformation of TTP is the most stable in agreement with experimental findings. For TMA and TMA(-) the cis-cis planar conformation is not the most stable. Contour plots of the five occupied -MO's of TTP show great similarity to those of naphthalene.Less detailed calculations were carried out for 3-hydroxy-prop-2-en-1-thione (HPT) and 3-mercapto-prop-2-en-1-thione (MPT). HPT was shown to be most stable in the cis planar hydrogen bonded conformation in agreement with the experimentally obtained results. For MPT the non-hydrogen bonded planar structure was found to be the most stable.     

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.     

Electronic states of the [2 n ]cyclophanes

Low energy singlet and triplet states for a series of [2 _n ]cyclophanes are discussed in terms of the results of a simple model calculation. Experimental trends can be explained under the assumption of significant - interaction involving the saturated bridges. This interaction destabilizes low energy excimer states, in contrast to the usual red shift observed for alkylbenzenes. The observed near-constancy of the onset of the absorption spectra can be explained by near-cancelation of through-bond and through-space contributions.From 01.01.1984, Royal Danish School of Educational Studies, Department of Chemistry, DK-2400 Copenhagen NV, Denmark     

A configuration interaction study of phosphine using bonded functions

A computational study is made of the effect of basis set upon the energy, properties and inversion barrier of the phosphine molecule. The calculations are performed at both the SCF and CI level. The flexibility of the double zeta basis is discussed in the light of the results.     

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.     

Theoretical study of the electronic structure of diazomethane

The least-energy dissociation path of the ground state of CH₂N₂ was determined fromab initio calculations using in a complementary way basis sets of minimal size (STO-3G) and double-zeta (DZ) quality. The results indicate that the least-energy point of attack of the N₂ molecule on CH₂ (¹ A ₁) is roughly perpendicular to the molecular plane (93 °), the C and N atoms being almost co-linear (angle C-N-N203 ° with outermost N atom pointing away from CH₂). The potential barrier of 1.2 eV found previously on theC _2v dissociation path, disappears completely along the least-energy dissociation path (point groupC _s (out-of-plane)). These findings corroborate the Woodward-Hoffman rules for this process since the outermost orbitals of the two intersecting states found in point groupC _2v (...2b ₁ and ...8a ₁) both correlate to the same irreducible representation (10á) in point groupC _s (out-of-plane).Larger basis set calculations (DZ + polarization functions on all centers, 3d _c and 3d _N developed here), were also carried out on CH₂N₂ (¹ A ₁,³ A ₂ and¹ A ₂) at the¹ A ₁ equilibrium geometry and on CH₂ (³ B ₁) and N₂ (¹ _g ⁺ ) at their respective equilibrium geometries. These calculations, together with consideration of correlation energy differences, yieldD ₀ ⁰ (CH₂N₂,¹ A ₁) = 19 kcal/mole and vertical excitation energies of 67 and 73 kcal/mole for the³ A ₂ and¹ A ₂ states respectively. The latter value is in good agreement with the measured experimental value: 72.4 kcal/mole corresponding to the maximum of intensity in the¹ A ₂¹ A ₁ absorption band.     

Theoretical assignments of the electronic spectrum of acetylene

Non empirical calculations of energies and properties of some excited states of acetylene are presented. A frozen core approximation is used and excitations to , and MO's are taken into account. Both valence and Rydberg states are considered. Assignments of the UV and electron impact spectra are proposed and some questions are raised.     

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.     

Theoretical approach of the conformations of benzylidene aniline using the PCILO method

The conformational properties of benzylidene aniline (BA) were investigated using the PCILO method. Evidence has been provided that the optimization of geometry plays a prominent part, in agreement with recent experimental work on photoelectron spectroscopy by Haselbach.     

Testing the arbitrariness and limits of a pseudopotential technique through calculations on the series of diatoms HF,AlH, HCl,AlF, AlCl,F2, Cl2

The pseudopotential techniques present some degrees of freedom, the influence of which on molecular calculations must be tested to assess the stability and accuracy of the results. The present work uses a semi-local pseudopotential extracted from near Hartree-Fock atomic calculations; the shape of the inner part of the pseudoorbital, the analytic form of the pseudopotential are shown to have less influence than the choice of the valence basis set which must be optimized. The calculated molecular constants perfectly agree with the large basis set all-electron calculations, even for polar molecules.Equipe de Recherche Associée au CNRS No. 821.     

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.     

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.