Isotropic band gaps and freeform waveguides observed in hyperuniform disordered photonic solids

Finite-size effects in hyperuniform vortex matter

Novel hyperuniform materials are emerging as an active field of applied and basic research since they can be designed to have exceptional physical properties. This ubiquitous state of matter presents a hidden order that is characterized by the density of constituents of the system being uniform at large scales, as in a perfect crystal, although they can be isotropic and disordered like a liquid. In the quest for synthesizing hyperuniform materials in experimental conditions, the impact of finite-size effects remains as an open question to be addressed. We use vortex matter in type-II superconductors as a toy model system to study this issue. We previously reported that vortex matter nucleated in samples with point disorder is effectively hyperuniform and thus presents the interesting physical properties inherent to hyperuniform systems. In this work we present experimental evidence that on decreasing the thickness of the vortex system its hyperuniform order is depleted. By means of hydrodynamic arguments we show that the experimentally observed depletion can be associated to two crossovers that we describe within a hydrodynamic approximation. The first crossover length is thickness-dependent and separates a class-II hyperuniform regime at intermediate lengthscales from a regime that can become asymptotically non-hyperuniform for large wavelengths in very thin samples. The second crossover takes place at smaller lengthscales and marks the onset of a faster increase of density fluctuations due to the dispersivity of the elastic constants. Our work points to a novel mechanism of emerging hyperuniformity controlled by the thickness of the host sample, an issue that has to be taken into account when growing hyperuniform structures for technological applications.

Vibrational properties of disordered stealthy hyperuniform 1D atomic chains

Disorder hyperuniformity is a recently discovered exotic state of many-body systems that possess a hidden order in between that of a perfect crystal and a completely disordered system. Recently, this novel disordered state has been observed in a number of quantum materials including amorphous 2D graphene and silica, which are endowed with unexpected electronic transport properties. Here, we numerically investigate 1D atomic chain models, including perfect crystalline, disordered stealthy hyperuniform (SHU) as well as randomly perturbed atom packing configurations to obtain a quantitative understanding of how the unique SHU disorder affects the vibrational properties of these low-dimensional materials. We find that the disordered SHU chains possess lower cohesive energies compared to the randomly perturbed chains, implying their potential reliability in experiments. Our inverse partition ratio (IPR) calculations indicate that the SHU chains can support fully delocalized states just like perfect crystalline chains over a wide range of frequencies, i.e.ω∈(0,100)cm-1, suggesting superior phonon transport behaviors within these frequencies, which was traditionally considered impossible in disordered systems. Interestingly, we observe the emergence of a group of highly localized states associated withω∼200cm-1, which is characterized by a significant peak in the IPR and a peak in phonon density of states at the corresponding frequency, and is potentially useful for decoupling electron and phonon degrees of freedom. These unique properties of disordered SHU chains have implications in the design and engineering of novel quantum materials for thermal and phononic applications.

Hyperuniformity in two-dimensional periodic and quasiperiodic point patterns

We study hyperuniform properties in various two-dimensional periodic and quasiperiodic point patterns. Using the histogram of the two-point distances, we develop an efficient method to calculate the hyperuniformity order metric, which quantifies the regularity of the hyperuniform point patterns. The results are compared with those calculated with the conventional running average method. To discuss how the lattice symmetry affects the order metric, we treat the trellis and Shastry-Sutherland lattices with the same point density as examples of periodic lattices, and Stampfli hexagonal and dodecagonal quasiperiodic tilings with the same point density as examples of quasiperiodic tilings. It is found that the order metric for the Shastry-Sutherland lattice (Stampfli dodecagonal tilings) is smaller than the other in the periodic (quasiperiodic) tiling, meaning that the order metric is deeply related to the lattice symmetry. Namely, the point pattern with higher symmetry is characterized by the smaller order metric when their point densities are identical. Order metrics for several other quasiperiodic tilings are also calculated.

Experimental evidence of isotropic transparency and complete band gap formation for ultrasound propagation in stealthy hyperuniform media

Following on recent experimental characterization of the transport properties of stealthy hyperuniform media for electromagnetic and acoustic waves, we report here measurements at ultrasonic frequencies of the multiple scattering of waves by 2D hyperuniform distributions of steel rods immersed in water. The transparency, for which the effective attenuation of the medium is canceled, is first evidenced by measuring the transmission of a plane wave propagating in a highly correlated and relatively dense medium. It is shown that a band gap occurs in the vicinity of the first Bragg frequency. The isotropy of both transparency and band gap are also evidenced for the case of waves generated by a point source in differently ordered and circular-shaped distributions. In other words, we thus obtain a representation of the Green's function. Our results demonstrate the huge potential of hyperuniform, as well as highly correlated, media for the design of functional materials.

How does a hyperuniform fluid freeze?

All phase transitions can be categorized into two different types: continuous and discontinuous phase transitions. Discontinuous phase transitions are normally accompanied with significant structural changes, and nearly all of them have the kinetic pathway of nucleation and growth, if the system does not suffer from glassy dynamics. Here, in a system of barrier-controlled reactive particles, we find that the discontinuous freezing transition of a nonequilibrium hyperuniform fluid into an absorbing state does not have the kinetic pathway of nucleation and growth, and the transition is triggered by long-wavelength fluctuations. The transition rate decreases with increasing the system size, which suggests that the metastable hyperuniform fluid could be kinetically stable in an infinitely large system. This challenges the common understanding of metastability in discontinuous phase transitions. Moreover, we find that the "metastable yet kinetically stable" hyperuniform fluid features a scaling in the structure factor [Formula: see text] in 2D, which is the third dynamic hyperuniform state in addition to the critical hyperuniform state with [Formula: see text] and the nonequilibrium hyperuniform fluid with [Formula: see text].

Computational design of anisotropic stealthy hyperuniform composites with engineered directional scattering properties

Disordered hyperuniform materials are an emerging class of exotic amorphous states of matter that endow them with singular physical properties, including large isotropic photonic band gaps, superior resistance to fracture, and nearly optimal electrical and thermal transport properties, to name but a few. Here we generalize the Fourier-space-based numerical construction procedure for designing and generating digital realizations of isotropic disordered hyperuniform two-phase heterogeneous materials (i.e., composites) developed by Chen and Torquato [Acta Mater. 142, 152 (2018)1359-645410.1016/j.actamat.2017.09.053] to anisotropic microstructures with targeted spectral densities. Our generalized construction procedure explicitly incorporates the vector-dependent spectral density function χ[over ̃]_{_{V}}(k) of arbitrary form that is realizable. We demonstrate the utility of the procedure by generating a wide spectrum of anisotropic stealthy hyperuniform microstructures with χ[over ̃]_{_{V}}(k)=0 for k∈Ω, i.e., complete suppression of scattering in an "exclusion" region Ω around the origin in Fourier space. We show how different exclusion-region shapes with various discrete symmetries, including circular-disk, elliptical-disk, square, rectangular, butterfly-shaped, and lemniscate-shaped regions of varying size, affect the resulting statistically anisotropic microstructures as a function of the phase volume fraction. The latter two cases of Ω lead to directionally hyperuniform composites, which are stealthy hyperuniform only along certain directions and are nonhyperuniform along others. We find that while the circular-disk exclusion regions give rise to isotropic hyperuniform composite microstructures, the directional hyperuniform behaviors imposed by the shape asymmetry (or anisotropy) of certain exclusion regions give rise to distinct anisotropic structures and degree of uniformity in the distribution of the phases on intermediate and large length scales along different directions. Moreover, while the anisotropic exclusion regions impose strong constraints on the global symmetry of the resulting media, they can still possess structures at a local level that are nearly isotropic. Both the isotropic and anisotropic hyperuniform microstructures associated with the elliptical-disk, square, and rectangular Ω possess phase-inversion symmetry over certain range of volume fractions and a percolation threshold ϕ_{c}≈0.5. On the other hand, the directionally hyperuniform microstructures associated with the butterfly-shaped and lemniscate-shaped Ω do not possess phase-inversion symmetry and percolate along certain directions at much lower volume fractions. We also apply our general procedure to construct stealthy nonhyperuniform systems. Our construction algorithm enables one to control the statistical anisotropy of composite microstructures via the shape, size, and symmetries of Ω, which is crucial to engineering directional optical, transport, and mechanical properties of two-phase composite media.

Computational approach for structure generation of anisotropic particles (CASGAP) with targeted distributions of particle design and orientational order

The macroscopic properties of materials are governed by their microscopic structure which depends on the materials' composition (i.e., building blocks) and processing conditions. In many classes of synthetic, bioinspired, or natural soft and/or nanomaterials, one can find structural anisotropy in the microscopic structure due to anisotropic building blocks and/or anisotropic domains formed through the processing conditions. Experimental characterization and complementary physics-based or data-driven modeling of materials' structural anisotropy are critical for understanding structure-property relationships and enabling targeted design of materials with desired macroscopic properties. In this pursuit, to interpret experimentally obtained characterization results (e.g., scattering profiles) of soft materials with structural anisotropy using data-driven computational approaches, there is a need for creating real space three-dimensional structures of the designer soft materials with realistic physical features (e.g., dispersity in building block sizes) and anisotropy (i.e., aspect ratios of the building blocks, their orientational and positional order). These real space structures can then be used to compute and complement experimentally obtained characterization results or be used as initial configurations for physics-based simulations/calculations that can then provide training data for machine learning models. To address this need, we present a new computational approach called CASGAP - Computational Approach for Structure Generation of Anisotropic Particles - for generating any desired three dimensional real-space structure of anisotropic building blocks (modeled as particles) adhering to target distributions of particle shape, size, and positional and orientational order.

Bandgap fluctuations and robustness in two-dimensional hyperuniform dielectric materials

We numerically study the statistical fluctuations of photonic band gaps in ensembles of stealthy hyperuniform disordered patterns. We find that at low stealthiness, where correlations are weak, band gaps of different system realizations appear over a wide frequency range, are narrow, and generally do not overlap. Interestingly, above a critical value of stealthiness χ≳0.35, the bandgaps become large and overlap significantly from realization to realization, while a second gap appears. These observations extend our understanding of photonic bandgaps in disordered systems and provide information on the robustness of gaps in practical applications.

Evolving scattering networks for engineering disorder

Network science provides a powerful tool for unraveling the complexities of social, technological and biological systems. Constructing networks using wave phenomena is also of great interest in devising advanced hardware for machine learning, as shown in optical neural networks. Although most wave-based networks have employed static network models, the impact of evolving models in network science provides strong motivation to apply dynamical network modeling to wave physics. Here the concept of evolving scattering networks for scattering phenomena is developed. The network is defined by links, node degrees and their evolution processes modeling multi-particle interferences, which directly determine scattering from disordered materials. I demonstrate the concept by examining network-based material classification, microstructure screening and preferential attachment in evolutions, which are applied to stealthy hyperuniformity. The results enable independent control of scattering from different length scales, revealing superdense material phases in short-range order. The proposed concept provides a bridge between wave physics and network science to resolve multiscale material complexities and open-system material design.

Transport through Amorphous Photonic Materials with Localization and Bandgap Regimes

We propose a framework that unifies the description of light transmission through three-dimensional amorphous dielectric materials that exhibit both localization and a photonic bandgap. We argue that direct, coherent reflection near and in the bandgap attenuates the generation of diffuse or localized photons. Using the self-consistent theory of localization and considering the density of states of photons, we can quantitatively describe the total transmission of light for all transport regimes: transparency, light diffusion, localization, and bandgap. Comparison with numerical simulations of light transport through hyperuniform networks supports our theoretical approach.

3D Tomographic Analysis of the Order-Disorder Interplay in the Pachyrhynchus congestus mirabilis Weevil

The bright colors of Pachyrhynchus weevils originate from complex dielectric nanostructures within their elytral scales. In contrast to previous work exhibiting highly ordered single-network diamond-type photonic crystals, here, it is shown by combining optical microscopy and spectroscopy measurements with 3D focused ion beam (FIB) tomography that the blue scales of P. congestus mirabilis differ from that of an ordered diamond structure. Through the use of FIB tomography on elytral scales filled with platinum (Pt) by electron beam-assisted deposition, it is revealed that the red scales of this weevil possess a periodic diamond structure, while the network morphology of the blue scales exhibit diamond morphology only on the single scattering unit level with disorder on longer length scales. Full wave simulations performed on the reconstructed volumes indicate that this local order is sufficient to open a partial photonic bandgap even at low dielectric constant contrast between chitin and air in the absence of long-range or translational order. The observation of disordered and ordered photonic crystals within a single organism opens up interesting questions on the cellular origin of coloration and studies on bio-inspired replication of angle-independent colors.

Growth rules for irregular architected materials with programmable properties

Biomaterials display microstructures that are geometrically irregular and functionally efficient. Understanding the role of irregularity in determining material properties offers a new path to engineer materials with superior functionalities, such as imperfection insensitivity, enhanced impact absorption, and stress redirection. We uncover fundamental, probabilistic structure-property relationships using a growth-inspired program that evokes the formation of stochastic architectures in natural systems. This virtual growth program imposes a set of local rules on a limited number of basic elements. It generates materials that exhibit a large variation in functional properties starting from very limited initial resources, which echoes the diversity of biological systems. We identify basic rules to control mechanical properties by independently varying the microstructure's topology and geometry in a general, graph-based representation of irregular materials.

Control of localization and optical properties with deep-subwavelength engineered disorder

The effect of deep subwavelength disorder in one-dimensional dichromic multilayer films on the optical transmission, localization length, and Goos-Hänchen shift around the critical angle is analyzed using sets of disordered multilayer films with different degrees of order metric τ. For each Gaussian-perturbed multilayer film designed by a Metropolis algorithm targeting the predetermined order metric τ, the numerically obtained localization length and transmission show excellent agreement with the recent theoretical analysis developed for disordered multilayer films, further revealing τ-dependence of the Goos-Hänchen shift across the critical angle. Emphasizing the role of deep subwavelength structures in disorder-induced transmission enhancement, our result thus paves the way toward the inverse design of a deep subwavelength disordered structural landscape for the targeted order metric τ or abnormal optical responses - including the Goos-Hänchen shift.

Hyperuniform Active Chiral Fluids with Tunable Internal Structure

Large density fluctuations observed in active systems and hyperuniformity are two seemingly incompatible phenomena. However, the formation of hyperuniform states has been recently predicted in nonequilibrium fluids formed by chiral particles performing circular motion with the same handedness. Here we report evidence of hyperuniformity realized in a chiral active fluid comprised of pear-shaped Quincke rollers of arbitrary handedness. We show that hyperuniformity and large density fluctuations, triggered by dynamic clustering, coexist in this system at different length scales. The system loses its hyperuniformity as the curvature of particles' motion increases, transforming them into localized spinners. Our results experimentally demonstrate a novel hyperuniform active fluid and provide new insights into an interplay between chirality, activity, and hyperuniformity.

Over 65% Sunlight Absorption in a 1 μm Si Slab with Hyperuniform Texture

Thin, flexible, and invisible solar cells will be a ubiquitous technology in the near future. Ultrathin crystalline silicon (c-Si) cells capitalize on the success of bulk silicon cells while being lightweight and mechanically flexible, but suffer from poor absorption and efficiency. Here we present a new family of surface texturing, based on correlated disordered hyperuniform patterns, capable of efficiently coupling the incident spectrum into the silicon slab optical modes. We experimentally demonstrate 66.5% solar light absorption in free-standing 1 μm c-Si layers by hyperuniform nanostructuring for the spectral range of 400 to 1050 nm. The absorption equivalent photocurrent derived from our measurements is 26.3 mA/cm², which is far above the highest found in literature for Si of similar thickness. Considering state-of-the-art Si PV technologies, we estimate that the enhanced light trapping can result in a cell efficiency above 15%. The light absorption can potentially be increased up to 33.8 mA/cm² by incorporating a back-reflector and improved antireflection, for which we estimate a photovoltaic efficiency above 21% for 1 μm thick Si cells.

Ultraefficient reconstruction of effectively hyperuniform disordered biphase materials via non-Gaussian random fields

Disordered hyperuniform systems are statistically isotropic and possess no Bragg peaks like liquids and glasses, yet they suppress large-scale density fluctuations in a similar manner as in perfect crystals. The unique hyperuniform long-range order in these systems endow them with nearly optimal transport, electronic, and mechanical properties. The concept of hyperuniformity was originally introduced for many-particle systems and has subsequently been generalized to biphase heterogeneous materials such as porous media, composites, polymers, and biological tissues for unconventional property discovery. Existing methods for rendering realizations of disordered hyperuniform biphase materials reconstruction typically employ stochastic optimization such as the simulated annealing approach, which requires many iterations. Here, we propose an explicit ultraefficient method for reconstructing effectively hyperuniform biphase materials, based on the second-order non-Gaussian random fields where no additional tuning step or iteration is needed. Both the effectively hyperuniform microstructure and the latent material property field can be simultaneously generated in a single reconstruction. Moreover, our method can also incorporate hierarchical uncertainties in the heterogeneous materials, including both uncertainties in the disordered material microstructure and material property variation within each phase. The efficiency and feasibility of the proposed reconstruction method are demonstrated via a wide spectrum of examples spanning from isotropic to anisotropic, effectively hyperuniform to nonhyperuniform, and antihyperuniform systems. Our ultraefficient reconstruction method can be readily incorporated into material design, probabilistic analysis, optimization, and discovery of novel disordered hyperuniform heterogeneous materials.

Hyperuniformity and phase enrichment in vortex and rotor assemblies

Ensembles of particles rotating in a two-dimensional fluid can exhibit chaotic dynamics yet develop signatures of hidden order. Such rotors are found in the natural world spanning vastly disparate length scales — from the rotor proteins in cellular membranes to models of atmospheric dynamics. Here we show that an initially random distribution of either driven rotors in a viscous membrane, or ideal vortices with minute perturbations, spontaneously self assemble into a distinct arrangement. Despite arising from drastically different physics, these systems share a Hamiltonian structure that sets geometrical conservation laws resulting in prominent structural states. We find that the rotationally invariant interactions isotropically suppress long-wavelength fluctuations — a hallmark of a disordered hyperuniform material. With increasing area fraction, the system orders into a hexagonal lattice. In mixtures of two co-rotating populations, the stronger population will gain order from the other and both will become phase enriched. Finally, we show that classical 2D point vortex systems arise as exact limits of the experimentally accessible microscopic membrane rotors, yielding a new system through which to study topological defects.

Rotor-like dynamics is observed in many natural systems, from the rotor proteins in cellular membranes to atmospheric models. Here, the authors uncover geometrical conservation laws that limit distribution of driven rotors in a membrane or a soap film and allow to predict their structural states.

Multiphase flow detection with photonic crystals and deep learning

Multiphase flows are ubiquitous in industrial settings. It is often necessary to characterize these fluid mixtures in support of process optimization. Unfortunately, existing commercial technologies often fail to provide frequent, accurate, and cost-efficient data necessary to enable process optimization. Here we show a new physics-based concept and testing with lab and field prototypes leveraging photonic crystals for real-time characterization of multiphase flows. In particular, low power (~1 mW) microwave transmission through photonic crystals filled with fluid mixtures may be interrogated by deep learning analysis techniques to provide a fast and accurate characterization of phase fraction and flow morphology. Moreover when these flow characteristics are known, the flow rate is accurately inferred from the differential pressure necessary for the flow to pass through the photonic crystal. This insight provides a basis to develop a unique class of inexpensive, accurate, and convenient techniques to characterize multiphase flows.

Critical pore radius and transport properties of disordered hard- and overlapping-sphere models

Transport properties of porous media are intimately linked to their pore-space microstructures. We quantify geometrical and topological descriptors of the pore space of certain disordered and ordered distributions of spheres, including pore-size functions and the critical pore radius δ_{c}. We focus on models of porous media derived from maximally random jammed sphere packings, overlapping spheres, equilibrium hard spheres, quantizer sphere packings, and crystalline sphere packings. For precise estimates of the percolation thresholds, we use a strict relation of the void percolation around sphere configurations to weighted bond percolation on the corresponding Voronoi networks. We use the Newman-Ziff algorithm to determine the percolation threshold using universal properties of the cluster size distribution. The critical pore radius δ_{c} is often used as the key characteristic length scale that determines the fluid permeability k. A recent study [Torquato, Adv. Wat. Resour. 140, 103565 (2020)10.1016/j.advwatres.2020.103565] suggested for porous media with a well-connected pore space an alternative estimate of k based on the second moment of the pore size 〈δ^{2}〉, which is easier to determine than δ_{c}. Here, we compare δ_{c} to the second moment of the pore size 〈δ^{2}〉, and indeed confirm that, for all porosities and all models considered, δ_{c}^{2} is to a good approximation proportional to 〈δ^{2}〉. However, unlike 〈δ^{2}〉, the permeability estimate based on δ_{c}^{2} does not predict the correct ranking of k for our models. Thus, we confirm 〈δ^{2}〉 to be a promising candidate for convenient and reliable estimates of the fluid permeability for porous media with a well-connected pore space. Moreover, we compare the fluid permeability of our models with varying degrees of order, as measured by the τ order metric. We find that (effectively) hyperuniform models tend to have lower values of k than their nonhyperuniform counterparts. Our findings could facilitate the design of porous media with desirable transport properties via targeted pore statistics.

Scalable Disordered Hyperuniform Architectures via Nanoimprint Lithography of Metal Oxides

Fabrication and scaling of disordered hyperuniform materials remain hampered by the difficulties in controlling the spontaneous phenomena leading to this novel kind of exotic arrangement of objects. Here, we demonstrate a hybrid top-down/bottom-up approach based on sol-gel dip-coating and nanoimprint lithography for the faithful reproduction of disordered hyperuniform metasurfaces in metal oxides. Nano- to microstructures made of silica and titania can be directly printed over several cm² on glass and on silicon substrates. First, we describe the polymer mold fabrication starting from a hard master obtained via spontaneous solid-state dewetting of SiGe and Ge thin layers on SiO₂. Then, we assess the effective disordered hyperuniform character of master and replica and the role of the thickness of the sol-gel layer on the metal oxide replicas and on the presence of a residual layer underneath. Finally, as a potential application, we show the antireflective character of titania structures on silicon. Our results are relevant for the realistic implementation over large scales of disordered hyperuniform nano- and microarchitectures made of metal oxides, thus opening their exploitation in the framework of wet chemical assembly.

Soft and disordered hyperuniform elastic metamaterials for highly efficient vibration concentration

Vibrations, which widely exist throughout the world, could be a nearly endless and locally obtained green energy source. It has been a long-standing challenge to efficiently utilize dispersed vibration energy, especially within the high-frequency range, since the amplitudes of high-frequency vibrations in local parts of objects are relatively weak. Here, for the first time, we propose a soft and disordered hyperuniform elastic metamaterial (DHEM), achieving a remarkable concentration of vibrations in broad frequency bands by a maximum enhancement factor of ∼4000 at 1930 Hz. The DHEM, with rational sizes from ∼1 cm to ∼1000 cm, covers a broad range of frequencies from ∼10 Hz to ∼10 kHz, which are emitted by many vibration sources including domestic appliances, factories and transportation systems, for example. Moreover, the performance of the soft DHEM under deformation is validated, enabling conformal attachments on uneven objects. Our findings lay the groundwork for reducing traditional energy consumption by recovering some of the energy dissipated by devices in the working world.

Dynamic evolution of hyperuniformity in a driven dissipative colloidal system

Hyperuniformity is evolving to become a unifying concept that can help classify and characterize equilibrium and nonequilibrium states of matter. Therefore, understanding the extent of hyperuniformity in dissipative systems is critical. Here, we study the dynamic evolution of hyperuniformity in a driven dissipative colloidal system. We experimentally show and numerically verify that the hyperuniformity of a colloidal crystal is robust against various lattice imperfections and environmental perturbations. This robustness even manifests during crystal disassembly as the system switches between strong (class I), logarithmic (class II), weak (class III), and non-hyperuniform states. To aid analyses, we developed a comprehensive computational toolbox, enabling real-time characterization of hyperuniformity in real- and reciprocal-spaces together with the evolution of several order metric features, and measurements showing the effect of external perturbations on the spatiotemporal distribution of the particles. Our findings provide a new framework to understand the basic principles that drive a dissipative system to a hyperuniform state.

From colloidal particles to photonic crystals: advances in self-assembly and their emerging applications

Over the last three decades, photonic crystals (PhCs) have attracted intense interests thanks to their broad potential applications in optics and photonics. Generally, these structures can be fabricated via either "top-down" lithographic or "bottom-up" self-assembly approaches. The self-assembly approaches have attracted particular attention due to their low cost, simple fabrication processes, relative convenience of scaling up, and the ease of creating complex structures with nanometer precision. The self-assembled colloidal crystals (CCs), which are good candidates for PhCs, have offered unprecedented opportunities for photonics, optics, optoelectronics, sensing, energy harvesting, environmental remediation, pigments, and many other applications. The creation of high-quality CCs and their mass fabrication over large areas are the critical limiting factors for real-world applications. This paper reviews the state-of-the-art techniques in the self-assembly of colloidal particles for the fabrication of large-area high-quality CCs and CCs with unique symmetries. The first part of this review summarizes the types of defects commonly encountered in the fabrication process and their effects on the optical properties of the resultant CCs. Next, the mechanisms of the formation of cracks/defects are discussed, and a range of versatile fabrication methods to create large-area crack/defect-free two-dimensional and three-dimensional CCs are described. Meanwhile, we also shed light on both the advantages and limitations of these advanced approaches developed to fabricate high-quality CCs. The self-assembly routes and achievements in the fabrication of CCs with the ability to open a complete photonic bandgap, such as cubic diamond and pyrochlore structure CCs, are discussed as well. Then emerging applications of large-area high-quality CCs and unique photonic structures enabled by the advanced self-assembly methods are illustrated. At the end of this review, we outlook the future approaches in the fabrication of perfect CCs and highlight their novel real-world applications.

Spatially Shaping Waves to Penetrate Deep inside a Forbidden Gap

It is well known that waves with frequencies within the forbidden gap inside a crystal are transported only over a limited distance-the Bragg length-before being reflected by Bragg interference. Here, we demonstrate how to send waves much deeper into crystals in an exemplary study of light in two-dimensional silicon photonic crystals. By spatially shaping the wave fronts, the internal energy density-probed via the laterally scattered intensity-is enhanced at a tunable distance away from the front surface. The intensity is up to 100× enhanced compared to random wave fronts, and extends as far as 8× the Bragg length, which agrees with an extended mesoscopic model. We thus report a novel control knob for mesoscopic wave transport that pertains to any kind of waves.

Stone-Wales defects preserve hyperuniformity in amorphous two-dimensional networks

Disordered hyperuniformity (DHU) is a recently discovered novel state of many-body systems that possesses vanishing normalized infinite-wavelength density fluctuations similar to a perfect crystal and an amorphous structure like a liquid or glass. Here, we discover a hyperuniformity-preserving topological transformation in two-dimensional (2D) network structures that involves continuous introduction of Stone-Wales (SW) defects. Specifically, the static structure factor [Formula: see text] of the resulting defected networks possesses the scaling [Formula: see text] for small wave number k, where [Formula: see text] monotonically decreases as the SW defect concentration p increases, reaches [Formula: see text] at [Formula: see text], and remains almost flat beyond this p. Our findings have important implications for amorphous 2D materials since the SW defects are well known to capture the salient feature of disorder in these materials. Verified by recently synthesized single-layer amorphous graphene, our network models reveal unique electronic transport mechanisms and mechanical behaviors associated with distinct classes of disorder in 2D materials.

Non-Gaussian tail in the force distribution: a hallmark of correlated disorder in the host media of elastic objects

Inferring the nature of disorder in the media where elastic objects are nucleated is of crucial importance for many applications but remains a challenging basic-science problem. Here we propose a method to discern whether weak-point or strong-correlated disorder dominates based on characterizing the distribution of the interaction forces between objects mapped in large fields-of-view. We illustrate our proposal with the case-study system of vortex structures nucleated in type-II superconductors with different pinning landscapes. Interaction force distributions are computed from individual vortex positions imaged in thousands-vortices fields-of-view in a two-orders-of-magnitude-wide vortex-density range. Vortex structures nucleated in point-disordered media present Gaussian distributions of the interaction force components. In contrast, if the media have dilute and randomly-distributed correlated disorder, these distributions present non-Gaussian algebraically-decaying tails for large force magnitudes. We propose that detecting this deviation from the Gaussian behavior is a fingerprint of strong disorder, in our case originated from a dilute distribution of correlated pinning centers.

A review of shaped colloidal particles in fluids: anisotropy and chirality

This review treats asymmetric colloidal particles moving through their host fluid under the action of some form of propulsion. The propulsion can come from an external body force or from external shear flow. It may also come from externally-induced stresses at the surface, arising from imposed chemical, thermal or electrical gradients. The resulting motion arises jointly from the driven particle and the displaced fluid. If the objects are asymmetric, every aspect of their motion and interaction depends on the orientation of the objects. This orientation in turn changes in response to the driving. The objects' shape can thus lead to a range of emergent anisotropic and chiral motion not possible with isotropic spherical particles. We first consider what aspects of a body's asymmetry can affect its drift through a fluid, especially chiral motion. We next discuss driving by injecting external force or torque into the particles. Then we consider driving without injecting force or torque. This includes driving by shear flow and driving by surface stresses, such as electrophoresis. We consider how time-dependent driving can induce collective orientational order and coherent motion. We show how a given particle shape can be represented using an assembly of point forces called a Stokeslet object. We next consider the interactions between anisotropic propelled particles, the symmetries governing the interactions, and the possibility of bound pairs of particles. Finally we show how the collective hydrodynamics of a suspension can be qualitatively altered by the particles' shapes. The asymmetric responses discussed here are broadly relevant also for swimming propulsion of active micron-scale objects such as microorganisms.

Impact of particle size and multiple scattering on the propagation of waves in stealthy-hyperuniform media

Propagation of waves in materials that exhibit stealthy-hyperuniform long-range correlations is investigated. By using a modal decomposition of the field that takes multiple scattering into account at all orders, we study the impact of the concentration of particles on the transparency of such materials at low frequency. An upper frequency limit for transparency is defined that include both the particle size and the degree of stealthiness. We show that the independent scattering approximation is not relevant to calculate elastic mean free paths when wavelength becomes comparable to the size of particles. We find that transparency is very robust with regard to the degree of heterogeneity of the host random medium and the polydispersity of particles. Finally, it is shown that resonances can be used as the frequency filter.

Hyperuniform Structures Formed by Shearing Colloidal Suspensions

In periodically sheared suspensions there is a dynamical phase transition, characterized by a critical strain amplitude γ_{c}, between an absorbing state where particle trajectories are reversible and an active state where trajectories are chaotic and diffusive. Repulsive nonhydrodynamic interactions between "colliding" particles' surfaces have been proposed as a source of this broken time reversal symmetry. A simple toy model called random organization qualitatively reproduces the dynamical features of this transition. Random organization and other absorbing state models exhibit hyperuniformity, a strong suppression of density fluctuations on long length scales quantified by a structure factor S(q→0)∼q^{α} with α>0, at criticality. Here we show experimentally that the particles in periodically sheared suspensions organize into structures with anisotropic short-range order but isotropic, long-range hyperuniform order when oscillatory shear amplitudes approach γ_{c}.

Machine learning identifies scale-free properties in disordered materials

The vast amount of design freedom in disordered systems expands the parameter space for signal processing. However, this large degree of freedom has hindered the deterministic design of disordered systems for target functionalities. Here, we employ a machine learning approach for predicting and designing wave-matter interactions in disordered structures, thereby identifying scale-free properties for waves. To abstract and map the features of wave behaviors and disordered structures, we develop disorder-to-localization and localization-to-disorder convolutional neural networks, each of which enables the instantaneous prediction of wave localization in disordered structures and the instantaneous generation of disordered structures from given localizations. We demonstrate that the structural properties of the network architectures lead to the identification of scale-free disordered structures having heavy-tailed distributions, thus achieving multiple orders of magnitude improvement in robustness to accidental defects. Our results verify the critical role of neural network structures in determining machine-learning-generated real-space structures and their defect immunity.

Hyperuniform Monocrystalline Structures by Spinodal Solid-State Dewetting

Materials featuring anomalous suppression of density fluctuations over large length scales are emerging systems known as disordered hyperuniform. The underlying hidden order renders them appealing for several applications, such as light management and topologically protected electronic states. These applications require scalable fabrication, which is hard to achieve with available top-down approaches. Theoretically, it is known that spinodal decomposition can lead to disordered hyperuniform architectures. Spontaneous formation of stable patterns could thus be a viable path for the bottom-up fabrication of these materials. Here, we show that monocrystalline semiconductor-based structures, in particular Si_{1-x}Ge_{x} layers deposited on silicon-on-insulator substrates, can undergo spinodal solid-state dewetting featuring correlated disorder with an effective hyperuniform character. Nano- to micrometric sized structures targeting specific morphologies and hyperuniform character can be obtained, proving the generality of the approach and paving the way for technological applications of disordered hyperuniform metamaterials. Phase-field simulations explain the underlying nonlinear dynamics and the physical origin of the emerging patterns.

Optimized Large Hyperuniform Binary Colloidal Suspensions in Two Dimensions

The creation of disordered hyperuniform materials with extraordinary optical properties (e.g., large complete photonic band gaps) requires a capacity to synthesize large samples that are effectively hyperuniform down to the nanoscale. Motivated by this challenge, we propose a feasible equilibrium fabrication protocol using binary paramagnetic colloidal particles confined in a 2D plane. The strong and long-ranged dipolar interaction induced by a tunable magnetic field is free from screening effects that attenuate long-ranged electrostatic interactions in charged colloidal systems. Specifically, we numerically find a family of optimal size ratios that makes the two-phase system effectively hyperuniform. We show that hyperuniformity is a general consequence of low isothermal compressibilities, which makes our protocol suitable to treat more general systems with other long-ranged interactions, dimensionalities, and/or polydispersity. Our methodology paves the way to synthesize large photonic hyperuniform materials that function in the visible to infrared range and hence may accelerate the discovery of novel photonic materials.

Topological Hyperbolic Lattices

Non-Euclidean geometry, discovered by negating Euclid's parallel postulate, has been of considerable interest in mathematics and related fields for the description of geographical coordinates, Internet infrastructures, and the general theory of relativity. Notably, an infinite number of regular tessellations in hyperbolic geometry-hyperbolic lattices-are expected to extend Euclidean Bravais lattices and the consequent wave phenomena to non-Euclidean geometry. However, topological states of matter in hyperbolic lattices have yet to be reported. Here we investigate topological phenomena in hyperbolic geometry, exploring how the quantized curvature and edge dominance of the geometry affect topological phases. We report a recipe for the construction of a Euclidean photonic platform that inherits the topological band properties of a hyperbolic lattice under a uniform, pseudospin-dependent magnetic field, realizing a non-Euclidean analog of the quantum spin Hall effect. For hyperbolic lattices with different quantized curvatures, we examine the topological protection of helical edge states and generalize Hofstadter's butterfly, by employing two empirical parameters that measure the edge confinement and defect immunity. We demonstrate that the proposed platforms exhibit the unique spectral-magnetic sensitivity of topological immunity in highly curved hyperbolic planes. Our approach is applicable to general non-Euclidean geometry and enables the exploitation of infinite lattice degrees of freedom for band theory.

Hyperuniformity and density fluctuations at a rigidity transition in a model of biological tissues

The suppression of density fluctuations at different length scales is the hallmark of hyperuniformity. Here, we explore the presence of this hidden order in a manybody interacting model of biological tissue, known to exhibit a transition, or sharp crossover, from a solid to a fluid like phase. We show that the density fluctuations in the rigid phase are only suppressed up to a finite lengthscale. This length scale monotonically increases and grows rapidly as we approach the fluid phase reminiscent to divergent behavior at a critical point, such that the system is effectively hyperuniform in the fluid phase. Furthermore, complementary behavior of the structure factor across the critical point also indicates that hyperuniformity found in the fluid phase is stealthy.

Observation of Protected Photonic Edge States Induced by Real-Space Topological Lattice Defects

Topological defects (TDs) in crystal lattices are elementary lattice imperfections that cannot be removed by local perturbations, due to their real-space topology. In the emerging field of topological photonics, photonic topological edge states arise from the nontrivial topology of the band structure defined in momentum space and are generally protected against defects. Here we show that adding TDs into a valley photonic crystal generates a lattice disclination that acts like a domain wall and hosts photonic topological edge states. Unlike previous topological waveguides, the disclination forms an open arc and functions as a free-form waveguide connecting a pair of TDs of opposite topological charge. This interplay between the real-space topology of lattice defects and momentum-space band topology provides a novel scheme to implement large-scale photonic structures with complex arrangements of robust topological waveguides and resonators.

Multifunctional composites for elastic and electromagnetic wave propagation

Composites are ideally suited to achieve desirable multifunctional effective properties since the best properties of different materials can be judiciously combined with designed microstructures. Here, we establish cross-property relations for two-phase composite media that link effective elastic and electromagnetic wave characteristics to one another, including the respective effective wave speeds and attenuation coefficients, which facilitate multifunctional material design. This is achieved by deriving accurate formulas for the effective electromagnetic and elastodynamic properties that depend on the wavelengths of the incident waves and the microstructure via the spectral density. Our formulas enable us to explore the wave characteristics of a broad class of disordered microstructures because they apply, unlike conventional formulas, to a wide range of incident wavelengths (i.e., well beyond the long-wavelength regime). This capability enables us to study the dynamic properties of exotic disordered "hyperuniform" composites that can have advantages over crystalline ones, such as nearly optimal, direction-independent properties and robustness against defects. We specifically show that disordered "stealthy" hyperuniform microstructures exhibit novel wave characteristics (e.g., low-pass filters that transmit waves "isotropically" up to a finite wavenumber). Our cross-property relations for the effective wave characteristics can be applied to design multifunctional composites via inverse techniques. Design examples include structural components that require high stiffness and electromagnetic absorption; heat sinks for central processing units and sound-absorbing housings for motors that have to efficiently emit thermal radiation and suppress mechanical vibrations; and nondestructive evaluation of the elastic moduli of materials from the effective dielectric response.

Hyperuniform disordered waveguides and devices for near infrared silicon photonics

We introduce a hyperuniform-disordered platform for the realization of near-infrared photonic devices on a silicon-on-insulator platform, demonstrating the functionality of these structures in a flexible silicon photonics integrated circuit platform unconstrained by crystalline symmetries. The designs proposed advantageously leverage the large, complete, and isotropic photonic band gaps provided by hyperuniform disordered structures. An integrated design for a compact, sub-volt, sub-fJ/bit, hyperuniform-clad, electrically controlled resonant optical modulator suitable for fabrication in the silicon photonics ecosystem is presented along with simulation results. We also report results for passive device elements, including waveguides and resonators, which are seamlessly integrated with conventional silicon-on-insulator strip waveguides and vertical couplers. We show that the hyperuniform-disordered platform enables improved compactness, enhanced energy efficiency, and better temperature stability compared to the silicon photonics devices based on rib and strip waveguides.

Phoamtonic designs yield sizeable 3D photonic band gaps

We show that it is possible to construct foam-based heterostructures with complete photonic band gaps. Three-dimensional foams are promising candidates for the self-organization of large photonic networks with combinations of physical characteristics that may be useful for applications. The largest band gap found is based on 3D Weaire-Phelan foam, a structure that was originally introduced as a solution to the Kelvin problem of finding the 3D tessellation composed of equal-volume cells that has the least surface area. The photonic band gap has a maximal size of 16.9% (at a volume fraction of 21.6% for a dielectric contrast [Formula: see text]) and a high degree of isotropy, properties that are advantageous in designing photonic waveguides and circuits. We also present results for 2 other foam-based heterostructures based on Kelvin and C15 foams that have somewhat smaller but still significant band gaps.

Hydrodynamics of random-organizing hyperuniform fluids

Disordered hyperuniform structures are locally random while uniform like crystals at large length scales. Recently, an exotic hyperuniform fluid state was found in several nonequilibrium systems, while the underlying physics remains unknown. In this work, we propose a nonequilibrium (driven-dissipative) hard-sphere model and formulate a hydrodynamic theory based on Navier-Stokes equations to uncover the general mechanism of the fluidic hyperuniformity (HU). At a fixed density, this model system undergoes a smooth transition from an absorbing state to an active hyperuniform fluid and then, to the equilibrium fluid by changing the dissipation strength. We study the criticality of the absorbing-phase transition. We find that the origin of fluidic HU can be understood as the damping of a stochastic harmonic oscillator in q space, which indicates that the suppressed long-wavelength density fluctuation in the hyperuniform fluid can exhibit as either acoustic (resonance) mode or diffusive (overdamped) mode. Importantly, our theory reveals that the damping dissipation and active reciprocal interaction (driving) are the two ingredients for fluidic HU. Based on this principle, we further demonstrate how to realize the fluidic HU in an experimentally accessible active spinner system and discuss the possible realization in other systems.

Hyperuniformity on spherical surfaces

We study and characterize local density fluctuations of ordered and disordered hyperuniform point distributions on spherical surfaces. In spite of the extensive literature on disordered hyperuniform systems in Euclidean geometries, to date few works have dealt with the problem of hyperuniformity in curved spaces. Indeed, some systems that display disordered hyperuniformity, like the spatial distribution of photoreceptors in avian retina, actually occur on curved surfaces. Here we will focus on the local particle number variance and its dependence on the size of the sampling window (which we take to be a spherical cap) for regular and uniform point distributions, as well as for equilibrium configurations of fluid particles interacting through Lennard-Jones, dipole-dipole, and charge-charge potentials. We show that the scaling of the local number variance as a function of the window size enables one to characterize hyperuniform and nonhyperuniform point patterns also on spherical surfaces.

Disordered hyperuniform obstacles enhance sorting of dynamically chiral microswimmers

Disordered hyperuniformity, a brand new type of arrangement with novel physical properties, provides various practical applications in extensive fields. To highlight the great potential of applying disordered hyperuniformity to active systems, a practical example is reported here by an optimal sorting of dynamically chiral microswimmers in disordered hyperuniform obstacle environments in comparison with regular or disordered ones. This optimal chirality sorting stems from a competition between advantageous microswimmer-obstacle collisions and disadvantageous trapping of microswimmers by obstacles. Based on this mechanism, optimal chirality sorting is also realized by tuning other parameters including the number density of obstacles, the strength of driven force and the noise intensity. Our findings may open a new perspective on both theoretical and experimental investigations for further applications of disordered hyperuniformity in active systems.

Methodology to construct large realizations of perfectly hyperuniform disordered packings

Disordered hyperuniform packings (or dispersions) are unusual amorphous two-phase materials that are endowed with exotic physical properties. Such hyperuniform systems are characterized by an anomalous suppression of volume-fraction fluctuations at infinitely long-wavelengths, compared to ordinary disordered materials. While there has been growing interest in such singular states of amorphous matter, a major obstacle has been an inability to produce large samples that are perfectly hyperuniform due to practical limitations of conventional numerical and experimental methods. To overcome these limitations, we introduce a general theoretical methodology to construct perfectly hyperuniform packings in d-dimensional Euclidean space R^{d}. Specifically, beginning with an initial general tessellation of space by disjoint cells that meets a "bounded-cell" condition, hard particles of general shape are placed inside each cell such that the local-cell particle packing fractions are identical to the global packing fraction. We prove that the constructed packings with a polydispersity in size are perfectly hyperuniform in the infinite-sample-size limit, regardless of particle shapes, positions, and numbers per cell. We use this theoretical formulation to devise an efficient and tunable algorithm to generate extremely large realizations of such packings. We employ two distinct initial tessellations: Voronoi as well as sphere tessellations. Beginning with Voronoi tessellations, we show that our algorithm can remarkably convert extremely large nonhyperuniform packings into hyperuniform ones in R^{2} and R^{3}. Implementing our theoretical methodology on sphere tessellations, we establish the hyperuniformity of the classical Hashin-Shtrikman multiscale coated-spheres structures, which are known to be two-phase media microstructures that possess optimal effective transport and elastic properties. A consequence of our work is a rigorous demonstration that packings that have identical tessellations can either be nonhyperuniform or hyperuniform by simply tuning local characteristics. It is noteworthy that our computationally designed hyperuniform two-phase systems can easily be fabricated via state-of-the-art methods, such as 2D photolithographic and 3D printing technologies. In addition, the tunability of our methodology offers a route for the discovery of novel disordered hyperuniform two-phase materials.

Foam as a self-assembling amorphous photonic band gap material

We show that slightly polydisperse disordered 2D foams can be used as a self-assembled template for isotropic photonic band gap (PBG) materials for transverse electric (TE) polarization. Calculations based on in-house experimental and simulated foam structures demonstrate that, at sufficient refractive index contrast, a dry foam organization with threefold nodes and long slender Plateau borders is especially advantageous to open a large PBG. A transition from dry to wet foam structure rapidly closes the PBG mainly by formation of bigger fourfold nodes, filling the PBG with defect modes. By tuning the foam area fraction, we find an optimal quantity of dielectric material, which maximizes the PBG in experimental systems. The obtained results have a potential to be extended to 3D foams to produce a next generation of self-assembled disordered PBG materials, enabling fabrication of cheap and scalable photonic devices.

Enhanced absorption of waves in stealth hyperuniform disordered media

We study the propagation of waves in a set of absorbing subwavelength scatterers positioned on a stealth hyperuniform point pattern. We show that spatial correlations in the disorder substantially enhance absorption compared to a fully disordered structure with the same density of scatterers. The non-resonant nature of the mechanism provides broad angular and spectral robustness. These results demonstrate the possibility to design low-density materials with blackbody-like absorption.

Hyperuniformity order metric of Barlow packings

The concept of hyperuniformity has been a useful tool in the study of density fluctuations at large length scales in systems ranging across the natural and mathematical sciences. One can rank a large class of hyperuniform systems by their ability to suppress long-range density fluctuations through the use of a hyperuniformity order metric Λ[over ¯]. We apply this order metric to the Barlow packings, which are the infinitely degenerate densest packings of identical rigid spheres that are distinguished by their stacking geometries and include the commonly known fcc lattice and hcp crystal. The "stealthy stacking" theorem implies that these packings are all stealthy hyperuniform, a strong type of hyperuniformity, which involves the suppression of scattering up to a wave vector K. We describe the geometry of three classes of Barlow packings, two disordered classes and small-period packings. In addition, we compute a lower bound on K for all Barlow packings. We compute Λ[over ¯] for the aforementioned three classes of Barlow packings and find that, to a very good approximation, it is linear in the fraction of fcc-like clusters, taking values between those of least-ordered hcp and most-ordered fcc. This implies that the value of Λ[over ¯] of all Barlow packings is primarily controlled by the local cluster geometry. These results highlight the special nature of anisotropic stacking disorder, which provides impetus for future research on the development of anisotropic order metrics and hyperuniformity properties.

Hyperuniformity and anti-hyperuniformity in one-dimensional substitution tilings

This work considers the scaling properties characterizing the hyperuniformity (or anti-hyperuniformity) of long-wavelength fluctuations in a broad class of one-dimensional substitution tilings. A simple argument is presented which predicts the exponent α governing the scaling of Fourier intensities at small wavenumbers, tilings with α > 0 being hyperuniform, and numerical computations confirm that the predictions are accurate for quasiperiodic tilings, tilings with singular continuous spectra and limit-periodic tilings. Quasiperiodic or singular continuous cases can be constructed with α arbitrarily close to any given value between -1 and 3. Limit-periodic tilings can be constructed with α between -1 and 1 or with Fourier intensities that approach zero faster than any power law.

Hidden Hyperuniformity in Soft Polymeric Materials

We investigate the nature of long-range density fluctuations in melts of model "soft" polymers, specifically stars and bottlebrushes, over a wide temperature range by molecular dynamics simulation. The cores of the stars and the backbones of bottlebrush polymers are found to have a hyperuniform distribution; i.e., they exhibit anomalously small density fluctuations over a wide temperature range above the glass transition temperature. The hyperuniformity of these substituent polymer subregions is hidden since the fluid as a whole does not exhibit this property. These findings offer a strategy for the practical design of hyperuniform polymeric materials.