It's the whole day that I try to figure out at what point exactly the phase factor appearing in both creation and annihilation operators, went away maybe absorbed inside occupation number states. Here the procedure I followed: sorry in advanced for the "heavy" notation, I'm dumb and I need to have all informations in front of my eyes all the time.
So everything starts from these commutation relations \begin{gather*} \hat{n}_\mu(t,\boldsymbol{k}) \doteq \hat{a}_\mu^\dagger(t,\boldsymbol{k}) \hat{a}_\mu(t,\boldsymbol{k}) \\ \begin{cases} \left[ \hat{a}_\mu(t,\boldsymbol{k}), \hat{n}_\nu\left(\tilde{t},\tilde{\boldsymbol{k}}\right) \right] = \delta_{\mu\nu} \delta\left(t,\tilde{t}\right) \delta\left(\boldsymbol{k},\tilde{\boldsymbol{k}}\right) \hat{a}_\nu\left(\tilde{t},\tilde{\boldsymbol{k}}\right) \\ \left[ \hat{a}_\mu^\dagger(t,\boldsymbol{k}), \hat{n}_\nu\left(\tilde{t},\tilde{\boldsymbol{k}}\right) \right] = - \delta_{\nu\mu} \delta\left(t,\tilde{t}\right) \delta\left(\boldsymbol{k},\tilde{\boldsymbol{k}}\right) \hat{a}_\nu^\dagger\left(\tilde{t},\tilde{\boldsymbol{k}}\right) \end{cases} \end{gather*} Because of that, after some calculations \begin{gather*} \begin{cases} \hat{a}_\mu(t,\boldsymbol{k}) |n_\mu(t,\boldsymbol{k})\rangle \equiv c^-_\mu(t,\boldsymbol{k}) |(n-1)_\mu(t,\boldsymbol{k})\rangle \\ \hat{a}_\mu^\dagger(t,\boldsymbol{k}) |n_\mu(t,\boldsymbol{k})\rangle \equiv c^+_\mu(t,\boldsymbol{k}) |(n+1)_\mu(t,\boldsymbol{k})\rangle \end{cases} \\ \begin{cases} |c^-_\mu(t,\boldsymbol{k})|^2 = n_\mu(t,\boldsymbol{k}) \\ |c^+_\mu(t,\boldsymbol{k})|^2 = n_\mu(t,\boldsymbol{k})+1 \end{cases} \\ \begin{cases} \hat{a}_\mu(t,\boldsymbol{k}) |n_\mu(t,\boldsymbol{k})\rangle = e^{i \alpha^-_\mu\left(n_\mu(t,\boldsymbol{k}),t,\boldsymbol{k}\right)} \sqrt{n_\mu(t,\boldsymbol{k})} |(n-1)_\mu(t,\boldsymbol{k})\rangle \\ \hat{a}_\mu^\dagger(t,\boldsymbol{k}) |n_\mu(t,\boldsymbol{k})\rangle = e^{i \alpha^+_\mu\left(n_\mu(t,\boldsymbol{k}),t,\boldsymbol{k}\right)} \sqrt{n_\mu(t,\boldsymbol{k})+1} |(n+1)_\mu(t,\boldsymbol{k})\rangle \end{cases} \end{gather*} where $\alpha^-,\alpha^+$ are dinstinct functions with $(n_\mu(t,\boldsymbol{k}),t,\boldsymbol{k})$ arguments (right?). Using the above commutation relations I'm only able to prove that, unless a multiple of $2\pi$, the following equality holds \begin{equation*} \alpha^+\left(n_\mu(t,\boldsymbol{k}),t,\boldsymbol{k}\right) = -\alpha^-\left((n+1)_\mu(t,\boldsymbol{k}),t,\boldsymbol{k}\right) \end{equation*} but what I desire to prove is also that \begin{equation*} \alpha^+\left(n_\mu(t,\boldsymbol{k}),t,\boldsymbol{k}\right) = \alpha^-\left(n_\mu(t,\boldsymbol{k}),t,\boldsymbol{k}\right) \end{equation*} Is there a way to prove it or is not possible? In this other case how can I make the phase factor go away?