Transport Phase Diagram and Anderson Localization in Hyperuniform Disordered Photonic Materials

Transport Phase Diagram and Anderson Localization in Hyperuniform Disordered Photonic Materials


Hyperuniform disordered photonic materials (HDPM) are spatially correlated dielectric structures with unconventional optical properties (1); (2). They can be transparent to long-wavelength radiation while at the same time have isotropic band gaps in another frequency range (3); (4). This phenomenon raises fundamental questions concerning photon transport through disordered media. While optical transparency is robust against recurrent multiple scattering, little is known about other transport regimes like diffusive multiple scattering or Anderson localization (5). Here we investigate band gaps, and we report Anderson localization in two-dimensional stealthy HDPM using numerical simulations of the density of states and optical transport statistics. To establish a unified view, we propose a transport phase diagram. Our results show that, depending only on the degree of correlation, a dielectric material can transition from localization behavior to a bandgap crossing an intermediate regime dominated by tunneling between weakly coupled states.

Light propagation through a dielectric medium is determined by the spatial distribution of the material. Photons scatter at local variations of the refractive index. For a periodically organized system, interference dominates light transport and is responsible for optical phenomena in opal gems and photonic crystals (6). In random media, transport becomes diffusive through successive scattering events. The characteristic length scale over which isotropic diffusion takes place is the transport mean free path. Material thicker than the mean free path appears cloudy or white. However, when scattering centers are locally correlated, diffraction effects can be significant. The description of light transport then becomes a challenging problem with many applications, such as the transparency of the cornea to visible light (7), the strong wavelength dependence of the optical thickness of colloidal suspensions (8) and amorphous photonic structures (9); (10), and structural colors in biology (11). Critical opalescence and the relatively large electrical conductivity of disordered liquid metals (12) are closely related phenomena.

In the weak scattering limit photon transport is diffusive and can be described by a local collective scattering approximation, which states that the mean free path is inversely proportional to the number density of scatterers and to the effective transport cross section (12); (13); (14)


Here is the differential cross section for an isolated scatterer and is the momentum transfer. This equation relates positional correlations of the optical medium to transport via the structure factor . In the past, the local collective scattering approximation has been applied to dense and strongly scattering media (14), however the validity of this approach is limited as it does not include near-field corrections and recurrent scattering (15); (16). Clearly, the recently discovered isotropic bandgaps (3); (17) in HDPM cannot be derived from Eq. 1.

For uncorrelated or fully random media with optical transport properties are well understood. According to the single parameter scaling (SPS) hypothesis one expects a transition from diffuse scattering to Anderson localization (18); (19). Statistical properties of transport are governed by a single parameter that can be expressed as the ratio of a characteristic (localization) length to system size . While for transport is diffusive, SPS predicts a crossover to the Anderson localization regime for both disordered wires (quasi-one-dimensional system) and two-dimensional systems with any amount of disorder.

For correlated disordered media, over the last years, attention focused on the emergence of optical transparency (1); (4) and on photonic bandgaps (PBGs) in two-dimensional high-refractive-index disordered materials (3); (17); (20). In particular, the concept of stealthy hyperuniformity as a measure for the hidden order in amorphous materials has drawn significant attention (1); (2). HDPM are disordered, but uniform without specific defects, and the density of states can be strictly zero. A key parameter controlling structural correlations of HDPM, and thus the band gap width, is the stealthiness parameter , defined as the ratio between the number of constrained degrees of freedom to the total number of degrees of freedom (1). It is currently not known what effect has on the statistics of wave transport in HDPM outside the gap and whether the transport properties in HDPM can be understood in terms of SPS. Neither is anything known about the transport properties for frequencies near the band edge, a regime where disordered crystalline materials exhibit a rich and complex transition towards Anderson localization (21). The main goal of our work is to address these fundamental questions and to propose a phase diagram for wave transport through ideal hyperuniform disordered dielectric materials in two dimensions.

A relatively simple example of HDPM is a collection of infinite parallel cylinders of high-refractive index that are distributed according to a two-dimensional stealthy hyperuniform point pattern (SHU, see the insets in Fig. 1). The surface area and the number of points define a characteristic length through the number density . When the electric field is parallel to the cylinder axes (TM polarization), the propagation of light with a wave vector perpendicular to the axes can be described by a two-dimensional Schrödinger-like scalar equation with identical scatterers. This setup provides an ideal laboratory to explore the statistical properties of wave transport in two-dimensional HDPM. The normalized density of states (NDOS) as a function of frequency for a two-dimensional SHU pattern of monodisperse high refractive index cylinders has been reported for different values of  (3); (17).

Based on the NDOS and the results for two-dimensional disordered crystals (22); (23), we first make a hypothesis about the transport phase diagram of two-dimensional HDPM in Fig. 1. For strongly correlated, nearly crystalline materials (high ), stealth (transparent) and gap intervals are adjacent to each other, reminiscent of the conduction band and stop gap for electrons and photons in crystals. Lowering the degree of hyperuniformity to , the gap and stealth intervals shrink. In the newly accessible intervals in-between the stealth and the gap regions, the material possesses a reduced density of states and displays significant scattering. We expect a tunneling-regime (non-SPS) for frequencies near the full band gap and a diffusive scattering regime further apart. The tunneling regime is intimately related to the pseudo gap in which the density of states is low, but finite. Once the density of states is sufficiently high, far away from the gap, diffusive transport sets in and a crossover to classical Anderson localization should occur as . For high frequencies we also expect diffusive transport, independent of , since in this spectral region the scattering strength decays below the threshold where a gap opens.

Although the calculation of the standard transport mean free path can be useful to describe transport processes in disordered correlated media in frequency ranges where the system is almost transparent or transport is diffusive (14); (4), it loses its meaning in a region where the density of states is zero. Identifying the different transport regimes requires the numerical solution of the full multiple scattering problem in a wide spectral range. We calculate the decay of the intensity of the wave fields along the propagation direction and the full statistics of wave transport. First, we generate a statistical ensemble of 1000 spatial patterns of monodisperse high refractive index cylinders for each value of as described in (17). Our spatial two-dimensional point patterns show stealthy hyperuniformity, which means that not only long-range density fluctuations are suppressed, , but the structure factor vanishes over a finite range, . The critical wave number sets the phase boundary of the stealth phase (4); (17). With these parameters a PBG opens at a central frequency , and the PBG expands to a maximum width of  (17). Following prior work on disordered crystals (22); (23) we introduce the length


as an effective measure of the mean distance between states. Fig. 2(a) shows a map of . We can clearly identify the gap with no states and infinite (white region bounded by dashed line) and a pseudo gap region where the NDOS is significantly reduced, .

To obtain the characteristic decay length and to sample the statistical properties of transport, we apply the well-known generalized scattering matrix (GSM) method (24) to a system of fixed width as a function of the thickness in propagation direction. Periodic boundary conditions (inset in Fig. 3) define a set of transversal propagation channels. The optical analog of electrical conductance can be computed from Landauer’s formula. If is the intensity transmitted from incoming channel to outgoing channel , then the conductance of the system is . The exponential decay length of the conductance is determined using a fit to  (25). Fluctuations are suppressed by performing an ensemble average over 1000 realizations for each frequency and stealthiness. Fig. 2(b) shows a map of . The solid line corresponds to the boundary . Based on studies of two-dimensional disordered crystals (23), we expect this line to indicate the boundary between non-SPS and SPS regions.

So far we identified regions with exponential decay of the conductance but did not reveal the physical mechanism responsible for this rapid decay. To this end we analyze the optical transport statistics (5). Fig. 3 shows a color map of fluctuations in the logarithmic conductance for transmission through the whole system. As expected, regions associated with a PBG do not show substantial fluctuations despite a small decay length signaling transport by direct tunneling through the whole sample. Outside the gap, large fluctuations are most pronounced in the pseudo gap regime, indicating tunneling-like transport (23). Eventually, further away from the gap, in the regime with , fluctuations decay rapidly.

To discriminate between tunneling-like transport and Anderson localization, we consider the statistical distribution of the conductance. A quantitative description for transport fluctuations in quasi-one-dimensional systems is given by the Dorokhov-Mello-Pereyra-Kumar (DMPK) equation (26); (27); (28). DMPK predicts a crossover between the diffusive regime and the localized regime. While is Gaussian in the diffusive regime, it has a peculiar shape at the onset of the localized regime corresponding to with a marked discontinuity in the first derivative of the distribution and a sharp cut-off beyond  (28); (29). As decreases, this cut-off eventually leads to one-sided log-normal distributions for the conductance. Although the DMPK results were derived for the quasi-one-dimensional case, the characteristic shape of at the onset of the localized regime was exactly reproduced by numerical simulations in two-dimensional disordered systems in the SPS regime (30). We apply a similarity analysis of the numerical results in our HDPM system and the DMPK result (28) to identify regions where Anderson localization occurs. A map of the similarity S is shown in Fig. 4. Outside the gap and pseudo gap regions and light transport evolves from diffusion, where the conductance distributions are Gaussian, to Anderson localization, where the distribution is one-side log-normal. The behaviour of in the pseudo-gap region indicates remarkable differences in the statistical properties, which suggests that SPS is violated in this region.

In conclusion, our numerical results fully support the proposed transport phase diagram in Fig. 1. We expect that our findings are not restricted to SHU correlations, encoded by the parameter , but equally apply to other types of uniform spatial correlations (17). A similar transport phase diagram can be established in three dimensions. In this case a sharp phase boundary, known as the mobility edge and set by the Ioffe-Regel criterion  (21), is expected between photon diffusion and Anderson localization (dash-dotted line in Fig. 1). An important direction for future work are the implications of the phase diagram for electronic transport. It has been argued that SHU plays a role in the formation of electronic band gaps, for example in amorphous silicon (31). However, the influence of structural correlations and hyperuniformity on electron transport and localization in two or higher dimensions is far from being understood (32).

I Methods

Generation of stealthy hyperuniform point pattern. We employ a simulated annealing relaxation scheme to generate disordered SHU patterns with for as described in our previous work, ref.  (17). Patterns below the critical parameter , above which quasi-long-range order gradually appears (2), already show significant short-range order. The point patterns are decorated with dielectric cylinders as described in the text.

Band structure calculation. We calculate the normalized photonic density of states (NDOS) using the supercell method (6) implemented in the open source code MIT Photonic Bands (33). The supercell is repeated periodically and the band structure calculated by following the path in reciprocal space. The procedure is described in detail in our reference (17). Nearly identical strategies have also been applied by others (3).

Generalized scattering matrix (GSM) method. We use the improved generalized scattering matrix method. A summary of the method is provided in  (24). In the GSM the systems is discretized in slices in the propagation direction. For each slice the wave equation is solved and the scattering matrix calculated. By sequentially combining the corresponding scattering matrices, the total scattering matrix of the system up to a length is obtained.

Dorokhov-Mello-Pereyra-Kumar (DMPK) equation. The Dorokhov-Mello-Pereyra-Kumar (DMPK) equation is a Fokker-Plank equation describing the evolution of coherent wave transport statistics as a function of the ratio of the system length to the transport mean free path. It was shown to describe quantitatively the transport properties of disordered quasi-one-dimensional systems. There is evidence that the DMPK distributions retain the main properties of the conductance distributions in the metallic, critical and localized regime also in higher dimensions (29); (19). The DMPK distributions were obtained as described in ref. (34), see also supplementary material.

Similarity analysis. The similarity between the conductance distributions of correlated systems and the ideal DMPK distribution is characterized by a similarity function. We first create numerically a finite sample of DMPK conductance histograms with a certain bin size. We quantify these differences between the ideal DMPK sample and the results obtained by the GSM by calculating a squared distance distribution function . The squared distance between two distributions and is given by


In other contexts the squared distance is denoted or chi-squared. We do not make use of this notation to avoid confusion with the stealthiness parameter . For a set of different, but otherwise statistically equivalent, finite-sampling distributions we can define the distribution of squared distances considering all distances corresponding to all pairs with . To assess the similarity between and the reference distribution, we compute the area below both distributions. Similarity is defined as


By this definition is bounded in (see also supplementary material).


This research was supported by the Swiss National Science Foundation through the National Center of Competence in Research Bio-Inspired Materials and through Project No. 149867 and 169074 (L.F. and F.S.), Spanish Ministerio de Economía y Competitividad (MICINN) and European Regional Development Fund (ERDF) through Project FIS2015-69295-C3-3-P and the Basque Dep. de Educación through Project PI-2016-1-0041 (J.J.S.). M.E. acknowledges funding by Deutsche Forschungsgemeinschaft through the Cluster of Excellence Engineering of Advanced Materials and support from the Central Institute for Scientific Computing (ZISC) and the Interdisciplinary Center for Functional Particle Systems (IZ-FPS) at FAU Erlangen-Nuremberg.


Figure 1: Photonic transport phase diagram for hyperuniform disordered dielectric materials. A box of size is filled with high refractive index () cylinders, radius , at a filling fraction of  (17). In the central region (I) a wide photonic gap opens inhibiting the propagation of TM polarized electromagnetic waves (3); (17). In this region the density of states is exactly zero and evanescent light waves decay exponentially over distances shorter than the characteristic structural length scale . In the stealth region (II), on the left side of the band gap, the material is transparent. In the vicinity of the gap, for smaller values of , the density of states is suppressed but non-zero (III), and evanescent waves can tunnel between isolated states with a decay length . Sufficiently far from the gap diffuse light transport may cross over to strong Anderson localization. The transition from diffusive to localizated transport (dotted line) is system-size-dependent and we expect Anderson localization (IV) in the limit of sample sizes . The degree of hyperuniformity is denoted . is the optical frequency, the vacuum speed of light, and is the mean distance between scatterers. Eventually, at high (not shown), we enter the crystal regime (3); (17).
Figure 2: Ratio of the characteristic size of the system to (a) the mean distance between states and (b) the decay length (color maps) as a function of frequency and stealthiness for SHU systems. The decay length can only be extracted if . Gray areas in (b) thus correspond to systems where the condition has not been reached. The dashed line delineates the full band gap and the solid line indicates in Figs. 2-4.
Figure 3: Fluctuations of the logarithm of the conductance after TM-polarized electromagnetic waves passed through the system, . The inset shows the simulation setup of the SHU pattern with dielectric cylinders using periodic boundary conditions in transverse directions.
Figure 4: (a) Similarity map S for the conductance distribution at . High values indicate that the conductance distribution is similar to the one predicted by DMPK. The function S is bounded in . (b-g) Conductance distribution for two typical situations, in the SPS regime, , and in the tunneling regime inside the pseudo gap, . Numerical results are compared to predictions of DMPK (black lines).


  1. Batten, R. D., Stillinger, F. H. & Torquato, S. Classical disordered ground states: Super-ideal gases and stealth and equi-luminous materials. Journal of Applied Physics 104, 033504 (2008). eprint 0803.1442.
  2. Torquato, S., Zhang, G. & Stillinger, F. H. Ensemble theory for stealthy hyperuniform disordered ground states. Physical Review X 5, 21020 (2015).
  3. Florescu, M., Torquato, S. & Steinhardt, P. J. Designer disordered materials with large, complete photonic band gaps. Proceedings of the National Academy of Sciences 106, 20658–20663 (2009).
  4. Leseur, O., Pierrat, R. & Carminati, R. High-density hyperuniform materials can be transparent. Optica 3, 763 (2016).
  5. Chabanov, A. A., Stoytchev, M. & Genack, A. Z. Statistical signatures of photon localization. Nature 404, 850–853 (2000).
  6. Joannopoulos, J. D., Johnson, S. G., Winn, J. N. & Meade, R. D. Photonic Crystals: Molding the Flow of Light (Princeton University Press, 2008).
  7. Benedek, G. B. Theory of transparency of the eye. Applied Optics 10, 459 (1971).
  8. Rojas-Ochoa, L. F., Mendez-Alcaraz, J. M., Sáenz, J. J., Schurtenberger, P. & Scheffold, F. Photonic properties of strongly correlated colloidal liquids. Physical Review Letters 93, 73903 (2004).
  9. Reufer, M., Rojas-Ochoa, L. F., Eiden, S., Sáenz, J. J. & Scheffold, F. Transport of light in amorphous photonic materials. Applied Physics Letters 91, 171904 (2007).
  10. García, P. D., Sapienza, R., Blanco, Á. & López, C. Photonic glass: A novel random material for light. Advanced Materials 19, 2597–2602 (2007).
  11. Prum, R. O., Torres, R. H., Williamson, S. & Dyck, J. Coherent light scattering by blue feather barbs. Nature 396, 28–29 (1998).
  12. Ashcroft, N. W. & Lekner, J. Structure and resistivity of liquid metals. Physical Review 145, 83–90 (1966).
  13. Fraden, S. & Maret, G. Multiple light scattering from concentrated, interacting suspensions. Physical Review Letters 65, 512–515 (1990).
  14. Conley, G. M., Burresi, M., Pratesi, F., Vynck, K. & Wiersma, D. S. Light transport and localization in two-dimensional correlated disorder. Physical Review Letters 112, 143901 (2014).
  15. Wiersma, D. S., van Albada, M. P., van Tiggelen, B. A. & Lagendijk, A. Experimental evidence for recurrent multiple scattering events of light in disordered media. Physical Review Letters 74, 4193–4196 (1995).
  16. Rezvani Naraghi, R., Sukhov, S., Sáenz, J. J. & Dogariu, A. Near-field effects in mesoscopic light transport. Physical Review Letters 115, 203903 (2015).
  17. Froufe-Pérez, L. S. et al. Role of short-range order and hyperuniformity in the formation of band gaps in disordered photonic materials. Physical Review Letters 117, 053902 (2016).
  18. Abrahams, E., Anderson, P. W., Licciardello, D. C. & Ramakrishnan, T. V. Scaling theory of localization: Absence of quantum diffusion in two dimensions. Physical Review Letters 42, 673–676 (1979).
  19. Markoš, P. Numerical analysis of the Anderson localization. Acta Physica Slovaca 56, 561–685 (2006).
  20. Edagawa, K., Kanoko, S. & Notomi, M. Photonic amorphous diamond structure with a 3D photonic band gap. Physical Review Letters 100, 013901 (2008).
  21. John, S. Localization of light: Theory of photonic band gap materials. In Photonic Band Gap Materials, 563–665 (Springer Netherlands, Dordrecht, 1996).
  22. Deych, L. I., Lisyansky, A. A. & Altshuler, B. L. Single-parameter scaling in one-dimensional Anderson localization: Exact analytical solution. Physical Review B 64, 224202 (2001).
  23. Prior, J., Somoza, A. M. & Ortuño, M. Conductance fluctuations and single-parameter scaling in two-dimensional disordered systems. Physical Review B 72, 024206 (2005).
  24. A. Torres, J. & J. Sáenz, J. Improved generalized scattering matrix method: Conduction through ballistic nanowires. Journal of the Physical Society of Japan 73, 2182–2193 (2004).
  25. Beenakker, C. W. J. Random-matrix theory of quantum transport. Reviews of Modern Physics 69, 731–808 (1997).
  26. Dorokhov, O. N. Transmission coefficient and the localization length of an electron in N bound disordered chains. ZhETF Pisma Redaktsiiu 36, 259 (1982).
  27. Mello, P., Pereyra, P. & Kumar, N. Macroscopic approach to multichannel disordered conductors. Annals of Physics 181, 290–317 (1988).
  28. Froufe-Pérez, L. S., García-Mochales, P., Serena, P. A., Mello, P. A. & Sáenz, J. J. Conductance Distributions in Quasi-One-Dimensional Disordered Wires. Physical Review Letters 89, 246403 (2002).
  29. Muttalib, K. A. & Wölfle, P. ?One-sided? log-normal distribution of conductances for a disordered quantum wire. Physical Review Letters 83, 3013–3016 (1999).
  30. Somoza, A. M., Prior, J., Ortuño, M. & Lerner, I. V. Crossover from diffusive to strongly localized regime in two-dimensional systems. Physical Review B 80, 212201 (2009).
  31. Xie, R. et al. Hyperuniformity in amorphous silicon based on the measurement of the infinite-wavelength limit of the structure factor. Proceedings of the National Academy of Sciences 110, 13250–13254 (2013).
  32. Izrailev, F., Krokhin, A. & Makarov, N. Anomalous localization in low-dimensional systems with correlated disorder. Physics Reports 512, 125–254 (2012).
  33. Johnson, S. & Joannopoulos, J. Block-iterative frequency-domain methods for Maxwell’s equations in a planewave basis. Optics Express 8, 173–190 (2001).
  34. Mello, P., Yépez, M., Froufe-Pérez, L. & Sáenz, J. Statistical scattering of waves in disordered waveguides: Universal properties. Physica A: Statistical Mechanics and its Applications 372, 203–209 (2006).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minumum 40 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description