Particle-reduced Hamiltonians

After the energy is expressed as a functional of the 2-RDM, a systematic hierarchy of V-representabihty constraints, known as p-positivity conditions, is derived [17]. We develop the details of the 2-positivity, 3-positivity, and partial 3-positivity conditions [21, 27, 34, 33]. In Section II.E the formal solution of V-representability for the 2-RDM is presented through a convex set of two-particle reduced Hamiltonian matrices [7, 21]. It is shown that the positivity conditions correspond to certain classes of reduced Hamiltonian matrices, and consequently, they are exact for certain classes of Hamiltonian operators at any interaction strength. In Section II.F the size of the 2-RDM is reduced through the use of spin and spatial symmetries [32, 34], and in Section II.G the variational 2-RDM method is extended to open-shell molecules [35]. 

The formal solution of Al-representability for the 2-RDM is developed in terms of a convex set of two-particle reduced Hamiltonian matrices. To complement the derivation of the positivity conditions from the metric matrices, we derive them from classes of these two-particle reduced Hamiltonian matrices. This interpretation allows us to demonstrate that the 2-positivity conditions are exact for certain classes of Hamiltonian operators for any interaction strength. In this section all of the ROMs are normalized to unity. Much of this discussion appeared originally in Refs. [21, 29]. 

Direct minimization of the energy as a functional of the p-RDM may be achieved if the p-particle density matrix is restricted to the set of Al-represen-table p-matrices, that is, p-matrices that derive from the contraction of at least one A-particle density matrix. The collection of ensemble Al-representable p-RDMs forms a convex set, which we denote as P. To define P, we first consider the convex set of p-particle reduced Hamiltonians, which are... 

When N = p, the set B simply contains the p-particle reduced Hamiltonians, which are positive semidefinite, but when N = p + 1, because the lifting process raises the lowest eigenvalue of the reduced Hamiltonian, the set also contains p-particle reduced Hamiltonians that are lifted to positive semidefinite matrices. Consequently, the number of Wrepresentability constraints must increase with N, that is, B C B. To constrain the p-RDMs, we do not actually need to consider all pB in B, but only the members of the convex set B, which are extreme A member of a convex set is extreme if and only if it cannot be expressed as a positively weighted ensemble of other members of the set (i.e., the extreme points of a square are the four corners while every point on the boundary of a circle is extreme). These extreme constraints form a necessary and sufficient set of A-representability conditions for the p-RDM [18, 41, 42], which we can formally express as... 

A significant class of p-particle reduced Hamiltonians in the set B includes... 

A quantum system of N particles may also be interpreted as a system of (r — N) holes, where r is the rank of the one-particle basis set. The complementary nature of these two perspectives is known as the particle-hole duality [13, 44, 45]. Even though we treated only the iV-representability for the particles in the formal solution, any p-hole RDM must also be derivable from an (r — A)-hole density matrix. While the development of the formal solution in the literature only considers the particle reduced Hamiltonian, both the particle and the hole representations for the reduced Hamiltonian are critical in the practical solution of N-representability problem for the 1-RDM [6, 7]. The hole definitions for the sets and are analogous to the definitions for particles except that the number (r — N) of holes is substituted for the number of particles. In defining the hole RDMs, we assume that the rank r of the one-particle basis set is finite, which is reasonable for practical calculations, but the case of infinite r may be considered through the limiting process as r —> oo. 

Any arbitrary one-particle reduced Hamiltonian shifted by its A-particle ground-state energy must be expressible by the extreme HamUtonian elements in the convex set As we showed in Eq. (52), keeping the 1-RDM positive semide-finite is equivalent to applying the Al-representability constraints in Eq. (50) for the class of extreme positive semidefinite which may be parameterized by... 

The ground-state energies of atoms and molecules where the A-particle Hamiltonian is defined by Eq. (48) may be expressed through three different representations of the 2-RDM and the two-particle reduced Hamiltonian ... 

The three complementary representations of the reduced Hamiltonian offer a framework for understanding the D-, the Q-, and the G-positivity conditions for the 2-RDM. Each positivity condition, like the conditions in the one-particle case, correspond to including a different class of two-particle reduced Hamiltonians in the A-representability constraints of Eq. (50). The positivity of arises from employing all positive semidefinite in Eq. (50) while the Q- and the G-conditions arise from positive semidefinite and B, respectively. To understand these positivity conditions in the particle (or D-matrix) representation, we define the D-form of the reduced Hamiltonian in terms of the Q- and the G-representations ... 

The Q- and the G-conditions are thus equivalent to the constraints in Eq. (50) with the two-particle reduced Hamiltonians in Eqs. (68) and (69), where B >0 and B > 0. Unlike the one-particle case, these reduced Hamiltonians do not exhaust all of the extreme constraints in Eq. (50), and yet the explicit forms of the Hamiltonians give us insight into the variety of correlated Hamiltonians that can be treated accurately. 

Many methods in chemistry for the correlation energy are based on a form of perturbation theory, but the positivity conditions are quite different. Traditional perturbation theory performs accurately for all kinds of two-particle reduced Hamiltonians, which are close enough to a mean-field (Hartree-Fock) reference. There are a myriad of chemical systems, however, where the correlated wave-function (or 2-RDM) is not sufficiently close to a statistical mean field. Different from perturbation theory, the positivity conditions function by increasing the number of extreme two-particle Hamiltonians in which are employed as constraints upon the 2-RDM in Eq. (50) and, hence, they exactly treat a certain convex set of reduced Hamiltonians to all orders of perturbation theory. For the... 

D-, the Q-, and the G-conditions we have If the two-particle reduced Hamiltonian shifted by its N-particle ground-state energy can be written as an ensemble of the reduced Hamiltonians in the set > 0 as well as the Q- and the G-reduced Hamiltonians parameterized in Eqs. (68) and (69), then the energy for an N-particle system may be computed exactly. 

Nakatsuji [37] in 1976 first proved that with the assumption of N-representability [3] a 2-RDM and a 4-RDM will satisfy the CSE if and only if they correspond to an A-particle wavefunction that satishes the corresponding Schrodinger equation. Just as the Schrodinger equation describes the relationship between the iV-particle Hamiltonian and its wavefunction (or density matrix D), the CSE connects the two-particle reduced Hamiltonian and the 2-RDM. However, because the CSE depends on not only the 2-RDM but also the 3- and 4-RDMs, it cannot be solved for the 2-RDM without additional constraints. Two additional types of constraints are required (i) formulas for building the 3- and 4-RDMs from the 2-RDM by a process known as reconstruction, and (ii) constraints on the A-representability of the 2-RDM, which are applied in a process known as purification. 

As in the previous section, by connected we mean all terms that scale linearly with N. Wedge products of cumulant RDMs can scale linearly if and only if they are connected by the indices of a matrix that scales linearly with N transvec-tion). In the previous section we only considered the indices of the one-particle identity matrix in the contraction (or number) operator. In the CSE we have the two-particle reduced Hamiltonian matrix, which is defined in Eqs. (2) and (3). Even though the one-electron part of scales as N, the division by A — 1 in Eq. (3) causes it to scale linearly with N. Hence, from our definition of connected, which only requires the matrix to scale linearly with N, the transvection... 

The unitary decomposition may be applied to any Hermitian, antisymmetric two-particle matrix including the 2-RDM, the two-hole RDM, and the two-particle reduced Hamiltonian. The decomposition is also readily generalized to treat p-particle matrices [80-82]. The trial 2-RDM to be purified may be written... 

Recall that in his Theorems 3 and 4 Hans Kummer [3] defined a contraction operator, L, which maps a linear operator on A-space onto an operator on p-space and an expansion operator, E, which maps an operator on p-space onto an operator on A-space. Note that the contraction and expansion operators are super operators in the sense that they act not on spaces of wavefunctions but on linear spaces consisting of linear operators on wavefunction spaces. If the two-particle reduced Hamiltonian is defined as... 

---