Intraband discrete breathers in disordered nonlinear systems. I. Delocalization

The existence of Discrete Breathers or DBs (also called Intrinsic Localized Modes or ILMs) and multibreathers, is investigated in a simple one-dimensional chain of random anharmonic oscillators with quartic potentials coupled by springs. When the breather frequency is outside and above the linearized (phonon) spectrum, the existence theorems and numerical methods previously used in periodic nonlinear models for finding time-periodic and spatially localized solutions, hold identically in random nonlinear systems. These solutions are extraband discrete breathers (EDBs). When the frequencies penetrate inside the linearized spectrum, the existence theorems do not hold. Our numerical investigations demonstrate that the strict continuation of (localized) EDBs as intraband discrete breathers (IDBs) is impossible because of cascades of bifurcations generating many discontinuities. A detailed analysis of these bifurcations for small systems with increasing sizes, shows that only a relatively small subset of the spatially extended multibreathers can be strictly continued while their frequency varies inside the phonon spectrum. We propose an ansatz for finding the coding sequences of these solutions and continuing safely these multibreathers in finite systems of any size. This continuation ends at a lower limit frequency where the solution annihilates through a bifurcation with another multibreather. A smaller subset of these multibreather solutions can be continued to amplitude zero and become linear localized modes at this limit. Conversely, any linear localized mode can be continued when increasing its frequency as an extended multibreather. Extrapolation of these results to infinite systems yields the main conclusion of this first part which is that nonlinearity in disordered systems (with localized eigenmodes only) restores their capability of energy transportation by generating infinitely many spatially extended time-periodic solutions. This approach yields mainly spatially extended solutions, except sometimes at their bifurcation points. In the second part of this work, which is presented in our next article, we develop an accurate method for calculating in situ localized IDBs.