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.     

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.     

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.     

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.     

Electronic structure of unstable intermediates

The geometry of the linear molecule HBO, has been investigated within the restricted Hartree-Fock LCAO-MO-SCF approximation. The calculated bond lengths for the near Hartree-Fock calculation were R(H-B)=2.1913 bohr, R(B-O)=2.2284 bohr. Several one electron properties have been calculated for the minimum energy configuration.     

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.     

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.     

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.     

Excited states of phenyl carbonyl compounds

In an attempt to clarify the origin of the dual phosphorescence in phenyl alkyl ketones, we have made some calculation (within the C.I.P.S.I. method in an excitonic scheme) to elucidate the conformation of both ground states and excited states of propiophenone. Our calculations have shown the presence of two stable isomers in the ground state, first n ^* state, and first ^* singlet and triplet states. So our work suggests that the origin of the dual phosphorescence of propiophenone could be related to the conformational change of the molecule in the n ^* state, because the most stable conformations in the n ^* state and in the ground state are different.     

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.     

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.     

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.     

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.     

Basis set dependence of localized orbitals

The effects of Gaussian basis set contraction and addition of polarization functions on H₂O localized orbitals have been studied at the experimental geometry. It is shown that the electric moments and moment features of localized orbitals are not influenced very much by basis set quality variations, as going from medium size to enlarged basis sets. The difference between bond pair and lone pair charge densities was found to be larger on approaching the Hartree-Fock limit. A minimal basis set, however, does not suitably characterize the localized charge distributions.     

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.     

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.     

Atom-surface elastic scattering in the presence of laser radiation

A qualitative discussion on atom-surface elastic scattering taking place in the presence of laser radiation is given. It is suggested that appreciable effects of laser radiation on diffraction patterns may be expected if the laser radiation is capable of inducing electronic transitions in atoms with a large probability.Research sponsored by the Air Force Office of Scientific Research (AFSC), United States Air Force, under Contract No. F49620-78-C-0005, and the National Science Foundation under Grant No. CHE77-27826. The United States Government is authorized to reproduce and distribute reprints for governmental purposes notwithstanding any copyright notation hereon.Alfred P. Sloan Research Fellow, 1976–1980; Camille and Henry Dreyfus Teacher-Scholar, 1975–1980.