Pauli Periodicity: Spectral and Atomic Perspectives

• Pauli periodicity is a multifaceted concept characterizing periodic spectral properties in Dirac systems, shell organization in atoms, and correlation patterns in many-fermion configurations.
• The spectral theory approach employs a Pauli matrix decomposition to reveal half-periodicity and σ2-symmetry, leading to precise criteria for gap closures and instability interval vanishing.
• Methodologies like gauge transformations, Thomas-Fermi scaling, and unbiased neighbor counting offer actionable insights for bridging mathematical formalism with experimental quantum and atomic physics.

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 (Currie et al., 2017). 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 (Thompson, 2022). 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” (Fremling et al., 2019). 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$ (Currie et al., 2017). 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 $\sigma_2$0 and encode the natural $\sigma_2$1-decomposition of $\sigma_2$2.

The potential is written as

$\sigma_2$3

where $\sigma_2$4 belongs to the subspace skew-commuting with $\sigma_2$5 and $\sigma_2$6 belongs to the subspace commuting with $\sigma_2$7 (Currie et al., 2017). In components,

$\sigma_2$8

Because $\sigma_2$9, one obtains the explicit gauge transformation

$$J Y'(z,\lambda) + Q(z) Y(z,\lambda) = \lambda Y(z,\lambda), \qquad 0\le z\le \pi,$$0

which conjugates the original system into one with new potential

$$J Y'(z,\lambda) + Q(z) Y(z,\lambda) = \lambda Y(z,\lambda), \qquad 0\le z\le \pi,$$1

while preserving the same boundary conditions at $$J Y'(z,\lambda) + Q(z) Y(z,\lambda) = \lambda Y(z,\lambda), \qquad 0\le z\le \pi,$$2 (Currie et al., 2017).

This Pauli-based decomposition is the structural origin of the periodicity theorems. The periodicity properties are not imposed directly on the original $$J Y'(z,\lambda) + Q(z) Y(z,\lambda) = \lambda Y(z,\lambda), \qquad 0\le z\le \pi,$$3; 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 $$J Y'(z,\lambda) + Q(z) Y(z,\lambda) = \lambda Y(z,\lambda), \qquad 0\le z\le \pi,$$4 and $$J Y'(z,\lambda) + Q(z) Y(z,\lambda) = \lambda Y(z,\lambda), \qquad 0\le z\le \pi,$$5 organize the gauge-invariant content of the potential.

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

The monodromy matrix is $$J Y'(z,\lambda) + Q(z) Y(z,\lambda) = \lambda Y(z,\lambda), \qquad 0\le z\le \pi,$$6, and the discriminant is defined by

$$J Y'(z,\lambda) + Q(z) Y(z,\lambda) = \lambda Y(z,\lambda), \qquad 0\le z\le \pi,$$7

Since $$J Y'(z,\lambda) + Q(z) Y(z,\lambda) = \lambda Y(z,\lambda), \qquad 0\le z\le \pi,$$8, the Floquet multipliers $$J Y'(z,\lambda) + Q(z) Y(z,\lambda) = \lambda Y(z,\lambda), \qquad 0\le z\le \pi,$$9 satisfy

$J=\begin{pmatrix}0&1\-1&0\end{pmatrix},$0

(Currie et al., 2017). The commuting part $J=\begin{pmatrix}0&1\-1&0\end{pmatrix},$1 contributes the overall factor

$J=\begin{pmatrix}0&1\-1&0\end{pmatrix},$2

so that the discriminant of the conjugated system is

$J=\begin{pmatrix}0&1\-1&0\end{pmatrix},$3

and the relevant entire functions are $J=\begin{pmatrix}0&1\-1&0\end{pmatrix},$4.

Theorem A identifies half-periodicity through the zero structure of $J=\begin{pmatrix}0&1\-1&0\end{pmatrix},$5: the zeros of

$J=\begin{pmatrix}0&1\-1&0\end{pmatrix},$6

are all double if and only if the original system is unitarily equivalent, via $J=\begin{pmatrix}0&1\-1&0\end{pmatrix},$7, to one whose Hermitian potential $J=\begin{pmatrix}0&1\-1&0\end{pmatrix},$8 is $J=\begin{pmatrix}0&1\-1&0\end{pmatrix},$9-periodic almost everywhere (Currie et al., 2017). In the forward direction, $Q(z)$0-periodicity implies that $Q(z)$1 and $Q(z)$2 satisfy the same Dirac equation and differ by the constant matrix $Q(z)$3; setting $Q(z)$4 gives

$Q(z)$5

In the reverse direction, assuming that $Q(z)$6 has only double zeros, a Hadamard factorization argument together with high-energy asymptotics shows that

$Q(z)$7

vanishes identically, and the same matrix identity $Q(z)$8 then forces $Q(z)$9 (Currie et al., 2017).

Theorem B gives the complementary $\pi$0-symmetry criterion. The zeros of

$\pi$1

are all double if and only if the system is unitarily equivalent, via the same $\pi$2, to one satisfying

$\pi$3

(Currie et al., 2017). If $\pi$4 is $\pi$5-symmetric, then $\pi$6 and $\pi$7 solve the same Dirac equation and differ by the constant $\pi$8, leading to

$\pi$9

Conversely, the double-zero structure of $Y(z,\lambda)$0 forces the same matrix equality and thereby the $Y(z,\lambda)$1-symmetry of $Y(z,\lambda)$2 (Currie et al., 2017).

These theorems constitute the core Pauli periodicity mechanism in the spectral setting. Half-periodicity is encoded through $Y(z,\lambda)$3, whereas $Y(z,\lambda)$4-conjugation symmetry is encoded through $Y(z,\lambda)$5.

3. Vanishing instability intervals and the canonical $Y(z,\lambda)$6 form

The same formalism yields an Ambarzumyan-type statement for Dirac systems. All spectral instability intervals, defined as the open gaps in $Y(z,\lambda)$7, vanish if and only if, after the $Y(z,\lambda)$8-gauge, the potential has the form

$Y(z,\lambda)$9

with $2\times2$0 and $2\times2$1 real and $2\times2$2-periodic (Currie et al., 2017). Equivalently, in the formulation given in the abstract, all instability intervals vanish if and only if

$2\times2$3

for some real-valued $2\times2$4-periodic functions $2\times2$5 and $2\times2$6 integrable on compact sets.

The proof combines Theorems A and B. Vanishing of all gaps means that for every real $2\times2$7 both $2\times2$8 and $2\times2$9 have only double zeros. This forces $Y(0,\lambda)=I$0 to be simultaneously $Y(0,\lambda)=I$1-periodic and $Y(0,\lambda)=I$2-symmetric. Periodicity together with $Y(0,\lambda)=I$3-anti-periodicity on the skew part $Y(0,\lambda)=I$4 implies $Y(0,\lambda)=I$5, so the potential reduces to the commuting $Y(0,\lambda)=I$6 sector (Currie et al., 2017).

In that case the system can be solved explicitly:

$Y(0,\lambda)=I$7

and

$Y(0,\lambda)=I$8

Hence $Y(0,\lambda)=I$9 for all real $\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}$0, so no gaps occur (Currie et al., 2017).

In this sense, Pauli periodicity is not merely a descriptive label but a complete spectral classification. The joint action of $\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}$1 and $\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}$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 $\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}$3, and a Pauli potential is postulated as a Flory-Huggins-type mean-field excluded-volume penalty between the four-dimensional contours of different electrons (Thompson, 2022). With pair densities $\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}$4, the total Pauli energy is

$\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}$5

with

$\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}$6

This yields the local form

$\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}$7

where $\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}$8.

The uniform-limit scaling argument connects this Pauli potential to Thomas-Fermi theory. In mean field, $\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}$9, which has the wrong density scaling compared to the Thomas-Fermi kinetic energy density

$\sigma_2$00

A polymer-scaling argument instead gives

$\sigma_2$01

in precise agreement with Thomas-Fermi scaling, and matching prefactors fixes

$\sigma_2$02

with no further adjustable parameters (Thompson, 2022).

The self-consistent-field equations follow from the free-energy functional

$\sigma_2$03

leading to the mean-field potentials

$\sigma_2$04

where the electron-electron Hartree term carries a Fermi-Amaldi factor $\sigma_2$05 to reduce self-interaction (Thompson, 2022). The auxiliary propagators satisfy

$\sigma_2$06

with

$\sigma_2$07

These equations are iterated to self-consistency.

Under spherical averaging, each $\sigma_2$08 and $\sigma_2$09 is expanded in spherical Bessel functions $\sigma_2$10, and the total radial density is reconstructed as

$\sigma_2$11

where $\sigma_2$12 labels grouped shell or sub-shell pairs such as $\sigma_2$13, $\sigma_2$14, and $\sigma_2$15 (Thompson, 2022). Numerical calculations at $\sigma_2$16 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 $\sigma_2$17 for the lightest elements and drop to $\sigma_2$18 or less for atoms heavier than nitrogen; beyond $\sigma_2$19 the binding-energy error stabilizes around $\sigma_2$20–$\sigma_2$21 (Thompson, 2022).

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” (Thompson, 2022).

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 $\sigma_2$22-particle snapshots to a chosen seed configuration $\sigma_2$23 by minimizing the configuration distance

$\sigma_2$24

and then averages the aligned one-body densities:

$\sigma_2$25

Peaks in $\sigma_2$26 had been taken as evidence of “Pauli crystals,” but the analysis in (Fremling et al., 2019) 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 $\sigma_2$27 and $\sigma_2$28 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 (Fremling et al., 2019). The configuration variance

$\sigma_2$29

further shows that the CD peaks do not imply rigid particle binding to seed sites: for random points $\sigma_2$30, for non-interacting fermions with $\sigma_2$31 one finds $\sigma_2$32–$\sigma_2$33, and for a strongly correlated liquid with $\sigma_2$34 it drops further but does not vanish (Fremling et al., 2019).

The alternative structure proposed in the paper is the “Pauli anti-crystal.” For a single Slater determinant with occupied orbitals $\sigma_2$35, the one- and two-body densities are

$\sigma_2$36

and

$\sigma_2$37

For flat one-body density $\sigma_2$38, the reduced pair correlator is

$\sigma_2$39

with

$\sigma_2$40

Since $\sigma_2$41 can only subtract one particle from the uniform background, $\sigma_2$42 must dip below $\sigma_2$43 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” (Fremling et al., 2019).

In the concrete one-dimensional example on a ring of length $\sigma_2$44, choosing plane waves

$\sigma_2$45

yields constructive interference in $\sigma_2$46 whenever $\sigma_2$47, producing equally spaced dips in $\sigma_2$48 (Fremling et al., 2019). 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, (Fremling et al., 2019) introduces neighbour counting statistics (NCS), defined as the full counting statistics of the number of neighbors within radius $\sigma_2$49 of a typical particle:

$\sigma_2$50

The generating function is

$\sigma_2$51

the mean

$\sigma_2$52

reproduces the usual $\sigma_2$53, and the variance

$\sigma_2$54

measures the rigidity of the $\sigma_2$55-shell (Fremling et al., 2019). 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 $\sigma_2$56, the $\sigma_2$57 state shows a featureless $\sigma_2$58 beyond the short exclusion hole, whereas the $\sigma_2$59 Wigner crystal exhibits plateaux in $\sigma_2$60 at $\sigma_2$61 and deep minima in $\sigma_2$62; for a classical Lennard-Jones plus Gaussian-bump model, lowering temperature sharpens the shell counts $\sigma_2$63 (Fremling et al., 2019). 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 $\sigma_2$64-periodicity or $\sigma_2$65-symmetry after a Pauli-based gauge reduction (Currie et al., 2017). In atomic density-functional modeling, it denotes shell periodicity emerging from a Pauli potential derived from excluded volume in thermal-space (Thompson, 2022). 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 (Fremling et al., 2019). 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.