arXiv:2108.11534

Lies My Teacher Told Me About Density Functional Theory: Seeing Through Them with the Hubbard Dimer

Kieron Burke, John Kozlowski

Published Aug 26, 2021• 3 citations

Abstract

Most realistic calculations of moderately correlated materials begin with a ground-state density functional theory (DFT) calculation. While Kohn-Sham DFT is used in about 40,000 scientific papers each year, the fundamental underpinnings are not widely appreciated. In this chapter, we analyze the inherent characteristics of DFT in their simplest form, using the asymmetric Hubbard dimer as an illustrative model. We begin by working through the core tenets of DFT, explaining what the exact ground-state density functional yields and does not yield. Given the relative simplicity of the system, almost all properties of the exact exchange-correlation functional are readily visualized and plotted. Key concepts include the Kohn-Sham scheme, the behavior of the XC potential as correlations become very strong, the derivative discontinuity and the difference between KS gaps and true charge gaps, and how to extract optical excitations using time-dependent DFT. By the end of this text and accompanying exercises, the reader will improve their ability to both explain and visualize the concepts of DFT, as well as better understand where others may go wrong.

DOI: 10.48550/arxiv.2108.11534

Related papers

The Hubbard dimer: a density functional case study of a many-body problem

D J Carrascal, J Ferrer, J C Smith +1 more

This review explains the relationship between density functional theory and strongly correlated models using the simplest possible example, the two-site Hubbard model. The relationship to traditional quantum chemistry is included. Even in this elementary example, where the exact ground-state energy and site occupations can be found analytically, there is much to be explained in terms of the underlying logic and aims of density functional theory. Although the usual solution is analytic, the density functional is given only implicitly. We overcome this difficulty using the Levy-Lieb construction to create a parametrization of the exact function with negligible errors. The symmetric case is most commonly studied, but we find a rich variation in behavior by including asymmetry, as strong correlation physics vies with charge-transfer effects. We explore the behavior of the gap and the many-body Green's function, demonstrating the 'failure' of the Kohn-Sham (KS) method to reproduce the fundamental gap. We perform benchmark calculations of the occupation and components of the KS potentials, the correlation kinetic energies, and the adiabatic connection. We test several approximate functionals (restricted and unrestricted Hartree-Fock and Bethe ansatz local density approximation) to show their successes and limitations. We also discuss and illustrate the concept of the derivative discontinuity. Useful appendices include analytic expressions for density functional energy components, several limits of the exact functional (weak- and strong-coupling, symmetric and asymmetric), various adiabatic connection results, proofs of exact conditions for this model, and the origin of the Hubbard model from a minimal basis model for stretched H2.

arXiv:1502.02194Journal of Physics Condensed MatterSep 17, 2015104 citations

The Hubbard dimmer: a case study crossing the bridge between DFT and strongly correlated systems

Carrascal Camino, J. M. Diego

The use of Density Functional Theory (DFT) has become a standard procedure in the recent years for the simulation of nanostructures. However it is well known that DFT suffers from severe limitations when covering the description of strongly correlated systems for which those approximations are not good.\n\nStrongly correlated systems can be described by means of tight-binding models. Such models are much simpler than the more realistic real-space models as they depend on few parameters and in some cases can be solved algebraically.\n\nThe first three chapters of this thesis are devoted to the use of the DFT machinery to the description of such lattice models in order to assess the true limitations when DFT is applied to strongly correlated systems.\n\nThe fourth and fifth chapters are oriented to the direct application of DFT to the study of two different and interesting problems: the search of molecular magnets with large magnetic anisotropy energy (MAE) and the study of molecular junctions using graphene leads.\n\nIn the first chapter we derive the exact functional and the exchange-correlation energies that let us building the Kohn-Sham (KS) potentials of the single and double site Anderson-Hubbard models in order to find out if the exact DFT gives the correct description of such systems. We find that this description is in general incomplete.\n\nThe second chapter presents an analysis of the double-site asymmetric Hubbard model with 2 electrons. The properties of its ground and excited states are discussed. We then use the Mean-Field (MF) and the Bethe-ansatz LDA (BALDA) approximations to assess their performance as a function of the external potential-interactions ratio. We then apply the concepts of DFT to this model, establishing its KS description and the corresponding adiabatic connection. A parametrized functional is also devised using a method that allows improving the accuracy of any approximated functional whenever the relation between the kinetic and the potential correlation terms is known.\n\nThe third chapter covers the DFT estimation of the fundamental gap in strongly correlated systems. We study the asymmetric Hubbard dimer. We find the exact quasi-particle excitation spectrum, including the ionization potential and electron affinity, which let us derive the real gap. We set a method to calculate the KS potentials for non-integer particle numbers as a necessary step to find the KS Green's function. We prove that the Koopmans' and the Janak's theorems are obeyed. We find that the KS spectrum is in general incorrect, except for the HOMO by definition, when strong correlations are dominant. The KS gap is always smaller than the true gap with interactions and we find the role of the derivative discontinuity on this.\n\nIn the fourth chapter we focus on the practical use of DFT and we test several organometallic molecules based with the objective of finding stable molecular magnets at room temperatures. The current technology uses thin-sheets for manufacturing magnetic bits for data storage. Different alternatives such as the use of single molecular magnets or clusters of transition metals to reduce the physical size of a bit. Heavy metal dimers may show huge magnetic anisotropies provided by a strong spin-orbit coupling (SOC). We analyze the magnetic properties of several organometallic molecules based on transition metals to determine if such huge anisotropies survive when the metal is encapsulated into the organic molecule.\n\nThe fifth chapter analyzes the use of graphene in the development of nanoelectronic circuits based on molecular junctions. The graphene has unique properties excellently matching the requirements for ideal electrodes. However the use of graphene nanoelectronics is difficulty because the shape of the graphene edges can change arbitrarily the coupling and switching functionality of the molecules attached. The aim of this chapter is to investigate the electronics and transport properties of different edge morphologies of graphene-wedge single-molecule junctions. It is shown that some graphene-wedge BDT and BPD junctions show negative differential resistance (NDR), spin-filtering (SP) behavior and Fano resonances.

Jul 17, 20152 citations

The Hubbard Dimer: A density functional case study of a many-body problem

D. J. Carrascal, J. Ferrer, J. C. Smith +1 more

This review explains the relationship between density functional theory and strongly correlated models using the simplest possible example, the two-site Hubbard model. The relationship to traditional quantum chemistry is included. Even in this elementary example, where the exact ground-state energy and site occupations can be found analytically, there is much to be explained in terms of the underlying logic and aims of Density Functional Theory. Although the usual solution is analytic, the density functional is given only implicitly. We overcome this difficulty using the Levy-Lieb construction to create a parametrization of the exact function with negligible errors. The symmetric case is most commonly studied, but we find a rich variation in behavior by including asymmetry, as strong correlation physics vies with charge-transfer effects. We explore the behavior of the gap and the many-body Green’s function, demonstrating the ‘failure ’ of the Kohn-Sham method to reproduce the fundamental gap. We perform benchmark calculations of the occupation and components of the KS potentials, the correlation kinetic energies, and the adiabatic connection. We test several approximate functionals (restricted and unrestricted Hartree-Fock and Bethe Ansatz Local Density Approximation) to show their successes and limitations. We also discuss and illustrate the concept of the derivative discontinuity. Useful appendices include analytic expressions for Density

Jan 1, 2015106 citations

Excited-state-specific Kohn–Sham formalism for the asymmetric Hubbard dimer

Pierre‐François Loos, Sara Giarrusso

Building on our recent study [Giarrusso and Loos, J. Phys. Chem. Lett. 14, 8780 (2023)], we explore the generalization of the ground-state Kohn-Sham (KS) formalism of density-functional theory (DFT) to the (singlet) excited states of the asymmetric Hubbard dimer at half-filling. While we found that the KS-DFT framework can be straightforwardly generalized to the highest-lying doubly excited state, the treatment of the first excited state presents significant challenges. In particular, using a density-fixed adiabatic connection, we show that the density of the first excited state lacks non-interacting v-representability. However, by employing an analytic continuation of the adiabatic path, we demonstrate that the density of the first excited state can be generated by a complex-valued external potential in the non-interacting case. More practically, by performing state-specific KS calculations with exact and approximate correlation functionals-each state possessing a distinct correlation functional-we observe that spurious stationary solutions of the KS equations may arise due to the approximate nature of the functional.

arXiv:2412.14945The Journal of Chemical PhysicsApr 8, 20256 citations

Exact maps in density-functional theory

Tanja Dimitrov

Providing sufficient accuracy at low-computational cost, not surprisingly, density-functional theory ( DFT ) is a standard tool in theoretical material science. Similarly, quantum-electrodynamical density-functional theory ( QEDFT ) provides a promising route to describe experiments at the interface of quantum optics and quantum chemistry. Formally, in the Kohn-Sham ( KS ) formulation of DFT and QEDFT the exact many-particle Schrödinger equation is replaced by a set of coupled non-linear one-particle differential equations. The electron-electron, the electron-ion and the electron-photon interactions of the interacting system are mimicked by their corresponding effective KS potential with an, in principle unknown exchange-correlation ( xc ) term. In practice, DFT and QEDFT require approximations for this xc potential. The quality of which fully determines how much the DFT and QEDFT observables deviate from the exact ones. Understanding the shortcomings and improving the quality of approximate functionals makes the study of the exact functional of model systems indispensable. Here, we employ exact diagonalization for the Hamiltonian of different types of real-space lattice models in one spatial dimension to access the ground-state wave function, and thus to explicitly construct the exact functionals of DFT and QEDFT . To study the functional behavior of purely electronic systems, we consider two singlet-state electrons moving in an external potential with two connected atomic wells of variable depth and distance. To study coupled electron-photon systems, we consider the Rabi Hamiltonian with variable electron-photon coupling strength. By varying the model parameters, we systematically construct the functionals for systems from weak to strong coupling between the building blocks of the global system. Using the concept of the effective density-to- potential map, we discuss the intrasystem steepening which is an exact condition on the xc functional that captures the static correlation present in the system and carve out how it relates to the well-known intersystem derivative discontinuity. We further construct arbitrary observables as functionals of the internal variables which are directly affected by such features of the effective density-to-potential map. By constructing the reduced dipole map in the four-particle case of coupled electron-photon system, we show that the density-to-potential map helps to identify and understand recently discovered physical effects in the field of strongly coupled electron-photon systems. This work provides a contribution to the development of approximate function- als required in practical density-functional calculations of arbitrary experimentally accessible observables. Both the intrasystem steepening and the intersystem derivative discontinuity play a crucial role for systems with highly-localized electrons. Including both features in approximate functionals will improve the accu- racy of DFT evaluation of observables for such systems. Those results provide novel insight into the properties of the exact xc potentials in QEDFT . This work opens a new route to construct new functionals using coarse-grained or embed- ding approaches within the general and exact frameworks of DFT and QEDFT.

Jan 1, 20171 citation

Exact thermal density functional theory for a model system: Correlation components and accuracy of the zero-temperature exchange-correlation approximation

Justin C. Smith, Aurora Pribram−Jones, Kieron Burke

Thermal density functional theory calculations often use the Mermin-Kohn-Sham scheme, but employ ground-state approximations to the exchange-correlation (XC) free energy. In the simplest solvable nontrivial model, an asymmetric Hubbard dimer, we calculate the exact many-body energies and the exact Mermin-Kohn-Sham functionals for this system and extract the exact XC free energy. For moderate temperatures and weak correlation, we find this approximation to be excellent. We extract various exact free-energy correlation components and the exact adiabatic connection formula.

arXiv:1509.03097Physical review. B./Physical review. BJun 14, 201642 citations