Wyart M, Silbert LE, Nagel SR, Witten TA. Effects of compression on the vibrational modes of marginally jammed solids.
PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2005;
72:051306. [PMID:
16383602 DOI:
10.1103/physreve.72.051306]
[Citation(s) in RCA: 194] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/26/2004] [Revised: 08/18/2005] [Indexed: 05/05/2023]
Abstract
Glasses have a large excess of low-frequency vibrational modes in comparison with most crystalline solids. We show that such a feature is a necessary consequence of the weak connectivity of the solid, and that the frequency of modes in excess is very sensitive to the pressure. We analyze, in particular, two systems whose density D(omega) of vibrational modes of angular frequency omega display scaling behaviors with the packing fraction: (i) simulations of jammed packings of particles interacting through finite-range, purely repulsive potentials, comprised of weakly compressed spheres at zero temperature and (ii) a system with the same network of contacts, but where the force between any particles in contact (and therefore the total pressure) is set to zero. We account in the two cases for the observed (a) convergence of D(omega) toward a nonzero constant as omega-->0, (b) appearance of a low-frequency cutoff omega*, and (c) power-law increase of omega* with compression. Differences between these two systems occur at a lower frequency. The density of states of the modified system displays an abrupt plateau that appears at omega*, below which we expect the system to behave as a normal, continuous, elastic body. In the unmodified system, the pressure lowers the frequency of the modes in excess. The requirement of stability despite the destabilizing effect of pressure yields a lower bound on the number of extra contact per particle deltaz:deltaz> or =p1/2, which generalizes the Maxwell criterion for rigidity when pressure is present. This scaling behavior is observed in the simulations. We finally discuss how the cooling procedure can affect the microscopic structure and the density of normal modes.
Collapse