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:2301.10193 · The Journal of Chemical Physics

Refining and relating fundamentals of functional theory

Julia Liebert, Adam Yanis Chaou, Christian Schilling

Published Jun 1, 2023• 27 citations
arXivHTMLPDF

Abstract

To advance the foundation of one-particle reduced density matrix functional theory (1RDMFT), we refine and relate some of its fundamental features and underlying concepts. We define by concise means the scope of a 1RDMFT, identify its possible natural variables, and explain how symmetries could be exploited. In particular, for systems with time-reversal symmetry, we explain why there exist six equivalent universal functionals, prove concise relations among them, and conclude that the important notion of v-representability is relative to the scope and choice of variable. All these fundamental concepts are then comprehensively discussed and illustrated for the Hubbard dimer and its generalization to arbitrary pair interactions W. For this, we derive by analytical means the pure and ensemble functionals with respect to both the real- and complex-valued Hilbert space. The comparison of various functionals allows us to solve the underlying v-representability problems analytically, and the dependence of its solution on the pair interaction is demonstrated. Intriguingly, the gradient of each universal functional is found to always diverge repulsively on the boundary of the domain. In that sense, this key finding emphasizes the universal character of the fermionic exchange force, recently discovered and proven in the context of translationally invariant one-band lattice models.

DOI
10.1063/5.0143657

Related papers

Diverging Exchange Force and Form of the Exact Density Matrix Functional

Christian Schilling, R. Schilling

For translationally invariant one-band lattice models, we exploit the ab initio knowledge of the natural orbitals to simplify reduced density matrix functional theory (RDMFT). Striking underlying features are discovered. First, within each symmetry sector, the interaction functional $\mathcal{F}$ depends only on the natural occupation numbers $\mathbit{n}$. The respective sets ${\mathcal{P}}_{N}^{1}$ and ${\mathcal{E}}_{N}^{1}$ of pure and ensemble $N$-representable one-matrices coincide. Second, and most importantly, the exact functional is strongly shaped by the geometry of the polytope ${\mathcal{E}}_{N}^{1}\ensuremath{\equiv}{\mathcal{P}}_{N}^{1}$, described by linear constraints ${D}^{(j)}(\mathbit{n})\ensuremath{\ge}0$. For smaller systems, it follows as $\mathcal{F}[\mathbit{n}]=\ensuremath{\sum}_{i,{i}^{\ensuremath{'}}}{\overline{V}}_{i,{i}^{\ensuremath{'}}}\sqrt{{D}^{(i)}(\mathbit{n}){D}^{({i}^{\ensuremath{'}})}(\mathbit{n})}$. This generalizes to systems of arbitrary size by replacing each ${D}^{(i)}$ by a linear combination of ${{D}^{(j)}(\mathbit{n})}$ and adding a nonanalytical term involving the interaction $\stackrel{^}{V}$. Third, the gradient $d\mathcal{F}/d\mathbit{n}$ is shown to diverge on the boundary $\ensuremath{\partial}{\mathcal{E}}_{N}^{1}$, suggesting that the fermionic exchange symmetry manifests itself within RDMFT in the form of an ``exchange force.'' All findings hold for systems with a nonfixed particle number as well and $\stackrel{^}{V}$ can be any $p$-particle interaction. As an illustration, we derive the exact functional for the Hubbard square.

arXiv:1901.01321•Physical Review Letters•Jan 4, 2019•58 citations

Geometry of Generalized Density Functional Theories

Chih-Chun Wang

Density functional theory (DFT) is an indispensable ab initio method in both quantum chemistry and condensed matter physics. Based on recent advancements in reduced density matrix functional theory (RDMFT), a variant of DFT that is believed to be better suited for strongly correlated systems, we construct a mathematical framework generalizing all ground state functional theories, which in particular applies to fermionic, bosonic, and spin systems. Within the special class of such functional theories where the space of external potentials forms the Lie algebra of a compact Lie group, the $N$-representability problem is readily solved by applying techniques from the study of momentum maps in symplectic geometry, an approach complementary to Klyachko's famous solution to the quantum marginal problem. The ``boundary force'', a diverging repulsive force from the boundary of the functional's domain observed in previous works but only qualitatively understood in isolated systems, is studied quantitatively and extensively in our work. Specifically, we present a formula capturing the exact behavior of the functional close to the boundary. In the case where the Lie algebra is abelian, a completely rigorous proof of the boundary force formula based on Levy-Lieb constrained search is given. Our formula is a first step towards developing more accurate functional approximations, with the potential of improving current RDMFT approximate functionals such as the Piris natural orbital functionals (NOFs). All key concepts and ideas of our work are demonstrated in translation invariant bosonic lattice systems.

quant-pharXiv:2511.14822•Nov 19, 2025•0 citations

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

Unified Framework for Functional Theories of Quantum Systems

Chih-Chun Wang, Julia Liebert, Markus Penz +1 more

We introduce and study a unified framework for density-functional theory and its variants for quantum systems on finite-dimensional Hilbert spaces. These theories seek to reduce the complexity inherent in the many-body quantum problem by describing ground states through reduced variables. The central ingredients of our unified framework are a generalized choice of basic observables, whose expectation values define precisely those reduced variables, and a fixed part of the Hamiltonian characterizing the class of quantum systems under consideration. It is this minimal structure, which we call the scope of a functional theory, that is necessary and sufficient for the formulation of a functional theory. In particular, it allows one to define the universal functionals, establish their convexity and differentiability properties, address representability questions, and prove a Hohenberg-Kohn-type uniqueness result. A purification construction also relates ensemble and weighted-ensemble functionals to the pure-state variant. Particular emphasis is placed on functional theories with Lie-algebra observable structures, connecting the variational framework to symplectic geometry. The result of this work is a systematic mathematical formulation in which structural results can be proved once and applied across a broad class of finite-dimensional functional theories.

math-phCross-listedarXiv:2606.06676•Jun 8, 2026•0 citations

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
To advance the foundation of one-particle reduced density matrix functional theory (1RDMFT), we refine and relate some of its fundamental features and underlying concepts. We define by concise means the scope of a 1RDMFT, identify its possible natural variables, and explain how symmetries could be exploited. In particular, for systems with time-reversal symmetry, we explain why there exist six equivalent universal functionals, prove concise relations among them, and conclude that the important notion of v-representability is relative to the scope and choice of variable. All these fundamental concepts are then comprehensively discussed and illustrated for the Hubbard dimer and its generalization to arbitrary pair interactions W. For this, we derive by analytical means the pure and ensemble functionals with respect to both the real- and complex-valued Hilbert space. The comparison of various functionals allows us to solve the underlying v-representability problems analytically, and the dependence of its solution on the pair interaction is demonstrated. Intriguingly, the gradient of each universal functional is found to always diverge repulsively on the boundary of the domain. In that sense, this key finding emphasizes the universal character of the fermionic exchange force, recently discovered and proven in the context of translationally invariant one-band lattice models.
For translationally invariant one-band lattice models, we exploit the ab initio knowledge of the natural orbitals to simplify reduced density matrix functional theory (RDMFT). Striking underlying features are discovered. First, within each symmetry sector, the interaction functional F\mathcal{F}F depends only on the natural occupation numbers $\mathbit{n}$. The respective sets PN1{\mathcal{P}}_{N}^{1}PN1​ and EN1{\mathcal{E}}_{N}^{1}EN1​ of pure and ensemble NNN-representable one-matrices coincide. Second, and most importantly, the exact functional is strongly shaped by the geometry of the polytope ${\mathcal{E}}_{N}^{1}\ensuremath{\equiv}{\mathcal{P}}_{N}^{1}$, described by linear constraints ${D}^{(j)}(\mathbit{n})\ensuremath{\ge}0$. For smaller systems, it follows as $\mathcal{F}[\mathbit{n}]=\ensuremath{\sum}_{i,{i}^{\ensuremath{'}}}{\overline{V}}_{i,{i}^{\ensuremath{'}}}\sqrt{{D}^{(i)}(\mathbit{n}){D}^{({i}^{\ensuremath{'}})}(\mathbit{n})}$. This generalizes to systems of arbitrary size by replacing each D(i){D}^{(i)}D(i) by a linear combination of ${{D}^{(j)}(\mathbit{n})}$ and adding a nonanalytical term involving the interaction $\stackrel{^}{V}$. Third, the gradient $d\mathcal{F}/d\mathbit{n}$ is shown to diverge on the boundary $\ensuremath{\partial}{\mathcal{E}}_{N}^{1}$, suggesting that the fermionic exchange symmetry manifests itself within RDMFT in the form of an ``exchange force.'' All findings hold for systems with a nonfixed particle number as well and $\stackrel{^}{V}$ can be any ppp-particle interaction. As an illustration, we derive the exact functional for the Hubbard square.
Density functional theory (DFT) is an indispensable ab initio method in both quantum chemistry and condensed matter physics. Based on recent advancements in reduced density matrix functional theory (RDMFT), a variant of DFT that is believed to be better suited for strongly correlated systems, we construct a mathematical framework generalizing all ground state functional theories, which in particular applies to fermionic, bosonic, and spin systems. Within the special class of such functional theories where the space of external potentials forms the Lie algebra of a compact Lie group, the NNN-representability problem is readily solved by applying techniques from the study of momentum maps in symplectic geometry, an approach complementary to Klyachko's famous solution to the quantum marginal problem. The ``boundary force'', a diverging repulsive force from the boundary of the functional's domain observed in previous works but only qualitatively understood in isolated systems, is studied quantitatively and extensively in our work. Specifically, we present a formula capturing the exact behavior of the functional close to the boundary. In the case where the Lie algebra is abelian, a completely rigorous proof of the boundary force formula based on Levy-Lieb constrained search is given. Our formula is a first step towards developing more accurate functional approximations, with the potential of improving current RDMFT approximate functionals such as the Piris natural orbital functionals (NOFs). All key concepts and ideas of our work are demonstrated in translation invariant bosonic lattice systems.
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.
We introduce and study a unified framework for density-functional theory and its variants for quantum systems on finite-dimensional Hilbert spaces. These theories seek to reduce the complexity inherent in the many-body quantum problem by describing ground states through reduced variables. The central ingredients of our unified framework are a generalized choice of basic observables, whose expectation values define precisely those reduced variables, and a fixed part of the Hamiltonian characterizing the class of quantum systems under consideration. It is this minimal structure, which we call the scope of a functional theory, that is necessary and sufficient for the formulation of a functional theory. In particular, it allows one to define the universal functionals, establish their convexity and differentiability properties, address representability questions, and prove a Hohenberg-Kohn-type uniqueness result. A purification construction also relates ensemble and weighted-ensemble functionals to the pure-state variant. Particular emphasis is placed on functional theories with Lie-algebra observable structures, connecting the variational framework to symplectic geometry. The result of this work is a systematic mathematical formulation in which structural results can be proved once and applied across a broad class of finite-dimensional functional theories.
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.