 Solution

Let $K=H-\mu N$, take $M$ slices of width $ϵ=\beta /M$, and insert the coherent-state resolution of identity between factors of $e^{-ϵK}$. If $H$ is in normal ordering, each short-time matrix element isCombining the overlap with the Gaussian weight leaves ${\sum }_{j,n}{\bar{\psi }}_{n,j}({\psi }_{n,j}-{\psi }_{n,j-1})$ in the action. The coherent-state time slicing prescription fixes which adjacent labels appear in $H$, rather than allowing an arbitrary ordering change after taking the continuum limit.

The twist $\zeta$ in the thermal trace closes the path with the coherent-state thermal boundary conditions: periodic for bosons, antiperiodic for fermions. Taking the regulated continuum limit givesThe action is dimensionless because $\tau$ has inverse-energy units. For fermions, $\bar{\psi }$ and $\psi$ remain independent Grassmann fields. If the original Hamiltonian is not normally ordered, first express it in normal order, retaining all constants.