The exact conditions for density functionals and density matrix functionals in terms of fractional charges and fractional spins are known, and their violation in commonly used functionals has been shown to be the root of many major failures in practical applications. However, approximate functionals are designed for physical systems with integer charges and spins, not in terms of the fractional variables. Here we develop a general framework for extending approximate density functionals and many-electron theory to fractional-charge and fractional-spin systems. Our development allows for the fractional extension of any approximate theory that is a functional of G(0), the one-electron Green's function of the non-interacting reference system. The extension to fractional charge and fractional spin systems is based on the ensemble average of the basic variable, G(0). We demonstrate the fractional extension for the following theories: (1) any explicit functional of the one-electron density, such as the local density approximation and generalized gradient approximations; (2) any explicit functional of the one-electron density matrix of the non-interacting reference system, such as the exact exchange functional (or Hartree-Fock theory) and hybrid functionals; (3) many-body perturbation theory; and (4) random-phase approximations. A general rule for such an extension has also been derived through scaling the orbitals and should be useful for functionals where the link to the Green's function is not obvious. The development thus enables the examination of approximate theories against known exact conditions on the fractional variables and the analysis of their failures in chemical and physical applications in terms of violations of exact conditions of the energy functionals. The present work should facilitate the calculation of chemical potentials and fundamental bandgaps with approximate functionals and many-electron theories through the energy derivatives with respect to the fractional charge. It should play an important role in developing accurate approximate density functionals and many-body theory.
The exact conditions for density functionals and density matrix functionals in terms of fractional charges and fractional spins are known, and their violation in commonly used functionals has been shown to be the root of many major failures in practical applications. However, approximate functionals are designed for physical systems with integer charges and spins, not in terms of the fractional variables. Here we develop a general framework for extending approximate density functionals and many-electron theory to fractional-charge and fractional-spin systems. Our development allows for the fractional extension of any approximate theory that is a functional of G(0), the one-electron Green's function of the non-interacting reference system. The extension to fractional charge and fractional spin systems is based on the ensemble average of the basic variable, G(0). We demonstrate the fractional extension for the following theories: (1) any explicit functional of the one-electron density, such as the local density approximation and generalized gradient approximations; (2) any explicit functional of the one-electron density matrix of the non-interacting reference system, such as the exact exchange functional (or Hartree-Fock theory) and hybrid functionals; (3) many-body perturbation theory; and (4) random-phase approximations. A general rule for such an extension has also been derived through scaling the orbitals and should be useful for functionals where the link to the Green's function is not obvious. The development thus enables the examination of approximate theories against known exact conditions on the fractional variables and the analysis of their failures in chemical and physical applications in terms of violations of exact conditions of the energy functionals. The present work should facilitate the calculation of chemical potentials and fundamental bandgaps with approximate functionals and many-electron theories through the energy derivatives with respect to the fractional charge. It should play an important role in developing accurate approximate density functionals and many-body theory.
Standard approximations for the exchange-correlation functional have been found to give big errors for the linearity condition of fractional charges, leading to delocalization error, and the constancy condition of fractional spins, leading to static correlation error. These two conditions are now unified and extended to states with both fractional charge and fractional spin to give a much more stringent condition: the exact energy functional is a plane, linear along the fractional charge coordinate and constant along the fractional spin coordinate with a line of discontinuity at the integer. Violation of this condition underlies the failure of all known approximate functionals to describe the gaps in strongly correlated systems. It is shown that explicitly discontinuous functionals of the density or orbitals that go beyond these currently used smooth approximations is the key for the application of density functional theory to strongly correlated systems.
A framework for developing new approximate electronic structure methods is presented, in which the correlation energy of a many-electron system in the ground state is computed as in the single-reference second-order many-body perturbation theory but with the reference one-body Hamiltonian modified by additional (generally non-local) potentials which are universal functionals of the Hartree-Fock density. The existence of such functionals that reproduce the exact correlation energy justifies the search for approximate models which may overcome some important deficiencies of the traditional density-functional methods and also perform better than the usual low-order wavefunction methods.
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.
The study of systems with fractional charges and spins has become an extremely important tool to understand errors in approximate electronic structure methods, particularly in the context of density functional theory. Meanwhile, similar studies with wavefunction (WF)-based methods beyond second-order perturbation theory have been lacking. In this contribution, we study the performance of different coupled cluster (CC) and many-body perturbation theory (MBPT)-based methods for fractional charges. The use of the conventional and renormalized formulations of fractional-charge MBPT is discussed. The fractional spin behavior of the coupled cluster doubles (CCD) method is also investigated. Overall, all tested WF methods show very promising performance for the fractional charge problem. CCD is also quite accurate for the fractional spin problem in He+ across most of the range, although it breaks down to near Hartree-Fock quality in the strongly correlated limit. Beyond the study of fractional charge and spin curves, the implementation of CC methods with fractional occupation numbers offers a promising route to treating problems with multi-reference character in a single-reference framework.
The description of many-electron systems with a fractional electron number, Ntot, and fractional z-projection of the spin, Mtot, is of great importance in physical chemistry, solid-state physics, and materials science. In this study, we analyze the fundamental question of what the ensemble ground state of a general, finite, many-electron system at zero temperature is, with a given Ntot and Mtot, distinguishing between low- and high-spin cases (separated by the boundary spin MB). For the low-spin case, the general form of the ensemble ground state has been rigorously derived in Goshen and Kraisler [J. Phys. Chem. Lett. 15, 2337 (2024)], generalizing the piecewise linearity and the flat-plane conditions for many-electron systems. Here, we provide an alternative proof for this case, discuss the ambiguity in the description of the ground state, and show that this ambiguity can be removed via maximization of the system's entropy. For the high-spin case, we find that the form of the ensemble ground state strongly depends on the system in question. We prove three general properties that characterize the ground state at high spins and narrow down the list of pure states it may consist of. We illustrate the aforementioned properties of high-spin cases by examining the ensemble ground state when Mtot approaches MB from above during the addition of (a fraction of) an up- and down-electron to a given system. Furthermore, we relate the frontier orbital energies of Kohn-Sham (KS) density functional theory (DFT) to total energy differences at high spin values, particularly the ionization potential (IP), the fundamental gap, and the spin flip energies. Analyzing the frontier energies on both sides of each boundary in the total energy profile, where the energy slope changes abruptly, we derive expressions for new derivative discontinuities, which are predicted to appear as jumps in the corresponding KS potentials. In this way, we generalize the well-known IP theorem of DFT to cases with fractional electron number and to cases with high spin. Our analytical results are supported by an extensive numerical analysis of the Atomic Spectra Database of the National Institute of Standards. The new exact conditions for many-electron systems derived in this study are instrumental for the development of advanced approximations in DFT and other many-electron methods.
Density-Functional Theory (DFT), although widely used and very successful in the calculation of several observables, fails to correctly describe strongly correlated materials. In the first part of this work we, therefore, introduce reduced-density-matrix-functional theory (RDMFT) which is one possible way to treat electron correlation beyond DFT. Within this theory the one-body reduced density matrix (1-RDM) is used as the basic variable. Our main interest is the calculation of the fundamental gap which proves very problematic within DFT. In order to calculate the fundamental gap we generalize RDMFT to fractional particle numbers M by describing the system as an ensemble of an N and an N+1 particle system (with N≤ M≤N+1). For each fixed particle number, M, the total energy is minimized with respect to the natural orbitals and their occupation numbers. This leads to the total energy as a function of M. The derivative of this function with respect to the particle number has a discontinuity at integer particle number which is identical to the gap. In addition, we investigate the necessary and sufficient conditions for the 1-RDM of a system with fractional particle number to be N-representable. Numerical results are presented for alkali atoms, small molecules, and periodic systems. Another problem within DFT is the description of non-relativistic many- electron systems in the presence of magnetic fields. It requires the paramagnetic current density and the spin magnetization to be used as basic variables besides the electron density. However, electron-gas-based functionals of current-spin-density-functional Theory (CSDFT) exhibit derivative discontinuities as a function of the magnetic field whenever a new Landau level is occupied, which makes them difficult to use in practice. Since the appearance of Landau levels is, intrinsically, an orbital effect it is appealing to use orbital-dependent functionals. We have developed a CSDFT version of the optimized effective potential (OEP) method which allows for the use of explicitly orbital-dependent functionals in the context of CSDFT. We present the derivation of corresponding equations and show results for a quantum dot in external magnetic fields in the second part of this thesis.