Mean-field density of states of a small-world model and a jammed soft spheres model
Abstract
We consider a class of random block matrix models in d dimensions, d 1, motivated by the study of the vibrational density of states (DOS) of soft spheres near the isostatic point. The contact networks of average degree Z = z0 + ζ are represented by random z0-regular graphs (only the circle graph in d=1 with z0=2) to which Erd\"os-Renyi graphs having a small average degree ζ are superimposed. In the case d=1, for ζ small the shifted Kesten-McKay DOS with parameter Z is a mean-field solution for the DOS. Numerical simulations in the z0=2 model, which is the k=1 Newman-Watts small-world model, and in the z0=3 model lead us to conjecture that for ζ 0 the cumulative function of the DOS converges uniformly to that of the shifted Kesten-McKay DOS, in an interval [0, ω0], with ω0 < z0-1 + 1. For 2 d 4, we introduce a cutoff parameter Kd 0.5 modeling sphere repulsion. The case Kd=0 is the random elastic network case, with the DOS close to the Marchenko-Pastur DOS with parameter t=Zd. For Kd large the DOS is close for small ω to the shifted Kesten-McKay DOS with parameter t=Zd; in the isostatic case the DOS has around ω=0 the expected plateau. The boson peak frequency in d=3 with K3 large is close to the one found in molecular dynamics simulations for Z=7 and 8.