Boundedness is demonstrated for the fifth derivative of the one-particle reduced density matrix for non-relativistic Coulombic wavefunctions in the vicinity of the diagonal. To prove this result, improved pointwise bounds are obtained for cluster derivatives of wavefunctions involving multiple clusters.
We obtain bounds for all derivatives of the non-relativistic Coulombic one-particle density matrix γ(x,y) near the diagonal x=y.
An explicit expression relating the magnitude of the fifth-order off-diagonal cusp in the real part of the one-electron reduced density matrix to the “on-top” two-electron density is derived in a rigorous manner from the behavior of the underlying electronic wavefunction at the electron–electron coalescence. The implications of the presence of this cusp upon electronic structure calculations of quantum chemistry and solid-state physics, including the limits imposed upon their accuracy, are elucidated. In particular, the power-law decay of the occupation numbers of the natural orbitals is demonstrated for 1S states of systems composed of arbitrary even numbers of electrons. The practical importance of analogous off-diagonal cusps in many-electron reduced density matrices is briefly discussed.
Abstract A partial trace over the occupation numbers of all but k states in the density matrix of an ensemble with an arbitrary number of single‐particle states is defined as the (reduced) k‐state density matrix. This matrix is used to obtain a complete, practical solution to the problem of determining the representability of the diagonal elements of the one‐ and two‐particle (reduced) density matrices. This solution is expressed as a series of linear inequalities involving the density‐matrix elements; the inequalities are identical with those derived previously by Davidson and McCrae by a different method. In addition, our method is used to obtain nonlinear, matrix inequalities on the off‐diagonal elements of the density matrices.
We tabulate, on the diagonal, the nth power of the primitive approximation for the pair Coulomb density matrix on a sphere, with up to n=5, and we discuss its behavior at particles contact as n grows. Our results show that even on a curved surface quantum statistical physics is fundamentally different from its classical counterpart. If in the latter one needs to artificially introduce a regularization of the divergent attractive Coulomb potential preventing particles contact, for example with a hard core, in the former this is not necessary.