Disordered harmonic chains with random masses and springs: A combinatorial approach
Phys. Rev. E 113, 014143 – Published 29 January, 2026
DOI: https://doi.org/10.1103/bc9p-fhyz
Abstract
We study harmonic chains with i.i.d. random spring constants and i.i.d. random masses . We introduce a combinatorial approach which allows us to derive a compact and general approximate formula for the complex Lyapunov exponent, in terms of the solutions of two transcendental equations involving the distributions of the spring constants and the masses. Our result makes easy the asymptotic analysis of the low-frequency properties of the eigenmodes (spectral density and localization) for arbitrary disorder distribution, as well as their high-frequency properties. We apply the method to the case of power-law distributions with and with (with ). At low frequency, the spectral density presents the power law , where the exponent exhibits first-order phase transitions on the line and on the line . The exponent of the nondisordered chain () is recovered when and are both finite, i.e., and . The Lyapunov exponent (inverse localization length) shows also a power-law behavior , where the exponent exhibits several phase transitions: the exponent is for or ( or infinite) and when and ( and both finite). In the intermediate region it is given by . On the transition lines, and receive logarithmic corrections. Finally, we illustrate the versatility of our combinatorial approach by considering the case of the Anderson model with random couplings (a model known to be mapped onto the random spring chain for “Dyson type I” disorder).