Deorthogonalization of atomic originals in the CNDO approach

The procedure for deorthogonalization (D) of atomic orbitals in the semiempirical CNDO approach is reviewed. For comparative studies, CNDO/2, CNDO/2D, and STO -3G calculations of molecular dipole moments and Mulliken populations are carried out on 35 prototype molecules containing H, C, N, O, and F atoms. The calculated values are assessed on the basis of how well they agree with experimental trends, chemical bonding theories, and ab initio molecular orbital (MO) values. Results of analyses indicate that the CNDO/2D values for dipole moments are in reasonable agreement with experimental values, and those for net atomic charges and electron populations bear greater resemblance to the ab initio (STO -3G and 6-31G**) values than the original CNDO/2 values. These findings, together with those of previous investigators, demonstrate unequivocally the advantages of incorporating deorthogonalization into routine CNDO/2 or INDO calculations as a means to obtain reasonable estimates of charge distributions.     

A CNDO/INDO molecular orbital formalism for the elements H to Br. applications

The CNDO/INDO molecular orbital formalism introduced in the preceding paper has been applied to a large number of atom combinations up to bromine under the inclusion of the first transition metal series. The results are compared with experimental data (geometries, ionization potentials, dipole moments) or with the results of sophisticatedab initio calculations (one electron energies, net charges, atomic populations). The semiempirical model reproduces for a wide range of molecules the experimental andab initio data with remarkable success.     

Calculations of the isotropic hyperfine coupling constants in free radicals

For a number of free radicals the results of non-empirical (ab initio) and semi-empirical (INDO, DEPAC, CNDO/SP) calculations of the isotropic hyperfine coupling constants are compared.     

Valence shell calculations on polyatomic molecules

Dipole moments and charge distributions for twenty molecules of widely different types have been calculated using (a) the CNDO/2 method and (b) a CNDO/2D method in which the orbitals from the CNDO/2 method are deorthogonalized by a Löwdin transformation and are then used to calculate the dipole moments in a rigorous manner. A statistical analysis of the results for the dipole moments calculated by the CNDO/2D method shows that they are in very slightly better agreement with experiment than those from the CNDO/2 method. The net charge distributions from the CNDO/2D method follow more closely the trends of ab initio calculations than do the CNDO/2 net charges.

---

Zusammenfassung Dipolmomente und Ladungsdichten von Molekülen unterschiedlichen Typs wurden mittels des CNDO/2- und CNDO/2D-Verfahrens (d. i. mit delokalisierten Löwdin-Orbitalen als AO's) berechnet. Eine statistische Analyse zeigt, daß die Resultate der zweiten Methode etwas besser als die der ersten den experimentellen Ergebnissen folgen. Das Analoge gilt für die Nettoladungsverteilungen in bezug auf die Trends bei ab initio-Rechnungen.

---

Résumé Les moments dipolaires et les distributions de charge pour vingt molécules de types divers ont été calculés par: a) la méthode CNDO/2; b) une méthode CNDO/2D où les orbitales de CNDO/2 sont déorthogonalisées par une transformation de Löwdin. Une analyse statistique montre que les moments dipolaires calculés par CNDO/2D sont légèrement en meilleur accord avec l'expérience que ceux calculés par CNDO/2. Les distributions de charge de CNDO/2D sont plus ressemblantes à celles de calculs ab-initio que ne le sont les distributions de CNDO/2.

---

---

This work represents part of the Ph.D. Dissertation submitted to the University of Virginia by D. D. S. and was supported by Grants No. 1-F01-GM41986-01 from the National Institutes of Health, Bethesda, Maryland, U.S.A., and No. AF-AFOSR-1184-67 from the Air Force Directorate of Scientific Research.

---

NASA Research Trainee.     

CNDO CI calculations on second-row molecules

The CNDO/S method is extended to second-row molecules and parametrized for phosphorus and sulphur. Both sp and spd basis sets are considered. The method is applied to the aromatic molecules phosphorin and thiophen. The uv transitions, ionization potentials and dipole moments of these molecules are satisfactorily explained.     

Chlorosilanes: Development and application of MM2 force field parameters

The geometries, relative conformational energies, and dipole moments of mono and polychlorosilanes have been calculated using ab initio molecular orbital (MO) theory. Calculations at the HF/3–21G(*) level, with the exception of dipole moments, give reasonable agreement with experimental data. A new MM2 force field for chlorosilanes, which includes terms for bond length shortening and bond angle compression due to the attachment of electronegative Cl atoms, has been developed on the basis of experimental and ab initio results. The new force field is generally successful in predicting structural parameters, but is unable to reproduce the dipole moments of several model systems. While dipole moment predictions are not the authors' main interest, this failure defines a shortcoming in the MM2 method. The new parameters have been applied to problems in the prediction of stereochemistries of cyclic systems, and compared with experimental results where data are available.     

A theoretical study on the structure and properties of phenothiazine derivatives and their radical cations

Semiempirical CNDO, AM1, PM3 and ab initio HF/STO-3G, HF/3-21G(d), and HF/6-31(d) methods were employed in the geometry optimization of the phenothiazine and the corresponding radical cation. The results obtained from the PM3 performances were as good as those from the ab initio calculations in the structure optimization of both phenothiazine and phenothiazine radical cation. The PM3 method was used to optimize the structures of a series of N-substituted phenothiazine derivatives and their radical cations. The PM3-optimized results were then analyzed with the ab initio calculation at the 6-311G(d,p) level, which yielded the total energy, frontier molecular orbitals, dipole moments, and charge and spin density distributions of the phenothiazine derivatives and their radical cations.     

Eine Variante des CNDO-Verfahrens unter Einbeziehung vond-Funktionen

A modified handling of thed-functions within the CNDO framework is presented. No change of the (CNDO/2) atomic parameters is necessary from sodium to iodine. One-centre kinetic energy contributions to the core-matrix elements are introduced explicitly. The calculation of dipole moments and equilibrium geometries (estimated with the aid of analytically calculated gradients of the molecular energy) of 45 molecules gives good results,sp- andspd-CNDO/2 have been applied for comparison. Some examples are presented, where the shortcomings of the conventional CNDO/2 are no longer present. A qualitatively accurate picture of the dipole moments of X-Y-molecules (X: Halogen; Y: Hydrogen, Halogen, Alkyl-) is achieved.

     

Application of the theory of symmetry of nonrigid molecules to the calculation of the conformational dependences of the torsional potential and dipole moment of acetone-related molecules

We have investigated the conformational dependences of the torsional potential and dipole moment of double-rotor molecules related to acetone, using semiempirical and ab initio calculations and expressing the results in terms of limited Fourier-series expansions. The use of the isodynamic operations of nonrigid molecules to obtain symmetry-adapted quasianalytic forms for the various properties helps compute the representative surfaces with a minimum number of points. Potential surfaces have been calculated for planar ground-state acetone (CNDO /2, STO /3G, and STO /4-31G) and both pyramidal excited-triplet acetone and ground-state dimethylamine (CNDO /2). For groundstate acetone STO /4-31G brings the results obtained with STO /3G closer to those from CNDO /2 and from experiment. The potential surface of excited-triplet acetone appears intermediate between those of ground-state acetone and dimethylamine. For dipole moments the convergence of the harmonic expansions of the vector components is slower than that of the torsional potential whereas that of the vector magnitude is faster.     

Conformational energies and rotational barriers in 3-methyl-1-butene and 1-butene: An ab initio and molecular mechanics study

The calculated result obtained with MM2(87) for the rotation of the isopropyl group in 3-methyl-1-butene is not in agreement with experimental data. In order to reparametrize the C_sp2-C_sp3-C_sp-C_sp3 torsional angle, 3-methyl-1-butene and 1-butene have been studied by molecular mechanics (MM2(87)) and ab initio (MP2/6-31G* and MP3/6-31G*) calculations. The reparametrization of the torsional angle gives calculated results from MM2(87) in agreement with experimental data and ab initio calculations for both 3-methyl-1-butene and 1-butene. The calculated barriers for the rotation of alkyl groups in alkylbenzenes are improved with these new parameters.     

Predicted vibrational frequencies and assignment of the spectra of uracil and its n,n-dideutero derivative by quantum chemical methods

A recent theoretical study on the vibrational spectrum and force field of maleimide (ref.1) lead to scale factors transferable to molecules of similar structure (ref.2). The theoretical vibrational spectrum of uracil (refs.2–3) calculated from scaled CNDO/2 force constants using the scale factors of ref.1 agrees better with the experimental results (ref.4) than former calculations by $\text{ab initio}$ STO-3G (ref.5) or MINDO/3 (ref.6) methods. The scaled CNDO/2 force field of uracil has also been used to predict the normal frequencies of a set of its deuterated derivatives.     

A CNDO/INDO crystal orbital model for transition metal polymers of the 3d series—basis equations

The crystal orbital formalism in the tight-binding approximation is combined with a recently developed CNDO/INDO model for transition metal species of the 3d series in order to allow band structure calculations on the Hartree-Fock (HF) SCF level for one-dimensional (1D) chains with organometallic unit cells. The band structure approach based on the CNDO and INDO approximation can be used for any atom combination up to bromine under the inclusion of the 3d series. The matrix elements for the tight-binding Hamiltonian are derived for an improved CNDO and INDO framework. The total energy of the 1D chain is partitioned into one-center contributions and into two-center increments of the intracell and intercell type. Semiempirical band structure calculations on simple model systems are compared with available ab initio data of high quality.     

Zur quantenchemischen Behandlung von Kation-Amid-Solvatkomplexen

A comparison ofab initio calculations employing different basis sets with corresponding CNDO/2 results for the Li⁺/HCONH₂ complex shows that these methods lead to completely different energy surfaces for this system. Reduction of the basis set, even to the minimal size, does not bring about serious changes in the results of theab initio calculations, whereas in the semiempirical treatment some methodical errors seem to occur. When using however, theab initio minimum geometry the CNDO/2 calculations also give a qualitatively correct picture for the influence of the cation on the amide modecule.     

Determination of one-centre core integrals from the average energies of atomic configurations

One-center core integrals for valence orbitals are determined from the experimental average energies of neutral atomic configurations from Li through Zn. These values are compared with those estimated from CNDO /1, “INDO /1”, CNDO /2, “INDO /2” and with theoretical values calculated from a pseudo-potential method. The agreement is good between values obtained from neutral atoms and from the psuedo-potential calculation except for the 3d orbitals of the transition elements where the theoretically calculated integrals over single ξ functions are not realistic. These two methods reproduce both term and average configuration energies for the first two rows of atoms; the semiempirical method reliably reproduces them for the third row. The CNDO /1 and INDO /1 methods underestimate atomic energies, while the CNDO /2 and INDO /2 procedures fail rather poorly. The propriety of using core integrals estimated semiempirically in molecular orbital calculations is discussed.     

Geometry optimization and electronic structures of molecules H_3AXAH_3

The geometries of molecules H_3AXAH_3(X=O,S,Se and A=C,Si)have been optimizedusing STO-3G ab initio calculations and gradient method and the results are in good agreement withreported experimental values.From the STO-3G optimized geometries,we have also calculated theelectronic structures of these molecules using 4-31G and 6-31G basis sets to obtain the MO energies.atomic net charges and dipole moments.The ionization potentials calculated by 6-31G basis set are ingood agreement with experimental values.     

Conformational analysis—CXVII : Methyl ethyl disulfide. Molecular mechanics and molecular orbital calculations

The conformation of methyl ethyl disulfide was investigated by molecular mechanics calculations using a recently developed force field for sulfur-containing alkanes. The results indicate that in the gas phase the molecule exists predominantly in two conformations, both with the CSSC dihedral angle gauche (84°), and the SSCC dihedral angle either gauche (72°) or trans (179°), and the methyl protons staggered. Ab initio molecular orbital calculations using an STO-3G basis set were employed to corroborate that these two conformations are of roughly equal stability, and that the next most stable conformation (by 0.6 kcal/mole) has the SSCC dihedral angle gauche (295°) with the terminal methyls proximal. In contrast to earlier CNDO/2 (spd) predictions, the SSCC cis conformer is the least stable, and no sizable attractive S?HC nonbonded interactions are discerned. Reasons for this are traced to a failure of the CNDO/2 method, which is especially serious when d orbitals are included in the basis set (spd) and the rigid rotor approximation is used. The present results are found to be consistent with recent electron diffraction, IR, Raman spectroscopic and X-ray diffraction data. The conformation of diethyl disulfide was also investigated by molecular mechanics calculations, and again gauche and trans SSCC arrangements are predicted to be preferred.     

On the molecular electrostatic potentials obtained with CNDO and INDO wave functions

Molecular electrostatic potentials computed with CNDO/2 and INDO wave functions are shown to present systematic differences with respect to ab initio potentials in the case of out-of-plane potentials and in-plane vicinal hetero atoms in planar hetero molecules.     

The electronic structure of the HCN dimer and trimer

The electronic structure and energy of dimerization and trimerization of HCN are computed with an STO-3G basis and the results found to be in good agreement with the experimental E. Unlike CNDO/2, this small ab initio basis predicts the correct geometry for the dimer of hydrogen cyanide. The charge redistribution effects found in this H-bond involving a C-H proton donor and sp hybridized acceptor are similar to those found in previous H-bonded studies.     

A charge-conserving approximation method for ab initio calculations

A new method is presented for approximate ab initio calculations in quantum chemistry. It is called CCAM (charge conserving approximation method). The calculation method does not include the use of empirical parameters. We use Slater type orbitals as basis set, replacing STO's by STO-2G functions to evaluate three- and four-center integrals and making the STO-2G two-orbital charge distributions have the same total charge as STO. The results are presented for test calculations on five molecules. In view of these results, CCAM is better than ab initio calculations over STO-6G in the results on total energies, kinetic energies and occupied orbital energies. In atomic populations, dipole moments and unoccupied orbital energies, CCAM is also satisfactory. We estimate that CCAM would be as fast as ab initio calculations over STO-2G in evaluating molecular integrals.