© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

Russian Journal of Physical Chemistry B

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

M. M. Mestechkin

Published Mar 1, 2017• 0 citations

Abstract

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).

DOI: 10.1134/s1990793117020087

Related papers

On the size-consistency of the reduced-density-matrix method and the unitary invariant diagonal N-representability conditions

Maho Nakata, James S. M. Anderson

A promising variational approach for determining the ground state energy and its properties is by using the second-order reduced density matrix (2-RDM). However, the leading obstacle with this approach is the N-representability problem. By employing a subset of conditions (typically the P, Q, G, T1 and T2′ conditions) results comparable to those of CCSD(T) can be achieved. However, these conditions do not guarantee size-consistency. In this work, we show that size-consistency can be satisfied if the 2-RDM satisfies the following conditions: (i) the 2-RDM is unitary invariant diagonal N-representable; (ii) the 2-RDM corresponding to each (unspecified) subsystem is the eigenstate of the number of corresponding electrons; and (iii) the 2-RDM satisfies at least one of the P, Q, G, T1 and T2′ conditions. This is the first time that a computationally feasible (though demanding) sufficient condition for the RDM method that guarantees size-consistency in all chemical systems has been published in the literature.

arXiv:1108.5665•AIP Advances•Jul 24, 2012•10 citations

Uncertainty relations and reduced density matrices: Mapping many-body quantum mechanics onto four particles

David A. Mazziotti, Robert Erdahl

For the description of ground-state correlation phenomena an accurate mapping of many-body quantum mechanics onto four particles is developed. The energy for a quantum system with no more than two-particle interactions may be expressed in terms of a two-particle reduced density matrix (2-RDM), but variational optimization of the 2-RDM requires that it corresponds to an N-particle wave function. We derive N-representability conditions on the 2-RDM that guarantee the validity of the uncertainty relations for all operators with two-particle interactions. One of these conditions is shown to be necessary and sufficient to make the RDM solutions of the dispersion condition equivalent to those from the contracted Schr\"odinger equation (CSE) [Mazziotti, Phys. Rev. A 57, 4219 (1998)]. In general, the CSE is a stronger N-representability condition than the dispersion condition because the CSE implies the dispersion condition as well as additional N-representability constraints from the Hellmann-Feynman theorem. Energy minimization subject to the representability constraints is performed for a boson model with 10, 30, and 75 particles. Even when traditional wave-function methods fail at large perturbations, the present method yields correlation energies within 2%.

Physical Review A•Mar 21, 2001•187 citations

Quantum Many-Body Theory from a Solution of the N -Representability Problem

David A. Mazziotti

Here we present a many-body theory based on a solution of the N-representability problem in which the ground-state two-particle reduced density matrix (2-RDM) is determined directly without the many-particle wave function. We derive an equation that re-expresses physical constraints on higher-order RDMs to generate direct constraints on the 2-RDM, which are required for its derivation from an N-particle density matrix, known as N-representability conditions. The approach produces a complete hierarchy of 2-RDM constraints that do not depend explicitly upon the higher RDMs or the wave function. By using the two-particle part of a unitary decomposition of higher-order constraint matrices, we can solve the energy minimization by semidefinite programming in a form where the low-rank structure of these matrices can be potentially exploited. We illustrate by computing the ground-state electronic energy and properties of the H_{8} ring.

arXiv:2304.08570•Physical Review Letters•Apr 12, 2023•41 citations

Geometry of density matrices. II. Reduced density matrices and N representability

John E. Harriman

The reduction of density matrices defined with respect to a finite basis set is considered. A basis set can be introduced into the vector space of Hermitian matrices acting on functions of the coordinates of $p$ particles for each value of $p$. An analogy with the construction of spin eigenfunctions is used to obtain a particular basis set called the reducing basis. Reduction of one of these basis elements is either one to one or maps into the origin. Any element in the preimage of a density matrix can thus be resolved into two components, one uniquely determined by the density matrix and the other arbitrary within a certin subset. A restatement of the $N$-representability problem is given, and two sufficient conditions and one necessary condition for $N$ representability are given in terms of distances and norms.

Physical Review A•Apr 1, 1978•84 citations

Construction of Approximately N -Representable Density Matrices

Jack Simons, John E. Harriman

A method is given for obtaining the density matrix or matrices, expressible in terms of a given set of spin geminals, which are as nearly $N$ representable as possible. A measure of $N$ representability is introduced as the norm of the antisymmetric component of a normalized wave function leading to the density matrix. An error estimate is given for the maximum deviation of an energy calculated with the density matrix from that calculated with an antisymmetric wave function. Linear variational parameters may occur in the density matrix, and their number can be increased at the expense of $N$ representability. Implementation of the method requires diagonalization of a large matrix, and the consequences of truncation are considered. The method is also related to the problem of the exact $N$ representability of a given density matrix.

Physical Review A•Sep 1, 1970•27 citations

Integral–geometrical consideration of density matrices

Boris P. Zapol

Abstract The ensemble N‐representability problem for the k‐th order reduced density matrix (k‐RDM) as well as the problem of reconstruction of the N‐particle system density matrices (N‐DM) from a given k‐RDM are studied. The spatial parts of the k‐RDM expansion in terms of spin tensorial operators Θ are represented using particular values (at specially chosen ) of the Radon transform of the N‐DM spatial parts (or their sums) 𝒟Nλ(x′ | x″) (here, is a d‐plane in the n‐space ℝn of x = (x′, x″)), with n = 6N, d = 3 (N − k), x′ ≡ (r′1, ⃛, r′N), x′ ≡ (r1″, ⃛, r″N ()). In this way, the problem is reduced to investigation of the properties of the functions . For a normalizable N – DM, it is proved that are bounded functions. The properties of implied by the N‐DM permutational symmetry, Hermiticity, and positive definiteness are found. A formal procedure of reconstruction of all N‐DM corresponding to a given k‐RDM is proposed. © 1995 John Wiley & Sons, Inc.

International Journal of Quantum Chemistry•Dec 5, 1995•2 citations

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).

A promising variational approach for determining the ground state energy and its properties is by using the second-order reduced density matrix (2-RDM). However, the leading obstacle with this approach is the N-representability problem. By employing a subset of conditions (typically the P, Q, G, T1 and T2′ conditions) results comparable to those of CCSD(T) can be achieved. However, these conditions do not guarantee size-consistency. In this work, we show that size-consistency can be satisfied if the 2-RDM satisfies the following conditions: (i) the 2-RDM is unitary invariant diagonal N-representable; (ii) the 2-RDM corresponding to each (unspecified) subsystem is the eigenstate of the number of corresponding electrons; and (iii) the 2-RDM satisfies at least one of the P, Q, G, T1 and T2′ conditions. This is the first time that a computationally feasible (though demanding) sufficient condition for the RDM method that guarantees size-consistency in all chemical systems has been published in the literature.

For the description of ground-state correlation phenomena an accurate mapping of many-body quantum mechanics onto four particles is developed. The energy for a quantum system with no more than two-particle interactions may be expressed in terms of a two-particle reduced density matrix (2-RDM), but variational optimization of the 2-RDM requires that it corresponds to an N-particle wave function. We derive N-representability conditions on the 2-RDM that guarantee the validity of the uncertainty relations for all operators with two-particle interactions. One of these conditions is shown to be necessary and sufficient to make the RDM solutions of the dispersion condition equivalent to those from the contracted Schr\"odinger equation (CSE) [Mazziotti, Phys. Rev. A 57, 4219 (1998)]. In general, the CSE is a stronger N-representability condition than the dispersion condition because the CSE implies the dispersion condition as well as additional N-representability constraints from the Hellmann-Feynman theorem. Energy minimization subject to the representability constraints is performed for a boson model with 10, 30, and 75 particles. Even when traditional wave-function methods fail at large perturbations, the present method yields correlation energies within 2%.

Here we present a many-body theory based on a solution of the N-representability problem in which the ground-state two-particle reduced density matrix (2-RDM) is determined directly without the many-particle wave function. We derive an equation that re-expresses physical constraints on higher-order RDMs to generate direct constraints on the 2-RDM, which are required for its derivation from an N-particle density matrix, known as N-representability conditions. The approach produces a complete hierarchy of 2-RDM constraints that do not depend explicitly upon the higher RDMs or the wave function. By using the two-particle part of a unitary decomposition of higher-order constraint matrices, we can solve the energy minimization by semidefinite programming in a form where the low-rank structure of these matrices can be potentially exploited. We illustrate by computing the ground-state electronic energy and properties of the H_{8} ring.

The reduction of density matrices defined with respect to a finite basis set is considered. A basis set can be introduced into the vector space of Hermitian matrices acting on functions of the coordinates of

p

particles for each value of

p

. An analogy with the construction of spin eigenfunctions is used to obtain a particular basis set called the reducing basis. Reduction of one of these basis elements is either one to one or maps into the origin. Any element in the preimage of a density matrix can thus be resolved into two components, one uniquely determined by the density matrix and the other arbitrary within a certin subset. A restatement of the

N

-representability problem is given, and two sufficient conditions and one necessary condition for

N

representability are given in terms of distances and norms.

A method is given for obtaining the density matrix or matrices, expressible in terms of a given set of spin geminals, which are as nearly

N

representable as possible. A measure of

N

representability is introduced as the norm of the antisymmetric component of a normalized wave function leading to the density matrix. An error estimate is given for the maximum deviation of an energy calculated with the density matrix from that calculated with an antisymmetric wave function. Linear variational parameters may occur in the density matrix, and their number can be increased at the expense of

N

representability. Implementation of the method requires diagonalization of a large matrix, and the consequences of truncation are considered. The method is also related to the problem of the exact

N

representability of a given density matrix.

Abstract The ensemble N‐representability problem for the k‐th order reduced density matrix (k‐RDM) as well as the problem of reconstruction of the N‐particle system density matrices (N‐DM) from a given k‐RDM are studied. The spatial parts of the k‐RDM expansion in terms of spin tensorial operators Θ are represented using particular values (at specially chosen ) of the Radon transform of the N‐DM spatial parts (or their sums) 𝒟Nλ(x′ | x″) (here, is a d‐plane in the n‐space ℝn of x = (x′, x″)), with n = 6N, d = 3 (N − k), x′ ≡ (r′1, ⃛, r′N), x′ ≡ (r1″, ⃛, r″N ()). In this way, the problem is reduced to investigation of the properties of the functions . For a normalizable N – DM, it is proved that are bounded functions. The properties of implied by the N‐DM permutational symmetry, Hermiticity, and positive definiteness are found. A formal procedure of reconstruction of all N‐DM corresponding to a given k‐RDM is proposed. © 1995 John Wiley & Sons, Inc.