Statistical mechanics of an integrable system

TL;DR

This paper demonstrates through numerical simulations that integrable systems like the Toda chain can reach thermal equilibrium, supporting the idea that chaos is not essential for thermalization in high-dimensional systems.

Contribution

It provides evidence that integrable Hamiltonian systems can thermalize without chaos, challenging traditional views on the necessity of chaos for equilibrium in statistical mechanics.

Findings

Fast thermalization in the Toda chain from atypical initial conditions.

Equilibrium fluctuations match predictions from the Gibbs ensemble.

No conflict with the Generalized Gibbs Ensemble for the Toda model.

Abstract

We provide here an explicit example of Khinchin's idea that the validity of equilibrium statistical mechanics in high dimensional systems does not depend on the details of the dynamics. This point of view is supported by extensive numerical simulation of the one-dimensional Toda chain, an integrable non-linear Hamiltonian system where all Lyapunov exponents are zero by definition. We study the relaxation to equilibrium starting from very atypical initial conditions and focusing on energy equipartion among Fourier modes, as done in the original Fermi-Pasta-Ulam-Tsingou numerical experiment. We find evidence that in the general case, i.e., not in the perturbative regime where Toda and Fourier modes are close to each other, there is a fast reaching of thermal equilibrium in terms of a single temperature. We also find that equilibrium fluctuations, in particular the behaviour of specific heat as function of temperature, are in agreement with analytic predictions drawn from the ordinary Gibbs ensemble, still having no conflict with the established validity of the Generalized Gibbs Ensemble for the Toda model. Our results suggest thus that even an integrable Hamiltonian system reaches thermalization on the constant energy hypersurface, provided that the considered observables do not strongly depend on one or few of the conserved quantities. This suggests that dynamical chaos is irrelevant for thermalization in the large-$N$ limit, where any macroscopic observable reads of as a collective variable with respect to the coordinate which diagonalize the Hamiltonian. The possibility for our results to be relevant for the problem of thermalization in generic quantum systems, i.e., non-integrable ones, is commented at the end.