Nonintegrability is a necessary condition for the thermalization of a generic Hamiltonian system. In practice, the integrability can be broken in various ways. As illustrating examples, we numerically studied the thermalization behaviors of two types of one-dimensional (1D) diatomic chains in the thermodynamic limit. One chain was the diatomic Toda chain whose nonintegrability was introduced by unequal masses. The other chain was the diatomic Fermi-Pasta-Ulam-Tsingou-β chain whose nonintegrability was introduced by quartic nonlinear interaction. We found that these two different methods of destroying the integrability led to qualitatively different routes to thermalization in the near-integrable region, but the thermalization time, T_{eq}, followed the same scaling law; T_{eq} was inversely proportional to the square of the perturbation strength. This law also agreed with the existing results of 1D monatomic lattices. All these results imply that there is a universal scaling law of thermalization that is independent of the method of breaking integrability.
Pressure plays a vital role in changing the transport properties of matter. To understand this phenomenon at a microscopic level, we here focus on a more fundamental problem, i.e., how pressure affects the thermalization properties of solids. As illustrating examples, we study the thermalization behavior of the monatomic chain and the mass-disordered chain of Fermi-Pasta-Ulam-Tsingou-β under different strains in the thermodynamic limit. It is found that the pressure-induced change in integrability results in qualitatively different thermalization processes for the two kinds of chains. However, for both cases, the thermalization time follows the same law-it is inversely proportional to the square of the nonintegrability strength. This result suggests that pressure can significantly change the integrability of a system, which provides a new perspective for understanding the pressure-dependent thermal transport behavior.
The Fermi--Pasta--Ulam--Tsingou (FPUT) problem provides a paradigmatic framework for understanding thermalization in weakly nonlinear many-body Hamiltonian systems. This focused review summarizes major developments in near-integrable dynamics, wave resonances, and phonon kinetic theory, emphasizing recent results on thermalization-time scaling in nonlinear lattices. A coherent picture emerges in the thermodynamic limit based on the eigenmode properties of an appropriate integrable reference system. If the reference system has extended eigenmodes and the leading resonant or quasi-resonant processes form a sufficiently connected network, the thermalization time follows
Teq∝g−2, where
g measures the effective deviation from integrability. This scaling is found broadly in ordered and weakly disordered lattices and is robust to dimensionality, interaction potential, integrability-breaking mechanism, and multimode initial conditions. Identifying the correct integrable reference is essential; in some one-dimensional FPUT-type lattices with cubic interactions, the nearby Toda lattice, rather than the harmonic chain, provides the proper reference. If all reference eigenmodes are localized, spatial-overlap constraints progressively fragment low-order resonance networks as
g decreases, and thermalization becomes controlled by higher-order processes. Numerical studies reveal successive regimes
Teq∝g−γ with
γ=2,4,6, together with weak system-size dependence. Whether this hierarchy persists asymptotically and whether a finite thermalization threshold exists remain open questions. We also discuss finite-size effects, strongly nonintegrable dynamics, heat transport, localization, and higher-dimensional lattices.
Most of the previous studies on classical thermalization focus on the wave-vector space, encountering limitations when extended beyond quasi-integrable regions. In this study, we propose a scheme to study the thermalization of the classical Hamiltonian chain of interacting oscillators in real space by developing a thermalization indicator proposed by Parisi [Europhys. Lett. 40, 357 (1997)0295-507510.1209/epl/i1997-00471-9], which approaches zero in the thermal state. Upon reaching the steady state characterized by the generalized Gibbs ensemble for a harmonic chain, a quench protocol is implemented to change the Hamiltonian to a nonintegrable form instantaneously, thereby preparing nonequilibrium initial states. This approach enables investigations of thermalization in real space, particularly valuable for exploring regions beyond quasi-integrability. For the FPUT-β lattice, we observe that the thermalization time as a function of the nonintegrable strength follows a -2 scaling law in the quasi-integrable region and -1/4 in the strongly integrable region. Moreover, numerical results reveal the thermalization time is proportional to the Lyapunov time, which bridges microscopic chaotic dynamics and the macroscopic thermalization process.
The Fermi-Pasta-Ulam-Tsingou problem, i.e., the problem of energy equipartition among normal modes in a weakly nonlinear lattice, is here studied in two types of two-dimensional (2D) lattices, more precisely in lattices with square cell and triangular cell. We apply the wave-turbulence approach to describe the dynamics and find multi-wave resonances play a major role in the transfer of energy among the normal modes. We show that, in general, the thermalization time in 2D systems is inversely proportional to the squared perturbation strength in the thermodynamic limit. Numerical simulations confirm that the results are consistent with the theoretical prediction no matter systems are translation-invariant or not. It leads to the conclusion that such systems can always be thermalized by arbitrarily weak many-body interactions. Moreover, the validity for disordered lattices implies that the localized states are unstable.
The Salerno model constitutes an intriguing interpolation between the integrable Ablowitz-Ladik (AL) model and the more standard (nonintegrable) discrete nonlinear Schrödinger (DNLS) one. The competition of local on-site nonlinearity and nonlinear dispersion governs the thermalization of this model. Here, we investigate the statistical mechanics of the Salerno one-dimensional lattice model in the nonintegrable case and illustrate the thermalization in the Gibbs regime. As the parameter interpolating between the two limits (from DNLS toward AL) is varied, the region in the space of initial energy and norm densities leading to thermalization expands. The thermalization in the non-Gibbs regime heavily depends on the finite system size; we explore this feature via direct numerical computations for different parametric regimes.