On the mathematical formulation and interpretation of LACO MO theory

The mathematical basis of LCAO MO theory is studied, both within the Hartree-Fock approximation and in more exact formulations. The basic LCAO expansion for molecular orbitals ¦> in terms of atomic orbitals ¦x> is conveniently written ¦> = ¦x> S ^–1 B where S is the overlap matrix for atomic orbitals and B is the matrix of atomic orbital-molecular orbital overlaps. It is suggested that matrices P and Q, defined by P=B B and Q=BnB where n is the matrix of molecular orbital occupation numbers, are appropriate to the interpretation of molecular calculations in terms of atomic orbital components, electronic populations and the degree of bonding. Implications for Hartree-Fock calculations are investigated.     

Simplified accurate approximation for the Kohn-Sham potential using the KLI method

The method of Krieger, Li, and Iafrate (KLI) [Phys. Rev. A46, 5453 (1992) and A47, 165 (1993)] is employed to calculate the Kohn-Sham (KS) potential, V_κσ, for the exchange-only case in which the electron-electron interaction between “core” electrons in the Hartree-Fock exchange energy functional is treated in the local-spin-density (LSD) approximation with and without self-interaction-correction (SIC). The resulting V_κσ(r) maintains the important analytic properties exhibited by the exact KS potential. When the core is taken to include all occupied states except those in the last two occupied subshells of the atom, we find that properties strongly dependent on the valence electron states continue to be accurately approximated. In particular, when the LSDSIC approximation is employed, we find the results of self-consistent calculations of the ionization potential and electron affinity are within 0.3 mRy of the exact KS results and that the energy eigenvalue corresponding to the highest energy occupied orbital and <r²> have an average error of a few tenths of 1% for both atoms and negative ions for Z ≤ 20. Similarly, slightly less accurate results are obtained when the LSD approximation is employed. These results suggest that the KLI method may be accurately and more easily applied to multiatom systems when this additional approximation is made. © 1997 John Wiley & Sons, Inc.     

Non-additivity in water-ion-water interactions

The three-body system Li⁺(H₂O)₂ was analyzed to study that non-additive part of the interaction potential which can be obtained by the Hartree-Fock approximation.For long and intermediate distances the three-body correction was found to be well represented by the induction energy, where bond dipoles are induced on each water molecule by point charges located on the (unpolarizable) lithium ion and on the other molecule respectively: for shorter distances this approximation was corrected by means of an exponential repulsive term. Such a potential model for non-additive interactions was extended to the more general situation Li⁺(H₂O)_n, and Monte-Carlo calculations were carried out on clusters containing up to six water molecules; comparison with other simulation results and with available data showed a significantly improved agreement with experiment. Tentative values for H are presented for n =7, 8,..., 20, where experimental data are not available.     

Improving numerical integration through basis set expansion

Calculations are presented to assess a theorem presented by S.F. Boys [(1969) Proc. R. Soc. A. 309:195], regarding the accuracy of numerical integration in quantum chemical calculations. The theorem states that the error due to numerical integration can be made proportional to the error due to basis set truncation, and thus goes to zero in the limit of a complete basis. We test this theorem on the hydrogen atom, showing that with a solution-spanning basis, the numerically exact orbital energy can indeed be calculated with a small number of integration points. Moreover, tests for H and H₂⁺ demonstrate that even when only a near-complete basis is employed, Boys Theorem can significantly reduce integration error. However, for other systems, like the oxygen atom and the CO₂ molecule, the theorem yields no advantage for some occupied orbitals. It is concluded that the theorem would be most useful for calculations that demand large basis sets.     

Natural transition orbitals for complex two-component excited state calculations

While the natural transition orbital (NTO) method has allowed electronic excitations from time-dependent Hartree-Fock and density functional theory to be viewed in a traditional orbital picture, the extension to multicomponent molecular orbitals such as those used in relativistic two-component methods or generalized Hartree-Fock (GHF) or generalized Kohn-Sham (GKS) is less straightforward due to mixing of spin-components and the inherent inclusion of spin-flip transitions in time-dependent GHF/GKS. An extension of single-component NTOs to the two-component framework is presented, in addition to a brief discussion of the practical aspects of visualizing two-component complex orbitals. Unlike the single-component analog, the method explicitly describes the spin and frequently obtains solutions with several significant orbital pairs. The method is presented using calculations on a mercury atom and a CrO₂Cl₂ complex.     

Local physical quantities for spin based on the relativistic quantum field theory in molecular systems

Local physical quantities for spin are investigated on the basis of the four‐ and two‐component relativistic quantum theory. In the quantum field theory, local physical quantities for spin such as the spin angular momentum density, spin torque density, zeta force density, and zeta potential play important roles in spin dynamics. We discuss how to calculate these local physical quantities based on the two‐component relativistic quantum theory. Some different types of relativistic numerical calculations of local physical quantities in Li atom and C₆H₆ are demonstrated and compared. Local physical quantities for each orbital are also discussed, and it is seen that a total local zeta potential is given as a result of some cancellation of large contributions from each orbital. © 2016 Wiley Periodicals, Inc.     

Ab initio Studies on polymers

Ab initio crystal orbital calculations have been performed on the infinite all-trans polyene. A structure with r _c=c= 1.346 Å, r _c-c = 1.446 Å, r _c-h = 1.08 Å, and CCC = 125.3 ° was found to be most stable. The most important force constants, the band structure and the density of states were determined as well.     

Additivity model for calculations of UHF spin densities in some aza-aromatic radical anions

The concept of a universal basis set for electronic structure calculation is illustrated by presenting results obtained when basis sets are transferred from one atom to another. A single Slater-orbital basis set, consisting of nine 1s and six 2p functions, produces Hartree-Fock total energies and orbital energies in good agreement with the most accurate calculations of these energies obtained using different basis sets individually optimized for each atom. Transferability of integrals is a natural consequence of the use of the same basis set for each atom in a molecule.     

Numerical integration of overlap and ligand contributions to the electric field gradient

The total electric field gradient (EFG) tensor V _pq is calculated by numerical integration of threedimensional integrals. Each of them is solved a) by integrating over one dimension analytically and b) by integrating over the remaining two dimensions on the basis of a Gauss-type integration rule. The use of 100 abscissas in the twodimensional numerical integration scheme yields satisfactory accuracy which was checked by evaluating overlap integrals; an increase to 400 abscissas does not increase the result drastically. Calculating quadrupole splittings E _Q from numerically integrated electric field gradient tensors V _pq we observe that depending a) on the amount of covalency and b) on the amount of deviation from octahedral or tetrahedral symmetry, involved in a molecular system, overlap and ligand contributions to V _pq play an important role. Especially for the sandwich compound ferrocene, Fe(C₅H₅)₂, we find a significant difference between E _Q ^{num. int.} which follows from the numerical integration method, and E _Q ^conventional which is derived from effective charges.Supported by Deutsche Forschungsgemeinschaft     

GRECP/5e‐MRD‐CI calculation of the electronic structure of PbH

The correlation calculation of the electronic structure of PbH is carried out with the generalized relativistic effective core potential (GRECP) and multireference single‐ and double‐excitation configuration interaction (MRD‐CI) methods. The 22‐electron GRECP for Pb is used and the outer core 5s, 5p, and 5d pseudospinors are frozen using the level‐shift technique, so only five external electrons of PbH are correlated. A new configuration selection scheme with respect to the relativistic multireference states is employed in the framework of the MRD‐CI method. The [6, 4, 3, 2] correlation spin–orbit basis set is optimized in the coupled cluster calculations on the Pb atom using a recently proposed procedure, in which functions in the spin–orbital basis set are generated from calculations of different ionic states of the Pb atom and those functions are considered optimal that provide the stationary point for some energy functional. Spectroscopic constants for the two lowest‐lying electronic states of PbH (²Π_1/2, ²Π_3/2) are found to be in good agreement with the experimental data. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002     

Accurate local spin-polarized exchange potential: Reconciliation of generalized Slater and Kohn–Sham methods

Using simple physical arguments, a local spin-polarized exchange potential, V_xσ, is constructed from the single-particle Hartree–Fock (HF ) potentials (generalized Slater method) that reduces to the usual Kohn–Sham (KS ) result in the uniform gas limit. Numerical results for 10 closed subshell atoms demonstrate that the total energy calculated employing this V_xσ is closer to the exact KS results than those of other standard exchange approximations with electron densities and highest occupied orbital eigenvalues that closely approximate the HF results.     

<Emphasis Type="Italic">Ab initio</Emphasis> structural study of M(mda)<Subscript>2</Subscript> complexes (M = Be,Mg, Ca; mda = C<Subscript>3</Subscript>O<Subscript>2</Subscript>H<Subscript>3</Subscript>)

The structure of M(mda)₂ (M = Be, Mg, Ca; mda = C₃O₂H₃) bis-complexes was investigated by the ab initio Hartree-Fock method and by including electron correlation in terms of second order Möller-Plesset perturbation theory; for calculations we used triple-zeta valence basis sets complemented with polarization functions. Two most probable geometrical nuclear configurations (D _2h and D _2d ) are considered for each molecule. The structure with two mutually orthogonal chelate ligands (D _2d symmetry) corresponds to the potential energy surface (PES) minimum. The planar D _2h configuration corresponds to the first order saddle point on PES; consequently, its relative energy determines the height of the barrier to the D _2d D _2h D _2d intramolecular rearrangement. Correlation equations that relate the calculated values of equilibrium internuclear distances, force constants, and rearrangement barrier heights to the value of the ionic radius of the metal atom have been obtained. These correlations were employed to evaluate the molecular constants for Sr(mda)₂ and Ba(mda)₂. The theoretical data are compared with the available experimental literature data.Original Russian Text Copyright © 2004 by V. V. Sliznev, S. B. Lapshina, and G. V. GirichevTranslated from Zhurnal Strukturnoi Khimii, Vol. 45, No. 4, pp. 611–623, July–August, 2004This revised version was published online in April 2005 with a corrected cover date.     

Some developments in the spin-extended Hartree-Fock method

All the second-order density matrix spin components for the spin-extended Hartree-Fock method are obtained. The coefficients in the final formulae are only ω_sM, ω_sM±1, ω_sM±2, where ω_sM are the weights of pure states of spin s in the initial unprojected determinant with spin projection M. The eigenvalue problem for the best electron density natural orbitals in the spin-extended method is formulated. All the second-order transition density matrix spin components between pure spin basis functions built of orthogonal orbitals and distinguished by different core choice are also found. This basis may be used on CI calculations.     

Semiempirical studies of core-electron binding energy shifts

The inner-core binding-energy shifts (BEs) of boron and carbon atoms in various chemical environments were studied by the semiempirical Self-Consistent Charge Molecular Orbital (SCC MO) method. The calculations are based on the initial ground state electrostatic potential model. The main feature of our approach is the empirical treatment of the coefficient relating BEs with the orbital populations of the host atom and the Madelung energy term. These adjustable parameters absorb a large portion of relaxation energy. The so obtained results are in good agreement with experimental data. They are better than earlier CNDO/2 results obtained by using either ground state or relaxation potential models. Present results indicate that semiempirical methods like SCC MO are able to account for changes in BE(1s) with a fair accuracy although the inner-shell electrons are not explicitly considered in the actual calculations.     

Organomineral sorbents based on clinoptilolite-containing tuffs

The anion-exchange properties of a new organomineral sorbent obtained by modification of clinoptilolite-containing tuffs by polyhexamethyleneguanidine have been studied after different periods of storage (time after synthesis) and numbers of sorption-regeneration cycles. The sorbent can be used as a cation- and an anion-exchanger simultaneously. Selectivity coefficients (exchange constants) for F^–, SO₄ ^2–, and HPO₄ ^2– ionsvs. Cl^– ions on modified clinoptilolite-containing tuffs have been determined. The modification improves the mechanical properties of clinoptilolite tuffs.For part 1, see Ref. 1.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 9, pp. 1554–1556, September, 1994.     

Radial dependence of exterior electron distributions of molecular orbitals

A molecular surface is introduced to divide interior electron densities from exterior electron densities (EED). The radial distribution of EED (RADEED) is defined for each molecular orbital as a function of the distance from the molecular surface. Logarithmic plots of RADEED for NH₃ using various basis sets in ab initio MO calculations revealed some important features: (i) the Hartree-Fock limit for the orbital function tail may be suggested and thus qualities of basis sets can be discussed, and (ii) the slope of the curve shows the decay rate of the orbital which can be compared with the curve derived from the theoretical behavior of the long-range asymptotic form involving either the lowest ionization potential or the orbital energy of the highest occupied orbital.Dedicated to Professor J. Koutecký on the occasion of his 65th birthday     

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.     

Self consistent modified extended Hückel (SC-MEH) calculations on heavy metal systems. I. Platinum(II) tetragonal planar complexes with and without relativistic effects

The Self Consistent Modified Extended Hückel molecular orbital method had been applied to several square planar complexes of platinum (II). Calculations including both the limited 5d, 6s, 6p and extended 5s, 5p, 5d, 6s, 6p starting bases for platinum were made. It is shown that in PtCl ₄ ^2– both the nuclear quadrupole moment and minimum total energy vs. bond distance are calculated to be in good agreement with experiment, only with the extended platinum AO basis.Specific inclusion of relativistic parameters via a pseudo-relativistic approximation are shown to have a significant effect on the energy molecular energy levels, however no meaningful rationalization can be made without the simultaneous inclusion of ligand field parameters as well.Supported in part by a grant made available through the Cancer Association of Greater New Orleans. Use of the facilities of the Computer Research Center of the University of New Orleans is gratefully acknowledged.     

Ab initio Calculation of the charge distribution and the ligand field splitting in the tetrahedral halo complexes CuCl 4 2− 4 and NiCl 4 2−

The charge distribution and the ligand field splitting in the tetrachloro complexes CuCl ₄ ^2– and NiCl ₄ ^2– have been investigated by means of the restricted Hartree-Fock method. A rather large basis set of contracted Gaussian type orbitals has been employed. The charge distributions have been analysed by means of Mulliken population analyses. The ligand field splitting 10Dq has been compared with literature results known for the octahedral cluster NiF ₆ ^4– occurring in KNiF₃. A detailed analysis has been carried out for CuCl ₄ ^2– . From calculations on a selected number of states of NiCl ₄ ^2– the Racah parameters B and C have been obtained.