Exchange correlation basis

Integrals involving the exchange-correlation potential r,c or the exchange-correlation energy density cannot be evaluated analytically so that further sets of auxiliary functions are introduced. (In practice and 6,<, behave similarly so that a common set is used to fit both functions.) The exchange-correlation basis (XCB) also consists of Hermite Gaussians... 

The density is computed as p(r) = 2. n i ). (/ )p. Often, p(r) is expanded in an AO basis, which need not be the same as the basis used for the and the expansion coefficients of p are computed in tenns of those of the It is also connnon to use an AO basis to expand p (r) which, together with p, is needed to evaluate the exchange-correlation fiinctionaTs contribution toCg. 

There is no systematic way in which the exchange correlation functional Vxc[F] can be systematically improved in standard HF-LCAO theory, we can improve on the model by increasing the accuracy of the basis set, doing configuration interaction or MPn calculations. What we have to do in density functional theory is to start from a model for which there is an exact solution, and this model is the uniform electron gas. Parr and Yang (1989) write... 

Again the set of fitting functions may or may not be the same as the orbital and/or the density basis functions. Once the potential has been fitted, the exchange—correlation energy may be evaluated from integrals involving three functions, analogously to eq. 

We have used the basis set of the Linear-Muffin-Tin-Orbital (LMTO) method in the atomic sphere approximation (ASA). The LMTO-ASA is based on the work of Andersen and co-workers and the combined technique allows us to treat all phases on equal footing. To treat itinerant magnetism we have employed the Vosko-Wilk-Nusair parametrization for the exchange-correlation energy density and potential. In conjunction with this we have treated the alloying effects for random and partially ordered phases with a multisublattice generalization of the coherent potential approximation (CPA). 

Note again the formal simplicity of equation (7-17) as compared to equation (7-18) in spite of the fact that the former is exact provided the correct Vxc is inserted, while the latter is inherently an approximation. The calculation of the formally L2/2 one-electron integrals contained in hllv, equation (7-13) is a fairly simple task compared to the determination of the classical Coulomb and the exchange-correlation contributions. However, before we turn to the question, how to deal with the Coulomb and Vxc integrals, we want to discuss what kind of basis functions are nowadays used in equation (7-4) to express the Kohn-Sham orbitals. 

In words, the integral of equation (7-33) for the exchange-correlation potential is approximated by a sum of P terms. Each of these is computed as the product of the numerical values of the basis functions and rp, with the exchange-correlation potential Vxc at each point rp on the grid. Each product is further weighted by the factor Wp, whose value depends on the actual numerical technique used. 

Rashin et al.45 obtained the dipole moment of 32 molecules of biological relevance by means of the DFT(SVWN) and DFT(B88/P86) calculations. The results showed a rather weak dependence of calculated dipole moments on the functional form of the exchange-correlation functional but a strong dependence on the basis set. 

---