---
title: Correlated-Disordered Metasurfaces
url: https://www.emergentmind.com/topics/correlated-disordered-metasurfaces
type: topic
---

# Correlated-Disordered Metasurfaces

Correlated-disordered metasurfaces are planar assemblies of meta-atoms or scatterers that lack global periodicity but retain engineered spatial correlations, so that pair separations, local exclusion rules, reciprocal-space constraints, or connectivity statistics become design variables rather than fabrication byproducts. In this setting, disorder is not restricted to Poisson-like randomness: it can be imposed through minimum inter-particle distances, Gaussian-correlated lattice perturbations, hyperuniform or stealthy-hyperuniform point patterns, shuffled lattices, or connected morphologies near a critical packing threshold. Across these realizations, correlated disorder is used to control specular and diffuse scattering, collective resonances, photon density of states, localization, directional emission, broadband absorption, and visual appearance in ways that are not accessible with either perfect periodicity or purely uncorrelated disorder [2306.13540], [2310.14349].

## 1. Definitions, disorder classes, and structural control

A central distinction is between uncorrelated and correlated positional disorder. In perturbed lattices, uncorrelated disorder is generated by independently displacing each lattice site by a random vector, whereas correlated disorder introduces spatial smoothing so that nearby sites move coherently [2306.13540]. In one generative model for metasurfaces of pitch $P$, each nanostructure at position $\mathbf r_i$ is displaced by $\Delta \mathbf r_i$ with components uniformly drawn in $[-S_d P, S_d P]$, and correlations are added through a Gaussian weight
$$
C_{ij}=\exp\!\left\{-\left[\frac{r_{ij}}{2L_cP}\right]^2\right\},
$$
with full width at half maximum $2\sqrt{2\ln 2}\,L_cP$, so that the total perturbation is
$$
\mathbf r_i'=\mathbf r_i+\Delta \mathbf r_i+\sum_{j\neq i}\Delta \mathbf r_j\,C_{ij}.
$$
Here $S_d$ controls displacement amplitude and $L_c$ the correlation length and smoothness of the distortion [2306.13540].

A second broad class arises from exclusion constraints. Random sequential addition, Poisson-disk sampling, and related hard-core processes enforce a minimum spacing and therefore generate short-range repulsion and short-range order rather than independent point placement. In the optics literature this correlated-disorder class is frequently used to suppress diffuse scattering near the specular direction or to regularize fabrication-robust layouts [2212.02541], [2310.14349]. In shuffled lattices, such as randomly jittered silicon nanopillar arrays of nominal period $P_0=310\,\mathrm{nm}$, short-range positional correlations remain centered around the original lattice spacing while long-range order and Bragg peaks are progressively smeared as the disorder parameter increases; at large jitter, close, far, and even overlapping arrangements can occur because no strict hard-core exclusion is imposed [2503.12438].

Hyperuniform and stealthy-hyperuniform patterns constitute a reciprocal-space definition of correlated disorder. With mean density $\rho$, pair correlation $g(r)$, and $h(r)=g(r)-1$, the static structure factor is
$$
S(\mathbf k)=1+\rho\int h(\mathbf r)e^{-i\mathbf k\cdot \mathbf r}\,d\mathbf r.
$$
Hyperuniformity requires $S(k\to 0)\to 0$, equivalently a number variance scaling slower than Poisson, while stealthy hyperuniformity imposes $S(\mathbf k)=0$ for $|\mathbf k|<K$ [2306.13540], [2110.12034], [1701.06799]. In metasurface design this suppresses long-wavelength density fluctuations and concentrates scattering into selected momentum channels.

A further structural regime appears when the notion of isolated meta-atoms begins to break down. In ultrathin dielectric metasurfaces patterned from square posts, a “critical packing regime” occurs when a significant fraction of metaatoms become physically connected; in the reported scanning electron micrographs this corresponds to about half of the metaatoms touching, with short chains and loops coexisting with isolated posts [2505.02244]. This regime is neither a dilute particulate topology nor a fully semi-continuous aggregate, and its optical response cannot be reduced to either limit.

## 2. Structural descriptors beyond conventional disorder parameters

For metasurfaces with weak or uncorrelated positional disorder, conventional descriptors such as $g(r)$, nearest-neighbor distributions, and Fourier-based metrics are often sufficient. Under strong correlations, however, these descriptors become less reliable as universal disorder coordinates. In strongly correlated Gaussian-distorted lattices, collective point displacements obscure the original lattice, broaden the accessible configuration space, and make the generative parameter $S_d$ a poor proxy for “how disordered” a sample actually is; finite-size windows and measurement noise further degrade statistical descriptors [2306.13540].

A topology-inspired alternative uses persistent homology on the point coordinates. A Vietoris–Rips filtration is built by drawing balls of radius $r$, connecting points whose balls intersect, and tracking topological features as $r$ grows. In two dimensions, the main features are connected components $H_0$ and loops $H_1$, summarized by Betti numbers $\beta_0(\epsilon)$ and $\beta_1(\epsilon)$ as functions of $\epsilon=2r$, or by the birth–death pairs $(b,d)$ forming a persistence diagram $D$ [2306.13540]. On this basis, two disorder descriptors were introduced:
$$
\mathrm{nSH}=\frac{1}{N}\sum_{(b,d)\in H_0\subset D}(d-b),
$$
and
$$
\mathrm{TD}=\sum_i\left[1+\frac{\sum_{(b,d)\in H_i\subset D}\left[\frac{d-b}{L_i}\right]\ln\!\left(\frac{d-b}{L_i}\right)}{\ln(\Omega_i)}\right],
$$
with $L_i=\sum_{(b,d)\in H_i\subset D}(d-b)$ and $\Omega_i$ the number of $H_i$ features. In the reported interpretation, $\mathrm{nSH}$ is a proxy for the typical nearest-neighbor scale, whereas $\mathrm{TD}$ is minimal for ordered lattices, invariant under global rescaling, independent of dataset size $N$, and orthogonal to $\mathrm{nSH}$ [2306.13540].

The computational pipeline used these descriptors on 1203 perturbed square-lattice patterns with $P\in\{500,600,700\}\,\mathrm{nm}$, $S_d\in[0,0.4]$, and $L_c\in\{0,2,8\}$, each containing $25\times25=625$ points. Vietoris–Rips filtrations were computed with Ripser; Wasserstein distances between persistence diagrams were computed with GUDHI and embedded by classical multidimensional scaling. The reported behavior was that uncorrelated disorder yielded clear clustering by $(P,S_d)$, intermediate correlation widened and overlapped clusters, and strong correlation merged them, making the generative parameters unreliable predictors. Coloring the same embeddings by $\mathrm{nSH}$ recovered the spacing scale across $L_c$, while coloring by $\mathrm{TD}$ recovered disorder strength across $P$ and $L_c$ [2306.13540].

Topological descriptors are not the only nontrivial structural measures in the field. In morphogenetically generated correlated media, a translational order metric
$$
\tau=\frac{1}{\tau_{\max}}\sum_{\mathbf k\neq 0}\left[S(\mathbf k)-S_p(\mathbf k)\right]^2
$$
was used to follow the evolution from correlated disorder to crystallization. The reported normalized values were $\tau_n\approx 0.03$ for early correlated disorder, $\tau_n\approx 0.13$ for a transitional regime, and $\tau_n\approx 0.98$ for a hexagonal crystal [2306.10931]. The broader methodological implication is that correlated-disordered metasurfaces are often better characterized by descriptors that remain reference-free under changes of pitch, correlation length, or morphology than by a single generative disorder amplitude.

## 3. Scattering theory, collective electrodynamics, and reciprocity-space design

The electromagnetic response of disordered metasurfaces is commonly decomposed into coherent and diffuse contributions. For the scattered field,
$$
\mathbf E_s=\langle\mathbf E_s\rangle+\delta \mathbf E_s,\qquad \langle\delta\mathbf E_s\rangle=0,
$$
and the averaged intensity separates as
$$
\langle|\mathbf E_s|^2\rangle=|\langle \mathbf E_s\rangle|^2+\langle|\delta\mathbf E_s|^2\rangle.
$$
Within an independent scattering approximation, the radiant intensity is often written as a single-particle form factor multiplied by a structure factor, while effective-field and quasi-crystalline closures incorporate collective interactions more self-consistently [2212.02541], [2310.14349]. For coherent specular reflection and transmission by a particle monolayer, the analytical expressions in the independent scattering approximation and the effective field approximation depend explicitly on the one-point density and the single-particle scattering amplitude; pair correlations enter the coherent channel only through more advanced closures such as QCA or through full-wave simulations [2212.02541].

In dense or resonant metasurfaces, discrete-element electrodynamics becomes essential. For asymmetrically split-ring arrays, each resonant arc is modeled as a damped oscillator and all recurrent scattering processes are retained through a many-body coupling matrix $\mathcal C$, so that the collective amplitudes obey
$$
\dot{\mathbf b}=\mathcal C\,\mathbf b+\mathbf F(t),
$$
with collective eigenvalues
$$
\lambda_j=-\gamma_j/2-i\delta\omega_j.
$$
This microscopic description captures radiative and non-radiative decay, retardation, electric–electric and magnetic–magnetic dipole interactions, and electric–magnetic cross-coupling [1812.10452]. The central physical point is that the response of a disordered metasurface is mapped from pair separations and multiple scattering sequences into collective linewidths, frequency shifts, mode localization, and far-field resonances.

For arrays supporting surface lattice resonances, periodic references remain useful. In a square lattice of pitch $a$ and effective refractive index $n_{\mathrm{eff}}$, the Rayleigh anomaly for diffraction order $(p,q)$ satisfies
$$
\lambda_{\mathrm{RA}}=\frac{n_{\mathrm{eff}}\,a}{\sqrt{p^2+q^2}}.
$$
Correlated disorder can emulate effective lattice vectors through short-range order and peaks in $S(k)$, thereby supporting SLR-like diffractive coupling without perfect periodicity [2306.13540].

A distinct analytical framework emerges when disorder itself is designed in reciprocal space. For lattices perturbed by correlated random displacements, the disorder statistics define three scattering components: a diffuse background, Bragg-like diffraction orders, and correlation halos. The halo term is absent for uncorrelated disorder, depends on the increment distributions $\Delta_m=\epsilon_{j+m}-\epsilon_j$, and can be positive or negative. In consequence, correlation halos are not broadened diffraction peaks; they are independent features whose positions depend on the correlation range and can remain visible after ordinary diffraction orders vanish [2602.17208]. This result expands the standard view in which disorder merely broadens reciprocal-lattice peaks into a broadband background.

## 4. Resonances, localization, and experimentally observed optical regimes

The reported optical consequences of correlated disorder are diverse because the relevant mechanisms differ across metasurface classes. In plasmonic nanoparticle lattices designed with topological descriptors, the strength of surface lattice resonances was correlated with $\mathrm{TD}$ while keeping $\mathrm{nSH}$ approximately fixed. In theory, for gold nanocylinders of height $50\,\mathrm{nm}$ and diameter $120\,\mathrm{nm}$ in a medium of refractive index $n=1.41$, three metasurfaces selected at the same $(L_c,S_d,P)=(8,0.3,500\,\mathrm{nm})$ but different $\mathrm{TD}$ showed progressively weaker resonances as $\mathrm{TD}$ increased, with reported quality factors $Q\approx 8.2$, $7.5$, and $6.5$ from lowest to highest $\mathrm{TD}$. In experiment, focused-ion-beam-fabricated elongated gold nanodisks coated with IC1-200 and measured at normal incidence showed that in 5 out of 6 correlated comparisons the metasurface with lower $\mathrm{TD}$, even when $S_d$ was larger, exhibited higher $Q$ [2306.13540].

| Pattern | TD | Reported $Q_{\parallel},Q_{\perp}$ |
|---|---:|---|
| periodic, $L_c=0,S_d=0$ | 0.000 | 10.1, 11.5 |
| $L_c=6,S_d=0.2$ | 0.030 | 4.0, 5.2 |
| $L_c=6,S_d=0.4$ | 0.012 | 9.3, 6.8 |
| $L_c=8,S_d=0.2$ | 0.025 | 6.7, 4.4 |
| $L_c=8,S_d=0.4$ | 0.005 | 7.8, 7.0 |
| $L_c=10,S_d=0.2$ | 0.026 | 8.0, 11.5 |
| $L_c=10,S_d=0.4$ | 0.002 | 10.1, 11.5 |

The main exception occurred for $L_c=10$ and perpendicular polarization, where both correlated samples had $Q_\perp\approx 11.5$ because large $L_c$ smooths short-range disorder while preserving long-scale distortions; $\mathrm{TD}$ detects the latter, whereas SLRs remain robust to it [2306.13540].

In metamaterial arrays with strong radiative interactions, positional disorder drives a different transition. In a regular $30\times 36$ array of asymmetrically split rings, approximately $70\%$ of the driven response concentrates into a single spatially extended many-body subradiant eigenmode with magnetic dipoles in phase and collective decay rate $\gamma\approx 0.21\Gamma$. A gradual increase of positional disorder rapidly localizes the mode and red-shifts the far-field transmission resonance through a cooperative Lamb shift; for one representative realization, $\gamma$ increased only weakly to $\approx 0.28\Gamma$ at $D\approx 0.33$, showing that localization and subradiance can coexist at moderate disorder [1812.10452].

Correlated disorder also modifies the density of states and localization in two-dimensional resonant media. For vector TE waves in stealth-hyperuniform point patterns, localization occurs at moderate density in the same window where the density of states exhibits a pseudo-gap; the reported localization island appears around $k_0a\approx 4.2$–$5.2$ and $\delta\approx 0$, whereas no signature of localization is found for white-noise disorder. For scalar TM waves, localization occurs at high density irrespective of correlations [2110.12034]. The proposed microscopic origin is destructive interference between independent scattering and recurrent loop scattering weighted by the short-range peak of $g_2(r)$ near $r\approx a$.

Near the critical packing threshold, connected dielectric metasurfaces display abrupt far-field changes tied to a redistribution of quasi-normal modes in the complex-frequency plane. In particulate arrays, the photon density of states contains clustered clouds and bandgap-like voids; at critical packing these voids begin to vanish and new collective resonances shift to lower frequencies; in aggregate regimes the PDoS becomes broad and nearly uniform across the visible, with $Q$ factors generally below $50$ [2505.02244]. In the far field, correlated particulate arrays show a blue shift of the diffuse BRDF maximum with increasing density, while transitions into critical and aggregate regimes produce broadband diffuse whitening. For $f=0.1$, $n=10\,\mu\mathrm m^{-2}$, and $l=170\,\mathrm{nm}$, the reported diffuse brightness approaches $\sim 80\%$ across the visible even though the structure is a single $145\,\mathrm{nm}$ poly-Si layer, and at critical packing a broadband specular “mirror” peak reaches $\sim 50\%$ near $\lambda\approx 670\,\mathrm{nm}$ with $R_s/T_s\approx 15$ [2505.02244].

Hyperuniform reciprocal-space engineering leads to yet another regime: isotropic annular scattering and emission. Gold metasurfaces derived from stealth-hyperuniform point sets with stealthiness $\chi=0.49$ and 4000 points exhibit a single broad isotropic diffraction maximum, and both scattering and fluorescence measurements show rotationally symmetric rings at in-plane momentum
$$
|k_{\mathrm{ring}}|\approx \kappa \frac{2\pi}{a}, \qquad \kappa\approx 1.03\text{–}1.09,
$$
so that $|k_{\mathrm{ring}}|/k_0\approx \kappa \lambda/a$. The opening of the ring therefore scales with the inverse correlation length parameter $a^{-1}$, and in fluorescence the ring closes near $a/\lambda\approx 0.9$ because of band folding into the light cone [1701.06799].

A separate application to infrared silicon photodetection uses a shuffled-lattice disordered metasurface integrated with upconversion nanoparticles. In that work the layout is generated by independent uniform jitter of Si nanopillars around a square lattice, and the paper does not report explicit correlation descriptors such as $g(r)$ or $S(k)$; nevertheless, the highly disordered configuration with $\sigma=0.7$ is reported to increase infrared absorption by $2.6$-fold and the near field by $3.9$-fold relative to the ordered structure, while the measured responsivity at $1550\,\mathrm{nm}$ reaches $0.22\,\mathrm{A/W}$ at room temperature, corresponding to an external quantum efficiency of $17.6\%$ [2503.12438]. The stated mechanism is disorder-induced mode packing of hybrid Mie–plasmonic cavities together with field enhancement that boosts Er$^{3+}$ upconversion and hot-electron generation.

## 5. Design methodologies, fabrication routes, and inverse design

A practical design workflow for correlated-disordered metasurfaces has been articulated most explicitly for topological learning. The sequence is: specify a target optical property, choose a disorder model such as minimum-distance-constrained patterns, Gaussian-correlated displacements, or hyperuniform and stealthy-hyperuniform patterns, generate candidate point sets over ranges of $S_d$, $L_c$, pitch, or minimum spacing, compute Vietoris–Rips persistence diagrams and the descriptors $\mathrm{nSH}$ and $\mathrm{TD}$, select patterns at matched $\mathrm{nSH}$ and desired $\mathrm{TD}$, predict the optical response with a dipole approximation or a learned regression, and then iterate until the target metrics are met before fabrication and validation [2306.13540]. In this formulation, $\mathrm{nSH}$ fixes local spacing while $\mathrm{TD}$ ranks positional disorder.

Morphogenetic design provides a different route to correlated disorder by replacing global optimization with local reaction–diffusion rules. In the Gray–Scott model,
$$
\partial_t A=d_A\nabla^2A-AB^2+f(1-A),\qquad
\partial_t B=d_B\nabla^2B+AB^2-(f+k)B,
$$
the reported parameter set $f=0.036$, $k=0.065$, $d_A=0.5$, $d_B=1$ generates self-replicating spots whose early-time patterns are stealthy-hyperuniform-like and isotropic in $S(\mathbf k)$, whose transitional states develop weak local crystallites, and whose late-time states converge toward a compact hexagonal crystal [2306.10931]. After thresholding at $0.45$, the reported supercells yield spot diameters of about $3\,\mathrm{mm}$ and average center-to-center spacing of about $5\,\mathrm{mm}$. Extruded into dielectric rods of $\varepsilon_r=3$ in air, these patterns produce isotropic microwave TM bandgaps with reported normalized widths of approximately $6\%$ in the correlated-disordered state, $7\%$ in the transitional state, and $8\%$ in the crystal [2306.10931]. The method is notable because the paper states that it eliminates cost-function minimization and is natively scalable to large domains.

Reciprocal-space engineering begins from a target $S(\mathbf k)$ rather than a real-space motif. In hyperuniform gold metasurfaces, a point pattern was chosen so that $S(\mathbf k)$ had a pronounced exclusion region near $k=0$ and a single broad isotropic maximum at $ak/2\pi\approx 1.03$; pillar-type and network-type layouts were then produced while preserving the dominant isotropic resonance [1701.06799]. In a more general correlated-noise framework, the averaged far-field intensity is expressed as the sum of a diffuse term, a Bragg term, and a correlation-halo term, with the displacement probability density function $\rho_\epsilon(\epsilon)$ and the correlation length $L_c$ as independent design knobs. The paper gives an explicit constructive example,
$$
\rho_\epsilon(\epsilon)=2\left(\frac{\sin(\pi \epsilon)}{\pi \epsilon}\right)^2\left[\cos(2\pi q\epsilon)+1\right],
$$
whose Fourier transform preserves only selected diffraction orders; for $q=3$, the retained orders are $k_x=0,\pm 3$ with reported intensities $1$ at specular and $0.25$ at $k_x=\pm 3$ [2602.17208]. The same paper describes the resulting framework as a practical method for inverse design, namely finding the disorder that produces desired scattering patterns.

Fabrication routes reflect the breadth of the field. Reported examples include focused-ion-beam fabrication of gold nanodisks on $\sim 12\times12\,\mu\mathrm m$ areas followed by spin-coating with IC1-200 [2306.13540], electron-beam lithography and reactive-ion etching of $145\,\mathrm{nm}$ poly-Si on fused silica with tunable connectivity [2505.02244], electron-beam lithography and ICP etching of Si nanopillars followed by $40\,\mathrm{nm}$ Al deposition and spin-coated NaYF$_4$:Er$^{3+}$ core–shell nanoparticles [2503.12438], electron-beam lithography of hyperuniform gold pillars and networks on glass followed by a $50\,\mathrm{nm}$ PMMA layer doped with DCM dye [1701.06799], and colloidal deposition of approximately $100\,\mathrm{nm}$ silver nanocubes above a SiO$_2$/Si reflector [2211.09520]. Large-area and bottom-up routes are likewise emphasized in the broader review literature because correlated disorder can be more fabrication-resilient than tightly periodic phase profiles [2310.14349].

## 6. Misconceptions, limitations, and research directions

A recurrent misconception is that disorder in metasurfaces is synonymous with uncontrollable whitish diffuse scattering. Several results contradict that reduction. Correlated disorder can support sharp or quasi-sharp collective resonances, isotropic momentum-space rings, specular suppression near $k\approx 0$, critical-packing mirror peaks, prescribed diffuse whitening, or Morpho-like correlation halos [2306.13540], [1701.06799], [2505.02244], [2602.17208]. Another misconception is that a single scalar disorder amplitude is always a sufficient disorder coordinate. Under strong positional correlations, $S_d$ can become ambiguous, and topological descriptors or reciprocal-space statistics become more informative than the generative parameter itself [2306.13540].

Model validity is strongly regime-dependent. In coherent specular theory, ISA is accurate for very dilute monolayers, and EFA remains quantitatively useful up to about $10\%$ coverage and large angles for high-index dielectric particles, but for plasmonic particles at comparable coverage strong near-field coupling and hot spots degrade both approximations [2212.02541]. In critical-packing dielectric metasurfaces, extended Maxwell–Garnett modeling matches particulate regimes but underestimates the strong reflection peak near the connectivity threshold because connected-cluster physics and collective electric–magnetic dipolar response are not captured by a local effective medium [2505.02244]. In SLR design, global $\mathrm{TD}$ can overestimate disorder that does not materially affect the resonance quality factor because long-scale distortions matter less than short-range coupling; a local $\mathrm{TD}$ computed in sliding windows was proposed as a better match to near-field interaction ranges [2306.13540].

The physical role of correlations is also not monotonic. Large $L_c$ can improve short-range order and preserve high-$Q$ resonances even while increasing global topological disorder; hyperuniformity can reduce low-angle scattering and broaden absorption; connected morphologies can whiten diffuse light but also erase spectral selectivity by filling PDoS voids [2306.13540], [2505.02244]. Polarization dependence remains important in several platforms, from elongated plasmonic nanodisks to TE/TM localization windows in resonant point patterns [2306.13540], [2110.12034]. In particular, the localization results in correlated two-dimensional media show that vector TE waves localize in a correlated moderate-density regime but exhibit no localization signature in white-noise disorder, whereas scalar TM waves localize at high density regardless of correlation class [2110.12034].

Current research directions accordingly emphasize descriptors and models that preserve structural universality while remaining tied to physically relevant coupling scales. Reported directions include localized topological descriptors aligned with coupling radii and integration with inverse design, topology optimization, and graph neural networks using topological features as priors or constraints [2306.13540]; nonlocal effective models and polarization control for critical-packing topologies [2505.02244]; and broader k-space design of hyperuniform, quasi-periodic, or multilayer correlated-disordered metasurfaces for broadband antireflection, structural color, transparent displays, chiral films, light trapping, and wavefront manipulation [2310.14349]. The field’s unifying idea is that disorder becomes most useful when it is specified statistically, topologically, or in reciprocal space rather than treated as an uncontrolled deviation from a crystal.

Source: https://www.emergentmind.com/topics/correlated-disordered-metasurfaces