WJCMP  Vol.4 No.3 , August 2014
The Pursuit of Fallacy in Density Functional Theory: The Quest for Exchange and Correlation, the Rigorous Treatment of Exchange in the Kohn-Sham Formalism and the Continuing Search for Correlation
Abstract: As pointed out in the paper preceding this one, in the case of functionals whose independent variable must obey conditions of integral normalization, conventional functional differentiation, defined in terms of an arbitrary test function, is generally inapplicable and functional derivatives with respect to the density must be evaluated through the alternative and widely used limiting procedure based on the Dirac delta function. This leads to the determination of the rate of change of the dependent variable with respect to its independent variable at each isolated pair, , that may not be part of a functional (a set of ordered pairs). This extends the concept of functional derivative to expectation values of operators with respect to wave functions leading to a density even if the wave functions (and expectation values) do not form functionals. This new formulation of functional differentiation forms the basis for the study of the mathematical integrity of a number of concepts in density functional theory (DFT) such as the existence of a universal functional of the density, of orbital-free density functional theory, the derivative discontinuity of the exchange and correlation functional and the extension of DFT to open systems characterized by densities with fractional normalization. It is shown that no universal functional exists but, rather, a universal process based only on the density and independent of the possible existence of a potential, leads to unique functionals of the density determined through the minimization procedure of the constrained search. The mathematical integrity of two methodologies proposed for the treatment of the Coulomb interaction, the self-interaction free method and the optimized effective potential method is examined and the methodologies are compared in terms of numerical calculations. As emerges from this analysis, the optimized effective potential method is found to be numerically approximate but formally invalid, contrary to the rigorously exact results of the self-interaction-free method.
Cite this paper: Gonis, A. (2014) The Pursuit of Fallacy in Density Functional Theory: The Quest for Exchange and Correlation, the Rigorous Treatment of Exchange in the Kohn-Sham Formalism and the Continuing Search for Correlation. World Journal of Condensed Matter Physics, 4, 200-225. doi: 10.4236/wjcmp.2014.43023.

[1]   Hohenberg, P. and Kohn, W. (1964) Inhomogeneous Electron Gas. Physical Review, 136, B864-B871.

[2]   Kohn, W. and Sham, L.J. (1965) Self-Consistent Equations Including Exchange and Correlation Effects. Physical Review, 140, A1133-A1138.

[3]   Parr, R.G. and Yang, C.Y. (1989) Density Functional Theory of Atoms and Molecules. Oxford University Press, Oxford.

[4]   Dreitzler, R.M. and Gross, E.K.U. (1990) Density Functional Theory. Springer Verlag, Berlin, New York.

[5]   Zumbach, G. and Maschke, K. (1983) New Approach to the Calculation of Density Functionals. Physical Review A, 28, 544-554.

[6]   The strict mathematical definition of a functional is that of a function of coordinates mapped onto a number, usually in terms of an integral. The definition used here includes the mapping of a function onto a function, or a set of functions.

[7]   Lieb, E.H. (1982) Density Functionals for Coulomb Systems. In: Shimony, A. and Feshbach, H., Eds., Physics as Natural Philosophy, Essays in Honor of Laszlo Tisza, MIT Press Cambridge, Massachusetts, 111-149.
See also, Lieb, E.H. (1985) Density Functional Methods in Physics. In: Dreizler, R.M., Ed., NATO ASI Series B 123, Plenum, New York, 31, and Lieb, E.H. (1983) in Int. J. Quantum Chem. 24, 243.

[8]   Levy, M. (1979) Universal Variational Functionals of Electron Densities, First-Order Density Matrices, and Natural Spin-Orbitals and Solution of the v-Representability Problem. Proceedings of the National Academy of Sciences of the United States of America, 76, 6062-6065.

[9]   It is occasionally argued [3] that the variational principle can be expressed in terms of the variational procedure, allowing, in principle, a fractionally normalized charge density. Here, the quantity, , is a Lagrange multiplier that guarantees the presence of a given normalization associated with wave functions not normalized to unity in terms of which the variational theorem takes the form. These considerations are both unnecessary as well as frought with formal danger. First, the wave functions considered with respect of the variational theorem can always be normalized to unity through division by the square root of its squared modulus. Second, a variational procedure based on a fractionally normalized density can easily misdirect developments along the treatment of ensembles of open systems in terms of formalism appropriate exclusively for pure states, as is discussed explicitly in the body of the paper. As a quick reminder, the wave functions used in the variational theorem must be of the form emerging as solutions to the Schr¨odinger equation, and those are normalized to unity.

[10]   Cioslowski, J. (1988) Density Functionals for the Energy of Electronic Systems: Explicit Variational Construction. Physical Review Letters, 60, 2141-2143.

[11]   Cioslowski, J. (1988) Density Driven Self-Consistent Field Method. I. Derivation and Basic Properties. The Journal of Chemical Physics, 89, 4871-4874.

[12]   Cioslowski, J. (1989) Density Driven Self-Consistent Field Method. II. Construction of All One-Particle Wave Functions that Are Orthonormal and Sum up to a Given Density. Quantum Chemistry Symposium. International Journal of Quantum Chemistry, 36, 255-262.

[13]   It is useful to use a function defined through a process to demonstrate the difficulties that can arise with respect to functionals defined by means of a procedure, rather than as functions of the independent variable. Consider the function, where is a real number, and is determined through the process of flipping a coin for any given, setting if the coin comes up heads, and if it comes up tales. Unlike a functional, however, this function of process has no parametric derivative (thus no derivative).

[14]   Some exceptions are known. The ground-state wave function of a single particle system is given in terms of the density, and the corresponding kinetic energy functional takes the von Weizsacker form. Also the kinetic energy and exchange energy of an infinite number of plane waves forming a density are given explicitly by means of the expression. These cases are the exceptions that prove the rule that functionals of the density defy representation as explicit functions of the density.

[15]   Becke, A.D. (1993) A New Mixing of Hartree-Fock and Local Density-Functional Theories. The Journal of Chemical Physics, 98, 1372.

[16]   Perdew, J.P., Ernzerhof, M. and Burke, K. (1996) Rationale for Mixing Exact Exchange with Density Functional Approximations. The Journal of Chemical Physics, 105, 9982.

[17]   Kim, K. and Jordan, K.D. (1994) Comparison of Density Functional and MP2 Calculations on the Water Monomer and dimer. The Journal of Physical Chemistry, 98, 10089-10094.

[18]   Becke, A.D. (1988) Density-Functional Exchange-Energy Approximation with Correct Asymptotic Behavior. Physical Review A, 38, 3098-3100.

[19]   Lee, C., Yang, W.T. and Parr, R.G. (1988) Development of the Colle-Salvetti Correlation-Energy Formula into a Functional of the Electron Density. Physical Review B, 37, 785-789.

[20]   Becke, A.D. (1993) Density-Functional Thermochemistry. III. The Role of Exact Exchange. The Journal of Chemical Physics, 98, 5648.

[21]   Heyd, J., Scuseria, G.E. and Ernzerhof, M. (2003) Hybrid Functionals Based on a Screened Coulomb Potential. The Journal of Chemical Physics, 118, 8207.

[22]   Zhao, Y. and Truhlar, D.G. (2008) The M06 Suite of Density Functionals for Main Group Thermochemistry, Thermochemical Kinetics, Noncovalent Interactions, Excited States, and Transition Elements: Two New Functionals and Systematic Testing of Four M06-Class Functionals and 12 Other Functionals. Theoretical Chemistry Accounts, 120, 215241.

[23]   Zhao, Y. and Truhlar, D.G. (2006) Density Functional for Spectroscopy: No Long-Range Self-Interaction Error, Good Performance for Rydberg and Chargetransfer States, and Better Performance on Average than B3LYP for Ground States. The Journal of Physical Chemistry A, 110, 13126-13130.

[24]   Kümmel, S. and Kronik, L. (2008) Orbital-Dependent Density Functionals: Theory and Applications. Reviews of Modern Physics, 80, 3-60.

[25]   Dane, M., Gonis, A., Nicholson, D.M. and Stocks, G.M. (2013) On the Solution of the Self-Interaction Problem in Kohn-Sham Density Functional Theory. Journal of Physics and Chemistry of Solids (Online Version), 75, 1160-1178.

[26]   This situation is not unlike that of the function defined over the set of natural numbers that exhibits a minimum but no derivative at the point of the minimum, e.g., defined over the integers and zero, has a minimum at 0, but no derivative there.

[27]   Trott, M. Functional Derivative. From MathWorld-A Wolfram Web Resource, Created by Eric W. Weisstein.

[28]   English, H. and English, R. (1984) Exact Density Functional for Ground-State Energies I. General Results. Physica Status Solidi (b), 123, 373.

[29]   English, H. and English, R. (1984) Exact Density Functional for Ground-State Energies II. Details and Remarks. Physica Status Solidi (b), 124, 711.

[30]   Perdew, J.P., Parr, R.G., Levy, M. and Balduz, J.L. (1982) Density-Functional Theory for Fractional Particle Number: Derivative Discontinuities of the Energy. Physical Review Letters, 49, 1691-1694.

[31]   Indeed, one is hard put to find a case in quantum chemistry or condensed matter physics where functional differentiation is carried out in terms other that based on the Dirac delta function used as the arbitrary test function.

[32]   Fock, V. (1930) Naherungsmethode zur Losung des quantenmechanischen Mehrkorperproblems. Zeitschrift für Physik, 61, 126-148.

[33]   Fermi, E. and Amaldi, E. (1934) Le orbite 8s degli elementi. Mem. Accad. d’Italia, 6, 119-149, number 82 in Refs. [34], [35].

[34]   Fermi, E. and Segre, E. (2011) Collected Papers, Note E Memorie, of Enrico Fermi V1: Italy, 1921-1938. Literary Licensing, LLC, Whitefish, MT.

[35]   Fermi, E. (1962) Enrico Fermi-Collected Papers (Note E Memorie), 1921-1938, Volume 1. The University of Chicago Press, Chicago and London; Accademia Nationale dei Lincei, Rome.

[36]   LibXC.

[37]   Perdew, J.P. and Zunger, A. (1981) Self-Interaction Correction to Density-Functional Approximations for Many-Electron Systems. Physical Review B, 23, 5048-5079.

[38]   Kümmel, S. and Perdew, J.P. (2003) Simple Iterative Construction of the Optimized Effective Potential for Orbital Functionals, Including Exact Exchange. Physical Review Letters, 90, Article ID: 043004.

[39]   Grabo, T., Kreibich, T., Kurth, S. and Gross, E. (2000) Strong Coulomb Correlations in Electronic Structure Calculations: Beyond the Local Density Approximation. Gordon and Breach Science Publishers, Amsterdam, 203-311.

[40]   Yang, W.T. and Wu, Q. (2002) Direct Method for Optimized Effective Potentials in Density-Functional Theory. Physical Review Letters, 89, Article ID: 143002.

[41]   The reader unfamiliar with the mathematics of the OEP method may benefit by reading standard treatises, e.g., [24].

[42]   Yang, W.T., Ayers, P.W. and Wu, Q. (2004) Potential Functionals: Dual to Density Functionals and Solution to the v-Representability Problem. Physical Review Letters, 92, 146404.

[43]   Gidopoulos, N.I. and Lathiotakis, N.N. (2012) Constraining Density Functional Approximations to Yield Self-Interaction Free Potentials. The Journal of Chemical Physics, 136, Article ID: 224109.

[44]   Heaton-Burgess, T., Bulat, F.A. and Yang, W.T. (2007) Optimized Effective Potentials in Finite Basis Sets. Physical Review Letters, 98, Article ID: 256401.

[45]   Harbola, M.K. and Sahni, V. (1989) Quantum-Mechanical Interpretation of the Exchange-Correlation Potential of Kohn-Sham Density Functional Theory. Physical Review Letters, 62, 489-492.

[46]   Ceperley, D.M. and Alder, B.J. (1980) Ground State of the Electron Gas by a Stochastic Method. Physical Review Letters, 45, 566-569.

[47]   Gal, T. and Nagy, A. (2000) A Method to Get an Analytical Expression for the Non-Interacting Kinetic Energy Density Functional. Journal of Molecular Structure, 501-502, 167-171.

[48]   Wang, Y. and Nagy, A. (1993) Construction of Exact Kohn-Sham Orbitals from a Given Electron Density. Physical Review A, 47, R1591-R1593.

[49]   Parr, R.G., Kugler, A.A. and Nagy, A. (1995) Some Identities in Density-Functional Theory. Physical Review A, 52, 969-976.

[50]   Levy, M., Perdew, J.P. and Sahni, V. (1984) Exact Differential Equation for the Density and Ionization Energy of a Many-Particle System. Physical Review A, 30, 2745-2748.

[51]   Wang, A. and Carter, E.A. (2000) Chapter 5 of “Theoretical Methods in Condensed Phase Chemistry”. In: Schwartz, S.D., Ed., A New Book Series of “Progress in Theoretical Chemistry and Physics”, Kluwer, Dordrecht, 117-184.

[52]   Garcia-Gonzales, P., Alvarellos, J.E. and Chacon, E. Kinetic-Energy Density Functional: Atoms and Shell Structure. Physical Review A, 54, 1897.

[53]   Mirtschink, A., Seidl, M. and Gori-Giorgi, P. (2013) Derivative Discontinuity in the Strong-Interaction Limit of Density-Functional Theory. Physical Review Letters, 111, Article ID: 126402.

[54]   Landau, L.D. and Lifschitz, E.M. (1977) Quantum Mechanics: Non-Relativistic Theory. 3rd Edition, Elsevier Butterworth-Heinemann Linacre House, Jordan Hill, Oxford, Burlington.