N-representability of the Jastrow wave function pair density of the lowest-order
Conditions for the N-representability of the pair density (PD) are needed for the development of the PD functional theory. We derive sufficient conditions for the N-representability of the PD that is calculated from the Jastrow wave function within the lowest order. These conditions are used as the constraints on the correlation function of the Jastrow wave function. A concrete procedure to search the suitable correlation function is also presented.
<!>Results<!>Note that the normalization constant<!>Discussions<!>Concluding remarks
<p>The pair density (PD), which is defined as the diagonal element of the second-order reduced density matrix, has more information about the electron correlation than the electron density 1-4 . Therefore, the PD functional theory has been expected to be one of the promising schemes beyond the density functional theory 5-28 . We have proposed the PD functional theory that yields the best PD within the set of the Jastrow wave function PDs of the lowest-order (LO-Jastrow PDs) 17 . The search region for the ground-state PD is substantially extended as compared with the previous theory 18, 19 . On the other hand, however, there remains a significant problem related to the N-representability of the LO-Jastrow PDs 17 . Concerning the N-representability of the PD or the second-order reduced density matrix, there are a lot of works 29-39 . But unfortunately, the necessary and sufficient conditions for the N-representability of the PD are not yet known in a practical form.</p><p>Let us revisit the problem here. We shall consider an N 0 -electron system. The Jastrow wave function is given by</p><p>where x i denotes the coordinates including the spatial coordinate r i and spin coordinate η i , and where</p><p>SSD N 1 0 are the normalization constant, correlation function and single Slater determinant (SSD), respectively. The LO-Jastrow PD is given by 17, 40, 41</p><p>0 where γ ′ ′ = rr rr ( ; )</p><p>0 is the PD calculated from the SSD. In the preceding paper 17 , we have confirmed that Eq. ( 2) meets four kinds of necessary conditions for the N-representability of the PD, and may become "approximately N-representable" 42-47 . However, the possibility of it being N-representable has not been discussed 17 . This is an arguable problem that is concerned with whether the reproduced PD is physically reasonable or not.</p><p>The aim of this paper is to discuss the N-representability of Eq. ( 2) and to show the way to search the suitable correlation function</p><p>The organization of this paper is as follows. For the convenience of the subsequent discussions, we first examine the properties of the LO-Jastrow PD in the next section. Then, the sufficient conditions for the N-representability of Eq. ( 2), which are imposed on the correlation function, will be derived recursively. Next, concrete steps for searching the correlation function that meets these conditions are discussed. Finally, concluding remarks are given in the last section.</p><!><p>Properties of the LO-Jastrow PD. In this section, we shall discuss the properties of the LO-Jastrow PD. To this aim, the properties of PDs that are calculated from SSDs are investigated. The cofactor expansion of</p><p>In what follows, suppose that these spin orbitals are given as the solutions of simultaneous equations of previous work 17 , and therefore they are orthonormal to each other. 3) denote (N 0 − 1) -electron SSDs that are defined as the minor determinants multiplied by</p><p>. The PD that is calculated from Eq. ( 3) is given by</p><p>0 can be expressed by the sum of γ ′ ′ = − rr rr ( ; )</p><p>. By the repetition of this procedure, we arrive at</p><p>are PDs calculated from the two-electron SSDs. These two-electron SSDs, which are denoted by Φ ⋅⋅⋅ x x ( , )</p><p>, are obtained by the above-mentioned successive cofactor expansions.</p><p>Multiplying both sides of Eq. ( 5) by</p><p>This relation is the starting point to examine the N-representability of the LO-Jastrow PD.</p><p>Sufficient conditions for the N-representablity of the LO-Jastrow PD. We shall start with consid-</p><p>2 that appears in the right-hand side of Eq. ( 6). Suppose that the two-electron wave function that yields</p><p>2 is given by the following Jastrow wave function;</p><p>where</p><p>is the normalization constant. From Eq. ( 7), the PD is calculated as</p><p>. Therefore, we get</p><p>2 . Hereafter, we assume that such conditions hold for all values of ⋅ ⋅ ⋅ i j p , , , and q, i.e.,</p><p>ij pq 2 , and check whether the conditions Eq. ( 8) are satisfied or not. As mentioned later, the conditions Eq. ( 8) are the parts of sufficient conditions for the N-representability of Eq. ( 2).</p><!><p>Under the conditions Eq. ( 8), we have</p><p>where γ ′ ′ = ˆrr rr ( ;</p><p>denotes the PD operator for an n-electron system. Substituting Eq. ( 9) into Eq. ( 6) and rearranging, we get</p><p>It should be noticed that the right-hand side of Eq. ( 10) has a characteristic form. Concerning the N-representability of this form, the following theorem holds: Theorem. If there exists the set of single-valued, continuous, smooth and finite functions {a α (x n + 1 )} (1 ≤ α ≤ n + 1) that satisfy the conditions;</p><p>then the following equations hold:</p><p>Here</p><p>(1 ≤ α ≤ n + 1) denote the n-electron wave functions, and σ is a permutation operator upon the electron coordinates, and σ is the number of interchanges in σ.</p><p>Proof. The left-hand side of Eq. ( 13) seems to be related to the average of γ ′ ′ = ˆrr rr ( ;</p><p>with respect to a density matrix for a mixed state. Indeed, if the density matrix for the mixed state is given by</p><p>is calculated as</p><p>ˆˆn rr rr rr rr Tr ( ; ) 1 1 ( ; ) ( 16)</p><p>(2)</p><p>By using Eq. ( 16), the left-hand side of Eq. ( 13) is rewritten as</p><p>. On the other hand, it is expected that the average</p><p>may be given as the expectation value of γ ′ ′ = ˆrr rr ( ;</p><p>with respect to a pure state for the whole system that includes the n-electron system as a subsystem 48, 49 . We shall take an (n + 1)-electron system as the whole system, and suppose that the wave function for such the (n + 1)-electron system is given by Eq. ( 14) with Eq. (11). Then we indeed get</p><p>1 is antisymmetric, i.e., if Eq. ( 12) holds, then the expectation value of</p><p>Consequently, Eq. ( 13) immediately follows from Eqs ( 16), ( 17) and (18). This means that Eq. ( 13) holds under the conditions that there exists the set of functions {a α (x n + 1 )} (1 ≤ α ≤ n + 1) that satisfy Eqs (11) and (12). Q.E.D. First, we apply the above theorem to the term</p><p>that appears in Eq. (10). According to the theorem, this term can be rewritten as</p><p>if there exists the set of functions</p><p>3 (1 ≤ q ≤ 3) that satisfy the following conditions:</p><p>q q ijk op N ijk opq q q ijk op N ijk opq</p><p>In addition to Eq. ( 8), the existence conditions for ⋅⋅⋅ a x { ( )} q ijk op 3 (1 ≤ q ≤ 3), i.e., Eqs ( 21) and (22), are also the parts of sufficient conditions for the N-representability of Eq. ( 2). We assume that the set of functions ⋅⋅⋅ { } a x ( ) q ijk op 3</p><p>(1 ≤ q ≤ 3) is obtained. Substitution of Eq. ( 20) into Eq. ( 10) leads to</p><p>that appears in Eq. ( 23). Similarly to the above, we obtain  The existence conditions for</p><p>4 (1 ≤ p ≤ 4), i.e., Eqs ( 26) and (27), are added to the set of sufficient conditions for the N-representability of Eq. (2). At this stage, Eq. ( 8) and the existence conditions for</p><p>4 (1 ≤ p ≤ 4) belong to the set of sufficient conditions. Thus the theorem is applied repeatedly, so that further conditions are added to the set of sufficient conditions. Continuing until the conditions for a x { ( )} i N 0 (1 ≤ i ≤ N 0 ), we eventually obtain sufficient conditions for the N-representability of Eq. ( 2).</p><!><p>Concrete steps for constructing the N-representable LO-Jastrow PD. In the preceding section, the sufficient conditions for the N-representability of the LO-Jastrow PD are derived. In this section, we consider the concrete steps for searching the correlation function that meets these conditions or checking its existence.</p><p>1. First, we give a trial form of the correlation function. Using this, simultaneous equations for the N 0 -electron system are solved in a self-consistent way 17 . 2. Let us consider the SSD that consists of the resultant spin orbitals for the simultaneous equations. The SSD can generally be expanded using the cofactor. The SSD for the N 0 -electron system is expanded along the N 0 th row, then we get the N 0 number of SSDs for the (N 0 − 1)-electron system, i.e., Φ ⋅ ⋅ ⋅ − x x ( , , )</p><p>. Successively, each of the SSDs for the (N 0 − 1)-electron system is expanded along the (N 0 − 1)th row, and then the (N 0 − 1) number of SSDs for the (N 0 − 2)-electron system, i.e., Φ ⋅ ⋅ ⋅</p><p>, can be obtained for each i. After that and later, the cofactor expansions are likewise repeated, and we finally arrive at the SSDs for the two-electron system, i.e.,Φ ⋅⋅⋅⋅ x x ( , )</p><p>, , and q are integers such that 1 for all cases of ⋅ ⋅ ⋅ ⋅ i j k o p , , , , , and q through Eq. ( 7). 4. By means of these</p><p>2 , three-electron wave functions are defined as Eq. ( 20). We will check whether there exists the set of</p><p>that are satisfied with Eqs. ( 21) and (22). Note that the check has to be performed for all cases of</p><p>3 cannot be found, we again return to the step 1 and modify the correlation function. Suppose the suitable</p><p>is found in this process, then three-electron antisymmetric wave functions Ψ = ⋅⋅⋅⋅ x x x ( , , )</p><p>3 for all cases of … i j k o , , , , and p can be constructed from Eq. (20). 5. We successively proceed the case of four-electron wave functions that are defined as Eq. (25). In a similar way to the step 4, we check whether the set of</p><p>meets the conditions of Eqs. ( 26) and (27). If no, we restart from the step 1 with the modified correlation function. Suppose that the suitable 4 can be obtained from Eq. ( 25). 6. Likewise, we further proceed the problem constructing the antisymmetric wave functions for more-electron systems. We search the set of</p><p>that are satisfied with Eqs. ( 11) and ( 12), together with modifying the correlation function. If we successfully find the correlation function that meets the conditions (11) and (12) for any n(≤N 0 − 1), the LO-Jastrow PD of the N 0 -electron system becomes N-representable.</p><p>As easily inferred, the above steps are feasible only for the small-electron systems from the practical viewpoint. However, it should be noted that the correlation function that makes the LO-Jastrow PD N-representable may, in principle, be found along the above steps, though there is a possibility that the suitable correlation function may not exist in some system.</p><p>Other possibilities to obtain antisymmetric wave functions. In the above-mentioned concrete steps, antisymmetric (n + 1)-electrons wave functions are built up from given antisymmetric n-electrons wave functions via Eq. ( 14). In this subsection, we show that there are other possibilities to obtain antisymmetric (n+1)-electrons wave functions without using Eq. ( 14).</p><p>Instead of Eq. ( 14), we can use the following expression for</p><p>where a α (x i ) and b α (x n + 1 ) denotes functions that should be determined. If these functions satisfy the following conditions;</p><p>1 becomes an antisymmetric wave function and yields the PD just given by the left-hand side of Eq. (13). Therefore, if we use the expression Eq. ( 28), then we have to search a set of functions a α (x i ) and b α (x n + 1 ) that satisfy Eqs. ( 29)-(31) instead of Eqs. (11) and (12) in the above-mentioned concrete steps.</p><p>Similarly to the above, we can choose the other expression for Ψ</p><p>1 instead of Eq. ( 14) or (28). For example, Ψ</p><p>SSD n 1 that are constructed by using given spin orbitals φ</p><p>where a α (x i , x n + 1 ) denote a function that should be determined. If a α (x i , x n + 1 ) satisfies the following conditions;  then Eq. ( 32) becomes an antisymmetric (n + 1)-electrons wave function and yields the PD just given by the left-hand side of Eq. (13). Therefore, if we use the expression Eq. (32), then what we have to do in the concrete steps is to search a set of functions a α (x i , x n+1 ) that satisfy Eqs (33) and (34) instead of searching a set of functions that satisfy Eqs (11) and (12). Thus, various expressions for Ψ ⋅ ⋅ ⋅ = + + x x ( , , )</p><p>N n n 1 1</p><p>1 can be adopted. This means that the present method can provide various prescriptions to construct the N-representable LO-Jastrow PD.</p><!><p>In this paper, the sufficient conditions for the N-representability of the LO-Jastrow PD are discussed. Using the properties of the LO-Jastrow PD, we derive the sufficient conditions that are imposed on the correlation function of the Jastrow wave function. As shown in the previous section, additional steps to search the suitable correlation function, which satisfies the sufficient conditions, are attached to the computational scheme proposed previously 17 . Although the number of steps rapidly increases with that of electrons, the concrete steps that are presented in the previous section are feasible for a small-electron system. Of course, there is a possibility that the suitable correlation function cannot be found out. In this case, as mention in the previous section, we can adopt other expressions for Ψ ⋅ ⋅ ⋅ = + + x x ( , , )</p><p>N n n 1 1</p><p>1 , so that we may possibly find out a suitable correlation function. Otherwise, as mentioned in the previous paper, LO-Jastrow PDs are approximately N-representable in a sense that they satisfy four kinds of necessary conditions 17, 42-47 .</p>