• Abstract and Figures
• Public Full-text

The structural, vibrational, and mechanical properties of jammed packings of deformable particles in three dimensions Dong Wang,1 John D. Treado,1, 2 Arman Boromand,1 Blake Norwick,3 Michael P. Murrell,4, 3, 5 Mark D. Shattuck,6 and Corey S. O’Hern1, 2, 3, 7, ∗ 1Department of Mechanical Engineering & Materials Science, Yale University, New Haven, Connecticut 06520, USA 2Integrated Graduate Program in Physical and Engineering Biology, Yale University, New Haven, Connecticut 06520, USA 3Department of Physics, Yale University, New Haven, Connecticut 06520, USA 4Department of Biomedical Engineering, Yale University, New Haven, Connecticut 06520, USA 5Systems Biology Institute, Yale University, West Haven, Connecticut 06516, USA 6Benjamin Levich Institute and Physics Department, The City College of New York, New York, New York 10031, USA 7Department of Applied Physics, Yale University, New Haven, Connecticut 06520, USA (Dated: August 24, 2021) We investigate the structural, vibrational, and mechanical properties of jammed packings of de- formable 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 hGi scales with pressure P as hGi∼ 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, spher- ical particles with fixed shape. These studies underscore the importance of incorporating particle deformability and shape change when modeling the properties of jammed soft materials. I. INTRODUCTION model in 2D that combines optimal features of both classes of models for soft particles [30, 31]. The DP Numerous physical systems are composed of discrete, model treats each particle as a discrete object, and thus soft particles that can change shape under applied stress. in contrast to vertex- or Voronoi-based models, the DP Examples include collections of emulsion droplets [1, 2], model can be used to study systems over a wide range colloids [3], bubbles [4], and hydrogel particles [5, 6]. of packing fractions—from isolated particles to confluent These systems display complex, spatio-temporal response systems. In addition, the DP model considers shape- to applied deformations, including shear jamming [7–11], energy functions for each particle individually (through the shape parameter = p2/4πa, where p and a are the shear banding [12, 13], aging [14–17], and memory for- A mation [18–20]. perimeter and area of the particle, and the bending en- Many of the physics-based, theoretical models that are ergy), and thus the DP model provides control over the used to investigate the mechanical and vibrational re- shape of each deformable particle separately. In previ- ous studies, we investigated the mechanical and vibra- arXiv:2108.09616v1 [cond-mat.soft] 22 Aug 2021 sponse of soft materials fall into one of two classes: 1) “soft-particle” models for which the interparticle forces tional properties of jammed packings of DP particles in are generated by overlaps between discrete particles of 2D with and without bending energy[32]. We showed fixed shape [21–23] and 2) vertex- or Voronoi-based mod- that packings of DP particles without bending energy are els [24–29] that treat the system as space-filling polygons hypostatic (with fewer contacts than the isostatic value) in two dimensions (2D) or polyhedra in 3D, with interpar- over the full range of shape parameters, and that the ticle forces determined by shape-energy functions written missing contacts are stabilized by low-frequency, quartic in terms of the vertices of all polygons in 2D or polyhedra modes of the dynamical matrix. When perturbing the in 3D. system along quartic modes with amplitude δ, the en- 4 2 Recently, we introduced the deformable particle (DP) ergy of the system increases as δ , rather than δ as for non-quartic modes. Particles with non-zero bending en- ergy undergo a buckling transition when > ∗, above which the minimal energy shape is not a regularA A polygon ∗ [email protected] and ∗ increases with the bending stiffness. Packings of A 2 unbuckled particles with < ∗ are isostatic with no particles, as well two types of point processes. In Appen- quartic modes. In contrast,A packingsA of buckled parti- dices C and D, we show the influence of adding thermal cles with > ∗ possess quartic modes, but we showed fluctuations to the compression protocol for generating that it is difficultA A to determine how many degrees of free- jammed packings of deformable particles on the proper- dom are associated with each buckled particle [33]. The ties of jammed packings of deformable particles in two shape degrees of freedom contribute significantly to the and three dimensions. vibrational response for packings of DP particles without bending energy, which gives rise to power-law scaling of the ensemble-averaged shear modulus with pressure that II. METHODS differs from that for jammed packings of particles with fixed shape. We model deformable particles in 3D as surface- In this article, we develop the DP model in 3D, which triangulated polyhedra with Nv vertices as shown in considers particles as surface-triangulated polyhedra, and Fig. 1. The vertices are connected via Delaunay trian- then use it to investigate the structural, mechanical, and gulation, resulting in Nf triangles and Ne edges on the vibrational response of jammed packings of deformable surface of each polyhedron. We characterize the shape of particles in 3D. The 3D DP model allows us to deter- 3D deformable particles using the non-dimensional shape mine whether the structural, vibrational, and mechanical parameter (or asphericity) = s3/2/(6√πv), where s A properties of jammed packings of deformable particles in and v are the total surface area and volume of the par- 3D are similar to those in 2D, which is important for ticle, respectively. = 1 when the particle is a sphere, A classifying the critical behavior [34] of the jamming tran- and > 1 for any non-spherical shape. A sition for deformable particle packings. We will show that The total potential energy U for a collection of N de- many of the results for jammed packings of deformable formable particles in 3D obeys the following: particles are similar in 2D and 3D. For example, packings N 2 of deformable particles with no bending energy possess N ǫ v 2 N f ǫ a U = v n 1 + a nf 1 low-frequency, quartic modes, whose number matches the 2 v − 2 f − n=1 0 n=1 a0n ! number of missing contacts from simple contact counting. X X fX=1 (1) N N Also, the pressure-dependent mechanical response varies e ǫ + b (θ )2 + U int, strongly with the particle deformability; the ensemble- 2 ne averaged shear modulus scales with pressure as G P β n=1 e=1 h i∼ X X with β 0.75 for truly deformable particles, whereas f β 0.5 for≈ particles with non-zero bending energy. How- where v0 is the preferred volume for each particle and a0n ever,≈ in contrast to the results for 2D, we show that for is the preferred area for the fth triangle on the surface all non-zero values of the bending energy (i.e. both un- of particle n. The bending angle θne is the angle between buckled and buckled particles), DP packings in 3D are the two unit normals to the triangles that share the eth isostatic at jamming onset and do not possess quartic edge on particle n. The three coefficients ǫv, ǫa, and ǫb modes. control fluctuations in the particle volume, surface trian- gle area, and curvature, respectively. To prevent overlap The remainder of the article is organized as follows. between deformable particles, we include the purely re- In Sec. II, we describe the shape-energy function for the pulsive, linear spring interaction potential between over- DP model in 3D and the computational methods used lapping spherical vertices on neighboring particles: to generate jammed packings of deformable particles in 3D and to calculate the dynamical matrix, density of vi- N N Nv Nv 2 ǫ r r brational modes, stress tensor, and shear modulus for U int = c 1 ni,mj Θ 1 ni,mj , 2 − σ − σ these packings. In Sec. III, we discuss the results in- n=1 m>n i=1 j=1 ni,mj ni,mj X X X X cluding calculations of the vibrational modes for a single (2) deformable particle (Sec. IIIA), and the packing frac- where rni,mj is the distance between the centers of spher- tion and contact number at jamming onset (Sec. IIIB), ical vertices i and j on separate particles n and m and the density of vibrational modes (Sec. IIIC), the con- σni,mj = (σni +σmj )/2 is the average diameter of the two tribution of the shape degrees of freedom to the vibra- vertices. The Heaviside step function Θ( ) enforces that tional modes (Sec.

Load more

  16

Copy Link