Jamming and replica symmetry breaking of weakly disordered crystals

Abstract: We discuss the physics of crystals with small polydispersity near the jamming transition point. For this purpose, we introduce an effective single-particle model taking into account the nearest neighbor structure of crystals. The model can be solved analytically by using the replica method in the limit of large dimensions. In the absence of polydispersity, the replica symmetric solution is stable until the jamming transition point, which leads to the standard scaling of perfect crystals. On the contrary, for finite polydispersity, the model undergoes the full replica symmetry breaking (RSB) transition before the jamming transition point. In the RSB phase, the model exhibits the same scaling as amorphous solids near the jamming transition point. These results are fully consistent with the recent numerical simulations of crystals with polydispersity. The simplicity of the model also allows us to derive the scaling behavior of the vibrational density of states that can be tested in future experiments and numerical simulations.