Mean-field density of states of a small-world model and a jammed soft spheres model

TL;DR

This paper studies the vibrational density of states in small-world and jammed sphere models using random matrix theory, revealing how network structure and sphere repulsion influence vibrational spectra.

Contribution

It introduces a mean-field approach to analyze the density of states in small-world and jammed sphere models, connecting network topology with vibrational properties.

Findings

The DOS converges to the shifted Kesten-McKay distribution as small-world randomness vanishes.

Sphere repulsion cutoff parameter affects the DOS, approaching the Marchenko-Pastur distribution.

The boson peak frequency in 3D models aligns with molecular dynamics simulations.

Abstract

We consider a class of random block matrix models in $d$ dimensions, $d \ge 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 = z_0 + \zeta$ are represented by random $z_0$-regular graphs (only the circle graph in $d=1$ with $z_0=2$) to which Erd\"os-Renyi graphs having a small average degree $\zeta$ are superimposed. In the case $d=1$, for $\zeta$ small the shifted Kesten-McKay DOS with parameter $Z$ is a mean-field solution for the DOS. Numerical simulations in the $z_0=2$ model, which is the $k=1$ Newman-Watts small-world model, and in the $z_0=3$ model lead us to conjecture that for $\zeta \to 0$ the cumulative function of the DOS converges uniformly to that of the shifted Kesten-McKay DOS, in an interval $[0, \omega_0]$, with $\omega_0 < \sqrt{z_0-1} + 1$. For $2 \le d \le 4$, we introduce a cutoff parameter $K_d \le 0.5$ modeling sphere repulsion. The case $K_d=0$ is the random elastic network case, with the DOS close to the Marchenko-Pastur DOS with parameter $t=\frac{Z}{d}$. For $K_d$ large the DOS is close for small $\omega$ to the shifted Kesten-McKay DOS with parameter $t=\frac{Z}{d}$; in the isostatic case the DOS has around $\omega=0$ the expected plateau. The boson peak frequency in $d=3$ with $K_3$ large is close to the one found in molecular dynamics simulations for $Z=7$ and $8$.