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

Disentangling the marginal problem in quantum chemistry

Carlos L. Benavides-Riveros, José M. Gracia Bondía

Published Jan 1, 2015• 4 citations
OpenAlex

Abstract

It is well known that determining the energy of molecules and other quantum many-body systems reduces in the standard approximation to optimizing a simple linear functional of a 12-variable object, the two-electron reduced density matrix (2-RDM). The difficulty is that the variation ensemble for that functional has never been satisfactorily determined. This is known as the N-representability problem of quantum chemistry (which to a large extent is a problem of quantum information theory). The situation has given rise to competing research programs, typically trading more complicated functionals for simpler representability conditions. Chief among them is density functional theory, based on a three-variable object for which the N-representability is trivial, whereas the exact functional is very strange indeed, and probably forever unknowable. An intermediate position is occupied by 1-RDM functional theory. Postulated by Pauli to explain the electronic structure of atoms and molecules, the exclusion principle establishes an upper bound of 1 for the fermionic natural occupation numbers, accordingly allowing no more than one electron in each quantum state. This is a necessary and sufficient condition for a 1-RDM to be the contraction of an ensemble N-body density matrix. The fermionic one-body quantum marginal problem asks whether given natural occupation numbers can arise from an antisymmetric (ensemble or pure) N-particle state. The configuration interaction method affords optimal descriptions of quantum states of atoms and molecules by expanding the wave function in terms of orbital-based configurations of Slater determinants. For these systems, the dimension of the Hilbert space grows binomially with m, the number of spin-orbitals in the basis of the one-particle Hilbert space, and N, the number of electrons of the system. It has been observed that for the rank-six approximation of a pure-state N=3 system, the occupation numbers satisfy some additional constraints, stronger than the Pauli exclusion principle. The recent analysis by Alexander Klyachko and coworkers of the one-body marginal problem of the pure N-fermion state establishes a systematic approach to this type of constraints. In fact, for a pure quantum system of N electrons arranged in m spin-orbitals the occupation numbers satisfy a set of linear inequalities, known as generalized Pauli constraints (GPC). These inequalities define a convex polytope of allowed states in R^m. They are conditions for a 1-RDM to be the contraction of a N-body state. When one of the GPC is completely saturated, the system is said to be pinned, and it lies on one of the facets of the polytope. The nature of those conditions has been explored till now only in a few systems: a model of three spinless fermions confined to a one-dimensional harmonic potential, the lithium isoelectronic series and ground and excited states of some three- and four-electron molecules for the rank being at least twice the number of electrons. For all these systems the inequalities are (quite often) nearly saturated. This is the so-called quasipinning phenomenon. In this PhD thesis we have taken up the challenge of using numerical and analytical methods to examine pinned and quasipinned states, for atoms and molecules, starting from scratch with configuration-interaction and multiconfiguration self-consistent methods. This procedure serves to study the occurrence of quasipinning in realistic systems. A second goal is to show how the subsets of pinned states defined by GPC give rise to the most efficient approach from a computational viewpoint, yielding the leading order of the electron-electron correlations. As a consequence, we underline in this thesis the theoretical and practical importance of Klyachko's approach to the quantum marginal problem and its impact on the competing research programs to determine feasible electronic densities, 1-RDM, 2-RDM, intracular distributions or Wigner density quasiprobabilities. In relation with the above, our research provides a new variational optimization method for few-fermion ground states. We quantitatively confirm its high accuracy for quasipinned systems and derive an upper bound on the error of the correlation energy given by the ratio of the numerical value of the Klyachko inequality and the distance to the Hartree-Fock point. Depending on the details of the algorithm, we are able to reach 98%-99% of the correlation energy for such systems.

Related papers

Generalized Pauli constraints in reduced density matrix functional theory

Iris Theophilou, Nektarios N. Lathiotakis, Miguel A. L. Marques +1 more

Functionals of the one-body reduced density matrix (1-RDM) are routinely minimized under Coleman’s ensemble N-representability conditions. Recently, the topic of pure-state N-representability conditions, also known as generalized Pauli constraints, received increased attention following the discovery of a systematic way to derive them for any number of electrons and any finite dimensionality of the Hilbert space. The target of this work is to assess the potential impact of the enforcement of the pure-state conditions on the results of reduced density-matrix functional theory calculations. In particular, we examine whether the standard minimization of typical 1-RDM functionals under the ensemble N-representability conditions violates the pure-state conditions for prototype 3-electron systems. We also enforce the pure-state conditions, in addition to the ensemble ones, for the same systems and functionals and compare the correlation energies and optimal occupation numbers with those obtained by the enforcement of the ensemble conditions alone.

•Jan 1, 2015•34 citations

Generalized Pauli constraints in reduced density matrix functional theory

Iris Theophilou, Nektarios N. Lathiotakis, Miguel A. L. Marques +1 more

Functionals of the one-body reduced density matrix (1-RDM) are routinely minimized under Coleman's ensemble N-representability conditions. Recently, the topic of pure-state N-representability conditions, also known as generalized Pauli constraints, received increased attention following the discovery of a systematic way to derive them for any number of electrons and any finite dimensionality of the Hilbert space. The target of this work is to assess the potential impact of the enforcement of the pure-state conditions on the results of reduced density-matrix functional theory calculations. In particular, we examine whether the standard minimization of typical 1-RDM functionals under the ensemble N-representability conditions violates the pure-state conditions for prototype 3-electron systems. We also enforce the pure-state conditions, in addition to the ensemble ones, for the same systems and functionals and compare the correlation energies and optimal occupation numbers with those obtained by the enforcement of the ensemble conditions alone.

arXiv:1503.00746•The Journal of Chemical Physics•Apr 16, 2015•52 citations

Sparsity of the wavefunction from the generalized Pauli exclusion principle

Romit Chakraborty, David A. Mazziotti

Electron occupations that arise from pure quantum states are restricted by a stringent set of conditions that are said to generalize the Pauli exclusion principle. These generalized Pauli constraints (GPCs) define the boundary of the set of one-electron reduced density matrices (1-RDMs) that are derivable from at least one N-electron wavefunction. In this paper, we investigate the sparsity of the Slater-determinant representation of the wavefunction that is a necessary, albeit not sufficient, condition for its 1-RDM to lie on the boundary of the set of pure N-representable 1-RDMs or in other words saturate one of the GPCs. The sparse wavefunction, we show, is exact not only for 3 electrons in 6 orbitals but also for 3 electrons in 8 orbitals. For larger numbers of electrons and/or orbitals in the lowest spin state, the exact wavefunction does not generally saturate one of the GPCs, and hence, the sparse representation is typically an approximation. Because the sparsity of the wavefunction is a necessary but not sufficient condition for saturation of one of the GPCs, optimization of the sparse wavefunction Ansatz to minimize the ground-state energy does not necessarily produce a wavefunction whose 1-RDM exactly saturates one of the GPCs. While the sparse Ansatz can be employed with arbitrary orbitals or optimized orbitals, in this paper, we explore the Ansatz with the natural orbitals from full configuration interaction, which yields an upper bound to the ground-state energy that equals the exact energy for a given basis set if the full-configuration-interaction wavefunction saturates the Ansatz’s GPC. With calculations on the boron isoelectronic sequence, the dinitrogen cation N2+, hydrogen chains, and cyclic conjugated π systems, we examine the quality of the sparse wavefunction Ansatz from the amount of correlation energy recovered.

The Journal of Chemical Physics•Feb 5, 2018•12 citations

Structure of the first order reduced density matrix in three electron systems: A generalized Pauli constraints assisted study

Iris Theophilou, Nektarios N. Lathiotakis, Nicole Helbig

We investigate the structure of the one-body reduced density matrix of three electron systems, i.e., doublet and quadruplet spin configurations, corresponding to the smallest interacting system with an open-shell ground state. To this end, we use configuration interaction (CI) expansions of the exact wave function in Slater determinants built from natural orbitals in a finite dimensional Hilbert space. With the exception of maximally polarized systems, the natural orbitals of spin eigenstates are generally spin dependent, i.e., the spatial parts of the up and down natural orbitals form two different sets. A measure to quantify this spin dependence is introduced and it is shown that it varies by several orders of magnitude depending on the system. We also study the ordering issue of the spin-dependent occupation numbers which has practical implications in reduced density matrix functional theory minimization schemes, when generalized Pauli constraints (GPCs) are imposed and in the form of the CI expansion in terms of the natural orbitals. Finally, we discuss the aforementioned CI expansion when there are GPCs that are almost "pinned."

arXiv:1712.08670•The Journal of Chemical Physics•Mar 20, 2018•4 citations

Solving one-body ensemble N-representability problems with spin

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

The Pauli exclusion principle is fundamental to understanding electronic quantum systems. It namely constrains the expected occupancies n i of orbitals φ i according to 0 ≤ n i ≤ 2 . In this work, we first refine the underlying one-body N -representability problem by taking into account simultaneously spin symmetries and a potential degree of mixedness w of the N -electron quantum state. We then derive a comprehensive solution to this problem by using basic tools from representation theory, convex analysis and discrete geometry. Specifically, we show that the set of admissible orbital one-body reduced density matrices is fully characterized by linear spectral constraints on the natural orbital occupation numbers, defining a convex polytope Σ N , S ( w ) ⊂ [ 0 , 2 ] d . These constraints are independent of M and the number d of orbitals, while their dependence on N , S is linear, and we can thus calculate them for arbitrary system sizes and spin quantum numbers. Our results provide a crucial missing cornerstone for ensemble density (matrix) functional theory.

arXiv:2412.01805•Quantum•Dec 2, 2025•0 citations

Structure of the first order reduced density matrix in three electron systems: A generalized Pauli constraints assisted study

Theophilou, I., Lathiotakis, N., Helbig, N.

We investigate the structure of the one-body reduced density matrix of three electron systems, i.e., doublet and quadruplet spin configurations, corresponding to the smallest interacting system with an open-shell ground state. To this end, we use configuration interaction (CI) expansions of the exact wave function in Slater determinants built from natural orbitals in a finite dimensional Hilbert space. With the exception of maximally polarized systems, the natural orbitals of spin eigenstates are generally spin dependent, i.e., the spatial parts of the up and down natural orbitals form two different sets. A measure to quantify this spin dependence is introduced and it is shown that it varies by several orders of magnitude depending on the system. We also study the ordering issue of the spin-dependent occupation numbers which has practical implications in reduced density matrix functional theory minimization schemes, when generalized Pauli constraints (GPCs) are imposed and in the form of the CI expansion in terms of the natural orbitals. Finally, we discuss the aforementioned CI expansion when there are GPCs that are almost “pinned.”

•Mar 21, 2018•4 citations
It is well known that determining the energy of molecules and other quantum many-body systems reduces in the standard approximation to optimizing a simple linear functional of a 12-variable object, the two-electron reduced density matrix (2-RDM). The difficulty is that the variation ensemble for that functional has never been satisfactorily determined. This is known as the N-representability problem of quantum chemistry (which to a large extent is a problem of quantum information theory). The situation has given rise to competing research programs, typically trading more complicated functionals for simpler representability conditions. Chief among them is density functional theory, based on a three-variable object for which the N-representability is trivial, whereas the exact functional is very strange indeed, and probably forever unknowable. An intermediate position is occupied by 1-RDM functional theory. Postulated by Pauli to explain the electronic structure of atoms and molecules, the exclusion principle establishes an upper bound of 1 for the fermionic natural occupation numbers, accordingly allowing no more than one electron in each quantum state. This is a necessary and sufficient condition for a 1-RDM to be the contraction of an ensemble N-body density matrix. The fermionic one-body quantum marginal problem asks whether given natural occupation numbers can arise from an antisymmetric (ensemble or pure) N-particle state. The configuration interaction method affords optimal descriptions of quantum states of atoms and molecules by expanding the wave function in terms of orbital-based configurations of Slater determinants. For these systems, the dimension of the Hilbert space grows binomially with m, the number of spin-orbitals in the basis of the one-particle Hilbert space, and N, the number of electrons of the system. It has been observed that for the rank-six approximation of a pure-state N=3 system, the occupation numbers satisfy some additional constraints, stronger than the Pauli exclusion principle. The recent analysis by Alexander Klyachko and coworkers of the one-body marginal problem of the pure N-fermion state establishes a systematic approach to this type of constraints. In fact, for a pure quantum system of N electrons arranged in m spin-orbitals the occupation numbers satisfy a set of linear inequalities, known as generalized Pauli constraints (GPC). These inequalities define a convex polytope of allowed states in R^m. They are conditions for a 1-RDM to be the contraction of a N-body state. When one of the GPC is completely saturated, the system is said to be pinned, and it lies on one of the facets of the polytope. The nature of those conditions has been explored till now only in a few systems: a model of three spinless fermions confined to a one-dimensional harmonic potential, the lithium isoelectronic series and ground and excited states of some three- and four-electron molecules for the rank being at least twice the number of electrons. For all these systems the inequalities are (quite often) nearly saturated. This is the so-called quasipinning phenomenon. In this PhD thesis we have taken up the challenge of using numerical and analytical methods to examine pinned and quasipinned states, for atoms and molecules, starting from scratch with configuration-interaction and multiconfiguration self-consistent methods. This procedure serves to study the occurrence of quasipinning in realistic systems. A second goal is to show how the subsets of pinned states defined by GPC give rise to the most efficient approach from a computational viewpoint, yielding the leading order of the electron-electron correlations. As a consequence, we underline in this thesis the theoretical and practical importance of Klyachko's approach to the quantum marginal problem and its impact on the competing research programs to determine feasible electronic densities, 1-RDM, 2-RDM, intracular distributions or Wigner density quasiprobabilities. In relation with the above, our research provides a new variational optimization method for few-fermion ground states. We quantitatively confirm its high accuracy for quasipinned systems and derive an upper bound on the error of the correlation energy given by the ratio of the numerical value of the Klyachko inequality and the distance to the Hartree-Fock point. Depending on the details of the algorithm, we are able to reach 98%-99% of the correlation energy for such systems.
Functionals of the one-body reduced density matrix (1-RDM) are routinely minimized under Coleman’s ensemble N-representability conditions. Recently, the topic of pure-state N-representability conditions, also known as generalized Pauli constraints, received increased attention following the discovery of a systematic way to derive them for any number of electrons and any finite dimensionality of the Hilbert space. The target of this work is to assess the potential impact of the enforcement of the pure-state conditions on the results of reduced density-matrix functional theory calculations. In particular, we examine whether the standard minimization of typical 1-RDM functionals under the ensemble N-representability conditions violates the pure-state conditions for prototype 3-electron systems. We also enforce the pure-state conditions, in addition to the ensemble ones, for the same systems and functionals and compare the correlation energies and optimal occupation numbers with those obtained by the enforcement of the ensemble conditions alone.
Functionals of the one-body reduced density matrix (1-RDM) are routinely minimized under Coleman's ensemble N-representability conditions. Recently, the topic of pure-state N-representability conditions, also known as generalized Pauli constraints, received increased attention following the discovery of a systematic way to derive them for any number of electrons and any finite dimensionality of the Hilbert space. The target of this work is to assess the potential impact of the enforcement of the pure-state conditions on the results of reduced density-matrix functional theory calculations. In particular, we examine whether the standard minimization of typical 1-RDM functionals under the ensemble N-representability conditions violates the pure-state conditions for prototype 3-electron systems. We also enforce the pure-state conditions, in addition to the ensemble ones, for the same systems and functionals and compare the correlation energies and optimal occupation numbers with those obtained by the enforcement of the ensemble conditions alone.
Electron occupations that arise from pure quantum states are restricted by a stringent set of conditions that are said to generalize the Pauli exclusion principle. These generalized Pauli constraints (GPCs) define the boundary of the set of one-electron reduced density matrices (1-RDMs) that are derivable from at least one N-electron wavefunction. In this paper, we investigate the sparsity of the Slater-determinant representation of the wavefunction that is a necessary, albeit not sufficient, condition for its 1-RDM to lie on the boundary of the set of pure N-representable 1-RDMs or in other words saturate one of the GPCs. The sparse wavefunction, we show, is exact not only for 3 electrons in 6 orbitals but also for 3 electrons in 8 orbitals. For larger numbers of electrons and/or orbitals in the lowest spin state, the exact wavefunction does not generally saturate one of the GPCs, and hence, the sparse representation is typically an approximation. Because the sparsity of the wavefunction is a necessary but not sufficient condition for saturation of one of the GPCs, optimization of the sparse wavefunction Ansatz to minimize the ground-state energy does not necessarily produce a wavefunction whose 1-RDM exactly saturates one of the GPCs. While the sparse Ansatz can be employed with arbitrary orbitals or optimized orbitals, in this paper, we explore the Ansatz with the natural orbitals from full configuration interaction, which yields an upper bound to the ground-state energy that equals the exact energy for a given basis set if the full-configuration-interaction wavefunction saturates the Ansatz’s GPC. With calculations on the boron isoelectronic sequence, the dinitrogen cation N2+, hydrogen chains, and cyclic conjugated π systems, we examine the quality of the sparse wavefunction Ansatz from the amount of correlation energy recovered.
We investigate the structure of the one-body reduced density matrix of three electron systems, i.e., doublet and quadruplet spin configurations, corresponding to the smallest interacting system with an open-shell ground state. To this end, we use configuration interaction (CI) expansions of the exact wave function in Slater determinants built from natural orbitals in a finite dimensional Hilbert space. With the exception of maximally polarized systems, the natural orbitals of spin eigenstates are generally spin dependent, i.e., the spatial parts of the up and down natural orbitals form two different sets. A measure to quantify this spin dependence is introduced and it is shown that it varies by several orders of magnitude depending on the system. We also study the ordering issue of the spin-dependent occupation numbers which has practical implications in reduced density matrix functional theory minimization schemes, when generalized Pauli constraints (GPCs) are imposed and in the form of the CI expansion in terms of the natural orbitals. Finally, we discuss the aforementioned CI expansion when there are GPCs that are almost "pinned."
The Pauli exclusion principle is fundamental to understanding electronic quantum systems. It namely constrains the expected occupancies n i of orbitals φ i according to 0 ≤ n i ≤ 2 . In this work, we first refine the underlying one-body N -representability problem by taking into account simultaneously spin symmetries and a potential degree of mixedness w of the N -electron quantum state. We then derive a comprehensive solution to this problem by using basic tools from representation theory, convex analysis and discrete geometry. Specifically, we show that the set of admissible orbital one-body reduced density matrices is fully characterized by linear spectral constraints on the natural orbital occupation numbers, defining a convex polytope Σ N , S ( w ) ⊂ [ 0 , 2 ] d . These constraints are independent of M and the number d of orbitals, while their dependence on N , S is linear, and we can thus calculate them for arbitrary system sizes and spin quantum numbers. Our results provide a crucial missing cornerstone for ensemble density (matrix) functional theory.
We investigate the structure of the one-body reduced density matrix of three electron systems, i.e., doublet and quadruplet spin configurations, corresponding to the smallest interacting system with an open-shell ground state. To this end, we use configuration interaction (CI) expansions of the exact wave function in Slater determinants built from natural orbitals in a finite dimensional Hilbert space. With the exception of maximally polarized systems, the natural orbitals of spin eigenstates are generally spin dependent, i.e., the spatial parts of the up and down natural orbitals form two different sets. A measure to quantify this spin dependence is introduced and it is shown that it varies by several orders of magnitude depending on the system. We also study the ordering issue of the spin-dependent occupation numbers which has practical implications in reduced density matrix functional theory minimization schemes, when generalized Pauli constraints (GPCs) are imposed and in the form of the CI expansion in terms of the natural orbitals. Finally, we discuss the aforementioned CI expansion when there are GPCs that are almost “pinned.”