Incompressibility of polydisperse random-close-packed colloidal particles

Formations of force network and softening of amorphous elastic materials from a coarsen-grained particle model

Amorphous materials, such as granular substances, glasses, emulsions, foams, and cells, play significant roles in various aspects of daily life, serving as building materials, plastics, food products, and agricultural items. Understanding the mechanical response of these materials to external forces is crucial for comprehending their deformation, toughness, and stiffness. Despite the recognition of the formation of force networks within amorphous materials, the mechanisms behind their formation and their impact on macroscopic physical properties remain elusive. In this study, we employ a coarse-grained particle model to investigate the mechanical response, wherein local physical properties are integrated into the softness of the particles. Our findings reveal the emergence of a chain-like force distribution, which correlates with the planar distribution of softness and heterogeneous density variations. Additionally, we observe that the amorphous material undergoes softening due to the heterogeneous distribution of softness, a phenomenon explicable through a simple theoretical framework. Moreover, we demonstrate that the ambiguity regarding the size ratio of the blob to the force network can be adjusted by the amplitude of planar fluctuations in softness, underscoring the robustness of the coarse-grained particle model.

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 of maximally random jammed packings of hyperspheres across spatial dimensions

The maximally random jammed (MRJ) state is the most random (i.e., disordered) configuration of strictly jammed (mechanically rigid) nonoverlapping objects. MRJ packings are hyperuniform, meaning their long-wavelength density fluctuations are anomalously suppressed compared to typical disordered systems, i.e., their structure factors S(k) tend to zero as the wave number |k| tends to zero. Here we show that generating high-quality strictly jammed states for Euclidean space dimensions d=3,4, and 5 is of paramount importance in ensuring hyperuniformity and extracting precise values of the hyperuniformity exponent α>0 for MRJ states, defined by the power-law behavior of S(k)∼|k|^{α} in the limit |k|→0. Moreover, we show that for fixed d it is more difficult to ensure jamming as the particle number N increases, which results in packings that are nonhyperuniform. Free-volume theory arguments suggest that the ideal MRJ state does not contain rattlers, which act as defects in numerically generated packings. As d increases, we find that the fraction of rattlers decreases substantially. Our analysis of the largest truly jammed packings suggests that the ideal MRJ packings for all dimensions d≥3 are hyperuniform with α=d-2, implying the packings become more hyperuniform as d increases. The differences in α between MRJ packings and the recently proposed Manna-class random close packed (RCP) states, which were reported to have α=0.25 in d=3 and be nonhyperuniform (α=0) for d=4 and d=5, demonstrate the vivid distinctions between the large-scale structure of RCP and MRJ states in these dimensions. Our paper clarifies the importance of the link between true jamming and hyperuniformity and motivates the development of an algorithm to produce rattler-free three-dimensional MRJ packings.

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.

Comparison of the static structure factor at long wavelengths for a dusty plasma liquid and other liquids

Especially small values of the static structure factor S(k) at long wavelengths, i.e., small k, were obtained in an analysis of experimental data, for a two-dimensional dusty plasma in its liquid state. For comparison, an analysis of S(k) data was carried out for many previously published experiments with other liquids. The latter analysis indicates that the magnitude of S(k) at small k is typically in a range 0.02-0.13. In contrast, the corresponding value for a dusty plasma liquid was found to be as small as 0.0139. Another basic finding for the dusty plasma liquid is that S(k) at small k generally increases with temperature, with its lowest value, noted above, occurring near the melting point. Simulations were carried out for the dusty plasma liquid, and their results are generally consistent with the experiment. Since a dusty plasma has a soft interparticle interaction, our findings support earlier theoretical suggestions that a useful design strategy for creating materials having exceptionally low values of S(0), so-called hyperuniform materials, is the use of a condensed material composed of particles that interact softly at their periphery.

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.

Quantifying the long-range structure of foams and other cellular patterns with hyperuniformity disorder length spectroscopy

We investigate the local and long-range structure of several space-filling cellular patterns: bubbles in a quasi-two-dimensional foam, and Voronoi constructions made around points that are uncorrelated (Poisson patterns), low discrepancy (Halton patterns), and displaced from a lattice by Gaussian noise (Einstein patterns). We study local structure with distributions of quantities including cell areas and side numbers. The former is the widest for the bubbles making foams the most locally disordered, while the latter show no major differences between the cellular patterns. To study long-range structure, we begin by representing the cellular systems as patterns of points, both unweighted and weighted by cell area. For this, foams are represented by their bubble centroids and the Voronoi constructions are represented by the centroids as well as the points from which they are created. Long-range structure is then quantified in two ways: by the spectral density, and by a real-space analog where the variance of density fluctuations for a set of measuring windows of diameter D is made more intuitive by conversion to the distance h(D) from the window boundary where these fluctuations effectively occur. The unweighted bubble centroids have h(D) that collapses for the different ages of the foam with random Poissonian fluctuations at long distances. The area-weighted bubble centroids and area-weighted Voronoi points all have constant h(D)=h_{e} for large D; the bubble centroids have the smallest value h_{e}=0.084sqrt[〈a〉], meaning they are the most uniform. Area-weighted Voronoi centroids exhibit collapse of h(D) to the same constant h_{e}=0.084sqrt[〈a〉] as for the bubble centroids. A similar analysis is performed on the edges of the cells and the spectra of h(D) for the foam edges show h(D)∼D^{1-ε} where ε=0.30±0.15.

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.

Composition fluctuations in polydisperse liquids: Glasslike effects well above the glass transition

We study a two-dimensional glass-forming system of slightly polydisperse (LJ) particles using molecular dynamics simulations and demonstrate that in the liquid regime (well above the vitrification temperature) this model shows a number of features typical of the glass transition: (i) the relation between compressibility and structure factor S(q) is strongly violated; (ii) the dynamical structure factor S(q,t) at low q shows a two-step relaxation; (iii) the time-dependent heat capacity c_{v}(t) shows a long-time power-law tail. We show that these phenomena can be rationalized with the idea of composition fluctuations and provide a quantitative theory for the effects (i) and (ii). It implies that such effects must be inherent in all polydisperse colloidal models, including binary LJ mixtures.

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.

Disordered hyperuniformity in two-dimensional amorphous silica

Disordered hyperuniformity (DHU) is a recently proposed new state of matter, which has been observed in a variety of classical and quantum many-body systems. DHU systems are characterized by vanishing infinite-wavelength normalized density fluctuations and are endowed with unique novel physical properties. Here, we report the discovery of disordered hyperuniformity in atomic-scale two-dimensional materials, i.e., amorphous silica composed of a single layer of atoms, based on spectral-density analysis of high-resolution transmission electron microscopy images. Moreover, we show via large-scale density functional theory calculations that DHU leads to almost complete closure of the electronic bandgap compared to the crystalline counterpart, making the material effectively a metal. This is in contrast to the conventional wisdom that disorder generally diminishes electronic transport and is due to the unique electron wave localization induced by the topological defects in the DHU state.

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.

Universal hidden order in amorphous cellular geometries

Partitioning space into cells with certain extreme geometrical properties is a central problem in many fields of science and technology. Here we investigate the Quantizer problem, defined as the optimisation of the moment of inertia of Voronoi cells, i.e., similarly-sized ‘sphere-like’ polyhedra that tile space are preferred. We employ Lloyd’s centroidal Voronoi diagram algorithm to solve this problem and find that it converges to disordered states associated with deep local minima. These states are universal in the sense that their structure factors are characterised by a complete independence of a wide class of initial conditions they evolved from. They moreover exhibit an anomalous suppression of long-wavelength density fluctuations and quickly become effectively hyperuniform. Our findings warrant the search for novel amorphous hyperuniform phases and cellular materials with unique physical properties.

Disordered hyperuniformity implies a hidden order on length scales that can be found in various amorphous materials. Klatt et al. analyse the evolution of random point patterns using Llyod’s algorithm and show that they converge to an effectively hyperuniform state regardless of the initial conditions.

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.

Two Diverging Length Scales in the Structure of Jammed Packings

At densities higher than the jamming transition for athermal, frictionless repulsive spheres we find two distinct length scales, both of which diverge as a power law as the transition is approached. The first, ξ_{Z}, is associated with the two-point correlation function for the number of contacts on two particles as a function of the particle separation. The second, ξ_{f}, is associated with contact-number fluctuations in subsystems of different sizes. On scales below ξ_{f}, the fluctuations are highly suppressed, similar to the phenomenon of hyperuniformity usually associated with density fluctuations. The exponents for the divergence of ξ_{Z} and ξ_{f} are different and appear to be different in two and three dimensions.

Rational design of stealthy hyperuniform two-phase media with tunable order

Disordered stealthy hyperuniform materials are exotic amorphous states of matter that have attracted recent attention because of their novel structural characteristics (hidden order at large length scales) and physical properties, including desirable photonic and transport properties. It is therefore useful to devise algorithms that enable one to design a wide class of such amorphous configurations at will. In this paper, we present several algorithms enabling the systematic identification and generation of discrete (digitized) stealthy hyperuniform patterns with a tunable degree of order, paving the way towards the rational design of disordered materials endowed with novel thermodynamic and physical properties. To quantify the degree of order or disorder of the stealthy systems, we utilize the discrete version of the τ order metric, which accounts for the underlying spatial correlations that exist across all relevant length scales in a given digitized two-phase (or, equivalently, a two-spin state) system of interest. Our results impinge on a myriad of fields, ranging from physics, materials science and engineering, visual perception, and information theory to modern data science.

Characterization of maximally random jammed sphere packings. III. Transport and electromagnetic properties via correlation functions

In the first two papers of this series, we characterized the structure of maximally random jammed (MRJ) sphere packings across length scales by computing a variety of different correlation functions, spectral functions, hole probabilities, and local density fluctuations. From the remarkable structural features of the MRJ packings, especially its disordered hyperuniformity, exceptional physical properties can be expected. Here we employ these structural descriptors to estimate effective transport and electromagnetic properties via rigorous bounds, exact expansions, and accurate analytical approximation formulas. These property formulas include interfacial bounds as well as universal scaling laws for the mean survival time and the fluid permeability. We also estimate the principal relaxation time associated with Brownian motion among perfectly absorbing traps. For the propagation of electromagnetic waves in the long-wavelength limit, we show that a dispersion of dielectric MRJ spheres within a matrix of another dielectric material forms, to a very good approximation, a dissipationless disordered and isotropic two-phase medium for any phase dielectric contrast ratio. We compare the effective properties of the MRJ sphere packings to those of overlapping spheres, equilibrium hard-sphere packings, and lattices of hard spheres. Moreover, we generalize results to micro- and macroscopically anisotropic packings of spheroids with tensorial effective properties. The analytic bounds predict the qualitative trend in the physical properties associated with these structures, which provides guidance to more time-consuming simulations and experiments. They especially provide impetus for experiments to design materials with unique bulk properties resulting from hyperuniformity, including structural-color and color-sensing applications.

Optimizing Hyperuniformity in Self-Assembled Bidisperse Emulsions

We study long range density fluctuations (hyperuniformity) in two-dimensional jammed packings of bidisperse droplets. Taking advantage of microfluidics, we systematically span a large range of size and concentration ratios of the two droplet populations. We identify various defects increasing long range density fluctuations mainly due to organization of local particle environment. By choosing an appropriate bidispersity, we fabricate materials with a high level of hyperuniformity. Interesting transparency properties of these optimized materials are established based on numerical simulations.

Microstructure and mechanical properties of hyperuniform heterogeneous materials

A hyperuniform random heterogeneous material is one in which the local volume fraction fluctuations in an observation window decay faster than the reciprocal window volume as the window size increases. Recent studies show that this class of materials are endowed with superior physical properties such as large isotropic photonic band gaps and optimal transport properties. Here we employ a stochastic optimization procedure to systematically generate realizations of hyperuniform heterogeneous materials with controllable short-range order, which is partially quantified using the two-point correlation function S_{2}(r) associated with the phase of interest. Specifically, our procedure generalizes the widely used Yeong-Torquato reconstruction procedure by including an additional constraint for hyperuniformity, i.e., the volume integral of the autocovariance function χ(r)=S_{2}(r)-ϕ^{2} over the whole space is zero. In addition, we only require the reconstructed S_{2} to match the target function up to a certain cutoff distance γ, in order to give the system sufficient degrees of freedom to satisfy the hyperuniform condition. By systematically increasing the γ value for a given S_{2}, one can produce a spectrum of hyperuniform heterogeneous materials with varying degrees of partial short-range order compatible with the specified S_{2}. The mechanical performance including both elastic and brittle fracture behaviors of the generated hyperuniform materials is analyzed using a volume-compensated lattice-particle method. For the purpose of comparison, the corresponding nonhyperuniform materials with the same short-range order (i.e., with S_{2} constrained up to the same γ value) are also constructed and their mechanical performance is analyzed. Here we consider two specific S_{2} including the positive exponential decay function and the correlation function associated with an equilibrium hard-sphere system. For the constructed systems associated with these two specific functions, we find that although the hyperuniform materials are softer than their nonhyperuniform counterparts, the former generally possess a significantly higher brittle fracture strength than the latter. This superior mechanical behavior is attributed to the lower degree of stress concentration in the material resulting from the hyperuniform microstructure, which is crucial to crack initiation and propagation.

Characterizing pixel and point patterns with a hyperuniformity disorder length

We introduce the concept of a "hyperuniformity disorder length" h that controls the variance of volume fraction fluctuations for randomly placed windows of fixed size. In particular, fluctuations are determined by the average number of particles within a distance h from the boundary of the window. We first compute special expectations and bounds in d dimensions, and then illustrate the range of behavior of h versus window size L by analyzing several different types of simulated two-dimensional pixel patterns-where particle positions are stored as a binary digital image in which pixels have value zero if empty and one if they contain a particle. The first are random binomial patterns, where pixels are randomly flipped from zero to one with probability equal to area fraction. These have long-ranged density fluctuations, and simulations confirm the exact result h=L/2. Next we consider vacancy patterns, where a fraction f of particles on a lattice are randomly removed. These also display long-range density fluctuations, but with h=(L/2)(f/d) for small f, and h=L/2 for f→1. And finally, for a hyperuniform system with no long-range density fluctuations, we consider "Einstein patterns," where each particle is independently displaced from a lattice site by a Gaussian-distributed amount. For these, at large L,h approaches a constant equal to about half the root-mean-square displacement in each dimension. Then we turn to gray-scale pixel patterns that represent simulated arrangements of polydisperse particles, where the volume of a particle is encoded in the value of its central pixel. And we discuss the continuum limit of point patterns, where pixel size vanishes. In general, we thus propose to quantify particle configurations not just by the scaling of the density fluctuation spectrum but rather by the real-space spectrum of h(L) versus L. We call this approach "hyperuniformity disorder length spectroscopy".

Hyperuniformity disorder length spectroscopy for extended particles

The concept of a hyperuniformity disorder length h was recently introduced for analyzing volume fraction fluctuations for a set of measuring windows [Chieco et al., Phys. Rev. E 96, 032909 (2017).PLEEE81539-375510.1103/PhysRevE.96.032909]. This length permits a direct connection to the nature of disorder in the spatial configuration of the particles and provides a way to diagnose the degree of hyperuniformity in terms of the scaling of h and its value in comparison with established bounds. Here, this approach is generalized for extended particles, which are larger than the image resolution and can lie partially inside and partially outside the measuring windows. The starting point is an expression for the relative volume fraction variance in terms of four distinct volumes: that of the particle, the measuring window, the mean-squared overlap between particle and region, and the region over which particles have nonzero overlap with the measuring window. After establishing limiting behaviors for the relative variance, computational methods are developed for both continuum and pixelated particles. Exact results are presented for particles of special shape and for measuring windows of special shape, for which the equations are tractable. Comparison is made for other particle shapes, using simulated Poisson patterns. And the effects of polydispersity and image errors are discussed. For small measuring windows, both particle shape and spatial arrangement affect the form of the variance. For large regions, the variance scaling depends only on arrangement but particle shape sets the numerical proportionality. The combined understanding permit the measured variance to be translated to the spectrum of hyperuniformity lengths versus region size, as the quantifier of spatial arrangement. This program is demonstrated for a system of nonoverlapping particles at a series of increasing packing fractions as well as for an Einstein pattern of particles with several different extended shapes.

Large-scale structure of randomly jammed spheres

We numerically analyze the density field of three-dimensional randomly jammed packings of monodisperse soft frictionless spherical particles, paying special attention to fluctuations occurring at large length scales. We study in detail the two-point static structure factor at low wave vectors in Fourier space. We also analyze the nature of the density field in real space by studying the large-distance behavior of the two-point pair correlation function, of density fluctuations in subsystems of increasing sizes, and of the direct correlation function. We show that such real space analysis can be greatly improved by introducing a coarse-grained density field to disentangle genuine large-scale correlations from purely local effects. Our results confirm that both Fourier and real space signatures of vanishing density fluctuations at large scale are absent, indicating that randomly jammed packings are not hyperuniform. In addition, we establish that the pair correlation function displays a surprisingly complex structure at large distances, which is however not compatible with the long-range negative correlation of hyperuniform systems but fully compatible with an analytic form for the structure factor. This implies that the direct correlation function is short ranged, as we also demonstrate directly. Our results reveal that density fluctuations in jammed packings do not follow the behavior expected for random hyperuniform materials, but display instead a more complex behavior.

The Perfect Glass Paradigm: Disordered Hyperuniform Glasses Down to Absolute Zero

Rapid cooling of liquids below a certain temperature range can result in a transition to glassy states. The traditional understanding of glasses includes their thermodynamic metastability with respect to crystals. However, here we present specific examples of interactions that eliminate the possibilities of crystalline and quasicrystalline phases, while creating mechanically stable amorphous glasses down to absolute zero temperature. We show that this can be accomplished by introducing a new ideal state of matter called a "perfect glass". A perfect glass represents a soft-interaction analog of the maximally random jammed (MRJ) packings of hard particles. These latter states can be regarded as the epitome of a glass since they are out of equilibrium, maximally disordered, hyperuniform, mechanically rigid with infinite bulk and shear moduli, and can never crystallize due to configuration-space trapping. Our model perfect glass utilizes two-, three-, and four-body soft interactions while simultaneously retaining the salient attributes of the MRJ state. These models constitute a theoretical proof of concept for perfect glasses and broaden our fundamental understanding of glass physics. A novel feature of equilibrium systems of identical particles interacting with the perfect-glass potential at positive temperature is that they have a non-relativistic speed of sound that is infinite.

Disordered hyperuniform heterogeneous materials

Disordered hyperuniform many-body systems are distinguishable states of matter that lie between a crystal and liquid: they are like perfect crystals in the way they suppress large-scale density fluctuations and yet are like liquids or glasses in that they are statistically isotropic with no Bragg peaks. These systems play a vital role in a number of fundamental and applied problems: glass formation, jamming, rigidity, photonic and electronic band structure, localization of waves and excitations, self-organization, fluid dynamics, quantum systems, and pure mathematics. Much of what we know theoretically about disordered hyperuniform states of matter involves many-particle systems. In this paper, we derive new rigorous criteria that disordered hyperuniform two-phase heterogeneous materials must obey and explore their consequences. Two-phase heterogeneous media are ubiquitous; examples include composites and porous media, biological media, foams, polymer blends, granular media, cellular solids, and colloids. We begin by obtaining some results that apply to hyperuniform two-phase media in which one phase is a sphere packing in d-dimensional Euclidean space [Formula: see text]. Among other results, we rigorously establish the requirements for packings of spheres of different sizes to be 'multihyperuniform'. We then consider hyperuniformity for general two-phase media in [Formula: see text]. Here we apply realizability conditions for an autocovariance function and its associated spectral density of a two-phase medium, and then incorporate hyperuniformity as a constraint in order to derive new conditions. We show that some functional forms can immediately be eliminated from consideration and identify other forms that are allowable. Specific examples and counterexamples are described. Contact is made with well-known microstructural models (e.g. overlapping spheres and checkerboards) as well as irregular phase-separation and Turing-type patterns. We also ascertain a family of autocovariance functions (or spectral densities) that are realizable by disordered hyperuniform two-phase media in any space dimension, and present select explicit constructions of realizations. These studies provide insight into the nature of disordered hyperuniformity in the context of heterogeneous materials and have implications for the design of such novel amorphous materials.

Hyperuniformity and its generalizations

Disordered many-particle hyperuniform systems are exotic amorphous states of matter that lie between crystal and liquid: They are like perfect crystals in the way they suppress large-scale density fluctuations and yet are like liquids or glasses in that they are statistically isotropic with no Bragg peaks. These exotic states of matter play a vital role in a number of problems across the physical, mathematical as well as biological sciences and, because they are endowed with novel physical properties, have technological importance. Given the fundamental as well as practical importance of disordered hyperuniform systems elucidated thus far, it is natural to explore the generalizations of the hyperuniformity notion and its consequences. In this paper, we substantially broaden the hyperuniformity concept along four different directions. This includes generalizations to treat fluctuations in the interfacial area (one of the Minkowski functionals) in heterogeneous media and surface-area driven evolving microstructures, random scalar fields, divergence-free random vector fields, and statistically anisotropic many-particle systems and two-phase media. In all cases, the relevant mathematical underpinnings are formulated and illustrative calculations are provided. Interfacial-area fluctuations play a major role in characterizing the microstructure of two-phase systems (e.g., fluid-saturated porous media), physical properties that intimately depend on the geometry of the interface, and evolving two-phase microstructures that depend on interfacial energies (e.g., spinodal decomposition). In the instances of random vector fields and statistically anisotropic structures, we show that the standard definition of hyperuniformity must be generalized such that it accounts for the dependence of the relevant spectral functions on the direction in which the origin in Fourier space is approached (nonanalyticities at the origin). Using this analysis, we place some well-known energy spectra from the theory of isotropic turbulence in the context of this generalization of hyperuniformity. Among other results, we show that there exist many-particle ground-state configurations in which directional hyperuniformity imparts exotic anisotropic physical properties (e.g., elastic, optical, and acoustic characteristics) to these states of matter. Such tunability could have technological relevance for manipulating light and sound waves in ways heretofore not thought possible. We show that disordered many-particle systems that respond to external fields (e.g., magnetic and electric fields) are a natural class of materials to look for directional hyperuniformity. The generalizations of hyperuniformity introduced here provide theoreticians and experimentalists new avenues to understand a very broad range of phenomena across a variety of fields through the hyperuniformity "lens."

Characterization of maximally random jammed sphere packings. II. Correlation functions and density fluctuations

In the first paper of this series, we introduced Voronoi correlation functions to characterize the structure of maximally random jammed (MRJ) sphere packings across length scales. In the present paper, we determine a variety of different correlation functions that arise in rigorous expressions for the effective physical properties of MRJ sphere packings and compare them to the corresponding statistical descriptors for overlapping spheres and equilibrium hard-sphere systems. Such structural descriptors arise in rigorous bounds and formulas for effective transport properties, diffusion and reactions constants, elastic moduli, and electromagnetic characteristics. First, we calculate the two-point, surface-void, and surface-surface correlation functions, for which we derive explicit analytical formulas for finite hard-sphere packings. We show analytically how the contact Dirac delta function contribution to the pair correlation function g_{2}(r) for MRJ packings translates into distinct functional behaviors of these two-point correlation functions that do not arise in the other two models examined here. Then we show how the spectral density distinguishes the MRJ packings from the other disordered systems in that the spectral density vanishes in the limit of infinite wavelengths; i.e., these packings are hyperuniform, which means that density fluctuations on large length scales are anomalously suppressed. Moreover, for all model systems, we study and compute exclusion probabilities and pore size distributions, as well as local density fluctuations. We conjecture that for general disordered hard-sphere packings, a central limit theorem holds for the number of points within an spherical observation window. Our analysis links problems of interest in material science, chemistry, physics, and mathematics. In the third paper of this series, we will evaluate bounds and estimates of a host of different physical properties of the MRJ sphere packings that are based on the structural characteristics analyzed in this paper.

Critical slowing down and hyperuniformity on approach to jamming

Hyperuniformity characterizes a state of matter that is poised at a critical point at which density or volume-fraction fluctuations are anomalously suppressed at infinite wavelengths. Recently, much attention has been given to the link between strict jamming (mechanical rigidity) and (effective or exact) hyperuniformity in frictionless hard-particle packings. However, in doing so, one must necessarily study very large packings in order to access the long-ranged behavior and to ensure that the packings are truly jammed. We modify the rigorous linear programming method of Donev et al. [J. Comput. Phys. 197, 139 (2004)JCTPAH0021-999110.1016/j.jcp.2003.11.022] in order to test for jamming in putatively collectively and strictly jammed packings of hard disks in two dimensions. We show that this rigorous jamming test is superior to standard ways to ascertain jamming, including the so-called "pressure-leak" test. We find that various standard packing protocols struggle to reliably create packings that are jammed for even modest system sizes of N≈10^{3} bidisperse disks in two dimensions; importantly, these packings have a high reduced pressure that persists over extended amounts of time, meaning that they appear to be jammed by conventional tests, though rigorous jamming tests reveal that they are not. We present evidence that suggests that deviations from hyperuniformity in putative maximally random jammed (MRJ) packings can in part be explained by a shortcoming of the numerical protocols to generate exactly jammed configurations as a result of a type of "critical slowing down" as the packing's collective rearrangements in configuration space become locally confined by high-dimensional "bottlenecks" from which escape is a rare event. Additionally, various protocols are able to produce packings exhibiting hyperuniformity to different extents, but this is because certain protocols are better able to approach exactly jammed configurations. Nonetheless, while one should not generally expect exact hyperuniformity for disordered packings with rattlers, we find that when jamming is ensured, our packings are very nearly hyperuniform, and deviations from hyperuniformity correlate with an inability to ensure jamming, suggesting that strict jamming and hyperuniformity are indeed linked. This raises the possibility that the ideal MRJ packings have no rattlers. Our work provides the impetus for the development of packing algorithms that produce large disordered strictly jammed packings that are rattler free, which is an outstanding, challenging task.

Experimental study of the relationship between local particle-size distributions and local ordering in random close packing

We experimentally study the structural properties of a sediment of size distributed colloids. By determining each particle size using a size estimation algorithm, we are able to investigate the relationship between local environment and local ordering. Our results show that ordered environments of particles tend to generate where the local particle-size distribution is within 5%. In addition, we show that particles whose size is close to the average size have 12 coordinate neighbors, which matches the coordination number of the fcc and hcp crystals. On the other hand, bcc structures are observed around larger particles. Our results represent experiments to show a size dependence of the specific ordering in colloidal systems.

Nanoscale engineering of two-dimensional disordered hyperuniform block-copolymer assemblies

Disordered hyperuniform (DH) media have been recognized as a new state of disordered matter that broadens our vision of material engineering. Here, long-range correlated disordered two-dimensional patterns are fabricated by self-assembling of spherical diblock-copolymer (BCP) micelles. Control of the self-assembling parameters leads to the formation of DH patterns of micelles that can host nanoscale material inclusions, therefore providing an effective strategy for fabricating multimaterial DH structures at molecular scale. Centroidal patterns are accurately determined by virtue of BCP micelles loaded with metal nanoparticles. Our analysis reveals the signature of nearly ideal DH BCP assemblies in the local density fluctuation and a dominant linear scaling in the local number fluctuation.

Ground states of stealthy hyperuniform potentials: I. Entropically favored configurations

Systems of particles interacting with "stealthy" pair potentials have been shown to possess infinitely degenerate disordered hyperuniform classical ground states with novel physical properties. Previous attempts to sample the infinitely degenerate ground states used energy minimization techniques, introducing algorithmic dependence that is artificial in nature. Recently, an ensemble theory of stealthy hyperuniform ground states was formulated to predict the structure and thermodynamics that was shown to be in excellent agreement with corresponding computer simulation results in the canonical ensemble (in the zero-temperature limit). In this paper, we provide details and justifications of the simulation procedure, which involves performing molecular dynamics simulations at sufficiently low temperatures and minimizing the energy of the snapshots for both the high-density disordered regime, where the theory applies, as well as lower densities. We also use numerical simulations to extend our study to the lower-density regime. We report results for the pair correlation functions, structure factors, and Voronoi cell statistics. In the high-density regime, we verify the theoretical ansatz that stealthy disordered ground states behave like "pseudo" disordered equilibrium hard-sphere systems in Fourier space. The pair statistics obey certain exact integral conditions with very high accuracy. These results show that as the density decreases from the high-density limit, the disordered ground states in the canonical ensemble are characterized by an increasing degree of short-range order and eventually the system undergoes a phase transition to crystalline ground states. In the crystalline regime (low densities), there exist aperiodic structures that are part of the ground-state manifold but yet are not entropically favored. We also provide numerical evidence suggesting that different forms of stealthy pair potentials produce the same ground-state ensemble in the zero-temperature limit. Our techniques may be applied to sample the zero-temperature limit of the canonical ensemble of other potentials with highly degenerate ground states.

Thermal fluctuations, mechanical response, and hyperuniformity in jammed solids

Jamming is a geometric phase transition occurring in dense particle systems in the absence of temperature. We use computer simulations to analyze the effect of thermal fluctuations on several signatures of the transition. We show that scaling laws for bulk and shear moduli only become relevant when thermal fluctuations are extremely small, and propose their relative ratio as a quantitative signature of jamming criticality. Despite the nonequilibrium nature of the transition, we find that thermally induced fluctuations and mechanical responses obey equilibrium fluctuation-dissipation relations near jamming, provided the appropriate fluctuating component of the particle displacements is analyzed. This shows that mechanical moduli can be directly measured from particle positions in mechanically unperturbed packings, and suggests that the definition of a "nonequilibrium index" is unnecessary for amorphous materials. We find that fluctuations of particle displacements are spatially correlated, and define a transverse and a longitudinal correlation length scale which both diverge as the jamming transition is approached. We analyze the frozen component of density fluctuations and find that it displays signatures of nearly hyperuniform behavior at large length scales. This demonstrates that hyperuniformity in jammed packings is unrelated to a vanishing compressibility and explains why it appears remarkably robust against temperature and density variations. Differently from jamming criticality, obstacles preventing the observation of hyperuniformity in colloidal systems do not originate from thermal fluctuations.

Diagnosing hyperuniformity in two-dimensional, disordered, jammed packings of soft spheres

Hyperuniformity characterizes a state of matter for which (scaled) density fluctuations diminish towards zero at the largest length scales. However, the task of determining whether or not an image of an experimental system is hyperuniform is experimentally challenging due to finite-resolution, noise, and sample-size effects that influence characterization measurements. Here we explore these issues, employing video optical microscopy to study hyperuniformity phenomena in disordered two-dimensional jammed packings of soft spheres. Using a combination of experiment and simulation we characterize the possible adverse effects of particle polydispersity, image noise, and finite-size effects on the assignment of hyperuniformity, and we develop a methodology that permits improved diagnosis of hyperuniformity from real-space measurements. The key to this improvement is a simple packing reconstruction algorithm that incorporates particle polydispersity to minimize the free volume. In addition, simulations show that hyperuniformity in finite-sized samples can be ascertained more accurately in direct space than in reciprocal space. Finally, our experimental colloidal packings of soft polymeric spheres are shown to be effectively hyperuniform.

Disordered strictly jammed binary sphere packings attain an anomalously large range of densities

Previous attempts to simulate disordered binary sphere packings have been limited in producing mechanically stable, isostatic packings across a broad spectrum of packing fractions. Here we report that disordered strictly jammed binary packings (packings that remain mechanically stable under general shear deformations and compressions) can be produced with an anomalously large range of average packing fractions 0.634≤φ≤0.829 for small to large sphere radius ratios α restricted to α≥0.100. Surprisingly, this range of average packing fractions is obtained for packings containing a subset of spheres (called the backbone) that are exactly strictly jammed, exactly isostatic, and also generated from random initial conditions. Additionally, the average packing fractions of these packings at certain α and small sphere relative number concentrations x approach those of the corresponding densest known ordered packings. These findings suggest for entropic reasons that these high-density disordered packings should be good glass formers and that they may be easy to prepare experimentally. We also identify an unusual feature of the packing fraction of jammed backbones (packings with rattlers excluded). The backbone packing fraction is about 0.624 over the majority of the α-x plane, even when large numbers of small spheres are present in the backbone. Over the (relatively small) area of the α-x plane where the backbone is not roughly constant, we find that backbone packing fractions range from about 0.606 to 0.829, with the volume of rattler spheres comprising between 1.6% and 26.9% of total sphere volume. To generate isostatic strictly jammed packings, we use an implementation of the Torquato-Jiao sequential linear programming algorithm [Phys. Rev. E 82, 061302 (2010)], which is an efficient producer of inherent structures (mechanically stable configurations at the local maxima in the density landscape). The identification and explicit construction of binary packings with such high packing fractions could have important practical implications for granular composites where density is critical both to material properties and fabrication cost, including for solid propellants, concrete, and ceramics. The densities and structures of jammed binary packings at various α and x are also relevant to the formation of a glass phase in multicomponent metallic systems.

The physics of the colloidal glass transition

As one increases the concentration of a colloidal suspension, the system exhibits a dramatic increase in viscosity. Beyond a certain concentration, the system is said to be a colloidal glass; structurally, the system resembles a liquid, yet motions within the suspension are slow enough that it can be considered essentially frozen. For several decades, colloids have served as a valuable model system for understanding the glass transition in molecular systems. The spatial and temporal scales involved allow these systems to be studied by a wide variety of experimental techniques. The focus of this review is the current state of understanding of the colloidal glass transition, with an emphasis on experimental observations. A brief introduction is given to important experimental techniques used to study the glass transition in colloids. We describe features of colloidal systems near and in glassy states, including increases in viscosity and relaxation times, dynamical heterogeneity and ageing, among others. We also compare and contrast the glass transition in colloids to that in molecular liquids. Other glassy systems are briefly discussed, as well as recently developed synthesis techniques that will keep these systems rich with interesting physics for years to come.

Maximally random jammed packings of Platonic solids: hyperuniform long-range correlations and isostaticity

We generate maximally random jammed (MRJ) packings of the four nontiling Platonic solids (tetrahedra, octahedra, dodecahedra, and icosahedra) using the adaptive-shrinking-cell method [S. Torquato and Y. Jiao, Phys. Rev. E 80, 041104 (2009)]. Such packings can be viewed as prototypical glasses in that they are maximally disordered while simultaneously being mechanically rigid. The MRJ packing fractions for tetrahedra, octahedra, dodecahedra, and icosahedra are, respectively, 0.763±0.005, 0.697±0.005, 0.716±0.002, and 0.707±0.002. We find that as the number of facets of the particles increases, the translational order in the packings increases while the orientational order decreases. Moreover, we show that the MRJ packings are hyperuniform (i.e., their infinite-wavelength local-number-density fluctuations vanish) and possess quasi-long-range pair correlations that decay asymptotically with scaling r(-4). This provides further evidence that hyperuniform quasi-long-range correlations are a universal feature of MRJ packings of frictionless particles of general shape. However, unlike MRJ packings of ellipsoids, superballs, and superellipsoids, which are hypostatic, MRJ packings of the nontiling Platonic solids are isostatic. We provide a rationale for the organizing principle that the MRJ packing fractions for nonspherical particles with sufficiently small asphericities exceed the corresponding value for spheres (∼0.64). We also discuss how the shape and symmetry of a polyhedron particle affects its MRJ packing fraction.