@ACuriousMind Just want to double confirm something here. Is what is happening in this formula
$$\gamma (\vec{x}_1,\vec{x}_2;\vec{x}_1',\vec{x}_2')=N_e(N_e-1)\int \Psi^*(\vec{x}_1',\vec{x}_2',\dots ,\vec{x}_{N_e}) \Psi(\vec{x},\vec{x}_2,\dots ,\vec{x}_{N_e})d\vec{x}_3\cdots d\vec{x}_{N_e}$$ we are basically taking the partial trace of $N_e-2$ of the $N_e$ electrons in the $N_e$ electron density matrix?