Existence of isostatic, maximally random jammed monodisperse hard-disk packings

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.

Exceptionally dense and resilient critically jammed polydisperse disk packings

Understanding the way disordered particle packings transition between jammed (rigid) and unjammed (fluid) states is of both great practical importance and strong fundamental interest. The values of critical packing fraction (and other state variables) at the jamming transition are protocol dependent. Here, we demonstrate that this variability can be systematically traced to structural measures of packing, as well as to energy measures inside the jammed regime. A novel generalized simultaneous particle swap algorithm constructs overjammed states of desired energy, which upon decompression lead to predictable critical packing fractions. Thus, for a given set of particle sizes, states with extraordinarily high critical packing fractions can be found efficiently, which sustain substantial shear strain and preserve their special structure over the entire jammed domain. The close relation revealed here between the energy landscape of overjammed soft-particle packings and the behavior near the jamming transition points towards new ways of understanding and constructing disordered materials with exceptional properties.

Estimating random close packing density from circle radius distributions

Circles of a single size can pack together densely in a hexagonal lattice, but adding in size variety disrupts the order of those packings. We conduct simulations which generate dense random packings of circles with specified size distributions and measure the area fraction in each case. While the size distributions can be arbitrary, we find that for a wide range of size distributions the random close-packing area fraction ϕ_{rcp} under this general protocol is determined to high accuracy by the polydispersity and skewness of the size distribution. At low skewness, all packings tend to a minimum packing fraction ϕ_{0}≈0.840 independent of polydispersity. In the limit of high skewness, ϕ_{rcp} becomes independent of skewness, asymptoting to a polydispersity-dependent limit. We show how these results can be predicted from the behavior of bidisperse or bi-Gaussian circle size distributions.

Mixed Nanosphere Assemblies at a Liquid-Liquid Interface

The in-plane packing of gold (Au), polystyrene (PS), and silica (SiO2) spherical nanoparticle (NP) mixtures at a water-oil interface is investigated in situ by UV-vis reflection spectroscopy. All NPs are functionalized with carboxylic acid such that they strongly interact with amine-functionalized ligands dissolved in an immiscible oil phase at the fluid interface. This interaction markedly increases the binding energy of these nanoparticle surfactants (NPSs). The separation distance between the Au NPSs and Au surface coverage are measured by the maximum plasmonic wavelength (λmax) and integrated intensities as the assemblies saturate for different concentrations of non-plasmonic (PS/SiO2) NPs. As the PS/SiO2 content increases, the time to reach intimate Au NP contact also increases, resulting from their hindered mobility. λmax changes within the first few minutes of adsorption due to weak attractive inter-NP forces. Additionally, a sharper peak in the reflection spectrum at NP saturation reveals tighter Au NP packing for assemblies with intermediate non-plasmonic NP content. Grazing incidence small angle X-ray scattering (GISAXS) and scanning electron microscopy (SEM) measurements confirm a decrease in Au NP domain size for mixtures with larger non-plasmonic NP content. The results demonstrate a simple means to probe interfacial phase separation behavior using in situ spectroscopy as interfacial structures densify into jammed, phase-separated NP films.

Dynamical Approach to the Jamming Problem

A simple dynamical model, biased random organization (BRO), appears to produce configurations known as random close packing (RCP) as BRO's densest critical point in dimension d=3. We conjecture that BRO likewise produces RCP in any dimension; if so, then RCP does not exist in d=1-2 (where BRO dynamics lead to crystalline order). In d=3-5, BRO produces isostatic configurations and previously estimated RCP volume fractions 0.64, 0.46, and 0.30, respectively. For all investigated dimensions (d=2-5), we find that BRO belongs to the Manna universality class of dynamical phase transitions by measuring critical exponents associated with the steady-state activity and the long-range density fluctuations. Additionally, BRO's distribution of near contacts (gaps) displays behavior consistent with the infinite-dimensional theoretical treatment of RCP when d≥4. The association of BRO's densest critical configurations with random close packing implies that RCP's upper-critical dimension is consistent with the Manna class d_{uc}=4.

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.

A geometric probabilistic approach to random packing of hard disks in a plane

In this paper, the random packing fraction of hard disks in a plane is analyzed, following a geometric probabilistic approach. First, the random close packing (RCP) of equally sized disks is modelled. Subsequently, following the same methodology, a simple, statistical geometric model is proposed for the random loose packing (RLP) of monodisperse disks. This very basic derivation of RLP leads to a packing value (≈0.66) that is in very good agreement with values that have been obtained previously for 2D disk packings. The present geometrical model also enables a closed-form expression for the contact (coordination) number as a function of the packing density at different states of compaction. These predictions are thoroughly compared with empirical and simulation results, among others the Rényi parking model, yielding good agreement.

Estimating random close packing in polydisperse and bidisperse hard spheres via an equilibrium model of crowding

We show that an analogy between crowding in fluid and jammed phases of hard spheres captures the density dependence of the kissing number for a family of numerically generated jammed states. We extend this analogy to jams of mixtures of hard spheres in d = 3 dimensions and, thus, obtain an estimate of the random close packing volume fraction, ϕ_RCP, as a function of size polydispersity. We first consider mixtures of particle sizes with discrete distributions. For binary systems, we show agreement between our predictions and simulations using both our own results and results reported in previous studies, as well as agreement with recent experiments from the literature. We then apply our approach to systems with continuous polydispersity using three different particle size distributions, namely, the log-normal, Gamma, and truncated power-law distributions. In all cases, we observe agreement between our theoretical findings and numerical results up to rather large polydispersities for all particle size distributions when using as reference our own simulations and results from the literature. In particular, we find ϕ_RCP to increase monotonically with the relative standard deviation, s_σ, of the distribution and to saturate at a value that always remains below 1. A perturbative expansion yields a closed-form expression for ϕ_RCP that quantitatively captures a distribution-independent regime for s_σ < 0.5. Beyond that regime, we show that the gradual loss in agreement is tied to the growth of the skewness of size distributions.

Structural Measures as Guides to Ultrastable States in Overjammed Packings

Jammed, disordered packings of given sets of particles possess a multitude of equilibrium states with different mechanical properties. Identifying and constructing desired states, e.g., of superior stability, is a complex task. Here, we show that in two-dimensional particle packings the energy of all metastable states (inherent structures) is reliably classified by simple scalar measures of local steric packing. These structural measures are insensitive to the particle interaction potential and so robust that they can be used to guide a modified swap algorithm that anneals polydisperse packings toward low-energy metastable states exceptionally fast. The low-energy states are extraordinarily stable against applied shear, so that the approach also efficiently identifies ultrastable packings.

Dense packings of geodesic hard ellipses on a sphere

Packing problems are abundant in nature and have been researched thoroughly both experimentally and in numerical models. In particular, packings of anisotropic, elliptical particles often emerge in models of liquid crystals, colloids, and granular and jammed matter. While most theoretical studies on anisotropic particles have thus far dealt with packings in Euclidean geometry, there are many experimental systems where anisotropically-shaped particles are confined to a curved surface, such as Pickering emulsions stabilized by ellipsoidal particles or protein adsorbates on lipid vesicles. Here, we study random close packing configurations in a two-dimensional model of spherical geodesic ellipses. We focus on the interplay between finite-size effects and curvature that is most prominent at smaller system sizes. We demonstrate that on a spherical surface, monodisperse ellipse packings are inherently disordered, with a non-monotonic dependence of both their packing fraction and the mean contact number on the ellipse aspect ratio, as has also been observed in packings of ellipsoids in both 2D and 3D flat space. We also point out some fundamental differences with previous Euclidean studies and discuss the effects of curvature on our results. Importantly, we show that the underlying spherical surface introduces frustration and results in disordered packing configurations even in systems of monodispersed particles, in contrast to the 2D Euclidean case of ellipse packing. This demonstrates that closed curved surfaces can be effective at introducing disorder in a system and could facilitate the study of monodispersed random packings.

Relaxation and Aging of Nanosphere Assemblies at a Water-Oil Interface

The relaxation and aging of an assembly of spherical nanoparticles (NPs) at a water-oil interface are characterized in situ by grazing incidence X-ray photon correlation spectroscopy. The dynamics of the interfacial assembly is measured while the interface saturates with NPs. Weak attractions between NPs lead to gel-like structures in the assembly, where the in-plane ordering is inhibited by the broad size distribution of the NPs. Structural rearrangements on the length scale of the NP-NP center-to-center distances proceed by intermittent fluctuations instead of continuous cooperative motions. The coexistence of rapid and slow NP populations is confirmed, as commonly observed in soft glass-forming materials. Dynamics are increasingly slowed as the NPs initially segregate to the locally clustered interface. The structural relaxation of the NPs in these localized clusters is 5 orders of magnitude slower than that of free particles in the bulk. When the interface is nearly saturated, the time for relaxation increases suddenly due to the onset of local jamming, and the dynamics slow exponentially afterward until the system reaches collective jamming by cooperative rearrangements. This investigation provides insights into structural relaxations near the glass transition and the evolution of the structure and dynamics of the assemblies as they transition from an isotropic liquid to a dense disordered film.

Descriptor-free unsupervised learning method for local structure identification in particle packings

Local structure identification is of great importance in many scientific and engineering fields. However, mathematical and supervised learning methods mostly rely on specific descriptors of local structure and can only be applied to particular packing configurations. In this work, we propose an improved unsupervised learning method, which is descriptor-free, for local structure identification in particle packing. The point cloud is used as the input of the improved method, which directly comes from spatial positions of particles and does not rely on specific descriptors. The improved method constructs an autoencoder based on the point cloud network combined with Gaussian mixture models for dimension reduction and clustering. Numerical examples show that the improved method performs well in local structure identification of quasicrystal disk and sphere packings, achieving comparable accuracy with previous methods. For disordered packings which have been considered nearly having no local structures, the improved method identifies a nontrivial 7-neighbor motif in the maximally dense random packing of disks, and finds acentric structural motifs in the random close packing of spheres, which demonstrate the ability on identification of new and unknown local structures. The improved unsupervised learning method would help mining information from massive simulation and experimental results, as well as devising new order parameters for particle packings.

On the origins of order

A cardinal feature common to embryonic development and tissue reorganization, as well as to wound healing and cancer cell invasion, is collective cellular migration. During collective migratory events the phenomena of cell jamming and unjamming are increasingly recognized, and underlying mechanical, genomic, transcriptional, and signaling events are increasingly coming to light. In this brief perspective I propose a synthesis that brings together in a new way two key concepts. On the one hand, it has been suggested that the unjammed phase of the cellular collective evolved under a selective pressure favoring fluid-like migratory dynamics as would be required so as to accommodate episodes of tissue evolution, development, plasticity, and repair. Being dynamic, such an unjammed migratory phase is expected to be energetically expensive compared with the jammed non-migratory phase, which is presumed to have evolved under a selective pressure favoring a solid-like homeostatic regime that, by comparison, is energetically economical and mechanically stable. On the other hand, well before the discovery of cell jamming and unjamming Kauffman proposed the general biological principle that living systems exist in a solid regime near the edge of chaos, and that natural selection achieves and sustains such a poised state. Here I propose that, in certain systems at least, this poised solid-like state as predicted in the abstract by Kauffman is realized in the particular by the jammed regime just at the brink of unjamming.

Revealing the characteristic length of random close packing via critical-like random pinning

By randomly pinning particles in fluidized states and finding the local energy minima, we form static packings of mono-disperse disks that resemble random close packing, when only n_c = 2.6% of the particles are pinned. The packings are isostatic and exhibit typical critical scalings of the jamming transition. The non-triviality of n_c is manifested mainly in two aspects. First, n_c acts as a critical point, leading to bifurcated critical scalings in its vicinity. The criticality of n_c is also demonstrated in the packings of weakly polydisperse disks. Second, n_c sets a length scale in agreement with the characteristic length of random close packing. With robust evidence, we show that this agreement is generally true for both mono- and poly-disperse particles and in both two and three dimensions. The randomness inherited from fluidized states by random pinning thus interprets the randomness of random close packing from a unique perspective.

Are cell jamming and unjamming essential in tissue development?

The last decade has seen a surge of evidence supporting the existence of the transition of the multicellular tissue from a collective material phase that is regarded as being jammed to a collective material phase that is regarded as being unjammed. The jammed phase is solid-like and effectively 'frozen', and therefore is associated with tissue homeostasis, rigidity, and mechanical stability. The unjammed phase, by contrast, is fluid-like and effectively 'melted', and therefore is associated with mechanical fluidity, plasticity and malleability that are required in dynamic multicellular processes that sculpt organ microstructure. Such multicellular sculpturing, for example, occurs during embryogenesis, growth and remodeling. Although unjamming and jamming events in the multicellular collective are reminiscent of those that occur in the inert granular collective, such as grain in a hopper that can flow or clog, the analogy is instructive but limited, and the implications for cell biology remain unclear. Here we ask, are the cellular jamming transition and its inverse --the unjamming transition-- mere epiphenomena? That is, are they dispensable downstream events that accompany but neither cause nor quench these core multicellular processes? Drawing from selected examples in developmental biology, here we suggest the hypothesis that, to the contrary, the graded departure from a jammed phase enables controlled degrees of malleability as might be required in developmental dynamics. We further suggest that the coordinated approach to a jammed phase progressively slows those dynamics and ultimately enables long-term mechanical stability as might be required in the mature homeostatic multicellular tissue.

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.

Random Close Packing as a Dynamical Phase Transition

Sphere packing is an ancient problem. The densest packing is known to be a face-centered cubic (FCC) crystal, with space-filling fraction ϕ_{FCC}=π/sqrt[18]≈0.74. The densest "random packing," random close packing (RCP), is yet ill defined, although many experiments and simulations agree on a value ϕ_{RCP}≈0.64. We introduce a simple absorbing-state model, biased random organization (BRO), which exhibits a Manna class dynamical phase transition between absorbing and active states that has as its densest critical point ϕ_{c_{max}}≈0.64≈ϕ_{RCP} and, like other Manna class models, is hyperuniform at criticality. The configurations we obtain from BRO appear to be structurally identical to RCP configurations from other protocols. This leads us to conjecture that the highest-density absorbing state for an isotropic biased random organization model produces an ensemble of configurations that characterizes the state conventionally known as RCP.

Does a central limit theorem hold for the k-skeleton of Poisson hyperplanes in hyperbolic space?

Poisson processes in the space of $$(d-1)$$ ( d - 1 ) -dimensional totally geodesic subspaces (hyperplanes) in a d-dimensional hyperbolic space of constant curvature $$-1$$ - 1 are studied. The k-dimensional Hausdorff measure of their k-skeleton is considered. Explicit formulas for first- and second-order quantities restricted to bounded observation windows are obtained. The central limit problem for the k-dimensional Hausdorff measure of the k-skeleton is approached in two different set-ups: (i) for a fixed window and growing intensities, and (ii) for fixed intensity and growing spherical windows. While in case (i) the central limit theorem is valid for all $$d\ge 2$$ d ≥ 2 , it is shown that in case (ii) the central limit theorem holds for $$d\in \{2,3\}$$ d ∈ { 2 , 3 } and fails if $$d\ge 4$$ d ≥ 4 and $$k=d-1$$ k = d - 1 or if $$d\ge 7$$ d ≥ 7 and for general k. Also rates of convergence are studied and multivariate central limit theorems are obtained. Moreover, the situation in which the intensity and the spherical window are growing simultaneously is discussed. In the background are the Malliavin–Stein method for normal approximation and the combinatorial moment structure of Poisson U-statistics as well as tools from hyperbolic integral geometry.

Emergent Field-Driven Robot Swarm States

We present an ecology-inspired form of active matter consisting of a robot swarm. Each robot moves over a planar dynamic resource environment represented by a large light-emitting diode array in search of maximum light intensity; the robots deplete (dim) locally by their presence the local light intensity and seek maximum light intensity. Their movement is directed along the steepest local light intensity gradient; we call this emergent symmetry breaking motion "field drive." We show there emerge dynamic and spatial transitions similar to gas, crystalline, liquid, glass, and jammed states as a function of robot density, resource consumption rates, and resource recovery rates. Paradoxically the nongas states emerge from smooth, flat resource landscapes, not rough ones, and each state can directly move to a glassy state if the resource recovery rate is slow enough, at any robot density.

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.

Critical point of jamming transition in two-dimensional monodisperse systems

The existence of amorphous packings in two-dimensional monodisperse system is a classical unsolved problem. We get the energy minimum state by the energy minimization method of enthalpy under constant pressure conditions. Firstly, we find that there are two peaks in the experiment, which demonstrate the interesting features of the coexistence of crystals and amorphous crystals. And then, we confirm the critical point of jamming transition of the two-dimensional monodisperse is [Formula: see text]. Finally, we prove that the jamming scaling is still satisfied in two-dimensional monodispersed system: [Formula: see text] and vanishes as [Formula: see text], and the boson peak shifts to lower frequencies for less compressed systems.

Bidisperse Nanospheres Jammed on a Liquid Surface

Jammed packings of bidisperse nanospheres were assembled on a nonvolatile liquid surface and visualized to the single-particle scale by using an in situ scanning electron microscopy method. The PEGylated silica nanospheres, mixed at different number fractions and size ratios, had large enough in-plane mobilities prior to jamming to form uniform monolayers reproducibly. From the collected nanometer-resolution images, local order and degree of mixing were assessed by standard metrics. For equimolar mixtures, a large-to-small size ratio of about 1.5 minimized correlated metrics for local orientational and positional order, as previously predicted in simulations of bidisperse disk jamming. Despite monolayer uniformity, structural and depletion interactions caused spheres of a similar size to cluster, a feature evident at size ratios above 2. Uniform nanoparticle monolayers of high packing disorder are sought in many liquid interface technologies, and these experiments outlined key design principles, buttressing extensive theory/simulation literature on the topic.

Stress-induced ordering of two-dimensional packings of elastic spheres

Packing of particles in confined environments is a common problem in multiple fields. Here, based on the two-dimensional Hertzian particle model, we study the packing of deformable spherical particles under compression and reveal the crucial role of stress as an ordering field in regulating particle arrangement. Specifically, under increasing compression, the squeezed particles spontaneously organize into regular packings in the sequence of triangular and square lattices, pentagonal tessellation, and the reentrant triangular lattice. The rich ordered patterns and complex structures revealed in this work suggest a fruitful organizational strategy based on the interplay of external stress and intrinsic elastic instability of particle arrays.

Unjamming and collective migration in MCF10A breast cancer cell lines

Each cell comprising an intact, healthy, confluent epithelial layer ordinarily remains sedentary, firmly adherent to and caged by its neighbors, and thus defines an elemental constituent of a solid-like cellular collective [1,2]. After malignant transformation, however, the cellular collective can become fluid-like and migratory, as evidenced by collective motions that arise in characteristic swirls, strands, ducts, sheets, or clusters [3,4]. To transition from a solid-like to a fluid-like phase and thereafter to migrate collectively, it has been recently argued that cells comprising the disordered but confluent epithelial collective can undergo changes of cell shape so as to overcome geometric constraints attributable to the newly discovered phenomenon of cell jamming and the associated unjamming transition (UJT) [1,2,5-9]. Relevance of the jamming concept to carcinoma cells lines of graded degrees of invasive potential has never been investigated, however. Using classical in vitro cultures of six breast cancer model systems, here we investigate structural and dynamical signatures of cell jamming, and the relationship between them [1,2,10,11]. In order of roughly increasing invasive potential as previously reported, model systems examined included MCF10A, MCF10A.Vector; MCF10A.14-3-3ζ; MCF10.ErbB2, MCF10AT; and MCF10CA1a [12-15]. Migratory speed depended on the particular cell line. Unsurprisingly, for example, the MCF10CA1a cell line exhibited much faster migratory speed relative to the others. But unexpectedly, across different cell lines higher speeds were associated with enhanced size of cooperative cell packs in a manner reminiscent of a peloton [9]. Nevertheless, within each of the cell lines evaluated, cell shape and shape variability from cell-to-cell conformed with predicted structural signatures of cell layer unjamming [1]. Moreover, both structure and migratory dynamics were compatible with previous theoretical descriptions of the cell jamming mechanism [2,10,11,16,17]. As such, these findings demonstrate the richness of the cell jamming mechanism, which is now seen to apply across these cancer cell lines but remains poorly understood.

Building Reconfigurable Devices Using Complex Liquid-Fluid Interfaces

Liquid-fluid interfaces provide a platform both for structuring liquids into complex shapes and assembling dimensionally confined, functional nanomaterials. Historically, attention in this area has focused on simple emulsions and foams, in which surface-active materials such as surfactants or colloids stabilize structures against coalescence and alter the mechanical properties of the interface. In recent decades, however, a growing body of work has begun to demonstrate the full potential of the assembly of nanomaterials at liquid-fluid interfaces to generate functionally advanced, biomimetic systems. Here, a broad overview is given, from fundamentals to applications, of the use of liquid-fluid interfaces to generate complex, all-liquid devices with a myriad of potential applications.

Universal Features of Metastable State Energies in Cellular Matter

Mechanical equilibrium states of cellular matter are overwhelmingly metastable and separated from each other by topology changes. Using theory and simulations, it is shown that for a wide class of energy functionals in 2D, including those describing tissue cell layers, local energy differences between neighboring metastable states as well as global energy differences between initial states and ground states are governed by simple, universal relations. Knowledge of instantaneous length of an edge undergoing a T1 transition is sufficient to predict local energy changes, while the initial edge length distribution yields a successful prediction for the global energy difference. An analytical understanding of the model parameters is provided.

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.

Response of jammed packings to thermal fluctuations

We focus on the response of mechanically stable (MS) packings of frictionless, bidisperse disks to thermal fluctuations, with the aim of quantifying how nonlinearities affect system properties at finite temperature. In contrast, numerous prior studies characterized the structural and mechanical properties of MS packings of frictionless spherical particles at zero temperature. Packings of disks with purely repulsive contact interactions possess two main types of nonlinearities, one from the form of the interaction potential (e.g., either linear or Hertzian spring interactions) and one from the breaking (or forming) of interparticle contacts. To identify the temperature regime at which the contact-breaking nonlinearities begin to contribute, we first calculated the minimum temperatures T_{cb} required to break a single contact in the MS packing for both single- and multiple-eigenmode perturbations of the T=0 MS packing. We find that the temperature required to break a single contact for equal velocity-amplitude perturbations involving all eigenmodes approaches the minimum value obtained for a perturbation in the direction connecting disk pairs with the smallest overlap. We then studied deviations in the constant volume specific heat C[over ¯]_{V} and deviations of the average disk positions Δr from their T=0 values in the temperature regime T_{C[over ¯]_{V}}<T<T_{r}, where T_{r} is the temperature beyond which the system samples the basin of a new MS packing. We find that the deviation in the specific heat per particle ΔC[over ¯]_{V}^{0}/C[over ¯]_{V}^{0} relative to the zero-temperature value C[over ¯]_{V}^{0} can grow rapidly above T_{cb}; however, the deviation ΔC[over ¯]_{V}^{0}/C[over ¯]_{V}^{0} decreases as N^{-1} with increasing system size. To characterize the relative strength of contact-breaking versus form nonlinearities, we measured the ratio of the average position deviations Δr^{ss}/Δr^{ds} for single- and double-sided linear and nonlinear spring interactions. We find that Δr^{ss}/Δr^{ds}>100 for linear spring interactions is independent of system size. This result emphasizes that contact-breaking nonlinearities are dominant over form nonlinearities in the low-temperature range T_{cb}<T<T_{r} for model jammed systems.

Hexagonal Patterning of the Insect Compound Eye: Facet Area Variation, Defects, and Disorder

The regular hexagonal array morphology of facets (ommatidia) in the Drosophila compound eye is accomplished by regulation of cell differentiation and planar cell polarity during development. Mutations in certain genes disrupt regulation, causing a breakdown of this perfect symmetry, so that the ommatidial pattern shows onset of disorder in the form of packing defects. We analyze a variety of such mutants and compare them to normal (wild-type), finding that mutants show increased local variation in ommatidial area, which is sufficient to induce a significant number of defects. A model formalism based on Voronoi construction is developed to predict the observed correlation between ommatidium size variation and the number of defects, and to study the onset of disorder in this system with statistical tools. The model uncovers a previously unknown large-scale systematic size variation of the ommatidia across the eye of both wild-type and mutant animals. Such systematic variation of area, as well as its statistical fluctuations, are found to have distinct effects on eye disorder that can both be quantitatively modeled. Furthermore, the topological order is also influenced by the internal structure of the ommatidia, with cells of greater relative mechanical stiffness providing constraints to ommatidial deformation and thus to defect generation. Without free parameters, the simulation predicts the size-topology correlation for both wild-type and mutant eyes. This work develops formalisms of size-topology correlation that are very general and can be potentially applied to other cellular structures near the onset of disorder.

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

The nonequilibrium process by which hard-particle systems may be compressed into disordered, jammed states has received much attention because of its wide utility in describing a broad class of amorphous materials. While dynamical signatures are known to precede jamming, the task of identifying static structural signatures indicating the onset of jamming have proven more elusive. The observation that compressing hard-particle packings towards jamming is accompanied by an anomalous suppression of density fluctuations (termed "hyperuniformity") has paved the way for the analysis of jamming as an "inverted critical point" in which the direct correlation function c(r), rather than the total correlation function h(r), diverges. We expand on the notion that c(r) provides both universal and protocol-specific information as packings approach jamming. By considering the degree and position of singularities (discontinuities in the nth derivative) as well as how they are changed by the convolutions found in the Ornstein-Zernike equation, we establish quantitative statements about the structure of c(r) with regards to singularities it inherits from h(r). These relations provide a concrete means of identifying features that must be expressed in c(r) if one hopes to reproduce various details in the pair correlation function accurately and provide stringent tests on the associated numerics. We also analyze the evolution of systems of three-dimensional monodisperse hard spheres of diameter D as they approach ordered and disordered jammed configurations. For the latter, we use the Lubachevsky-Stillinger (LS) molecular dynamics and Torquato-Jiao (TJ) sequential linear programming algorithms, which both generate disordered packings, but can show perceptible structural differences. We identify a short-ranged scaling c(r)∝-1/r as r→0 that accompanies the formation of the delta function at c(D) that indicates the formation of contacts in all cases, and show that this scaling behavior is, in this case, a consequence of the growing long rangedness in c(r), e.g., c∝-1/r^{2} as r→∞ for disordered packings. At densities in the vicinity of the freezing density, we find striking qualitative differences in the structure factor S(k) as well as c(r) between TJ- and LS-generated configurations, including the early formation of a delta function at c(D) in the TJ algorithm's packings, indicating the early formation of clusters of particles in near contact. Both algorithms yield structure factors that tend towards zero in the low-wave-number limit as jamming is approached. Correspondingly, we observe the expected power-law decay in c(r) for large r, in agreement with previous theoretical work. Our work advances the notion that static signatures are exhibited by hard-particle packings as they approach jamming and underscores the utility of the direct correlation function as a sensitive means of monitoring for the appearance of an incipient rigid network.

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.

Assembling semiconducting molecules by covalent attachment to a lamellar crystalline polymer substrate

We have investigated the potential of polymers containing precisely spaced side-branches for thin film applications, particularly in the context of organic electronics. Upon crystallization, the side-branches were excluded from the crystalline core of a lamellar crystal. Thus, the surfaces of these crystals were covered by side-branches. By using carboxyl groups as side-branches, which allow for chemical reactions, we could functionalize the crystal with semiconducting molecules. Here, we compare properties of crystals differing in size: small nanocrystals and large single crystals. By assembling nanocrystals on a Langmuir trough, large areas could be covered by monolayers consisting of randomly arranged nanocrystals. Alternatively, we used a method based on local supersaturation to grow large area single crystals of the precisely side-branched polymer from solution. Attachment of the semiconducting molecules to the lamellar surface of large single crystals was possible, however, only after an appropriate annealing procedure. As a function of the duration of the grafting process, the morphology of the resulting layer of semiconducting molecules changed from patchy to compact.

Pattern formations and optimal packing

Patterns of different symmetries may arise after solution to reaction-diffusion equations. Hexagonal arrays, layers and their perturbations are observed in different models after numerical solution to the corresponding initial-boundary value problems. We demonstrate an intimate connection between pattern formations and optimal random packing on the plane. The main study is based on the following two points. First, the diffusive flux in reaction-diffusion systems is approximated by piecewise linear functions in the framework of structural approximations. This leads to a discrete network approximation of the considered continuous problem. Second, the discrete energy minimization yields optimal random packing of the domains (disks) in the representative cell. Therefore, the general problem of pattern formations based on the reaction-diffusion equations is reduced to the geometric problem of random packing. It is demonstrated that all random packings can be divided onto classes associated with classes of isomorphic graphs obtained from the Delaunay triangulation. The unique optimal solution is constructed in each class of the random packings. If the number of disks per representative cell is finite, the number of classes of isomorphic graphs, hence, the number of optimal packings is also finite.

Confined disordered strictly jammed binary sphere packings

Disordered jammed packings under confinement have received considerably less attention than their bulk counterparts and yet arise in a variety of practical situations. In this work, we study binary sphere packings that are confined between two parallel hard planes and generalize the Torquato-Jiao (TJ) sequential linear programming algorithm [Phys. Rev. E 82, 061302 (2010)] to obtain putative maximally random jammed (MRJ) packings that are exactly isostatic with high fidelity over a large range of plane separation distances H, small to large sphere radius ratio α, and small sphere relative concentration x. We find that packing characteristics can be substantially different from their bulk analogs, which is due to what we term "confinement frustration." Rattlers in confined packings are generally more prevalent than those in their bulk counterparts. We observe that packing fraction, rattler fraction, and degree of disorder of MRJ packings generally increase with H, though exceptions exist. Discontinuities in the packing characteristics as H varies in the vicinity of certain values of H are due to associated discontinuous transitions between different jammed states. When the plane separation distance is on the order of two large-sphere diameters or less, the packings exhibit salient two-dimensional features; when the plane separation distance exceeds about 30 large-sphere diameters, the packings approach three-dimensional bulk packings. As the size contrast increases (as α decreases), the rattler fraction dramatically increases due to what we call "size-disparity" frustration. We find that at intermediate α and when x is about 0.5 (50-50 mixture), the disorder of packings is maximized, as measured by an order metric ψ that is based on the number density fluctuations in the direction perpendicular to the hard walls. We also apply the local volume-fraction variance σ(τ)(2)(R) to characterize confined packings and find that these packings possess essentially the same level of hyperuniformity as their bulk counterparts. Our findings are generally relevant to confined packings that arise in biology (e.g., structural color in birds and insects) and may have implications for the creation of high-density powders and improved battery designs.

A Geometric-Structure Theory for Maximally Random Jammed Packings

Maximally random jammed (MRJ) particle packings can be viewed as prototypical glasses in that they are maximally disordered while simultaneously being mechanically rigid. The prediction of the MRJ packing density ϕMRJ, among other packing properties of frictionless particles, still poses many theoretical challenges, even for congruent spheres or disks. Using the geometric-structure approach, we derive for the first time a highly accurate formula for MRJ densities for a very wide class of two-dimensional frictionless packings, namely, binary convex superdisks, with shapes that continuously interpolate between circles and squares. By incorporating specific attributes of MRJ states and a novel organizing principle, our formula yields predictions of ϕMRJ that are in excellent agreement with corresponding computer-simulation estimates in almost the entire α-x plane with semi-axis ratio α and small-particle relative number concentration x. Importantly, in the monodisperse circle limit, the predicted ϕMRJ = 0.834 agrees very well with the very recently numerically discovered MRJ density of 0.827, which distinguishes it from high-density "random-close packing" polycrystalline states and hence provides a stringent test on the theory. Similarly, for non-circular monodisperse superdisks, we predict MRJ states with densities that are appreciably smaller than is conventionally thought to be achievable by standard packing protocols.