Thermalization of strongly disordered nonlinear chains.

Thermalization of systems described by the discrete nonlinear Schrödinger equation, in the strong disorder limit, is investigated both theoretically and numerically. We show that introducing correlations in the disorder potential, while keeping the "effective" disorder fixed (as measured by the localization properties of wave-packet dynamics), strongly facilitates the thermalization process and leads to a standard grand canonical distribution of the probability norms associated with each site.