--- title: Independent Unitary Pairing url: https://www.emergentmind.com/topics/independent-unitary-pairing type: topic --- # Independent Unitary Pairing Independent Unitary Pairing describes a set of occupation-number bounds and structural constraints for fermionic quantum systems, arising from the fundamental antisymmetry of multi-particle wavefunctions. In contrast to the classical Pauli exclusion principle, which restricts individual occupation numbers between 0 and 1, the phenomenon encompasses additional linear inequalities—known as generalized Pauli constraints (GPCs)—that nontrivially correlate the occupations of different orbitals. These constraints manifest geometrically as the definition of a convex polytope in occupation-number space, whose facets correspond to the saturation of specific linear combinations. Saturation (or near-saturation—“pinning” or “quasipinning”) of GPCs produces selection rules for the many-body wave function, dramatically simplifying its Slater determinant expansion and modifying physical observables, correlation energies, and response properties. The operational relevance is substantiated in models such as harmonically trapped fermions (Harmonium), the Hubbard model, and variants with explicit symmetry, with both theoretical and experimental probes possible. ## 1. Generalized Pauli Constraints and Occupation-Number Polytopes For a pure state of $N$ fermions in a $d$-dimensional single-particle Hilbert space, the natural occupation numbers (NONs) $\vec\lambda = (\lambda_1 \geq \lambda_2 \geq ... \geq \lambda_d \geq 0)$—eigenvalues of the one-body reduced density operator—must obey \[ 0 \leq \lambda_i \leq 1 \,,\quad \sum_{i=1}^d \lambda_i = N \] (classical Pauli bounds). Klyachko and collaborators established that complete antisymmetry of the wavefunction induces further, finitely many linear inequalities (GPCs), represented canonically as \[ D_j(\lambda) \equiv \kappa_j^{(0)} + \sum_{i=1}^d \kappa_j^{(i)}\lambda_i \geq 0 \,,\quad j = 1, ...\, r_{N,d} \] with integer coefficients. The intersection of these inequalities defines a convex polytope $\mathcal{P} \subset \mathbb{R}^d$, the set of all physical NON vectors for pure fermionic states [1602.05198][1507.00299]. For instance, in the $(N,d) = (3,6)$ Borland–Dennis case: \[ \lambda_1 + \lambda_6 = \lambda_2 + \lambda_5 = \lambda_3 + \lambda_4 = 1 \] with the single nontrivial GPC: \[ D^{(3,6)}(\lambda) = 2 - (\lambda_1 + \lambda_2 + \lambda_4) \geq 0. \] ## 2. Geometric Interpretation: Polytope Facets and Pinning The GPCs define proper facets ("faces") within the Pauli simplex $\Sigma = \{\vec\lambda \mid 1 \geq \lambda_1 \geq ... \geq \lambda_d \geq 0,\, \sum_i\lambda_i = N\}$. The physical occupation numbers $\vec\lambda$ are confined to the strictly smaller region $\mathcal{P} \subset \Sigma$. When $\vec\lambda$ falls exactly on a facet $F_{D_j}$ (i.e., $D_j(\lambda) = 0$), this state is "pinned" and the wave function exhibits additional symmetry; quantitative proximity is measured via $D_{\min}(\lambda) = \min_{j} D_j(\lambda)$ ("quasipinning" when $D_{\min}$ is small but not zero) [1602.05198][1210.5531][1506.08833]. In harmonium and small-site Hubbard models, pinning and quasipinning structure transitions can be induced and tracked by scanning system parameters (e.g., interaction strength). These transitions have direct consequences for wavefunction entanglement and the active configuration space. ## 3. Quasipinning Metrics, Scaling Laws, and Physical Content The degree of quasipinning is explicitly characterized: for three spinless fermions in one-dimensional harmonium, $D_{\min}(\kappa) \sim c \cdot \kappa^8$ in the weak-coupling regime, where $\kappa$ is the dimensionless coupling ($\kappa \to 0$) [1602.05198][1210.5531]. For $N$ particles ($N \leq 8$), the quasipinning exponent becomes $2N$, i.e., \[ D_{\min}(\kappa) \sim d_N\,\kappa^{2N} \] These exponents persist over intermediate couplings before vanishing for strong interactions. Comparison to the distance from the boundary of the Pauli simplex yields the Q-parameter, which quantifies the "strength" of GPC-induced quasipinning relative to Pauli constraints: \[ Q_j(\lambda)= \log_{10} \left[\frac{\textrm{dist}_1(\lambda, \text{facet of }\Sigma)}{\textrm{dist}_1(\lambda, F_{D_j})}\right] \] For modest couplings, $Q(\kappa) \sim 2|\log_{10}\kappa|$, signifying a major enhancement [1602.05198]. ## 4. Selection Rules and Simplification of the Wave Function Pinning leads to strict selection rules for the expansion of the many-body wavefunction over Slater determinants. Specifically, only those configurations $|i_1,\ldots,i_N\rangle$ which reside in the zero-eigenspace of the corresponding operator $\hat{D}_j$ are permitted: \[ D_j(\lambda) = 0 \implies \hat{D}_j|\Psi\rangle = 0 \] This truncates the active CI (Configuration Interaction) space, typically to include only a handful of determinants. For the Borland–Dennis polytope, pinning enforces \[ |\Psi_{\textrm{BD}}\rangle = \alpha\,|1,2,3\rangle + \beta\,|1,4,5\rangle + \gamma\,|2,4,6\rangle \] with all other determinants excluded [1711.09129][1210.5531]. When only quasipinning occurs, stability results guarantee the wavefunction lies close (in Hilbert space norm) to the pinned subspace, making the construction robust against small deviations [1210.5531]. ## 5. Implications for Variational Methods and Correlation Energy The structural simplification enables a "natural extension of Hartree–Fock" ("NEHF", *Editor's term*) that variationally optimizes over only the permitted configurations selected by the saturated GPC. This approach gives universal geometric bounds on the missing correlation energy for pinned/quasipinned states \[ \Delta E \leq C\, D_j(\vec\lambda) \] \[ \frac{\Delta E}{E_{\rm corr}} \leq K \frac{D_j(\lambda)}{S(\lambda)} \] where $S(\lambda)$ is the $l^1$ distance from the Hartree–Fock corner, and $C,K$ depend on the spectral structure of the Hamiltonian [1711.09129]. In the lithium atom example, three Slater determinants recover $\sim$87% of the correlation energy, with more than 95% possible upon basis refinement. ## 6. Occupancy Bounds in Composite and Correlated Systems In models of composite fermion formation—e.g., the extended Lipkin and Faddeev approaches for two- and three-body bound states—the algebraic structure of the pairing operator and resulting in-medium occupation introduces further fermion-like bounds. Pauli blocking constrains the maximal density of two-fermion composites ("bosons") to the "Mott density" $n_2^{\rm Mott} = 1/\alpha$ (with form-factor normalization $\alpha$), while Bose enhancement in three-fermion compounds partially offsets Pauli blocking, yielding a higher critical Mott density $n_3^{\rm Mott} = 1/\beta > n_2^{\rm Mott}$ [1406.0396]. This physical interplay yields Borromean binding phenomena: three-body bound states may persist even where all two-body substructures are dissolved, governed by unitary pairing and exclusion principles. | System | GPC polytope non-trivial? | Pinning regime(s) | |----------------------|--------------------------|-----------------------------| | 1D Harmonium (N=3) | Yes | Quasipinning ($\kappa^8$ law) | | Hubbard (few sites) | Yes | Exact pinning for $U