---
title: 'Pauli Periodicity: Spectral and Atomic Perspectives'
url: https://www.emergentmind.com/topics/pauli-periodicity
type: topic
---

# Pauli Periodicity: Spectral and Atomic Perspectives

Searching arXiv for the cited papers to ground the article in the source literature.
“Pauli periodicity” is not a single standardized term; rather, it plausibly denotes several technically distinct periodicity phenomena whose formulation is governed by Pauli structure. In one usage, it refers to Borg-type periodicity theorems for first-order self-adjoint systems in which the Pauli matrices provide the natural decomposition of a Hermitian $\pi$-periodic potential and sharply characterize when half-periodicity, $\sigma_2$-symmetry, or gap closure occurs [1707.00982]. In a second usage, it refers to the shell periodicity of atoms obtained from a Pauli-exclusion “potential” formulated as a classical excluded-volume interaction in a four-dimensional thermal-space [2211.15764]. In a third, more polemical context, it concerns alleged “Pauli crystals,” where apparent periodic patterns in aligned single-shot configurations are argued to be artifacts, while exclusion instead produces anti-correlation structures termed “Pauli anti-crystals” [1903.02067]. Taken together, these uses show that “periodicity” enters through spectral multiplicity, shell organization, or correlation geometry, depending on the underlying formalism.

## 1. Pauli-matrix periodicity in first-order self-adjoint systems

A precise spectral-theoretic instance of Pauli periodicity arises for the first-order self-adjoint Dirac-type system
$$
J Y'(z,\lambda) + Q(z) Y(z,\lambda) = \lambda Y(z,\lambda), \qquad 0\le z\le \pi,
$$
with
$$
J=\begin{pmatrix}0&1\\-1&0\end{pmatrix},
$$
where $Q(z)$ is Hermitian, $\pi$-periodic, and integrable on compact sets, and $Y(z,\lambda)$ is a fundamental $2\times2$ matrix solution satisfying $Y(0,\lambda)=I$ [1707.00982]. The Pauli matrices
$$
\sigma_0=I,\qquad
\sigma_1=\begin{pmatrix}0&1\\1&0\end{pmatrix},\qquad
\sigma_2=\begin{pmatrix}0&-i\\ i&0\end{pmatrix},\qquad
\sigma_3=\begin{pmatrix}1&0\\0&-1\end{pmatrix}
$$
play a double role: they form an orthonormal basis of $M_2(\mathbb C)$ and encode the natural $J$-decomposition of $Q$.

The potential is written as
$$
Q(z)=Q_1(z)+Q_2(z),
$$
where $Q_1$ belongs to the subspace skew-commuting with $J$ and $Q_2$ belongs to the subspace commuting with $J$ [1707.00982]. In components,
$$
Q_2(z)=(q_1(z)+q_2(z))\sigma_0+\operatorname{Im}q(z)\,\sigma_2.
$$
Because $[J,Q_2]=0$, one obtains the explicit gauge transformation
$$
R(z)=\exp\!\Bigl(J\!\int_0^z Q_2(t)\,dt\Bigr),
$$
which conjugates the original system into one with new potential
$$
\widetilde Q(z)=R(z)^{-1}Q_1(z)R(z)\in S_J,
$$
while preserving the same boundary conditions at $z=0,\pi$ [1707.00982].

This Pauli-based decomposition is the structural origin of the periodicity theorems. The periodicity properties are not imposed directly on the original $Q$; rather, they emerge after passing to the gauge-equivalent “canonical” system in which the commuting part has been removed. This suggests that, in this setting, Pauli periodicity is fundamentally a statement about how $\sigma_0$ and $\sigma_2$ organize the gauge-invariant content of the potential.

## 2. Discriminant, Floquet multipliers, and double-zero criteria

The monodromy matrix is $Y(\pi,\lambda)$, and the discriminant is defined by
$$
\Delta(\lambda)=\operatorname{Tr}Y(\pi,\lambda)=y_{11}(\pi,\lambda)+y_{22}(\pi,\lambda).
$$
Since $\det Y\equiv1$, the Floquet multipliers $p$ satisfy
$$
p^2-\Delta(\lambda)p+1=0
$$
[1707.00982]. The commuting part $Q_2$ contributes the overall factor
$$
e^{-i\beta},\qquad \beta:=\int_0^\pi \operatorname{Im}q(t)\,dt,
$$
so that the discriminant of the conjugated system is
$$
\widetilde\Delta(\lambda)=e^{i\beta}\Delta(\lambda),
$$
and the relevant entire functions are $\widetilde\Delta(\lambda)\pm2$.

Theorem A identifies half-periodicity through the zero structure of $\Delta(\lambda)+2e^{-i\beta}$: the zeros of
$$
\Delta(\lambda)+2e^{-i\beta}
$$
are all double if and only if the original system is unitarily equivalent, via $R(z)$, to one whose Hermitian potential $\widetilde Q(z)$ is $\pi/2$-periodic almost everywhere [1707.00982]. In the forward direction, $\pi/2$-periodicity implies that $Y(z+\pi)$ and $Y(z+\pi/2)$ satisfy the same Dirac equation and differ by the constant matrix $Y(\pi/2)$; setting $z=0$ gives
$$
Y(\pi)=Y(\pi/2)^2.
$$
In the reverse direction, assuming that $\widetilde\Delta(\lambda)+2$ has only double zeros, a Hadamard factorization argument together with high-energy asymptotics shows that
$$
F(z,\lambda)=Y(z+\pi)+Y(z)
$$
vanishes identically, and the same matrix identity $Y(\pi)=Y(\pi/2)^2$ then forces $\widetilde Q(z+\pi/2)=\widetilde Q(z)$ [1707.00982].

Theorem B gives the complementary $\sigma_2$-symmetry criterion. The zeros of
$$
\Delta(\lambda)-2e^{-i\beta}
$$
are all double if and only if the system is unitarily equivalent, via the same $R(z)$, to one satisfying
$$
\widetilde Q(z)=\sigma_2\,\widetilde Q(z)\,\sigma_2
$$
[1707.00982]. If $\widetilde Q$ is $\sigma_2$-symmetric, then $Y(z+\pi)$ and $\sigma_2Y(z+\pi/2)$ solve the same Dirac equation and differ by the constant $\sigma_2Y(\pi/2)$, leading to
$$
Y(\pi)=(\sigma_2Y(\pi/2))^2.
$$
Conversely, the double-zero structure of $\widetilde\Delta(\lambda)-2$ forces the same matrix equality and thereby the $\sigma_2$-symmetry of $\widetilde Q$ [1707.00982].

These theorems constitute the core Pauli periodicity mechanism in the spectral setting. Half-periodicity is encoded through $\widetilde\Delta+2$, whereas $\sigma_2$-conjugation symmetry is encoded through $\widetilde\Delta-2$.

## 3. Vanishing instability intervals and the canonical $\sigma_0,\sigma_2$ form

The same formalism yields an Ambarzumyan-type statement for Dirac systems. All spectral instability intervals, defined as the open gaps in $\{|\Delta(\lambda)|\le 2\}$, vanish if and only if, after the $R$-gauge, the potential has the form
$$
\widetilde Q(z)=r(z)\sigma_0+q(z)\sigma_2,
$$
with $r$ and $q$ real and $\pi$-periodic [1707.00982]. Equivalently, in the formulation given in the abstract, all instability intervals vanish if and only if
$$
Q=r\sigma_0+q\sigma_2
$$
for some real-valued $\pi$-periodic functions $r$ and $q$ integrable on compact sets.

The proof combines Theorems A and B. Vanishing of all gaps means that for every real $\lambda$ both $\widetilde\Delta(\lambda)+2$ and $\widetilde\Delta(\lambda)-2$ have only double zeros. This forces $\widetilde Q$ to be simultaneously $\pi/2$-periodic and $\sigma_2$-symmetric. Periodicity together with $\sigma_2$-anti-periodicity on the skew part $Q_1$ implies $Q_1\equiv0$, so the potential reduces to the commuting $\sigma_0,\sigma_2$ sector [1707.00982].

In that case the system can be solved explicitly:
$$
Y(z,\lambda)=\exp\bigl(J\lambda z-i\sigma_2\int_0^z q(t)\,dt-iI\int_0^z r(t)\,dt\bigr),
$$
and
$$
\Delta(\lambda)=2\cos\!\Bigl(\lambda\pi-\int_0^\pi q\Bigr)e^{-i\int_0^\pi r}.
$$
Hence $|\Delta(\lambda)|\le2$ for all real $\lambda$, so no gaps occur [1707.00982].

In this sense, Pauli periodicity is not merely a descriptive label but a complete spectral classification. The joint action of $\sigma_0$ and $\sigma_2$ eliminates the instability intervals by restricting the potential to the commuting subspace.

## 4. Shell periodicity from a Pauli-exclusion potential

A distinct use of Pauli periodicity appears in atomic structure, where the periodicity in question is the shell structure of the periodic table. In the ring-polymer representation, each quantum particle is modeled as a Gaussian thread parameterized by an imaginary-time coordinate $s\in[0,\beta]$, and a Pauli potential is postulated as a Flory-Huggins-type mean-field excluded-volume penalty between the four-dimensional contours of different electrons [2211.15764]. With pair densities $n_i(\mathbf r,\beta)$, the total Pauli energy is
$$
U_P[\{n\}]
=\frac12\sum_{i\ne j}\iint d\mathbf r\,d\mathbf r'\,
n_i(\mathbf r)\,V_{\rm xv}(\mathbf r-\mathbf r')\,n_j(\mathbf r'),
$$
with
$$
V_{\rm xv}(\mathbf r-\mathbf r')=g_0^{-1}\delta(\mathbf r-\mathbf r').
$$
This yields the local form
$$
U_P[\{n\}]
=\frac12\sum_{i,j}g_{ij}^{-1}\int d\mathbf r\,n_i(\mathbf r)\,n_j(\mathbf r),
$$
where $g_{ij}^{-1}=(1-\delta_{ij})g_0^{-1}$.

The uniform-limit scaling argument connects this Pauli potential to Thomas-Fermi theory. In mean field, $U_P/V\to \frac{g_0^{-1}}{2}n_0^2$, which has the wrong density scaling compared to the Thomas-Fermi kinetic energy density
$$
T_{\rm TF}[n]=C_{\rm TF}\int n^{5/3}\,d\mathbf r,\qquad
C_{\rm TF}=\tfrac{3}{10}(3\pi^2)^{2/3}.
$$
A polymer-scaling argument instead gives
$$
\frac{U_P}{N}\sim D^{-2}\sim n_0^{2/3}
\quad\Longrightarrow\quad
\frac{U_P}{V}\sim n_0^{5/3},
$$
in precise agreement with Thomas-Fermi scaling, and matching prefactors fixes
$$
g_0^{-1}=2C_{\rm TF}=\frac{3}{5}(3\pi^2)^{2/3}\approx5.7425
$$
with no further adjustable parameters [2211.15764].

The self-consistent-field equations follow from the free-energy functional
$$
F[\{n\},\{w\}]
= -\frac1\beta\sum_i N_i\ln Q_i
-\sum_i\int w_i\,n_i
+\bigl[U_{\rm ext}[\rho]+U_{\rm ee}[\rho]\bigr]
+U_P[\{n\}],
$$
leading to the mean-field potentials
$$
w_i(\mathbf r)=w_{0,i}(\mathbf r)+w_{P,i}(\mathbf r),
\qquad
w_{P,i}(\mathbf r)=g_0^{-1}\sum_{j\ne i}n_j(\mathbf r),
$$
where the electron-electron Hartree term carries a Fermi-Amaldi factor $(N-1)/N$ to reduce self-interaction [2211.15764]. The auxiliary propagators satisfy
$$
\frac{\partial q_i}{\partial s}
=\frac{\hbar^2}{2m}\nabla^2 q_i-w_i(\mathbf r)\,q_i,
\qquad
q_i(\mathbf r_0,\mathbf r,0)=V\,\delta(\mathbf r-\mathbf r_0),
$$
with
$$
Q_i=\frac1V\int d\mathbf r\,q_i(\mathbf r,\mathbf r,\beta),
\qquad
n_i(\mathbf r)=\frac{N_i}{Q_i}q_i(\mathbf r,\mathbf r,\beta).
$$
These equations are iterated to self-consistency.

Under spherical averaging, each $n_i(r)$ and $w_i(r)$ is expanded in spherical Bessel functions $j_0(\alpha_k r/R)$, and the total radial density is reconstructed as
$$
\rho(r)=\sum_i n_i(r)=\sum_l \rho_l(r),
$$
where $l$ labels grouped shell or sub-shell pairs such as $1s$, $2s$, and $2p$ [2211.15764]. Numerical calculations at $\beta=80$ a.u. for H through Ar show that radial electron densities display the correct number of peaks (“shells”) for each element. The errors on the total binding energies compared to known binding energies are less than $9\%$ for the lightest elements and drop to $3\%$ or less for atoms heavier than nitrogen; beyond $Z\approx7$ the binding-energy error stabilizes around $2$–$3\%$ [2211.15764].

Here the periodicity is chemical rather than Floquet-theoretic. The periodic table’s shell organization is traced to an effective Pauli potential, and the paper explicitly states that “the periodicity of chemical elements emerges naturally from classical-statistical correlations (excluded volume) among non-local Gaussian threads in an extra (imaginary-time) dimension” [2211.15764].

## 5. Claimed Pauli crystals, anti-crystals, and periodic anti-correlation patterns

A third context concerns whether Pauli exclusion alone can generate crystalline periodicity in single-shot many-fermion configurations. The configuration density (CD) method aligns Monte Carlo–sampled $N$-particle snapshots to a chosen seed configuration $\mathbf r'$ by minimizing the configuration distance
$$
g(\mathbf r|\mathbf r')
=\sum_{i,j=1}^N \exp\!\Bigl[-\frac1{G^2}|r_i-r'_j|^2\Bigr],
$$
and then averages the aligned one-body densities:
$$
\rho_{\rm CD}(x)=\frac1M\sum_{m=1}^M\sum_{i=1}^N
\delta^{(2)}\!\bigl(x-[R_m r_i^{(m)}]\bigr).
$$
Peaks in $\rho_{\rm CD}$ had been taken as evidence of “Pauli crystals,” but the analysis in [1903.02067] argues that such peaks can be imposed by the method itself.

The crucial numerical observation is that the same CD procedure applied to non-interacting fermions in a two-dimensional trap, Laughlin $q=1$ and $q=3$ states, and uniform random points produces the same seed-configuration pattern. This is presented as proof that the peaks can be artifacts of the alignment procedure rather than evidence of true underlying order [1903.02067]. The configuration variance
$$
C_V(r_\alpha)=\langle n^2\rangle_{r_\alpha}-\langle n\rangle_{r_\alpha}^2
$$
further shows that the CD peaks do not imply rigid particle binding to seed sites: for random points $C_V\approx1$, for non-interacting fermions with $q=1$ one finds $C_V\approx0.8$–$0.9$, and for a strongly correlated liquid with $q=3$ it drops further but does not vanish [1903.02067].

The alternative structure proposed in the paper is the “Pauli anti-crystal.” For a single Slater determinant with occupied orbitals $\{\psi_k(x)\}_{k\in K}$, the one- and two-body densities are
$$
\rho(x)=\sum_{k\in K}|\psi_k(x)|^2,
$$
and
$$
\rho(x,y)=\rho(x)\rho(y)-|\rho(x\to y)|^2+\delta(x-y)\rho(x),
\qquad
\rho(x\to y)=\sum_{k\in K}\psi_k(x)\psi_k^*(y).
$$
For flat one-body density $\rho(x)=\rho_0$, the reduced pair correlator is
$$
g(r)=\frac1N\int dR\,\rho(R+r,R)=\rho_0-f(r)+\delta(r),
$$
with
$$
f(r)=\frac1N\int dR\,|\rho(R\to R+r)|^2|,
\qquad
\int d^dr\,f(r)=1.
$$
Since $f(r)$ can only subtract one particle from the uniform background, $g(r)$ must dip below $\rho_0$ in a set whose total hole-volume is exactly one particle. By suitable choice of orbitals, that single “Pauli hole” can be split into multiple smaller holes at a periodic array of separation vectors; that pattern of anti-correlation dips is termed a “Pauli anti-crystal” [1903.02067].

In the concrete one-dimensional example on a ring of length $L$, choosing plane waves
$$
\psi_k(x)=\frac1{\sqrt L}e^{i2\pi kx/L},\qquad k=\Delta_k n,\ n=1,\dots,N,
$$
yields constructive interference in $|\rho(x\to y)|^2$ whenever $x-y=m(L/\Delta_k)$, producing equally spaced dips in $g(r)$ [1903.02067]. Thus, exclusion generates periodic deficits rather than periodic density peaks.

## 6. Neighbour counting statistics and the scope of “Pauli periodicity”

To characterize genuine local order without seed bias, [1903.02067] introduces neighbour counting statistics (NCS), defined as the full counting statistics of the number of neighbors within radius $r$ of a typical particle:
$$
p_n(r)=\frac1N\sum_{i=1}^N
\Big\langle
\delta_{\,n,\sum_{j\ne i}\Theta(r-|r_j-r_i|)}
\Big\rangle.
$$
The generating function is
$$
\chi(\lambda;r)=\sum_{n=0}^N p_n(r)e^{i\lambda n},
$$
the mean
$$
n(r)=\sum_n n\,p_n(r)
$$
reproduces the usual $g(r)$, and the variance
$$
\sigma^2(r)=\sum_n n^2p_n(r)-n(r)^2
$$
measures the rigidity of the $r$-shell [1903.02067]. The method is presented as more reliable because it requires no seed configuration or alignment, fully captures multipoint correlations around a particle, and directly reveals shell occupancies and stiffness.

The numerical applications illustrate the distinction between incipient order and artifact. For Laughlin states on the sphere with $N=200$, the $q=1$ state shows a featureless $\sigma^2(r)$ beyond the short exclusion hole, whereas the $q=90$ Wigner crystal exhibits plateaux in $n(r)$ at $6,18,36,\dots$ and deep minima in $\sigma^2(r)$; for a classical Lennard-Jones plus Gaussian-bump model, lowering temperature sharpens the shell counts $4,8,12,\dots$ [1903.02067]. These are bona fide shell signatures, but they are not identified with non-interacting “Pauli crystals.”

Across these literatures, “Pauli periodicity” therefore spans several non-equivalent notions. In spectral theory, it refers to exact equivalences between double zeros of discriminants and either $\pi/2$-periodicity or $\sigma_2$-symmetry after a Pauli-based gauge reduction [1707.00982]. In atomic density-functional modeling, it denotes shell periodicity emerging from a Pauli potential derived from excluded volume in thermal-space [2211.15764]. In many-body imaging, the term must be handled cautiously: non-interacting Pauli crystals are ruled out, while periodicity can survive only as anti-correlation holes or as genuine shell structure diagnosed by unbiased counting statistics [1903.02067]. A plausible implication is that the common element is not a single periodic object, but the way Pauli structure constrains admissible spectra, densities, and correlations.

Source: https://www.emergentmind.com/topics/pauli-periodicity