Fermionic elimination [fdip-009I]
✍️sourceAGENTDRAFTED
- November 1, 2026
Fermionic elimination [fdip-009I]
✍️sourceAGENTDRAFTED
- November 1, 2026
What remains after some fermion variables are eliminated? The paired Grassmann integral replaces a Gaussian normalizer by a determinant, but the effective quadratic form is still a Schur complement. The eliminated block must be invertible; the full matrix need not be positive or invertible for this algebraic identity.
A commuting complex field supplies the inverse determinant when the Hermitian part of its matrix is positive. This is an ordinary convergent integral on twice as many real coordinates. Its cancellation with the Grassmann determinant uses the same matrix in both quadratic forms. A non-Hermitian accretive matrix usually gives a complex weight, so this cancellation does not construct a positive joint probability law.
The logarithm of the determinant needs an additional choice. On the accretive cone there is a unique analytic logarithm normalized at the identity; it need not be the principal scalar logarithm of the determinant. The four-cycle makes these distinctions concrete with the same effective matrix as in real Gaussian elimination.
Theorem 1. Partial Grassmann integration and the Schur complement [fdip-009J]AGENTDRAFTED
- November 1, 2026
Theorem 1. Partial Grassmann integration and the Schur complement [fdip-009J]AGENTDRAFTED
- November 1, 2026
Let \(R,U\) be disjoint ordered finite sets. Order the independent paired generators with all retained \(R\)-pairs first and all eliminated \(U\)-pairs second. For \(n=|R|+|U|\), let \[K=\begin {pmatrix}A&B\\C&D\end {pmatrix} \in M_n(\mathbb C).\] Require \(D=K_{UU}\) to be invertible and put \(S=A-BD^{-1}C\). No adjoint relation between \(B,C\) is assumed. Define \(J_U\) by taking the coefficient of the ordered \(U\)-pair top word after moving retained generators to the left. The result is an element of the retained exterior algebra. Define \(J_R\) by its paired top word, so \(J_RJ_U=J_n\). The eliminated top word has even degree, so moving it past a retained coefficient gives no further sign. Then \[J_U\left (e^{-\bar \psi K\psi }\right ) =\det D\,e^{-\bar \psi _R S\psi _R},\] and consequently \[\det K=\det D\det S.\] The empty-block conventions are determinant one and the identity coefficient map. Invertibility of \(K\) or \(S\) is unnecessary.
Proof.
- November 1, 2026
Proof.
- November 1, 2026
Introduce the odd linear combinations \[\begin {aligned} \chi _U&=\psi _U+D^{-1}C\psi _R,\\ \bar \chi _U&=\bar \psi _U+\bar \psi _RBD^{-1}. \end {aligned}\] The bars remain labels of independent generators. Expanding in the displayed multiplication order gives \[\bar \psi K\psi =\bar \psi _R S\psi _R+\bar \chi _U D\chi _U.\] For example, the last term of the expansion of \(\bar \chi _U D\chi _U\) is \(\bar \psi _RBD^{-1}C\psi _R\), which cancels the subtraction in \(S\). All these quadratic elements are even. They commute, so their finite exponentials multiply.
Replacing an eliminated generator by itself plus a retained odd linear combination preserves its partial top coefficient. Indeed, expand each exterior monomial in the new symbols: every term using a retained replacement has fewer eliminated generators, while the term using no replacement has coefficient one and its original order. Terms of lower eliminated degree cannot contribute to the top degree. Thus the coefficient in the shifted quadratic exponential equals the coefficient for independent \(U\)-generators. The paired determinant identity evaluates it as \(\det D\). The retained exponential factors on the left, proving the partial formula. Applying \(J_R\) and the same determinant identity proves the full factorization.
The effective matrix is the same block expression as in real positive Gaussian elimination. The scalar here is \(\det D\); the real Gaussian fiber integral instead has \((2\pi )^{|U|/2}(\det D)^{-1/2}\) under its real positive hypotheses and volume convention. Neither positivity nor an ordinary integration measure was used in the Grassmann proof.
The finite block identity is also supplied by
Matrix.det_fromBlocks₂₂, for finite index types over a
commutative ring with an invertible eliminated block. This
formal statement supports the determinant factorization, not
the partial Grassmann coefficient operation.
Theorem 2. Accretive complex Gaussian integrals and cancellation [fdip-009K]AGENTDRAFTED
- November 1, 2026
Theorem 2. Accretive complex Gaussian integrals and cancellation [fdip-009K]AGENTDRAFTED
- November 1, 2026
Call \(K\in M_n(\mathbb C)\) strictly accretive when its Hermitian part \(H=(K+K^*)/2\) is positive definite. For \(z_j=x_j+iy_j\), use ordinary Lebesgue volume \(d^{2n}z=\prod _j dx_j\,dy_j\). Then \(K\) is invertible, the following integral converges absolutely, and \[I(K):=\int _{\mathbb C^n}e^{-z^*Kz}\,d^{2n}z =\frac {\pi ^n}{\det K}.\] Consequently, with \(Z_B(K)=\pi ^{-n}I(K)\) and the paired Grassmann value \(Z_F(K)=J_n(e^{-\bar \psi K\psi })\), \[Z_B(K)Z_F(K)=1.\] For \(n=0\) both factors are one.
Proof.
- November 1, 2026
Proof.
- November 1, 2026
For \(n>0\), choose \(h>0\) with \(H\ge hI\). Then \[|e^{-z^*Kz}|=e^{-z^*Hz}\le e^{-h\|z\|^2}.\] The right side is integrable by the product of scalar Gaussian integrals. If \(Kz=0\), its real quadratic value is \(z^*Hz=0\), so \(z=0\); finite dimension makes \(K\) invertible.
First suppose \(K=K^*>0\). A unitary eigenbasis gives positive eigenvalues \(\lambda _j\). The unitary change of variables preserves real volume, since its real Jacobian has absolute determinant \(|\det _{\mathbb C}U|^2=1\). In each complex coordinate the scalar Gaussian formula gives \[\int _{\mathbb R^2}e^{-\lambda _j(x^2+y^2)}dx\,dy =\frac \pi {\lambda _j}.\] Multiplication gives \(I(K)=\pi ^n/\det K\).
For general \(K\), write \(K=H+A\), where \(A^*=-A\), and set \(K(w)=H+wA\) for a complex parameter \(w\). Its Hermitian part is \(H+i(\operatorname {Im}w)A\). Choose \(\epsilon >0\) such that \(\epsilon \|A\|<h/2\), with any positive choice if \(A=0\). On the connected strip \(|\operatorname {Im}w|<\epsilon \), this part is at least \(hI/2\), so \(K(w)\) is invertible.
The integral \(I(K(w))\) is holomorphic on that strip. Differentiation of the integrand brings down powers of \(-z^*Az\); on every smaller neighborhood these are bounded by a polynomial in \(\|z\|\) times \(e^{-h\|z\|^2/2}\). This integrable bound justifies differentiation under the integral. The function \(\pi ^n/\det K(w)\) is holomorphic there too. For \(w=it\) with real \(|t|<\epsilon \), the matrix \(H+itA\) is Hermitian positive, so the preceding calculation identifies the two functions. These points accumulate at the interior point zero. The one-variable identity theorem therefore identifies the functions throughout the strip, in particular at \(w=1\). The Grassmann determinant then proves cancellation.
For Hermitian positive \(K\), the normalized density \((\det K/\pi ^n)e^{-z^*Kz}\) is a genuine Gaussian probability law on \(2n\) real coordinates. For general accretive \(K\), the weight can be complex: \(K=(1+i)\) gives \(Z_F=1+i\) and \(Z_B=(1+i)^{-1}\). Normalization and determinant cancellation still hold, but there is no positive Gaussian law with this density. A single real Gaussian has only a determinant power \(-1/2\); Grassmann algebra by itself supplies no convergence for a commuting integral.
Accretivity also survives the Schur operation. Its principal block \(D\) is accretive by testing \(K\) on \((0,u)\), hence invertible. For \(v\ne 0\), put \(x=(v,-D^{-1}Cv)\). Then \(Kx=(Sv,0)\), and \[\operatorname {Re}(v^*Sv) =\operatorname {Re}(x^*Kx)>0.\] Thus \(S\) is accretive. This extra conclusion uses accretivity of the full matrix, not just invertibility of \(D\).
[wingate2025lattice, Section 9.1, equations (9.3)-(9.5), p. 19] describes integrating a fermion matrix and re-expressing its determinant with commuting complex pseudofermions. Its equation (9.5) uses \(K^{-1}\) in the commuting quadratic form, yielding a determinant rather than its inverse when convergence holds. The cancellation here uses \(K\) in both forms and the specified factor \(\pi ^{-n}\); the proof above states its convergence hypotheses explicitly.
Proposition 3. The logarithm on the accretive cone [fdip-009L]AGENTDRAFTED
- November 1, 2026
Proposition 3. The logarithm on the accretive cone [fdip-009L]AGENTDRAFTED
- November 1, 2026
On the open convex cone \(\mathcal C_n\) of strictly accretive complex matrices, there is a unique holomorphic function \(\mathcal L\) such that \[e^{\mathcal L(K)}=\det K,\qquad \mathcal L(I)=0.\] For \(K_t=I+t(K-I)\), it is given by \[\mathcal L(K)=\int _0^1 \operatorname {tr}\left [K_t^{-1}(K-I)\right ]dt.\] On Hermitian positive matrices it equals the real \(\log \det K\). Thus the fermionic logarithmic contribution \(W_F(K)=-\mathcal L(K)\) is unambiguous on this cone.
Proof.
- November 1, 2026
Proof.
- November 1, 2026
Write \(K_t=I+t(K-I)\). Convexity makes \(K_t\) accretive and invertible for every \(0\le t\le 1\). In a neighborhood of each \(K\in \mathcal C_n\), the Hermitian parts of all these matrices have a common positive lower bound. Their inverses and derivatives in the matrix entries are therefore locally bounded uniformly in \(t\). The integral of the holomorphic rational integrand is holomorphic in \(K\).
Put \(g(t)=\det K_t\). The determinant derivative, as in the real log-determinant calculation, gives \[\frac {g'(t)}{g(t)} =\operatorname {tr}(K_t^{-1}(K-I)).\] This complex version follows directly from \(\det (K_t+hB)=\det K_t\det (I+hK_t^{-1}B)\) and \(\det (I+hM)=1+h\operatorname {tr}M+O(h^2)\), obtained by the permutation expansion. Therefore the derivative of \(g(t)\exp (-\int _0^t g'(s)/g(s)\,ds)\) is zero. Its value at zero is one, proving \(\exp \mathcal L(K)=\det K\); the identity matrix gives \(\mathcal L(I)=0\).
Two such logarithms differ by a continuous function with values in \(2\pi i\mathbb Z\). The cone is connected, so this function is constant, and normalization makes it zero. On Hermitian positive \(K\), all \(K_t\) are positive and \(g(t)>0\); ordinary real logarithmic differentiation identifies the integral with \(\log \det K\).
This branch is not always the principal scalar logarithm of \(\det K\). For \(K=e^{i\pi /3}I_4\), the Hermitian part is \(I_4/2>0\). Along \(K(s)=e^{is}I_4\), \(0\le s\le \pi /3\), continuity and normalization give \(\mathcal L(K(s))=4is\). Hence \[\mathcal L(e^{i\pi /3}I_4)=\frac {4\pi i}{3},\] whereas the principal scalar logarithm of its determinant is \(-2\pi i/3\). The determinant can even reach the scalar principal branch cut while the matrix stays accretive.
For a general invertible complex matrix, writing \(-\log \det K\) requires a declared branch on the domain under consideration; there is no single-valued continuous logarithm on all invertible matrices when \(n\ge 1\). The loop \(\operatorname {diag}(e^{it},1,\ldots ,1)\), \(0\le t\le 2\pi \), already has determinant winding once. If \(\det K=0\), no logarithm exists, although the Grassmann determinant identity still does. The boson-fermion cancellation itself uses no logarithm branch.
Example 4. The same four-cycle Schur matrix with fermions [fdip-009M]AGENTDRAFTED
- November 1, 2026
Example 4. The same four-cycle Schur matrix with fermions [fdip-009M]AGENTDRAFTED
- November 1, 2026
Take the alternating four-cycle with \(t,m>0\), set \(a=m^2\), \(s=1+t+a\), and use \(K=Q_{t,m}=L_t+aI\) now as a complex matrix. Order retained vertices \(R=(0,2)\) before eliminated vertices \(U=(1,3)\). Its blocks are \[A=D=sI_2,\qquad B=C=-\begin {pmatrix}1&t\\t&1\end {pmatrix}.\] The fermionic Schur formula gives \[J_U(e^{-\bar \psi K\psi }) =s^2e^{-\bar \psi _R S\psi _R},\] where \[\begin {aligned} S&=\begin {pmatrix}\alpha &-\beta \\-\beta &\alpha \end {pmatrix},\\ \alpha &=s-\frac {1+t^2}{s},\qquad \beta =\frac {2t}{s}. \end {aligned}\] This is exactly the real Gaussian effective precision. The fermionic scalar is \(s^2\); the real Gaussian eliminated factor in the cited example is \(2\pi /s\).
The two eigenvalues of \(S\) are \[\frac {a(a+2+2t)}s,\qquad \frac {(a+2)(a+2t)}s.\] Writing \(\mathcal D=a(a+2)(a+2t)(a+2+2t)\), it follows that \[\det S=\frac {\mathcal D}{s^2},\qquad Z_F(K)=s^2\det S=\mathcal D.\] Positivity of \(K\) also allows the ordinary complex integral with normalized volume \(\pi ^{-4}d^8z\). The complex Gaussian theorem gives \[Z_B(K)=\mathcal D^{-1},\qquad Z_F(K)Z_B(K)=1.\] In contrast, the original real Gaussian, with exponent \(-\phi ^TK\phi /2\) and volume \(d^4\phi \), has partition function \((2\pi )^2/\sqrt {\mathcal D}\). These powers and constants refer to different field spaces and measures.
At \(t=m=1\), \(s=3\), \(\mathcal D=45\), and \[S=\begin {pmatrix}7/3&-2/3\\-2/3&7/3\end {pmatrix}, \qquad \det S=5.\] The Grassmann eliminated factor \(9\) times the retained factor \(5\) gives \(45\), while the normalized complex-boson value is \(1/45\). At \(t=2,m=1\), \(s=4\), \(\det S=105/16\); the product is \(16\cdot 105/16=105\). These are the same two effective matrices compared in the weighted-cycle example.
At \(m=0,t>0\), the eliminated block \((1+t)I_2\) is still invertible and the partial Grassmann formula is valid. The retained Schur matrix has a constant zero mode, so \(\det K=0\): the full Grassmann integral vanishes. The commuting full Gaussian instead diverges in its constant mode. There is no finite inverse determinant, cancellation identity or logarithmic contribution at this massless point.