Metadisorder for designer light in random systems

See allHide authors and affiliations

Science Advances  14 Oct 2016:
Vol. 2, no. 10, e1501851
DOI: 10.1126/sciadv.1501851


Disorder plays a critical role in signal transport by controlling the correlation of a system, as demonstrated in various complex networks. In wave physics, disordered potentials suppress wave transport, because of their localized eigenstates, from the interference between multiple scattering paths. Although the variation of localization with tunable disorder has been intensively studied as a bridge between ordered and disordered media, the general trend of disorder-enhanced localization has remained unchanged, and the existence of complete delocalization in highly disordered potentials has not been explored. We propose the concept of “metadisorder”: randomly coupled optical systems in which eigenstates can be engineered to achieve unusual localization. We demonstrate that one of the eigenstates in a randomly coupled system can always be arbitrarily molded, regardless of the degree of disorder, by adjusting the self-energy of each element. Ordered waves with the desired form are then achieved in randomly coupled systems, including plane waves and globally collective resonances. We also devise counterintuitive functionalities in disordered systems, such as “small-world–like” transport from non–Anderson-type localization, phase-conserving disorder, and phase-controlled beam steering.

  • Disorder
  • Anderson localization
  • wave transport
  • metamaterial
  • network
  • random-walk potential
  • invisibility
  • beam steering
  • chirality
  • orbital angular momentum


Network modeling has provided an intuitive picture to understand various complex systems in nature and society. In these system networks, disorder is a crucial factor in signal transport over the system. For example, in the pioneering work in graph theory (1), D. J. Watts discovered “small-world disordered networks” with randomly rewired connections, which can model various elemental systems in physics, biology, and sociology: seismic networks (2), Caenorhabditis elegans neurons (3), brain connectome (4), and affinity groups in social networks (1).

The role of disorder is also evident in wave physics, especially when compared to the “order” in potential energy, such as periodic or quasiperiodic potentials. Periodic potentials allow ballistic transport through extended Bloch eigenstates (5), whereas broken correlation in disordered potentials leads to Anderson-localized eigenstates (68), which significantly suppress wave transport. Although serious attempts have been made (812) to fill the gap between the order and disorder in wave systems, the increase of disorder has only led to a monotonic change from extended to localized eigenstates (9), prohibiting the existence of eigenstates with unusual forms: for example, completely delocalized eigenstates in highly disordered potentials. Thus, a clear distinction has been maintained between the applications of order and disorder in wave systems, in accordance with the contrast between their eigenstates: transporting devices using ordered potentials (5, 13) and focusing devices using disordered potentials (1416). In disordered systems, the design of nontrivial waves such as globally collective and scattering-free propagations remains a challenge.

Here, we demonstrate the existence of globally collective and delocalized waves in randomly coupled optical systems, by introducing the concept of “metadisorder.” By using platforms of coupled waveguides (or resonators), we prove that designer mode excitation and wave propagation can be achieved in the random network by controlling the self-energy of the optical elements (for example, size and loss and gain of waveguides) (Fig. 1). Distinct from random scattering and Anderson localization in classical random systems (Fig. 1, A and B), unusual phenomena are demonstrated in the “artificial” disordered optical potentials (Fig. 1, C and D): perfect plane waves without any scattering or phase distortion, designer guided waves, globally collective resonances, conservative waves in complex potentials, and the invisible disorder of phase conservation, all of which provide ordered waves with disorder-like energy bands. By creating artificial disorder to achieve non–Anderson-type localization, we also reveal the separate control of eigenstate localization and wave transport, analogous to the separation of transport and clustering in the “small-world network” (1). Our method, paving the way toward disorder-robust small-world optics, is also applied to functional wave devices in randomly coupled optical systems, including tunable focusing, in-phase spatial oscillation, parity converters, point-source excitation of plane waves, and modal filtering from incoherent light sources.

Fig. 1 The concept of metadisorder systems.

(A and B) Random scattering and Anderson localization in randomly coupled systems composed of identical waveguides (constant widths w) (A) or resonators (constant radii r) (B). (C and D) Delocalized propagations of designer modal profiles in metadisorder systems. The control of self-energy distributions can be obtained by manipulating the waveguide widths wp,q (C) or the resonator radii rp,q (D). The positions of optical elements in coupled systems are assumed to be completely random. Blue dotted arrows denote propagation directions.


Wave transport and localization in randomly coupled systems

We start with random systems composed of weakly coupled optical elements, such as waveguides (5, 17) or resonators (18, 19). On the basis of discrete models of coupled-mode theory (CMT) (17, 20) or tight-binding (TB) analysis (19), an N-body system is governed by the eigenvalue equation Hψ = γψ (see Materials and Methods). The Hamiltonian H can be decomposed into H = D + K, where D is the diagonal matrix for the self-energy of each optical element [for example, wavevector β of each waveguide or resonant frequency f of each resonator, both of which define the phase evolution of the field in uncoupled optical elements (20)], and K is the off-diagonal matrix for the interaction energy between elements [for example, coupling coefficient κ between coupled waveguides or between coupled resonators (20)], which represents the network of the system. Figure 2A shows a one-dimensional (1D) randomly coupled system with off-diagonal disorder (21), having the identical self-energy β and disordered interactions {[D]p = γo0 and [K]pq = κpq = κ0 + Δκ·u(−1,1), where κpq is the coupling between the pth and qth elements for 1 ≤ p, qN, κ0 and Δκ represent the averaged and disordered coupling, respectively, and u is the uniform probability density function}.

Fig. 2 Effect of disorder in optical systems.

(A) A schematic of a random optical system composed of coupled waveguides, analogous to the randomly coupled pendulums with identical oscillating features. The oscillation of each pendulum describes the phase evolution during propagation. (B to D) The first three eigenstates for ordered (B), weakly disordered (C), and Anderson (D) potentials, calculated by using CMT (see Materials and Methods). The potential n denotes the effective waveguide index of a single waveguide. (E to G) Corresponding wave transport calculated by using the transfer-matrix method (TMM) (see Materials and Methods). (H and I) The variations of modal size w (H) and diffusion exponent α (I) are shown as a function of the disorder Δκ. (J) Eigenbands (neff) for the disorder in (B) (black dotted line) and (D) (green points and line). The points in (H) to (J) represent each statistical ensemble, and the solid lines are the averages of 200 ensembles. The green dotted line in (I) denotes the diffusion state (α = 1). γo0 = 1.6·k0, κ0 = 0.01·k0, and N = 51 for (B) to (J), where k0 = 2π/λ0 is the free-space wave number. The practical waveguide design and the distance between waveguides in the mid-infrared regime (λ0 = 3 μm) are calculated in note S1 using COMSOL Multiphysics.

Figure 2 (B to D) presents a few eigenstates of the optical systems at different degrees of off-diagonal disorder (21) (see note S1 for its practical realization in the mid-infrared regime). As the strength of the disorder Δκ increases, Bloch eigenstates (Fig. 2B) begin to be localized (Fig. 2C), eventually exhibiting the wavelength-scale Anderson localization of exponential decay envelopes (Fig. 2D). The localization naturally results in the suppression of wave transport within the network, from ballistic (Fig. 2E; α = 2) to diffusive (Fig. 2F; 1 < α < 2) and even subdiffusive (Fig. 2G; α < 1) transport [α, diffusion exponent (22); see Materials and Methods]. Not restricted to eigenstate localization (Fig. 2H; w, modal width; see Materials and Methods and note S2) and the following forms of suppressed wave transport (Fig. 2I), the increase of the disorder in the system also alters the spreading of its eigenvalues, linearizing the eigenband (Fig. 2J and note S3). In these results, a continuous transition between the regimes of order and disorder is evident, confirming the classical relationship (612) between localization, transport, and disorder. Motivated by the generic form of the Hamiltonian H = D + K, we now demonstrate that the design of unconventional eigenstates in highly disordered potentials can be achieved by using the degree of freedom on the self-energy of each element in D, which was neglected in Fig. 2.

Molding of an eigenstate in randomly coupled systems

Suppose that we desire to “mold” an eigenstate ψ with the spatial distribution of vm = [vm1, vm2,…, vmN]T with the eigenvalue γm while preserving the randomly coupled network of the system. We design such a system in the reciprocal space basis; one of its basis vectors has the spatial distribution of vm. For this purpose, we develop the eigen-decomposition matrix V = [vm, v2,…, vN] using the Gram-Schmidt process, where vm and the set of column vectors vs = [vs1, vs2,…, vsN]T (s = 2, 3,…, N) together compose the orthonormal basis set. Because of the orthonormality (VV = I), the eigenvalue equation in the V-reciprocal space becomes VHV(Vψ) = γ(Vψ), or Hrψr = γψr, where Hr = VHV and ψr = Vψ. The randomly coupled network of the system is represented by a complex matrix Kr = VKV in the V-reciprocal space, and then, the V-reciprocal Hamiltonian Hr = VDV + Kr has the following componentsEmbedded Image(1)where γop = [D]p and Krpq = [Kr]pq. To design the eigenstate ψ of spatial representation vm, one of the reciprocal eigenstates should be ψr = [1, 0,…, 0]T. This condition is uniquely fulfilled when the first column of Hr has only one nonzero component of [Hr]11, and its eigenvalue can always be set as desired by assigning [Hr]11 = γm. We then achieve the self-energy of each element γop deterministically as follows (see Materials and Methods), from [[Hr]11, [Hr]21,…, [Hr]N1]T = [γm, 0,…, 0]TEmbedded Image(2)or, simply, γo = [γo1, γo2,…, γoN]T = [diag(vm)]−1VmKr11, −Kr21,…, −KrN1]T.

Equation 2 indicates that if the inverse of the matrix diag(vm) exists, that is, vmi ≠ 0 for all i, the self-energy vector γo can always be found, determining the optical potential of each element from γop (see fig. S1C for γo = n). Therefore, for any network K regardless of the degree of disorder, a single eigenstate can always be molded into the desired shape of vm with the eigenvalue γm by adjusting the potential of each element (Fig. 2A versus Fig. 3A) while preserving the network of the system, which we herein call metadisorder to describe the artificial change (engineered D) of the disordered system (random K). The proposed metadisorder system allows for the nontrivial form of an eigenstate in all regimes of disordered networks K, for example, the globally collective eigenstate, in contrast to the case of identical elements (Fig. 2, C and D). We note that such a degree of freedom for an eigenstate in an N-body system originates from the N number of the adjustable design parameters, the self-energy of N elements, in an N × N Hamiltonian matrix.

Fig. 3 1D metadisorder systems.

(A) A schematic of a 1D metadisorder system composed of coupled waveguides, analogous to the randomly coupled pendulums with different self-oscillating features, such as oscillating period (rod length) and gain or loss parameters (color). Each waveguide has different real parts of self-energy due to changing the width of the waveguide (note S1). The colors of the waveguides represent the imaginary part of self-energy: gain and loss (treated in note S4). (B to I) Designed eigenstates and optical potentials, eigenstate propagations (B to E), and eigenvalues (neff) (F to I) of 1D metadisorder systems are calculated by using the CMT for plane wave (vm(x) = 1) (B and F), Gaussian wave (vm(x) = exp[−x2/(2·σ2)]) (C and G), interface wave (vm(x) = exp[−|x|/(2·σ)]) (D and H), and surface wave (vm(x) = exp[−|xxL|/(2·σ)]) (E and I) eigenstates, where the spatial bandwidth σ = Lst/16 in (C) to (E) and (G) to (I), the left boundary xL = −Lst/2, and Lst is the overall potential length. Δκ = κ0 in (B) to (E) for the extreme degree of disorder. Blue symbols represent Δκ = κ0, green symbols represent Δκ = 0.53·κ0, and black dotted lines represent Δκ = 0 in (F) to (I). All of the other parameters are the same as those in Fig. 2, based on note S1.

1D metadisorder systems for designer waves and non-Anderson localization

Figure 3 shows examples of designer waves in 1D metadisorder systems, where the optical potential of each waveguide is calculated using Eq. 2. Compared to the highly disordered system composed of identical waveguides in Fig. 2, we achieve various examples of wave systems having a designer eigenstate vm: plane wave (Fig. 3B), Gaussian-enveloped guided wave with nonexponential decay (Fig. 3C), and interface (Fig. 3D) and surface (Fig. 3E) waves both with Anderson-like (6, 9, 10, 21) exponential decay. We note that the real-valued vm corresponds to the designer ground state in the disordered eigenband (red arrows in Fig. 3, F to I; see note S4 for designer excited states, which allow conservative waves in complex potentials).

Our design method allows for a scattering-free plane wave in highly disordered systems (Fig. 3B; Δκ = κ0), in stark contrast to conventional disordered systems (Fig. 2G), which only lead to strong localization from random scattering. The plane wave eigenstate in highly disordered systems, which has a modal size equal to the overall system size, is more extended than that in finite-N ordered systems with the field discontinuity at the boundary due to the broken translational symmetry (Fig. 2B; see also note S5 for the stability of the delocalization in metadisorder systems, against self-energy and coupling errors). Furthermore, unconventional localization forms, such as the non-Anderson Gaussian localization (Fig. 3C) or the designed Anderson-like exponential localization at the interface (Fig. 3D) or surface (Fig. 3E), can also be achieved as a general extension to the accidental emergence of classical Anderson localization (Fig. 2D). Because a potential with a globally extended eigenstate should have the reduced random scattering in the overall system, other eigenstates of similar eigenvalues also tend to have wider spatial bandwidth. We note that the designed eigenstates in metadisorder systems can be excited through the evanescent coupling (23) from a single waveguide (note S6).

The concept of metadisorder also creates a new class of disordered potentials, which support the counterintuitive relation between eigenstate localization (w) and transport (α) by imposing the designer eigenstate vm with unconventional localizations. In Fig. 4, we consider the localized designer eigenstate of the specific form vm(x) = exp[−|x|g/(2·σg)] (Fig. 4A), where g = 1 for Anderson-like exponential localization and g = 2, 4, and 6 for the convenient examples of non-Anderson localizations. Figure 4 (F to N) presents the localization-transport (w-α) relation of 1D non-Anderson metadisorder systems compared to classical Anderson disorder (Fig. 4B) and Anderson-like metadisorder (Fig. 4, C to E). Although Anderson-like metadisorders (g = 1) provide a similar w-α relation to that of Anderson disorder (Fig. 4B versus Fig. 4, C to E), the non-Anderson metadisorder (g > 1) enables more “localized” waves yet achieves ballistic transport (for example, Fig. 4B versus Fig. 4F; a factor of ~2 decrease for w, and α ≈ 2, for Δκ < κ0/10), analogous to the separate control of localization and transport in “small-world” networks (1). Such a nonclassical wave transport even enables “disorder-induced” wave transport [increase in α for smaller values of w(Δκ) or σ; Fig. 4H], showing the reversed relation between w and α compared to that of conventional disordered systems (710, 17). This anomalous relation is more apparent for metadisorder systems with larger g (Fig. 4, I to N; for g = 4 and 6), allowing not only separate control of localization and wave transport with g and σ but also robustness of wave transport to the disorder Δκ, analogous to the difference between clustering and characteristic path length in small-world networks (1).

Fig. 4 Non-Anderson metadisorder systems with disorder-induced wave transport.

(A) Shapes of designer eigenstates vm(x) = exp[−|x|g/(2·σg)] with different g. (B to N) Eigenstate localization (w) and wave transport (α) for Anderson disorder with identical elements (B), g = 1 Anderson-like metadisorders (C to E), and g = 2 (F to H), g = 4 (I to K), and g = 6 (L to N) non-Anderson metadisorders . In metadisorder systems, the bandwidths of designer eigenstates are σ = Lst/16 (C, F, I, and L), σ = Lst/24 (D, G, J, and M), and σ = Lst/32 in (E, H, K, and N). Error bars denote the SD of 200 ensembles in (B) to (N). All of the other parameters are the same as those in Fig. 2, based on note S1.

Extension to 2D metadisorder systems for designer wave functionalities

The eigenstate design in a V-reciprocal space allows for its extension to multidimensional problems in a straightforward manner by including all of the coupling coefficients in a multidimension in the network matrix K. Here, we consider 2D disordered systems, obtained by the random deformation of the periodic lattice (Fig. 5A). Figure 5B shows an example of disordered coupled resonator systems from the random deformation of a 17 × 17 square lattice (see note S7 for the practical realization of 2D coupled resonator systems in the terahertz regime).

Fig. 5 2D metadisorder systems.

(A) Schematics of a coupled resonator lattice with identical elements (left) and its metadisorder-transformed structure (right). The self-energy of each element in the metadisorder system is adjusted by controlling the size of the resonator or using gain or loss materials. Nearest-neighbor (yellow arrows) and next-nearest-neighbor (green arrows) couplings are presented. (B to M) 2D metadisorder systems. The practical weak-coupling design for the CMT analysis is calculated in note S7, using COMSOL Multiphysics, for the terahertz regime (λ0 = c/f0 = 265.3 μm). For the periodicity of ax = ay = 0.16·λ0, the position of each resonator is randomly shifted by Δx = Δy = 0.03·λ0·u(−1, 1) in (B) to (M), and the self-energy of each element is adjusted by f0·(1 + Δf) following the metadisorder design from Eq. 2. (B) A sample of the obtained resonator distribution for the 17 × 17 lattice in (C) to (H). (C to E) Standing-wave collective resonances using real potentials for uniform (C), quadrupole (D), and exponentially localized (E) distributions. (F to H) Traveling-wave chiral collective resonances, using complex potentials, for dipole (F), quadrupole (G), and octopole (H) distributions. (I to M) Metadisorder-based functionalities in the 17 × 17 lattice for invisible disorder (I), steered focusing (J), spatial oscillation (K), parity-converted beam splitter (L), and the point-source excitation of oblique plane waves (M). The fields in input and output waveguides (I to M) are multiplied by 40 for better visualization. The designed resonance is set to γm = f0 for all cases. See movies S1 to S11 for the temporal dynamics of each case in (C) to (M).

Figure 5 (C to H) and movies S1 to S6 show collective and ordered light behaviors in highly disordered systems, which support strongly correlated phase information over the entire system. By adjusting each resonant frequency of constituent resonators following Eq. 2, we design free-form standing-wave resonances with a perfectly uniform field distribution (Fig. 5C), quadrupole phase distribution (Fig. 5D), and designed localization (Fig. 5E), despite randomly deformed interactions between resonators. Furthermore, the introduction of complex potentials allows for one-way traveling-wave resonances, for example, by imposing the form of exp(−iρθ) on vm (Fig. 5, F to H, for the azimuth θ). Such a “chiral” rotation of the phase in collective resonances derives the orbital angular momentum (OAM) (24) of resonant light. Notably, although the chiral feature of light in the proposed metadisorder systems also requires complex optical potentials with gain and loss, our method involves neither parity-time symmetry nor periodicity (25). Propagating light with nonzero OAM can also be achieved by using waveguide-based disordered systems.

Finally, having demonstrated collective resonance modes in highly disordered systems, we present the excitation of designer eigenstates (Fig. 5, I to M, and movies S7 to S11) with the external coupling of conventional waveforms. When the inner connection of the system is sufficiently strong, the separation of eigenvalues (that is, free spectral range) becomes sufficiently large to achieve wave dynamics dependent solely on a single eigenstate. Figure 5I and movie S7 show the wave flow through the perfectly uniform collective eigenstate over the entire system. Following the property of the eigenstate, the disordered system becomes “invisible” for incident plane waves, prohibiting any alteration of phase and amplitude [transmission (T), ~100%; zero effective index]. With eigenstate-based metadisorder design, we also implement high-level functionalities with excellent throughputs, including tunable light focusing (Fig. 5J; T, ~96%), phase-conserved spatial oscillation (Fig. 5K; T, ~98%), parity converters (even to odd) (Fig. 5L; T, ~99%) of real potentials, and the point-source excitation of oblique plane waves (5.6°) (Fig. 5M and movie S11; T, ~97%) using complex potentials. Because of the spectral selectivity of the designed eigenstate, metadisorder systems also allow the modal filtering function, or the excitation of target waveforms from spatially incoherent incidences (note S8).


To summarize, we revealed a new class of random wave systems, that is, the metadisorder, which can be globally collective and deliberately controlled. Exploiting metadisorders of non–Anderson-type localization, we first derive the counterintuitive relation between localization and transport, including small-world–like (1) or disorder-induced transport. As demonstrated in collective wave dynamics and functionalities, our eigenstate-based approach also provides a powerful means to control wave flow while preserving or manipulating the phase information. Although we controlled only the self-energy D for the Hamiltonian H = D + K, our method can also be easily extended to determine the “network” K of disordered self-energy distributions for the designer eigenstate, revealing the unexplored regime between order and disorder from the correlation between diagonal (or self-energy) and off-diagonal (or interaction energy) disorders. We emphasize that our method is distinct from other approaches to handling the flow of waves: the supersymmetric technique (23, 2628) controls eigenspectra but transforms eigenstates in a fixed manner, and the transformation optics technique (29) treats a “continuous” potential landscape, lacking the degree of freedom in interaction energy. From their small-world–inherited disorder robustness and globally collective features, we also envisage the application of metadisorder systems to many other nontrivial physics, such as hyperuniformity (30, 31), topological networks (18), or quasiparticles in disordered potentials.


Analysis of randomly coupled discrete optical systems

Consider the N-body system composed of weakly coupled optical elements. In CMT (17) or TB (19) methods, the governing equation, including self-energy and interaction energy, becomesEmbedded Image(3)where p = (1, 2,…, N) is the element number, ψp is the field at the pth element, ξ is the wave evolution axis [time t for coupled resonators (13) and space x, y, or z for coupled waveguides (17)], γop is the self-energy of the pth element, and κpq is the coupling coefficient between the pth and qth elements. For the steady-state solution (∂ξiγ), Eq. 3 becomes the matrix eigenvalue problem Hψ = γψ, whereEmbedded Image(4)

ψ = [ψ1, ψ2,…, ψN]T, D is the diagonal self-energy matrix, and K is the off-diagonal network matrix. The randomly coupled system can then be described by assigning random numbers to the components of the K matrix. Because each element number p corresponds to the physical location of the pth element Xp [for example, pxp for 1D and p → (xp, yp) for 2D problems], the obtained eigenstate ψ can be reexpressed in the spatial domain ψ = ψ(X).

Calculation of the diffusion exponent α

Consider the 1D system of Fig. 2 (ξ = y), which has an eigenstate of ψk(x) and corresponding eigenvalues γk (k = 1, 2,…, N; γk is the effective wavevector of ψk). For the incidence of φi(x) = Σak·ψk, the propagating field can be obtained through the TMM as φ(x,y) = Σak·ψk·exp(iγk·y). To analyze the transporting feature of the system without boundary effects, the incident wave was excited at the center waveguide [x = xm, where m = (N + 1)/2 for odd N]; we then calculated the spatially varying mean-square displacement (MSD) (22) M(y) as followsEmbedded Image(5)

When the MSD M(y) was fitted for y exponentially as M(y) ≈ cα·yα, we achieved the diffusion exponent α: α = 2 for ballistic transport and α = 1 for diffusive transport (22). The calculated results are shown in Figs. 2 (E to G and I) and 4. Note that m does not have to be the center waveguide precisely when N is sufficiently large, and thus, the boundary effect can be neglected.

Calculation of modal size

For the 1D system (ξ = y) with ψk(x) and γk, the modal size for each eigenstate is defined as (9)Embedded Image(6)where Δxp is the size of the pth element, obtained from the distance between waveguides for each value of the coupling coefficient (see note S1). See also note S2 for the eigenstate-dependent localizations.

The derivation of Eq. 2

The relation [[Hr]11, [Hr]21,…, [Hr]N1]T = [γm, 0,…, 0]T in the main text, for the design of ψ of spatial representation vm, has the following explicit formEmbedded Image(7)

Equation 7 can be recast into the form ofEmbedded Image(8)

Solving Eq. 8 for γo = [γo1, γo2,…, γoN]T, we arrived at Eq. 2 in the main text, with the condition of vmi ≠ 0 for all i.


Supplementary material for this article is available at

note S1. Design of waveguide-based 1D random networks.

note S2. Eigenstate-dependent localization.

note S3. The variation of the eigenband by disorders.

note S4. Designer excited states in 1D metadisorder systems.

note S5. Stability of the delocalization in metadisorder systems.

note S6. The excitation of the designer eigenstate in metadisorder systems.

note S7. Design of resonator-based 2D random optical networks.

note S8. Modal filtering through 2D metadisorder structures.

fig. S1. The design of waveguide-based elements for 1D random networks.

fig. S2. The localization of each eigenstate.

fig. S3. The variation of the eigenband for different degrees of disorder Δκ.

fig. S4. 1D metadisorder systems with complex potentials.

fig. S5. The variation of the ground-state modal size.

fig. S6. The excitation of designer eigenstates through the evanescent coupling.

fig. S7. The design of resonator-based elements for 2D random networks.

fig. S8. Modal filtering through 2D metadisorder systems, from incoherent source.

movie S1. Temporal wave dynamics of Fig. 5C.

movie S2. Temporal wave dynamics of Fig. 5D.

movie S3. Temporal wave dynamics of Fig. 5E.

movie S4. Temporal wave dynamics of Fig. 5F.

movie S5. Temporal wave dynamics of Fig. 5G.

movie S6. Temporal wave dynamics of Fig. 5H.

movie S7. Temporal wave dynamics of Fig. 5I.

movie S8. Temporal wave dynamics of Fig. 5J.

movie S9. Temporal wave dynamics of Fig. 5K.

movie S10. Temporal wave dynamics of Fig. 5L.

movie S11. Temporal wave dynamics of Fig. 5M.

References (3246)

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: We thank S. Torquato for his encouragement on our research on disordered systems, and we thank the anonymous reviewers for their valuable comments and insights. Funding: This work was supported by the National Research Foundation of Korea (NRF) through the Global Frontier Program NRF-2014M3A6B3063708 and the Global Research Laboratory Program K20815000003, which are all funded by the Ministry of Science, ICT and Future Planning of the Korean government. S.Y. was supported by the Basic Science Research Program (2016R1A6A3A04009723) through the NRF, funded by the Ministry of Education of the Korean government. Author contributions: S.Y. conceived the presented idea, developed the theory, and performed the computations. X.P. contributed to the analysis of wave transport. X.P. and J.H. verified the analytical methods. N.P. encouraged S.Y. to investigate the inverse design and small-world nature of eigensystems while supervising the findings of this work. All authors discussed the results and contributed to the final manuscript. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

Stay Connected to Science Advances

Navigate This Article