Athermal jamming versus thermalized glassiness in sheared frictionless particles

Considerations on the relaxation time in shear-driven jamming

We study the jamming transition in a model of elastic particles under shear at zero temperature, with a focus on the relaxation time τ_{1}. This relaxation time is from two-step simulations where the first step is the ordinary shearing simulation and the second step is the relaxation of the energy after stopping the shearing. τ_{1} is determined from the final exponential decay of the energy. Such relaxations are done with many different starting configurations generated by a long shearing simulation in which the shear variable γ slowly increases. We study the correlations of both τ_{1}, determined from the decay, and the pressure, p_{1}, from the starting configurations as a function of the difference in γ. We find that the correlations of p_{1} are longer lived than the ones of τ_{1} and find that the reason for this is that the individual τ_{1} is controlled both by p_{1} of the starting configuration and a random contribution which depends on the relaxation path length-the average distance moved by the particles during the relaxation. We further conclude that it is γ_{τ}, determined from the correlations of τ_{1}, which is the relevant one when the aim is to generate data that may be used for determining the critical exponent that characterizes the jamming transition.

Comparison of compression versus shearing near jamming, for a simple model of athermal frictionless disks in suspension

Using a simplified model for a non-Brownian suspension, we numerically study the response of athermal, overdamped, frictionless disks in two dimensions to isotropic and uniaxial compression, as well as to pure and simple shearing, all at finite constant strain rates ε[over ̇]. We show that isotropic and uniaxial compression result in the same jamming packing fraction ϕ_{J}, while pure-shear- and simple-shear-induced jamming occurs at a slightly higher ϕ_{J}^{*}, consistent with that found previously for simple shearing. A critical scaling analysis of pure shearing gives critical exponents consistent with those previously found for both isotropic compression and simple shearing. Using orientational order parameters for contact bond directions, we compare the anisotropy of the force and contact networks at both lowest nematic order, as well as higher 2n-fold order.

Universal Relaxation Dynamics of Sphere Packings below Jamming

We show that non-Brownian suspensions of repulsive spheres below jamming display a slow relaxational dynamics with a characteristic timescale that diverges at jamming. This slow timescale is fully encoded in the structure of the unjammed packing and can be readily measured via the vibrational density of states. We show that the corresponding dynamic critical exponent is the same for randomly generated and sheared packings. Our results show that a wide variety of physical situations, from suspension rheology to algorithmic studies of the jamming transition are controlled by a unique diverging timescale, with a universal critical exponent.

Dimensionality and Viscosity Exponent in Shear-driven Jamming

Collections of bidisperse frictionless particles at zero temperature in three dimensions are simulated with a shear-driven dynamics with the aim to compare with the behavior in two dimensions. Contrary to the prevailing picture, and in contrast to results from isotropic jamming from compression or quench, we find that the critical exponents in three dimensions are different from those in two dimensions and conclude that shear-driven jamming in two and three dimensions belong to different universality classes.

The physics of jamming for granular materials: a review

Granular materials consist of macroscopic grains, interacting via contact forces, and unaffected by thermal fluctuations. They are one of a class systems that undergo jamming, i.e. a transition between fluid-like and disordered solid-like states. Roughly twenty years ago, proposals by Cates et al for the shear response of colloidal systems and by Liu and Nagel, for a universal jamming diagram in a parameter space of packing fraction, ϕ, shear stress, τ, and temperature, T raised key questions. Contemporaneously, experiments by Howell et al and numerical simulations by Radjai et al and by Luding et al helped provide a starting point to explore key insights into jamming for dry, cohesionless, granular materials. A recent experimental observation by Bi et al is that frictional granular materials have a a re-entrant region in their jamming diagram. In a range of ϕ, applying shear strain, γ, from an initially force/stress free state leads to fragile (in the sense of Cates et al), then anisotropic shear jammed states. Shear jamming at fixed ϕ is presumably conjugate to Reynolds dilatancy, involving dilation under shear against deformable boundaries. Numerical studies by Radjai and Roux showed that Reynolds dilatancy does not occur for frictionless systems. Recent numerical studies by several groups show that shear jamming occurs for finite, but not infinite, systems of frictionless grains. Shear jamming does not lead to known ordering in position space, but Sarkar et al showed that ordering occurs in a space of force tiles. Experimental studies seeking to understand random loose and random close packings (rlp and rcp) and dating back to Bernal have probed granular packings and their response to shear and intruder motion. These studies suggest that rlp's are anisotropic and shear-jammed-like, whereas rcp's are likely isotropically jammed states. Jammed states are inherently static, but the jamming diagram may provide a context for understanding rheology, i.e. dynamic shear in a variety of systems that include granular materials and suspensions.

Shear banding, discontinuous shear thickening, and rheological phase transitions in athermally sheared frictionless disks

We report on numerical simulations of simple models of athermal, bidisperse, soft-core, massive disks in two dimensions, as a function of packing fraction ϕ, inelasticity of collisions as measured by a parameter Q, and applied uniform shear strain rate γ[over ̇]. Our particles have contact interactions consisting of normally directed elastic repulsion and viscous dissipation, as well as tangentially directed viscous dissipation, but no interparticle Coulombic friction. Mapping the phase diagram in the (ϕ,Q) plane for small γ[over ̇], we find a sharp first-order rheological phase transition from a region with Bagnoldian rheology to a region with Newtonian rheology, and show that the system is always Newtonian at jamming. We consider the rotational motion of particles and demonstrate the crucial importance that the coupling between rotational and translational degrees of freedom has on the phase structure at small Q (strongly inelastic collisions). At small Q, we show that, upon increasing γ[over ̇], the sharp Bagnoldian-to-Newtonian transition becomes a coexistence region of finite width in the (ϕ,γ[over ̇]) plane, with coexisting Bagnoldian and Newtonian shear bands. Crossing this coexistence region by increasing γ[over ̇] at fixed ϕ, we find that discontinuous shear thickening can result if γ[over ̇] is varied too rapidly for the system to relax to the shear-banded steady state corresponding to the instantaneous value of γ[over ̇].

Motility-driven glass and jamming transitions in biological tissues

Cell motion inside dense tissues governs many biological processes, including embryonic development and cancer metastasis, and recent experiments suggest that these tissues exhibit collective glassy behavior. To make quantitative predictions about glass transitions in tissues, we study a self-propelled Voronoi (SPV) model that simultaneously captures polarized cell motility and multi-body cell-cell interactions in a confluent tissue, where there are no gaps between cells. We demonstrate that the model exhibits a jamming transition from a solid-like state to a fluid-like state that is controlled by three parameters: the single-cell motile speed, the persistence time of single-cell tracks, and a target shape index that characterizes the competition between cell-cell adhesion and cortical tension. In contrast to traditional particulate glasses, we are able to identify an experimentally accessible structural order parameter that specifies the entire jamming surface as a function of model parameters. We demonstrate that a continuum Soft Glassy Rheology model precisely captures this transition in the limit of small persistence times, and explain how it fails in the limit of large persistence times. These results provide a framework for understanding the collective solid-to-liquid transitions that have been observed in embryonic development and cancer progression, which may be associated with Epithelial-to-Mesenchymal transition in these tissues.

Constant Stress and Pressure Rheology of Colloidal Suspensions

We study the constant stress and pressure rheology of dense hard-sphere colloidal suspensions using Brownian dynamics simulation. Expressing the flow behavior in terms of the friction coefficient-the ratio of shear to normal stress-reveals a shear arrest point from the collapse of the rheological data in the non-Brownian limit. The flow curves agree quantitatively (when scaled) with the experiments of Boyer et al. [Phys. Rev. Lett. 107, 188301 (2011)]. Near suspension arrest, both the shear and the incremental normal viscosities display a universal power law divergence. This work shows the important role of jamming on the arrest of colloidal suspensions and illustrates the care needed when conducting and analyzing experiments and simulations near the flow-arrest transition.

Divergence of Viscosity in Jammed Granular Materials: A Theoretical Approach

A theory for jammed granular materials is developed with the aid of a nonequilibrium steady-state distribution function. The approximate nonequilibrium steady-state distribution function is explicitly given in the weak dissipation regime by means of the relaxation time. The theory quantitatively agrees with the results of the molecular dynamics simulation on the critical behavior of the viscosity below the jamming point without introducing any fitting parameter.

Universal rescaling of flow curves for yield-stress fluids close to jamming

The experimental flow curves of four different yield-stress fluids with different interparticle interactions are studied near the jamming concentration. By appropriate scaling with the distance to jamming all rheology data can be collapsed onto master curves below and above jamming that meet in the shear-thinning regime and satisfy the Herschel-Bulkley and Cross equations, respectively. In spite of differing interactions in the different systems, master curves characterized by universal scaling exponents are found for the four systems. A two-state microscopic theory of heterogeneous dynamics is presented to rationalize the observed transition from Herschel-Bulkley to Cross behavior and to connect the rheological exponents to microscopic exponents for the divergence of the length and time scales of the heterogeneous dynamics. The experimental data and the microscopic theory are compared with much of the available literature data for yield-stress systems.

Aspects of jamming in two-dimensional athermal frictionless systems

In this work we provide an overview of jamming transitions in two dimensional systems focusing on the limit of frictionless particle interactions in the absence of thermal fluctuations. We first discuss jamming in systems with short range repulsive interactions, where the onset of jamming occurs at a critical packing density and where certain quantities show a divergence indicative of critical behavior. We describe how aspects of the dynamics change as the jamming density is approached and how these dynamics can be explored using externally driven probes. Different particle shapes can produce jamming densities much lower than those observed for disk-shaped particles, and we show how jamming exhibits fragility for some shapes while for other shapes this is absent. Next we describe the effects of long range interactions and jamming behavior in systems such as charged colloids, vortices in type-II superconductors, and dislocations. We consider the effect of adding obstacles to frictionless jamming systems and discuss connections between this type of jamming and systems that exhibit depinning transitions. Finally, we discuss open questions such as whether the jamming transition in all these different systems can be described by the same or a small subset of universal behaviors, as well as future directions for studies of jamming transitions in two dimensional systems, such as jamming in self-driven or active matter systems.

Finite size analysis of zero-temperature jamming transition under applied shear stress by minimizing a thermodynamic-like potential

By finding local minima of a thermodynamic-like potential, we generate jammed packings of frictionless spheres under constant shear stress σ and obtain the yield stress σy by sampling the potential energy landscape. For three-dimensional systems with harmonic repulsion, σy satisfies the finite size scaling with the limiting scaling relation σy∼ϕ-ϕc,∞, where ϕc,∞ is the critical volume fraction of the jamming transition at σ=0 in the thermodynamic limit. The finite size scaling implies a length ξ∼(ϕ-ϕc,∞)-ν with ν=0.81±0.05, which turns out to be a robust and universal length scale exhibited as well in the finite size scaling of multiple quantities measured without shear and independent of particle interaction. Moreover, comparison between our new approach and quasistatic shear reveals that quasistatic shear tends to explore low-energy states.

Connecting the random organization transition and jamming within a unifying model system

While the random organization transition describes the change from reversible to irreversible dynamics in a nonequilibrium system, the athermal jamming transition at zero shear rate occurs when particles can no longer avoid overlaps. Despite the obvious differences between these two transitions, we show that they both occur within the same model packing problem. In this unifying model system the particles are first randomly distributed and then displaced in each step if they overlap. For random displacements we obtain a random organization transition, while jamming occurs in the case of deterministic shifts. We also analyze the critical behavior of random organization. Our results show that random organization and jamming are opposite limits of random sphere packings, and we expect that various equilibrium and nonequilibrium transitions can be formulated as related intermediate packing problems.

Origin of rigidity in dry granular solids

Solids are distinguished from fluids by their ability to resist shear. In traditional solids, the resistance to shear is associated with the emergence of broken translational symmetry as exhibited by a nonuniform density pattern. In this work, we focus on the emergence of shear rigidity in a class of solids where this paradigm is challenged. Dry granular materials have no energetically or entropically preferred density modulations. We show that, in contrast to traditional solids, the emergence of shear rigidity in these granular solids is a collective process, which is controlled solely by boundary forces, the constraints of force and torque balance, and the positivity of the contact forces. We develop a theoretical framework based on these constraints, which connects rigidity to broken translational symmetry in the space of forces, not positions of grains. We apply our theory to experimentally generated shear-jammed states and show that these states are indeed characterized by a persistent, non-uniform density modulation in force space, which emerges at the shear-jamming transition.