Single-particle energies and density of states in density functional theory

Time-dependent density functional theory (TD-DFT) is commonly used as the foundation to obtain neutral excited states and transition weights in DFT, but does not allow direct access to density of states and single-particle energies, i.e. ionisation energies and electron affinities. Here we show that by extending TD-DFT to a superfluid formulation, which involves operators that break particle-number symmetry, we can obtain the density of states and single-particle energies from the poles of an appropriate superfluid response function. The standard Kohn– Sham eigenvalues emerge as the adiabatic limit of the superfluid response under the assumption that the exchange– correlation functional has no dependence on the superfluid density. The Kohn– Sham eigenvalues can thus be interpreted as approximations to the ionisation energies and electron affinities. Beyond this approximation, the formalism provides an incentive for creating a new class of density functionals specifically targeted at accurate single-particle eigenvalues and bandgaps.     

An efficient implementation of two-component relativistic density functional theory with torque-free auxiliary variables

An integration and assembly strategy for efficient evaluation of the exchange correlation term in relativistic density functional theory within two-component Kohn–Sham framework is presented. Working equations that both take into account all the components of the spin magnetization and can exploit parallelism, optimized cache utilization, and micro-architecture specific-floating point operations are discussed in detail in this work. The presented assembly of the exchange correlation potential, suitable for both open and closed shell systems, uses spinor density and a set of auxiliary variables, ensuring easy retrofitting of existing density functionals designed for collinear density. The used auxiliary variables in this paper, based on the scalar and non-collinear density, can preserve non-zero exchange correlation magnetic field local torque, without violating the required overall zero torque, even for GGA functionals. This is mandatory to obtain accurate spin dynamics and proper time evolution of the magnetization. Spin frustrated hydrogen rings are used to validate the current implementation and phenoxy radicals of different sizes are used to monitor the performance. This approach is a step towards extending the applicability of relativistic two-component DFT to systems of large size (>100 atoms).     

Wavefunction stability analysis without analytical electronic Hessians: application to orbital-optimised second-order Møller–Plesset theory and VV10-containing density functionals

Wavefunction stability analysis is commonly applied to converged self-consistent field (SCF) solutions to verify whether the electronic energy is a local minimum with respect to second-order variations in the orbitals. By iterative diagonalisation, the procedure calculates the lowest eigenvalue of the stability matrix or electronic Hessian. However, analytical expressions for the electronic Hessian are unavailable for most advanced post-Hartree–Fock (HF) wave function methods and even some Kohn–Sham (KS) density functionals. To address such cases, we formulate the Hessian-vector product within the iterative diagonalisation procedure as a finite difference of the electronic gradient with respect to orbital perturbations in the direction of the vector. As a model application, following the lowest eigenvalue of the orbital-optimised second-order Møller–Plesset perturbation theory (OOMP2) Hessian during H₂ dissociation reveals the surprising stability of the spin-restricted solution at all separations, with a second independent unrestricted solution. We show that a single stable solution can be recovered by using the regularised OOMP2 method (δ-OOMP2), which contains a level shift. Internal and external stability analyses are also performed for SCF solutions of a recently developed range-separated hybrid density functional, ωB97X-V, for which the analytical Hessian is not yet available due to the complexity of its long-range non-local VV10 correlation functional.     

Random-phase approximation correlation methods for molecules and solids

Random-phase approximation (RPA) correlation methods based on Kohn–Sham density-functional theory and Hartree–Fock are derived using the adiabatic-connection fluctuation dissipation theorem. It is shown that the correlation energy within the adiabatic-connection fluctuation-dissipation theorem is exact in a Kohn–Sham framework while for Hartree–Fock reference states this is not the case. This shows that Kohn–Sham reference states are probably better suited to describe electron correlation for use in RPA methods than Hartree–Fock reference states. Both, Kohn–Sham and Hartree–Fock RPA methods are related to each other both by comparing the underlying correlation functionals and numerically through the comparison of total energies and reaction energies for a set of small organic molecules.     

Time-dependent generalized Kohn–Sham theory

Generalized Kohn–Sham (GKS) theory extends the realm of density functional theory (DFT) by providing a rigorous basis for non-multiplicative potentials, the use of which is outside original Kohn–Sham theory. GKS theory is of increasing importance as it underlies commonly used approximations, notably (conventional or range-separated) hybrid functionals and meta-generalized-gradient-approximation (meta-GGA) functionals. While this approach is often extended in practice to time-dependent DFT (TDDFT), the theoretical foundation for this extension has been lacking, because the Runge–Gross theorem and the van Leeuwen theorem that serve as the basis of TDDFT have not been generalized to non-multiplicative potentials. Here, we provide the necessary generalization. Specifically, we show that with one simple but non-trivial additional caveat – upholding the continuity equation in the GKS electron gas – the Runge–Gross and van Leeuwen theorems apply to time-dependent GKS theory. We also discuss how this is manifested in common GKS-based approximations.     

Thirty years of density functional theory in computational chemistry: an overview and extensive assessment of 200 density functionals

In the past 30 years, Kohn–Sham density functional theory has emerged as the most popular electronic structure method in computational chemistry. To assess the ever-increasing number of approximate exchange-correlation functionals, this review benchmarks a total of 200 density functionals on a molecular database (MGCDB84) of nearly 5000 data points. The database employed, provided as Supplemental Data, is comprised of 84 data-sets and contains non-covalent interactions, isomerisation energies, thermochemistry, and barrier heights. In addition, the evolution of non-empirical and semi-empirical density functional design is reviewed, and guidelines are provided for the proper and effective use of density functionals. The most promising functional considered is ωB97M-V, a range-separated hybrid meta-GGA with VV10 nonlocal correlation, designed using a combinatorial approach. From the local GGAs, B97-D3, revPBE-D3, and BLYP-D3 are recommended, while from the local meta-GGAs, B97M-rV is the leading choice, followed by MS1-D3 and M06-L-D3. The best hybrid GGAs are ωB97X-V, ωB97X-D3, and ωB97X-D, while useful hybrid meta-GGAs (besides ωB97M-V) include ωM05-D, M06-2X-D3, and MN15. Ultimately, today's state-of-the-art functionals are close to achieving the level of accuracy desired for a broad range of chemical applications, and the principal remaining limitations are associated with systems that exhibit significant self-interaction/delocalisation errors and/or strong correlation effects.     

Density functional. Theory and application to atoms and molecules

The density functional theory is one of the most efficient and promising methods of quantum physics and chemistry. It is a theory of electronic structure formulated in terms of the electron density as the basic unknown function instead of the electron wave function. According to the fundamental theorems of Hohenberg and Kohn, the electron density carries all the information one might need to determine any property of the electron system. However, the way of obtaining it, is not at all trivial. In this report, the recent advances are summarized. After a review of the Hohenberg–Kohn theorems, the method of constrained search and the Kohn–Sham scheme, exact theorems, relations and inequalities are discussed. There are several important concepts of chemistry (e.g. electronegativity, hardness, softness) that have recently obtained a firm foundation in the density functional theory. The optimized potential method and the methods that generate the potential from the electron density are reviewed. The local and nonlocal approximate functionals are compared. Extensions of the ground-state density functional theory (excited states, time-dependent, relativistic and finite temperature) are summarized. A review of the applications to atoms and molecules is presented.     

Confinement effects on the spin potential of first row transition metal cations

Abstract

By using a spherical confinement method, the behaviour of spin potential and pairing energy is studied and compared to the free ion limit for a representative sample of first row transition metal cations. The study was carried out using three approximations within the Kohn–Sham model; exchange-only, exchange plus correlation contribution and correcting the self-interaction error. For the three approaches, the spin potential shows a close connection with the capability of a system to perform a spin-flip process. Namely, in accordance with Hund’s rule, the spin potential increases from low d occupation up to maximum for the half filled configurations; and it decreases from that point on, as d occupation grows. Such a conclusion is reached for confined and non-confined cations, even under extreme confinement conditions. In addition, two important observations are obtained: (a) In contrast to the neutral atoms situation, in the case of cations no eigenvalue crossings are observed under confinement conditions for the whole sample of ions tested. (b) The self-interaction error found in many exchange–correlation functionals does not affect the pairing energy over confined atoms, even when this error has an important contribution on a single eigenvalue. Therefore, pairing energy predicted by exchange–correlation functionals non-corrected by the self-interaction error can be made safely on transition metal cations under high pressures.     

Atomic C 6 dispersion coefficients: a four-component relativistic Kohn–Sham study

We present C ₆ homo- and heteroatomic dispersion coefficients for all closed-shell atoms of the periodic table based on dipole–dipole polarizabilities at imaginary frequencies calculated using our recent extension of the complex polarization propagator approach to the four-component relativistic Kohn–Sham approach. Lack of proper reference data bars definite conclusions as to which density functional shows the overall best performance, and we therefore call for state-of-the-art wave function-based correlated calculations of dispersion coefficients. Scalar relativistic effects are significant already for elements as light as zinc, whereas spin–orbit effects must be taken into account only for very heavy elements.     

Reformulation of density functional theory for N-representable densities and the resolution of the v-representability problem

Density functional theory for the case of general, N-representable densities is reformulated in terms of density functional derivatives of expectation values of operators evaluated with wave functions leading to a density, making no reference to the concept of potential. The developments provide proof of existence of a mathematical procedure that determines whether a density is v-representable and in the case of an affirmative answer determines the potential (within an additive constant) as a derivative with respect to the density of a constrained search functional. It also establishes the existence of an energy functional of the density that, for v-representable densities, assumes its minimum value at the density describing the ground state of an interacting many-particle system. The theorems of Hohenberg and Kohn emerge as special cases of the formalism. Numerical results for one-dimensional non-interacting systems illustrate the formalism. Some direct formal and practical implications of the present reformulation of DFT are also discussed.     

Sum-rules of the response potential in the strongly-interacting limit of DFT

The response part of the exchange-correlation potential of Kohn–Sham density functional theory plays a very important role, for example for the calculation of accurate band gaps and excitation energies. Here we analyze this part of the potential in the limit of infinite interaction in density functional theory, showing that in the one-dimensional case it satisfies a very simple sum rule.     

Density scaling and virial theorem

The virial theorem, the Levy–Perdew relation and the differential virial theorem are derived for density–scaled Kohn–Sham systems. Earlier it was shown that there exists a value of the scaling factor for which the correlation energy disappears and we should treat only exchange for which a simple approximation was proposed. The new Levy–Perdew relation is applied to judge the quality of this approximation.     

On the exact formulation of multi-configuration density-functional theory: electron density versus orbitals occupation

The exact formulation of multi-configuration density-functional theory is discussed in this work. As an alternative to range-separated methods, where electron correlation effects are split in the coordinate space, the combination of configuration interaction methods with orbital occupation functionals is explored at the formal level through the separation of correlation effects in the orbital space. When applied to model Hamiltonians, this approach leads to an exact site-occupation embedding theory (SOET). An adiabatic connection expression is derived for the complementary bath functional and a comparison with density matrix embedding theory is made. Illustrative results are given for the simple two-site Hubbard model. SOET is then applied to a quantum chemical Hamiltonian, thus leading to an exact complete active space site-occupation functional theory (CASSOFT) where active electrons are correlated explicitly within the CAS and the remaining contributions to the correlation energy are described with an orbital occupation functional. The computational implementation of SOET and CASSOFT as well as the development of approximate functionals are left for future work.     

Adaptive local basis set for Kohn–Sham density functional theory in a discontinuous Galerkin framework I: Total energy calculation

Kohn–Sham density functional theory is one of the most widely used electronic structure theories. In the pseudopotential framework, uniform discretization of the Kohn–Sham Hamiltonian generally results in a large number of basis functions per atom in order to resolve the rapid oscillations of the Kohn–Sham orbitals around the nuclei. Previous attempts to reduce the number of basis functions per atom include the usage of atomic orbitals and similar objects, but the atomic orbitals generally require fine tuning in order to reach high accuracy. We present a novel discretization scheme that adaptively and systematically builds the rapid oscillations of the Kohn–Sham orbitals around the nuclei as well as environmental effects into the basis functions. The resulting basis functions are localized in the real space, and are discontinuous in the global domain. The continuous Kohn–Sham orbitals and the electron density are evaluated from the discontinuous basis functions using the discontinuous Galerkin (DG) framework. Our method is implemented in parallel and the current implementation is able to handle systems with at least thousands of atoms. Numerical examples indicate that our method can reach very high accuracy (less than 1 meV) with a very small number (4–40) of basis functions per atom.     

Partition density functional theory and its extension to the spin-polarized case

Partition density functional theory (PDFT) [P. Elliott, K. Burke, M.H. Cohen, and A. Wasserman, Phys. Rev. A 82 (2), 024501 (2010)] is a formally exact method for obtaining molecular properties from Kohn–Sham calculations on isolated fragments. Here, we express the partition energy of PDFT as an implicit functional of the molecular spin-densities for a given choice of fragmentation, and use the principle of von Barth and Hedin to formulate the spin-decomposed version of PDFT. We introduce a partition energy functional of the spin-up and spin-down electronic densities and derive the associated polarized partition potentials, which are found to be global quantities that influence every fragment in the molecule. Along with the formal theory, we propose a simplified approach to computing the spin-partition potentials, and illustrate its utility and accuracy with two simple examples. Finally, we propose a viable approach to including external electric and magnetic fields in the framework of spin-PDFT.     

Multiple exciton generation in silicon quantum dot arrays: density functional perturbation theory computation

We use the density functional theory (DFT) combined with the many-body perturbation theory to derive expressions for the rates of the optical photon→exciton and photon→bi-exciton processes in nanoparticles, and for quantum efficiency, all to the leading order in the screened Coulomb interaction between Kohn–Sham quasiparticles. Also, we calculate exciton→bi-exciton rates due to the impact ionisation (II) mechanism in Si₂₉H₃₆ quantum dots (QDs) with both crystalline and amorphous core structures, and in quasi-one dimensional (1-D) arrays constructed from these QDs. We observe significant dependence of the carrier multiplication rates on the structure’s morphology and structural disorder. Amorphous silicon QD arrays are predicted to have more efficient bi-exciton generation rates as a function of exciton energy compared to their crystalline counterparts, and the isolated QDs of both kinds.