Critical Scaling at the Jamming Transition

1. Critical Scaling at the Jamming Transition Peter Olsson, Umeå University Stephen Teitel, University of Rochester Supported by: US Department of Energy Swedish High Performance Computing Center North

2. outline • introduction - jamming phase diagram • our model for a granular material • simulations in 2D at T= 0 • scaling collapse for shear viscosity • dynamic correlation length • conclusions

3. flowing ➝ rigid but disordered granular materials large grains ⇒ T= 0 upon increasing the volume density of particles above a critical value the sudden appearance of a finite shear stiffness signals a transition from a flowing state to a rigid but disordered state - this is the jamming transition “point J” sheared foamspolydisperse densely packed gas bubbles upon decreasing the applied shear stress below a critical yield stress, the foam ceases to flow and behaves like an elastic solid structural glass upon decreasing the temperature, the viscosity of a liquid grows rapidly and the liquid freezes into a disordered rigid solid animations from Leiden granular group website

4. T glass  jamming J 1/ surface below which states are jammed conjecture by Liu and Nagel (Nature 1998) jamming “point J” is a special critical point in a larger phase diagram with the three axes:  volume density T temperature  applied shear stress (nonequilibrium axis) understanding T = 0 jamming at “point J” in granular materials may have implications for understanding the structural glass transition at finite T here we consider the  plane at T = 0 in 2D

5. ⇒ shear flow in fluid state shear stress  velocity gradient shear viscosity below jamming above jamming shear viscosity of a flowing granular material if jamming is like a critical point we expect

6. for N disks in area LxLy the volume density is non-overlapping ⇒ non-interacting overlapping ⇒ harmonic repulsion r (Durian, PRL 1995 (foams); O’Hern, Silbert, Liu, Nagel, PRE 2003) model granular material bidisperse mixture of soft disks in two dimensionsat T = 0 equal numbers of disks with diameters d1 = 1, d2 = 1.4 interaction V(r)(frictionless)

7. Lees-Edwards boundary conditions create a uniform shear strain  Ly Ly Lx diffusively moving particles (particles in a viscous liquid) interactions strain rate strain  driven by uniform applied shear stress  dynamics

8. simulation parameters Lx = Ly N = 1024 for < 0.844 N = 2048 for ≥ 0.844 t ~ 1/N, integrate with Heun’s method (ttotal) ~ 10, ranging from 1 to 200 depending onNand finite size effects negligible (can’t get too close to c) animation at: = 0.830  0.838 c ≃ 0.8415 = 10-5

9. results for small = 10-5(represents  → 0 limit, “point J”) as N increases, -1() vanishes continuously at c≃ 0.8415 smaller systems jam below c

10. c c results for finite shear stress 

11.  J  c Scaling hypothesis (in analogy to 2nd order equilibrium phase transitions): if we rescale lengths by a factor b,  will scale with  and  according to choose length rescaling factor crossover scaling function crossover scaling variable crossover scaling exponent  scaling about “point J” for finite shear stress  scaling variables  and 

12. scaling collapse of viscosity  stress  is a relevant variable point J is a true 2nd order critical point unclear if jamming remains at finite 

13. transverse velocity correlation function (average shear flow along x)   distance to minimum gives correlation length  correlation length regions separated by  are anti-correlated motion is by rotation of regions of size 

14. scaling collapse of correlation length   diverges at point J

15. ' flowing shear stress  ' jammed  cz c 0 c volume density  “point J” phase diagram in  plane

16. conclusions • point J is a true 2nd order critical point • correlation length diverges at point J • critical scaling extends to non-equilibrium driven steady states at finite shear stress   in agreement with proposal by Liu and Nagel • shear stress  is a relevant variable that changes the critical behavior at point J • jamming transition at finite  remains to be clarified • finite temperature?