I can't solve this apparent paradox; I have the free lagrangian for massive fermions $\mathscr L = i\bar\Psi\gamma^\mu\partial_\mu\Psi - m\bar\Psi\Psi$ which is invariant under the global phase transformation
$$ \Psi \rightarrow e^{i\lambda}\Psi \hspace{40pt} \bar\Psi \rightarrow e^{-i\lambda}\bar\Psi$$
Therefore I can compute the Noether conserved current
$$J^\mu = \frac{\partial\mathscr L}{\partial(\partial_\mu\Psi)}\frac{\delta\Psi}{\delta\lambda} +\frac{\partial\mathscr L}{\partial(\partial_\mu\bar\Psi)}\frac{\delta\bar\Psi}{\delta\lambda} = \bar\Psi\gamma^\mu\Psi$$
Now I see in the notes of my Professor the following calculation:
Since usually we want a real current, we can rewrite $J^\mu = i\bar\Psi\gamma^\mu\Psi$ or with explicit spinor indices $J^\mu = i\bar\Psi^\alpha(\gamma^\mu)_\alpha^{\ \ \ \beta}\Psi_\beta = i\Psi^{\dagger\sigma}(\gamma^0)_\sigma^{\ \ \ \alpha}(\gamma^\mu)_\alpha^{\ \ \ \beta}\Psi_\beta$ where $\alpha, \beta, \sigma = 1, 2, 3, 4$. We can prove this is real in fact
$$(J^\mu)^\dagger = -i(\Psi^\dagger)_\beta(\gamma^{\mu\dagger})^\beta_{\ \ \ \alpha}(\gamma^{0\dagger})^\alpha_{\ \ \ \sigma}(\Psi^\sigma)$$
but we know that for any $\psi^\alpha$ and $\chi^\alpha$ we can write
$$\psi^\alpha\chi_\alpha = \varepsilon^{\alpha\beta}\psi_\beta\chi_\alpha = - \psi_\beta\varepsilon^{\beta\alpha}\chi_\alpha = -\psi_\beta\chi^\beta$$
this means that when we raise and lower a contracted spinor index, we get a minus sign in front of our expression, hence
$$(J^\mu)^\dagger = -i(\Psi^\dagger)_\beta(\gamma^{\mu\dagger})^\beta_{\ \ \ \alpha}(\gamma^{0\dagger})^\alpha_{\ \ \ \sigma}(\Psi^\sigma)$$
$\color{red}{\mathrm{raise \ and \ lower} \ \beta}$
$$(J^\mu)^\dagger = -(\color{red}-) i(\Psi^\dagger)^\beta(\gamma^{\mu\dagger})_{\beta\alpha}(\gamma^{0\dagger})^\alpha_{\ \ \ \sigma}(\Psi^\sigma) $$
raise and lower $\color{green}\alpha$ then $\color{blue}\sigma$
$$(J^\mu)^\dagger = -(\color{red}-)(\color{green}-)(\color{blue}-) i(\Psi^\dagger)^\beta(\gamma^{\mu\dagger})_\beta^{\ \ \ \alpha}(\gamma^{0\dagger})_\alpha^{\ \ \ \sigma}(\Psi_\sigma) = + i\Psi^\dagger\gamma^{\mu\dagger}\gamma^{0\dagger}\Psi$$
now we can insert the identity $\gamma^0\gamma^0 = \Bbb I_4$ and since $\gamma^{0\dagger} = \gamma^0$ and $\gamma^0\gamma^{\mu\dagger}\gamma^0 = \bar\gamma^\mu = \gamma^\mu$ we eventually get
$$ (J^\mu)^\dagger = i\color{red}{\Psi^\dagger\gamma^0}\color{\green}{\gamma^0\gamma^{\mu\dagger}\gamma^0}\Psi = i\color{red}{\bar\Psi}\color{green}{\gamma^\mu}\Psi = J^\mu$$
we conclude that $J^\mu$ is real.
Well I can't understand this procedure since, without adding that $i$, I would have simply done in this way
\begin{aligned} (J^\mu)^\dagger & = \Psi^\dagger\gamma^{\mu\dagger}\gamma^{0\dagger}\Psi = \ \ \ \ \ \ \ \ \ \ \ \ \mathrm{inserting \ identity} \gamma^0\gamma^0 = \Bbb I_4 \\ & = \color{red}{\Psi^\dagger\gamma^0}\color{green}{\gamma^0\gamma^{\mu\dagger}\gamma^0}\Psi = \\ & = \color{red}{\bar\Psi}\color{green}{\gamma^\mu}\Psi = J^\mu \end{aligned}
and I would have done in this way also because if I were to apply the same procedure used by my professor to the kinetic term of the lagrangian I would get
$$[i\bar\Psi\gamma^\mu\partial_\mu\Psi]^\dagger = [i\Psi^{\dagger\sigma}(\gamma^0)_\sigma^{\ \ \ \alpha}(\gamma^\mu)_\alpha^{\ \ \ \beta}(\partial_\mu\Psi)_\beta]^\dagger = - i(\partial_\mu\Psi^\dagger)_\beta(\gamma^{\mu\dagger})^\beta_{\ \ \ \alpha}(\gamma^{0\dagger})^\alpha_{\ \ \ \sigma}\Psi^\sigma$$
now I raise and lower $\color{red}\beta$, $\color{green}\alpha$ and $\color{blue}\sigma$ and i get
$$-(\color{red}-)(\color{green}-)(\color{blue}-) i(\partial_\mu\Psi^\dagger)^\beta(\gamma^{\mu\dagger})_\beta^{\ \ \ \alpha}(\gamma^{0\dagger})_\alpha^{\ \ \ \sigma}\Psi_\sigma = +i\partial_\mu\Psi^\dagger\gamma^{\mu\dagger}\gamma^{0\dagger}\Psi$$
Now I can insert the identity $\gamma^0\gamma^0 = \Bbb I_4$ and then integrate by parts
$$i\partial_\mu\Psi^\dagger\gamma^{\mu\dagger}\gamma^{0\dagger}\Psi = -i\Psi^\dagger\gamma^0\gamma^0\gamma^{\mu\dagger}\gamma^0\partial_\mu\Psi + i\partial_\mu(\Psi^\dagger\gamma^0\gamma^0\gamma^{\mu\dagger}\gamma^0\Psi) = -i\bar\Psi\gamma^\mu\partial_\mu\Psi + i\partial_\mu(\bar\Psi\gamma^\mu\Psi)$$
i.e. up to a total derivative I get
$$(i\bar\Psi\gamma^\mu\partial_\mu\Psi)^\dagger = - i\bar\Psi\gamma^\mu\partial_\mu\Psi $$
Therefore with this procedure the lagrangian (i.e. the action) would no longer be real. What is wrong with this?