Bidisperse Disks Research Articles

Designing the pressure-dependent shear modulus using tessellated granular metamaterials.

Jammed packings of granular materials display complex mechanical response. For example, the ensemble-averaged shear modulus 〈G〉 increases as a power law in pressure p for static packings of soft spherical particles that can rearrange during compression. We seek to design granular materials with shear moduli that can either increase or decrease with pressure without particle rearrangements even in the large-system limit. To do this, we construct tessellated granular metamaterials by joining multiple particle-filled cells together. We focus on cells that contain a small number of bidisperse disks in two dimensions. We first study the mechanical properties of individual disk-filled cells with three types of boundaries: periodic boundary conditions (PBC), fixed-length walls (FXW), and flexible walls (FLW). Hypostatic jammed packings are found for cells with FLW, but not in cells with PBC and FXW, and they are stabilized by quartic modes of the dynamical matrix. The shear modulus of a single cell depends linearly on p. We find that the slope of the shear modulus with pressure λ_{c}<0 for all packings in single cells with PBC where the number of particles per cell N≥6. In contrast, single cells with FXW and FLW can possess λ_{c}>0, as well as λ_{c}<0, for N≤16. We show that we can force the mechanical properties of multicell granular metamaterials to possess those of single cells by constraining the end points of the outer walls and enforcing an affine shear response. These studies demonstrate that tessellated granular metamaterials provide a platform for the design of soft materials with specified mechanical properties.

The Edwards volume ensemble in cyclically sheared granular experiments.

We experimentally investigate the Edwards volume ensemble in cyclically sheared bidisperse disks of two friction coefficients (μ ≈ 0.3 and μ → ∞) subjected to a range of shear amplitudes γm. Despite the local and global anisotropy, hysteresis, and the potential long-range correlation of the free volume, the Edwards volume ensemble surprisingly provides an excellent statistical description of disk packings in cyclically sheared systems. Our finding can be better understood from the comprehensive analysis of the geometric and statistical properties of Voronoi cells of individual particles. First, the average degrees of anisotropy of Voronoi cells are weak at both the microscopic and macroscopic scales within a range of shear amplitudes γm of up to γm = 12% regardless of the inter-particle friction coefficients μ even though the azimuthal distributions of the Voronoi cell depend on μ. Second, there is only negligible hysteresis of global compactivity and volume fluctuations. Finally, the spatial correlations of the free volume and the orientation are weakly anisotropic and short ranged for practical purposes. Both results are independent of μ. Interestingly, our free-volume statistical results are consistent with the simple physical picture that the free volume is directly proportional to the compactivity.

Force chain structure in a rod-withdrawn granular layer

When a rod is vertically withdrawn from a granular layer, oblique force chains can be developed by effective shearing. In this study, the force-chain structure in a rod-withdrawn granular layer was experimentally investigated using a photoelastic technique. The rod is vertically withdrawn from a two-dimensional granular layer consisting of bidisperse photoelastic disks. During the withdrawal, the development process of force chains is visualized by the photoelastic effect. By systematic analysis of photoelastic images, force chain structures newly developed by the rod withdrawing are identified and analyzed. In particular, the relation between the rod-withdrawing force [Formula: see text], total force-chains force [Formula: see text], and their average orientation [Formula: see text] are discussed. We find that the oblique force chains are newly developed by withdrawing. The force-chain angle [Formula: see text] is almost constant (approximately [Formula: see text] from the horizontal), and the total force [Formula: see text] gradually increases by the withdrawal. In addition, [Formula: see text] shows a clear correlation with [Formula: see text].

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.

Structured randomness: jamming of soft discs and pins.

Simulations are used to find the zero temperature jamming threshold, φj, for soft, bidisperse disks in the presence of small fixed particles, or "pins", arranged in a lattice. The presence of pins leads, as one expects, to a decrease in φj. Structural properties of the system near the jamming threshold are calculated as a function of the pin density. While the correlation length exponent remains ν = 1/2 at low pin densities, the system is mechanically stable with more bonds, yet fewer contacts than the Maxwell criterion implies in the absence of pins. In addition, as pin density increases, novel bond orientational order and long-range spatial order appear, which are correlated with the square symmetry of the pin lattice.

Dynamic phases, stratification, laning, and pattern formation for driven bidisperse disk systems in the presence of quenched disorder

Using numerical simulations, we examine the dynamics of driven two-dimensional bidisperse disks flowing over quenched disorder. The system exhibits a series of distinct dynamical phases as a function of applied driving force and packing fraction, including a phase-separated state as well as a smectic state with liquid-like or polycrystalline features. At low driving forces, we find a clogged phase with an isotropic density distribution, while at intermediate driving forces the disks separate into bands of high and low densities with either liquid-like or polycrystalline structure in the high-density bands. In addition to the density phase separation, we find that in some cases there is a fractionation of the disk species, particularly when the disk size ratio is large. The species phase-separated regimes form a variety of patterns such as large disks separated by chains of smaller disks. Our results show that the formation of laning states can be enhanced by tuning the ratio of disk radius of the two species, due to the clumping of small disks in the interstitial regions between the large disks. This system could be experimentally realized using sterically interacting colloidal particles suspended in a viscous fluid driven over random pinning arrays or granular matter suspended in fluid moving over a random landscape.

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.

Local and global avalanches in a two-dimensional sheared granular medium.

We present the experimental and numerical studies of a two-dimensional sheared amorphous material composed of bidisperse photoelastic disks. We analyze the statistics of avalanches during shear including the local and global fluctuations in energy and changes in particle positions and orientations. We find scale-free distributions for these global and local avalanches denoted by power laws whose cutoffs vary with interparticle friction and packing fraction. Different exponents are found for these power laws depending on the quantity from which variations are extracted. An asymmetry in time of the avalanche shapes is evidenced along with the fact that avalanches are mainly triggered by the shear bands. A simple relation independent of the intensity is found between the number of local avalanches and the global avalanches they form. We also compare these experimental and numerical results for both local and global fluctuations to predictions from mean-field and depinning theories.

Anomalous stress fluctuations in athermal two-dimensional amorphous solids.

We numerically study the local stress distribution within athermal, isotropically stressed, mechanically stable, packings of bidisperse frictionless disks above the jamming transition in two dimensions. Considering the Fourier transform of the local stress, we find evidence for algebraically increasing fluctuations in both isotropic and anisotropic components of the stress tensor at small wave numbers, contrary to recent theoretical predictions. Such increasing fluctuations imply a lack of self-averaging of the stress on large length scales. The crossover to these increasing fluctuations defines a length scale ℓ_{0}, however, it appears that ℓ_{0} does not vary much with packing fraction ϕ, nor does ℓ_{0} seem to be diverging as ϕ approaches the jamming ϕ_{J}. We also find similar large length scale fluctuations of stress in the inherent states of a quenched Lennard-Jones liquid, leading us to speculate that such fluctuations may be a general property of amorphous solids in two dimensions.

Clogging and jamming transitions in periodic obstacle arrays.

We numerically examine clogging transitions for bidisperse disks flowing through a two-dimensional periodic obstacle array. We show that clogging is a probabilistic event that occurs through a transition from a homogeneous flowing state to a heterogeneous or phase-separated jammed state where the disks form dense connected clusters. The probability for clogging to occur during a fixed time increases with increasing particle packing and obstacle number. For driving at different angles with respect to the symmetry direction of the obstacle array, we show that certain directions have a higher clogging susceptibility. It is also possible to have a size-specific clogging transition in which one disk size becomes completely immobile while the other disk size continues to flow.

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 etal. [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.

Search for hyperuniformity in mechanically stable packings of frictionless disks above jamming.

We numerically simulate mechanically stable packings of soft-core, frictionless, bidisperse disks in two dimensions, above the jamming packing fraction ϕ(J). For configurations with a fixed isotropic global stress tensor, we investigate the fluctuations of the local packing fraction ϕ(r) to test whether such configurations display the hyperuniformity that has been claimed to exist exactly at ϕ(J). For our configurations, generated by a rapid quench protocol, we find that hyperuniformity persists only out to a finite length scale and that this length scale appears to remain finite as the system stress decreases towards zero, i.e., towards the jamming transition. Our result suggests that the presence of hyperuniformity at jamming may be sensitive to the specific protocol used to construct the jammed configurations.

Stress Relaxation for Granular Materials near Jamming under Cyclic Compression.

We have explored isotropically jammed states of semi-2D granular materials through cyclic compression. In each compression cycle, systems of either identical ellipses or bidisperse disks transition between jammed and unjammed states. We determine the evolution of the average pressure P and structure through consecutive jammed states. We observe a transition point ϕ_{m} above which P persists over many cycles; below ϕ_{m}, P relaxes slowly. The relaxation time scale associated with P increases with packing fraction, while the relaxation time scale for collective particle motion remains constant. The collective motion of the ellipses is hindered compared to disks because of the rotational constraints on elliptical particles.

Statistics of conserved quantities in mechanically stable packings of frictionless disks above jamming

We numerically simulate mechanically stable packings of soft-core, frictionless, bidisperse disks in two dimensions, above the jamming packing fraction ϕ(J). For configurations with a fixed isotropic global stress tensor, we compute the averages, variances, and correlations of conserved quantities (stress Γ(C), force-tile area A(C), Voronoi volume V(C), number of particles N(C), and number of small particles N(sC)) on compact subclusters of particles C, as a function of the cluster size and the global system stress. We find several significant differences depending on whether the cluster C is defined by a fixed radius R or a fixed number of particles M. We comment on the implications of our findings for maximum entropy models of jammed packings.

Spatiotemporal chaotic unjamming and jamming in granular avalanches.

We have investigated the spatiotemporal chaotic dynamics of unjamming and jamming of particles in a model experiment – a rotating drum partially filled with bidisperse disks to create avalanches. The magnitudes of the first Lyapunov vector δu(t) and velocity v(t) of particles are directly measured for the first time to yield insights into their spatial correlation Cδu,v, which is on statistical average slightly larger near the unjamming than the value near the jamming transition. These results are consistent with the recent work of Banigan et al (Nature Phys. 2013), and it is for the first time to validate their theoretical models in a real scenario. v(t) shows rich dynamics: it grows exponentially for unstable particles and keeps increasing despite stochastic interactions; after the maximum, it decays with large fluctuations. Hence the spatiotemporal chaotic dynamics of avalanche particles are entangled, causing temporal correlations of macroscopic quantities of the system. We propose a simple model for these observations.

Aggregation and segregation of confined active particles.

We simulate a model of self-propelled disks with soft repulsive interactions confined to a box in two dimensions. For small rotational diffusion rates, monodisperse disks spontaneously accumulate at the walls. At low densities, interaction forces between particles are strongly inhomogeneous, and a simple model predicts how these inhomogeneities alter the equation of state. At higher densities, collective effects become important. We observe signatures of a jamming transition at a packing fraction ϕ ∼ 0.88, which is also the jamming point for non-active athermal monodisperse disks. At this ϕ, the system develops a critical finite active speed necessary for wall aggregation. At packing fractions above ϕ ∼ 0.6, the pressure decreases with increasing density, suggesting that strong interactions between particles are affecting the equation of state well below the jamming transition. A mixture of bidisperse disks segregates in the absence of any adhesion, identifying a new mechanism that could contribute to cell sorting in embryonic development.

Vibrations of jammed disk packings with Hertzian interactions

Contact breaking and Hertzian interactions between grains can both give rise to nonlinear vibrational response of static granular packings. We perform molecular dynamics simulations at constant energy in 2D of frictionless bidisperse disks that interact via Hertzian spring potentials as a function of energy and measure directly the vibrational response from the Fourier transform of the velocity autocorrelation function. We compare the measured vibrational response of static packings near jamming onset to that obtained from the eigenvalues of the dynamical matrix to determine the temperature above which the harmonic approximation breaks down. We compare packings that interact via single-sided (purely repulsive) and double-sided Hertzian spring interactions to disentangle the effects of the shape of the potential from contact breaking. Our studies show that while Hertzian interactions lead to weak nonlinearities in the vibrational behavior (e.g. the generation of harmonics of the eigenfrequencies of the dynamical matrix), the vibrational response of static packings with Hertzian contact interactions is dominated near jamming by contact breaking as found for systems with repulsive linear spring interactions.

Size-topology correlations in disk packings: terminal bidispersity in order–disorder transitions

In random packings or tilings, the size distribution of individual elements (domains) and the statistics of numbers of neighbours of those domains are strongly correlated. In the case of circular disks forming a random packing in the plane, it has long been known empirically that a certain critical amount of bidispersity avoids crystallization of the packing. We demonstrate how the formalism of a simplified granocentric model allows for an analytical computation of the size-topology correlation as a function of both size ratio and frequency of small disks. The results, obtained without free parameters, are in excellent agreement with the empirical findings of packing simulations concerning critical (terminal) bidispersity. It is also shown that, at equal size variance, the discrete (bidisperse) disk size distributions induce stronger disorder than continuously polydisperse disks.

Size-Topology Relations in Packings of Grains, Emulsions, Foams, and Biological Cells

Particulate packings in 3D are used to study the effects of compression and polydispersity on the geometry of the tiling in these systems. We find that the dependence of the neighbor number on cell size is quasilinear in the monodisperse case and becomes nonlinear above a threshold polydispersity, independent of the method of creation of the tiling. These size-topology relations can be described by a simple analytical theory, which quantifies the effects of positional disorder in the monodisperse case and those of size disorder in the polydisperse case and is applicable in two and three dimensions. The theory thus gives a unifying framework for a wide range of amorphous systems, ranging from biological tissues, foams, and bidisperse disks to compressed emulsions and granular matter.

Tuning jammed frictionless disk packings from isostatic to hyperstatic

We perform extensive computational studies of two-dimensional static bidisperse disk packings using two distinct packing-generation protocols. The first involves thermally quenching equilibrated liquid configurations to zero temperature over a range of thermal quench rates r and initial packing fractions followed by compression and decompression in small steps to reach packing fractions φ(J) at jamming onset. For the second, we seed the system with initial configurations that promote micro- and macrophase-separated packings followed by compression and decompression to φ(J). Using these protocols, we generate more than 10(4) static packings over a wide range of packing fraction, contact number, and compositional and positional order. We find that disordered, isostatic packings exist over a finite range of packing fractions in the large-system limit. In agreement with previous calculations, the most dilute mechanically stable packings with φ min ≈ 0.84 are obtained for r > r*, where r* is the rate above which φ(J) is insensitive to rate. We further compare the structural and mechanical properties of isostatic versus hyperstatic packings. The structural characterizations include the contact number, several order parameters, and mixing ratios of the large and small particles. We find that the isostatic packings are positionally and compositionally disordered (with only small changes in a number of order parameters), whereas bond-orientational and compositional order increase strongly with contact number for hyperstatic packings. In addition, we calculate the static shear modulus and normal mode frequencies (in the harmonic approximation) of the static packings to understand the extent to which the mechanical properties of disordered, isostatic packings differ from partially ordered packings. We find that the mechanical properties of the packings change continuously as the contact number increases from isostatic to hyperstatic.