Instability dynamics of nonlinear normal modes in the Fermi–Pasta–Ulam–Tsingou chains

Nonlinear normal modes are periodic orbits that survive in nonlinear chains, whose instability plays a crucial role in the dynamics of many-body Hamiltonian systems toward thermalization. Here we focus on how the stability of nonlinear modes depends on the perturbation strength and the system size to observe whether they have the same behavior in different models. To this end, as illustrating examples, the instability dynamics of the N/2 mode in both the Fermi–Pasta–Ulam–Tsingou-α and -β chains under fixed boundary conditions are studied systematically. Applying the Floquet theory, we show that for both models the stability time T as a function of the perturbation strength λ follows the same behavior; i.e., T ∝ ( λ − λ c ) − 1 2 , where λ c is the instability threshold. The dependence of λ c on N is also obtained. The results of T and λ c agree well with those obtained by the direct molecular dynamics simulations. Finally, the effect of instability dynamics on the thermalization properties of a system is briefly discussed.