Elastic collapse in disordered isostatic networks

TL;DR

This paper investigates how elastic moduli in disordered isostatic networks diminish with system size, revealing exponential decay in directed networks and finite moduli when overconstraints are present, with implications for sphere packings and network stability.

Contribution

It provides a comprehensive analysis of elastic collapse in disordered isostatic networks, introducing models for stress growth and the effects of overconstraints on elastic moduli.

Findings

Elastic moduli decrease as inverse power-laws or exponentially with system size.

Directed isostatic networks exhibit exponential decay of elastic moduli.

Presence of overconstraints leads to finite asymptotic elastic moduli.

Abstract

Isostatic networks are minimally rigid and therefore have, generically, nonzero elastic moduli. Regular isostatic networks have finite moduli in the limit of large sizes. However, numerical simulations show that all elastic moduli of geometrically disordered isostatic networks go to zero with system size. This holds true for positional as well as for topological disorder. In most cases, elastic moduli decrease as inverse power-laws of system size. On directed isostatic networks, however, of which the square and cubic lattices are particular cases, the decrease of the moduli is exponential with size. For these, the observed elastic weakening can be quantitatively described in terms of the multiplicative growth of stresses with system size, giving rise to bulk and shear moduli of order exp{-bL}. The case of sphere packings, which only accept compressive contact forces, is considered separately. It is argued that these have a finite bulk modulus because of specific correlations in contact disorder, introduced by the constraint of compressivity. We discuss why their shear modulus, nevertheless, is again zero for large sizes. A quantitative model is proposed that describes the numerically measured shear modulus, both as a function of the loading angle and system size. In all cases, if a density p>0 of overconstraints is present, as when a packing is deformed by compression, or when a glass is outside its isostatic composition window, all asymptotic moduli become finite. For square networks with periodic boundary conditions, these are of order sqrt{p}. For directed networks, elastic moduli are of order exp{-c/p}, indicating the existence of an "isostatic length scale" of order 1/p.