Exchange-correlation potential method

A common problem for both methods lies in the use of potentials that do not possess the correct net attractiveness. This can have the consequence that continuum feamres appear shifted in energy. In particular, there is evidence that the LB94 exchange-correlation potential currently used for the B-spline calculations, although possessing the correct asymptotic behavior for ion plus electron, is too attractive, and near threshold features can then disappear below the ionization threshold. An empirical correction can be made, offsetting the energy scale, but this can mean that dynamics within a few electronvolts of threshold get an inadequate description or are lost. There is limited scope to tune the Xa potential, principally by adjustment of the assumed a parameter, but for the B-spline method a preferable alternative for the future may well be use of the SAOP functional that also has correct asymptotic behavior, but appears to be better calibrated for such problems [79]. 

FIGURE 7.1 Exchange-correlation potential Vxc (in atomic units) for neon as a function of distance r (in atomic units) from the nucleus. The potential is obtained from the ground-state density by employing the ZP method. 

To perform excited-state calculations, one has to approximate the exchange-correlation potential. Local self-interaction-free approximate exchange-correlation potentials have been proposed for this purpose [73]. We can try to construct these functionals as orbital-dependent functionals. There are different exchange-correlation functionals for the different excited states, and we suppose that the difference between the excited-state functionals can be adequately modeled through the occupation numbers (i.e., the electron configuration). Both the OPM and the KLI methods have been generalized for degenerate excited states [37,40]. 

Table 9.1 presents excitation energies for a few atoms and ions. Calculations were performed with the generalized KLI approximation [69,74], For comparison, experimental data and the results obtained with the local-spin-density (LSD) exchange-correlation potential [75] are shown. The KLI method contains only the exchange. 

The long-range asymptotic form of the exchange potential, a constituent of the exchange-correlation potential (56), will be discussed in this and the remaining subsections. Since the known derivations of this form are very complicated and involve an analysis of an integral equation of the OP method (see, e.g. [17], [27],... 

The line-integral expression (189), for the exchange-correlation potential of the HF-KS approach in Sect. 2.4, offers an interesting way to reconstruct the exchange potential for a given system, from the known HF solution for this system. Being alternative to schemes discussed in Sect. 3, this method provides an expression for the exchange potential solely in terms of the HF orbitals in the form... 

The basic concepts of the one-electron Kohn-Sham theory have been presented and the structure, properties and approximations of the Kohn-Sham exchange-correlation potential have been overviewed. The discussion has been focused on the most recent developments in the theory, such as the construction of from the correlated densities, the methods to obtain total energy and energy differences from the potential, and the orbital dependent approximations to v. The recent achievements in analysis of the atomic shell and molecular bond midpoint structure of have been... 

The difficulty of this problem can be appreciated by noticing that in order to solve the Kohn-Sham equations exactly, one must have the exact exchange-correlation potential which, moreover, must be obtained from the exact exchange-correlation functional c[p( )] given by Eq. (160). As this functional is not known, the attempts to obtain a direct solution to the Kohn-Sham equations have had to rely on the use of approximate exchange-correlation functionals. This approximate direct method, however, does not satisfy the requirement of functional iV-representability,... 

In this paper we present preliminary results of an ab-initio study of quantum diffusion in the crystalline a-AlMnSi phase. The number of atoms in the unit cell (138) is sufficiently small to permit computation with the ab-initio Linearized Muffin Tin Orbitals (LMTO) method and provides us a good starting model. Within the Density Functional Theory (DFT) [15,16], this approach has still limitations due to the Local Density Approximation (LDA) for the exchange-correlation potential treatment of electron correlations and due to the approximation in the solution of the Schrodinger equation as explained in next section. However, we believe that this starting point is much better than simplified parametrized tight-binding like s-band models. 

In spite of the current popularity of density functional methods and the many efforts to construct functionals that accurately describe molecular electronic properties, an exact exchange-correlation potential expressed only in terms of the electron density is elusive. This paper presents some thoughts around this problem that try to reflect some of the concerns about density functionals expressed by Per-Olov Lowdin on many occasions. 

The construction of exchange correlation potentials and energies becomes a task for which not much guidance can be obtained from fundamental theory. The form of dependence on the electron density is generally not known and can only to a limited extent be obtained from theoretical considerations. The best one can do is to assume some functional dependence on the density with parameters to satisfy some consistency criteria and to fit calculated results to some model systems for which applications of proper quantum mechanical theory can be used as comparisons. At best this results in some form of ad-hoc semi-empirical method, which may be used with success for simulations of molecular ground state properties, but is certainly not universal. 

Clever use of this approach and choice of exchange correlation potential has produced methods for calculating molecular ground state properties and energetics in reasonable agreement with experimental numbers for an impressive array of systems. The fact that these methods are computationally inexpensive make them applicable to larger systems than can routinely be studied with proper quantum mechanical approaches. This makes them useful simulation tools. 

---