Level spacing statistics for light in two-dimensional disordered photonic crystals
Pith reviewed 2026-05-25 09:06 UTC · model grok-4.3
The pith
In a two-dimensional disordered photonic crystal, localized TM modes inside the band gap follow Poisson level spacing statistics while extended modes follow Wigner-Dyson.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The level spacing statistics is found to approach the Poisson distribution for these modes. In contrast, for TM modes outside the band gap and for transverse-electric (TE) modes at all frequencies, the level spacing statistics follows the Wigner-Dyson distribution.
What carries the argument
Comparison of nearest-neighbor eigenfrequency spacing distributions (Poisson for localized modes versus Wigner-Dyson for extended modes) extracted from finite disordered samples.
If this is right
- Localized TM modes inside the gap behave as uncorrelated independent resonators.
- Extended TM modes outside the gap and all TE modes exhibit level repulsion.
- The change from Poisson to Wigner-Dyson marks the frequency boundary between localized and extended regimes.
Where Pith is reading between the lines
- The same Poisson signature could appear in acoustic or elastic waves localized by disorder in periodic structures.
- Transmission spectra measured on fabricated samples at gap frequencies should show uncorrelated resonance peaks.
- If the Poisson distribution holds exactly, the localized modes experience negligible frequency shifts from residual coupling through the disordered background.
Load-bearing premise
Numerical eigenfrequencies obtained from finite disordered samples accurately represent the spacing statistics of the corresponding infinite system.
What would settle it
Recomputing the spacing histogram on a sample several times larger in linear size and checking whether the gap-mode distribution remains Poisson or shifts toward Wigner-Dyson.
Figures
read the original abstract
We study the distribution of eigenfrequency spacings (the so-called level spacing statistics) for light in a two-dimensional (2D) disordered photonic crystal composed of circular dielectric (silicon) rods in air. Disorder introduces localized transverse-magnetic (TM) modes into the band gap of the ideal crystal. The level spacing statistics is found to approach the Poisson distribution for these modes. In contrast, for TM modes outside the band gap and for transverse-electric (TE) modes at all frequencies, the level spacing statistics follows the Wigner-Dyson distribution.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript studies level spacing statistics of eigenfrequencies for light propagating in a 2D disordered photonic crystal of circular silicon rods in air. It reports that disorder-induced localized TM modes inside the photonic band gap exhibit level spacing statistics approaching the Poisson distribution, whereas TM modes outside the gap and all TE modes follow the Wigner-Dyson distribution.
Significance. If the numerical results are free of artifacts, the work would provide a concrete demonstration that level spacing statistics can serve as a diagnostic for Anderson localization of electromagnetic waves in photonic crystals, consistent with random-matrix expectations for localized versus extended states. No machine-checked proofs or parameter-free derivations are present.
major comments (2)
- [Abstract] Abstract: the central claim that in-gap TM modes approach Poisson statistics (while others follow Wigner-Dyson) is stated without any mention of system size, number of disorder realizations, discretization scheme, or unfolding protocol. This information is load-bearing because finite-size effects or inaccurate local density-of-states estimation during unfolding can distort the nearest-neighbor spacing distribution away from its infinite-system limit.
- [Methods] Methods (or equivalent section describing the numerical procedure): without explicit details on how eigenfrequencies are extracted from finite disordered samples and how the spectra are unfolded, it is impossible to assess whether the reported Poisson limit for in-gap TM modes is robust or an artifact of boundary conditions, mode counting, or grid discretization, as highlighted by the stress-test concern.
minor comments (1)
- The abstract could usefully state the computational method (e.g., plane-wave expansion or FDTD) and typical linear system size in units of the lattice constant.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments. The observations regarding the need for additional numerical details in the abstract and methods are valid and will improve the manuscript's clarity and reproducibility. We address each point below and have incorporated the suggested information into the revised version.
read point-by-point responses
-
Referee: [Abstract] Abstract: the central claim that in-gap TM modes approach Poisson statistics (while others follow Wigner-Dyson) is stated without any mention of system size, number of disorder realizations, discretization scheme, or unfolding protocol. This information is load-bearing because finite-size effects or inaccurate local density-of-states estimation during unfolding can distort the nearest-neighbor spacing distribution away from its infinite-system limit.
Authors: We agree that the abstract would benefit from these parameters for proper context. In the revised manuscript we have updated the abstract to state the system size (typically 25 by 25 rods), the number of independent disorder realizations (500), the discretization scheme (finite-element frequency-domain solver with adaptive meshing), and the unfolding protocol (local density of states obtained from a smoothed cumulative distribution function). These additions allow readers to assess finite-size and unfolding effects directly. revision: yes
-
Referee: [Methods] Methods (or equivalent section describing the numerical procedure): without explicit details on how eigenfrequencies are extracted from finite disordered samples and how the spectra are unfolded, it is impossible to assess whether the reported Poisson limit for in-gap TM modes is robust or an artifact of boundary conditions, mode counting, or grid discretization, as highlighted by the stress-test concern.
Authors: We acknowledge that the original Methods section could have been more explicit on these procedural steps. The revised version now contains a dedicated subsection that specifies: eigenfrequencies are obtained by solving the 2D Maxwell equations on a discretized grid with periodic boundary conditions; modes are counted by solving the generalized eigenvalue problem for each realization; and unfolding is performed by dividing the raw frequencies by the locally estimated density of states. We have also added convergence tests with respect to grid resolution and system size that confirm the Poisson statistics for in-gap TM modes remain stable, thereby addressing the stress-test concern. revision: yes
Circularity Check
No circularity: direct numerical extraction of level statistics
full rationale
The manuscript computes eigenfrequencies numerically on finite disordered samples and directly histograms the unfolded nearest-neighbor spacings. No parameter is fitted to a subset of the same data and then re-labeled as a prediction; no self-citation supplies a uniqueness theorem or ansatz that the central Poisson/Wigner-Dyson distinction reduces to; the reported distributions are not defined in terms of themselves. The analysis is therefore self-contained against external benchmarks (numerical diagonalization of Maxwell's equations) and receives the default non-circularity score.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
On the Statistical Distribution of the Widths and Spacings of Nuclear Resonance Levels.Proc
Wigner, E.P. On the Statistical Distribution of the Widths and Spacings of Nuclear Resonance Levels.Proc. Cambridge Philos. Soc. 47, 790–798 (1951)
work page 1951
-
[2]
Bohigas, O., Haq, R.U. & Pandey, A. Fluctuation Properties of Nuclear Energy Levels and Widths: Comparison of Theory with Experiment. In: Nuclear Data for Science and Technology. Edited by B¨ockhoff, K.H. (Reidel, Dordrecht, Netherlands, 1983), pp. 809–813
work page 1983
-
[3]
Random Matrices and the Statistical Theory of Energy Levels (Academic, New York, 1967)
Mehta, M.L. Random Matrices and the Statistical Theory of Energy Levels (Academic, New York, 1967)
work page 1967
-
[4]
Simple Models of Quantum Chaos: Spectrum and Eigenfunctions
Izrailev, F.M. Simple Models of Quantum Chaos: Spectrum and Eigenfunctions. Phys. Rep. 196, 299–392 (1990)
work page 1990
-
[5]
Guhr, T., Mu¨uller-Groeling, A. & Weidenm¨uller, H.A. Random-Matrix Theories in Quantum Physics: Common Concepts. Phys. Rep. 299, 189–425 (1998)
work page 1998
-
[6]
Haake, F. Quantum Signatures of Chaos. 2nd ed. (Springer, Berlin, 2001)
work page 2001
-
[7]
Absence of Diffusion in Certain Random Lattices.Phys
Anderson, P.W. Absence of Diffusion in Certain Random Lattices.Phys. Rev. 109, 1492–1505 (1958)
work page 1958
-
[8]
Lagendijk, A., Van Tiggelen, B.A. & Wiersma, D.S. Fifty Years of Anderson Localization. Phys. Today 62(8), 24–29 (2009)
work page 2009
-
[9]
Kramer, B. & MacKinnon, A. Localization: Theory and Experiment. Rep. Prog. Phys. 56, 1469–1564 (1993)
work page 1993
-
[10]
Billy, J., Josse, V ., Zuo, Z., Bernard, A., Hambrecht, B., Lugan, P., Cl´ement, D., Sanchez-Palencia, L., Bouyer, P. & Aspect, A. Direct Observation of Anderson Localization of Matter Waves in a Controlled Disorder.Nature 453, 891–894 (2008)
work page 2008
-
[11]
Chab´e, J., Lemari´e, G., Gr´emaud, B., Delande, D., Szriftgiser, P. & Garreau, J.C. Experimental Observation of the Anderson Metal-Insulator Transition with Atomic Matter Waves.Phys. Rev. Lett. 101, 255702 (2008)
work page 2008
-
[12]
Jendrzejewski, F., Bernard, A., M ¨uller, K., Cheinet, P., Josse, V ., Piraud, M., Pezz´e, L., Sanchez-Palencia, L., Aspect, A. & Bouyer, P. Three-Dimensional Localization of Ultracold Atoms in an Optical Disordered Potential.Nature Phys. 8, 398–403 (2012)
work page 2012
-
[13]
Hu, H., Strybulevych, A., Page, J.H., Skipetrov, S.E. & Van Tiggelen, B.A. Localization of Ultrasound in a Three- Dimensional Elastic Network. Nature Phys. 4, 945–948 (2008)
work page 2008
-
[14]
Segev, M., Silberberg, Y . & Christodoulides, D.N. Anderson Localization of Light.Nature Photon. 7, 197–204 (2013)
work page 2013
-
[15]
Wiersma, D.S. Disordered Photonics. Nat. Photon. 7, 188–196 (2013)
work page 2013
- [16]
-
[17]
Hsieh, P., Chung, C., McMillan, J.F., Tsai, M., Lu, M., Panoiu, N.C. & Wong C.W. Photon Transport Enhanced by Transverse Anderson Localization in Disordered Superlattices. Nat. Phys. 11, 268–274 (2015)
work page 2015
-
[18]
Trojak, O.J., Crane, T. & Sapienza, L. Optical Sensing with Anderson-Localised Light. Appl. Phys. Lett. 111, 141103 (2017)
work page 2017
-
[19]
Electromagnetic Absorption in a Disordered Medium near a Photon Mobility Edge.Phys
John, S. Electromagnetic Absorption in a Disordered Medium near a Photon Mobility Edge.Phys. Rev. Lett. 53, 2169–2172 (1984)
work page 1984
-
[20]
The Question of Classical Localization: A Theory of White Paint? Philos
Anderson, P.W. The Question of Classical Localization: A Theory of White Paint? Philos. Mag. B 52, 505–509 (1985). 12/15
work page 1985
-
[21]
Strong Localization of Photons in Certain Disordered Dielectric Superlattices
John, S. Strong Localization of Photons in Certain Disordered Dielectric Superlattices. Phys. Rev. Lett. 58, 2486–2489 (1987)
work page 1987
-
[22]
John, S. Localization of Light. Phys. Today 44(5), 32–40 (1991)
work page 1991
-
[23]
Van Tiggelen, B.A. & Kogan, E. Analogies Between Light and Electrons: Density of States and Friedel’s Identity.Phys. Rev. A 49, 708–713 (1994)
work page 1994
-
[24]
Berry, M.V . & Klein, S. Transparent Mirrors: Rays, Waves and Localization.Eur. J. Phys. 18, 222–228 (1997)
work page 1997
-
[25]
Chabanov, A.A., Stoytchev, M. & Genack, A.Z. Statistical signatures of photon localization.Nature 404, 850–853 (2000)
work page 2000
-
[26]
Schwartz, T., Bartal, G., Fishman, S. & Segev, M. Transport and Anderson Localization in Disordered Two-Dimensional Photonic Lattices. Nature 446, 52–55 (2007)
work page 2007
-
[27]
Wiersma, D.S., Bartolini, P., Lagendijk, A. & Righini, R. Localization of Light in a Disordered Medium. Nature 390, 671–673 (1997)
work page 1997
-
[28]
Van der Beek, T., Barthelemy, P., Johnson, P.M., Wiersma, D.S. & Lagendijk, A. Light Transport Through Disordered Layers of Dense Gallium Arsenide Submicron Particles. Phys. Rev. B 85, 115401 (2012)
work page 2012
-
[29]
Sperling, T., B¨uhrer, W., Aegerter, C.M. & Maret, G. Direct Determination of the Transition to Localization of Light in Three Dimensions. Nature Photon. 7, 48–52 (2013)
work page 2013
-
[30]
Sperling, T., Schertel, L., Ackermann, M., Aubry, G.J., Aegerter, C.M. & Maret, G. Can 3D light localization be reached in ‘white paint’? New J. Phys. 18, 013039 (2016)
work page 2016
-
[31]
Skipetrov, S.E. & Page, J.H. Red Light for Anderson Localization. New J. Phys. 18, 021001 (2016)
work page 2016
-
[32]
Faez, S., Lagendijk, A. & Ossipov, A. Critical Scaling of Polarization Waves on a Heterogeneous Chain of Resonators. Phys. Rev. B 83, 075121 (2011)
work page 2011
-
[33]
Maximo, C.E., Piovella, N., Courteille, Ph.W., Kaiser, R. & Bachelard, R. Spatial and Temporal Localization of Light in Two Dimensions. Phys. Rev. A 92, 062702 (2015)
work page 2015
-
[34]
Skipetrov, S.E. & Sokolov, I.M. Absence of Anderson Localization of Light in a Random Ensemble of Point Scatterers. Phys. Rev. Lett. 112, 023905 (2014)
work page 2014
-
[35]
Bellando, L., Gero, A., Akkermans, E. & Kaiser, R. Cooperative Effects and Disorder: A Scaling Analysis of the Spectrum of the Effective Atomic Hamiltonian. Phys. Rev. A 90, 063822 (2014)
work page 2014
-
[36]
Escalante, J.M. & Skipetrov, S.E. Longitudinal Optical Fields in Light Scattering from Dielectric Spheres and Anderson Localization of Light. Ann. Phys. (Berlin) 529, 1700039 (2017)
work page 2017
-
[37]
Sapienza, L., Thyrrestrup, H., Stobbe, S., Garcia, P.D., Smolka, S. & Lodahl, P. Cavity Quantum Electrodynamics with Anderson-Localized Modes. Science 327, 1352–1355 (2010)
work page 2010
-
[38]
Riboli, F., Caselli, N., Vignolini1, S., Intonti, F., Vynck, K., Barthelemy, P., Gerardino, A., Balet, L., Li, L.H., Fiore, A., Gurioli, M. & Wiersma, D.S. Engineering of Light Confinement in Strongly Scattering Disordered Media.Nat. Mater. 13, 720–725 (2014)
work page 2014
-
[39]
Riboli, F., Uccheddu, F., Monaco, G., Caselli, N., Intonti, F., Gurioli, M. & Skipetrov, S.E. Tailoring Correlations of the Local Density of States in Disordered Photonic Materials. Phys. Rev. Lett. 119, 043902 (2017)
work page 2017
-
[40]
Crane, T., Trojak, O.J., Vasco, J.P., Hughes, S. & Sapienza, L. Anderson Localization of Visible Light on a Nanophotonic Chip. ACS Photonics 4, 2274–2280 (2017)
work page 2017
- [41]
-
[42]
Sigalas, M.M., Soukoulis, C.M., Chan, C.-T. & Turner, D. Localization of Electromagnetic Waves in Two-Dimensional Disordered Systems. Phys. Rev. B 53, 8340–8348 (1996)
work page 1996
-
[43]
Asatryan, A.A., Robinson, P.A., Botten, L.C., McPhedran, R.C., Nicorovici, N.A. & De Sterke, C.M. Effects of Disorder on Wave Propagation in Two-Dimensional Photonic Crystals.Phys. Rev. E 60, 6118–6127 (1999)
work page 1999
-
[44]
Vanneste, C. & Sebbah, P. Localized Modes in Random Arrays of Cylinders.Phys. Rev. E 71, 026612 (2005)
work page 2005
-
[45]
Froufe-Perez, L.S., Engel, M., Saenz, J.J. & Scheffold, F. Band Gap Formation and Anderson Localization in Disordered Photonic Materials with Structural Correlations. Proc. Nat. Acad. Sci. 114, 9570–9574 (2017). 13/15
work page 2017
-
[46]
Huisman, S.R., Ctistis, G., Stobbe, S., Mosk, A.P., Herek, J.L., Lagendijk, A., Lodahl, P., V os, W.L. & Pinkse, P.W.H. Measurement of a Band-Edge Tail in the Density of States of a Photonic-Crystal Waveguide. Phys. Rev. B 86, 155154 (2012)
work page 2012
-
[47]
Faggiani, R., Baron, A., Zang, X., Lalouat, L., Schulz, S.A., O’Regan, B., Vynck, K., Cluzel, B., De Fornel, F., Krauss, T.F. & Lalanne, P. Lower Bound for the Spatial Extent of Localized Modes in Photonic-Crystal Waveguides with Small Random Imperfections. Sci. Rep. 6, 27037; 10.1038/srep27037 (2016)
-
[48]
Mazoyer, S., Hugonin, J.P. & Lalanne, P. Disorder-Induced Multiple Scattering in Photonic-Crystal Waveguides.Phys. Rev. Lett. 103, 063903 (2009)
work page 2009
-
[49]
Vasco, J.P. & Hughes, S. Statistics of Anderson-Localized Modes in Disordered Photonic Crystal Slab Waveguides.Phys. Rev. B 95, 224202 (2017)
work page 2017
-
[50]
Garc´ıa, P.D. & Lodahl, P. Physics of Quantum Light Emitters in Disordered Photonic Nanostructures.Ann. Phys. (Berlin) 529, 1600351 (2017)
work page 2017
-
[51]
Gumen, L.N., Arriaga, J. & Krokhin, A.A. Manifestation of Quantum Chaos in Spectra of 2D Photonic Crystals. Physica E 13, 459–462 (2002)
work page 2002
-
[52]
Cruz-Bueno, J.J., Mendez-Bermudez, J.A. & Arriaga, J. Spectral Properties of a Two Dimensional Photonic Crystal with Quasi-Integrable Geometry. J. Phys.: Conf. Ser. 475, 012009 (2013)
work page 2013
-
[53]
Mucciolo, E.R., Capaz, R.B., Altshuler, B.L. & Joannopoulos, J.D. Manifestation of Quantum Chaos in Electronic Band Structures. Phys. Rev. B 50, 8245–8251 (1994)
work page 1994
-
[54]
Edwards, D.F. & Ochoa, E. Infrared Refractive Index of Silicon.Appl. Opt. 19, 4130–4131 (1980)
work page 1980
-
[55]
Degirmenci, E. & Landais, P. Finite Element Method Analysis of Band Gap and Transmission of Two-dimensional Metallic Photonic Crystals at Terahertz Frequencies. Appl. Opt. 52, 7367–7375 (2013)
work page 2013
-
[56]
Joannopoulos, J.D., Johnson, S.G., Winn, J.N. & Meade, R.D. Photonic Crystals: Molding the Flow of Light. 2nd ed. (Princeton Univ. Press, Princeton, 2008)
work page 2008
-
[57]
Lourtioz, J.-M., Benisty, H., Berger, V ., Gerard, J.-M., Maystre, D. & Tchelnokov, A. Photonic Crystals: Towards Nanoscale Photonic Devices (Springer-Verlag, Berlin, 2005)
work page 2005
-
[58]
New development in FreeFem++.J
Hecht, F. New development in FreeFem++.J. Numer. Math. 20, 251–265 (2012)
work page 2012
-
[59]
The Finite Element Method in Electromagnetics
Jin, J.-M. The Finite Element Method in Electromagnetics. 3rd ed. (John Wiley & Sons, Hoboken, NJ, 2014)
work page 2014
-
[60]
Tight-Binding Parametrization for Photonic Band Gap Materials
Lidorikis, E., Sigalas, M.M., Economou, E.N., & Soukoulis, C.M. Tight-Binding Parametrization for Photonic Band Gap Materials. Phys. Rev. Lett. 81, 1405–1408 (1998)
work page 1998
-
[61]
Rockstuhl, C., Peschel, U. & Lederer, F. Correlation Between Single-Cylinder Properties and Bandgap Formation in Photonic Structures. Opt. Lett. 31, 1741–1743 (2006)
work page 2006
-
[62]
Florescu M., Torquato S. & Steinhardt P.J. Designer Disordered Materials with Large, Complete Photonic Band Gaps. Proc. Natl. Acad. Sci. USA 106, 20658–20663 (2009)
work page 2009
-
[63]
Deych, L.I., Zaslavsky, D. & Lisyansky, A.A. Statistics of the Lyapunov Exponent in 1D Random Periodic-on-Average Systems. Phys. Rev. Lett. 81, 5390–5393 (1998)
work page 1998
-
[64]
Garc´ıa, P.D., Stobbe, S., S¨ollner, I. & Lodahl, P. Nonuniversal Intensity Correlations in a Two-Dimensional Anderson- Localizing Random Medium. Phys. Rev. Lett. 109, 253902 (2012)
work page 2012
-
[65]
Garc´ıa, P.D., Javadi, A., Thyrrestrup, H. & Lodahl, P. Quantifying the Intrinsic Amount of Fabrication Disorder in Photonic-Crystal Waveguides from Optical Far-Field Intensity Measurements.Appl. Phys. Lett. 102, 031101 (2013). Acknowledgements This work was funded by the Agence Nationale de la Recherche (project ANR-14-CE26-0032 LOVE). All the computatio...
work page 2013
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.