Search

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 $\langle G \rangle$ scales with pressure $P$ as $\langle G \rangle \sim P^{\beta}$, with $\beta \approx 0.75$ for jammed packings of deformable particles with zero bending energy. In contrast, $\beta \approx 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.

Soft, amorphous solids such as tissues, foams, and emulsions are composed of deformable particles. However, the effect of single-particle deformability on the collective behavior of soft solids is still poorly understood. We perform numerical simulations of two-dimensional jammed packings of explicitly deformable particles to study the mechanical response of model soft solids. We find that jammed packings of deformable particles with excess shape degrees of freedom possess low-frequency quartic vibrational modes that stabilize the packings even though they possess fewer interparticle contacts than the nominal isostatic value. Adding intra-particle constraints can rigidify the particles, but these particles undergo a buckling transition and gain an effective shape degree of freedom when their preferred perimeter is above a threshold value. We find that the mechanical response of jammed packings of deformable particles with shape degrees of freedom differs significantly from that of jammed packings of rigid-shape particles, which emphasizes the importance of particle deformability in modelling soft solids., Comment: 12 pages total, 6 main figures, 4 appendices, 5 appendix figures

The mechanical response of packings of purely repulsive, spherical particles to athermal, quasistatic simple shear near jamming onset is highly nonlinear. Previous studies have shown that, at small pressure $p$, the ensemble-averaged static shear modulus $\langle G-G_0 \rangle$ scales with $p^\alpha$, where $\alpha \approx 1$, but above a characteristic pressure $p^{**}$, $\langle G-G_0 \rangle \sim p^\beta$, where $\beta \approx 0.5$. However, we find that the shear modulus $G^i$ for an individual packing typically decreases linearly with $p$ along a geometrical family where the contact network does not change. We resolve this discrepancy by showing that, while the shear modulus does decrease linearly within geometrical families, $\langle G \rangle$ also depends on a contribution from discontinuous jumps in $\langle G \rangle$ that occur at the transitions between geometrical families. For $p > p^{**}$, geometrical-family and rearrangement contributions to $\langle G \rangle$ are of opposite signs and remain comparable for all system sizes. $\langle G \rangle$ can be described by a scaling function that smoothly transitions between the two power-law exponents $\alpha$ and $\beta$. We also demonstrate the phenomenon of {\it compression unjamming}, where a jammed packing can unjam via isotropic compression., Comment: 6 pages, 4 pages

We perform computational studies of jammed particle packings in two dimensions undergoing isotropic compression using the well-characterized soft particle (SP) model and the deformable particle (DP) model that we developed for compressed bubbles and emulsions. In the SP model, circular particles are allowed to overlap, generating purely repulsive forces. In the DP model, particles minimize their perimeter, while deforming at fixed area to avoid overlap during compression. We directly compare the structural and mechanical properties of jammed particle packings generated using the SP and DP models as a function of the true packing fraction $\rho$, instead of the reduced number density $\phi$. We show that near jamming onset the excess contact number $\Delta z=z-z_J$ and shear modulus ${\cal G}$ scale as $\Delta \rho^{0.5}$ in the large system limit for both the SP and DP models, where $\Delta \rho = \rho-\rho_J$ and $z_J \approx 4$ and $\rho_J \approx 0.842$ are the values at jamming onset. $\Delta z$ and ${\cal G}$ for the SP and DP models begin to differ for $\rho \gtrsim 0.88$. In this regime, $\Delta z \sim {\cal G}$ can be described by a sum of two power-laws in $\Delta \rho$, i.e. $\Delta z \sim {\cal G} \sim C_0\Delta \rho^{0.5} +C_1\Delta \rho^{1.0}$ to lowest order. We show that the ratio $C_1/C_0$ is much larger for the DP model compared to to that for the SP model. We also characterize the void space in jammed packings as a function of $\rho$. We find that, unlike the SP model, the DP model is able to describe the formation of Plateau borders as the system approaches $\rho = 1$. We further show that the results for $z$ and the shape factor ${\cal A}$ versus $\rho$ for the DP model agree with recent experimental studies of compressed foams and emulsions., Comment: 12 pages, 11 figures, 9 main figures, and 2 supplemental figures

There are two main classes of physics-based models for two-dimensional cellular materials: packings of repulsive disks and the vertex model. These models have several disadvantages. For example, disk interactions are typically a function of particle overlap, yet the model assumes that the disks remain circular during overlap. The shapes of the cells can vary in the vertex model, however, the packing fraction is fixed at $\phi=1$. Here, we describe the deformable particle model (DPM), where each particle is a polygon composed of a large number of vertices. The total energy includes three terms: two quadratic terms to penalize deviations from the preferred particle area $a_0$ and perimeter $p_0$ and a repulsive interaction between DPM polygons that penalizes overlaps. We performed simulations to study the onset of jamming in packings of DPM polygons as a function of asphericity, ${\cal A} = p_0^2/4\pi a_0$. We show that the packing fraction at jamming onset $\phi_J({\cal A})$ grows with increasing ${\cal A}$, reaching confluence at ${\cal A} \approx 1.16$. ${\cal A}^*$ corresponds to the value at which DPM polygons completely fill the cells obtained from a surface-Voronoi tessellation. Further, we show that DPM polygons develop invaginations for ${\cal A} > {\cal A}^*$ with excess perimeter that grows linearly with ${\cal A}-{\cal A}^*$. We confirm that packings of DPM polygons are solid-like over the full range of ${\cal A}$ by showing that the shear modulus is nonzero., Comment: 6 pages, 5 figures