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.     

Recursion formulae for calculation of overlap integrals

Multi-ζ Slater-type orbitals are frequently used in molecular orbital calculations. Master formulae and numerical tables are available in literature for overlap integrals between s, p, and d atomic orbitals up to principal quantum number (n) = 3 and for some other selected quantum numbers. However, no master formula or numerical table is available for quantum numbers n = 5 and above and involving ? orbitals. In this article recursion formulae have been presented for the calculation of the overlap integral between any two s, p, d, and ? atomic orbitals formed by a linear combination of Slater-type orbitals. These formulae, when expanded, would give rise to all the master formulae reported in the literature as well as formulae hitherto unreported.     

Multicenter point charge model for high-quality molecular electrostatic potentials from AM1 calculations

The natural atomic orbital/point charge (NAO-PC) model based upon the AM1 wave function has been developed to calculate molecular electrostatic potentials (MEPs). Up to nine point charges (including the core charge) are used to represent heavy atoms. The positions and magnitudes of the eight charges that represent the atomic electron cloud are calculated from the natural atomic orbitals (NAOs) and their occupations. Each hybrid NAO is represented by two point charges situated at the centroid of each lobe. The positions of the centroids and the magnitudes of the charges were obtained by numerical integration of the Slater-type hybrids and the results used to set up polynomials and look-up tables that replace the integration step in the actual MEP calculation. The MEPs calculated using this method are found to be in better agreement with those obtained using RHF/6-31G* than those obtained from the AM1 wave function using Coulson charges or with MOPAC-ESP. The MEP calculations are extremely fast and have, for instance, been incorporated into an interactive graphics package. © 1993 John Wiley & Sons, Inc.     

STOP: A slater-type orbital package for molecular electronic structure determination

A general ab initio package using Slater-type atomic orbitals is presented. This package, called STOP, uses the one-center two-range expansion method to evaluate the multicenter electronic integrals. Thoroughly optimized numerical techniques, in particular, convergence accelerators and suitable Gauss quadratures, are used in the algorithms which provide accurate numerical values for all these integrals. STOP thus provides wavefunctions for general molecular structures at the self-consistent field level for the first time over a Slater-type orbital basis. © 1996 John Wiley & Sons, Inc.     

On the use of spatial symmetry in ab initio calculations. Transformation of the two-electron integrals from atomic orbital to localized molecular orbital basis

Anm ⁵-dependent integral transformation procedure from atomic orbital basis to localized molecular orbitals is described for spatially extended systems with some Abelian symmetry groups. It is shown that exploiting spatial symmetry, the number of non-redundant integrals for normal saturated hydrocarbons can be reduced by a factor of 2.5-3.5, depending on the size of the system and on the basis. Starting from a list of integrals over basis functions in canonical order, the number of multiplications of the four-index transformation is reduced by a factor of 2.8-3.5 as compared to that of Diercksen's algorithm. It is pointed out that even larger reduction can be achieved if negligible integrals over localized molecular orbitals are omitted from the transformation in advance.     

EHMO方法的分子轨道定域化

用Foster-Boys的定域化准则讨论了EHMO方法的分子轨道定域化问题,提出用双中心重叠积分近似计算双中心轨道偶极矩积分方法,得到的EHMO定域分子轨道与严格定域化结果接近,与从头计算方法的定域化结果定性一致。     

A unified scheme for the calculation of differentiated and undifferentiated molecular integrals over solid-harmonic Gaussians

Utilizing the fact that solid-harmonic combinations of Cartesian and Hermite Gaussian atomic orbitals are identical, a new scheme for the evaluation of molecular integrals over solid-harmonic atomic orbitals is presented, where the integration is carried out over Hermite rather than Cartesian atomic orbitals. Since Hermite Gaussians are defined as derivatives of spherical Gaussians, the corresponding molecular integrals become the derivatives of integrals over spherical Gaussians, whose transformation to the solid-harmonic basis is performed in the same manner as for integrals over Cartesian Gaussians, using the same expansion coefficients. The presented solid-harmonic Hermite scheme simplifies the evaluation of derivative molecular integrals, since differentiation by nuclear coordinates merely increments the Hermite quantum numbers, thereby providing a unified scheme for undifferentiated and differentiated four-center molecular integrals. For two- and three-center two-electron integrals, the solid-harmonic Hermite scheme is particularly efficient, significantly reducing the cost relative to the Cartesian scheme.     

Finite jellium models. I. Restricted Hartree-Fock calculations

Restricted Hartree-Fock calculations have been performed on the Fermi configurations of n electrons confined within a cube. The self-consistent-field orbitals have been expanded in a basis of N particle-in-a-box wave functions. The difficult one- and two-electron integrals have been reduced to a small set of canonical integrals that are calculated accurately using quadrature. The total energy and exchange energy per particle converge smoothly toward their limiting values as n increases; the highest occupied molecular orbital-lowest unoccupied molecular orbital gap and Dirac coefficient converge erratically. However, the convergence in all cases is slow.     

Analytical evaluation of three-center nuclear-attraction integrals over s-Slater orbitals for asymmetrical conformations of the centers

Analytical formulas for three-center nuclear-attraction integrals over Slater orbitals are given for any location of the three atomic centers. In the mathematical derivations the Neumann expansion has been used and new general auxiliary integrals which depend on the elliptical coordinates of one of the centers are defined. The orbital exponents within the integrals may be different.     

Theoretical investigation of the magnetic interactions of Ni9 complexes

On the basis of density-functional theory (DFT) calculations, a theoretical analysis of the exchange interactions in Ni9L2(O2CMe)8{(2-py)2CO2}4, was performed, where L is a bridging ligand, OH- (1) or N3- (2). Each magnetic interaction between the Ni spin centers is analyzed for 1 and 2 in terms of exchange integrals (J values), orbital overlap integrals (T values) and natural orbitals. It was found that a J3 interaction, which is a magnetic interaction via the bridging ligand orbitals, mainly controls the whole magnetic properties, and the dominant interaction is a sigma-type orbital interaction between Ni dz2 orbitals. Further investigations on the magnetostructural correlations are performed on the J3 interactions using simplest Ni-L-Ni models. These models reproduced the magnetic interactions qualitatively well not only for the Ni9 complexes but also for other inorganic complexes. Strong correlations have been found between the magnetic orbital overlaps (T values) and the Ni-L-Ni angle. These results revealed that the difference of the magnetic properties between OH- and N3- is caused by the orbital overlap integral (T values) of the sigma-type J3 interaction pathway. The magnetic interactions are also discussed from a Hubbard model by evaluating the transfer integral (t) and on-site Coulomb integrals (U), in relation to the Heisenberg picture.     

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.     

Ab initioLCAO-MO-SCF Calculation of chlorpromazine and promazine

Ab initio LCAO -MO -SCF Gaussian basis function calculations have been performed for chlorpromazine and promazine. By prescreening for the size of the integrals before calculation, it was only necessary to calculate 12 million out of the possible 38.5 million integrals for chlorpromazine. By a novel procedure of processing the integral tapes for the SCF it was possible to cut down significantly on the amount of time for the SCF . The SCF calculations converged smoothly for both promazine and chlorpromazine. There is a sizeable energy gap between the energy of the highest occupied molecular orbitals in these molecules (which is of the order of ?0.3 a.u.) and the lowest unoccupied molecular orbital (which is of the order of + 0.15 a.u.). The gross atomic populations of chlorpromazine and promazine resemble each other and differ only somewhat on the carbon atom to which the substituent is attached and the carbons and their hydrogens adjacent to it.     

New method for the direct calculation of electron density in many-electron systems. I. Application to closed-shell atoms

Analytical properties of hydrogen-like atomic orbitals (HAO ) that are used in the MOLCAO approach to the quantum theory of molecules have been studied. Addition and expansion theorems for HAO have been proved, both in coordinate and momentum representations. A close relation has been established between HAO and the reduced Bessel functions of half-integer indices. New methods are suggested to calculate integrals for atomic and molecular form factors, and multicenter integrals, for the HAO basis in the MO LCAO theory.     

Evaluation of Orbital- and Ground State Energies of Some Open- and Closed-Shell Atoms over Integer and Noninteger Slater Type Orbitals

Ab initio calculations of the orbital and the ground state energies of some open- and closed-shell atoms over Slater type orbitals with quantum numbers integer and Slater type orbitals with quantum numbers noninteger have been performed. In order to increase the efficiency of these calculations the atomic two-electron integrals were expressed in terms of incomplete beta function. Results were observed to be in good agreement with the literature.     

Synthesis,Characterization and Electronic Properties of a Symmetric Bidentate Schiff-Base Ligand and Its Complex with CadmIum(II)

A symmetric bidentate Schiff-base ligand and its complex with Cd(II) is described. The ligand and its complex were characterized by microanalysis, UV-Vis, GC-Mass and FT-IR spectroscopic methods. The present work also studied the structures and electronic properties of the bidentate Schiff-base ligand by using ab initio and AM1 molecular orbital methods. Ab initio and AM1 geometrical predictions have been compared. The analysis of molecular orbitals (MOs) indicates that the N(3) and N(6) atoms could be the coordination sites in this ligand.     

On obtaining localized bonding schemes from delocalized molecular-orbital wave functions

A straightforward procedure is proposed for expanding a molecular orbital determinantal wave function into a set of determinantal wave functions composed of atomic orbitals localized at the atoms of a molecule. By employing this method, atomic orbital determinants and their weights can be derived for a molecule from the computed molecular-orbital wave function. The procedure permits the interpretation of a molecular orbital determinantal wave function in terms of bonding schemes related to the classic resonance structures used by organic chemists. By using the unrestricted molecular orbital determinant, bonding schemes and their weights are obtained for butadiene, the butadiene radical cation and the acrylonitrile radical anion. Their dominant bonding schemes are in accord with the relevant resonance structures for these molecules. For the butadiene radical cation and the acrylonitrile anion they are shown to be compatible with the accepted mechanisms of the electrochemical coupling reactions of butadiene and acrylonitrile. Received: 7 August 1996 / Accepted: 18 March 1997     

Localized orbitals based on the fermi hole

The Fermi hole provides a direct (non-iterative) method for tansforming canonical SCF molecular orbitals into localized orbitals. Except for simple overlap integrals required to maintain orthogonality, this method requires no integrals over orbitals or basis functions. This method is demonstrated by application to a furanone (C₄H₄O₂), methylacetylene, and boron trifluoride. The results of these calculations are compared to those determined by the orbital centroid criterion of localization.     

Study of localized molecular orbitals using group theory methods and its approach to the many-electron correlation problem. IV. The symmetry-adaptation of many-center integrals and hamiltonian matrix elements in MCSCF calculations

Group theoretic methods are presented for the transformations of integrals and the evaluation of matrix elements encountered in multiconfigurational self-consistent field (MCSCF) and configuration interaction (CI) calculations. The method has the advantages of needing only to deal with a symmetry unique set of atomic orbitals (AO) integrals and transformation from unique atomic integrals to unique molecular integrals rather than with all of them. Hamiltonian matrix element is expressed by a linear combination of product terms of many-center unique integrals and geometric factors. The group symmetry localized orbitals as atomic and molecular orbitals are a key feature of this algorithm. The method provides an alternative to traditional method that requires a table of coupling coefficients for products of the irreducible representations of the molecular point group. Geometric factors effectively eliminate these coupling coefficients. The saving of time and space in integral computations and transformations is analyzed. © 1994 by John Wiley & Sons, Inc.     

Spin-orbit interaction in polyatomic molecules: Ab initio computations with Gaussian orbitals

Standard sets of Gaussian atomic orbitals (STO -3G , STO -4.31G ) are used to evaluate spin-orbit coupling constants in linear molecules (CO, NNN) and spin-orbit effects on singlet–triplet transition intensities in formaldehyde. All spin-other orbit effects have been included. In all cases spin-other orbit interactions form a large fraction of the matrix elements. Simple formulae to evaluate spin-orbit one- and two-electron integrals over atomic orbitals are presented. Standard molecular integral programs can be used for the computation of spin-orbit integrals.     

等性扩展基杂化轨道的构造

杂化轨道理论近来有了长足的发展,杂化轨道的构造方法主要有群论方法最大重叠方法,自然杂化轨道法和其它以分子轨道为基础构造杂化轨道的方法等。其中最大重叠杂化轨道不仅满足正交化条件而且能较定量地考虑到配体轨道的作用,因而已经得到广泛应用。本文在最大重叠原理的基础上得到了扩展基杂化轨道的解析形式。扩展基杂化轨道对一给定几何构型的分子M(X_1X_2…X_n),其中心原子M的n个杂化轨道与诸配体{X_i}形成一组方向键。M的杂化轨道(HO)和原子轨道(AO)可分别作为该分子对称操作群的表示之基。这两种不同基的表示进行约化之后,属于同一不可约表示的HO和AO是线性相