Expansion of linear combinations of slater-type orbitals in eigenfunctions of the three-dimensional isotropic harmonic oscillator

Several classes of functions related to the Gaussian have been used with success as basis sets for the representation of atomic and molecular orbitals. We have compared the representation of a hydrogen 1s orbital by a sum of Gaussian lobe functions with its expansion in eigenfunctions of the three-dimensional isotropic harmonic oscillator. The lobe functions are shown to achieve better expectation values of the energy, with fewer terms. The lobe functions have the further computational advantage of not containing high powers of the radius. It is concluded that the lobe functions are a superior basis set for use in calculations of the electronic structure of atoms and molecules.     

Angular dependence of Gaussian-Lobe orbitals. II. Set of axial Gaussian-Lobe orbitals

The topological properties of real spherical harmonic representations on the unit sphere have been found to provide a convenient tool to infer the lobe edifices which mimic these orbitals. The prohibitive number of lobes required in such an approach for l > 2, can be avoided in using only axial Gaussian-lobe orbitals (AGLO ). It is proved that 2l + 1 independent Y_lo-like functions correctly span the relevant Y_lm (m = ?l,l) subspace. The multipolar component analysis of any spatial arrangement of lobes is derived, and allows the optimization of the angular dependence of AGLO S. The cases of d- and f-orbitals are studied in detail and accurate optimized functions are proposed. This method can be easily extended to obtain the atomic orbitals of any azimuthal quantum number l-subspace.     

Angular dependence of Gaussian-Lobe orbitals. I. Analysis of standard p- and d-orbitals

Multipole expansions of Gaussian-lobe atomic orbitals around their centers are theoretically investigated in order to study the exact angular dependence of such functions. Analytical expressions of the multipole coefficients are derived for standard lobe orbitals. It is shown that the average-square values of multipole components are related to a unique orbital parameter λ. The numerical values of p- and d-components are given for selected λ and the choice of this parameter is discussed on the basis of symmetry and computational arguments. The transferability of optimized atomic exponents from harmonic (or Cartesian) functions to lobe functions is established so that the possibility of applying the Gaussian-lobe orbital approach in chemical studies is greatly extended.     

Quality of correlating functions generated from commonly used basis sets

Tests have been performed on the quality of correlating functions generated from commonly used Gaussian basis sets, such as the 4-31G and MIDI-4 sets. The atoms tested were carbon, nitrogen, and oxygen. Self-consistent field and configuration interaction (CI) calculations were performed for the ground and lower excited states of neutral atoms as well as for positive and negative ions, using the original sets. Next, after adding (1) one d, and (2) two d and one f primitive Gaussian-type functions (GTFs) to the original sets, the CI calculations were repeated. In order to investigate the quality of the correlating orbitals generated from the GTF sets, parallel calculations to those for the GTF sets were carried out with an extended set of Slater-type functions. It was found that the excitation energies change in a stepwise manner as the basis sets changed from the original sets to the original set + 1d and the original set +2d1f. The improvements in excitation energies and ionization energies were almost independent of the original sets and were found to be strongly dependent on the augmented correlation functions. © 1996 by John Wiley & Sons, Inc.     

Notes on hydrogenic Gaussian lobe functions

The Gaussian lobe representation of 2p and 3d_xy functions of a one-electron atom has been reconsidered. Contrary to earlier results, it has been found that the minimum energy associated with the simplest case of these functions is reached when the expansion centers coincide with the nucleus. This also seems to be true for representations with larger expansion lengths. Under similar conditions, Whitten's octahedral representation of the 2p function gives a better energy value than Poshusta's tetrahedral one.     

Gaussian expansions of hydrogenic functions involving few variational parameters

Fully optimized lengthy Gaussian expansions of hydrogenic 1s and 2p_z functions have been practically reproduced by generating the variational parameters by means of suitably chosen functions involving only few parameters. Hydrogenic 3d_xy and 4f_xyz functions have been approximated in a similar way. The usefulness of such generating functions is briefly discussed.     

Picture change error correction in the radial distributions of canonical orbital densities and total electron density of radon atom: the effect of the size of nucleus and the basis set limit

The 2nd order Douglas-Kroll-Hess (DKH2) and the Infinite Order Two Component (IOTC) radial distributions of electron density of canonical Hartree-Fock (HF) orbitals of radon atom are presented. Furthermore, the total electron density is revisited. The picture change error (PCE) correction is investigated by analytical means. The point charge model of nucleus and the Gaussian nucleus model are employed. The basis set is extrapolated by means of including tight s and also p Gaussians within the original triple zeta basis set. It is found that the DKH1 PCE corrected DKH2 total electron and s orbital contact densities are negative for the point charge model of nucleus if tight enough s Gaussians are included in the basis set. It is shown that this failure is caused due to the missing terms of the second order Douglas-Kroll transformation for the DKH2 electron density. PCE is found the most striking in the DKH2/IOTC electron density of s orbitals close to the nucleus. The radial distributions of the 2-component p _1/2 orbital densities are considerably affected by PCE at the nucleus as well. Furthermore, the PCE corrected DKH2/IOTC scalar p orbital densities have a non-zero value of electron density at nucleus and can be considered as an spin-orbit (SO) average of the p _1/2 and p _3/2 orbitals. The d and f orbitals are affected by PCE in the vicinity of the nucleus only little. The PCE corrected DKH2 and IOTC radial distributions of orbital densities are nodeless, which is completely in agreement with the radial distribution of the analytic or numeric DCH orbital densities.     

Localized orbital calculations of the bonding in SO,SO2F2, ClO3F and SOCL2

Localized valence molecular orbitals have been obtained for SO, SO₂F₂, ClO₃F and SOCl₂ by the method due to Boys and Foster. The bonding in these molecules, in which the second row atom is exhibiting an excess valency, is discussed in terms of the form of these localized orbitals. The bonding of the second row atom to an oxygen atom is described by three bent bond orbitals, whilst bonding to a halogen atom is described by a single bond orbital. The participation of 3d functions in the various bonding and nonbonding orbitals is analysed in this localized orbital framework.     

Unrestricted CNDO-MO calculations. II. MnO 4 2− and CrO 4 3−

The electronic structures of the tetrahedral molecule ions MnO ₄ ^2– and CrO ₄ ^3– have been investigated within an unrestricted CNDO-MO approximation [Theoret. Chim. Acta (Berl.)20, 317 (1971)]. Calculations assuming the unpaired electron occupies the 3a ₁, 2e, and 4t2 molecular orbitals indicate that the 3a ₁ and2e orbitals have similar orbital energies and that the 4t ₂ orbital is at a higher energy. The experimentally indicated2e orbital for the unpaired electron is obtained with expanded O^1– type atomic orbitals for oxygen and valence metal orbitals of the expanded 3d and plus one ion 4p types. The metal 4s orbitals must be held to the neutral atom type. The optimum valence orbitals above with a slightly contracted 4s type metal orbitals yield the minimum total energy and places the unpaired electron in the 3a ₁ orbital. Since the contracted 4s metal orbital produces results that are not in agreement with experimental data, the method used apparently does not adequately take into account the increased electron-electron repulsions that contracted 4s orbitals produce.     

Multipole expansion of Cartesian Gaussian orbitals about a new origin

Cartesian Gaussian orbitals (CGO) about one center are expanded as a series of spherical harmonics around a new origin. Radial coefficients are given in analytical form. Developments of solid harmonic Gaussian orbitals (SHGO) are also obtained and used to study the convergence properties of expansions of s, p, d, and f Gaussians in terms of a single dimensionless parameter λ. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 64 : 647–653, 1997     

Ab initio calculations,using a small Gaussian basis set,of the electronic structure of the sulphate ion

Ab initio molecular orbital calculations of the electronic structure of the sulphate ion have been performed in which three Gaussian-type functions are used to simulate each member of a minimal basis of Slater-type orbitals. Comparative calculations on H₂S show that such a basis excellently reproduces the properties of the valence electrons given by calculations in a Slater basis. The expansion of the basis by the addition of sulphur 3d orbitals results in a large decrease in the molecular energy (1 a.u.) and has a pronounced effect on the ordering and energy of the molecular orbitals. The results of a number of semiempirical schemes are discussed in the light of these results.     

Contracted polarization functions for B to Ar

Using optimal exponents for B through Ne given by Dunning and those for Al through Ar by Woon and Dunning, d-type contracted polarization functions (2d/1d), (3d/1d), and (3d/2d) are generated from natural orbitals of atomic single and double excitation configuration interaction (SDCI) calculations, where the numbers before and after the slash are those of the primitive and contracted Gaussian type functions. The resulting contracted functions are tested on N₂ and P₂ molecules by self-consistent field and SDCI calculations, which clarify characteristics of the present polarization functions. Received: 5 June 1997 / Accepted: 20 August 1997     

Electronic structure of UH,UF, and their ions

A relativistic effective core potential (REP) has been generated for the uranium atom and used in self-consistent-field calculations of the A states of UH, UF, and their ions. Energy curves were calculated at the base configuration level which ensures the dissociating atoms are described by Hartree–Fock wavefunctions. The electronic bonding of these molecules is found to be similar to that of comparable alkaline–earth hydrides and fluorides. The uranium 6p, 6d, and 5f orbitals retain their atomic character but the orbitals extend into the bonding region and are distorted by overlap repulsion and electrostatic effects. Nonetheless, the atomic energetic coupling determines that low energy states will have the maximum spin multiplicity and maximum orbital angular momentum projection consonant with the charge-transfer bonding.     

Electronic states of NiFe. Anab initio HF-CI study

Ab initio Hartree-Fock and configuration interaction methods have been employed in describing the interaction between a Ni and an Fe atom. The chemical bond between the atoms is due to a 4sσ molecular orbital. The 3d orbitals merely cause small splittings between the potential energy curves. Equilibrium distance, dissociation energy and vibrational frequency are predicted for the ground state of the molecule.     

Configuration interaction calculations on the 2P ground state of boron atom and C+ using Slater orbitals

Configuration Interaction (CI) calculations on the ground ²P state of boron atom are presented using a wave function expansion constructed with L‐S eigenfunction configurations of s‐, p‐, and d‐Slater orbitals. Two procedures of optimization of the orbital exponents have been investigated. First, CI(SD) calculations including few types of configurations and full optimization of the orbital exponents led to the energy ?24.63704575 a.u. Second, full‐CI (FCI) calculations including a large number of configuration types using a fixed set of orbital exponents for all configurations gave ?24.63405222 a.u. using the basis [4s3p2d] and 2157 configurations, and to an improved result of ?24.64013999 a.u. for 3957 configurations and a [5s4p3d] basis. This last result is better than earlier calculations of Schaefer and Harris (Phys Rev 1968, 167, 67), and compares well with the recent ones from Froese Fischer and Bunge (personal communication). In addition, using the same wave functions, CI calculations of the boron isoelectronic ion C⁺ have been performed obtaining an energy of ?37.41027598 a.u. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Accurate analytical self-consistent field wave functions for Tm2+ and Tm3+

The analytical expansion self-consistent field method was employed to perform ab initio calculations for the ground states of the rare-earth ions, Tm²⁺, 4f¹³, ²F, and Tm³⁺, 4f¹², ³H, (Z = 69). In each case the total number of basis functions used in the analytical expansions was 29, distributed as follows: 10, 8, 5, and 6, for the symmetries s, p, d, and f, respectively. All of the orbital exponents of the basis functions were optimized repeatedly, to the extent of the single-precision computer representation. Values of 〈rⁿ〉 for the 4f orbital of both ions are also presented, for the convenience of experimentalists.     

Analytical optimization of exponent values in protonic and deuteronic Gaussian‐type functions by elimination of translational and rotational motions from multi‐component molecular orbital scheme

To optimize the exponent values in protonic and deuteronic Gaussian‐type functions (GTF) by the elimination of translational and rotational motions, we have proposed the new scheme of an analytical gradient formula with respect to the exponent values in the multi‐component molecular orbital scheme, which can take into account the quantum effects of protons and deuterons, under the Hartree‐Fock level of theory. Numerical assessment of H₂ and D₂ molecules confirms that there is a clear difference between distributions of protonic and deuteronic orbitals following the elimination of translational and rotational motions. In particular, the d‐type GTF in the protonic orbital drastically improves the total energy. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008