In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method...

In this article, we examine Hohenberg–Kohn theorems for Current Density Functional Theory, that is, generalizations of the classical Hohenberg–Kohn theorem that includes both electric and magnetic fields. In the Vignale and Rasolt formulation (Vignale and Rasolt, Phys. Rev. Lett. 1987, 59, 2360), which uses the paramagnetic current density, we address the issue of degenerate ground states and prove that the ensemble‐representable particle and paramagnetic current density determine the degenerate ground states. For the formulation that uses the total current density, we note that the proof suggested by Diener (Diener, J. Phys.: Condens. Matter. 1991, 3, 9417) is unfortunately not correct. Furthermore, we give a proof that the magnetic field and the ensemble‐representable particle density determine the scalar and vector potentials up to a gauge transformation. This generalizes the result of Grayce and Harris (Grayce and Harris, Phys. Rev. A 1994, 50, 3089) to the case of degenerate ground states. We moreover prove the existence of a positive wavefunction that is the ground state of infinitely many different Hamiltonians. © 2014 Wiley Periodicals, Inc.

The mathematical foundation of the so-called extended coupled-cluster method for the solution of the many-fermion Schrödinger equation is here developed. We prove an existence and uniqueness result, both in the full infinite-dimensional amplitude space as well as for discretized versions of it. The extended coupled-cluster method is formulated as a critical point of an energy function using a generalization of the Rayleigh-Ritz principle: the bivariational principle. This gives a quadratic bound for the energy error in the discretized case. The existence and uniqueness results are proved using a type of monotonicity property for the flipped gradient of the energy function.Comparisons to the analysis of the standard coupled-cluster method is made, and it is argued that the bivariational principle is a useful tool, both for studying coupled-cluster type methods, and for developing new computational schemes in general.

A detailed account of the Kohn-Sham algorithm from quantum chemistry, formulated rigorously in the very general setting of convex analysis on Banach spaces, is given here. Starting from a Levy-Lieb-type functional, its convex and lower semi-continuous extension is regularized to obtain differentiability. This extra layer allows to rigorously introduce, in contrast to the common unregularized approach, a well-defined Kohn-Sham iteration scheme. Convergence in a weak sense is then proven. This generalized formulation is applicable to a wide range of different density-functional theories and possibly even to models outside of quantum mechanics.

In this article, we investigate the numerical and theoretical aspects of the coupled-cluster method tailored by matrix-product states. We investigate formal properties of the used method, such as energy size consistency and the equivalence of linked and unlinked formulation. The existing mathematical analysis is here elaborated in a quantum chemical framework. In particular, we highlight the use of what we have defined as a complete active space-external space gap describing the basis splitting between the complete active space and the external part generalizing the concept of a HOMO−LUMO gap. Furthermore, the behavior of the energy error for an optimal basis splitting, i.e., an active space choice minimizing the density matrix renormalization group-tailored coupled-cluster singles doubles error, is discussed. We show numerical investigations on the robustness with respect to the bond dimensions of the single orbital entropy and the mutual information, which are quantities that are used to choose a complete active space. Moreover, the dependence of the groundstate energy error on the complete active space has been analyzed numerically in order to find an optimal split between the complete active space and external space by minimizing the density matrix renormalization group-tailored coupled-cluster error.