I am wondering if there exists a method to compute path integrals using the stationary phase method for theories with both bosons and fermions. (I am aware of such a method for theories with bosons only.)
One of the most basic problems which bothers me is how to interpret stationary points of the action (or, in other words, solutions of the equations of motion). For examples let us consider the Lagrangian of QED: $$\mathcal{L}=-\frac{1}{4} F_{\mu\nu}F^{\mu\nu}+\bar\psi(i\gamma^\mu D_\mu-m)\psi.\tag{1}$$ The equations of motions are \begin{eqnarray*} \partial_\nu F^{\nu\mu}=e\bar\psi\gamma^\mu\psi,\tag{2}\\ (i\gamma^\mu D_\mu-m)\psi=0.\tag{3} \end{eqnarray*}
I am not sure how to interpret the first equation: the left hand side is bosonic, while the right hand side is fermionic (however it is even as a product of two fermions).