Search

We introduce a modified version of the disordered Klein-Gordon lattice model, having two parameters for controlling the disorder strength: $D$, which determines the range of the coefficients of the on-site potentials, and $W$, which defines the strength of the nearest-neighbor interactions. We fix $W=4$ and investigate how the properties of the system's normal modes change as we approach its ordered version, i.e. $D\rightarrow 0$. We show that the probability density distribution of the normal modes' frequencies takes a `U'-shaped profile as $D$ decreases. Furthermore, we use two quantities for estimating the modes' spatial extent, the so-called localization volume $V$ (which is related to the mode's second moment) and the mode's participation number $P$. We show that both quantities scale as $\propto D^{-2}$ when $D$ approaches zero and we numerically verify a proportionality relation between them as $V/P \approx 2.6$., Comment: 6 pages, 7 figures

We numerically investigate the characteristics of chaos evolution during wave packet spreading in two typical one-dimensional nonlinear disordered lattices: the Klein-Gordon system and the discrete nonlinear Schr\"{o}dinger equation model. Completing previous investigations \cite{SGF13} we verify that chaotic dynamics is slowing down both for the so-called `weak' and `strong chaos' dynamical regimes encountered in these systems, without showing any signs of a crossover to regular dynamics. The value of the finite-time maximum Lyapunov exponent $\Lambda$ decays in time $t$ as $\Lambda \propto t^{\alpha_{\Lambda}}$, with $\alpha_{\Lambda}$ being different from the $\alpha_{\Lambda}=-1$ value observed in cases of regular motion. In particular, $\alpha_{\Lambda}\approx -0.25$ (weak chaos) and $\alpha_{\Lambda}\approx -0.3$ (strong chaos) for both models, indicating the dynamical differences of the two regimes and the generality of the underlying chaotic mechanisms. The spatiotemporal evolution of the deviation vector associated with $\Lambda$ reveals the meandering of chaotic seeds inside the wave packet, which is needed for obtaining the chaotization of the lattice's excited part., Comment: 11 pages, 10 figures

We implement several symplectic integrators, which are based on two part splitting, for studying the chaotic behavior of one- and two-dimensional disordered Klein-Gordon lattices with many degrees of freedom and investigate their numerical performance. For this purpose, we perform extensive numerical simulations by considering many different initial energy excitations and following the evolution of the created wave packets in the various dynamical regimes exhibited by these models. We compare the efficiency of the considered integrators by checking their ability to correctly reproduce several features of the wave packets propagation, like the characteristics of the created energy distribution and the time evolution of the systems' maximum Lyapunov exponent estimator. Among the tested integrators the fourth order $ABA864$ scheme \cite{BCFLMM13} showed the best performance as it needed the least CPU time for capturing the correct dynamical behavior of all considered cases when a moderate accuracy in conserving the systems' total energy value was required. Among the higher order schemes used to achieve a better accuracy in the energy conservation, the sixth order scheme $s11ABA82\_6$ exhibited the best performance., Comment: 19 pages, 4 figures