Stealthy-Hyperuniformity Overview

• Stealthy-hyperuniformity is defined by the complete suppression of scattering over a finite wavevector range, setting it apart from generic hyperuniform systems.
• The collective-coordinate formulation enforces S(k)=0 on constrained modes, enabling precise control over disorder, isotropy, and tunable short-range order.
• This approach yields practical benefits such as enhanced photonic band gaps, accelerated transport, and improved connectivity in both ordered and disordered media.

Stealthy-hyperuniformity denotes a stronger-than-hyperuniform suppression of density fluctuations in statistically homogeneous point configurations and two-phase media: the structure factor satisfies $S(\mathbf{k})=0$ throughout a finite exclusion region $0<|\mathbf{k}|\le K$, rather than merely vanishing in the limit $|\mathbf{k}| \to 0$. Because this imposes an exact spectral gap around the origin, stealthy systems are Class I hyperuniform, exhibit anomalously small number-variance growth, and completely suppress single scattering over a finite band of wavevectors (Zhang et al., 2015). In current usage, the concept spans point patterns, particle packings, spin chains, networks, and composite media whose microstructure is constrained in Fourier space so that long-range homogeneity coexists with disorder, isotropy, and tunable short-range order.

1. Formal definition and statistical descriptors

For a statistically homogeneous point configuration of number density $\rho$, the pair correlation function $g_2(r)$ and total correlation function $h(r)=g_2(r)-1$ determine the structure factor through

$S(\mathbf{k})=1+\rho \,\tilde h(\mathbf{k}).$

Hyperuniformity is the condition

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

which is equivalent to anomalously suppressed long-wavelength density fluctuations. In real space, the number variance $\sigma_N^2(R)$ in a spherical observation window of radius $R$ grows more slowly than the window volume; for Class I systems, including stealthy ones, $0<|\mathbf{k}|\le K$0 as $0<|\mathbf{k}|\le K$1 (Zhang et al., 2015).

Stealthy-hyperuniformity strengthens this requirement by enforcing

$0<|\mathbf{k}|\le K$2

or, more generally for two-phase media,

$0<|\mathbf{k}|\le K$3

where $0<|\mathbf{k}|\le K$4 is the spectral density, i.e. the Fourier transform of the autocovariance. In packings of identical spheres or disks, these descriptions are directly linked through

$0<|\mathbf{k}|\le K$5

so stealthiness at the point-pattern level induces stealthiness in the corresponding two-phase medium (Kim et al., 2024).

This finite exclusion region in reciprocal space is the defining distinction from generic hyperuniformity. It produces an exact low-$0<|\mathbf{k}|\le K$6 gap in scattering, not merely a small-$0<|\mathbf{k}|\le K$7 asymptotic decay, and this sharper constraint is responsible for many of the optical, transport, and rigidity properties associated with the subject.

2. Collective-coordinate formulation and model interactions

A standard realization of stealthy-hyperuniformity uses collective coordinates in a periodic fundamental cell of volume $0<|\mathbf{k}|\le K$8. For $0<|\mathbf{k}|\le K$9 points at positions $|\mathbf{k}| \to 0$0, the collective density variable is

$|\mathbf{k}| \to 0$1

with single-configuration structure factor $|\mathbf{k}| \to 0$2 for $|\mathbf{k}| \to 0$3. If a pair interaction has Fourier transform $|\mathbf{k}| \to 0$4 that is nonnegative and supported only on $|\mathbf{k}| \to 0$5, then the total energy can be written as

$|\mathbf{k}| \to 0$6

Equivalently, over symmetry-independent constrained modes,

$|\mathbf{k}| \to 0$7

Since $|\mathbf{k}| \to 0$8 inside the stealthy window, minimizing $|\mathbf{k}| \to 0$9 enforces $\rho$0, hence $\rho$1, for all constrained wavevectors (Zhang et al., 2015).

The canonical model family takes

$\rho$2

With $\rho$3 fixed, the inverse Fourier transform yields bounded, long-ranged, oscillatory real-space pair potentials. The force on particle $\rho$4 follows directly: $\rho$5

This formulation also clarifies why numerical generation is subtle. Pure energy minimization from random initial conditions samples the degenerate ground-state manifold in an algorithm-dependent way. Canonical $\rho$6 sampling instead uses low-temperature molecular dynamics followed by minimization, which removes that artifact and, when done at sufficiently low equilibration temperature, yields ensembles that are insensitive to the detailed choice of $\rho$7 within the stealthy window (Zhang et al., 2015). High-accuracy collective-coordinate optimization later pushed the residual “distance to stealthiness” to $\rho$8 for systems up to $\rho$9 in two dimensions, showing that ultra-accurate enforcement of the constraints can be essential for diagnosing the energy landscape and wave-transport behavior (Morse et al., 2023).

3. Constraint fraction, phase behavior, and entropic selection

The natural control parameter is the fraction of constrained degrees of freedom,

$g_2(r)$0

where $g_2(r)$1 is the number of independent constrained wavevectors with $g_2(r)$2. For fixed $g_2(r)$3, decreasing the number density increases $g_2(r)$4, while increasing density decreases $g_2(r)$5. In the thermodynamic limit, $g_2(r)$6 is proportional to the volume of the exclusion ball in $g_2(r)$7-space and inversely proportional to density (Zhang et al., 2015).

For $g_2(r)$8, the canonical $g_2(r)$9 ensemble selects disordered stealthy ground states. These states have no Bragg peaks, satisfy $h(r)=g_2(r)-1$0 for $h(r)=g_2(r)-1$1, and develop increasing short-range order as $h(r)=g_2(r)-1$2 increases. Their structure factor just above $h(r)=g_2(r)-1$3 obeys the “pseudo” hard-sphere ansatz in Fourier space at small $h(r)=g_2(r)-1$4: the first peak of $h(r)=g_2(r)-1$5 above $h(r)=g_2(r)-1$6 plays the role of the direct-space first-neighbor peak of an equilibrium hard-sphere fluid. This correspondence is accurate up to $h(r)=g_2(r)-1$7 values on the order of $h(r)=g_2(r)-1$8–$h(r)=g_2(r)-1$9, depending on dimension, and then progressively breaks down as higher-order structure grows (Zhang et al., 2015).

At $S(\mathbf{k})=1+\rho \,\tilde h(\mathbf{k}).$0, the entropically favored states are crystalline in the infinite-system-size limit. In two dimensions, triangular lattices are favored over square lattices in Wang–Landau sampling; in three dimensions, evidence indicates that face-centered cubic order is favored near $S(\mathbf{k})=1+\rho \,\tilde h(\mathbf{k}).$1, with face-centered cubic and body-centered cubic entropy scalings becoming nearly indistinguishable at higher $S(\mathbf{k})=1+\rho \,\tilde h(\mathbf{k}).$2 (Zhang et al., 2015). Stacked-slider phases also belong to the exact ground-state manifold above $S(\mathbf{k})=1+\rho \,\tilde h(\mathbf{k}).$3, but they are not entropically favored in the canonical zero-temperature ensemble. They are nonperiodic, statistically anisotropic, possess long-range orientational order, and have zero shear modulus because exact sliding motions preserve all stealthy constraints (Chertkov et al., 2015).

A complementary geometric result sharpens this phase picture. Under the Fourier-space hard-sphere analogy, the densest Fourier-space hard-sphere configuration is a Bravais lattice, specifically the reciprocal lattice of the densest real-space Bravais lattice. Thus, while disordered stealthy ground states exist for $S(\mathbf{k})=1+\rho \,\tilde h(\mathbf{k}).$4 in $S(\mathbf{k})=1+\rho \,\tilde h(\mathbf{k}).$5, the densest possible stealthy configurations are ordered in every dimension (Morse et al., 2024).

4. Structural realizations: glasses, packings, polycrystals, and stackings

As $S(\mathbf{k})=1+\rho \,\tilde h(\mathbf{k}).$6, the degeneracy of the stealthy ground-state manifold becomes maximal. A modified collective-coordinate optimization that adds a short-ranged soft-core repulsion to the usual stealthy energy exploits this degeneracy to construct ultradense packings. In three dimensions, the resulting $S(\mathbf{k})=1+\rho \,\tilde h(\mathbf{k}).$7 packings are configurationally extremely close to maximally random jammed sphere packings: they have packing fraction $S(\mathbf{k})=1+\rho \,\tilde h(\mathbf{k}).$8, mean contact number $S(\mathbf{k})=1+\rho \,\tilde h(\mathbf{k}).$9, gap exponent $\lim_{|\mathbf{k}|\to 0} S(\mathbf{k})=0,$0, and small-$\lim_{|\mathbf{k}|\to 0} S(\mathbf{k})=0,$1 scaling $\lim_{|\mathbf{k}|\to 0} S(\mathbf{k})=0,$2, i.e. $\lim_{|\mathbf{k}|\to 0} S(\mathbf{k})=0,$3, matching the maximally random jammed state (Torquato et al., 14 Apr 2025).

The same soft-core construction yields ultradense disordered stealthy-hyperuniform packings throughout the disordered regime $\lim_{|\mathbf{k}|\to 0} S(\mathbf{k})=0,$4. In two dimensions the attainable packing fractions span $\lim_{|\mathbf{k}|\to 0} S(\mathbf{k})=0,$5, and in three dimensions $\lim_{|\mathbf{k}|\to 0} S(\mathbf{k})=0,$6, values far exceeding those obtained without the extra repulsion. Without soft-core repulsion, the maximal achievable packing fraction decreases to zero on average as $\lim_{|\mathbf{k}|\to 0} S(\mathbf{k})=0,$7 increases for $\lim_{|\mathbf{k}|\to 0} S(\mathbf{k})=0,$8; with it, the maximal packing fraction remains large and essentially $\lim_{|\mathbf{k}|\to 0} S(\mathbf{k})=0,$9-independent (Vanoni et al., 31 Mar 2025). In this regime, increasing $\sigma_N^2(R)$0 thins the contact network, eventually producing linear polymer-like chains of contacting particles with progressively shorter mean chain lengths.

Polycrystalline stealthy-hyperuniform media constitute a distinct realization. In two dimensions, collective-coordinate optimization generates ultradense polycrystalline disk packings with $\sigma_N^2(R)$1 between $\sigma_N^2(R)$2 and $\sigma_N^2(R)$3, local triangular order, and grain boundaries, while still maintaining $\sigma_N^2(R)$4 for $\sigma_N^2(R)$5. These polycrystals differ sharply from Lubachevsky–Stillinger reference polycrystals: their grain-size distribution develops a pronounced peak, and the grain–boundary two-phase autocovariance exhibits long-range oscillations, indicating inter-grain correlations induced by the stealthy constraint (Vanoni et al., 22 Jun 2026).

Stealthiness also occurs in ordered stacking families. All Barlow packings, including stacking-disordered ones, are stealthy hyperuniform by the stealthy-stacking theorem, with the common lower bound

$\sigma_N^2(R)$6

Within that family, the hyperuniformity order metric $\sigma_N^2(R)$7 is, to a very good approximation, linear in the fraction of fcc-like clusters, showing that large-scale fluctuation suppression is governed primarily by local three-layer stacking geometry (Middlemas et al., 2018). This illustrates a broader point: stealthiness does not imply either full periodicity or full amorphousness, but rather a reciprocal-space constraint that can coexist with crystals, polycrystals, and highly degenerate disordered states.

5. Wave, transport, and network consequences

The most direct physical consequence of stealthy-hyperuniformity is the exact suppression of scattering for $\sigma_N^2(R)$8. In two-phase media, nonlocal strong-contrast theory shows that layered and transversely isotropic stealthy-hyperuniform media possess finite transparency windows in wavevector space, and that for layered transverse polarization and transversely isotropic TM polarization this perfect transparency persists through third order in the expansion. Within the predicted interval there is, in practice, no Anderson localization because the localization length is much larger than any practically large sample size (Kim et al., 2024). The same low-$\sigma_N^2(R)$9 gap also accelerates interphase diffusion: the excess spreadability decays exponentially in stealthy media, whereas nonstealthy media exhibit algebraic long-time tails.

These principles have been realized experimentally in photonics. In three dimensions, an optimized disordered stealthy-hyperuniform dielectric network fabricated for microwaves exhibits the first experimental isotropic photonic band gap of this type. For refractive index $R$0, the complete isotropic gap is centered at $R$1 with width $R$2; at $R$3, not fabricated in that work, the calculated gap widens to $R$4 and shifts to $R$5 (Siedentop et al., 2024). A periodic diamond network of similar node density shows anisotropic gaps, while an amorphous non-stealthy reference shows no complete gap, demonstrating that isotropy alone is insufficient.

A related semiconductor realization uses a hole-based stealthy-hyperuniform pattern on a quantum-cascade platform for the mid-infrared. The measured TE reflection band is centered near $R$6 with a gap-midgap ratio of $R$7, and the reflection spectrum remains unchanged under in-plane rotation, directly confirming spatial isotropy in the optical response (Gallego et al., 2024).

Transport and effective-medium advantages extend beyond photonic gaps. Two-phase media derived from ultradense stealthy-hyperuniform packings exhibit faster diffusion spreadability than media generated without soft-core repulsion, and the imaginary part of the effective dynamic dielectric constant vanishes at a small wavevector, implying perfect transparency for the corresponding modes. Cross-property relations link the transparency-window size to the long-time spreadability rate, since both are governed by the same spectral density $R$8 (Vanoni et al., 31 Mar 2025).

Network models inherit related benefits. Delaunay triangulation networks built from disordered stealthy-hyperuniform point sets have lower percolation thresholds than Poisson networks under distance-dependent bond occupation, and the threshold decreases as $R$9 increases. In two dimensions, systems with large $0<|\mathbf{k}|\le K$00 fall into the same universality class as lattice percolation, whereas Poisson and low-$0<|\mathbf{k}|\le K$01 systems exhibit shifted critical exponents (Wang et al., 17 Mar 2026). This suggests that suppressing long-wavelength density fluctuations systematically improves global connectivity in statistically homogeneous disordered networks.

6. Generalizations, rigorous results, and current points of interpretation

Stealthy-hyperuniformity is not restricted to isotropic Hermitian point sets. In two-phase composites, inverse design can impose arbitrary realizable exclusion regions $0<|\mathbf{k}|\le K$02 in Fourier space, producing anisotropic stealthy-hyperuniform microstructures with circular, elliptical, square, rectangular, butterfly-shaped, or lemniscate-shaped spectral gaps. These media may be directionally hyperuniform—stealthy along selected directions but not others—while still possessing nearly isotropic local morphology (Shi et al., 2023). In one-dimensional deep-subwavelength multilayers, a stealthy gap in the spectrum of $0<|\mathbf{k}|\le K$03 classifies intermediate phases between crystals and uncorrelated disorder and enables angle-selective control of localization through the breakdown of effective-medium theory by Goos–Hänchen-dominated interface physics (Park et al., 2024).

Open-system wave physics has introduced a further extension. For complex material potentials $0<|\mathbf{k}|\le K$04, the scattering intensity splits into an even part

$0<|\mathbf{k}|\le K$05

and an odd real–imaginary cross term

$0<|\mathbf{k}|\le K$06

In this setting, hyperuniformity requires only that the real and imaginary parts are each hyperuniform, whereas stealthiness depends essentially on the cross-correlations because cancellation between $0<|\mathbf{k}|\le K$07 and $0<|\mathbf{k}|\le K$08 can create unidirectional stealthy phases unavailable in Hermitian media and in non-Hermitian crystals (Lee et al., 10 Feb 2026). Experimentally, residual single scattering inside the nominal stealthy region has also been observed in two-dimensional photonic-crystal slabs with stealthy-hyperuniform disorder, where a complex effective mass induced by radiative loss makes the dynamics intrinsically non-Hermitian and leaves a finite linewidth even below the nominal $0<|\mathbf{k}|\le K$09 threshold for Hermitian single scattering (Barsukova et al., 7 Jul 2025).

The concept also has rigorous probabilistic consequences. For translation-invariant stochastic processes whose structure function vanishes on an open set $0<|\mathbf{k}|\le K$10, the exterior configuration outside any bounded domain determines the exact interior configuration with probability $0<|\mathbf{k}|\le K$11; such systems are maximally rigid. When $0<|\mathbf{k}|\le K$12, stealthy-hyperuniform point processes additionally have bounded holes, with the maximal hole size bounded above by a constant proportional to the inverse radius of the spectral gap (Ghosh et al., 2017). These results supply mathematically precise versions of two empirically central features: suppressed voids and unusually strong long-range constraints.

Several apparent counterexamples have therefore been re-evaluated. Gyromorphs, whose structure factors contain rings of Bragg-like peaks with discrete rotational symmetry, were proposed as disordered media reproducing stealthy-hyperuniform advantages without stealthiness or hyperuniformity. Subsequent analysis showed instead that, in the large-$0<|\mathbf{k}|\le K$13 nearly isotropic limit, gyromorphs are hyperuniform of Class III, not stealthy, and display size-dependent pseudogaps densely populated by localized states rather than the smooth gaps of highly stealthy hyperuniform materials (Skolnick et al., 21 Jun 2026). The broader implication is that isotropy, disorder, or even hyperuniformity alone do not reproduce the full phenomenology of stealthy-hyperuniformity; the finite exclusion region $0<|\mathbf{k}|\le K$14 remains the distinctive structural ingredient.