Sparse Approximation of the Lieb Functional with Moment Constraints

Abstract. The aim of this paper is to present new sparsity results about the so-called Lieb functional, which is a key quantity in Density Functional Theory for electronic structure calculations of molecules. The Lieb functional was shown by Lieb to be a convexification of the so-called Lévy–Lieb functional. Given an electronic density for a system of [Formula: see text] electrons, which may be seen as a probability density on [Formula: see text], the value of the Lieb functional for this density is defined as the solution of a quantum multimarginal optimal transport problem, which reads as a minimization problem defined on the set of trace-class operators acting on the space of electronic wave-functions that are antisymmetric [Formula: see text] functions of [Formula: see text], with partial trace equal to the prescribed electronic density. We introduce a relaxation of this quantum optimal transport problem where the full partial trace constraint is replaced by a finite number of moment constraints on the partial trace of the set of operators. We show that, under mild assumptions on the electronic density, there exist sparse minimizers to the resulting moment-constrained approximation of the Lieb (MCAL) functional that read as operators with rank at most equal to the number of moment constraints. We also prove under appropriate assumptions on the set of moment functions that the value of the MCAL functional converges to the value of the exact Lieb functional as the number of moments go to infinity. We also prove some rates of convergence on the associated approximation of the ground state energy. We finally study the mathematical properties of the associated dual problem and introduce a suitable numerical algorithm, named Quantum Column Generation Algorithm (QGCA), in order to solve some simple toy models.