Generalized Pauli Principle (GPP)

Updated 15 March 2026

Generalized Pauli Principle is a framework of additional linear constraints on the natural occupation numbers that go beyond the standard Pauli exclusion principle.
It defines a convex polytope, using facets and geometric structures (e.g., the Borland–Dennis case) to precisely characterize allowed fermionic states.
GPP has significant practical implications, such as inducing pinning/quasipinning phenomena that simplify many-body wave function representations and computational approaches.

The Generalized Pauli Principle (GPP) formalizes and quantifies all constraints on the occupation numbers of the one-particle reduced density matrix (1-RDM) for a pure N-fermion state beyond the standard Pauli exclusion principle, expressing the full content of fermionic exchange symmetry. These constraints—known as Generalized Pauli Constraints (GPCs)—impose a set of homogeneous linear inequalities that sharply restrict the allowed region of natural occupation numbers and underpin new structural, computational, and physical phenomena in quantum many-body systems.

1. Standard and Generalized Pauli Constraints

The traditional Pauli exclusion principle asserts that, for the eigenvalues $\{\lambda_i\}$ of a 1-RDM of an $N$ -fermion pure state (the natural occupation numbers, or NONs), one has

$0 \le \lambda_i \le 1, \quad \sum_{i=1}^d \lambda_i = N.$

This defines the "Pauli simplex" $\Sigma = \{\vec{\lambda}\in \mathbb{R}^d\,|\, \sum_i \lambda_i = N,\, 1 \ge \lambda_1 \ge \cdots \ge \lambda_d \ge 0\}$ . However, the antisymmetry requirement of fermions leads to strictly stronger, additional linear inequalities—the GPCs—so that the set of allowed spectra is a convex polytope $\mathcal{P}\subset\Sigma$ cut out by finitely many constraints

$D_j(\vec{\lambda}) \equiv \kappa_j^{(0)} + \sum_{i=1}^d \kappa_{j,i}\lambda_i \ge 0, \quad j=1,\ldots, r_{N,d},$

with integer coefficients and minimal representation, as established by Klyachko and Altunbulak (Tennie et al., 2016, Smart et al., 2020, Reuvers, 2019). Only $\vec{\lambda}$ lying in $\mathcal{P}$ are compatible with pure $N$ -fermion states.

For $N=3$ , $d=6$ (the Borland–Dennis case), the explicit constraints are:

Three pair-sum equalities: $\lambda_1+\lambda_6 = \lambda_2+\lambda_5 = \lambda_3+\lambda_4 = 1$ ,
One genuine GPC: $2 - (\lambda_1+\lambda_2+\lambda_4) \ge 0$ .

2. Mathematical Derivation and Geometric Structure

GPCs arise from the pure-state quantum marginal (N-representability) problem: determining which 1-RDM spectra can be realized by pure, antisymmetric $N$ -body states. The solution, via Schur–Weyl duality and invariant theory, characterizes permissible $\vec{\lambda}$ by the intersection of the hyperplane $\sum_i \lambda_i=N$ , the descending ordering, and all GPCs, yielding the polytope $\mathcal{P}$ .

Each GPC $D_j$ defines a facet $F_j$ of $\mathcal{P}$ . Geometrically, $\mathcal{P}$ is a high-dimensional convex polytope, embedded in $\Sigma$ , with facets corresponding to both standard exclusion and deeper antisymmetry constraints (Tennie et al., 2016, Schilling, 2015).

For $N=2$ and all $d$ , the set is fully characterized by paired eigenvalues (Yang–Youla, no new inequalities). For $N=3, d=6$ , the Borland–Dennis polytope is three-dimensional with the above constraints, reflecting the first non-trivial GPP structure (Reuvers, 2019).

3. Physical Consequences: Pinning, Quasipinning, and Wave Function Simplification

Pinning occurs if the NON vector $\vec{\lambda}$ lies exactly on a GPC facet: $D_j(\vec{\lambda})=0$ . More generally, quasipinning refers to $\vec{\lambda}$ being extremely close to a facet, $D_j(\vec{\lambda}) \ll 1$ but nonzero. The minimal GPC distance,

$D_{\min} \equiv \min_j D_j(\vec{\lambda}),$

is a canonical measure of quasipinning.

Pinning triggers a superselection rule: the $N$ -fermion wave function must be a linear combination of only those Slater determinants whose occupation patterns also saturate the pinned GPC(s), drastically reducing the effective number of determinants in the expansion (Maciążek et al., 2019). For the Borland–Dennis case, pinning to the nontrivial facet restricts the wave function to a four-dimensional subspace of the 20 possible determinants.

The physical mechanism driving quasipinning is the competition between energy minimization (favoring occupation spread near the Fermi level) and exchange antisymmetry (forbidding repeated occupation), resulting in high but non-exact saturation of GPCs in highly correlated, low-dimensional states (Tennie et al., 2016).

4. Quantification: The Q-Parameter and Trivial vs Non-Trivial Quasipinning

To distinguish true GPP-induced quasipinning from situations where closeness to a facet simply reflects nearly integer NONs (close to the Pauli simplex boundary), one introduces the $Q$ -parameter:

$Q_j(\vec{\lambda}) = -\log_{10}\left[\frac{D_j(\vec{\lambda})}{c_j\,\mathrm{dist}_1(\vec{\lambda},\Sigma_{r,s})}\right],$

where $c_j$ is the optimal constant bounding $D_j$ in terms of the distance to the nearest Pauli facet $\Sigma_{r,s}$ . Large $Q_j$ indicates $\vec{\lambda}$ is much closer to a GPC facet than standard exclusion would require: truly non-trivial "quasipinning" (Tennie et al., 2015, Tennie et al., 2016).

This distinction is decisive in applications: in large, high-dimensional systems ( $d \gg N$ ), GPCs remove only an exponentially small fraction of the allowed volume of $\Sigma$ ; in low-dimensional (few-orbital) cases, GPC effects are prominent and non-trivial quasipinning can be detected and quantified (Reuvers, 2019). Empirically, this has significant consequences for static/dynamic correlation diagnostics and for efficient many-body ansätze (Tennie et al., 2015).

5. Physical Realizations and Empirical Tests

GPP and its consequences have been rigorously explored in exactly and numerically solvable models, including Harmonic models (Harmonium), small-lattice Hubbard models, and small atoms (e.g., Li, Be) (Tennie et al., 2016, Schilling et al., 2017, Tennie et al., 2016). Consistent findings include:

Weak-coupling scaling laws for quasipinning, e.g., $D_{\min}\sim\kappa^8$ for N=3, d=6 Harmonium,
Absence of exact pinning but strong non-trivial quasipinning in small atomic systems,
Drastic wave function simplification in pinning or strong quasipinning regimes.

Experimental tests have directly verified the GPP using quantum computing platforms. Random pure states of three fermions in six orbitals were prepared and tomographically analyzed, confirming the Borland–Dennis constraint to within one part in $10^{18}$ (Smart et al., 2020). The absence of any violations, compared to the high violation probability under mere Pauli exclusion, provides conclusive evidence for the physical reality of GPP constraints.

6. Extensions: Spin, Reduced Symmetry, and Generalized Statistics

Spin degrees of freedom and spatial dimension further partition the GPP polytope. Spin-adapted GPCs define stricter domains for sectors of fixed spin quantum numbers, and pinning/quasipinning effects typically weaken with reduced symmetry or increased dimensionality (Liebert et al., 21 Feb 2025, Tennie et al., 2016). For small, spin-polarized systems, spin-adapted GPCs can be precisely computed and yield stronger superselection rules for wave function structure.

Beyond standard fermions, the GPP concept extends to systems with partial distinguishability, exclusion statistics, and generalized particle statistics. In the context of "WH Statistics" (Hao et al., 19 Jan 2026), a unified exclusion rule

$n_i \leq g_i/\kappa$

interpolates between Fermi–Dirac, Bose–Einstein, anyons, and classical hard-core limits, providing a physical foundation for GPP across indistinguishability and interaction regimes. In systems with multiple species or internal structure, highly nontrivial combinatorial GPPs emerge, governing mutual exclusion and occupation constraints (Liu et al., 2011).

7. Computational and Methodological Frameworks

The presence of GPCs has profound impact on computational methods:

Active space selection in configuration interaction (CI) or multiconfigurational self-consistent field (MCSCF) schemes: Pinning/quasipinning restricts the relevant determinants/sectors for variational optimization (Maciążek et al., 2019, Theophilou et al., 2017).
Reduced-density-matrix functional theory (RDMFT): GPCs define the true boundary of the physical domain for functionals and introduce “exchange-force” divergences at the polytope boundary (Liebert et al., 21 Feb 2025).
Efficient large- $N$ expansions and group-theoretical treatments: By encoding GPP at the level of normal-mode quanta, one implements antisymmetry non-projectively for arbitrary $N$ , enabling scalable perturbation approaches in, e.g., the unitary Fermi gas (Watson, 2015).

Quasipinning analysis relies on precise computation of occupation numbers (truncating near-filled or near-empty orbitals), computation of all relevant $D_j$ , and assessment of $Q$ to identify nontrivial GPC effects. These tools, initially feasible only for $N\leq5, d\leq11$ , have been extended to larger (spin-adapted) systems using symmetry and group theory (Liebert et al., 21 Feb 2025).

The Generalized Pauli Principle, by systematically extending Pauli’s original exclusion rule, provides a mathematically rigorous, physically verified, and computationally indispensable framework for understanding correlation, symmetry, and structural simplicity in few- and many-fermion quantum systems. Its constraints manifest in both fundamental structural selection rules and emergent properties of correlated quantum matter (Tennie et al., 2016, Schilling, 2015, Reuvers, 2019, Tennie et al., 2016, Smart et al., 2020, Maciążek et al., 2019, Liebert et al., 21 Feb 2025, Hao et al., 19 Jan 2026, Theophilou et al., 2017, Schilling et al., 2017, Chakraborty et al., 2014).