When working with path integrals of both bosonic and fermionic field variables, I'm a bit unsure of how to do the usual complete the square trick when an interaction between the two is concerned. Say you have a generic partition function like \begin{equation}Z=\int D\phi D\bar\psi D\psi \,e^{iS[\phi,\bar\psi,\psi]}\end{equation} where the action has your standard quadratic part, but also an interaction term like $\mathcal{L}_{int}\sim \bar\psi \phi\psi$. If you wanted to integrate out the Bose fields first, you would complete the square but that would involve making a substitution like $\phi\rightarrow \phi'= \phi-\bar\psi\psi $ and integrating over $\phi'$ (I may be wrong here). I'm confused on what it means to:
- Subtract a product of Grassmann numbers from an ordinary number (they aren't the same kind of number)
- Integrate over the new Bose field, which is technically a function of the Grassmann variables (are there any subtleties with that process?)
Any help would be greatly appreciated.