Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2106.02560 · Physical Review Letters

Ensemble Reduced Density Matrix Functional Theory for Excited States and Hierarchical Generalization of Pauli’s Exclusion Principle

Christian Schilling, Stefano Pittalis

Published Jul 7, 2021• 48 citations
arXivHTMLPDF

Abstract

We propose and work out a reduced density matrix functional theory (RDMFT) for calculating energies of eigenstates of interacting many-electron systems beyond the ground state. Various obstacles which historically have doomed such an approach to be unfeasible are overcome. First, we resort to a generalization of the Ritz variational principle to ensemble states with fixed weights. This in combination with the constrained search formalism allows us to establish a universal functional of the one-particle reduced density matrix. Second, we employ tools from convex analysis to circumvent the too involved N-representability constraints. Remarkably, this identifies Valone's pioneering work on RDMFT as a special case of convex relaxation and reveals that crucial information about the excitation structure is contained in the functional's domain. Third, to determine the crucial latter object, a methodology is developed which eventually leads to a generalized exclusion principle. The corresponding linear constraints are calculated for systems of arbitrary size.

DOI
10.1103/physrevlett.127.023001

Related papers

Foundation of One-Particle Reduced Density Matrix Functional Theory for Excited States

Julia Liebert, Federico Castillo, Jean‐Philippe Labbé +1 more

In Phys. Rev. Lett . 2021, 127, 023001 a reduced density matrix functional theory (RDMFT) was proposed for calculating energies of selected eigenstates of interacting many-Fermion systems. Here, we develop a solid foundation for this so-called w -RDMFT and present the details of various derivations. First, we explain how a generalization of the Ritz variational principle to ensemble states with fixed weights w in combination with the constrained search would lead to a universal functional of the one-particle reduced density matrix. To turn this into a viable functional theory, however, we also need to implement an exact convex relaxation. This general procedure includes Valone’s pioneering work on ground state RDMFT as the special case w = (1,0, ···). Then, we work out in a comprehensive manner a methodology for deriving a compact description of the functional’s domain. This leads to a hierarchy of generalized exclusion principle constraints which we illustrate in great detail. By anticipating their future pivotal role in functional theories and to keep our work self-contained, several required concepts from convex analysis are introduced and discussed.

arXiv:2106.03918•Journal of Chemical Theory and Computation•Dec 21, 2021•38 citations

Reduced Density Matrix Functional Theory for Bosons: Foundations and Applications

Julia Liebert

Density functional theory constitutes the workhorse of modern electronic structure calculations due to its favourable computational cost despite the fact that it usually fails to describe strongly correlated systems. A particularly promising approach to overcome those difficulties is reduced density matrix functional theory (RDMFT): It abandons the complexity of the $N$-particle wave function and at the same time explicitly allows for fractional occupation numbers. It is the goal of this thesis to initiate and establish a bosonic RDMFT for both ground state and excited state energy calculations. Motivated by the Onsager and Penrose criterion which identifies RDMFT as a particularly suitable approach to describe Bose-Einstein condensates (BECs), we derive the universal functional for a homogeneous BEC in the Bogoliubov regime. Remarkably, the gradient of the universal functional is found to diverge repulsively in the regime of complete condensation. This introduces the new concept of a BEC force, which provides a universal explanation for quantum depletion since it is merely based on the geometry of quantum states. In the second part of the thesis, we propose and work out an ensemble RDMFT targeting excitations in bosonic quantum systems. This endeavour further highlights the potential of convex analysis for the development of functional theories in the future. Indeed by resorting to several concepts from convex analysis, we succeeded to provide a comprehensive foundation of $\boldsymbol{w}$-ensemble RDMFT for bosons which is further based on a generalization of the Ritz variational principle and a constrained search formalism. In particular, we solve the emerging $N$-representability problem leading to a generalization of Pauli's famous exclusion principle to bosonic mixed states.

arXiv:2205.02635•May 5, 2022•1 citation

Deriving density-matrix functionals for excited states

Julia Liebert, Christian Schilling

We initiate the recently proposed w w -ensemble one-particle reduced density matrix functional theory ( w w -RDMFT) by deriving the first functional approximations and illustrate how excitation energies can be calculated in practice. For this endeavour, we first study the symmetric Hubbard dimer, constituting the building block of the Hubbard model, for which we execute the Levy-Lieb constrained search. Second, due to the particular suitability of w w -RDMFT for describing Bose-Einstein condensates, we demonstrate three conceptually different approaches for deriving the universal functional in a homogeneous Bose gas for arbitrary pair interaction in the Bogoliubov regime. Remarkably, in both systems the gradient of the functional is found to diverge repulsively at the boundary of the functional’s domain, extending the recently discovered Bose-Einstein condensation force to excited states. Our findings highlight the physical relevance of the generalized exclusion principle for fermionic and bosonic mixed states and the curse of universality in functional theories.

arXiv:2210.00964•SciPost Physics•May 19, 2023•15 citations

Communication: Relating the pure and ensemble density matrix functional

Christian Schilling

A crucial theorem in Reduced Density Matrix Functional Theory (RDMFT) suggests that the universal pure and ensemble functionals coincide on their common domain of pure N-representable one-matrices. We refute this by a comprehensive analysis of the geometric picture underlying Levy’s constrained search. Moreover, we then show that the ensemble functional follows instead as the lower convex envelop of the pure functional. It is particularly remarkable that the pure functional determines the ensemble functional even outside its own domain of pure N-representable one-matrices. From a general perspective, this demonstrates that relaxing pure RDMFT to ensemble RDMFT does not necessarily circumvent the complexity of the one-body pure N-representability conditions (generalized Pauli constraints). Instead, the complexity may simply be transferred from the underlying space of pure N-representable one-matrices to the structure of the universal one-matrix functional.

arXiv:1812.09291•The Journal of Chemical Physics•Dec 21, 2018•62 citations

Refining ensemble N -representability of one-body density matrices from partial information

Julia Liebert, Anna O. Schouten, Irma Avdic +2 more

Abstract The N -representability problem places fundamental constraints on reduced density matrices (RDMs) that originate from physical many-fermion quantum states. Motivated by recent developments in functional theories, we introduce a hierarchy of ensemble one-body N -representability problems that incorporate partial knowledge of the one-body RDMs (1RDMs) within an ensemble of N -fermion states with fixed weights w i . Specifically, we propose a systematic relaxation that reduces the refined problem—where full 1RDMs are fixed for certain ensemble elements—to a more tractable form involving only natural occupation number vectors. Remarkably, we show that this relaxed problem is related to a generalization of Horn’s problem, enabling an explicit solution by combining its constraints with those of the weighted ensemble N -representability conditions. An additional convex relaxation yields a convex polytope that provides physically meaningful restrictions on lattice site occupations in ensemble density functional theory for excited states.

arXiv:2506.09960•New Journal of Physics•Dec 1, 2025•0 citations

An exact one-particle theory of bosonic excitations: from a generalized Hohenberg–Kohn theorem to convexified N-representability

Julia Liebert, Christian Schilling

Abstract Motivated by the Penrose–Onsager criterion for Bose–Einstein condensation we propose a functional theory for targeting low-lying excitation energies of bosonic quantum systems through the one-particle picture. For this, we employ an extension of the Rayleigh–Ritz variational principle to ensemble states with spectrum w and prove a corresponding generalization of the Hohenberg–Kohn theorem: the underlying one-particle reduced density matrix determines all properties of systems of N identical particles in their w -ensemble states. Then, to circumvent the v-representability problem common to functional theories, and to deal with energetic degeneracies, we resort to the Levy–Lieb constrained search formalism in combination with an exact convex relaxation. The corresponding bosonic one-body w -ensemble N-representability problem is solved comprehensively. Remarkably, this reveals a complete hierarchy of bosonic exclusion principle constraints in conceptual analogy to Pauli’s exclusion principle for fermions and recently discovered generalizations thereof.

arXiv:2204.12715•New Journal of Physics•Jan 1, 2023•18 citations
We propose and work out a reduced density matrix functional theory (RDMFT) for calculating energies of eigenstates of interacting many-electron systems beyond the ground state. Various obstacles which historically have doomed such an approach to be unfeasible are overcome. First, we resort to a generalization of the Ritz variational principle to ensemble states with fixed weights. This in combination with the constrained search formalism allows us to establish a universal functional of the one-particle reduced density matrix. Second, we employ tools from convex analysis to circumvent the too involved N-representability constraints. Remarkably, this identifies Valone's pioneering work on RDMFT as a special case of convex relaxation and reveals that crucial information about the excitation structure is contained in the functional's domain. Third, to determine the crucial latter object, a methodology is developed which eventually leads to a generalized exclusion principle. The corresponding linear constraints are calculated for systems of arbitrary size.
In Phys. Rev. Lett . 2021, 127, 023001 a reduced density matrix functional theory (RDMFT) was proposed for calculating energies of selected eigenstates of interacting many-Fermion systems. Here, we develop a solid foundation for this so-called w -RDMFT and present the details of various derivations. First, we explain how a generalization of the Ritz variational principle to ensemble states with fixed weights w in combination with the constrained search would lead to a universal functional of the one-particle reduced density matrix. To turn this into a viable functional theory, however, we also need to implement an exact convex relaxation. This general procedure includes Valone’s pioneering work on ground state RDMFT as the special case w = (1,0, ···). Then, we work out in a comprehensive manner a methodology for deriving a compact description of the functional’s domain. This leads to a hierarchy of generalized exclusion principle constraints which we illustrate in great detail. By anticipating their future pivotal role in functional theories and to keep our work self-contained, several required concepts from convex analysis are introduced and discussed.
Density functional theory constitutes the workhorse of modern electronic structure calculations due to its favourable computational cost despite the fact that it usually fails to describe strongly correlated systems. A particularly promising approach to overcome those difficulties is reduced density matrix functional theory (RDMFT): It abandons the complexity of the NNN-particle wave function and at the same time explicitly allows for fractional occupation numbers. It is the goal of this thesis to initiate and establish a bosonic RDMFT for both ground state and excited state energy calculations. Motivated by the Onsager and Penrose criterion which identifies RDMFT as a particularly suitable approach to describe Bose-Einstein condensates (BECs), we derive the universal functional for a homogeneous BEC in the Bogoliubov regime. Remarkably, the gradient of the universal functional is found to diverge repulsively in the regime of complete condensation. This introduces the new concept of a BEC force, which provides a universal explanation for quantum depletion since it is merely based on the geometry of quantum states. In the second part of the thesis, we propose and work out an ensemble RDMFT targeting excitations in bosonic quantum systems. This endeavour further highlights the potential of convex analysis for the development of functional theories in the future. Indeed by resorting to several concepts from convex analysis, we succeeded to provide a comprehensive foundation of w\boldsymbol{w}w-ensemble RDMFT for bosons which is further based on a generalization of the Ritz variational principle and a constrained search formalism. In particular, we solve the emerging NNN-representability problem leading to a generalization of Pauli's famous exclusion principle to bosonic mixed states.
We initiate the recently proposed w w -ensemble one-particle reduced density matrix functional theory ( w w -RDMFT) by deriving the first functional approximations and illustrate how excitation energies can be calculated in practice. For this endeavour, we first study the symmetric Hubbard dimer, constituting the building block of the Hubbard model, for which we execute the Levy-Lieb constrained search. Second, due to the particular suitability of w w -RDMFT for describing Bose-Einstein condensates, we demonstrate three conceptually different approaches for deriving the universal functional in a homogeneous Bose gas for arbitrary pair interaction in the Bogoliubov regime. Remarkably, in both systems the gradient of the functional is found to diverge repulsively at the boundary of the functional’s domain, extending the recently discovered Bose-Einstein condensation force to excited states. Our findings highlight the physical relevance of the generalized exclusion principle for fermionic and bosonic mixed states and the curse of universality in functional theories.
A crucial theorem in Reduced Density Matrix Functional Theory (RDMFT) suggests that the universal pure and ensemble functionals coincide on their common domain of pure N-representable one-matrices. We refute this by a comprehensive analysis of the geometric picture underlying Levy’s constrained search. Moreover, we then show that the ensemble functional follows instead as the lower convex envelop of the pure functional. It is particularly remarkable that the pure functional determines the ensemble functional even outside its own domain of pure N-representable one-matrices. From a general perspective, this demonstrates that relaxing pure RDMFT to ensemble RDMFT does not necessarily circumvent the complexity of the one-body pure N-representability conditions (generalized Pauli constraints). Instead, the complexity may simply be transferred from the underlying space of pure N-representable one-matrices to the structure of the universal one-matrix functional.
Abstract The N -representability problem places fundamental constraints on reduced density matrices (RDMs) that originate from physical many-fermion quantum states. Motivated by recent developments in functional theories, we introduce a hierarchy of ensemble one-body N -representability problems that incorporate partial knowledge of the one-body RDMs (1RDMs) within an ensemble of N -fermion states with fixed weights w i . Specifically, we propose a systematic relaxation that reduces the refined problem—where full 1RDMs are fixed for certain ensemble elements—to a more tractable form involving only natural occupation number vectors. Remarkably, we show that this relaxed problem is related to a generalization of Horn’s problem, enabling an explicit solution by combining its constraints with those of the weighted ensemble N -representability conditions. An additional convex relaxation yields a convex polytope that provides physically meaningful restrictions on lattice site occupations in ensemble density functional theory for excited states.
Abstract Motivated by the Penrose–Onsager criterion for Bose–Einstein condensation we propose a functional theory for targeting low-lying excitation energies of bosonic quantum systems through the one-particle picture. For this, we employ an extension of the Rayleigh–Ritz variational principle to ensemble states with spectrum w and prove a corresponding generalization of the Hohenberg–Kohn theorem: the underlying one-particle reduced density matrix determines all properties of systems of N identical particles in their w -ensemble states. Then, to circumvent the v-representability problem common to functional theories, and to deal with energetic degeneracies, we resort to the Levy–Lieb constrained search formalism in combination with an exact convex relaxation. The corresponding bosonic one-body w -ensemble N-representability problem is solved comprehensively. Remarkably, this reveals a complete hierarchy of bosonic exclusion principle constraints in conceptual analogy to Pauli’s exclusion principle for fermions and recently discovered generalizations thereof.