Local potential function, kinetic energy

This argument shows that the locality hypothesis fails for more than two electrons because the assumed Frechet derivative must be generalized to a Gateaux derivative, equivalent in the context of OEL equations to a linear operator that acts on orbital wave functions. The conclusion is that the use by Kohn and Sham of Schrodinger s operator t is variationally correct, but no equivalent Thomas-Fermi theory exists for more than two electrons. Empirical evidence (atomic shell structure, chemical binding) supports the Kohn-Sham choice of the nonlocal kinetic energy operator, in comparison with Thomas-Fermi theory [288]. A further implication is that if an explicit approximate local density functional Exc is postulated, as in the local-density approximation (LDA) [205], the resulting Kohn-Sham theory is variation-ally correct. Typically, for Exc = f exc(p)p d3r, the density functional derivative is a Frechet derivative, the local potential function vxc = exc + p dexc/dp. 

Calculations were done with a full-potential version of the LMTO method with nonoverlapping spheres. The contributions from the interstitial region were accounted for by expanding the products of Hankel functions in a series of atom-ce- -ered Hankels of three different kinetic energies. The corrected tetrahedron method was used for Brillouin zone integration. Electronic exchange and correlation contributions to the total energy were obtained from the local-density functional calculated by Ceperley and Alder " and parametrized by Vosko, Wilk, and Nusair. ... 

However, one feature of the HF potential is that it is not a local potential. In the case of perfect data (i.e. zero experimental error), the fitted orbitals obtained are no longer Kohn-Sham orbitals, as they would have been if a local potential (for example, the local exchange approximation [27]) had been used. Since the fitted orbitals can be described as orbitals which minimise the HF energy and are constrained produce the real density , they are obviously quite closely related to the Kohn-Sham orbitals, which are orbitals which minimise the kinetic energy and produce the real density . In fact, Levy [16] has already considered these kind of orbitals within the context of hybrid density functional theories. 

The obstacle to simultaneous quantum chemistry and quantum nuclear dynamics is apparent in Eqs. (2.16a)-(2.16c). At each time step, the propagation of the complex coefficients, Eq. (2.11), requires the calculation of diagonal and off-diagonal matrix elements of the Hamiltonian. These matrix elements are to be calculated for each pair of nuclear basis functions. In the case of ab initio quantum dynamics, the potential energy surfaces are known only locally, and therefore the calculation of these matrix elements (even for a single pair of basis functions) poses a numerical difficulty, and severe approximations have to be made. These approximations are discussed in detail in Section II.D. In the case of analytic PESs it is sometimes possible to evaluate these multidimensional integrals analytically. In either case (analytic or ab initio) the matrix elements of the nuclear kinetic energy... 

Since this is a contradiction, it follows that p(r) p (r). Thus, there is a one-to-one correspondence between the local potential F(r) and the electron density p(r). This implies that V, T, and are uniquely determined by the electron density, and therefore are functionals of the electron density. If is a unique functional of the electron density p, the kinetic and exchange-correlation energies, T and Exc must be functionals of p also. 

Here Ho is the kinetic energy operator of valence electrons Vps is the pseudopotential [40,41] which defines the atomic core. V = eUn(r) is the Hartree energy which satisfies the Poisson equation ArUn(r) = —4nep(r) with proper boundary conditions as discussed in the previous subsection. The last term is the exchange-correlation potential Vxc [p which is a functional of the density. Many forms of 14c exist and we use the simplest one which is the local density approximation [42] (LDA). One may also consider the generalized gradient approximation (GGA) [43,44] which can be implemented for transport calculations without too much difficulty [45]. Importantly a self-consistent solution of Eq. (2) is necessary because Hks is a functional of the charge density p. One constructs p from the KS states Ts, p(r) = (r p r) = ns Fs(r) 2, where p is the density matrix,... 

For direct Af-electron variational methods, the computational effort increases so rapidly with increasing N that alternative simplified methods must be used for calculations of the electronic structure of large molecules and solids. Especially for calculations of the electronic energy levels of solids (energy-band structure), the methodology of choice is that of independent-electron models, usually in the framework of density functional theory [189, 321, 90], When restricted to local potentials, as in the local-density approximation (LDA), this is a valid variational theory for any A-electron system. It can readily be applied to heavy atoms by relativistic or semirelativistic modification of the kinetic energy operator in the orbital Kohn-Sham equations [229, 384],... 

There is still a very important area of contact between LS-DFT and other approaches having to do with the direct evaluation of the Kohn-Sham potentials [6,62-66] from known densities. In this respect, the work, of Zhao, Morrison and Parr [66] is of particular interest as it provides an alternative to local-scaling transformations for a fixed-density variation of the kinetic energy functional. 

---