We study the chaotic behavior of multidimensional Hamiltonian systems in the presence of nonlinearity and disorder. It is known that any localized initial excitation in a large enough linear disordered system spreads for a finite amount of time and then halts forever. This phenomenon is called Anderson localization (AL). What happens to AL when nonlinearity is introduced is an interesting question which has been considered in several studies over the past decades. However, the characteristics and the asymptotic fate of such evolutions still remain an issue of intense debate due to their computational difficulty, especially in systems of more than one spatial dimension. As the spreading of initially localized wave packets is a non-equilibrium thermalization process related to the ergodic and chaotic properties of the system, in our work we investigate the properties of chaos studying the behavior of observables related to the systems tangent dynamics. In particular, we consider the disordered discrete nonlinear Schrodinger (DDNLS) equation of one (1D) and two (2D) spatial dimensions. We present detailed computations of the time evolution of the systems maximum Lyapunov exponent (MLE--$Lambda$), and the related deviation vector distribution (DVD). We find that although the systems MLE decreases in time following a power law $t^{alpha_Lambda}$ with $alpha_Lambda <0$ for both the weak and strong chaos regimes, no crossover to the behavior $Lambda propto t^{-1}$ (which is indicative of regular motion) is observed. In addition, the analysis of the DVDs reveals the existence of random fluctuations of chaotic hotspots with increasing amplitudes inside the excited part of the wave packet, which assist in homogenizing chaos and contribute to the thermalization of more lattice sites.
Download