Static structural signatures of nearly jammed disordered and ordered hard-sphere packings: Direct correlation function

Medium-Range Structural Order as the Driver of Activated Dynamics and Complexity Reduction in Glass-Forming Liquids

We analyze in depth the Elastically Collective Nonlinear Langevin Equation theory of activated dynamics in metastable liquids to establish that the predicted inter-relationships between the alpha relaxation time, local cage and collective elastic barriers, dynamic localization length, and shear modulus are causally related within the theory to the medium range order (MRO) static correlation length. The latter grows exponentially with density for metastable hard sphere fluids and as a nonuniversal inverse power law with temperature for supercooled liquids under isobaric conditions. The physical origin of predicted connections between the alpha time and other metrics of cage order and the thermodynamic inverse dimensionless compressibility is fully established. It is discovered that although kinetic constraints from the real space first coordination shell are important for the alpha time, they are of secondary importance compared to the consequences of the more universal MRO correlations in both the modestly and deeply metastable regimes. This understanding sheds new light on the theoretical basis for, and prior successes of, the predictive mapping of chemically complex thermal liquids to effective hard sphere fluids based on matching their dimensionless compressibilities, a scheme we call "complexity reduction". In essence, the latter is equivalent to the physical requirement that the thermal liquid MRO correlation equals that of its effective hard sphere analog. The mapping alone is shown to provide a remarkable level of quantitative predictive power for the glass transition temperature Tg of 21 molecular and polymer liquids. Predictions for the chemically specific absolute magnitude and growth with cooling of the MRO correlation length are obtained and lie in the window of 2-6 nm at Tg. Dynamic heterogeneity, elastic facilitation, and beyond pair structure issues are briefly discussed. Future opportunities to theoretically analyze the equilibrated deep glass regime are outlined.

Theoretical analysis of the structure, thermodynamics, and shear elasticity of deeply metastable hard sphere fluids

The structure, thermodynamics, and slow activated dynamics of the equilibrated metastable regime of glass-forming fluids remain a poorly understood problem of high theoretical and experimental interest. We apply a highly accurate microscopic equilibrium liquid state integral equation theory, in conjunction with naïve mode coupling theory of particle localization, to study in a unified manner the structural correlations, thermodynamic properties, and dynamic elastic shear modulus in deeply metastable hard sphere fluids. Distinctive behaviors are predicted including divergent inverse critical power laws for the contact value of the pair correlation function, pressure, and inverse dimensionless compressibility, and a splitting of the second peak and large suppression of interstitial configurations of the pair correlation function. The dynamic elastic modulus is predicted to exhibit two distinct exponential growth regimes with packing fraction that have strongly different slopes. These thermodynamic, structural, and elastic modulus results are consistent with simulations and experiments. Perhaps most unexpectedly, connections between the amplitude of long wavelength density fluctuations, dimensionless compressibility, local structure, and the dynamic elastic shear modulus have been theoretically elucidated. These connections are more broadly relevant to understanding the slow activated relaxation and mechanical response of colloidal suspensions in the ultradense metastable region and deeply supercooled thermal liquids in equilibrium.

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.

Resolving entropy contributions in nonequilibrium transitions

We derive a functional for the entropy contributed by any microscopic degrees of freedom as arising from their measurable pair correlations. Applicable both in and out of equilibrium, this functional yields the maximum entropy which a system can have given a certain correlation function. When applied to different correlations, the method allows us to identify the degrees of freedom governing a certain physical regime, thus capturing and characterizing dynamic transitions. The formalism applies also to systems whose translational invariance is broken by external forces and whose number of particles may vary. We apply it to experimental results for jammed bidisperse emulsions, capturing the crossover of this nonequilibrium system from crystalline to disordered hyperuniform structures as a function of mixture composition. We discover that the cross-correlations between the positions and sizes of droplets in the emulsion play the central role in the formation of the disordered hyperuniform states. We discuss implications of the approach for entropy estimation out of equilibrium and for characterizing transitions in disordered systems.

Non-hyperuniform metastable states around a disordered hyperuniform state of densely packed spheres: stochastic density functional theory at strong coupling

The disordered and hyperuniform structures of densely packed spheres near and at jamming are characterized by vanishing of long-wavelength density fluctuations, or equivalently by long-range power-law decay of the direct correlation function (DCF). We focus on previous simulation results that exhibit the degradation of hyperuniformity in jammed structures while maintaining the long-range nature of the DCF to a certain length scale. Here we demonstrate that the field-theoretic formulation of stochastic density functional theory is relevant to explore the degradation mechanism. The strong-coupling expansion method of stochastic density functional theory is developed to obtain the metastable chemical potential considering the intermittent fluctuations in dense packings. The metastable chemical potential yields the analytical form of the metastable DCF that has a short-range cutoff inside the sphere while retaining the long-range power-law behavior. It is confirmed that the metastable DCF provides the zero-wavevector limit of the structure factor in quantitative agreement with the previous simulation results of degraded hyperuniformity. We can also predict the emergence of soft modes localized at the particle scale by plugging this metastable DCF into the linearized Dean-Kawasaki equation, a stochastic density functional equation.

Cloaking the underlying long-range order of randomly perturbed lattices

Random, uncorrelated displacements of particles on a lattice preserve the hyperuniformity of the original lattice, that is, normalized density fluctuations vanish in the limit of infinite wavelengths. In addition to a diffuse contribution, the scattering intensity from the the resulting point pattern typically inherits the Bragg peaks (long-range order) of the original lattice. Here we demonstrate how these Bragg peaks can be hidden in the effective diffraction pattern of independent and identically distributed perturbations. All Bragg peaks vanish if and only if the sum of all probability densities of the positions of the shifted lattice points is a constant at all positions. The underlying long-range order is then "cloaked" in the sense that it cannot be reconstructed from the pair correlation function alone. On the one hand, density fluctuations increase monotonically with the strength of perturbations a, as measured by the hyperuniformity order metric Λ[over ¯]. On the other hand, the disappearance and reemergence of long-range order, depending on whether the system is cloaked as the perturbation strength increases, is manifestly captured by the τ order metric. Therefore, while the perturbation strength a may seem to be a natural choice for an order metric of perturbed lattices, the τ order metric is a superior choice. It is noteworthy that cloaked perturbed lattices allow one to easily simulate very large samples (with at least 10^{6} particles) of disordered hyperuniform point patterns without Bragg peaks.

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.

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.

Can exotic disordered "stealthy" particle configurations tolerate arbitrarily large holes?

The probability of finding a spherical cavity or "hole" of arbitrarily large size in typical disordered many-particle systems in the infinite-system-size limit (e.g., equilibrium liquid states) is non-zero. Such "hole" statistics are intimately linked to the thermodynamic and nonequilibrium physical properties of the system. Disordered "stealthy" many-particle configurations in d-dimensional Euclidean space [Doublestruck R]^d are exotic amorphous states of matter that lie between a liquid and crystal that prohibit single-scattering events for a range of wave vectors and possess no Bragg peaks [Torquato et al., Phys. Rev. X, 2015, 5, 021020]. In this paper, we provide strong numerical evidence that disordered stealthy configurations across the first three space dimensions cannot tolerate arbitrarily large holes in the infinite-system-size limit, i.e., the hole probability has compact support. This structural "rigidity" property apparently endows disordered stealthy systems with novel thermodynamic and physical properties, including desirable band-gap, optical and transport characteristics. We also determine the maximum hole size that any stealthy system can possess across the first three space dimensions.

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.