The Fermi-Pasta-Ulam problem and its underlying integrable dynamics: an approach through Lyapunov Exponents

TL;DR

This paper investigates the dynamics of Fermi-Pasta-Ulam models by numerically computing Lyapunov exponents, analyzing their asymptotic behavior for large systems and energies, and extending theoretical understanding of their chaotic properties.

Contribution

It compares different FPU models near linear and Toda limits by calculating Lyapunov exponents and extends existing theory to a broader class of models.

Findings

Lyapunov exponents follow power-law asymptotics with respect to energy.

Asymptotic behavior is slow, requiring significant computational effort.

Analytic methods extend successfully to linear hierarchy models, but not near Toda.

Abstract

FPU models, in dimension one, are perturbations either of the linear model or of the Toda model; perturbations of the linear model include the usual $\beta$-model, perturbations of Toda include the usual $\alpha+\beta$ model. In this paper we explore and compare two families, or hierarchies, of FPU models, closer and closer to either the linear or the Toda model, by computing numerically, for each model, the maximal Lyapunov exponent $\chi$. We study the asymptotics of $\chi$ for large $N$ (the number of particles) and small $\epsilon$ (the specific energy $E/N$), and find, for all models, asymptotic power laws $\chi\simeq C\epsilon^a$, $C$ and $a$ depending on the model. The asymptotics turns out to be, in general, rather slow, and producing accurate results requires a great computational effort. We also revisit and extend the analytic computation of $\chi$ introduced by Casetti, Livi and Pettini, originally formulated for the $\beta$-model. With great evidence the theory extends successfully to all models of the linear hierarchy, but not to models close to Toda.