One-body reduced density matrix

Question

Assume there is a N-particle state denoted as $|\Psi_N\rangle$, the density operator from its definition reads $\gamma_N = |\Psi_N\rangle\langle\Psi_N|$ and the density matrix elements take the form like \begin{equation} \langle\alpha_1,\alpha_2,\cdots,\alpha_N|\gamma_N|\beta_1,\beta_2,\cdots,\beta_N\rangle. \end{equation} The one-body reduced density matrix is defined as \begin{equation} \gamma_1(\alpha_1,\beta_1)=N\int\cdots\int \langle\alpha_1,\alpha_2,\cdots,\alpha_N|\gamma_N|\beta_1,\alpha_2,\cdots,\alpha_N\rangle d\alpha_2\cdots d\alpha_N. \end{equation} However, in some other books, someone writes down the one-body reduced density matrix using second quantization like \begin{equation} \gamma_1(\alpha_1,\beta_1) = \langle\Psi_N| \mathcal{c}^\dagger_{\alpha_1}\mathcal{c}_{\beta_1}|\Psi_N\rangle. \end{equation} So my question is do these two definitions of the one-body reduced density matrix equivalent? If they are, how to prove? Thank you!

Answer 1

First, $\alpha_1$ and $\beta_1$ in your second-quantized definition of the one-body density matrix should be interchanged (we will see below this is correct): $$ \gamma_1(\alpha_1,\beta_1)=\langle\Psi_N|c_{\beta_1}^+c_{\alpha_1}|\Psi_N\rangle $$ Then let's insert the unit operator $I=\sum_i|i\rangle\langle i|$, where $|i\rangle$ are all many-particle states, between $c_{\beta_1}^+$ and $c_{\alpha_1}$: $$ \gamma_1(\alpha_1,\beta_1)=\sum_i\langle\Psi_N|c_{\beta_1}^+|i\rangle\langle i|c_{\alpha_1}|\Psi_N\rangle. $$ Only the $(N-1)$-particle states $|i\rangle=|\alpha_2\ldots\alpha_N\rangle$ will survive in this sum: $$ \gamma_1(\alpha_1,\beta_1)=\int d\alpha_2\ldots d\alpha_N\langle\Psi_N|c_{\beta_1}^+|\alpha_2\ldots\alpha_N\rangle\langle \alpha_2\ldots\alpha_N|c_{\alpha_1}|\Psi_N\rangle=\\ =\int d\alpha_2\ldots d\alpha_N\langle \alpha_2\ldots\alpha_N|c_{\alpha_1}\gamma_Nc_{\beta_1}^+|\alpha_2\ldots\alpha_N\rangle. $$ Then, using the property of creation operators $$ c_{\beta_1}^+|\alpha_2\ldots\alpha_N\rangle=\sqrt{N}|\beta_1\alpha_2\ldots\alpha_N\rangle,\qquad\langle\alpha_2\ldots\alpha_N|c_{\alpha_1}=\sqrt{N}\langle\alpha_1\alpha_2\ldots\alpha_N|, $$ we get the first-quantized definition: $$ \gamma_1(\alpha_1,\beta_1)=N\int d\alpha_2\ldots d\alpha_N\langle\alpha_1\alpha_2\ldots\alpha_N|\gamma_N|\beta_1\alpha_2\ldots\alpha_N\rangle. $$