No Arabic abstract
Under shear, a system of particles changes its contact network and becomes unstable as it transitions between mechanically stable states. For hard spheres at zero pressure, contact breaking events necessarily generate an instability, but this is not the case at finite pressure, where we identify two types of contact changes: network events that do not correspond to instabilities and rearrangement events that do. The relative fraction of such events is constant as a function of system size, pressure and interaction potential, consistent with our observation that both nonlinearities obey the same finite-size scaling. Thus, the zero-pressure limit of the nonlinear response is highly singular.
Different from previous modelings of self-propelled particles, we develop a method to propel the particles with a constant average velocity instead of a constant force. This constant propulsion velocity (CPV) approach is validated by its agreement with the conventional constant propulsion force (CPF) approach in the flowing regime. However, the CPV approach shows its advantage of accessing quasistatic flows of yield stress fluids with a vanishing propulsion velocity, while the CPF approach is usually unable to because of finite system size. Taking this advantage, we realize the cyclic self-propulsion and study the evolution of the propulsion force with propelled particle displacement, both in the quasistatic flow regime. By mapping shear stress and shear rate to propulsion force and propulsion velocity, we find similar rheological behaviors of self-propelled systems to sheared systems, including the yield force gap between the CPF and CPV approaches, propulsion force overshoot, reversible-irreversible transition under cyclic propulsion, and propulsion bands in plastic flows. These similarities suggest the underlying connections between self-propulsion and shear, although they act on systems in different ways.
The interactions between particles in particulate systems are organized in `force networks, mesoscale features that bridge between the particle scale and the scale of the system as a whole. While such networks are known to be crucial in determining the system wide response, extracting their properties, particularly from experimental systems, is difficult due to the need to measure the interparticle forces. In this work, we show by analysis of the data extracted from simulations that such detailed information about interparticle forces may not be necessary, as long as the focus is on extracting the most dominant features of these networks. The main finding is that a reasonable understanding of the time evolution of force networks can be obtained from incomplete information such as total force on the particles. To compare the evolution of the networks based on the completely known particle interactions and the networks based on incomplete information (total force each grain) we use tools of algebraic topology. In particular we will compare simple measures defined on persistence diagrams that provide useful summaries of the force network features.
Using dynamical density functional theory (DDFT) methods we investigate the laning instability of a sheared colloidal suspension. The nonequilibrium ordering at the laning transition is driven by non-affine particle motion arising from interparticle interactions. Starting from a DDFT which incorporates the non-affine motion, we perform a linear stability analysis that enables identification of the regions of parameter space where lanes form. We illustrate our general approach by applying it to a simple one-component fluid of soft penetrable particles.
The similarity in mechanical properties of dense active matter and sheared amorphous solids has been noted in recent years without a rigorous examination of the underlying mechanism. We develop a mean-field model that predicts that their critical behavior should be equivalent in infinite dimensions, up to a rescaling factor that depends on the correlation length of the applied field. We test these predictions in 2d using a new numerical protocol, termed `athermal quasi-static random displacement, and find that these mean-field predictions are surprisingly accurate in low dimensions. We identify a general class of perturbations that smoothly interpolate between the uncorrelated localized forces that occur in the high-persistence limit of dense active matter, and system-spanning correlated displacements that occur under applied shear. These results suggest a universal framework for predicting flow, deformation, and failure in active and sheared disordered materials.
In a recent paper [S. Mandal et al., Phys. Rev. E 88, 022129 (2013)] the nature of spatial correlations of plasticity in hard sphere glasses was addressed both via computer simulations and in experiments. It was found that the experimentally obtained correlations obey a power law whereas the correlations from simulations are better fitted by an exponential decay. We here provide direct evidence--- via simulations of a hard sphere glass in 2D---that this discrepancy is a consequence of the finite system size in the 3D simulations. By extending the study to a 2D soft disk model at zero temperature, the robustness of the power-law decay in sheared amorphous solids is underlined. Deviations from a power law occur when either reducing the packing fraction towards the supercooled regime in the case of hard spheres or changing the dissipation mechanism from contact dissipation to a mean-field type drag for the case of soft disks.