Consider the Lagrangian of a massive pseudo-scalar and a fermion field, defined as $$\mathcal{L}=\bar{\psi}(i\gamma^\mu\partial_\mu-m)\psi-1/2(\partial_\mu\theta)(\partial^\mu\theta)-M^2\theta^2-ig\bar{\psi}\gamma^5\psi\theta,$$
where $g$ is the interaction constant between the fermion field and the pseudo-scalar field. The question is the following: writing the equations of motion for the two fields following the Euler-Lagrange equation gives us a system of three coupled equations which should (I guess) be solved perturbatively. Does anyone know of a source where I can follow the solution of those equations, or a source with a similar problem solved?
What I am also curious about is that if we calculate the Feynman diagrams with two and four external pseudo-scalar lines provides us with an effective potential $V(\theta)$ where the effective coupling constant takes into account the interaction with the pseudo-scalar field, so the effectively the two fields decouple. Then we can solve the equations of motion easily. Is that a valid method for finding those solutions?