For the set of all fermionic field operators $\Psi(x) | x \in \mathbb{R}^{3 +1}$, we won't find a $|\phi \rangle$ that is an eigenstate to the complete set of field operators, unless we make use of Grassmann numbers and construct a fermionic coherent state. However, we can do canonical quantization or nonrelativistic many-body quantum mechanics with fermionic field operators just fine, without ever hearing anything about Grassmann numbers.
Now I am wondering: Do I NOT encounter Grassmann numbers because I don't dig deep enough, or is it that a usual Hilbert space over complex numbers is fully fine and sufficient to also describe the set of $\Psi(x) | x \in \mathbb{R}^{3+1}$, with the $\Psi(x)$ all anticommuting?
I tried to specify ladder operators using only complex numbers, and for a finite number of states, this works pretty well (see for example this answer), but I don't know if I run into problems when I try the same with an uncountably infinite number of base states (one for every position, for example).
To make the question more specific: If I try to write down matrix elements, would I encounter Grassmann numbers? Is there another way that Grassmann numbers might appear?
In this question one answer is that the operators should already be operators on a supervectorspace - Is that a contradiction to my guess that we don't need them in general, or is it that due to the very specific question? (The OP of the other question specifically asked wether it's possible to construct creation and annihilation operators that also anticommute with grassmann numbers)?