Finite size analysis of zero-temperature jamming transition under applied shear stress

By finding local minima of an enthalpy-like energy, we can generate jammed packings of frictionless spheres under constant shear stress $\sigma$ and obtain the yield stress $\sigma_y$ by sampling the potential energy landscape. For three-dimensional systems with harmonic repulsion, $\sigma_y$ satisfies the finite size scaling with the limiting scaling relation $\sigma_y\sim\phi - \phi_{_{c,\infty}}$, where $\phi_{_{c,\infty}}$ is the critical volume fraction of the jamming transition at $\sigma=0$ in the thermodynamic limit. The width or uncertainty of the yield stress decreases with decreasing $\phi$ and decays to zero in the thermodynamic limit. The finite size scaling implies a length $\xi\sim (\phi-\phi_{_{c,\infty}})^{-\nu}$ with $\nu=0.81\pm 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 quasi-static shear reveals that quasi-static shear tends to explore low-energy states.