 Axial-current divergence for scalar and pseudoscalar backgrounds

For $\mathcal{L}=\bar{\psi }(i/\partial -S+iP{\gamma }^5)\psi$, with real scalar backgrounds $S,P$, the Dirac equation is $(i/\partial -S+iP{\gamma }^5)\psi =0$ and the adjoint Dirac equation is $i({\partial }_{\mu }\bar{\psi }){\gamma }^{\mu }+S\bar{\psi }-iP\bar{\psi }{\gamma }^5=0$. Use $\{{\gamma }^5,{\gamma }^{\mu }\}=0$ and $({\gamma }^5{)}^2=1$ in the derivative of $j_5^{\mu }=\bar{\psi }{\gamma }^{\mu }{\gamma }^5\psi$ to obtain the displayed classical identity. Taking $S=m+g\varphi$ and $P=G\Phi$ includes both scalar and Hermitian pseudoscalar Yukawa couplings.