One-body reduced density matrix 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!
 A: 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.
$$
