Numerical atomic basis sets generation

Linear dependencies of Gaussian-type orbital basis sets employed in the framework of the HF SCF method for periodic structures, which occur when diffuse basis functions are included in a basis set in an uncontrolled manner, were investigated [468]. The basis sets constructed avoid numerical linear dependences and were optimized for a number of periodic structures. The numerical AO basis sets for solids were generated in [469] by confining atoms within spheres and smoothing the orbitals so that the first and second derivatives go to zero at the boundary. This forms small atomic-like basis sets that can be applied to solid-state problems and are efficient for treating large systems. 

The band-structure code, called BAND, also uses STO basis sets with STO fit functions or numerical atomic orbitals. Periodicity can be included in one, two, or three dimensions. No geometry optimization is available for band-structure calculations. The wave function can be decomposed into Mulliken, DOS, PDOS, and COOP plots. Form factors and charge analysis may also be generated. 

Variational one-center restoration. In the variational technique of one-center restoration (VOCR) [79, 80], the proper behavior of the four-component molecular spinors in the core regions of heavy atoms can be restored as an expansion in spherical harmonics inside the sphere with a restoration radius, Rvoa, that should not be smaller than the matching radius, Rc, used at the RECP generation. The outer parts of spinors are treated as frozen after the RECP calculation of a considered molecule. This method enables one to combine the advantages of two well-developed approaches, molecular RECP calculation in a gaussian basis set and atomic-type one-center calculation in numerical basis functions, in the most optimal way. This technique is considered theoretically in [80] and some results concerning the efficiency of the one-center reexpansion of orbitals on another atom can be found in [75]. 

Within the density functional theory (DFT), several schemes for generation of pseudopotentials were developed. Some of them construct pseudopotentials for pseudoorbitals derived from atomic calculations [29] - [31], while the others make use [32] - [36] of parameterized analytical pseudopotentials. In a specific implementation of the numerical integration for solving the DFT one-electron equations, named Discrete-Variational Method (DVM) [37]- [41], one does not need to fit pseudoorbitals or pseudopotentials by any analytical functions, because the matrix elements of an effective Hamiltonian can be computed directly with either analytical or numerical basis set (or a mixed one). 

We carried out non-spin-polarized molecular orbital calculations using a DV-Xa code called SCAT (7). In the DV-Xa method, matrix elements in the secular equation are derived from the weighted sum of integrand values at sampling points. Hence numerical atomic orbitals can be used as basis sets. These atomic orbitals were generated on each iteration in atomic-like potentials derived from spherically averaging the molecular potential around each nucleus (7). The numerical wave functions thus obtained can be used efficiently for the molecular orbitals and have a practical character making them suitable for... 

Although the availability of numerical solutions of HF equations is still restricted to at most two-center (or linear) systems, the development of suitable basis sets enabled the computation of SCF solutions within the Roothaan linear combination of atomic orbitals (LCAO) SCF formalism [9], Generation of such solutions, even for systems with several hundreds of electrons, is no-... 

The generator coordinate method (GCM), as initially formulated in nuclear physics, is briefly described. Emphasis is then given to mathematical aspects and applications to atomic systems. The hydrogen atom Schrodinger equation with a Gaussian trial function is used as a model for former and new analytical, formal and numerical derivations. The discretization technique for the solution of the Hill-Wheeler equation is presented and the generator coordinate Hartree-Fock method and its applications for atoms, molecules, natural orbitals and universal basis sets are reviewed. A connection between the GCM and density functional theory is commented and some initial applications are presented. 

Let us now find out whether these classical enthalpies may be reproduced by electronic-structure calculations (VASP) on Sn/Zn supercells using ultra-soft pseudopotentials, plane-wave basis sets and the GGA. We therefore have to theoretically determine the total energies of all crystal structure types under consideration (a-Sn, j6-Sn, Zn) as a function of the composition SnxZni x by a variation of the available atomic sites in terms of Sn and Zn occupation, just as for the preceding oxynitrides (CoOi- N ). In the present case, supercells with a total of 16 atoms were generated, and nine different compositions per structure were numerically evaluated. Because this amounts to a significant computational task, the use of pseudopotentials is mandatory, and this also allows the rapid calculation of interatomic forces and stresses for structural... 

---