Grassmann Variables and Berezin Integration

Grassmann variables are anticommuting generators (

θ_{i} θ_{j} = - θ_{j} θ_{i}

) of an exterior (Grassmann) algebra, and Berezin integration is the unique linear functional on that algebra that behaves like both differentiation and integration. Key result:

\int d θ θ = 1

\int d θ 1 = 0

, and for a quadratic form

\int d^{n} \overset{ˉ}{θ} d^{n} θ e^{- \overset{ˉ}{θ}^{T} A θ} = det A

(note: determinant, not its inverse). Applies to: fermionic path integrals, supersymmetry, supermanifolds, random matrix theory, and the algebraic theory of determinants. Who/when: Hermann Grassmann (1844); formal integration by Felix Berezin (1966).

Grassmann (anticommuting) variables provide the algebraic framework for describing fermionic degrees of freedom, where the Pauli exclusion principle is encoded directly in the algebra $θ^{2} = 0$ . A Grassmann algebra over a field is the exterior algebra of a vector space; its elements are polynomials in nilpotent, mutually anticommuting generators that terminate after finitely many terms. Berezin integration is the calculus invented to integrate over these variables — remarkably, integration and differentiation coincide. The single most surprising fact: whereas ordinary Gaussian integration over commuting variables yields $(det A)^{- 1/2}$ , the Berezin–Gaussian integral over Grassmann variables yields $(det A)^{+ 1}$ , the inverse power. This sign-reversal in the exponent of the determinant is the mathematical heart of the boson–fermion distinction.

The Core Insight

A Grassmann algebra is finite-dimensional as a vector space because the generators square to zero: any function of $n$ Grassmann variables is a polynomial whose highest term is $θ_{1} θ_{2} \dots θ_{n}$ . Consequently the only meaningful linear operation respecting translation invariance is to pick off the top coefficient — integration = extracting the coefficient of the top-form. The same operation differentiates. This collapse of analysis into algebra is what makes fermionic path integrals exact and combinatorial.

Formal Statement

Grassmann Algebra

Let $V$ be an $n$ -dimensional vector space over a field $K$ (usually $R$ or $C$ ) with basis ${θ_{1}, \dots, θ_{n}}$ . The Grassmann algebra (exterior algebra) $Λ (V) = ⨁_{k = 0}^{n} Λ^{k} (V)$ is the associative unital $K$ -algebra generated by the $θ_{i}$ subject to:

θ_{i} θ_{j} + θ_{j} θ_{i} = 0 for all i, j, in particular θ_{i}^{2} = 0.

⟳Impose the anticommutation relation as the defining quotient of the free tensor algebra

T (V)

by the two-sided ideal generated by

θ_{i} \otimes θ_{j} + θ_{j} \otimes θ_{i}

A general element ("superfunction") has the finite expansion

f (θ) = k = 0 \sum n i_{1} < \dots < i_{k} \sum c_{i_{1} \dots i_{k}} θ_{i_{1}} \dots θ_{i_{k}}, c \in K,

⟳Reorder any monomial using [1]; repeated indices vanish by nilpotency, so distinct ordered index sets form a basis of

dim Λ (V) = 2^{n}

The algebra is $Z_{2}$ -graded: $f = f_{even} + f_{odd}$ , where even (commuting) elements form the body+nilpotent center and odd elements anticommute among themselves.

Grassmann Differentiation

The left derivative $\partial / \partial θ_{i}$ is the odd linear operator defined on monomials by anticommuting $θ_{i}$ to the far left and stripping it:

\frac{\partial}{\partial θ _{i}} (θ_{j_{1}} \dots θ_{j_{k}}) = a = 1 \sum k (- 1)^{a - 1} δ_{i j_{a}} θ_{j_{1}} \dots θ_{j_{a}} \dots θ_{j_{k}},

⟳The sign

(- 1)^{a - 1}

counts the transpositions needed to move

θ_{j_{a}}

past

θ_{j_{1}}, \dots, θ_{j_{a - 1}}

; the hat denotes omission.

These satisfy the graded Leibniz rule and the nilpotency $(\partial / \partial θ_{i})^{2} = 0$ and anticommutation ${\partial_{θ_{i}}, \partial_{θ_{j}}} = 0$ , plus ${\partial_{θ_{i}}, θ_{j}} = δ_{ij}$ .

Berezin Integral

The Berezin integral $\int d θ_{i}$ over a single Grassmann variable is the unique (up to normalization) linear functional satisfying translation invariance $\int d θ f (θ + η) = \int d θ f (θ)$ . It is fixed by:

\int d θ_{i} 1 = 0, \int d θ_{i} θ_{i} = 1.

⟳Translation invariance forces

\int d θ 1 = 0

; the normalization

\int d θ θ = 1

is convention. Hence Berezin integration coincides with left differentiation:

\int d θ_{i} = \partial / \partial θ_{i}

acting from the left.

For multiple variables one orders the measures and integrates from right to left:

\int d^{n} θ θ_{1} θ_{2} \dots θ_{n} = 1, d^{n} θ := d θ_{n} d θ_{n - 1} \dots d θ_{1},

⟳Apply [4] iteratively,

\int d θ_{n} \dots \int d θ_{1}

, each acting on the leftmost remaining factor; the prescribed ordering gives

+ 1

with no extra sign.

Crucially, the measure carries odd parity: $d θ_{i}$ anticommutes with all Grassmann generators and with other measures, $d θ_{i} d θ_{j} = - d θ_{j} d θ_{i}$ .

Berezin Gaussian Integral

Let $A$ be an $n \times n$ matrix and ${θ_{i}, \overset{ˉ}{θ}_{i}}$ be $2 n$ independent Grassmann generators. Then

\int i = 1 \prod n d \overset{ˉ}{θ}_{i} d θ_{i} exp (- i, j \sum \overset{ˉ}{θ}_{i} A_{ij} θ_{j}) = det A .

⟳Expand the exponential — it terminates because each

θ_{i}, \overset{ˉ}{θ}_{i}

appears at most once. Only the top term

\frac{1}{n !} (- \overset{ˉ}{θ} A θ)^{n}

survives integration; collecting it gives

\sum_{σ} sgn (σ) \prod_{i} A_{iσ (i)} = det A

For a single set of $n$ Grassmann variables with antisymmetric $A = - A^{T}$ ,

\int d^{n} θ exp (\frac{1}{2} θ^{T} A θ) = Pf (A), Pf (A)^{2} = det A,

⟳The top term is

\frac{1}{2 ^{n} n !} (θ^{T} A θ)^{n}

; matching against

θ_{1} \dots θ_{2 n}

reproduces the Pfaffian's signed sum over perfect matchings.

Why It Holds / Derivation

Berezin Gaussian = det A

Step 1. Write the bilinear exponent and expand. Because

\overset{ˉ}{θ}_{i}^{2} = θ_{j}^{2} = 0

, the Taylor series of the exponential truncates:

exp (- \overset{ˉ}{θ} A θ) = m = 0 \sum n \frac{( - 1 ) ^{m}}{m !} (ij \sum \overset{ˉ}{θ}_{i} A_{ij} θ_{j})^{m} .

Step 2. The Berezin measure

\prod_{i} d \overset{ˉ}{θ}_{i} d θ_{i}

returns zero on any term not containing every

θ_{i}

and every

\overset{ˉ}{θ}_{i}

exactly once (by [4]). Only the

m = n

term qualifies. Step 3. Expand

(\sum_{ij} \overset{ˉ}{θ}_{i} A_{ij} θ_{j})^{n}

. Surviving monomials use each

\overset{ˉ}{θ}

and each

θ

once, so the

θ

-indices and

\overset{ˉ}{θ}

-indices each range over a permutation. The combinatorial factor

n!

cancels the

1/ n!

. Step 4. Reordering the Grassmann factors into the canonical order

\overset{ˉ}{θ}_{1} θ_{1} \dots \overset{ˉ}{θ}_{n} θ_{n}

produces the signature

sgn (σ)

of the permutation relating

θ

-indices to

\overset{ˉ}{θ}

-indices. The

(- 1)^{n}

from the prefactor combines with the measure ordering to leave the standard determinant sign. Step 5. Collecting,

\int = \sum_{σ \in S_{n}} sgn (σ) \prod_{i} A_{i, σ (i)} = det A .

■

Change of Variables / Jacobian Inversion

Let $θ_{i} = \sum_{j} M_{ij} θ_{j}^{'}$ be an invertible linear change of Grassmann variables. The top-form transforms as $θ_{1} \dots θ_{n} = (det M) θ_{1}^{'} \dots θ_{n}^{'}$ .

Step 1. Since integration extracts the coefficient of the top monomial,

\int d^{n} θ f =

(coeff of

θ_{1} \dots θ_{n}

). Step 2. For

\int d^{n} θ f = \int d^{n} θ^{'} f

to give consistent values, the measure must transform inversely to the coordinates:

d^{n} θ = (det M)^{- 1} d^{n} θ^{'} .

⟳Demand

\int d^{n} θ (θ_{1} \dots θ_{n}) = \int d^{n} θ^{'} (θ_{1}^{'} \dots θ_{n}^{'}) = 1

; substitute the top-form scaling by

det M

and require cancellation.

This is the opposite of the ordinary Jacobian rule ( $d^{n} x = ∣ det M ∣ d^{n} x^{'}$ ), and is why the Berezin–Gaussian integral yields $det A$ rather than $(det A)^{- 1}$ . $■$

Felix Alexandrovich Berezin | 1931–1980

A Soviet mathematical physicist at Moscow State University, Berezin systematized the calculus of anticommuting variables in his 1966 monograph The Method of Second Quantization. He introduced what is now called the Berezin integral, the supertrace, the superdeterminant (Berezinian), and laid foundations for the geometry of supermanifolds and geometric quantization. His insistence that integration over fermions be defined purely algebraically — abandoning any measure-theoretic picture — proved exactly right for quantum field theory. Berezin drowned in a kayaking accident on the Kolyma River in 1980; the integral, the Berezinian, and Berezin–Toeplitz quantization carry his name.

Building Intuition

Think of a Grassmann variable as a "coordinate that has only two states": absent or present, with no notion of magnitude (since $θ^{2} = 0$ kills any second power). A function of $θ$ is therefore linear: $f (θ) = a + b θ$ . There is nothing to integrate in the Riemann sense — there is no interval, no area under a curve. The only translation-invariant linear readout of $f$ is its slope $b$ . So " $\int d θ$ " is just "read off the coefficient of $θ$ ," which is the same as differentiating. Integration and differentiation being the same operation is the signature weirdness of Grassmann calculus. The exclusion principle ( $θ^{2} = 0$ : you cannot put two fermions in the same state) is built into the algebra, not imposed afterward.

Boson vs. Fermion Gaussian integration: same quadratic form, inverse powers of the determinant

x

)

Worked Examples

One-variable integral with explicit coefficients

Given:

f (θ) = 3 + 5 θ

, a function of a single Grassmann variable. Step 1. By linearity,

\int d θ f (θ) = 3 \int d θ 1 + 5 \int d θ θ

. Step 2. Apply rules [4]:

\int d θ 1 = 0

and

\int d θ θ = 1

. Step 3. Therefore

\int d θ f (θ) = 3 \cdot 0 + 5 \cdot 1 = 5

\int d θ (3 + 5 θ) = 5 (dimensionless)

Note this equals $f^{'} (θ)$ — integration picked off the coefficient of $θ$ , identical to differentiation.

$2\times 2$ Berezin Gaussian = determinant

Given: Two Grassmann pairs

(θ_{1}, θ_{2}, \overset{ˉ}{θ}_{1}, \overset{ˉ}{θ}_{2})

and matrix

A = (2314), S = \overset{ˉ}{θ}_{1} (2 θ_{1} + θ_{2}) + \overset{ˉ}{θ}_{2} (3 θ_{1} + 4 θ_{2}) .

Step 1. Compute

\int d \overset{ˉ}{θ}_{1} d θ_{1} d \overset{ˉ}{θ}_{2} d θ_{2} e^{- S}

. Expand

e^{- S} = 1 - S + \frac{1}{2} S^{2}

(higher terms vanish, only

4

generators). Step 2. Only the term with all four generators survives. In

\frac{1}{2} S^{2}

the surviving cross-terms are

\overset{ˉ}{θ}_{1} θ_{1} \overset{ˉ}{θ}_{2} θ_{2}

-type. Collecting signs gives

det A

. Step 3. Directly:

det A = (2) (4) - (1) (3) = 8 - 3 = 5

\int d \overset{ˉ}{θ}_{1} d θ_{1} d \overset{ˉ}{θ}_{2} d θ_{2} e^{- \overset{ˉ}{θ} A θ} = det A = 5 (dimensionless)

Contrast: the bosonic Gaussian with the symmetric part of $A$ would give $(2 π)^{1} / det A = 2 π / 5$ , the inverse power.

Common Misconceptions

Integration is NOT the inverse of differentiation here

Although $\int d θ = \partial / \partial θ$ as operators, integration does not "undo" differentiation in the analytic sense. There is no constant of integration and no antiderivative. Likewise, do not import the absolute-value Jacobian $∣ det M ∣$ from real analysis: the Grassmann measure transforms by $(det M)^{- 1}$ with sign (equation [8]), because the measure is itself an odd object that transforms inversely to the coordinates. Forgetting this inversion is the single most common error in fermionic path-integral calculations.

The "obvious" measure scaling is wrong

Naively one might guess that rescaling $θ \to λ θ$ scales the measure the same way, $d θ \to λ d θ$ , as in ordinary calculus. This is false. Demand consistency: $1 = \int d θ θ = \int d (λ θ) (λ θ)$ . If $d (λ θ) = λ^{- 1} d θ$ , then $λ^{- 1} \cdot λ \int d θ θ = 1$ , consistent. The measure scales as $λ^{- 1}$ , inversely to the variable. A student who writes $d (λ θ) = λ d θ$ will obtain $λ^{2}$ and a wrong answer — and in the multivariate case will get $det M$ instead of $(det M)^{- 1}$ , flipping the determinant to the bosonic side.

Connections

Grassmann calculus is the algebraic core of several neighboring structures. The exterior algebra it defines is exactly the home of Differential Forms, where $d θ$ plays the role of $d x^{i}$ and the wedge product realizes [1]. The Berezin–Gaussian formula [6] is the fermionic analogue of the Gaussian Integral and underlies the Path Integral Formulation of fermionic quantum fields. The superdeterminant (Berezinian) generalizes Determinant to graded matrices, and the whole apparatus is the local model for Supermanifolds and the foundation of Supersymmetry. The Pfaffian formula [7] connects directly to the combinatorics of perfect matchings and to Random Matrix Theory.

Real-World Applications

Quantum Field Theory of fermions — The QED/QCD partition function integrates out Dirac fields via $\int D \overset{ˉ}{ψ} D ψ e^{- \overset{ˉ}{ψ} (\neq D + m) ψ} = det (\neq D + m)$ , the fermion determinant driving vacuum polarization and the chiral anomaly.
Lattice QCD — Pseudofermion algorithms and the Hybrid Monte Carlo method handle the fermion determinant arising from Berezin integration; the sign problem stems directly from $det$ becoming complex.
Supersymmetry and the MSSM — Superspace formalism builds Lagrangians as Berezin integrals $\int d^{2} θ W (Φ)$ over Grassmann superspace coordinates.
Random Matrix Theory — The supersymmetric method (Efetov) replaces ensemble averages of $1/ det$ by Grassmann–bosonic integrals, computing level correlations in disordered metals and quantum chaos.
Anderson localization & mesoscopic transport — Efetov's nonlinear sigma model is derived entirely through Berezin integration over fermionic and bosonic replicas.
Topological field theory & BRST quantization — Faddeev–Popov ghosts are Grassmann fields; the gauge-fixing determinant is a Berezin integral.
Combinatorics of graphs — The Pfaffian/dimer model (Kasteleyn) for counting perfect matchings on planar lattices is literally a Berezin–Gaussian integral [7].

Historical Development

1844 — Hermann Grassmann publishes Die lineale Ausdehnungslehre, introducing the exterior (anticommutative) algebra of "extensive quantities," largely ignored for decades.
1888 — Giuseppe Peano and later Élie Cartan revive and formalize the exterior algebra and differential forms, giving $θ_{i} \land θ_{j} = - θ_{j} \land θ_{i}$ geometric meaning.
1925–1928 — The Pauli exclusion principle and Dirac/Jordan–Wigner anticommutation relations ${a_{i}, a_{j}^{†}} = δ_{ij}$ make anticommuting objects central to physics.
1933 — Eddington and others use Pfaffians; the determinant–Pfaffian relation $Pf (A)^{2} = det A$ becomes standard.
1966 — Felix Berezin's The Method of Second Quantization defines integration over anticommuting variables, equations [4]–[6], inventing the modern Berezin integral.
1971–1975 — Salam–Strathdee, Wess–Zumino, and Akulov–Volkov develop superspace; Grassmann coordinates become physical superspace dimensions.
1975–1977 — Faddeev, Popov, and others recast gauge-theory quantization with Grassmann ghost fields; Berezin integration becomes routine in QFT.
1979–1983 — Efetov's supersymmetric method applies Berezin integration to disordered electronic systems and random matrices.
1987 — DeWitt's Supermanifolds and Rogers' work give a rigorous geometric framework, embedding Berezin integration in the theory of $Z_{2}$ -graded manifolds.

Subtleties & Limits

Infinite-dimensional rigor. For finitely many generators, Berezin integration is purely algebraic and unproblematic. In quantum field theory one needs infinitely many Grassmann variables, and the resulting fermionic functional integral

\int D \overset{ˉ}{ψ} D ψ

has no underlying measure-theoretic foundation comparable to the (already delicate) bosonic case; it is defined by its formal expansion and the determinant identity, with regularization required to make

det (\neq D + m)

finite. Curved Grassmann measures. On a general supermanifold, defining a global, invariant Berezin volume element requires the Berezinian and a choice of "integral forms"; the theory of integration on supermanifolds (Voronov, Witten's recent work on superstring measures and integral forms) involves genuine subtleties absent in ordinary geometry — e.g., the failure of naive integration-by-parts and the role of the "picture-changing" structure. The sign problem. When

det A

is complex or negative, Monte Carlo methods built on Berezin integration lose their probabilistic interpretation; whether a general efficient algorithm exists is a deep open question linked to computational complexity.