No Arabic abstract
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.
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.
At low volume fraction, disordered arrangements of frictionless spheres are found in un--jammed states unable to support applied stresses, while at high volume fraction they are found in jammed states with mechanical strength. Here we show, focusing on the hard sphere zero pressure limit, that the transition between un-jammed and jammed states does not occur at a single value of the volume fraction, but in a whole volume fraction range. This result is obtained via the direct numerical construction of disordered jammed states with a volume fraction varying between two limits, $0.636$ and $0.646$. We identify these limits with the random loose packing volume fraction $rl$ and the random close packing volume fraction $rc$ of frictionless spheres, respectively.
We focus on the response of mechanically stable (MS) packings of frictionless, bidisperse disks to thermal fluctuations, with the aim of quantifying how nonlinearities affect system properties at finite temperature. Packings of disks with purely repulsive contact interactions possess two main types of nonlinearities, one from the form of the interaction potential and one from the breaking (or forming) of interparticle contacts. To identify the temperature regime at which the contact-breaking nonlinearities begin to contribute, we first calculated the minimum temperatures $T_{cb}$ required to break a single contact in the MS packing for both single and multiple eigenmode perturbations of the $T=0$ MS packing. We then studied deviations in the constant volume specific heat $C_V$ and deviations of the average disk positions $Delta r$ from their $T=0$ values in the temperature regime $T_{cb} < T < T_{r}$, where $T_r$ is the temperature beyond which the system samples the basin of a new MS packing. We find that the deviation in the specific heat per particle $Delta {overline C}_V^0/{overline C}_V^0$ relative to the zero temperature value ${overline C}_V^0$ can grow rapidly above $T_{cb}$, however, the deviation $Delta {overline C}_V^0/{overline C}_V^0$ decreases as $N^{-1}$ with increasing system size. To characterize the relative strength of contact-breaking versus form nonlinearities, we measured the ratio of the average position deviations $Delta r^{ss}/Delta r^{ds}$ for single- and double-sided linear and nonlinear spring interactions. We find that $Delta r^{ss}/Delta r^{ds} > 100$ for linear spring interactions and is independent of system size.
We demonstrate the existence of unconventional rheological and memory properties in systems of soft-deformable particles whose energy depends on their shape, via numerical simulations. At large strains, these systems experience an unconventional shear weakening transition characterized by an increase in the mechanical energy and a drastic drop in shear stress, which stems from the emergence of short-ranged tetratic order. In these weakened states, the contact network evolves reversibly under strain reversal, keeping memory of its initial state, while the microscopic dynamics is irreversible.
We investigate the mechanical response of jammed packings of circulo-lines, interacting via purely repulsive, linear spring forces, as a function of pressure $P$ during athermal, quasistatic isotropic compression. Prior work has shown that the ensemble-averaged shear modulus for jammed disk packings scales as a power-law, $langle G(P) rangle sim P^{beta}$, with $beta sim 0.5$, over a wide range of pressure. For packings of circulo-lines, we also find robust power-law scaling of $langle G(P)rangle$ over the same range of pressure for aspect ratios ${cal R} gtrsim 1.2$. However, the power-law scaling exponent $beta sim 0.8$-$0.9$ is much larger than that for jammed disk packings. To understand the origin of this behavior, we decompose $langle Grangle$ into separate contributions from geometrical families, $G_f$, and from changes in the interparticle contact network, $G_r$, such that $langle G rangle = langle G_frangle + langle G_r rangle$. We show that the shear modulus for low-pressure geometrical families for jammed packings of circulo-lines can both increase {it 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 $langle G_f rangle$ 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.