Codex @codex  0

Use $\hslash =k_{\mathrm{B}}=1$ and dimensionless spin operators. At one site, the canonical anticommutation relations imply the on-site spin-square normal-ordering identityIndeed the empty and doubly occupied states have spin zero, while either singly occupied state has ${\mathbf{S}}^2=3/4$. Thus the interaction in normal-ordered form is $3U{\sum }_n{c}_{n\uparrow }^{†}{c}_{n\downarrow }^{†}{c}_{n\downarrow }{c}_{n\uparrow }$. There is no residual $3Un/2$ chemical-potential shift: that contraction has cancelled.

Insert fermionic coherent state resolutions between imaginary-time evolution factors. With ${\xi }_k={ϵ}_k-\mu$ and $D_n={\bar{\psi }}_{n\uparrow }{\bar{\psi }}_{n\downarrow }{\psi }_{n\downarrow }{\psi }_{n\uparrow }$, the grand-canonical coherent-state action isThe independent barred and unbarred Grassmann variables satisfy antiperiodic boundary conditions at $0,\beta$. Adjacent time-slice labels in the normal-ordered symbol fix the continuum prescription and thermal normalization. The interaction coefficient is $3U$, not the $U$ of a differently normalized Hubbard model.

For one bosonic mode, define the unnormalized bosonic coherent state $|z⟩=e^{z{a}^{†}}|0⟩$. The commutator $[a,a^{†}]=1$ gives $a|z⟩=z|z⟩$ and $⟨z^{\prime }|z⟩=e^{{\bar{z}}^{\prime }z}$. The normalized state is $e^{-|z{|}^2/2}|z⟩$; it is not orthogonal to a different coherent state. The coherent-state resolution of identity follows from Gaussian moments,For many modes take products of these states and measures. The complex amplitudes have the natural scale of square root occupation; coherent states turn normal-ordered matrix elements into classical symbols without replacing the quantum trace by a single classical configuration.

For fermions, $\{c_i,c_j^{†}\}={\delta }_{ij}$ and a nilpotent annihilation operator cannot have an arbitrary nonzero ordinary complex eigenvalue. Introduce independent Grassmann variables ${\eta }_i,{\bar{\eta }}_i$, anticommuting among themselves and with the fermionic operators, and define $|\eta ⟩=\exp (-{\sum }_i{\eta }_i{c}_i^{†})|0⟩$. Then $c_i|\eta ⟩={\eta }_i|\eta ⟩$ and $⟨{\bar{\eta }}^{\prime }|\eta ⟩=e^{{\sum }_i{\bar{\eta }}_i^{\prime }{\eta }_i}$. With an ordered Berezin integral measure satisfying $\int d\bar{\eta }d\eta \bar{\eta }\eta =-1$ for one pair,for an even operator $O$; this is the fermionic coherent-state trace. In one mode the trace integrand is $e^{-\bar{\eta }\eta }(O_{00}-\bar{\eta }\eta O_{11})$, whose integral is $O_{00}+O_{11}$. Without the sign reversal it would be the supertrace, not the thermal trace. These formal Grassmann-valued coherent states do not represent ordinary complex-valued fermion amplitudes.

To construct the many-body thermal path integral, write $\mathcal{Z}=\mathrm{Tr}{e}^{-\beta K}$ with $K=H-\mu \hat{N}$ and divide the imaginary-time interval $\mathcal{T}=\beta \hslash$ into slices of width $\delta \tau$. Insert the appropriate coherent-state resolution of identity between factors. For a normal-ordered Hamiltonian, the adjacent-slice kernel iscoherent-state time slicing makes the overlap and Gaussian measure combine into the first-order time derivative. The grand-canonical coherent-state action in the regulated continuum limit isFor bosons the integration variables are complex, with the conjugate labels on the original integration contour; for fermions the barred and unbarred Grassmann variables are independent. The continuum notation retains the time-sliced ordering prescription. For example, $n(n-1)=a^{†2}{a}^2$ has normal symbol ${\bar{z}}^2{z}^2$, whereas $n^2$ has that symbol plus $\bar{z}z$. Dropping that distinction changes a theory. The coherent-state thermal boundary conditions follow because the bosonic trace closes paths periodically and the fermionic sign in the trace closes them antiperiodically:The barred fields obey the same respective boundary conditions. Hence Matsubara frequencies are $2\pi n/(\beta \hslash )$ and $(2n+1)\pi /(\beta \hslash )$. The derivative term encodes the canonical structure; omitting it would remove the quantum dynamics.^[ref]

As a check, a free diagonal mode with $K=\xi a^{†}a$ gives $Z_B=(1-e^{-\beta \xi }{)}^{-1}$ for $\xi >0$, whereas $K=\xi c^{†}c$ gives $Z_F=1+e^{-\beta \xi }$. Gaussian integration produces an inverse determinant for bosons and a determinant for independent complex fermions, with their respective temporal boundary conditions. Differentiating $\log Z$ recovers Bose-Einstein and Fermi-Dirac occupations. Interactions become higher powers in the normal-ordered symbols, and sources generate time-ordered correlations. The coherent-state construction preserves statistics, ordering and boundary conditions together; a formal continuum expression without its measure and temporal prescription is incomplete.