I tried to give a QFT description of the hydrogen atom (at least conceptually and without going into mathematical details). The system is made by a spin$\frac{1}{2}$ fermion (electron) field $\phi_s$ coupled to the external electromagnetic field of the proton.
I can build a generic state of the electron field by acting on the vacuum with creation operators. At time $t=0$ a generic 1-electron state should be something like
$$\sum_{s} \int d^{3} \pmb{k} \, f(\pmb{k}, s) \hat{a} _{\pmb{k}, s}^{\dagger} |0\rangle$$
By choosing an appropriate $f$ it is possible to minimize the energy of the 1-electron state. I shall work with this particular state from now on and I will refer to it as $|gs\rangle$.
What if I want to predict (with a certain probability of course) where the electron is after a position measurement? Since $\pmb{x}$ is just a continuous index the ket $|\pmb{x}\rangle$ doesn't make much sense therefore I can't use it. My idea was to do the following procedure:
$|gs\rangle = \sum_{s} \int d^3 \pmb{x} \langle\phi_{s}(\pmb{x}) |gs\rangle |\phi_{s}(\pmb{x})\rangle $, where $|\phi_{s}(\pmb{x})\rangle \equiv \hat{\phi}_{s}(\pmb{x}) |0\rangle $. Since $|\phi_{s}(\pmb{x})\rangle$ describes a localized excitation of the field I thought I could interpret $\sum_{s} \lvert \langle\phi_{s}(\pmb{x}) |gs\rangle \rvert^{2} $ as the probability distribution of the electron in space.
Does what I did even make sense? Will $\langle\phi_{s}(\pmb{x}) |gs\rangle$ be related to the lowest energy solution of the Dirac equation (interpreted as the wave function of the electron)? Would lamb shift appear if I were to do the math?