The jamming transition of soft particles characterized by narrow size distributions has been well studied by physicists. However, polydispersed systems are more relevant to engineering, and the influence of polydispersity on jamming phenomena is still unexplored. Here, we numerically investigate jamming transitions of polydispersed soft particles in two dimensions. We find that polydispersity strongly influences contact forces, local coordination, and the jamming transition density. In contrast, the critical scaling of pressure and elastic moduli is not affected by the particle size distribution. Consistent with this observation, we find that the vibrational density of states is also insensitive to the polydispersity. Our results suggest that, regardless of particle size distributions, both mechanical and vibrational properties of soft particle packings near jamming are governed by the distance to jamming.@en

We numerically investigate stress relaxation in soft athermal disks to reveal critical slowing down when the system approaches the jamming point. The exponents describing the divergence of the relaxation time differ dramatically depending on whether the transition is approached from the jammed or unjammed phase. This contrasts sharply with conventional dynamic critical scaling scenarios, where a single exponent characterizes both sides. We explain this surprising difference in terms of the vibrational density of states, which is a key ingredient of linear viscoelastic theory. The vibrational density of states exhibits an extra slow mode that emerges below jamming, which we utilize to demonstrate the anomalous exponent below jamming.

@en

We numerically investigate nonlocal effects on inhomogeneous flows of soft athermal disks close to but below their jamming transition. We employ molecular dynamics to simulate Kolmogorov flows, in which a sinusoidal flow profile with fixed wave number is externally imposed, resulting in a spatially inhomogeneous shear rate. We find that the resulting rheology is strongly wave-number-dependent, and that particle migration, while present, is not sufficient to describe the resulting stress profiles within a conventional local model. We show that, instead, stress profiles can be captured with nonlocal constitutive relations that account for gradients to fourth order. Unlike nonlocal flow in yield stress fluids, we find no evidence of a diverging length scale.

@en