Ordered Grassmann Gaussian contractions ($⟨{\theta }_i{\bar{\theta }}_j{\theta }_k{\bar{\theta }}_l⟩=C_{ij}{C}_{kl}-C_{il}{C}_{kj},C=B^{-1}$)

Use the Berezin integral orientation $\int {\prod }_{i}d{\bar{\theta }}_{i}d{\theta }_i{\theta }_1{\bar{\theta }}_1\cdots {\theta }_n{\bar{\theta }}_n=1$. For an invertible matrix $B$, the normalized weight $e^{-\bar{\theta }B\theta }$ has $⟨{\theta }_i{\bar{\theta }}_j⟩=(B^{-1}{)}_{ij}$. Differentiate its determinant with respect to $B_{ji}$ to obtain this formula; two derivatives give the displayed four-point formula. The second pairing carries a fermionic sign. These are ordered contractions, so exchanging two odd factors changes their sign.

For $C=B^{-1}$, ordered Grassmann Gaussian contractions give $⟨\bar{\theta }M\theta ⟩=-\mathrm{tr}(MC)$ and $⟨(\bar{\theta }M\theta )(\bar{\theta }N\theta )⟩=\mathrm{tr}(MC)\mathrm{tr}(NC)-\mathrm{tr}(MCNC)$. The disconnected contribution has two fermion loops; the connected contribution has one, giving its minus sign. This covariance is an algebraic normalized Berezin integral, rather than a positivity statement about a probability measure.