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.     

Theoretical study of the potential energy curves of the diatomic radicals MeIIX. III. Application to MgCl,CaF and CaCl radicals and some preliminary remarks

The series of calculations of the potential energy curves of the diatomic radicals Me_IIX (Me_II = second group metal, X = halogen), has been extended to MgCl, CaF and CaCl. The calculations have been performed according to a stepwise procedure, outlined in previous works. The presently available results allow comparisons for the first members of the series.     

Geometry optimization inab initio SCF calculations

The floating orbital geometry optimization (FOGO) described previously [1, 2] for atoms without polarized inner-shell electrons, is extended to the general case. Instead of the Hellmann-Feynman force a special gradient is calculated analytically and utilized in a variable metric procedure simultaneously with the ordinary energy gradient. Test calculations on a sample of 12 molecules were performed to check the efficiency of the method. The geometries obtained are better than those obtained with the corresponding double-zeta basis set. The most striking results, however, are excellent dipole moments.     

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.     

Molecular mechanics calculations on the differentiation of diastereomeric complexes ofcis-decalin withβ-cyclodextrin

Molecular mechanics calculations with the latest available version of Allinger's MM2 force field (MM2(91)) on the diastereomeric complexes of both enantiomeric conformations ofcis-decalin with -cyclodextrin show a small preference (1.67 kJ mol^–1) for one of them, in agreement with the available¹³C-NMR results. Calculations were found to be sensitive to the procedure used.     

Ab initio Calculations of the potential energy surface of the reaction of singlet methylene with the hydrogen molecule

The potential energy surface for the insertion of singlet methylene into H₂ has been computed on theab initio SCF level as well as with inclusion of electron correlation by means of the CEPA method. The results are compared with those of previous semiempirical,ab initio SCF and CI calculations. The system is a prototype of a reaction where an allowed and a symmetry-forbidden path can compete. The electron correlation energy was found to be very different for different regions of the surface, but did not have much influence on the optimum reaction path. From the computed heat of the reaction, the heat of formation of singlet methylene was estimated to be 101.5 kcal/mol. According to the calculations the reaction does not need any activation energy.     

Inclusion of dynamic σ-π polarization in π-electronab initio calculations

Summary A method is presented, whereby dynamic - polarization, i.e. the correlation effect expressed by simultaneous (-*, -*) excitations, can be approximately included in a multi-reference configuration interaction (MRCI) or multi-configurational self-consistent field (MC-SCF) calculation, without need to explicitly correlate the sigma orbitals. The method, which we call the capacitance matrix method, is based on the use of conventional one-electron integrals, from which a polarization potential (SPP) contribution is computed and added to the one- and two-electron Hamiltonian. In the present form, the method requires one parameter for each type of atom, and one for each type of bond. These parameters were adjusted to reproduce the dynamic - polarization energy, computed by restricted multi-reference CI calculations, of a number of states of different hydrocarbons, and the agreement was within a few percent. Using the same parameters in CAS (Complete Active Space) SCF calculations of various states of benzene gives excitation energies, when SPP is included, which is comparable to those obtained by much more elaborate MRCI calculations.     

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.     

Spectroscopic and quantum chemical investigation of enol-enol tautomerism of 2-acetyltetronic acid

Vibratonal spectra of 2-acetyltetronic acid (ATA) (3-acetyloxolane-2,4-dione) were investigated in different aggregate states.Ab initio quantum chemical calculations (3-21G basis set) of 2-formyltetronic acid as analogue of ATA and semi empirical calculation of ATA were carried out. The tautomeric forms of ATA can be arranged in the following series with respect to thermodynamic stability (in decresing order):1a,1b,2a,2b.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 6, pp. 1043–1048, June, 1995.     

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

An optimum minimal basis of exponential atomic orbitals

An analysis of some of the energetic properties of the conventional minimal STO basis is used to suggest a new optimum set of exponential functions for use in molecular calculations.     

An examination of UHF methods for atoms in the first three rows

The reliability and usefulness of UHF methods of calculating hyperfine parameters has been examined. Comparisons with CI and SEHF methods are made.A detailed analysis of OPHF calculations for first and second row atoms has shown that various one-electron properties are accurately described by physically reasonable functions of the atomic number Z. In addition there is a strong correlation between these properties. This has led to a method whereby UHF spin densities across a row of atoms can be obtained from UHF calculations of only two atoms in that row.A strong correlation between experimental and UHF spin densities is shown to exist for atoms of the first three rows. This is used to predict experimental spin densities for atoms in these rows which have not yet been measured experimentally.     

Gas chromatography of halogenated adamantanes

The main thermodynamic characteristics of the adsorption of haloadamantanes on graphitized thermal carbon black were determined by experimental measurements and by calculations in terms of molecular-statistic theory of adsorption. Using experimental data, the Kovac retention indices were calculated and the optimum conditions for gas-chromatographic separation of haloadamantanes on stationary phases with different polarity were elucidated. The influence of the cage effect in the adamantane unit on the chromatographic properties of haloadamantanes was established.     

Application of the independent electron pair approach to the calculation of excitation energies,ionization potentials,and electron affinities of first row atoms

Excitation energies, first ionization potentials and electron affinities of first row atoms are calculated with a spin-adapted independent electron pair approximation (IEPA) combined with the direct determination of pair natural orbitals (PNOs). To enable comparison with molecular calculations Gaussian basis sets are used which are small enough to be also applicable to molecules. IEPA results for the above mentioned properties are accurate to 0.1–0.3 eV which is almost one order of magnitude better than the corresponding SCF-results. The same accuracy can be expected for molecules in which a localization of the doubly and singly occupied orbitals is possible, for instance for small hydrides. This is supported by the results of calculations on carbon hydrides.     

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.     

All-electron SCF LCAO MO calculations on various conformations of cis- and trans-stilbene

Ab initio SCF calculations of cis- and trans-stilbene at different conformations were performed using two program systems. Minimal energy is obtained for cis-stilbene when the phenyl rings are rotated by 52 ° out of the molecular plane. The deviation from planarity due to steric hindrance is smaller for the trans isomer yielding a rotational angle of 19 °. The trans isomer is calculated to be more stable by 5.7 kcal/mole than the cis isomer, confirming the experimental estimate according to which the energy of isomerization is about 3 kcal/mole. This is an improvement over semiempirical calculations which predict a lower energy for the trans configuration.     

A simulation of the exchange potential in unrestricted and restricted Hartree-Fock calculations studied on atoms

The spherical average of the Hartree-Fock exchange potential depending on each spin orbital is compared with Slater's exchange potential, V _xs, as demonstrated for the phosphorus atom. It is shown that the former potential can be simulated by (a + br)V _xs, where r is the radius and the constants a and b are calculated for each spin orbital. This simulation is tested for the iron atom and it is found that the results agree well with those obtained from unrestricted and restricted Hartree-Fock calculations, respectively. The applicability of this new method in energy band structure calculations is briefly discussed.Dedicated to Professor H. Hartmann on the occasion of his 65th birthday.     

Asymmetric reduction of acetophenone and propiophenone by NaAl(IPTOLate)H2 combined with enantiomeric enrichment of the reaction product as an inclusion complex with IPTOL

An improved procedure was developed for asymmetric reduction of acetophenone and propiophenone by the chiral reagent NaAl(IPTOLate)H₂. This procedure is based on isolation of the chiral alcohol that formed as a crystalline host—guest complex with the IPTOL ligand. The enantiomeric enrichment of the product was as high as 97% ee. The ability of IPTOL and its analogs to form host—guest complexes with a number of ether-type solvents, 1-phenylethanol, and 1-phenylpropan-1-ol as well as thermal stabilities of IPTOL-containing complexes with these alcohols were studied.