Independent Unitary Pairing

• Independent Unitary Pairing is a phenomenon in fermionic quantum systems defined by generalized Pauli constraints that form convex occupation-number polytopes.
• It introduces additional linear inequalities beyond classical Pauli bounds, which simplify the wavefunction expansion and impact correlation energy computations.
• Experimental models like harmonium and Hubbard demonstrate how pinning and quasipinning provide actionable insights into state selection and variational method improvements.

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 “^^^^2^^^^”) 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 (Tennie et al., 2016, Schilling, 2015).

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) (Tennie et al., 2016, Schilling et al., 2012, Schilling, 2015).

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$) (Tennie et al., 2016, Schilling et al., 2012). 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 (Tennie et al., 2016).

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 (Benavides-Riveros et al., 2017, Schilling et al., 2012).

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 (Schilling et al., 2012).

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 (Benavides-Riveros et al., 2017). 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}$ (Liebing et al., 2014).

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<U_c$
Composite Fermion	Yes (via pairing)	Mott-bound phase densities

7. Experimental Probes, Extensions, and Theoretical Outlook

Recent quantum-gas microscopy and quantum-dot platforms can resolve occupations and measure pinning/nonpinning transitions by tracking one-body spectra under variable interaction strengths (Hackl et al., 2021). Preparation of highly entangled fermionic pure states permits direct access to GPC facets, with state tomography and population measurements reconstructing $\vec\lambda$ for each system. High global purity is essential; otherwise, GPCs relax and only traditional Pauli bounds apply.

Extensions to general Hubbard, Kondo, and systems with longer-range interactions maintain the rigorous polytope bounds, with further momentum-occupation-number inequalities applicable in the thermodynamic limit (Lapa, 2021). Open questions include the constructive characterization of all GPC facets for higher $(N,d)$, and the development of algorithms exploiting unitary pairing-induced simplification for quantum simulation and reduced-density-matrix functional theories.

In sum, independent unitary pairing, as embodied in GPC-induced occupancy bounds and wave function structure, provides a rigorous framework for analyzing antisymmetric quantum matter beyond traditional exclusion-principle physics, with deep consequences for theory, simulation, and experiment.