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 developers and practitioners.

Benchmarking studies in computational chemistry use reference datasets to assess the accuracy of a method through error statistics. The commonly used error statistics, such as the mean signed and mean unsigned errors, do not inform end-users on the expected amplitude of prediction errors attached to these methods. We show that, the distributions of model errors being neither normal nor zero-centered, these error statistics cannot be used to infer prediction error probabilities. To overcome this limitation, we advocate for the use of more informative statistics, based on the empirical cumulative distribution function of unsigned errors, namely, (1) the probability for a new calculation to have an absolute error below a chosen threshold and (2) the maximal amplitude of errors one can expect with a chosen high confidence level. Those statistics are also shown to be well suited for benchmarking and ranking studies. Moreover, the standard error on all benchmarking statistics depends on the size of the reference dataset. Systematic publication of these standard errors would be very helpful to assess the statistical reliability of benchmarking conclusions.

It was found by Perdew, Parr, Levy, and Balduz [Phys. Rev. Lett. {\bf 49}, 1691 (1982)] and by Sham and Schl\"uter [Phys. Rev. Lett. {\bf 51}, 1884 (1983)] that the exact Kohn-Sham exchange-correlation potential of an open system may jump discontinuosly as the particle number crosses an integer, with important physical consequences. Recently, Sagvolden and Perdew [Phys. Rev. A {\bf 77}, 012517 (2008)] have analyzed the discontinuity of the exchange-correlation potential as the particle number crosses one, with an illustration that uses a model density for the H$^-$ ion. In this work, we extend their analysis to the case in which the external potential is the simple harmonic confinement, choosing spring-constant values for which the two-electron hamiltonian has an analytic solution. This way, we can obtain the exact, analytic, exchange and correlation potentials for particle number fluctuating between zero and two, illustrating the discontinuity as the particle number crosses one without introducing any model or approximation. We also discuss exchange and correlation separately.Comment: Submitted to Int. J. Quantum Chem., special issue honoring Prof. Mayer. New version, where an important error has been correcte