Hyperuniformity Disorder Length

• Hyperuniformity disorder length is a model-independent metric that quantifies the spatial scale over which density fluctuations are suppressed relative to Poisson randomness.
• Operational methods like scaled-variance thresholding, crossover analysis, and boundary-layer approaches enable extraction of this characteristic length in both simulations and experiments.
• Its variation across different material systems, such as crystals and jammed media, provides practical insights into designing materials with targeted fluctuation suppression.

A hyperuniformity disorder length quantifies the characteristic scale over which a many-particle, pixel, or cellular system exhibits suppressed density or field fluctuations compared to a reference random (e.g., Poisson) configuration. This length scale provides a physically meaningful, model-independent metric for the spatial extent of hyperuniform (or nearly hyperuniform) correlations in disordered and ordered materials. Various mathematical definitions and operational procedures connect real-space fluctuation statistics, spectral structure factors, and boundary-driven fluctuation analyses, yielding convergent interpretations of the hyperuniformity disorder length in computational, theoretical, and experimental contexts.

1. Mathematical Foundations and Definitions

Hyperuniformity is fundamentally characterized by the suppression of local density fluctuations at long wavelengths or large length scales. For a point configuration in $d$ dimensions, the number variance $\sigma_N^2(R)$ in a spherical window of radius $R$ is

$$\sigma_N^2(R) = \langle N(R)^2 \rangle - \langle N(R) \rangle^2.$$

A system is hyperuniform if

$\lim_{R\to\infty} \frac{\sigma_N^2(R)}{R^d}=0$

or, equivalently, in reciprocal space, if the structure factor $S(\mathbf k)$ satisfies

$\lim_{|\mathbf k|\to 0} S(\mathbf k) = 0$

(Torquato, 2018, Vanoni et al., 28 Jul 2025, Maher et al., 2024). In strictly hyperuniform systems, density fluctuations are suppressed so that the leading scaling of $\sigma_N^2(R)$ is governed by the surface-area term, i.e., $\sigma_N^2(R) \sim B_N R^{d-1}$.

The hyperuniformity disorder length—hereafter "$L_D$" (Editor's term)—furnishes a length scale quantifying how far hyperuniform-like order extends before fluctuations revert to non-hyperuniform (typically Poisson-like) behavior.

Several frameworks provide operational definitions:

• Variance-based approaches: $\sigma_N^2(R)$0 is the minimal radius $\sigma_N^2(R)$1 beyond which the scaled variance, suitably normalized, remains within a specified tolerance of its asymptotic hyperuniform scaling (Vanoni et al., 28 Jul 2025, Maher et al., 2024).
• Spectral approaches: $\sigma_N^2(R)$2 is connected to the small-$\sigma_N^2(R)$3 scaling of $\sigma_N^2(R)$4 as $\sigma_N^2(R)$5 and can be related to the leading nonzero coefficient in its expansion (Torquato, 2018).
• Boundary-layer perspective: $\sigma_N^2(R)$6 appears as the effective width of a region near the observation window's boundary where "residual" fluctuations are generated (Chieco et al., 2017, Chieco et al., 2021).

2. Operational Determination and Formulas

The most widely used procedures to extract $\sigma_N^2(R)$7 or related scales from data or simulation are:

• Scaled variance thresholding: In practical, finite systems, one forms the scaled variance

$\sigma_N^2(R)$8

with $\sigma_N^2(R)$9 an asymptotic prefactor and $R$0 the expected scaling exponent for the hyperuniform class (Vanoni et al., 28 Jul 2025). One then defines $R$1 as the minimum value of $R$2 for which $R$3 for a chosen tolerance $R$4.

• Crossover analysis: For instances where the number variance transitions between two scaling forms ($R$5 at small $R$6, $R$7 at large $R$8 in 2D), the crossover radius where the two terms in a fitted ansatz become comparable is taken as the disorder length (e.g., $R$9 for density, $$\sigma_N^2(R) = \langle N(R)^2 \rangle - \langle N(R) \rangle^2.$$0 for geometric observables) (Zheng et al., 2018).
• Real-space boundary approach: The hyperuniformity disorder length $$\sigma_N^2(R) = \langle N(R)^2 \rangle - \langle N(R) \rangle^2.$$1 is defined via the scaling of local volume-fraction variance $$\sigma_N^2(R) = \langle N(R)^2 \rangle - \langle N(R) \rangle^2.$$2 considering the fraction of fluctuations arising from a shell of thickness $$\sigma_N^2(R) = \langle N(R)^2 \rangle - \langle N(R) \rangle^2.$$3 around a window (Chieco et al., 2017):

$$\sigma_N^2(R) = \langle N(R)^2 \rangle - \langle N(R) \rangle^2.$$4

leading to an explicit inversion for $$\sigma_N^2(R) = \langle N(R)^2 \rangle - \langle N(R) \rangle^2.$$5

$$\sigma_N^2(R) = \langle N(R)^2 \rangle - \langle N(R) \rangle^2.$$6

where $$\sigma_N^2(R) = \langle N(R)^2 \rangle - \langle N(R) \rangle^2.$$7 is the variance ratio to random expectations.

• Spectroscopy approach: By direct computation or inversion of the measured variance or the spectral density, one can generate $$\sigma_N^2(R) = \langle N(R)^2 \rangle - \langle N(R) \rangle^2.$$8 or $$\sigma_N^2(R) = \langle N(R)^2 \rangle - \langle N(R) \rangle^2.$$9 spectra for further analysis (Chieco et al., 2021).

3. Disorder Lengths in Hyperuniformity Classes

The behavior and interpretation of $\lim_{R\to\infty} \frac{\sigma_N^2(R)}{R^d}=0$0 depend on the underlying hyperuniformity class, defined by small-$\lim_{R\to\infty} \frac{\sigma_N^2(R)}{R^d}=0$1 scaling of the structure factor $\lim_{R\to\infty} \frac{\sigma_N^2(R)}{R^d}=0$2. The large-$\lim_{R\to\infty} \frac{\sigma_N^2(R)}{R^d}=0$3 asymptotics for the number variance are:

• Class I ($\lim_{R\to\infty} \frac{\sigma_N^2(R)}{R^d}=0$4): $\lim_{R\to\infty} \frac{\sigma_N^2(R)}{R^d}=0$5
• Class II ($\lim_{R\to\infty} \frac{\sigma_N^2(R)}{R^d}=0$6): $\lim_{R\to\infty} \frac{\sigma_N^2(R)}{R^d}=0$7
• Class III ($\lim_{R\to\infty} \frac{\sigma_N^2(R)}{R^d}=0$8): $\lim_{R\to\infty} \frac{\sigma_N^2(R)}{R^d}=0$9

The disorder length $S(\mathbf k)$0 can be extracted analytically in terms of correction amplitudes and the asymptotic scaling, for example:

• Class I: $S(\mathbf k)$1 where $S(\mathbf k)$2 is the leading correction amplitude and $S(\mathbf k)$3 the power of the subleading scaling.
• Class III: $S(\mathbf k)$4 (Vanoni et al., 28 Jul 2025).

Empirically, Class I systems (crystals, quasicrystals, stealthy hyperuniform configurations) achieve hyperuniform suppression down to short scales (small $S(\mathbf k)$5). Class II and III systems manifest extended crossover regions, requiring larger observation windows before reaching their asymptotic fluctuations regime.

4. Application to Disordered and Jammed Systems

In disordered or jammed systems, the hyperuniformity disorder length provides a means to quantify the spatial reach of fluctuation suppression. For instance, in Voronoi jamming models, density fluctuations are suppressed only up to a finite length $S(\mathbf k)$6, which diverges at the rigidity transition ($S(\mathbf k)$7): $S(\mathbf k)$8 Beyond $S(\mathbf k)$9, density fluctuations revert to Poisson scaling. Perimeter fluctuations (a geometric observable) are suppressed over an even greater (but still finite) range $\lim_{|\mathbf k|\to 0} S(\mathbf k) = 0$0, diverging with a different exponent: $\lim_{|\mathbf k|\to 0} S(\mathbf k) = 0$1 Throughout the rigid phase, $\lim_{|\mathbf k|\to 0} S(\mathbf k) = 0$2, demonstrating that such systems suppress geometric fluctuations more efficiently than density ones (Zheng et al., 2018).

A similar paradigm holds in absorbing-state models, where the hyperuniform length $\lim_{|\mathbf k|\to 0} S(\mathbf k) = 0$3 separates scales with anomalous density fluctuations from those with ordinary central limit behavior. Higher moments of coarse-grained densities define an "extended" correlation length $\lim_{|\mathbf k|\to 0} S(\mathbf k) = 0$4, marking scales over which many-body (non-Gaussian) correlations persist (Zheng et al., 2020).

5. Hyperuniformity Disorder Length Spectroscopy and Computational Protocols

Hyperuniformity disorder length spectroscopy (HUDLS) emphasizes the extraction and interpretation of $\lim_{|\mathbf k|\to 0} S(\mathbf k) = 0$5 or $\lim_{|\mathbf k|\to 0} S(\mathbf k) = 0$6 as a function of window size, thereby connecting the real-space statistics to the underlying spatial organization. The practical workflow is:

1. Compute the local variance (number or volume fraction) as a function of window size.
2. Normalize to the random or ideal-hyperuniform expectation.
3. Invert the analytical variance formula for $\lim_{|\mathbf k|\to 0} S(\mathbf k) = 0$7 or fit the scaling behavior of the variance to extract $\lim_{|\mathbf k|\to 0} S(\mathbf k) = 0$8.
4. Identify the crossover scales and interpret $\lim_{|\mathbf k|\to 0} S(\mathbf k) = 0$9 asymptotics:
   • $\sigma_N^2(R)$0: Poisson/disordered.
   • $\sigma_N^2(R)$1: class-III or weakly hyperuniform.
   • $\sigma_N^2(R)$2 (constant): class-I or strong hyperuniformity.

This procedure provides a direct measure of the physical thickness of the "fluctuation-generating shell" and enables comparison across different material realizations, particle shapes, dimensionalities, and underlying order (Chieco et al., 2017, Chieco et al., 2021, Durian, 2017).

6. Generalizations, Geometrical and Field Fluctuations

The concept of a disorder length extends naturally to generalized observables and shapes:

• Geometric fields: The variance analysis can be applied to cell perimeters, areas, or edge networks in cellular or foam systems, yielding respective disorder lengths (e.g., $\sigma_N^2(R)$3 for area-weighted foams, $\sigma_N^2(R)$4 for perimeters) (Zheng et al., 2018, Chieco et al., 2021).
• Extended particles: The framework generalizes to polydisperse and nonpoint objects by evaluating specific overlap integrals and using tailored variance normalization (Durian, 2017).
• Polydispersity and discretization: The procedure accommodates mixtures and is robust across pixelated or continuum representations, provided variance formulas are properly matched to the system's microstructure.

7. Interpretation and Physical Significance

The hyperuniformity disorder length provides a unifying length scale for the spatial reach of suppressed fluctuations in random, hyperuniform, and nearly hyperuniform systems. It is model-independent and ties directly to measurable statistical minima, crossover phenomena, and fluctuation sources. A small $\sigma_N^2(R)$5 (or $\sigma_N^2(R)$6) implies nearly complete suppression of fluctuations over short scales, as in perfect crystals or optimal disordered hyperuniform states, whereas a large $\sigma_N^2(R)$7 signals a long crossover region with residual disorder. The distinct scaling of $\sigma_N^2(R)$8 across universality classes facilitates material comparison, the design of functional random media, and elucidates the limits of uniformity achievable in finite, noisy, or jammed systems (Vanoni et al., 28 Jul 2025, Chieco et al., 2021, Maher et al., 2024, Zheng et al., 2018, Chieco et al., 2017, Durian, 2017, Torquato, 2018).

---

Key References:

• (Torquato, 2018) Hyperuniform States of Matter
• (Zheng et al., 2018) Hyperuniformity and generalized fluctuations at Jamming
• (Chieco et al., 2017) Characterizing Pixel and Point Patterns with a Hyperuniformity Disorder Length
• (Durian, 2017) Hyperuniformity Disorder Length Spectroscopy for Extended Particles
• (Vanoni et al., 28 Jul 2025) When does hyperuniformity lead to uniformity across length scales?
• (Maher et al., 2024) Local order metrics for many-particle systems across length scales
• (Chieco et al., 2021) Quantifying the Long-Range Structure of Foams and Other Cellular Patterns with Hyperuniformity Disorder Length Spectroscopy
• (Zheng et al., 2020) Hidden Order Beyond Hyperuniformity in Critical Absorbing States