Jamming in systems composed of frictionless ellipse-shaped particles

Jamming on convex deformable surfaces

Jamming is a fundamental transition that governs the behavior of particulate media, including sand, foams and dense suspensions. Upon compression, such media change from freely flowing to a disordered, marginally stable solid that exhibits non-Hookean elasticity. While the jamming process is well established for fixed geometries, the nature and dynamics of jamming for a diverse class of soft materials and deformable substrates, including emulsions and biological matter, remains unknown. Here we propose a new scenario, metric jamming, where rigidification occurs on a surface that has been deformed from its ground state. Unlike classical jamming processes that exhibit discrete mechanical transitions, surprisingly we find that metric jammed states possess mechanical properties continuously tunable between those of classically jammed and conventional elastic media. The compact and curved geometry significantly alters the vibrational spectra of the structures relative to jamming in flat Euclidean space, and metric jammed systems also possess new types of vibrational mode that couple particle and shape degrees of freedom. Our work provides a theoretical framework that unifies our understanding of solidification processes that take place on deformable media and lays the groundwork to exploit jamming for the control and stabilization of shape in self-assembly processes.

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.

Theory of rigidity and numerical analysis of density of states of two-dimensional amorphous solids with dispersed frictional grains in the linear response regime

Using the Jacobian matrix, we obtain a theoretical expression of rigidity and the density of states of two-dimensional amorphous solids consisting of frictional grains in the linear response to an infinitesimal strain, in which we ignore the dynamical friction caused by the slip processes of contact points. The theoretical rigidity agrees with that obtained by molecular dynamics simulations. We confirm that the rigidity is smoothly connected to the value in the frictionless limit. We find that there are two modes in the density of states for sufficiently small k_{T}/k_{N}, which is the ratio of the tangential to normal stiffness. Rotational modes exist at low frequencies or small eigenvalues, whereas translational modes exist at high frequencies or large eigenvalues. The location of the rotational band shifts to the high-frequency region with an increase in k_{T}/k_{N} and becomes indistinguishable from the translational band for large k_{T}/k_{N}.

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.

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.

Measure of overlap between two arbitrary ellipses on a sphere

Various packing problems and simulations of hard and soft interacting particles, such as microscopic models of nematic liquid crystals, reduce to calculations of intersections and pair interactions between ellipsoids. When constrained to a spherical surface, curvature and compactness lead to non-trivial behaviour that finds uses in physics, computer science and geometry. A well-known idealized isotropic example is the Tammes problem of finding optimal non-intersecting packings of equal hard disks. The anisotropic case of elliptic particles remains, on the other hand, comparatively unexplored. We develop an algorithm to detect collisions between ellipses constrained to the two-dimensional surface of a sphere based on a solution of an eigenvalue problem. We investigate and discuss topologically distinct ways two ellipses may touch or intersect on a sphere, and define a contact function that can be used for construction of short- and long-range pair potentials.

Energetic rigidity. II. Applications in examples of biological and underconstrained materials

This is the second paper devoted to energetic rigidity, in which we apply our formalism to examples in two dimensions: Underconstrained random regular spring networks, vertex models, and jammed packings of soft particles. Spring networks and vertex models are both highly underconstrained, and first-order constraint counting does not predict their rigidity, but second-order rigidity does. In contrast, spherical jammed packings are overconstrained and thus first-order rigid, meaning that constraint counting is equivalent to energetic rigidity as long as prestresses in the system are sufficiently small. Aspherical jammed packings on the other hand have been shown to be jammed at hypostaticity, which we use to argue for a modified constraint counting for systems that are energetically rigid at quartic order.

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.

The structural, vibrational, and mechanical properties of jammed packings of deformable particles in three dimensions

We investigate the structural, vibrational, and mechanical properties of jammed packings of deformable particles with shape degrees of freedom in three dimensions (3D). Each 3D deformable particle is modeled as a surface-triangulated polyhedron, with spherical vertices whose positions are determined by a shape-energy function with terms that constrain the particle surface area, volume, and curvature, and prevent interparticle overlap. We show that jammed packings of deformable particles without bending energy possess low-frequency, quartic vibrational modes, whose number decreases with increasing asphericity and matches the number of missing contacts relative to the isostatic value. In contrast, jammed packings of deformable particles with non-zero bending energy are isostatic in 3D, with no quartic modes. We find that the contributions to the eigenmodes of the dynamical matrix from the shape degrees of freedom are significant over the full range of frequency and shape parameters for particles with zero bending energy. We further show that the ensemble-averaged shear modulus 〈G〉 scales with pressure P as 〈G〉 ∼ P^β, with β ≈ 0.75 for jammed packings of deformable particles with zero bending energy. In contrast, β ≈ 0.5 for packings of deformable particles with non-zero bending energy, which matches the value for jammed packings of soft, spherical particles with fixed shape. These studies underscore the importance of incorporating particle deformability and shape change when modeling the properties of jammed soft materials.

Testing mean-field theory for jamming of non-spherical particles: contact number, gap distribution, and vibrational density of states

We perform numerical simulations of the jamming transition of non-spherical particles in two dimensions. In particular, we systematically investigate how the physical quantities at the jamming transition point behave when the shapes of the particle deviate slightly from the perfect disks. For efficient numerical simulation, we first derive an analytical expression of the gap function, using the perturbation theory around the reference disks. Starting from disks, we observe the effects of the deformation of the shapes of particles by the n-th-order term of the Fourier series [Formula: see text]. We show that the several physical quantities, such as the number of contacts, gap distribution, and characteristic frequencies of the vibrational density of states, show the power-law behaviors with respect to the linear deviation from the reference disks. The power-law behaviors do not depend on n and are fully consistent with the mean-field theory of the jamming of non-spherical particles. This result suggests that the mean-field theory holds very generally for nearly spherical particles.

Mechanical response of packings of nonspherical particles: A case study of two-dimensional packings of circulo-lines

We investigate the mechanical response of jammed packings of circulo-lines in two spatial dimensions, interacting via purely repulsive, linear spring forces, as a function of pressure P during athermal, quasistatic isotropic compression. The surface of a circulo-line is defined as the collection of points that is equidistant to a line; circulo-lines are composed of a rectangular central shaft with two semicircular end caps. Prior work has shown that the ensemble-averaged shear modulus for jammed disk packings scales as a power law, 〈G(P)〉∼P^{β}, with β∼0.5, over a wide range of pressure. For packings of circulo-lines, we also find robust power-law scaling of 〈G(P)〉 over the same range of pressure for aspect ratios R≳1.2. However, the power-law scaling exponent β∼0.8-0.9 is much larger than that for jammed disk packings. To understand the origin of this behavior, we decompose 〈G〉 into separate contributions from geometrical families, G_{f}, and from changes in the interparticle contact network, G_{r}, such that 〈G〉=〈G_{f}〉+〈G_{r}〉. We show that the shear modulus for low-pressure geometrical families for jammed packings of circulo-lines can both increase and decrease with pressure, whereas the shear modulus for low-pressure geometrical families for jammed disk packings only decreases with pressure. For this reason, the geometrical family contribution 〈G_{f}〉 is much larger for jammed packings of circulo-lines than for jammed disk packings at finite pressure, causing the increase in the power-law scaling exponent for 〈G(P)〉.

Shear response of granular packings compressed above jamming onset

We investigate the mechanical response of jammed packings of repulsive, frictionless spherical particles undergoing isotropic compression. Prior simulations of the soft-particle model, where the repulsive interactions scale as a power law in the interparticle overlap with exponent α, have found that the ensemble-averaged shear modulus 〈G(P)〉 increases with pressure P as ∼P^{(α-3/2)/(α-1)} at large pressures. 〈G〉 has two key contributions: (1) continuous variations as a function of pressure along geometrical families, for which the interparticle contact network does not change, and (2) discontinuous jumps during compression that arise from changes in the contact network. Using numerical simulations, we show that the form of the shear modulus G^{f} for jammed packings within near-isostatic geometrical families is largely determined by the affine response G^{f}∼G_{a}^{f}, where G_{a}^{f}/G_{a0}=(P/P_{0})^{(α-2)/(α-1)}-P/P_{0}, P_{0}∼N^{-2(α-1)} is the characteristic pressure at which G_{a}^{f}=0, G_{a0} is a constant that sets the scale of the shear modulus, and N is the number of particles. For near-isostatic geometrical families that persist to large pressures, deviations from this form are caused by significant nonaffine particle motion. We further show that the ensemble-averaged shear modulus 〈G(P)〉 is not simply a sum of two power laws, but 〈G(P)〉∼(P/P_{c})^{a}, where a≈(α-2)/(α-1) in the P→0 limit and 〈G(P)〉∼(P/P_{c})^{b}, where b≳(α-3/2)/(α-1), above a characteristic pressure that scales as P_{c}∼N^{-2(α-1)}.

Translational and rotational critical-like behaviors in the glass transition of colloidal ellipsoid monolayers

Critical-like behaviors have been found in translational degrees of freedom near the glass transition of spherical particle systems mainly with local polycrystalline structures, but it is not clear if criticality exists in more general glassy systems composed of nonspherical particles without crystalline structures. Here, through experiments and simulations, we show critical-like behaviors in both translational and rotational degrees of freedom in monolayers of monodisperse colloidal ellipsoids in the absence of crystalline orders. We find rich features of the Ising-like criticality in structure and slow dynamics at the ideal glass transition point ϕ₀, showing the thermodynamic nature of glass transition at ϕ₀ A dynamic criticality is found at the mode-coupling critical point ϕ_c for the fast-moving clusters whose critical exponents increase linearly with fragility, reflecting a dynamic glass transition. These results cast light on the glass transition and explain the mystery that the dynamic correlation lengths diverge at two different temperatures.

Translational and rotational velocities in shear-driven jamming of ellipsoidal particles

We study shear-driven jamming of ellipsoidal particles at zero temperature with a focus on the microscopic dynamics. We find that a change from spherical particles to ellipsoids with aspect ratio α=1.02 gives dramatic changes of the microscopic dynamics with much lower translational velocities and a new role for the rotations. Whereas the velocity difference at contacts-and thereby the dissipation-in collections of spheres is dominated by the translational velocities and reduced by the rotations, the same quantity is in collections of ellipsoids instead totally dominated by the rotational velocities. By also examining the effect of different aspect ratios we find that the examined quantities show either a peak or a change in slope at α≈1.2, which thus gives evidence for a crossover between different regions of low and high aspect ratio.

Packing and dynamics of a protein solution approaching the jammed state

The packing of proteins and their collective behavior in crowded media is crucial for the understanding of biological processes. Here we study the structural and dynamical evolution of solutions of the globular protein bovine serum albumin with increasing concentration via drying using small angle X-ray scattering and dynamic light scattering. We probe an evolving correlation peak on the scattering profile, corresponding to the inter-protein distance, ξ, which decreases following a power law of the protein volume fraction, φ. The rate of decrease in ξ becomes faster above a protein concentration of ∼200 mg ml-1 (φ = 0.15). The power law exponent changes from 0.33, which is typical of colloidal or protein solutions, to 0.41. During the entire drying process, we observe the development and the growth of two-step relaxation dynamics with increasing φ as revealed by dynamic light scattering. We find three different regimes of the dependence of ξ as a function of φ. In the dilute regime (φ < 0.22), protein molecules are far apart from each other compared to their size. In this case, the dynamics mainly corresponds to Brownian motion. At an intermediate concentration (0.22 < φ < 0.47), inter-protein distances become comparable to the size of protein molecules, leading to a preferential orientation of the ellipsoidal protein molecules along with a possible deformation. In this regime, the dynamics shows two distinct relaxation times. At a very high concentration (φ > 0.47), the system reaches a jammed state. Subsequently, the secondary relaxation time in this state becomes extremely slow. In this state, the protein molecules have approximately one hydration layer. This study contributes to the understanding of protein molecular packing in crowded environments and the phenomenon of density-driven jamming for soft matter systems.

Packings of frictionless spherocylinders

We present simulation results on the properties of packings of frictionless spherocylindrical particles. Starting from a random distribution of particles in space, a packing is produced by minimizing the potential energy of interparticle contacts until a force-equilibrated state is reached. For different particle aspect ratios α=10⋯40, we calculate contacts z, pressure as well as bulk and shear modulus. Most important is the fraction f_{0}(α) of spherocylinders with contacts at both ends, as it governs the jamming threshold z_{c}(α)=8+2f_{0}(α). These results highlight the important role of the axial "sliding" degree of freedom of a spherocylinder, which is a zero-energy mode but only if no end contacts are present.

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.

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.

Jammed packings of 3D superellipsoids with tunable packing fraction, coordination number, and ordering

We carry out numerical studies of static packings of frictionless superellipsoidal particles in three spatial dimensions. We consider more than 200 different particle shapes by varying the three shape parameters that define superellipsoids. We characterize the structural and mechanical properties of both disordered and ordered packings using two packing-generation protocols. We perform athermal quasi-static compression simulations starting from either random, dilute configurations (Protocol 1) or thermalized, dense configurations (Protocol 2), which allows us to tune the orientational order of the packings. In general, we find that superellipsoid packings are hypostatic, with coordination number zJ < ziso, where ziso = 2df and df = 5 or 6 depending on whether the particles are axi-symmetric or not. Over the full range of orientational order, we find that the number of quartic modes of the dynamical matrix for the packings always matches the number of missing contacts relative to the isostatic value. This result suggests that there are no mechanically redundant contacts for ordered, yet hypostatic packings of superellipsoidal particles. Additionally, we find that the packing fraction at jamming onset for disordered packings of superellipsoidal depends on at least two particle shape parameters, e.g. the asphericity A and reduced aspect ratio β of the particles.

Vibrational properties of two-dimensional dimer packings near the jamming transition

Jammed particulate systems composed of various shapes of particles undergo the jamming transition as they are compressed or decompressed. To date, sphere packings have been extensively studied in many previous works, where isostaticity at the transition and scaling laws with the pressure of various quantities, including the contact number and the vibrational density of states, have been established. Additionally, much attention has been paid to nonspherical packings, and particularly recent work has made progress in understanding ellipsoidal packings. In this work, we study the dimer packings in two dimensions, which have been much less understood than systems of spheres and ellipsoids. We first study the contact number of dimers near the jamming transition. It turns out that packings of dimers have "rotational rattlers," each of which still has a free rotational motion. After correcting this effect, we show that dimers become isostatic at the jamming, and the excess contact number obeys the same critical law and finite-size scaling law as those of spheres. We next study the vibrational properties of dimers near the transition. We find that the vibrational density of states of dimers exhibits two characteristic plateaus that are separated by a peak. The high-frequency plateau is dominated by the translational degree of freedom, while the low-frequency plateau is dominated by the rotational degree of freedom. We establish the critical scaling laws of the characteristic frequencies of the plateaus and the peak near the transition. In addition, we present detailed characterizations of the real space displacement fields of vibrational modes in the translational and rotational plateaus.

Shear-driven flow of athermal, frictionless, spherocylinder suspensions in two dimensions: Stress, jamming, and contacts

We use numerical simulations to study the flow of a bidisperse mixture of athermal, frictionless, soft-core two-dimensional spherocylinders driven by a uniform steady-state shear strain applied at a fixed finite rate. Energy dissipation occurs via a viscous drag with respect to a uniformly sheared host fluid, giving a simple model for flow in a non-Brownian suspension and resulting in a Newtonian rheology. We study the resulting pressure p and deviatoric shear stress σ of the interacting spherocylinders as a function of packing fraction ϕ, strain rate γ[over ̇], and a parameter α that measures the asphericity of the particles; α is varied to consider the range from nearly circular disks to elongated rods. We consider the direction of anisotropy of the stress tensor, the macroscopic friction μ=σ/p, and the divergence of the transport coefficient η_{p}=p/γ[over ̇] as ϕ is increased to the jamming transition ϕ_{J}. From a phenomenological analysis of Herschel-Bulkley rheology above jamming, we estimate ϕ_{J} as a function of asphericity α and show that the variation of ϕ_{J} with α is the main cause for differences in rheology as α is varied; when plotted as ϕ/ϕ_{J}, rheological curves for different α qualitatively agree. However, a detailed scaling analysis of the divergence of η_{p} for our most elongated particles suggests that the jamming transition of spherocylinders may be in a different universality class than that of circular disks. We also compute the number of contacts per particle Z in the system and show that the value at jamming Z_{J} is a nonmonotonic function of α that is always smaller than the isostatic value. We measure the probability distribution of contacts per unit surface length P(ϑ) at polar angle ϑ with respect to the spherocylinder spine and find that as α→0 this distribution seems to diverge at ϑ=π/2, giving a finite limiting probability for contacts on the vanishingly small flat sides of the spherocylinder. Finally, we consider the variation of the average contact force as a function of location on the particle surface.

Densest versus jammed packings of bent-core trimers

We identify putatively maximally dense packings of tangent-sphere trimers with fixed bond angles (θ=θ_{0}), and contrast them to the disordered jammed states they form under quasistatic and dynamic athermal compression. Incommensurability of θ_{0} with three-dimensional (3D) close packing does not by itself inhibit formation of dense 3D crystals; all θ_{0} allow formation of crystals with ϕ_{max}(θ_{0})>0.97ϕ_{cp}. Trimers are always able to arrange into periodic structures composed of close-packed bilayers or trilayers of triangular-lattice planes, separated by "gap layers" that accommodate the incommensurability. All systems have ϕ_{J} significantly below the monomeric value, indicating that trimers' quenched bond-length and bond-angle constraints always act to promote jamming. ϕ_{J} varies strongly with θ_{0}; straight (θ_{0}=0) trimers minimize ϕ_{J} while closed (θ_{0}=120^{∘}) trimers maximize it. Marginally jammed states of trimers with lower ϕ_{J}(θ_{0}) exhibit quantifiably greater disorder, and the lower ϕ_{J} for small θ_{0} is apparently caused by trimers' decreasing effective configurational freedom as they approach linearity.

Jamming transition as a paradigm to understand the loss landscape of deep neural networks

Deep learning has been immensely successful at a variety of tasks, ranging from classification to artificial intelligence. Learning corresponds to fitting training data, which is implemented by descending a very high-dimensional loss function. Understanding under which conditions neural networks do not get stuck in poor minima of the loss, and how the landscape of that loss evolves as depth is increased, remains a challenge. Here we predict, and test empirically, an analogy between this landscape and the energy landscape of repulsive ellipses. We argue that in fully connected deep networks a phase transition delimits the over- and underparametrized regimes where fitting can or cannot be achieved. In the vicinity of this transition, properties of the curvature of the minima of the loss (the spectrum of the Hessian) are critical. This transition shares direct similarities with the jamming transition by which particles form a disordered solid as the density is increased, which also occurs in certain classes of computational optimization and learning problems such as the perceptron. Our analysis gives a simple explanation as to why poor minima of the loss cannot be encountered in the overparametrized regime. Interestingly, we observe that the ability of fully connected networks to fit random data is independent of their depth, an independence that appears to also hold for real data. We also study a quantity Δ which characterizes how well (Δ<0) or badly (Δ>0) a datum is learned. At the critical point it is power-law distributed on several decades, P_{+}(Δ)∼Δ^{θ} for Δ>0 and P_{-}(Δ)∼(-Δ)^{-γ} for Δ<0, with exponents that depend on the choice of activation function. This observation suggests that near the transition the loss landscape has a hierarchical structure and that the learning dynamics is prone to avalanche-like dynamics, with abrupt changes in the set of patterns that are learned.

Effects of spherical confinement and backbone stiffness on flexible polymer jamming

We use molecular simulations to study jamming of a crumpled bead-spring model polymer in a finite container and compare to jamming of repulsive spheres. After proper constraint counting, the onset of rigidity is seen to occur isostatically as in the case of repulsive spheres. Despite this commonality, the presence of the curved container wall and polymer backbone bonds introduce new mechanical properties. Notably, these include additional bands in the vibrational density of states that reflect the material structure as well as oscillations in local contact number and density near the wall but with lower amplitude for polymers. Polymers have fewer boundary contacts, and this low-density surface layer strongly reduces the global bulk modulus. We further show that bulk-modulus dependence on backbone stiffness can be described by a model of stiffnesses in series and discuss potential experimental and biological applications.

Biomechanical Feedback Strengthens Jammed Cellular Packings

Growth in confined spaces can drive cellular populations through a jamming transition from a fluidlike state to a solidlike state. Experiments have found that jammed budding yeast populations can build up extreme compressive pressures (over 1 MPa), which in turn feed back onto cellular physiology by slowing or even stalling cell growth. Using numerical simulations, we investigate how this feedback impacts the mechanical properties of model jammed cell populations. We find that feedback directs growth toward poorly coordinated regions, resulting in an excess number of cell-cell contacts that rigidify cell packings. Cell packings possess anomalously large shear and bulk moduli that depend sensitively on the strength of feedback. These results demonstrate that mechanical feedback on the single-cell level is a simple mechanism by which living systems may tune their population-level mechanical properties.

Machine learning characterization of structural defects in amorphous packings of dimers and ellipses

Structural defects within amorphous packings of symmetric particles can be characterized using a machine learning approach that incorporates structure functions of radial distances and angular arrangement. This yields a scalar field, softness, that correlates with the probability that a particle is about to be rearranged. However, when particle shapes are elongated, as in the case of dimers and ellipses, we find that the standard structure functions produce imprecise softness measurements. Moreover, ellipses exhibit deformation profiles in stark contrast to circular particles. In order to account for the effects of orientation and alignment, we introduce structure functions to recover the predictive performance of softness, as well as provide physical insight into local and extended dynamics. We study a model disordered solid, a bidisperse two-dimensional granular pillar, driven by uniaxial compression and composed entirely of monomers, dimers, or ellipses. We demonstrate how the computation of softness via a support vector machine extends to dimers and ellipses with the introduction of orientational structure functions. Then we highlight the spatial extent of rearrangements and defects, as well as their cross correlation, for each particle shape. Finally, we demonstrate how an additional machine learning algorithm, recursive feature elimination, provides an avenue to better understand how softness arises from particular structural aspects. We identify the most crucial structure functions in determining softness and discuss their physical implications.

The physics of jamming for granular materials: a review

Granular materials consist of macroscopic grains, interacting via contact forces, and unaffected by thermal fluctuations. They are one of a class systems that undergo jamming, i.e. a transition between fluid-like and disordered solid-like states. Roughly twenty years ago, proposals by Cates et al for the shear response of colloidal systems and by Liu and Nagel, for a universal jamming diagram in a parameter space of packing fraction, ϕ, shear stress, τ, and temperature, T raised key questions. Contemporaneously, experiments by Howell et al and numerical simulations by Radjai et al and by Luding et al helped provide a starting point to explore key insights into jamming for dry, cohesionless, granular materials. A recent experimental observation by Bi et al is that frictional granular materials have a a re-entrant region in their jamming diagram. In a range of ϕ, applying shear strain, γ, from an initially force/stress free state leads to fragile (in the sense of Cates et al), then anisotropic shear jammed states. Shear jamming at fixed ϕ is presumably conjugate to Reynolds dilatancy, involving dilation under shear against deformable boundaries. Numerical studies by Radjai and Roux showed that Reynolds dilatancy does not occur for frictionless systems. Recent numerical studies by several groups show that shear jamming occurs for finite, but not infinite, systems of frictionless grains. Shear jamming does not lead to known ordering in position space, but Sarkar et al showed that ordering occurs in a space of force tiles. Experimental studies seeking to understand random loose and random close packings (rlp and rcp) and dating back to Bernal have probed granular packings and their response to shear and intruder motion. These studies suggest that rlp's are anisotropic and shear-jammed-like, whereas rcp's are likely isotropically jammed states. Jammed states are inherently static, but the jamming diagram may provide a context for understanding rheology, i.e. dynamic shear in a variety of systems that include granular materials and suspensions.

Vibrational Modes and Dynamic Heterogeneity in a Near-Equilibrium 2D Glass of Colloidal Kites

Using video microscopy and particle-tracking techniques developed for dense Brownian systems of polygons, we study the structure-dynamics relationship in a near-equilibrium 2D glass consisting of anisotropic Penrose kite-shaped colloids. Detailed vibrational properties of kite glasses, both translational and rotational, are obtained using covariance matrix techniques. Different from other colloidal glasses of spheres and ellipsoids, the vibrational modes of kite glasses at low frequencies show a strong translational character with spatially localized rotational modes and extended translational modes. Low-frequency quasilocalized soft modes commonly found in sphere glasses are absent in the translational phonon modes of kite glasses. Soft modes are observed predominantly in the rotational vibrations and correlate well with the spatial distribution of Debye-Waller factors. The local structural entropy field shows a strong correlation with the observed dynamic heterogeneity.

Universality of jamming of nonspherical particles

Amorphous packings of nonspherical particles such as ellipsoids and spherocylinders are known to be hypostatic: The number of mechanical contacts between particles is smaller than the number of degrees of freedom, thus violating Maxwell's mechanical stability criterion. In this work, we propose a general theory of hypostatic amorphous packings and the associated jamming transition. First, we show that many systems fall into a same universality class. As an example, we explicitly map ellipsoids into a system of "breathing" particles. We show by using a marginal stability argument that in both cases jammed packings are hypostatic and that the critical exponents related to the contact number and the vibrational density of states are the same. Furthermore, we introduce a generalized perceptron model which can be solved analytically by the replica method. The analytical solution predicts critical exponents in the same hypostatic jamming universality class. Our analysis further reveals that the force and gap distributions of hypostatic jamming do not show power-law behavior, in marked contrast to the isostatic jamming of spherical particles. Finally, we confirm our theoretical predictions by numerical simulations.

Compression-driven jamming of athermal frictionless spherocylinders in two dimensions

We simulate numerically the compression-driven jamming of athermal, frictionless, soft-core spherocylinders in two dimensions, for a range of particle aspect ratios α. We find the critical packing fraction ϕ_{J}(α) for the jamming transition and the average number of contacts per particle z_{J}(α) at jamming. We find that both are nonmonotonic, with a peak at α≈1. We find that configurations at the compression-driven jamming point are always hypostatic for all α, with z_{J}<z_{iso}=2d_{f}=6 the isostatic value. We show that, for moderately elongated spherocylinders, there is no orientational ordering upon athermal compression through jamming. We analyze in detail the eigenmodes of the dynamical matrix close to the jamming point for a few different values of the aspect ratio, from nearly circular to moderately elongated. We find that there are low frequency bands containing N(z_{iso}-z_{J})/2 modes, such that the frequencies of these modes vanish as ϕ→ϕ_{J}. We consider the extended versus localized nature of these low frequency modes, and the extent to which they involve translational or rotational motion, and find many low frequency sliding modes where particles can move with little rotation. We highlight the importance of treating side-to-side contacts, along flat sides of the spherocylinder, properly for the correct determination of z_{J}. We note the singular nature of taking the α→0 limit. We discuss the similarities and differences with previous work on jammed ellipses and ellipsoids, to illustrate the effects that different particle shapes have on configurations at jamming.

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.

Anisotropic particles strengthen granular pillars under compression

We probe the effects of particle shape on the global and local behavior of a two-dimensional granular pillar, acting as a proxy for a disordered solid, under uniaxial compression. This geometry allows for direct measurement of global material response, as well as tracking of all individual particle trajectories. In general, drawing connections between local structure and local dynamics can be challenging in amorphous materials due to lower precision of atomic positions, so this study aims to elucidate such connections. We vary local interactions by using three different particle shapes: discrete circular grains (monomers), pairs of grains bonded together (dimers), and groups of three bonded in a triangle (trimers). We find that dimers substantially strengthen the pillar and the degree of this effect is determined by orientational order in the initial condition. In addition, while the three particle shapes form void regions at distinct rates, we find that anisotropies in the local amorphous structure remain robust through the definition of a metric that quantifies packing anisotropy. Finally, we highlight connections between local deformation rates and local structure.

Isostaticity and the solidification of semiflexible polymer melts

Using molecular dynamics simulations of a tangent-soft-sphere bead-spring polymer model, we examine the degree to which semiflexible polymer melts solidify at isostaticity. Flexible and stiff chains crystallize when they are isostatic as defined by appropriate degree-of-freedom-counting arguments. Semiflexible chains also solidify when isostatic if a generalized isostaticity criterion that accounts for the slow freezing out of configurational freedom as chain stiffness increases is employed. The configurational freedom associated with bond angles (θ) can be associated with the characteristic ratio C_∞ = (1 + 〈cos(θ)〉)/(1 - 〈cos(θ)〉). We find that the dependence of the average coordination number at solidification [Z(T_s)] on chains' characteristic ratio C_∞ has the same functional form [Z ≃ a - b ln(C_∞)] as the dependence of the average coordination number at jamming [Z(ϕ_J)] on C_∞ in athermal systems, suggesting that jamming-related phenomena play a significant role in thermal polymer solidification.

Slow creep in soft granular packings

Transient creep mechanisms in soft granular packings are studied numerically using a constant pressure and constant stress simulation method. Rapid compression followed by slow dilation is predicted on the basis of a logarithmic creep phenomenon. Characteristic scales of creep strain and time exhibit a power-law dependence on jamming pressure, and they diverge at the jamming point. Microscopic analysis indicates the existence of a correlation between rheology and nonaffine fluctuations. Localized regions of large strain appear during creep and grow in magnitude and size at short times. At long times, the spatial structure of highly correlated local deformation becomes time-invariant. Finally, a microscale connection between local rheology and local fluctuations is demonstrated in the form of a linear scaling between granular fluidity and nonaffine velocity.

Scaling Theory for the Frictionless Unjamming Transition

We develop a scaling theory of the unjamming transition of soft frictionless disks in two dimensions by defining local areas, which can be uniquely assigned to each contact. These serve to define local order parameters, whose distribution exhibits divergences as the unjamming transition is approached. We derive scaling forms for these divergences from a mean-field approach that treats the local areas as noninteracting entities, and demonstrate that these results agree remarkably well with numerical simulations. We find that the asymptotic behavior of the scaling functions arises from the geometrical structure of the packing while the overall scaling with the compression energy depends on the force law. We use the scaling forms of the distributions to determine the scaling of the total grain area A_{G} and the total number of contacts N_{C}.

Jamming of Semiflexible Polymers

We study jamming in model freely rotating polymers as a function of chain length N and bond angle θ_{0}. The volume fraction at jamming ϕ_{J}(θ_{0}) is minimal for rigid-rodlike chains (θ_{0}=0), and increases monotonically with increasing θ_{0}≤π/2. In contrast to flexible polymers, marginally jammed states of freely rotating polymers are highly hypostatic, even when bond and angle constraints are accounted for. Large-aspect-ratio (small θ_{0}) chains behave comparably to stiff fibers: resistance to large-scale bending plays a major role in their jamming phenomenology. Low-aspect-ratio (large θ_{0}) chains behave more like flexible polymers, but still jam at much lower densities due to the presence of frozen-in three-body correlations corresponding to the fixed bond angles. Long-chain systems jam at lower ϕ and are more hypostatic at jamming than short-chain systems. Implications of these findings for polymer solidification are discussed.

Elastic moduli and vibrational modes in jammed particulate packings

When we elastically impose a homogeneous, affine deformation on amorphous solids, they also undergo an inhomogeneous, nonaffine deformation, which can have a crucial impact on the overall elastic response. To correctly understand the elastic modulus M, it is therefore necessary to take into account not only the affine modulus M_{A}, but also the nonaffine modulus M_{N} that arises from the nonaffine deformation. In the present work, we study the bulk (M=K) and shear (M=G) moduli in static jammed particulate packings over a range of packing fractions φ. The affine M_{A} is determined essentially by the static structural arrangement of particles, whereas the nonaffine M_{N} is related to the vibrational eigenmodes. We elucidate the contribution of each vibrational mode to the nonaffine M_{N} through a modal decomposition of the displacement and force fields. In the vicinity of the (un)jamming transition φ_{c}, the vibrational density of states g(ω) shows a plateau in the intermediate-frequency regime above a characteristic frequency ω^{*}. We illustrate that this unusual feature apparent in g(ω) is reflected in the behavior of M_{N}: As φ→φ_{c}, where ω^{*}→0, those modes for ω<ω^{*} contribute less and less, while contributions from those for ω>ω^{*} approach a constant value which results in M_{N} to approach a critical value M_{Nc}, as M_{N}-M_{Nc}∼ω^{*}. At φ_{c} itself, the bulk modulus attains a finite value K_{c}=K_{Ac}-K_{Nc}>0, such that K_{Nc} has a value that remains below K_{Ac}. In contrast, for the critical shear modulus G_{c}, G_{Nc} and G_{Ac} approach the same value so that the total value becomes exactly zero, G_{c}=G_{Ac}-G_{Nc}=0. We explore what features of the configurational and vibrational properties cause such a distinction between K and G, allowing us to validate analytical expressions for their critical values.

Structural Transition in a Fluid of Spheroids: A Low-Density Vestige of Jamming

A thermodynamically equilibrated fluid of hard spheroids is a simple model of liquid matter. In this model, the coupling between the rotational degrees of freedom of the constituent particles and their translations may be switched off by a continuous deformation of a spheroid of aspect ratio t into a sphere (t=1). We demonstrate, by experiments, theory, and computer simulations, that dramatic nonanalytic changes in structure and thermodynamics of the fluids take place, as the coupling between rotations and translations is made to vanish. This nonanalyticity, reminiscent of a second-order liquid-liquid phase transition, is not a trivial consequence of the shape of an individual particle. Rather, free volume considerations relate the observed transition to a similar nonanalyticity at t=1 in structural properties of jammed granular ellipsoids. This observation suggests a deep connection to exist between the physics of jamming and the thermodynamics of simple fluids.

Structural characterization of the packings of granular regular polygons

By using a recently developed method for discrete modeling of nonspherical particles, we simulate the random packings of granular regular polygons with three to 11 edges under gravity. The effects of shape and friction on the packing structures are investigated by various structural parameters, including packing fraction, the radial distribution function, coordination number, Voronoi tessellation, and bond-orientational order. We find that packing fraction is generally higher for geometrically nonfrustrated regular polygons, and can be increased by the increase of edge number and decrease of friction. The changes of packing fraction are linked with those of the microstructures, such as the variations of the translational and orientational orders and local configurations. In particular, the free areas of Voronoi tessellations (which are related to local packing fractions) can be described by log-normal distributions for all polygons. The quantitative analyses establish a clearer picture for the packings of regular polygons.

Local origin of global contact numbers in frictional ellipsoid packings

In particulate soft matter systems the average number of contacts Z of a particle is an important predictor of the mechanical properties of the system. Using x-ray tomography, we analyze packings of frictional, oblate ellipsoids of various aspect ratios α, prepared at different global volume fractions ϕg. We find that Z is a monotonically increasing function of ϕg for all α. We demonstrate that this functional dependence can be explained by a local analysis where each particle is described by its local volume fraction ϕl computed from a Voronoi tessellation. Z can be expressed as an integral over all values of ϕl: Z(ϕg,α,X)=∫Zl(ϕl,α,X)P(ϕl|ϕg)dϕl. The local contact number function Zl(ϕl,α,X) describes the relevant physics in term of locally defined variables only, including possible higher order terms X. The conditional probability P(ϕl|ϕg) to find a specific value of ϕl given a global packing fraction ϕg is found to be independent of α and X. Our results demonstrate that for frictional particles a local approach is not only a theoretical requirement but also feasible.

Extraction of force-chain network architecture in granular materials using community detection

Force chains form heterogeneous physical structures that can constrain the mechanical stability and acoustic transmission of granular media. However, despite their relevance for predicting bulk properties of materials, there is no agreement on a quantitative description of force chains. Consequently, it is difficult to compare the force-chain structures in different materials or experimental conditions. To address this challenge, we treat granular materials as spatially-embedded networks in which the nodes (particles) are connected by weighted edges that represent contact forces. We use techniques from community detection, which is a type of clustering, to find sets of closely connected particles. By using a geographical null model that is constrained by the particles' contact network, we extract chain-like structures that are reminiscent of force chains. We propose three diagnostics to measure these chain-like structures, and we demonstrate the utility of these diagnostics for identifying and characterizing classes of force-chain network architectures in various materials. To illustrate our methods, we describe how force-chain architecture depends on pressure for two very different types of packings: (1) ones derived from laboratory experiments and (2) ones derived from idealized, numerically-generated frictionless packings. By resolving individual force chains, we quantify statistical properties of force-chain shape and strength, which are potentially crucial diagnostics of bulk properties (including material stability). These methods facilitate quantitative comparisons between different particulate systems, regardless of whether they are measured experimentally or numerically.

Relaxation dynamics in a binary hard-ellipse liquid

Structural relaxation in binary hard spherical particles has been shown recently to exhibit a wealth of remarkable features when size disparity or mixture composition is varied. In this paper, we test whether or not similar dynamical phenomena occur in glassy systems composed of binary hard ellipses. We demonstrate via event-driven molecular dynamics simulation that a binary hard-ellipse mixture with an aspect ratio of two and moderate size disparity displays characteristic glassy dynamics upon increasing density in both the translational and the rotational degrees of freedom. The rotational glass transition density is found to be close to the translational one for the binary mixtures investigated. More importantly, we assess the influence of size disparity and mixture composition on the relaxation dynamics. We find that an increase of size disparity leads, both translationally and rotationally, to a speed up of the long-time dynamics in the supercooled regime so that both the translational and the rotational glass transition shift to higher densities. By increasing the number concentration of the small particles, the time evolution of both translational and rotational relaxation dynamics at high densities displays two qualitatively different scenarios, i.e., both the initial and the final part of the structural relaxation slow down for small size disparity, while the short-time dynamics still slows down but the final decay speeds up in the binary mixture with large size disparity. These findings are reminiscent of those observed in binary hard spherical particles. Therefore, our results suggest a universal mechanism for the influence of size disparity and mixture composition on the structural relaxation in both isotropic and anisotropic particle systems.

Jamming in finite systems: stability, anisotropy, fluctuations, and scaling

Athermal packings of soft repulsive spheres exhibit a sharp jamming transition in the thermodynamic limit. Upon further compression, various structural and mechanical properties display clean power-law behavior over many decades in pressure. As with any phase transition, the rounding of such behavior in finite systems close to the transition plays an important role in understanding the nature of the transition itself. The situation for jamming is surprisingly rich: the assumption that jammed packings are isotropic is only strictly true in the large-size limit, and finite-size has a profound effect on the very meaning of jamming. Here, we provide a comprehensive numerical study of finite-size effects in sphere packings above the jamming transition, focusing on stability as well as the scaling of the contact number and the elastic response.

Effect of size polydispersity versus particle shape in dense granular media

We present a detailed analysis of the morphology of granular systems composed of frictionless pentagonal particles by varying systematically both the size span and particle shape irregularity, which represent two polydispersity parameters of the system. The microstructure is characterized in terms of various statistical descriptors such as global and local packing fractions, radial distribution functions, coordination number, and fraction of floating particles. We find that the packing fraction increases with the two parameters of polydispersity, but the effect of shape polydispersity for all the investigated structural properties is significant only at low size polydispersity where the positional and/or orientational ordering of the particles prevail. We focus in more detail on the class of side/side contacts, which is the interesting feature of our system as compared to a packing of disks. We show that the proportion of such contacts has weak dependence on the polydispersity parameters. The side- side contacts do not percolate but they define clusters of increasing size as a function of size polydispersity and decreasing size as a function of shape polydispersity. The clusters have anisotropic shapes but with a decreasing aspect ratio as polydispersity increases. This feature is argued to be a consequence of strong force chains (forces above the mean), which are mainly captured by side-side contacts. Finally, the force transmission is intrinsically multiscale, with a mean force increasing linearly with particle size.

Aspects of jamming in two-dimensional athermal frictionless systems

In this work we provide an overview of jamming transitions in two dimensional systems focusing on the limit of frictionless particle interactions in the absence of thermal fluctuations. We first discuss jamming in systems with short range repulsive interactions, where the onset of jamming occurs at a critical packing density and where certain quantities show a divergence indicative of critical behavior. We describe how aspects of the dynamics change as the jamming density is approached and how these dynamics can be explored using externally driven probes. Different particle shapes can produce jamming densities much lower than those observed for disk-shaped particles, and we show how jamming exhibits fragility for some shapes while for other shapes this is absent. Next we describe the effects of long range interactions and jamming behavior in systems such as charged colloids, vortices in type-II superconductors, and dislocations. We consider the effect of adding obstacles to frictionless jamming systems and discuss connections between this type of jamming and systems that exhibit depinning transitions. Finally, we discuss open questions such as whether the jamming transition in all these different systems can be described by the same or a small subset of universal behaviors, as well as future directions for studies of jamming transitions in two dimensional systems, such as jamming in self-driven or active matter systems.

Physics in ordered and disordered colloidal matter composed of poly(N-isopropylacrylamide) microgel particles

This review collects and describes experiments that employ colloidal suspensions to probe physics in ordered and disordered solids and related complex fluids. The unifying feature of this body of work is its clever usage of poly(N-isopropylacrylamide) (PNIPAM) microgel particles. These temperature-sensitive colloidal particles provide experimenters with a 'knob' for in situ control of particle size, particle interaction and particle packing fraction that, in turn, influence the structural and dynamical behavior of the complex fluids and solids. A brief summary of PNIPAM particle synthesis and properties is given, followed by a synopsis of current activity in the field. The latter discussion describes a variety of soft matter investigations including those that explore formation and melting of crystals and clusters, and those that probe structure, rearrangement and rheology of disordered (jammed/glassy) and partially ordered matter. The review, therefore, provides a snapshot of a broad range of physics phenomenology which benefits from the unique properties of responsive microgel particles.

Dynamics of sheared ellipses and circular disks: effects of particle shape

Much recent effort has focused on glassy and jamming properties of spherical particles. Very little is known about such phenomena for nonspherical particles, and we take a first step by studying ellipses. We find important differences between the dynamical and structural properties of disks and two-dimensional ellipses subject to continuous Couette shear. In particular, ellipses show slow dynamical evolution, without a counterpart in disks, in the mean velocity, local density, orientational order, and local stress. Starting from an unjammed state, ellipses can first jam under shear, and then slowly unjam. The slow unjamming process is understood as a result of gradual changes in their orientations, leading to a denser packing. For disks, the rotation of particles only contributes to the relaxation of frictional forces, and hence, does not significantly cause structural changes. For the shear-jammed states, the global building up and relaxation of stress, which occurs in the form of stress avalanches, is qualitatively different for disks and ellipses, and is manifested by different forms of rate dependence for ellipses versus disks. Unlike the weak rate dependence typical for many granular systems, ellipses show power-law dependence on the shearing rate Ω.

Variable-cell method for stress-controlled jamming of athermal, frictionless grains

A method is introduced to simulate jamming of polyhedral grains under controlled stress that incorporates global degrees of freedom through the metric tensor of a periodic cell containing grains. Jamming under hydrostatic (isotropic) stress and athermal conditions leads to a precise definition of the ideal jamming point at zero shear stress. The structures of tetrahedra jammed hydrostatically exhibit less translational order and lower jamming-point density than previously described maximally random jammed hard tetrahedra. Under the same conditions, cubes jam with negligible nematic order. Grains with octahedral symmetry having s>0.5 (where s interpolates from octahedra [s=0] to cubes [s=1]) jam with an abundance of face-face contacts in the absence of nematic order. For sufficiently large face-face contact number, percolating clusters form that span the entire simulation box. The response of hydrostatically jammed tetrahedra and cubes to shear-stress perturbation is also demonstrated with the variable-cell method.