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:1310.1246 · Physical Review A

Fermion N -representability for prescribed density and paramagnetic current density

Erik I. Tellgren, Simen Kvaal, Trygve Ulf Helgaker

Published Jan 30, 2014• 21 citations
arXivHTMLPDF

Abstract

The $N$-representability problem is the problem of determining whether there exists $N$-particle states with some prescribed property. Here we report an affirmative solution to the fermion $N$-representability problem when both the density and the paramagnetic current density are prescribed. This problem arises in current-density functional theory and is a generalization of the well-studied corresponding problem (only the density prescribed) in density functional theory. Given any density and paramagnetic current density satisfying a minimal regularity condition (essentially that a von Weiz\"acker--like canonical kinetic energy density is locally integrable), we prove that there exists a corresponding $N$-particle state. We prove this by constructing an explicit one-particle reduced density matrix in the form of a position-space kernel, i.e., a function of two continuous-position variables. In order to make minimal assumptions, we also address mathematical subtleties regarding the diagonal of, and how to rigorously extract paramagnetic current densities from, one-particle reduced density matrices in kernel form.

DOI
10.1103/physreva.89.012515

Related papers

N -Representability Problem for Fermion Density Matrices. II. The First-Order Density Matrix with N Even

Darwin W. Smith

For a quantum-mechanical system of $N$ identical fermions, the $N$-representability problem is the problem of recognizing whether, for a given pth-order reduced density matrix ${\ensuremath{\Gamma}}^{(p)}(12\ensuremath{\cdots}p|{1}^{\ensuremath{'}}{2}^{\ensuremath{'}}\ensuremath{\cdots}{p}^{\ensuremath{'}})$, there exists an antisymmetric $N$-particle wave function $\ensuremath{\Psi}(12\ensuremath{\cdots}N)$ such that ${\ensuremath{\Gamma}}^{(p)}(12\ensuremath{\cdots}p|{1}^{\ensuremath{'}}{2}^{\ensuremath{'}}\ensuremath{\cdots}{p}^{\ensuremath{'}})=({N}{p})\ensuremath{\int}\ensuremath{\Psi}(12\ensuremath{\cdots}N)\ifmmode\times\else\texttimes\fi{}{\ensuremath{\Psi}}^{*}({1}^{\ensuremath{'}}{2}^{\ensuremath{'}}\ensuremath{\cdots}{p}^{\ensuremath{'}},p+1\ensuremath{\cdots}N){d}_{{\ensuremath{\tau}}_{p+1}}\ensuremath{\cdots}{d}_{{\ensuremath{\tau}}_{N}}$. It is shown that if the Hamiltonian of a system is time-reversal invariant, and the number of particles, $N$, is even, the necessary and sufficient condition that an approximate first-order density matrix corresponding to a nondegenerate energy eigenstate be $N$-representable is that its natural spin-orbital occupation numbers be equal in pairs.

Physical Review•Jul 29, 1966•67 citations

N-representability of diagonal elements of second-order reduced density matrices

S.Kh. Samvelyan

The problem of pure-state N-representability of the two-particle spin-dependent density function ρ(x1, x2) is considered for an N-electron system, and a procedure for finding an N-representable ρ(x1, x2) is advanced. The problem is formulated in the framework of a family of N × N matrices formed from integrals of auxiliary two-particle functions θn(x1, x2) converging at n → ∞ to ρ(x1, x2)/[N(N−1)]. The simple requirement of positive definiteness of these matrices is shown to play a decisive role in finding an N-representable ρ(x1, x2). The results obtained may open new possibilities for using ρ(x1, x2) in the density-functional theory. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 65: 127–142, 1997

International Journal of Quantum Chemistry•Jan 1, 1997•32 citations

Diagonal N-representability as a method for solving the representability problem

M. M. Mestechkin

The latest results of density matrix theory are summarized in terms of analysis of off-diagonal components. For a given reduced density matrix (RDM) of arbitrary order, the full density matrix is explicitly constructed, from which the RDM is derived by means of the corresponding reduction procedure. Both the matrices satisfy all the necessary conditions, such antisymmetry, normalization, hermiticity, except that the non-negative definiteness of the full DM is not guaranteed. If an RDM satisfies the necessary inequalities with the same number of indices as it has itself, then, as proved here, there is an additive correction to the full DM matrix that turns it into a non-negative definite matrix leading to the same RDM and satisfying the other requirements. This establishes a fundamental fact that these inequalities are sufficient as well, providing the N-representability of the given RDM. The construction is based on the classification, introduced in the article, of the multi-index elements of the RDM and full DM according to their degree of off-diagonality. RDM elements with a predetermined degree of off-diagonality are expressed through full DM elements with the same off-diagonality. This shows that the conventional variational approach to calculating the energy of the system by convex programming methods under the conditions of N-representability (before and even after reaching the minimum energy) generally does not lead to constructing a mixed quantum state composed of solutions to the Schrödinger equation. Therefore, except for the energy, it does not guarantee a correct description of the other properties so long as the conditions of correctness for each particular property are not expressed in terms of the RDM and not incorporated into the specific conditions of N-representability (regardless of the requirement of energy minimization).

Russian Journal of Physical Chemistry B•Mar 1, 2017•0 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

Density matrices and density functionals : proceedings of the A. John Coleman symposium

Alan Coleman, Robert Erdahl, Vedene H. Smith

A Tribute to A. John Coleman - The Tame Mathematician.- Reduced Density Matrices: 1929-1989.- Some Aspects on the Development of the Theory of Reduced Density Matrices and the Representability Problem.- Representability Conditions.- On the Diagonal N-Representability Problem.- Fermion N-Representability Conditions Generated by a Decomposition of the 1-Particle Identity Operator onto Mutually Orthogonal Projection Operators.- The Unitarily Invariant Decomposition of Hermitian Operators.- Building Up N-Electron States with Symplectic Symmetry.- Time Dependent Antisymmetrized Geminal Power Theory Using a Coherent State Formulation.- Griffiths Inequalities for Fermion Systems.- Entropy of Reduced Density Matrices.- A Lower Bound to the Ground State Energy of a Boson System With Fermion Source.- Reduced Density Operators, Their Related von Neumann Density Operators, Close Cousins of These, and Their Physical Interpretation.- Theory and Practice of the Spin-Adapted Reduced Hamiltonians (SRH).- Variational Principle with Built-In Pure State N-Representability Conditions. The N-Electron Case.- Wigner Distributions as Representations of the Density Matrix.- Inter-Relationships Between Various Representations of One-Matrices and Related Densities: A Road Map and an Example ..- Current Problems in Density Functional Theory.- The Interface Between Reduced Density Matrices and Density Functional Theory.- The Physics Underlying the Langreth-Mehl Scheme for Non-Uniform Systems.- Understanding Energy Differences in Density Functional Theory.- Density Functional Calculations of Molecular Bond Energies.- Non-Local Effects on Atomic and Molecular Correlation Energies Studied with a Gradient-Corrected Density Functional.- An Evaluation of Local Electron Correlation Corrections and Non-Local Exchange Corrections to the Hartree-Fock-Slater Method from Calculations on Bond Energies and Electronic Spectra of Molecular Systems.- Correlation Energy Functionals of One-Matrices and Hartree-Fock Densities.- Some Remarks on Scaling Relations in Density Functional Theory.- Deduction of Semiempirical MO Methods from Density Functional Theory.- Charge and Spin Densities in Molecular Solids: Local Density Functional Calculations Versus Experiment.- A Functional of the Two-Particle Density Matrix for the Approximate Calculation of the Electronic Correlation Energy.- Extracules, Intracules, Correlation Holes, Potentials, Coefficients and All That.- The Exact Schroedinger Equation for the Electron Density.- Adiabatic Separation, Broken Symmetries and Geometry Optimization.- Asymptotic Results for Density Matrices and Electron Density in Atoms and Nearly Spherical Molecules.- An Algorithm for Calculating Isoelectronic Changes in Energies, Densities, and One-Matrices.- Atoms and Ions in the Limit of Large Nuclear Charge.- Improved Thomas-Fermi Theory for Atoms.- A Bond Energy from Quantum Mechanics.- Measured Electron Densities and Band Structure Calculations.- X-ray Orthonormal Orbital Model for Crystallography.

D. Reidel Pub. Co. eBooks•Jan 1, 1987•19 citations

n-Representability Problem for Reduced Density Matrices

Hans Kummer

In this paper we prove some theorems about the n-representability problem for reduced density operators. The first theorem (Theorem 6) sharpens a theorem proved by Garrod and Percus. Let Tnp be the set of all n-representable p-density operators. Then a density operator Dp belongs to Tnp¯ (the bar indicates the closure with respect to a certain topology) if and only if Tr (DpBp) ≥ 0 for all bounded self-adjoint p-particle operators Bp, such that their n-expansion (pn)ΓpnBp≡ ∑ i1<…<ipBp(i1…ip)is a positive operator in n-particle space. Moreover, it is shown that Tnp¯ is the closed convex hull of the exposed points of Tnp of finite one-rank (Theorem 9). A more practical version of this theorem may be formulated in the following manner (cf. Theorem 8). Consider the set γp of subspaces of the n-particle space, occurring as an eigenspace to the deepest eigenvalue of a bounded n-particle operator which is the n expansion of some p-particle operator. Choose from every element of γp one (and only one) vector (function) and form the corresponding reduced p-particle operator. Tnp¯ is the closed convex hull of all these p-density operators (cf. Theorem 9). For p = 1, this theorem reduces to Coleman's theorem about the n representability of the 1 matrix.

Journal of Mathematical Physics•Oct 1, 1967•207 citations
The NNN-representability problem is the problem of determining whether there exists NNN-particle states with some prescribed property. Here we report an affirmative solution to the fermion NNN-representability problem when both the density and the paramagnetic current density are prescribed. This problem arises in current-density functional theory and is a generalization of the well-studied corresponding problem (only the density prescribed) in density functional theory. Given any density and paramagnetic current density satisfying a minimal regularity condition (essentially that a von Weiz"acker--like canonical kinetic energy density is locally integrable), we prove that there exists a corresponding NNN-particle state. We prove this by constructing an explicit one-particle reduced density matrix in the form of a position-space kernel, i.e., a function of two continuous-position variables. In order to make minimal assumptions, we also address mathematical subtleties regarding the diagonal of, and how to rigorously extract paramagnetic current densities from, one-particle reduced density matrices in kernel form.
For a quantum-mechanical system of NNN identical fermions, the NNN-representability problem is the problem of recognizing whether, for a given pth-order reduced density matrix ${\ensuremath{\Gamma}}^{(p)}(12\ensuremath{\cdots}p|{1}^{\ensuremath{'}}{2}^{\ensuremath{'}}\ensuremath{\cdots}{p}^{\ensuremath{'}})$, there exists an antisymmetric NNN-particle wave function $\ensuremath{\Psi}(12\ensuremath{\cdots}N)$ such that ${\ensuremath{\Gamma}}^{(p)}(12\ensuremath{\cdots}p|{1}^{\ensuremath{'}}{2}^{\ensuremath{'}}\ensuremath{\cdots}{p}^{\ensuremath{'}})=({N}{p})\ensuremath{\int}\ensuremath{\Psi}(12\ensuremath{\cdots}N)\ifmmode\times\else\texttimes\fi{}{\ensuremath{\Psi}}^{*}({1}^{\ensuremath{'}}{2}^{\ensuremath{'}}\ensuremath{\cdots}{p}^{\ensuremath{'}},p+1\ensuremath{\cdots}N){d}_{{\ensuremath{\tau}}_{p+1}}\ensuremath{\cdots}{d}_{{\ensuremath{\tau}}_{N}}$. It is shown that if the Hamiltonian of a system is time-reversal invariant, and the number of particles, NNN, is even, the necessary and sufficient condition that an approximate first-order density matrix corresponding to a nondegenerate energy eigenstate be NNN-representable is that its natural spin-orbital occupation numbers be equal in pairs.
The problem of pure-state N-representability of the two-particle spin-dependent density function ρ(x1, x2) is considered for an N-electron system, and a procedure for finding an N-representable ρ(x1, x2) is advanced. The problem is formulated in the framework of a family of N × N matrices formed from integrals of auxiliary two-particle functions θn(x1, x2) converging at n → ∞ to ρ(x1, x2)/[N(N−1)]. The simple requirement of positive definiteness of these matrices is shown to play a decisive role in finding an N-representable ρ(x1, x2). The results obtained may open new possibilities for using ρ(x1, x2) in the density-functional theory. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 65: 127–142, 1997
The latest results of density matrix theory are summarized in terms of analysis of off-diagonal components. For a given reduced density matrix (RDM) of arbitrary order, the full density matrix is explicitly constructed, from which the RDM is derived by means of the corresponding reduction procedure. Both the matrices satisfy all the necessary conditions, such antisymmetry, normalization, hermiticity, except that the non-negative definiteness of the full DM is not guaranteed. If an RDM satisfies the necessary inequalities with the same number of indices as it has itself, then, as proved here, there is an additive correction to the full DM matrix that turns it into a non-negative definite matrix leading to the same RDM and satisfying the other requirements. This establishes a fundamental fact that these inequalities are sufficient as well, providing the N-representability of the given RDM. The construction is based on the classification, introduced in the article, of the multi-index elements of the RDM and full DM according to their degree of off-diagonality. RDM elements with a predetermined degree of off-diagonality are expressed through full DM elements with the same off-diagonality. This shows that the conventional variational approach to calculating the energy of the system by convex programming methods under the conditions of N-representability (before and even after reaching the minimum energy) generally does not lead to constructing a mixed quantum state composed of solutions to the Schrödinger equation. Therefore, except for the energy, it does not guarantee a correct description of the other properties so long as the conditions of correctness for each particular property are not expressed in terms of the RDM and not incorporated into the specific conditions of N-representability (regardless of the requirement of energy minimization).
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.
A Tribute to A. John Coleman - The Tame Mathematician.- Reduced Density Matrices: 1929-1989.- Some Aspects on the Development of the Theory of Reduced Density Matrices and the Representability Problem.- Representability Conditions.- On the Diagonal N-Representability Problem.- Fermion N-Representability Conditions Generated by a Decomposition of the 1-Particle Identity Operator onto Mutually Orthogonal Projection Operators.- The Unitarily Invariant Decomposition of Hermitian Operators.- Building Up N-Electron States with Symplectic Symmetry.- Time Dependent Antisymmetrized Geminal Power Theory Using a Coherent State Formulation.- Griffiths Inequalities for Fermion Systems.- Entropy of Reduced Density Matrices.- A Lower Bound to the Ground State Energy of a Boson System With Fermion Source.- Reduced Density Operators, Their Related von Neumann Density Operators, Close Cousins of These, and Their Physical Interpretation.- Theory and Practice of the Spin-Adapted Reduced Hamiltonians (SRH).- Variational Principle with Built-In Pure State N-Representability Conditions. The N-Electron Case.- Wigner Distributions as Representations of the Density Matrix.- Inter-Relationships Between Various Representations of One-Matrices and Related Densities: A Road Map and an Example ..- Current Problems in Density Functional Theory.- The Interface Between Reduced Density Matrices and Density Functional Theory.- The Physics Underlying the Langreth-Mehl Scheme for Non-Uniform Systems.- Understanding Energy Differences in Density Functional Theory.- Density Functional Calculations of Molecular Bond Energies.- Non-Local Effects on Atomic and Molecular Correlation Energies Studied with a Gradient-Corrected Density Functional.- An Evaluation of Local Electron Correlation Corrections and Non-Local Exchange Corrections to the Hartree-Fock-Slater Method from Calculations on Bond Energies and Electronic Spectra of Molecular Systems.- Correlation Energy Functionals of One-Matrices and Hartree-Fock Densities.- Some Remarks on Scaling Relations in Density Functional Theory.- Deduction of Semiempirical MO Methods from Density Functional Theory.- Charge and Spin Densities in Molecular Solids: Local Density Functional Calculations Versus Experiment.- A Functional of the Two-Particle Density Matrix for the Approximate Calculation of the Electronic Correlation Energy.- Extracules, Intracules, Correlation Holes, Potentials, Coefficients and All That.- The Exact Schroedinger Equation for the Electron Density.- Adiabatic Separation, Broken Symmetries and Geometry Optimization.- Asymptotic Results for Density Matrices and Electron Density in Atoms and Nearly Spherical Molecules.- An Algorithm for Calculating Isoelectronic Changes in Energies, Densities, and One-Matrices.- Atoms and Ions in the Limit of Large Nuclear Charge.- Improved Thomas-Fermi Theory for Atoms.- A Bond Energy from Quantum Mechanics.- Measured Electron Densities and Band Structure Calculations.- X-ray Orthonormal Orbital Model for Crystallography.
In this paper we prove some theorems about the n-representability problem for reduced density operators. The first theorem (Theorem 6) sharpens a theorem proved by Garrod and Percus. Let Tnp be the set of all n-representable p-density operators. Then a density operator Dp belongs to Tnp¯ (the bar indicates the closure with respect to a certain topology) if and only if Tr (DpBp) ≥ 0 for all bounded self-adjoint p-particle operators Bp, such that their n-expansion (pn)ΓpnBp≡ ∑ i1<…<ipBp(i1…ip)is a positive operator in n-particle space. Moreover, it is shown that Tnp¯ is the closed convex hull of the exposed points of Tnp of finite one-rank (Theorem 9). A more practical version of this theorem may be formulated in the following manner (cf. Theorem 8). Consider the set γp of subspaces of the n-particle space, occurring as an eigenspace to the deepest eigenvalue of a bounded n-particle operator which is the n expansion of some p-particle operator. Choose from every element of γp one (and only one) vector (function) and form the corresponding reduced p-particle operator. Tnp¯ is the closed convex hull of all these p-density operators (cf. Theorem 9). For p = 1, this theorem reduces to Coleman's theorem about the n representability of the 1 matrix.