Pressure Dependent Shear Response of Jammed Packings of Frictionless Spherical Particles

Unified Understanding of Nonlinear Rheology near the Jamming Transition Point

When slowly sheared, jammed packings respond elastically before yielding. This linear elastic regime becomes progressively narrower as the jamming transition point is approached, and rich nonlinear rheologies such as shear softening and hardening or melting emerge. However, the physical mechanism of these nonlinear rheologies remains elusive. To clarify this, we numerically study jammed packings of athermal frictionless soft particles under quasistatic shear γ. We find the universal scaling behavior for the ratio of the shear stress σ and the pressure P, independent of the preparation protocol of the initial configurations. In particular, we reveal shear softening σ/P∼γ^{1/2} over an unprecedentedly wide range of strain up to the yielding point, which a simple scaling argument can rationalize.

Protein folding as a jamming transition

Proteins fold to a specific functional conformation with a densely packed hydrophobic core that controls their stability. We develop a geometric, yet all-atom model for proteins that explains the universal core packing fraction of ϕ c = 0.55 found in experimental measurements. We show that as the hydrophobic interactions increase relative to the temperature, a novel jamming transition occurs when the core packing fraction exceeds ϕ c . The model also recapitulates the global structure of proteins since it can accurately refold to native-like structures from partially unfolded states.

How P. aeruginosa cells with diverse stator composition collectively swarm

Swarming is a macroscopic phenomenon in which surface bacteria organize into a motile population. The flagellar motor that drives swarming in Pseudomonas aeruginosa is powered by stators MotAB and MotCD. Deletion of the MotCD stator eliminates swarming, whereas deletion of the MotAB stator enhances swarming. Interestingly, we measured a strongly asymmetric stator availability in the wild-type (WT) strain, with MotAB stators produced at an approximately 40-fold higher level than MotCD stators. However, utilization of MotCD stators in free swimming cells requires higher liquid viscosities, while MotAB stators are readily utilized at low viscosities. Importantly, we find that cells with MotCD stators are ~10× more likely to have an active motor compared to cells uses the MotAB stators. The spectrum of motility intermittency can either cooperatively shut down or promote flagellum motility in WT populations. In P. aeruginosa, transition from a static solid-like biofilm to a dynamic liquid-like swarm is not achieved at a single critical value of flagellum torque or stator fraction but is collectively controlled by diverse combinations of flagellum activities and motor intermittencies via dynamic stator utilization. Experimental and computational results indicate that the initiation or arrest of flagellum-driven swarming motility does not occur from individual fitness or motility performance but rather related to concepts from the "jamming transition" in active granular matter.IMPORTANCEIt is now known that there exist multifactorial influences on swarming motility for P. aeruginosa, but it is not clear precisely why stator selection in the flagellum motor is so important. We show differential production and utilization of the stators. Moreover, we find the unanticipated result that the two motor configurations have significantly different motor intermittencies: the fraction of flagellum-active cells in a population on average with MotCD is active ~10× more often than with MotAB. What emerges from this complex landscape of stator utilization and resultant motor output is an intrinsically heterogeneous population of motile cells. We show how consequences of stator recruitment led to swarming motility and how the stators potentially relate to surface sensing circuitry.

Vibrational spectrum of Granular packings with random matrices

The vibrational spectrum of granular packings can be used as a signature of the jamming transition, with the density of states at zero frequency becoming nonzero at the transition. It has been proposed previously that the vibrational spectrum of granular packings can be approximately obtained from random matrix theory. Here, we show that the autocorrelation function of the density of states shows good agreement between dynamical numerical simulations of frictionless bead packs near the jamming point and the analytic predictions of the Laguerre orthogonal ensemble of random matrices; there is clear disagreement with the Gaussian orthogonal ensemble, establishing that the Laguerre ensemble correctly reproduces the universal statistical properties of jammed granular matter and excluding the Gaussian orthogonal ensemble. We also present a random lattice model which is a physically motivated variant of the random matrix ensemble. Numerical calculations reveal that this model reproduces the known features of the vibrational density of states of frictionless granular matter, while also retaining the correlation structure seen in the Laguerre random matrix theory.

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.

Local and global measures of the shear moduli of jammed disk packings

Strain-controlled isotropic compression gives rise to jammed packings of repulsive, frictionless disks with either positive or negative global shear moduli. We carry out computational studies to understand the contributions of the negative shear moduli to the mechanical response of jammed disk packings. We first decompose the ensemble-averaged, global shear modulus as 〈G〉=(1-F_{-})〈G_{+}〉+F_{-}〈G_{-}〉, where F_{-} is the fraction of jammed packings with negative shear moduli and 〈G_{+}〉 and 〈G_{-}〉 are the average values from packings with positive and negative moduli, respectively. We show that 〈G_{+}〉 and 〈|G_{-}|〉 obey different power-law scaling relations above and below pN^{2}∼1. For pN^{2}>1, both 〈G_{+}〉N and 〈|G_{-}|〉N∼(pN^{2})^{β}, where β∼0.5 for repulsive linear spring interactions. Despite this, 〈G〉N∼(pN^{2})^{β^{'}} with β^{'}≳0.5 due to the contributions from packings with negative shear moduli. We show further that the probability distribution of global shear moduli P(G) collapses at fixed pN^{2} and different values of p and N. We calculate analytically that P(G) is a Γ distribution in the pN^{2}≪1 limit. As pN^{2} increases, the skewness of P(G) decreases and P(G) becomes a skew-normal distribution with negative skewness in the pN^{2}≫1 limit. We also partition jammed disk packings into subsystems using Delaunay triangulation of the disk centers to calculate local shear moduli. We show that the local shear moduli defined from groups of adjacent triangles can be negative even when G>0. The spatial correlation function of local shear moduli C(r[over ⃗]) displays weak correlations for pn_{sub}^{2}<10^{-2}, where n_{sub} is the number of particles within each subsystem. However, C(r[over ⃗]) begins to develop long-ranged spatial correlations with fourfold angular symmetry for pn_{sub}^{2}≳10^{-2}.

How individual P. aeruginosa cells with diverse stator distributions collectively form a heterogeneous macroscopic swarming population

Swarming is a macroscopic phenomenon in which surface bacteria organize into a motile population. The flagellar motor that drives swarming in Pseudomonas aeruginosa is powered by stators MotAB and MotCD. Deletion of the MotCD stator eliminates swarming, whereas deletion of the MotAB stator enhances swarming. Interestingly, we measured a strongly asymmetric stator availability in the WT strain, with MotAB stators produced ∼40-fold more than MotCD stators. However, recruitment of MotCD stators in free swimming cells requires higher liquid viscosities, while MotAB stators are readily recruited at low viscosities. Importantly, we find that cells with MotCD stators are ∼10x more likely to have an active motor compared to cells without, so wild-type, WT, populations are intrinsically heterogeneous and not reducible to MotAB-dominant or MotCD-dominant behavior. The spectrum of motility intermittency can either cooperatively shut down or promote flagellum motility in WT populations. In P. aeruginosa , transition from a static solid-like biofilm to a dynamic liquid-like swarm is not achieved at a single critical value of flagellum torque or stator fraction but is collectively controlled by diverse combinations of flagellum activities and motor intermittencies via dynamic stator recruitment. Experimental and computational results indicate that the initiation or arrest of flagellum-driven swarming motility does not occur from individual fitness or motility performance but rather related to concepts from the 'jamming transition' in active granular matter.

Importance

After extensive study, it is now known that there exist multifactorial influences on swarming motility in P. aeruginosa , but it is not clear precisely why stator selection in the flagellum motor is so important or how this process is collectively initiated or arrested. Here, we show that for P. aeruginosa PA14, MotAB stators are produced ∼40-fold more than MotCD stators, but recruitment of MotCD over MotAB stators requires higher liquid viscosities. Moreover, we find the unanticipated result that the two motor configurations have significantly different motor intermittencies, the fraction of flagellum-active cells in a population on average, with MotCD active ∼10x more often than MotAB. What emerges from this complex landscape of stator recruitment and resultant motor output is an intrinsically heterogeneous population of motile cells. We show how consequences of stator recruitment led to swarming motility, and how they potentially relate to surface sensing circuitry.

Shear-Driven Solidification and Nonlinear Elasticity in Epithelial Tissues

Biological processes, from morphogenesis to tumor invasion, spontaneously generate shear stresses inside living tissue. The mechanisms that govern the transmission of mechanical forces in epithelia and the collective response of the tissue to bulk shear deformations remain, however, poorly understood. Using a minimal cell-based computational model, we investigate the constitutive relation of confluent tissues under simple shear deformation. We show that an initially undeformed fluidlike tissue acquires finite rigidity above a critical applied strain. This is akin to the shear-driven rigidity observed in other soft matter systems. Interestingly, shear-driven rigidity can be understood by a critical scaling analysis in the vicinity of the second order critical point that governs the liquid-solid transition of the undeformed system. We further show that a solidlike tissue responds linearly only to small strains and but then switches to a nonlinear response at larger stains, with substantial stiffening. Finally, we propose a mean-field formulation for cells under shear that offers a simple physical explanation of shear-driven rigidity and nonlinear response in a tissue.

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.

Disorder perturbation expansion for athermal crystals

We introduce a perturbation expansion for athermal systems that allows an exact determination of displacement fields away from the crystalline state as a response to disorder. We show that the displacement fields in energy-minimized configurations of particles interacting through central potentials with microscopic disorder can be obtained as a series expansion in the strength of the disorder. We introduce a hierarchy of force-balance equations that allows an order-by-order determination of the displacement fields, with the solutions at lower orders providing sources for the higher-order solutions. This allows the simultaneous force-balance equations to be solved, within a hierarchical perturbation expansion to arbitrary accuracy. We present exact results for an isotropic defect introduced into the crystalline ground state at linear order and second order in our expansion. We show that the displacement fields produced by the defect display interesting self-similar properties at every order. We derive a |δr|∼1/r and |δf|∼1/r^{2} decay for the displacement fields and excess interparticle forces at large distances r away from the defect. Finally, we derive nonlinear corrections introduced by the interactions between defects at second order in our expansion. We verify our exact results with displacement fields obtained from energy-minimized configurations of soft disks.

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)}.

Bulk modulus of soft particle assemblies under compression

Using a numerical approach based on the coupling of the discrete and finite element methods, we explore the variation of the bulk modulus K of soft particle assemblies undergoing isotropic compression. As the assemblies densify under pressure-controlled boundary conditions, we show that the non-linearities of K rapidly deviate from predictions standing on a small-strain framework or the, so-called, Equivalent Medium Theory (EMT). Using the granular stress tensor and extracting the bulk properties of single representative grains under compression, we propose a model to predict the evolution of K as a function of the sample’s solid fraction and a reference state as the applied pressure P→0. The model closely reproduces the trends observed in our numerical experiments confirming the behavior scalability of soft particle assemblies from the individual particle scale. Finally, we present the effect of the interparticle friction on K’s evolution and how our model easily adapts to such a mechanical constraint.

Contact network changes in ordered and disordered disk packings

We investigate the mechanical response of packings of purely repulsive, frictionless disks to quasistatic deformations. The deformations include simple shear strain at constant packing fraction and at constant pressure, "polydispersity" strain (in which we change the particle size distribution) at constant packing fraction and at constant pressure, and isotropic compression. For each deformation, we show that there are two classes of changes in the interparticle contact networks: jump changes and point changes. Jump changes occur when a contact network becomes mechanically unstable, particles "rearrange", and the potential energy (when the strain is applied at constant packing fraction) or enthalpy (when the strain is applied at constant pressure) and all derivatives are discontinuous. During point changes, a single contact is either added to or removed from the contact network. For repulsive linear spring interactions, second- and higher-order derivatives of the potential energy/enthalpy are discontinuous at a point change, while for Hertzian interactions, third- and higher-order derivatives of the potential energy/enthalpy are discontinuous. We illustrate the importance of point changes by studying the transition from a hexagonal crystal to a disordered crystal induced by applying polydispersity strain. During this transition, the system only undergoes point changes, with no jump changes. We emphasize that one must understand point changes, as well as jump changes, to predict the mechanical properties of jammed packings.

Jamming Below Upper Critical Dimension

Extensive numerical simulations in the past decades proved that the critical exponents of the jamming of frictionless spherical particles are the same in two and three dimensions. This implies that the upper critical dimension is d_{u}=2 or lower. In this Letter, we study the jamming transition below the upper critical dimension. We investigate a quasi-one-dimensional system: disks confined in a narrow channel. We show that the system is isostatic at the jamming transition point as in the case of standard jamming transition of the bulk systems in two and three dimensions. Nevertheless, the scaling of the excess contact number shows the linear scaling. Furthermore, the gap distribution remains finite even at the jamming transition point. These results are qualitatively different from those of the bulk systems in two and three dimensions.