## Abstract

Hyperuniform structures possess the ability to confine and drive light, although their fabrication is extremely challenging. Here we demonstrate that speckle patterns obtained by a superposition of randomly arranged sources of Bessel beams can be used to generate hyperunifrom scalar fields. By exploiting laser light tailored with a spatial filter, we experimentally produce (without requiring any computational power) a speckle pattern possessing maxima at locations corresponding to a hyperuniform distribution. By properly filtering out intensity fluctuation from the same speckle pattern, it is possible to retrieve an intensity profile satisfying the hyperuniformity requirements. Our findings are supported by extensive numerical simulations.

© 2018 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

## 1. Introduction

Disorder [1] is generally considered detrimental for imaging or optical information transport, nevertheless, complex systems can disclose extraordinary physical properties, allowing the confinement of light [2], and enabling high resolution imaging techniques [3].

The task of driving light at the nanoscale is usually fulfilled by homogeneous waveguides or photonic crystals supporting a band gap [4]: completely ordered structures which allow to tailor the light flow. On the other hand, recent studies on complex materials and photonic glasses [5, 6] demonstrated that complete gaps can be found also in case of designed disorder [7, 8], which is today inspiring new platforms to overcome the present limitations encountered in photonic applications [9] and enabling extraordinary optical phenomena [10–13].

Light propagating in disordered media generates speckle patterns that are also complex light features, they are originated by random scattering through or by a rough dielectric slab. These light structures can be tailored to produce semi-organized light patterns [14–16] able to constitute random potential-energy landscapes [17, 18] in use for particles [19] or ultracold atoms [20] manipulation and with striking properties if adopted for structured illuminations [21, 22]. The Amorphous Speckle Patterns (ASP) are a kind of light structures with a short-range correlated light arrangement constructed from random components residing on a fixed ring in momentum space, they correspond to a superposition of Bessel modes that lends the pattern the self-healing nature [14, 22–25], an advantageous property for enhancing optical imaging in thick tissues [22, 24]. Besides that, those specular features are labeled as “amorphous” for exhibiting a non-crystalline order, a property that has made them usable for the formation of amorphous photonic lattices in photorefractive crystals [23, 25] or photolithography [26]. Therefore, ASP can be exploited to fabricate static structures with disordered refractive index microstructure, as required to observe Anderson localization [2, 23, 25, 27] and are recurring features in non-linear phenomena like rogue waves [28] and random lasing [29, 30]. For these reasons ASPs are attracting large interest in the field of Optics and Photonics.

On the other hand, Anderson localization is often associated with hyperuniform dielectric structures and their direct correlation is currently examined [31, 32]. Hyperuniform systems [33] are exotic states of matter fabricated exploiting designed disorder laying between a crystal and a liquid: they behave similarly to a perfect crystal since they suppress large scale density fluctuation, but are also statistically isotropic with no Braggs peaks as in liquid and glasses [33–35]. They have been observed in different physical systems ranging from disordered ground state [36], to jammed particle packing [37], in processes of ultracold atoms [38] and in driven non-equilibrium systems [39]. A statistically homogeneous hyperuniform points configuration in *d*-dimensions is one in which the number variance ${\sigma}_{N}^{2}(R)$ of points within a spherical observation window of radius $R$ grows more slowly than ${R}^{d}$

The hyperuniform structures are currently subject of study in the regulation of light confinement associated with their photonic band gaps [31,32, 42] and are exploited for photonic applications at the micro-scale [42–46]. Nevertheless, the creation of hyperuniform heterogeneous material [40] requires very strict fabrication protocols resulting very challenging and costly, even though modern lithographic methods are able to produce super defined hyperuniform structures [47–49].

Hyperuniformity has been also generalized to heterogeneous media [35, 40] and random scalar fields [41]. In the latter case the hyperuniformity concept can be generalized by replacing the $S\left(k\right)$ with the spectral density $\tilde{\psi}\left(\text{k}\right)$ [41]. Given a continuous speckle intensity distribution C (a random scalar field), it is possible to calculate the auto-correlation function, $\chi \left(\text{r}\right)$ [50], as the inverse Fourier transform of the speckle energy-spectrum:

With the study of these properties, we demonstrate that it is possible to obtain a continuous hyperuniform field distribution from a standard laser speckle by applying an opportune intensity filtering in Fourier domain.

The hyperuniform condition from Eq. (2) can be simply extended to scalar fields that respect the following Eq.:

Currently, the challenge is on the identification of analogue real system capable to replicate hyperuniform continuous field distribution theorized with the recent works [33, 35, 41].In this paper we analyze ASP and investigate on their connection with hyperuniform continuous fields supporting our study with rigorous numerical simulations.

We demonstrate that the hyperuniform distribution of the ASP intensity peaks satisfy the conditions in Eqs. (1) and (6). To do so we examine normalized auto-correlation function, $\tilde{\chi}\left(\text{r}\right)$, and the spectral density $\tilde{\psi}\left(\text{k}\right)$ of tailored light fields distributions. Moreover, we demonstrate that it is possible to obtain a continuous variable hyper uniformity (by applying an opportune intensity filtering in Fourier domain) from a standard laser speckle.

## 2. Looking for hyperuniform continuous fields

In order to identify an optical equivalent that best replicate a continuous hyperuniform field, we simulate tailored distribution of non-random light fields and we examine synthetic speckle patterns directly designed from a given hyperuniform distribution. This preliminary study serves to extrapolate the parameters that characterize a hyperuniform speckle pattern and we intend to use those as starting point to build a more realistic hyperuniform continuous field distribution in our experiments.

The easiest way to simulate a hyperuniform continuous field is by adopting dephasors/modes with Gaussian spatial profile [50] placed in a [X,Y] grid according to the coordinates of a given hyperuniform matrix (the matrix has been provided by Authors of Ref [36].). A dephasor is a complex field with given amplitude and phase that, when is superimposed with other fields on a observation plane, generates an interference pattern. We analyse the speckle patterns generated from the superposition of 800 (number of points in the given hyperuniform matrix) Gauss modes (${\text{f}}_{\text{i}=1,\dots ,800}$) distributed in the grid according to a given $\text{\chi}=0.4$ stealthy hyperuniform map. The dimensionless parameter $\text{\chi}$ represents the ratio between the number of constrained degrees of freedom and the total number of degrees of freedom that have been set to build the hyperuniform map [33,36]. In practice, we replace the points of a hyperuniform distribution with interfering Gaussian dephasors/modes with the aim to form a continuous filed distribution. Taking in account that the average distance between neighbouring modes in the hyperuniform map is $\tilde{d}=368\text{pixels}$ (calculated as the distance between the peaks in $\tilde{\chi}\left(\text{r}\right)$ from the hyperuniform map), we compared the speckle patterns resulting from superposition of modes with $\frac{\text{\sigma}}{\tilde{d}}=0.054\text{}\left[\text{Fig}.\text{}1\left(\text{a}\right)\right],0.24\text{}\left[\text{Fig}.\text{}1\left(\text{b}\right)\right],0.43\text{}\left[\text{Fig}.\text{}1\left(\text{c}\right)\right],0.86\text{}\left[\text{Fig}.\text{}1\left(\text{d}\right)\right]$, with $\text{\sigma}$ the mode width (see Methods).

We first study the variation in $\tilde{\chi}\left(\text{r}\right)$ [line *II* of Fig. 1] and $\tilde{\psi}\left(\text{k}\right)$ [line *III* of Fig. 1] for the cases of $\text{\sigma}<\tilde{d}$ corresponding to Figs. 1(a), 1(b) and 1(c). Those patterns are secured (by the map) to maintain the hyperuniform geometry. Patterns (a), (b) and (c) from Fig. 1 are indeed three continuous fields exhibiting hyperuniform distributions as proved with the observation of their $\tilde{\psi}\left(\text{k}\right)$ reported in the line *III* of Figs. 1(a), 1(b) and 1(c). For the three cases, we have that $\tilde{\psi}\left(\text{k}\right)\text{}\to 0$ when $\text{k}\to 0$, meaning that the hyperuniform condition of Eq. (6) is respected.

When Gaussian sources are much smaller than their average distance, interference is negligible, on the contrary when they start consistently overlapping the interference between the fields becomes relevant and the hyperuniform structure is lost, this is the case of $\text{\sigma}\ge \tilde{d}\text{}$ in Fig. 1(d). Here, the hyperuniform from the map information is completely lost, i.e. the final speckle assumes the typical random distribution as demonstrated with the analysis of $\tilde{\chi}\left(\text{r}\right)$ and $\tilde{\psi}\left(\text{k}\right)$ shown in lines *II* and *III* of Fig. 1(d).

Therefore, the condition $\text{\sigma}\ll \tilde{d}$ must be respected to maintain hyperuniformity from points to continuous field distributions. Although this first intelligible analysis might result trivial, it discloses a very important piece of information: we have the evidence that when the $\text{\sigma}$ is larger (in the regime $\text{\sigma}\ll \tilde{d}$), the $\tilde{\chi}\left(\text{r}\right)$ tends to a Bessel function, suggesting that hyperuniformity for continuous fields requires a Bessel shaped auto-correlation function of the speckle distribution. This is shown by the statistical analysis illustrated in lines *II* and *III* of Fig. 1(c): the auto-correlations [blue curves in Figs. 1(*II*b) and (*II*c)] are directly compared with a Bessel function [dashed yellow curves in Figs. 1(*II*b) and (*II*c)] and spectral density takes the typical annular shape [Figs. 1(*III*b) and (*III*c)]. This is a statement that finds also theoretical support in recent studies demonstrating that continuous hyperunformity is well described with Bessel-like auto-correlation function [41]. From this fundamental concept we base our strategy to build hyperuniform continuous fields in our simulations and experiments.

#### 2.1 Superposition of Bessel modes

Standard speckle patterns result from interference of scrambled light fields, they represent a random distribution of maxima and minima in intensity. Their Fourier analysis is very well established [50]: correlation function with a sharp peak at short radius (around $r=0$) and decaying spectral density for larger wavenumbers (see Methods). Differently, ASPs exhibit Bessel-shaped auto-correlation function and a discontinuous spectral density: very bright around $k=0$ and null at larger wavenumbers unless those residing in a very tight interval around $k={k}^{\text{*}}\text{with}{k}^{\text{*}}\gg 0$ [14]. In this section we report numerical and experimental results demonstrating that the statistic of a superposition of Bessel beams is independent on their particular arrangement differently from what has resulted with the previous case illustrated in Fig. 1.

In our experiments, in order to generate ASPs, we apply a spatial frequency selection to the light transmitted through a scattering barrier, a filtering operation performed by an annular aperture (see Methods). Exploiting the annular filter we generate a superposition of plane waves with fixed *k* and random phases, which is known to give rise to an ASP.

On the other hand, in order to obtain a numerical comparison with our experiments, we simulate synthetic ASPs with superposition of modes/dephasors with Bessel spatial profile as described below:

*s*is the scaling factor that controls the diameter of the Bessel mode in a [X,Y] = 5000x5000 pixels grid, ${\text{R}}_{\text{i}}(\text{x}-{\text{x}}_{\text{i}},\text{y}-{\text{y}}_{\text{i}})$, ${\text{\phi}}_{\text{i}}\text{}$ represent respectively the position of the i-th mode in pixel unit and the phase associated. In the following we will always consider pixel as unit for the length, so that

*s*remains a dimensionless factor. Then, the resulting pattern is given by the following formula:

#### 2.2 Analysis of the synthetic speckle grains distribution: a comparison between random speckle pattern and amorphous speckle pattern

Here we describe the process of retrieving a hyperuniform pattern starting from the position of the intensity maxima from an amorphous pattern as the one shown in Fig. 2(c). Using a custom data analysis software we are able to find the local maxima of our speckle pattern and to obtain the position and the number of the speckle grains in each pattern. Similar analysis has been already adopted in previous works [52]. In such a way we are able to obtain a points’ map that provides the coordinate/position of the grains from the speckle patterns and then analyse the properties of their distribution. We extract the speckle grains coordinates from the two different conditions:

The software pinpoints each intensity peak (speckle grain) with yellow crosses as shown in Figs. 3(a) and 3(b) and extract their coordinates. Then we produce a synthetic frame replacing speckle intensity with a sum of Gaussians with equal intensity centred at the coordinates of the relative maxima. The synthetic frames are shown in Figs. 3(c) and 3(d). If two maxima are very close, then intensity of the two correspondent Gaussians is summed thus producing super-intense pixels in the synthetic frame. Note that these super-intense pixels are appearing only in Fig. 3(c) (standard speckles) highlighting that maxima are placed so close to overlap. On the other hand, in Fig. 3(d) (ASP) all maxima are sparse maintaining a sort of order. The structure factors $S(k)$ of the speckles of Figs. 3(a) and 3(b), are reported in Fig. 3(e) and 3(f). While for the case of conventional speckle pattern the $S(k)\ne 0$ for $k\to 0$ [Fig. 3(e)], in the case of ASP [Fig. 3(f)] the curve $S(k)$ tends to 0 when $k\to 0$, proving that the ASPs exhibit hyperuniform information intrinsically contained in the speckle grain distribution.

In other words, the sufficient condition to obtain intensity maxima with a hyperuniform distribution is to use Bessel beams as sources, therefore they will be the first candidates to produce *hyperuniform continuous fields*. Gaussian sources provided the same result but only by binding their locations accordingly with a hyperuniform map.

## 3. The generation of the hyperuniform continuous fields from experimental speckle patterns

From the previous section results that ASP maxima are hyperuniform, whereas the original pattern does not. This is due to strong intensity fluctuations (that depend upon local energy deposition) between different regions in the pattern. Normalizing the ASP with an appropriate low-pass we neglect the local average intensity fluctuations, thus we can examine only the speckle grains distribution. The procedure is described in the following section where we demonstrate that it is possible to obtain the hyperuniform continuous field starting from experimental amorphous speckle and with a non-perturbative data treatment. We define *A* as the 2-dimensional matrix with ${\text{n}}_{1}{\text{xn}}_{2}$ elements representing the ASP [Fig. 4(a)] obtained with our experimental apparatus. We generate a two-dimensional filter $h$ of the Gaussian type: a rotationally symmetric Gaussian low-pass filter of size ${h}_{size}$ with positive standard deviation ($\text{\sigma})$. We used ${h}_{size}$ $[{\text{n}}_{1}{\text{n}}_{2}]$ with same size as the speckle pattern matrix and $\text{\sigma}=20\text{pixels}$. We created the Gaussian filters using the following Eqs.:

*A*(the original speckle pattern) adopting the filter h in such a way to generate the envelope

*B*as follows:

*A*[Fig. 4(a)] dividing it with its envelope

*B*and we obtain the intensity corrected speckle pattern

*C*[Fig. 4(b)]:

*A*[panel (a) in Fig. 4] and

*C*[panel (b) in Fig. 4] we notice that they differ in their intensity distribution only as shown in the insets (c) and (d) of Fig. 4, zoom-in from patterns (a) and (b). In

*C*we observe intensity fluctuations (extraordinary peaks) as shown with the white circles in the Figs. 4(c) and 4(d). We retrieved the speckle grains positions in

*A*and

*C*and we calculated their cross correlation between the two maps obtaining the correlation coefficient 0.88 indicating that speckle grains position are not modified by the filtering process as expected. In practice, adopting the filtering process we eliminated the intensity fluctuations from

*A*and we obtain the “homogeneous” pattern shown in Fig. 4(b), the intensity corrected amorphous speckle pattern. If we analyse the statistics of the corrected pattern

*C*we have that $\tilde{\chi}\left(\text{k}\right)\to 0$ when $\text{k}\to 0$ as shown in Fig. 4(e), i.e. it is hyperuniform.

#### 3.1 Local field variance

In order to further confirm the hyperuniform distribution of ASPs we calculate the local scalar field fluctuations ${\text{R}}^{2}{\text{\sigma}}_{\text{F}}^{2}\left(\text{R}\right)$ of speckles within a circular observation window of radius $\text{R}$ [41, 42]. In this case, ${\text{\sigma}}_{\text{F}}^{}\left(\text{R}\right)$ is the natural extension of ${\text{\sigma}}_{\text{N}}$ for the analysis of continuous fields: it takes into account the total intensity variance within an observation window of size R as described in the work of Prof. Torquato [33, 35, 41].

The speckle patterns analysed are rasterized into a window with linear size $\text{W}$, where $\text{W}=400\text{pixels}$ according to the original image linear size.

Thereafter, we define a sub-window with variable radius R ranging from 2 pixels up to the maximum window radius ${\text{R}}_{\text{max}}$, where ${\text{R}}_{\text{max}}=\frac{\text{W}}{4}=100\text{pixels}$. We count the total intensity found within the region of interest (the sub-window) and we store this data while scanning with the sub-window the whole mask. The scan is performed in step of 10 pixels in both directions x and y in the mask. Boundary conditions are taken into account as the mask would have been replicated in each direction, in this way if the sub-window exits the mask on one side it emerges from the opposite. This allows to estimate the average intensity (the volume) residing in each sub-window size (from 2 pixels to ${\text{R}}_{\text{max}}$) and calculate the corresponding volume fraction variance. The results for each of the speckle pattern analysed is shown in the plot in Fig. 5.

In the left panel of Fig. 5 we compare the volume fraction variance calculated for synthetic ASP (red dashed curve), the synthetic random speckle pattern (blue dotted curve) and the_{intensity corrected ASP (yellow curve). We observe that while for the random distribution ${\text{R}}^{2}{\text{\sigma}}_{\text{F}}^{2}\left(\text{R}\right)$ tends to diverge with R, it tends to decay for the ASP. On the contrary, in case of theintensity corrected ASP ${\text{R}}^{2}{\text{\sigma}}_{\text{F}}^{2}\left(\text{R}\right)$ decays to 0 with R, as required for hyperuniform distributions. Indeed, the amorphous speckle pattern is characterized by the hyperuniform distribution of the speckle grains, in practice the number variance ${\text{\sigma}}_{\text{N}}^{2}\left(\text{R}\right)$ of the ASP’s grains within an observation window of size $\text{R}$ grows more slowly than the sub-windows volume in the large-$\text{R}$ limit, i.e., slower than ${\text{R}}^{\text{d}}$ (in our case $\text{d}=2$). On the contrary, for the standard speckle pattern grains distribution ${\text{\sigma}}_{\text{N}}^{2}\left(\text{R}\right)\text{}$ grows quadratic with R (i.e., non-hyperuniform).}

The same comparison is shown in the right panel of Fig. 5 for the experimental patterns, which appear in agreement with the numerical simulation.

## 4. Summary and outlook

By filtering with an annular aperture the spatial frequencies of a random distribution (in our case a standard speckle pattern) we obtain an amorphous pattern. The ASP emerges from the superposition of Bessel beams and always exhibits hyperuniformity in its speckle grain distribution. In essence, we bring to the scientific community the two following results:

- I) we demonstrate that the maxima of the ASP distribution are arranged in hyperuniform fashion
- II) We demonstrated that starting from the ASP it is possible to produce (after proper intensity filtering) an intensity distribution fulfilling requirements for the continuous variable hyperuniformity and consisting of a random superposition of Bessel beams.

The first result may be directly applied with multi-photon fabrication methods in which polymerization only happens above a threshold [53, 54] thus in principle enabling to produce a hyperuniform distribution with defects located at the high intensity locations of the speckle pattern, the speckle grains. It follows that hyperuniform light fields can be used to form engineered photonic architectures capable of strong light confinements [25–27, 30] and to tailor the transport of matter waves in disordered optical potentials [17]. Continuous hyperuniform fields instead would require an addition experimental step to be directly employed: a non-homogeneous filter which may be in principle realized employing a spatial light modulator.

We believe that the presented results will inspire new fabrications strategies based on complex light patterns with the potential also for other applications ranging from ultra-cold atoms to optical imaging.

## 5. Methods and validations

#### 5.1 Description of the experimental systems

Several methods are today available to generate amorphous speckle patterns exploiting expensive digital devices such as Spatial light Modulators (SLM). In our experiments speckle patterns are produced with the setup depicted in Fig. 6(a) resorting to cheap passive optical elements: scattering medium, annular rings and lenses. The beam emitted at 633nm wavelength by a He-Ne laser is reduced to 1mm in waist by the iris (F) and impinges on to a ground diffuser (D). Trespassing D, the light experiences multiple scattering and results completely scrambled in phases and directions. A 4f system composed of lenses L1 (focal length f1 = 25,4mm) and L2 (focal length f2 = 250mm) magnifies 10X the plane at 1mm from the back surface of D and projects it onto the CCD camera. It results the speckle pattern ${\text{C}}_{\text{experiment}}$ reported in the right panel of Fig. 7(b).

In order to generate ASPs, we apply a spatial frequency selection to the light transmitted through D, a filtering operation performed by an annular aperture. As depicted in Fig. 6(b), we center on the Fourier plane of the 4f system (25,4mm at the back of L1) a commercial ring filter (RF). The RF is the R1CA2000 annular aperture from Thorlabs with $1700\text{\mu m}$ obstruction diameter and $2000\text{\mu m}$ pinhole diameter. Exploiting a ring shaped filter we generate a superposition of plane waves with fixed k and random phase, which is known to give rise to an ASP as the one shown in Fig. 2(c).

#### 5.2 The numerical method and its analysis

The speckle phenomenon results from the superposition of a multitude of randomly phased elementary components (dephasor) and can be simulated implementing a simple numerical model that sums the dephasors on a hypothetical observation plane [50]. Within this section we validate our numerical model, and we aim to compare the speckle pattern resulting from our simulations with those from the experiments considering their well known statistics [50]. First, we simulate the light speckles formation, the intensity pattern on the camera plane (the observation plane), may be described as a random superposition of Gaussian components. Following this assumption, in our simulation we impose a superposition of dephasors with Gaussian spatial profile (Gaussian modes). Given a [X,Y] = 5000x5000 pixels grid, a single mode can be described as follow:

## Funding

FP7 Marie Curie ITN “OILTEBIA” (PITN-GA-2012-317526); H2020 “Laserlab Europe” (EC-GA 654148).

## Acknowledgment

We wish to acknowledge Prof. Salvatore Torquato and His graduate students Duyu Chen, Jaeuk Kim and Zheng Ma from Princeton University (New Jersey, USA) for providing us with hyeruniform matrix from Ref. [36] and their valuable feedbacks during the development of this work.

## References and links

**1. **M.J. Ziman, “Models of disorder: the theoretical physics of homogeneously disordered systems,” CUP Archive, 1979.

**2. **G. Ruocco, B. Abaie, W. Schirmacher, A. Mafi, and M. Leonetti, “Disorder-induced single-mode transmission,” Nat. Commun. **8**(14571), 14571 (2017). [CrossRef] [PubMed]

**3. **E. G. van Putten, D. Akbulut, J. Bertolotti, W. L. Vos, A. Lagendijk, and A. P. Mosk, “Scattering lens resolves sub-100 nm structures with visible light,” Phys. Rev. Lett. **106**(19), 193905 (2011). [CrossRef] [PubMed]

**4. **Y. Akahane, T. Asano, B. S. Song, and S. Noda, “High-Q photonic nanocavity in a two-dimensional photonic crystal,” Nature **425**(6961), 944–947 (2003). [CrossRef] [PubMed]

**5. **P. D. García, R. Sapienza, and C. López, “Photonic glasses: a step beyond white paint,” Adv. Mater. **22**(1), 12–19 (2010). [CrossRef] [PubMed]

**6. **D. S. Wiersma, “Disordered photonics,” Nat. Photonics **7**(3), 188–196 (2013). [CrossRef]

**7. **S. John, “Strong localization of photons in certain disordered dielectric superlattices,” Phys. Rev. Lett. **58**(23), 2486–2489 (1987). [CrossRef] [PubMed]

**8. **M. Florescu, S. Torquato, and P. J. Steinhardt, “Designer disordered materials with large, complete photonic band gaps,” Proc. Natl. Acad. Sci. U.S.A. **106**(49), 20658–20663 (2009). [CrossRef] [PubMed]

**9. **S. Yu, X. Piao, J. Hong, and N. Park, “Metadisorder for designer light in random systems,” Sci. Adv. **2**(10), e1501851 (2016). [CrossRef] [PubMed]

**10. **F. Riboli, N. Caselli, S. Vignolini, F. Intonti, K. Vynck, P. Barthelemy, A. Gerardino, L. Balet, L. H. Li, A. Fiore, M. Gurioli, and D. S. Wiersma, “Engineering of light confinement in strongly scattering disordered media,” Nat. Mater. **13**(7), 720–725 (2014). [CrossRef] [PubMed]

**11. **M. Jang, Y. Horie, A. Shibukawa, J. Brake, Y. Liu, S. M. Kamali, A. Arbabi, H. Ruan, A. Faraon, and C. Yang, “Wavefront shaping with disorder-engineered metasurfaces,” Nat. Photonics **12**(2), 84–90 (2018). [CrossRef] [PubMed]

**12. **P. Hsieh, C. Chung, J. F. McMillan, M. Tsai, M. Lu, N. C. Panoiu, and C. W. Wong, “Photon transport enhanced by transverse Anderson localization in disordered superlattices,” Nat. Phys. **11**(3), 268–274 (2015). [CrossRef]

**13. **R. Fickler, M. Ginoya, and R. W. Boyd, “Custom-tailored spatial mode sorting by controlled random scattering,” Phys. Rev. B **95**(16), 161108 (2017). [CrossRef]

**14. **D. Di Battista, D. Ancora, M. Leonetti, and G. Zacharakis, “Tailoring non-diffractive beams from amorphous light speckles,” Appl. Phys. Lett. **109**(12), 121110 (2016). [CrossRef]

**15. **R. Fischer, I. Vidal, D. Gilboa, R. R. B. Correia, A. C. Ribeiro-Teixeira, S. D. Prado, J. Hickman, and Y. Silberberg, “Light with tunable non-Markovian phase imprint,” Phys. Rev. Lett. **115**(7), 073901 (2015). [CrossRef] [PubMed]

**16. **Y. Bromberg and H. Cao, “Generating non-Rayleigh speckles with tailored intensity statistics,” Phys. Rev. Lett. **112**(21), 213904 (2014). [CrossRef]

**17. **L. Sanchez-Palencia and M. Lewenstein, “Disordered quantum gases under control,” Nat. Phys. **6**(2), 87–95 (2010). [CrossRef]

**18. **G. Volpe, G. Volpe, and S. Gigan, “Brownian motion in a speckle light field: tunable anomalous diffusion and selective optical manipulation,” Sci. Rep. **4**(3936), 3936 (2014). [PubMed]

**19. **G. Volpe, L. Kurz, A. Callegari, G. Volpe, and S. Gigan, “Speckle optical tweezers: micromanipulation with random light fields,” Opt. Express **22**(15), 18159–18167 (2014). [CrossRef] [PubMed]

**20. **F. Jendrzejewski, A. Bernard, K. Mueller, P. Cheinet, V. Josse, M. Piraud, L. Pezzé, L. Sanchez-Palencia, A. Aspect, and P. Bouyer, “Three-dimensional localization of ultracold atoms in an optical disordered potential,” Nat. Phys. **8**(5), 398–403 (2012). [CrossRef]

**21. **E. Mudry, K. Belkebir, J. Girard, J. Savatier, E. Le Moal, C. Nicoletti, M. Allain, and A. Sentenac, “Structured illumination microscopy using unknown speckle patterns,” Nat. Photonics **6**(5), 312–315 (2012). [CrossRef]

**22. **D. B. Phillips, R. He, Q. Chen, G. M. Gibson, and M. J. Padgett, “Non-diffractive computational ghost imaging,” Opt. Express **24**(13), 14172–14182 (2016). [CrossRef] [PubMed]

**23. **P. Ni, P. Zhang, X. Qi, J. Yang, Z. Chen, and W. Man, “Light localization and nonlinear beam transmission in specular amorphous photonic lattices,” Opt. Express **24**(3), 2420–2426 (2016). [CrossRef] [PubMed]

**24. **D. G. Pires, A. F. Sonsin, A. J. Jesus-Silva, and E. J. S. Fonseca, “Three-Dimensional Speckle Light Self-Healing-Based Imaging System,” Sci. Rep. **8**(1), 563 (2018). [CrossRef] [PubMed]

**25. **M. Boguslawski, S. Brake, J. Armijo, F. Diebel, P. Rose, and C. Denz, “Analysis of transverse Anderson localization in refractive index structures with customized random potential,” Opt. Express **21**(26), 31713–31724 (2013). [CrossRef] [PubMed]

**26. **J. Bingi and V. M. Murukeshan, “Speckle lithography for fabricating Gaussian, quasi-random 2D structures and black silicon structures,” Sci. Rep. **5**(1), 18452 (2016). [CrossRef] [PubMed]

**27. **G. M. Conley, M. Burresi, F. Pratesi, K. Vynck, and D. S. Wiersma, “Light transport and localization in two-dimensional correlated disorder,” Phys. Rev. Lett. **112**(14), 143901 (2014). [CrossRef] [PubMed]

**28. **C. Liu, R. E. Van Der Wel, N. Rotenberg, L. Kuipers, T. F. Krauss, A. Di Falco, and A. Fratalocchi, “Triggering extreme events at the nanoscale in photonic seas,” Nat. Phys. **11**(4), 358–363 (2015). [CrossRef]

**29. **C. Vanneste, P. Sebbah, and H. Cao, “Lasing with resonant feedback in weakly scattering random systems,” Phys. Rev. Lett. **98**(14), 143902 (2007). [CrossRef] [PubMed]

**30. **M. Leonetti, C. Conti, and C. Lopez, “The mode-locking transition of random lasers,” Nat. Photonics **5**(10), 615–617 (2011). [CrossRef]

**31. **T. Amoah and M. Florescu, “High-Q optical cavities in hyperuniform disordered materials,” Phys. Rev. B **91**(2), 020201 (2015). [CrossRef]

**32. **M. Florescu, P. J. Steinhardt, and S. Torquato, “Optical cavities and waveguides in hyperuniform disordered photonic solids,” Phys. Rev. B **87**(16), 165116 (2013). [CrossRef]

**33. **S. Torquato, “Hyperuniform states of matter,” Phys. Rep. in press (2018).

**34. **E. C. Zachary and S. Torquato, “Hyperuniformity in point patterns and two-phase random heterogeneous media,” J. Stat. Mech. **12**(12), 12015 (2009). [CrossRef]

**35. **S. Torquato, “Hyperuniformity and its generalizations,” Phys. Rev. E **94**(2), 022122 (2016). [CrossRef] [PubMed]

**36. **S. Torquato, G. Zhang, and F. H. Stillinger, “Ensemble theory for stealthy hyperuniform disordered ground states,” Phys. Rev. X **5**(2), 021020 (2015). [CrossRef]

**37. **S. Atkinson, G. Zhang, A. B. Hopkins, and S. Torquato, “Critical slowing down and hyperuniformity on approach to jamming,” Phys Rev E **94**(1), 012902 (2016). [CrossRef] [PubMed]

**38. **I. Lesanovsky and J. P. Garrahan, “Out-of-equilibrium structures in strongly interacting Rydberg gases with dissipation,” Phys. Rev. A **90**(1), 011603 (2014). [CrossRef]

**39. **J. H. Weijs, R. Jeanneret, R. Dreyfus, and D. Bartolo, “Emergent Hyperuniformity in Periodically Driven Emulsions,” Phys. Rev. Lett. **115**(10), 108301 (2015). [CrossRef] [PubMed]

**40. **S. Torquato, “Disordered hyperuniform heterogeneous materials,” J. Phys. Condens. Matter **28**(41), 414012 (2016). [CrossRef] [PubMed]

**41. **Z. Ma and S. Torquato, “Random Scalar Fields and Hyperuniformity,” J. Appl. Phys. **121**(24), 244904 (2017). [CrossRef]

**42. **W. Man, M. Florescu, E. P. Williamson, Y. He, S. R. Hashemizad, B. Y. C. Leung, D. R. Liner, S. Torquato, P. M. Chaikin, and P. J. Steinhardt, “Isotropic band gaps and freeform waveguides observed in hyperuniform disordered photonic solids,” Proc. Natl. Acad. Sci. U.S.A. **110**(40), 15886–15891 (2013). [CrossRef] [PubMed]

**43. **N. Muller, J. Haberko, C. Marichy, and F. Scheffold, “Silicon hyperuniform disordered photonic materials with a pronounced gap in the shortwave infrared,” Advanced Optical Materials **2**(2), 115–119 (2014). [CrossRef]

**44. **C. De Rosa, F. Auriemma, C. Diletto, R. Di Girolamo, A. Malafronte, P. Morvillo, G. Zito, G. Rusciano, G. Pesce, and A. Sasso, “Toward hyperuniform disordered plasmonic nanostructures for reproducible surface-enhanced Raman spectroscopy,” Phys. Chem. Chem. Phys. **17**(12), 8061–8069 (2015). [CrossRef] [PubMed]

**45. **R. Degl’Innocenti, Y. D. Shah, L. Masini, A. Ronzani, A. Pitanti, Y. Ren, D. S. Jessop, A. Tredicucci, H. E. Beere, and D. A. Ritchie, “Hyperuniform disordered terahertz quantum cascade laser,” Sci. Rep. **6**(1), 19325 (2016). [CrossRef] [PubMed]

**46. **O. Leseur, R. Pierrat, and R. Carminati, “High-density hyperuniform materials can be transparent,” Optica **3**(7), 763 (2016). [CrossRef]

**47. **L. S. Froufe-Pérez, M. Engel, P. F. Damasceno, N. Muller, J. Haberko, S. C. Glotzer, and F. Scheffold, “Role of short-range order and hyperuniformity in the formation of band gaps in disordered photonic materials,” Phys. Rev. Lett. **117**(5), 053902 (2016). [CrossRef] [PubMed]

**48. **S. R. Sellers, W. Man, S. Sahba, and M. Florescu, “Local self-uniformity in photonic networks,” Nat. Commun. **8**, 14439 (2017). [CrossRef] [PubMed]

**49. **N. Muller, J. Haberko, C. Marichy, and F. Scheffold, “Photonic hyperuniform networks obtained by silicon double inversion of polymer templates,” Optica **4**(3), 361–366 (2017). [CrossRef]

**50. **J. Goodman, *Speckle Phenomena in Optics: Theory and Applications* (Englewood: Roberts and Company Publishers, 2007).

**51. **O. Katz, P. Heidmann, M. Fink, and S. Gigan, “Non-invasive single-shot imaging through scattering layers and around corners via speckle correlations,” Nat. Photonics **8**(10), 784–790 (2014). [CrossRef]

**52. **J. Gateau, H. Rigneault, and M. Guillon, “Complementary Speckle Patterns: deterministic interchange of intrinsic vortices and maxima through Scattering Media,” Phys. Rev. Lett. **118**(4), 043903 (2017). [CrossRef] [PubMed]

**53. **S. Maruo, O. Nakamura, and S. Kawata, “Three-dimensional microfabrication with two-photon-absorbed photopolymerization,” Opt. Lett. **22**(2), 132–134 (1997). [CrossRef] [PubMed]

**54. **A. P. Joglekar, H. H. Liu, E. Meyhöfer, G. Mourou, and A. J. Hunt, “Optics at critical intensity: Applications to nanomorphing,” Proc. Natl. Acad. Sci. U.S.A. **101**(16), 5856–5861 (2004). [CrossRef] [PubMed]