Si tengo un sistema de $N$ fermiones no interactuantes, puedo escribir la función de onda del estado fundamental del sistema utilizando un determinante de Slater
$$ \Phi_{0}(\textbf{r}_{1}, ..., \textbf{r}_{N}) = \frac{1}{\sqrt{N!}} \begin{vmatrix} \phi_{k_{1}}(\textbf{r}_{1}) & \cdots & \phi_{k_{N}}(\textbf{r}_{1})\\ \vdots & \ddots & \vdots \\ \phi_{k_{1}}(\textbf{r}_{N}) & \cdots & \phi_{k_{N}}(\textbf{r}_{N}) \end{vmatrix}, $$
donde $\phi_{k}$ es la función de onda de una sola partícula de un fermión. Puedo definir los operadores $a_{k}$ y $a^{\dagger}_{k}$ que crean y aniquilan fermiones en el estado $\phi_{k}$ cuando actúan en el estado del sistema no interactuante.
Pero, si ahora estoy interesado en el estado de N fermiones interactuando, puedo escribir la función de onda total como una combinación lineal de todos los posibles determinantes de Slater
$$ \Psi_{0}(\textbf{r}_{1}, ..., \textbf{r}_{N}) = \sum_{n} C_{n} \Phi_{n} (\textbf{r}_{1}, ..., \textbf{r}_{N}).$$
Puedo interpretar $a^{\dagger}\Phi_{0}$ como el estado de un sistema de $N$ fermiones con un fermión añadido en estado $\phi_{k}$ . Mi pregunta es: ¿Puedo hacer la misma interpretación cuando aplico uno de los operadores de creación o aniquilación a $\Psi$ ? ¿Puedo decir que $a^{\dagger}_{k}\Psi_{0}$ es un estado de $N$ fermiones que interactúan con un fermión añadido en estado $\phi_{k}$ ?