The Dirac Lagrangian reads: \begin{equation*} \mathcal{L} = \bar\psi(i\not\partial-m)\psi.\tag{1} \end{equation*} It's invariant under the transformation $\psi(x)\rightarrow e^{i\alpha}\psi(x)$. Now suppose we perform the transformation $\psi(x)\rightarrow e^{i\alpha(x)}\psi(x)$ instead. We can find \begin{equation*} \mathcal{L}\rightarrow \bar\psi e^{-i\alpha(x)}(i\not\partial-m) e^{i\alpha(x)}\psi\\ = \mathcal{L}+\bar\psi e^{-i\alpha(x)}(i\not\partial e^{i\alpha(x)})\psi\\=\mathcal{L}-(\bar\psi\gamma^\mu\psi)\partial_\mu \alpha.\tag{2} \end{equation*} I'm trying to compare the above derivation to equation 2.10 on Peskin and Schroeder: $$ \mathcal{L(x)}\rightarrow\mathcal{L(x)}+\alpha\partial_\mu\mathcal{J^\mu(x)}.\tag{2.10} $$ Thus it seems like we have the relation
$$ \mathcal{J^\mu(x)} = -\bar\psi(x)\gamma^\mu\psi(x)\tag{3} $$ Is that correct? How do I find $$j^\mu(x) = \bar\psi(x)\gamma^\mu\psi(x)\tag{4}$$ for Dirac equation? Is there a way we can show that $$\partial_\mu j^\mu(x) = 0~?\tag{5}$$ I think I got confused about how to transform from $\mathcal{J^\mu(x)}$ to $j^\mu(x)$.