In this dissertation, I introduce both ground-state and thermal density functional theory. Throughout I use the asymmetric two-site Hubbard model, called the Hubbard dimer for short, to better understand and/or develop these theories. This model is used because it can be solved analytically and it contains all the necessary physics while still being conceptually simple enough to tease apart the various aspects of density functional theory. Ground-statedensity functional theory has seen broad use in many disciplines including physics, chemistry, geology, and material science and has led to a number of important physical and technological successes. In the first two chapters I elucidate the behavior of the ground-state theory using the Hubbard dimer. The simplicity of the model allows me to showcase aspects of the theory that are common points of confusion within the electronic structure community, e.g. the fundamental gap problem. The next two chapters focus on thermal density functional theory which has been coming to prominence as the study of warm dense matter has becomea growing interest at the national laboratories and in the astronomical body community. The Hubbard dimer allows me to do the first ever exact thermal density functional theory calculation. In this work I am better able to understand the approximations used in thermal density functional theory and can point to why they succeed and fail. This also allows me to illustrate old conditions and derive new ones. I conclude with an overview of the work and a few different directions in which the asymmetric Hubbard dimer could be used further.
In this dissertation, I introduce both ground-state and thermal density functional theory. Throughout I use the asymmetric two-site Hubbard model, called the Hubbard dimer for short, to better understand and/or develop these theories. This model is used because it can be solved analytically and it contains all the necessary physics while still being conceptually simple enough to tease apart the various aspects of density functional theory. Ground-statedensity functional theory has seen broad use in many disciplines including physics, chemistry, geology, and material science and has led to a number of important physical and technological successes. In the first two chapters I elucidate the behavior of the ground-state theory using the Hubbard dimer. The simplicity of the model allows me to showcase aspects of the theory that are common points of confusion within the electronic structure community, e.g. the fundamental gap problem. The next two chapters focus on thermal density functional theory which has been coming to prominence as the study of warm dense matter has becomea growing interest at the national laboratories and in the astronomical body community. The Hubbard dimer allows me to do the first ever exact thermal density functional theory calculation. In this work I am better able to understand the approximations used in thermal density functional theory and can point to why they succeed and fail. This also allows me to illustrate old conditions and derive new ones. I conclude with an overview of the work and a few different directions in which the asymmetric Hubbard dimer could be used further.
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.
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
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.
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.
In this Thesis, we develop and apply an accurate method to determine the ground-state properties of strongly-correlated electrons in the framework of lattice model Hamiltonians. In lattice density-functional theory (LDFT) the basic variable is the one-particle density matrix γ. From a generalized HK theorem, the ground state energy Egs[γgs] = minγ E[γ] is obtained by minimizing the energy-functional E[γ] over all the physical or representable γ. The energy functional can be divided into two contributions: the kinetic-energy functional T [γ], which linear dependence on γ is exactly known, and the correlation-energy functional W [γ], which dependence on γ is not explicitly known. Finding accurate approximations for W [γ] constitutes the actual challenge of this thesis. Part of this work is based on previous studies which derive an approximation of W [γ] for the Hubbard Hamiltonian based on scaling hypothesis and an exact analytical results for the Hubbard dimer. However, this approach is limited to spin independent and homogeneous systems. In order to extend the scope of LDFT, we have developed three different approaches to derive W [γ] in order to study broken symmetry problems. First we extend the previous scaling functional to system having charge transfer. By a systematic study of the dependence of W [γ] on the charge distribution we find similar scaling properties as for the homogeneous case. An extension to the Hubbard model on bipartite lattices is then derived and applied to finite and infinite systems with repulsive and attractive interactions. The high accuracy of this functional is demonstrated. Nevertheless, this approach is difficult to transfer to more complex systems, since the calculation of W [γ] applies to the system as a whole. To overcome this problem, we derive another approximation based on local scaling properties. This functional is formulated on-site so that it can be applied to any ordered or disordered Hamiltonian with on-site interaction. As applications we study the metal-insulator transition in the ionic Hubbard model in oneand two-dimensions, as well as the one-dimensional Hubbard chain with first and second nearest neighbor hoppings. Finally, we develop a numerical approach of W [γ], on the basis on exact diagonalizations of an effective many-body Hamiltonian corresponding to a cluster surrounded by an effective medium. This effective Hamiltonian depends on the density matrix γ and allows to derive approximation to W [γ] which improve systematically with increasing cluster size. The formulation is spin-dependent and allows a straightforward generalization to multi-orbital correlated systems such spd-Hamiltonians. In addition it take into account the effects of short-range charge and spin fluctuations in the functional. The accuracy of the method is demonstrated for the Hubbard model by comparison with the Bethe-Ansatz solution (1D) and quantum Monte Carlo simulations (2D). Finally, an outlook or relevant future developments of the theory is provided.
It is common wisdom that a better understanding of how to quantify electron correlation should lead to more efficient treatments of the correlation problem. Here, we not only substantiate this view but turn it into a concrete, systematic framework for ground states. Since prominent correlation measures such as the nonfreeness depend only on the one-particle reduced density matrix, modeling the correlation energy as a function of such a measure naturally places the framework within reduced density matrix functional theory, while keeping the cost at mean-field level. Strikingly, the variable conjugate to the measure emerges as a ``correlation temperature'', thereby connecting correlation measures, functional theory, and effective temperatures. The framework rationalizes methods such as i-DMFT, and we illustrate its key concepts analytically for the Hubbard dimer and numerically for small molecular systems. This brings a fundamental question within reach: how much of the correlation problem can be solved at mean-field cost?