 Axial-current divergence for a pseudoscalar Yukawa interaction

For ${\mathcal{L}}_{\mathrm{int}}=-g\varphi \bar{\psi }{\gamma }^5\psi$, the variational Dirac equation and adjoint Dirac equation are $(i/\partial -m-g\varphi {\gamma }^5)\psi =0$ and $i({\partial }_{\mu }\bar{\psi }){\gamma }^{\mu }+m\bar{\psi }+g\varphi \bar{\psi }{\gamma }^5=0$. Substitute these into the derivative of $j_5^{\mu }=\bar{\psi }{\gamma }^{\mu }{\gamma }^5\psi$. Anticommuting the chirality matrix past ${\gamma }^{\mu }$ gives two identical mass terms and two identical interaction terms, proving the displayed identity. These are classical field equations; the formula does not assert an unrenormalized quantum composite-operator identity. With Hermitian coupling $g=i{g}_P$, the second term is $-2g_P\varphi \bar{\psi }\psi$.