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

Structural properties of hyperuniform Voronoi networks

Disordered hyperuniform many-particle systems are recently discovered exotic states of matter, characterized by the complete suppression of normalized infinite-wavelength density fluctuations, as in perfect crystals, while lacking conventional long-range order, as in liquids and glasses. In this work, we begin a program to quantify the structural properties of nonhyperuniform and hyperuniform networks. In particular, large two-dimensional (2D) Voronoi networks (graphs) containing approximately 10,000 nodes are created from a variety of different point configurations, including the antihyperuniform hyperplane intersection process (HIP), nonhyperuniform Poisson process, nonhyperuniform random sequential addition (RSA) saturated packing, and both non-stealthy and stealthy hyperuniform point processes. We carry out an extensive study of the Voronoi-cell area distribution of each of the networks by determining multiple metrics that characterize the distribution, including their average areas and corresponding variances as well as higher-order cumulants (i.e., skewness γ_{1} and excess kurtosis γ_{2}). We show that the HIP distribution is far from Gaussian, as evidenced by a high skewness (γ_{1}=3.16) and large positive excess kurtosis (γ_{2}=16.2). The Poisson (with γ_{1}=1.07 and γ_{2}=1.79) and non-stealthy hyperuniform (with γ_{1}=0.257 and γ_{2}=0.0217) distributions are Gaussian-like distributions, since they exhibit a small but positive skewness and excess kurtosis. The RSA (with γ_{1}=0.450 and γ_{2}=-0.0384) and the highest stealthy hyperuniform distributions (with γ_{1}=0.0272 and γ_{2}=-0.0626) are also non-Gaussian because of their low skewness and negative excess kurtosis, which is diametrically opposite of the non-Gaussian behavior of the HIP. The fact that the cell-area distributions of large, finite-sized RSA and stealthy hyperuniform networks (e.g., with N≈10,000 nodes) are narrower, have larger peaks, and smaller tails than a Gaussian distribution implies that in the thermodynamic limit the distributions should exhibit compact support, consistent with previous theoretical considerations. Moreover, we compute the Voronoi-area correlation functions C_{00}(r) for the networks, which describe the correlations between the area of two Voronoi cells separated by a given distance r [M. A. Klatt and S. Torquato, Phys. Rev. E 90, 052120 (2014)1539-375510.1103/PhysRevE.90.052120]. We show that the correlation functions C_{00}(r) qualitatively distinguish the antihyperuniform, nonhyperuniform, and hyperuniform Voronoi networks considered here. Specifically, the antihyperuniform HIP networks possess a slowly decaying C_{00}(r) with large positive values, indicating large fluctuations of Voronoi cell areas across scales. While the nonhyperuniform Poisson and RSA network possess positive and fast decaying C_{00}(r), we find strong anticorrelations in C_{00}(r) (i.e., negative values) for the hyperuniform networks. The latter indicates that the large-scale area fluctuations are suppressed by accompanying large Voronoi cells with small cells (and vice versa) in the systems in order to achieve hyperuniformity. In summary, we have shown that cell-area distributions and pair correlation functions of Voronoi networks enable one to distinguish quantitatively antihyperuniform, standard nonhyperuniform, and hyperuniform networks from one another.

Three-dimensional construction of hyperuniform, nonhyperuniform, and antihyperuniform disordered heterogeneous materials and their transport properties via spectral density functions

Rigorous theories connecting physical properties of a heterogeneous material to its microstructure offer a promising avenue to guide the computational material design and optimization. The spectral density function χ[over ̃]_{_{V}}(k), which can be obtained experimentally from scattering data, enables accurate determination of various transport and wave propagation characteristics, including the time-dependent diffusion spreadability S(t) and effective dynamic dielectric constant ε_{e} for electromagnetic wave propagation. Moreover, χ[over ̃]_{_{V}}(k) determines rigorous upper bounds on the fluid permeability K. Given the importance of χ[over ̃]_{_{V}}(k), we present here an efficient Fourier-space based computational framework to construct three-dimensional (3D) statistically isotropic two-phase heterogeneous materials corresponding to targeted spectral density functions. In particular, we employ a variety of analytical functional forms for χ[over ̃]_{_{V}}(k) that satisfy all known necessary conditions to construct disordered stealthy hyperuniform, standard hyperuniform, nonhyperuniform, and antihyperuniform two-phase heterogeneous material systems at varying phase volume fractions. We show that by tuning the correlations in the system across length scales via the targeted functions, one can generate a rich spectrum of distinct structures within each of the above classes of materials. Importantly, we present the first realization of antihyperuniform two-phase heterogeneous materials in 3D, which are characterized by autocovariance function χ_{_{V}}(r) with a power-law tail, resulting in microstructures that contain clusters of dramatically different sizes and morphologies. We also determine the diffusion spreadability S(t) and estimate the fluid permeability K associated with all of the constructed materials directly from the corresponding spectral densities. Although it is well established that the long-time asymptotic scaling behavior of S(t) only depends on the functional form of χ[over ̃]_{_{V}}(k), with the stealthy hyperuniform and antihyperuniform media, respectively, achieving the most and least efficient transport, we find that varying the length-scale parameter within each class of χ[over ̃]_{_{V}}(k) functions can also lead to orders of magnitude variation of S(t) at intermediate and long time scales. Moreover, we find that increasing the solid volume fraction ϕ_{1} and correlation length a in the constructed media generally leads to a decrease in the dimensionless fluid permeability K/a^{2}, while the antihyperuniform media possess the largest K/a^{2} among the four classes of materials with the same ϕ_{1} and a. These results indicate the feasibility of employing parameterized spectral densities for designing composites with targeted transport properties.

Generating Ultradense Jammed Ellipse Packings Using Biased SWAP

Using a Lubachevsky-Stillinger-like growth algorithm combined with biased SWAP Monte Carlo and transient degrees of freedom, we generate ultradense disordered jammed ellipse packings. For all aspect ratios α, these packings exhibit significantly smaller intermediate-wavelength density fluctuations and greater local nematic order than their less-dense counterparts. The densest packings are disordered despite having packing fractions ϕJ(α) that are within less than 0.5% of that of the monodisperse-ellipse crystal [ϕxtal = π/(2√3) ≃ 0.9069] over the range 1.2 ≲ α ≲ 1.45 and coordination numbers ZJ(α) that are within less than 0.5% of isostaticity [Ziso = 6] over the range 1.2 ≲ α ≲ 2.6. Lower-α packings are strongly fractionated and consist of polycrystals of intermediate-size particles, with the largest and smallest particles isolated at the grain boundaries. Higher-α packings are also fractionated, but in a qualitatively different fashion; they are composed─of increasingly large locally nematic domains─reminiscent of liquid glasses.

Microstructural characterization of DEM-based random packings of monodisperse and polydisperse non-convex particles

Understanding of hard particles in morphologies and sizes on microstructures of particle random packings is of significance to evaluate physical and mechanical properties of many discrete media, such as granular materials, colloids, porous ceramics, active cells, and concrete. The majority of previous lines of research mainly dedicated microstructure analysis of convex particles, such as spheres, ellipsoids, spherocylinders, cylinders, and convex-polyhedra, whereas little is known about non-convex particles that are more close to practical discrete objects in nature. In this study, the non-convex morphology of a three-dimensional particle is devised by using a mathematical-controllable parameterized method, which contains two construction modes, namely, the uniformly distributed contraction centers and the randomly distributed contraction centers. Accordingly, three shape parameters are conceived to regulate the particle geometrical morphology from a perfect sphere to arbitrary non-convexities. Random packing models of hard non-convex particles with mono-/poly-dispersity in sizes are then established using the discrete element modeling Diverse microstructural indicators are utilized to characterize configurations of non-convex particle random packings. The compactness of non-convex particles in packings is characterized by the random close packing fraction fd and the corresponding average coordination number Z. In addition, four statistical descriptors, encompassing the radial distribution function g(r), two-point probability function S2(i)(r), lineal-path function L(i)(r), and cumulative pore size distribution function F(δ), are exploited to demonstrate the high-order microstructure information of non-convex particle random packings. The results demonstrate that the particle shape and size distribution have significant effects on Z and fd; the construction mode of the randomly distributed contraction centers can yield higher fd than that of the uniformly distributed contraction centers, in which the upper limit of fd approaches to 0.632 for monodisperse sphere packings. Moreover, non-convex particles of sizes following the famous Fuller distribution of the power-law distribution of the exponent q = 2.5, have the highest fd (≈0.761) with respect to other q. In contrast, the particle shapes have an almost negligible effect on the four statistical descriptors, but they are remarkably sensitive to particle packing fraction fp and size distribution. The results can provide sound guidance for custom-design of granular media by tailoring specific microstructures of particles.

Fast generation of spectrally shaped disorder

Media with correlated disorder display unexpected transport properties, but it is still a challenge to design structures with desired spectral features at scale. In this work, we introduce an optimal formulation of this inverse problem by means of the nonuniform fast Fourier transform, thus arriving at an algorithm capable of generating systems with arbitrary spectral properties, with a computational cost that scales O(NlogN) with system size. The method is extended to accommodate arbitrary real-space interactions, such as short-range repulsion, to simultaneously control short- and long-range correlations. We thus generate the largest-ever stealthy hyperuniform configurations in 2d (N=10^{9}) and 3d (N>10^{7}) and demonstrate the flexibility of the approach by generating structures with designed spectral features at scale. By an Ewald sphere construction we link the spectral and optical properties at the single-scattering level and show that stealthy hyperuniform structures generically display transmission gaps, providing a concrete example of fine-tuning of a physical property. We also show that large 3d power-law hyperuniformity in particle packings leads to single-scattering properties nearly identical to those of simple hard spheres. Finally, we demonstrate generalizations of the approach to impose features in either continuous or discrete real space, using constraints in either continuous or discrete reciprocal space. In particular, enforcing large spectral power at peaks with the right symmetry leads to the nondeterministic generation of quasicrystalline structures in 2d and 3d. This technique should become an essential tool to embed, and understand the role of, long-range correlations in disordered metamaterials.

Random close packing of semi-flexible polymers in two dimensions: Emergence of local and global order

Through extensive Monte Carlo simulations, we systematically study the effect of chain stiffness on the packing ability of linear polymers composed of hard spheres in extremely confined monolayers, corresponding effectively to 2D films. First, we explore the limit of random close packing as a function of the equilibrium bending angle and then quantify the local and global order by the degree of crystallinity and the nematic or tetratic orientational order parameter, respectively. A multi-scale wealth of structural behavior is observed, which is inherently absent in the case of athermal individual monomers and is surprisingly richer than its 3D counterpart under bulk conditions. As a general trend, an isotropic to nematic transition is observed at sufficiently high surface coverages, which is followed by the establishment of the tetratic state, which in turn marks the onset of the random close packing. For chains with right-angle bonds, the incompatibility of the imposed bending angle with the neighbor geometry of the triangular crystal leads to a singular intra- and inter-polymer tiling pattern made of squares and triangles with optimal local filling at high surface concentrations. The present study could serve as a first step toward the design of hard colloidal polymers with a tunable structural behavior for 2D applications.

Optimal three-dimensional particle shapes for maximally dense saturated packing

Saturated packing is a random packing state of particles widely applied in investigating the physicochemical properties of granular materials. Optimizing particle shape to maximize packing density is a crucial challenge in saturated packing research. The known optimal three-dimensional shape is an ellipsoid with a saturated packing density of 0.437 72(51). In this work, we generate saturated packings of three-dimensional asymmetric shapes, including spherocylinders, cones, and tetrahedra, via the random sequential adsorption algorithm and investigate their packing properties. Results show that the optimal shape of asymmetric spherocylinders gives the maximum density of 0.4338(9), while cones achieve a higher value of 0.4398(10). Interestingly, tetrahedra exhibit two distinct optimal shapes with significantly high densities of 0.4789(19) and 0.4769(18), which surpass all previous results in saturated packing. The study of adsorption kinetics reveals that the two optimal shapes of tetrahedra demonstrate notably higher degrees of freedom and faster growth rates of the particle number. The analysis of packing structures via the density pair-correlation function shows that the two optimal shapes of tetrahedra possess faster transitions from local to global packing densities.

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.

Machine learning approaches for the optimization of packing densities in granular matter

The fundamental question of how densely granular matter can pack and how this density depends on the shape of the constituent particles has been a longstanding scientific problem. Previous work has mainly focused on empirical approaches based on simulations or mean-field theory to investigate the effect of shape variation on the resulting packing densities, focusing on a small set of pre-defined shapes like dimers, ellipsoids, and spherocylinders. Here we discuss how machine learning methods can support the search for optimally dense packing shapes in a high-dimensional shape space. We apply dimensional reduction and regression techniques based on random forests and neural networks to find novel dense packing shapes by numerical optimization. Moreover, an investigation of the regression function in the dimensionally reduced shape representation allows us to identify directions in the packing density landscape that lead to a strongly non-monotonic variation of the packing density. The predictions obtained by machine learning are compared with packing simulations. Our approach can be more widely applied to optimize the properties of granular matter by varying the shape of its constituent particles.

Dense Disordered Jammed Packings of Hard Spherocylinders with a Low Aspect Ratio: A Characterization of Their Structure

This work numerically investigates dense disordered (maximally random) jammed packings of hard spherocylinders of cylinder length L and diameter D by focusing on L/D ∈ [0,2]. It is within this interval that one expects that the packing fraction of these dense disordered jammed packings ϕMRJ hsc attains a maximum. This work confirms the form of the graph ϕMRJ hsc versus L/D: here, comparably to certain previous investigations, it is found that the maximal ϕMRJ hsc = 0.721 ± 0.001 occurs at L/D = 0.45 ± 0.05. Furthermore, this work meticulously characterizes the structure of these dense disordered jammed packings via the special pair-correlation function of the interparticle distance scaled by the contact distance and the ensuing analysis of the statistics of the hard spherocylinders in contact: here, distinctly from all previous investigations, it is found that the dense disordered jammed packings of hard spherocylinders with 0.45 ≲ L/D ≤ 2 are isostatic.

Optimal shapes of disk assembly in saturated random packings

Particle morphology is one of the most significant factors influencing the packing structures of granular materials. With certain targeted properties or optimization criteria, inverse packing problems have drawn extensive attention in terms of their adaptability to many material design tasks. An important question hard to answer is which particle shape, especially within given shape families, forms the densest (loosest) random packing? In this paper, we address this issue for the disk assembly model in two dimensions with an infinite variety of shapes, which are simulated in the random sequential adsorption process to suppress crystallization. Via a unique shape representation method, particle shapes are transformed into genotype sequences in the continuous shape space where we utilize the genetic algorithm as an efficient shape optimizer. Specifically, we consider three representative species of disk assembly, i.e., congruent tangent disks, incongruent tangent disks, and congruent overlapping disks, and carry out shape optimization on their packing densities in the saturated random state. We numerically search optimal shapes in the three species with a variable number of constituent disks which yield the maximal and minimal packing densities. We obtain an isosceles circulo-triangle and an unclosed ring for the maximal and minimal packing density in saturated random packings, respectively. The perfect sno-cone and isosceles circulo-triangle are also specifically investigated which give remarkably high packing densities of around 0.6, much denser than those of ellipses. This study is beneficial for guiding the design of particle shapes as well as the inverse design of granular materials.

Mixing properties of bi-disperse ellipsoid assemblies: mean-field behaviour in a granular matter experiment

The structure and spatial statistical properties of amorphous ellipsoid assemblies have profound scientific and industrial significance in many systems, from cell assays to granular materials. This paper uses a fundamental theoretical relationship for mixture distributions to explain the observations of an extensive X-ray computed tomography study of granular ellipsoidal packings. We study a size-bi-disperse mixture of two types of ellipsoids of revolutions that have the same aspect ratio of α ≈ 0.57 and differ in size, by about 10% in linear dimension, and compare these to mono-disperse systems of ellipsoids with the same aspect ratio. Jammed configurations with a range of packing densities are achieved by employing different tapping protocols. We numerically interrogate the final packing configurations by analyses of the local packing fraction distributions calculated from the Voronoi diagrams. Our main finding is that the bi-disperse ellipsoidal packings studied here can be interpreted as a mixture of two uncorrelated mono-disperse packings, insensitive to the compaction protocol. Our results are consolidated by showing that the local packing fraction shows no correlation beyond their first shell of neighbours in the binary mixtures. We propose a model of uncorrelated binary mixture distribution that describes the observed experimental data with high accuracy. This analysis framework will enable future studies to test whether the observed mean-field behaviour is specific to the particular granular system or the specific parameter values studied here or if it is observed more broadly in other bi-disperse non-spherical particle systems.

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.

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.

Colloidal Shear-Thickening Fluids Using Variable Functional Star-Shaped Particles: A Molecular Dynamics Study

Complex colloidal fluids, depending on constituent shapes and packing fractions, may have a wide range of shear-thinning and/or shear-thickening behaviors. An interesting way to transition between different types of such behavior is by infusing complex functional particles that can be manufactured using modern techniques such as 3D printing. In this paper, we perform 2D molecular dynamics simulations of such fluids with infused star-shaped functional particles, with a variable leg length and number of legs, as they are infused in a non-interacting fluid. We vary the packing fraction (ϕ) of the system, and for each different system, we apply shear at various strain rates, turning the fluid into a shear-thickened fluid and then, in jammed state, rising the apparent viscosity of the fluid and incipient stresses. We demonstrate the dependence of viscosity on the functional particles’ packing fraction and we show the role of shape and design dependence of the functional particles towards the transition to a shear-thickening fluid.

Structural analysis of disordered dimer packings

Jammed disordered packings of non-spherical particles show significant variation in the packing density as a function of particle shape for a given packing protocol. Rotationally symmetric elongated shapes such as ellipsoids, spherocylinders, and dimers, e.g., pack significantly denser than spheres over a narrow range of aspect ratios, exhibiting a characteristic peak at aspect ratios of α_max ≈ 1.4-1.5. However, the structural features that underlie this non-monotonic behaviour in the packing density are unknown. Here, we study disordered packings of frictionless dimers in three dimensions generated by a gravitational pouring protocol in LAMMPS. Focusing on the characteristics of contacts as well as orientational and translational order metrics, we identify a number of structural features that accompany the formation of maximally dense packings as the dimer aspect ratio α is varied from the spherical limit. Our results highlight that dimer packings undergo significant structural changes as α increases up to α_max manifest in the reorganisation of the contact configurations between neighbouring dimers, increasing nematic order, and decreasing local translational order. Remarkably, for α > α_max our metrics remain largely unchanged, indicating that the peak in the packing density is related to the interplay of structural rearrangements for α < α_max and subsequent excluded volume effects with unchanged structure for α > α_max.

Structural characterization of many-particle systems on approach to hyperuniform states

The study of hyperuniform states of matter is an emerging multidisciplinary field, impinging on topics in the physical sciences, mathematics, and biology. The focus of this work is the exploration of quantitative descriptors that herald when a many-particle system in d-dimensional Euclidean space R^{d} approaches a hyperuniform state as a function of the relevant control parameter. We establish quantitative criteria to ascertain the extent of hyperuniform and nonhyperuniform distance-scaling regimes as well as the crossover point between them in terms of the "volume" coefficient A and "surface-area" coefficient B associated with the local number variance σ^{2}(R) for a spherical window of radius R. The larger the ratio B/A, the larger the hyperuniform scaling regime, which becomes of infinite extent in the limit B/A→∞. To complement the known direct-space representation of the coefficient B in terms of the total correlation function h(r), we derive its corresponding Fourier representation in terms of the structure factor S(k), which is especially useful when scattering information is available experimentally or theoretically. We also demonstrate that the free-volume theory of the pressure of equilibrium packings of identical hard spheres that approach a strictly jammed state either along the stable crystal or metastable disordered branch dictates that such end states be exactly hyperuniform. Using the ratio B/A, as well as other diagnostic measures of hyperuniformity, including the hyperuniformity index H and the direct-correlation function length scale ξ_{c}, we study three different exactly solvable models as a function of the relevant control parameter, either density or temperature, with end states that are perfectly hyperuniform. Specifically, we analyze equilibrium systems of hard rods and "sticky" hard-sphere systems in arbitrary space dimension d as a function of density. We also examine low-temperature excited states of many-particle systems interacting with "stealthy" long-ranged pair interactions as the temperature tends to zero, where the ground states are disordered, hyperuniform, and infinitely degenerate. We demonstrate that our various diagnostic hyperuniformity measures are positively correlated with one another. The same diagnostic measures can be used to detect the degree to which imperfections in nearly hyperuniform systems cause deviations from perfect hyperuniformity. Moreover, the capacity to identify hyperuniform scaling regimes should be particularly useful in analyzing experimentally or computationally generated samples that are necessarily of finite size.

Kinetic Frustration Effects on Dense Two-Dimensional Packings of Convex Particles and Their Structural Characteristics

The study of hard-particle packings is of fundamental importance in physics, chemistry, cell biology, and discrete geometry. Much of the previous work on hard-particle packings concerns their densest possible arrangements. By contrast, we examine kinetic effects inevitably present in both numerical and experimental packing protocols. Specifically, we determine how changing the compression/shear rate of a two-dimensional packing of noncircular particles causes it to deviate from its densest possible configuration, which is always periodic. The adaptive shrinking cell (ASC) optimization scheme maximizes the packing fraction of a hard-particle packing by first applying random translations and rotations to the particles and then isotropically compressing and shearing the simulation box repeatedly until a possibly jammed state is reached. We use a stochastic implementation of the ASC optimization scheme to mimic different effective time scales by varying the number of particle moves between compressions/shears. We generate dense, effectively jammed, monodisperse, two-dimensional packings of obtuse scalene triangle, rhombus, curved triangle, lens, and "ice cream cone" (a semicircle grafted onto an isosceles triangle) shaped particles, with a wide range of packing fractions and degrees of order. To quantify these kinetic effects, we introduce the kinetic frustration index K, which measures the deviation of a packing from its maximum possible packing fraction. To investigate how kinetics affect short- and long-range ordering in these packings, we compute their spectral densities χ̃_V(k) and characterize their contact networks. We find that kinetic effects are most significant when the particles have greater asphericity, less curvature, and less rotational symmetry. This work may be relevant to the design of laboratory packing protocols.

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.

Structural universality in disordered packings with size and shape polydispersity

We numerically investigate disordered jammed packings with both size and shape polydispersity, using frictionless superellipsoidal particles. We implement the set Voronoi tessellation technique to evaluate the local specific volume, i.e., the ratio of cell volume over particle volume, for each individual particle. We focus on the average structural properties for different types of particles binned by their sizes and shapes. We generalize the basic observation that the larger particles are locally packed more densely than the smaller ones in a polydisperse-sized packing into systems with coupled particle shape dispersity. For this purpose, we define the normalized free volume vf to measure the local compactness of a particle and study its dependency on the normalized particle size A. The definition of vf relies on the calibrated monodisperse specific volume for a certain particle shape. For packings with shape dispersity, we apply the previously introduced concept of equivalent diameter for a non-spherical particle to define A properly. We consider three systems: (A) linear superposition states of mixed-shape packings, (B) merely polydisperse-sized packings, and (C) packings with coupled size and shape polydispersity. For (A), the packing is simply considered as a mixture of different subsystems corresponding to monodisperse packings for different shape components, leading to A = 1, and vf = 1 by definition. We propose a concise model to estimate the shape-dependent factor αc, which defines the equivalent diameter for a certain particle. For (B), vf collapses as a function of A, independent of specific particle shape and size polydispersity. Such structural universality is further validated by a mean-field approximation. For (C), we find that the master curve vf(A) is preserved when particles possess similar αc in a packing. Otherwise, the dispersity of αc among different particles causes the deviation from vf(A). These findings show that a polydisperse packing can be estimated as the combination of various building blocks, i.e., bin components, with a universal relation vf(A).

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.

Hard convex lens-shaped particles: Characterization of dense disordered packings

Among the family of hard convex lens-shaped particles (lenses), the one with aspect ratio equal to 2/3 is "optimal" in the sense that the maximally random jammed (MRJ) packings of such lenses achieve the highest packing fraction ϕ_{MRJ}≃0.73 [G. Cinacchi and S. Torquato, Soft Matter 14, 8205 (2018)1744-683X10.1039/C8SM01519H]. This value is only a few percent lower than ϕ_{DKP}=0.76210⋯, the packing fraction of the corresponding densest-known crystalline (degenerate) packings [G. Cinacchi and S. Torquato, J. Chem. Phys. 143, 224506 (2015)JCPSA60021-960610.1063/1.4936938]. By exploiting the appreciably reduced propensity that a system of such optimal lenses has to positionally and orientationally order, disordered packings of them are progressively generated by a Monte Carlo method-based procedure from the dilute equilibrium isotropic fluid phase to the dense nonequilibrium MRJ state. This allows us to closely monitor how the (micro)structure of these packings changes in the process of formation of the MRJ state. The gradual changes undergone by the many structural descriptors calculated here can coherently and consistently be traced back to the gradual increase in contacts between the hard particles until the isostatic mean value of ten contact neighbors per lens is reached at the effectively hyperuniform MRJ state. Compared to the MRJ state of hard spheres, the MRJ state of such optimal lenses is denser (less porous), more disordered, and rattler-free. This set of characteristics makes them good glass formers. It is possible that this conclusion may also hold for other hard convex uniaxial particles with a correspondingly similar aspect ratio, be they oblate or prolate, and that, by using suitable biaxial variants of them, that set of characteristics might further improve.

Predicting maximally random jammed packing density of non-spherical hard particles via analytical continuation of fluid equation of state

We developed a formalism for accurately predicting the density of MRJ packing state of a wide spectrum of congruent non-spherical hard particles in 3D via analytical fluid EOS.

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.

DEM simulation of the local ordering of tetrahedral granular matter

The formation and growth of local order clusters in a tetrahedral granular assembly driven by 3D mechanical vibration were captured in DEM (discrete element method) dynamic simulation using a multi-sphere model. Two important kinds of clusters, dimer and wagon wheel structures, were observed based on which the growth behavior and mechanism of each local cluster with different orientations/structures were investigated. The results show that during vibration, dimer clusters are formed first and then most of them grow into linear trimers and tetramers. Wagon wheel clusters are also frequently observed that grow into hexamers and, further, octamer and nonamer local clusters. Coordination number (CN) evolution indicates that the decrease of local mean CN can be regarded as the signal for the formation of local clusters in the tetrahedral particle packing system. Nematic order metric analysis shows that although the two basic structures (dimer and wagon wheel structures) grow into complex local clusters during packing densification, these local clusters are randomly distributed in the tetrahedral particle packing system. Stress analysis indicates that the dimer-based local clusters are mostly formed in the compaction state of the tetrahedral particle packing system during the vibrated packing densification process. In comparison, the wagon wheel-based local clusters need much stronger interaction forces from tetrahedral particles during vibrated packing densification.

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.

Hard convex lens-shaped particles: metastable, glassy and jammed states

We generate and study dense positionally and/or orientationally disordered, including jammed, monodisperse packings of hard convex lens-shaped particles (lenses). Relatively dense isotropic fluid configurations of lenses of various aspect ratios are slowly compressed via a Monte Carlo method based procedure. Under this compression protocol, while 'flat' lenses form a nematic fluid phase (where particles are positionally disordered but orientationally ordered) and 'globular' lenses form a plastic solid phase (where particles are positionally ordered but orientationally disordered), 'intermediate', neither 'flat' nor 'globular', lenses do not form either mesophase. In general, a crystal solid phase (where particles are both positionally and orientationally ordered) does not spontaneously form during lengthy numerical simulation runs. In correspondence to those volume fractions at which a transition to the crystal solid phase would occur in equilibrium, a 'downturn' is observed in the inverse compressibility factor versus volume fraction curve beyond which this curve behaves essentially linearly. This allows us to estimate the volume fraction at jamming of the dense non-crystalline packings so generated. These packings are nematic for 'flat' lenses and plastic for 'globular' lenses, while they are robustly isotropic for 'intermediate' lenses, as confirmed by the calculation of the τ order metric, among other quantities. The structure factors S(k) of the corresponding jammed states tend to zero as the wavenumber k goes to zero, indicating they are effectively hyperuniform (i.e., their infinite-wavelength density fluctuations are anomalously suppressed). Among all possible lens shapes, 'intermediate' lenses with aspect ratio around 2/3 are special because they are those that reach the highest volume fractions at jamming while being positionally and orientationally disordered and these volume fractions are as high as those reached by nematic jammed states of 'flat' lenses and plastic jammed states of 'globular' lenses. All of their attributes, taken together, make such 'intermediate' lens packings particularly good glass-forming materials.

Evolutions of packing properties of perfect cylinders under densification and crystallization

Cylindrical particles are ubiquitous in nature and industry, and a cylinder is a representative shape of rod-like particles. However, the disordered packing results of cylinders in previous studies are quite inconsistent with each other. In this work, we obtain the MRJ (maximally random jammed) packings and the MDRPs (maximally dense random packings) of perfect cylinders with the aspect ratio (height/diameter) 0.2 ≤ w ≤ 6.0 using the ASC (adaptive shrinking cell) algorithm and the IMC (inverse Monte Carlo) method, respectively. The optimal aspect ratio corresponding to the maximal packing density is w = 0.9 in the MRJ state, while the value is w = 1.2 in the MDRP state. Then we investigate the evolutions of packing properties of perfect cylinders under densification and crystallization. We compare the different final packing states generated via the two methods with different compression rates and order constraints. In the densification procedure, we generate jammed and random packings of cylinders with various compression rates via the ASC and IMC method, respectively. When decreasing the compression rate, we find that the packing density increases but the optimal w remains the same in both methods. In the crystallization procedure, the order constraint in the IMC method is gradually released which means the degree of order in the packings is allowed to increase, and we find that the optimal w shifts from 1.2 to 0.9 while the packing density increases as well. Meanwhile, the random packings evolve into the jammed packings in the crystallization procedure which reflects the competition mechanism between the randomness and jamming. These results also indicate that the optimal w is solely related to the degree of order in the cylinder packings but not determined by the protocol or packing density. Furthermore, a uniform shape elongation effect on the random-packing densities of various shaped particles is found via a new proposed definition of the scaled aspect ratio. Finally, a rough linear relationship between the mean and standard deviation of the reduced Voronoi cell volumes is obtained only for the random packings. Our findings should lead to a better understanding toward the jammed and random packings and are helpful in guiding the granular material design.

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.

Hypostatic jammed packings of frictionless nonspherical particles

We perform computational studies of static packings of a variety of nonspherical particles including circulo-lines, circulo-polygons, ellipses, asymmetric dimers, dumbbells, and others to determine which shapes form packings with fewer contacts than degrees of freedom (hypostatic packings) and which have equal numbers of contacts and degrees of freedom (isostatic packings), and to understand why hypostatic packings of nonspherical particles can be mechanically stable despite having fewer contacts than that predicted from naive constraint counting. To generate highly accurate force- and torque-balanced packings of circulo-lines and cir-polygons, we developed an interparticle potential that gives continuous forces and torques as a function of the particle coordinates. We show that the packing fraction and coordination number at jamming onset obey a masterlike form for all of the nonspherical particle packings we studied when plotted versus the particle asphericity A, which is proportional to the ratio of the squared perimeter to the area of the particle. Further, the eigenvalue spectra of the dynamical matrix for packings of different particle shapes collapse when plotted at the same A. For hypostatic packings of nonspherical particles, we verify that the number of "quartic" modes along which the potential energy increases as the fourth power of the perturbation amplitude matches the number of missing contacts relative to the isostatic value. We show that the fourth derivatives of the total potential energy in the directions of the quartic modes remain nonzero as the pressure of the packings is decreased to zero. In addition, we calculate the principal curvatures of the inequality constraints for each contact in circulo-line packings and identify specific types of contacts with inequality constraints that possess convex curvature. These contacts can constrain multiple degrees of freedom and allow hypostatic packings of nonspherical particles to be mechanically stable.

Shape-dependent effective diffusivity in packings of hard cubes and cuboids compared with spheres and ellipsoids

We performed computational screening of effective diffusivity in different configurations of cubes and cuboids compared with spheres and ellipsoids. In total, more than 1500 structures are generated and screened for effective diffusivity. We studied simple cubic and face-centered cubic lattices of spheres and cubes, random configurations of cubes and spheres as a function of volume fraction and polydispersity, and finally random configurations of ellipsoids and cuboids as a function of shape. We found some interesting shape-dependent differences in behavior, elucidating the impact of shape on the targeted design of granular materials.

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.