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 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.
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.
We study numerically time evolution in classical lattices with weak or moderate nonlinearity which leads to interactions between linear modes. Our results show that in a certain strength range a moderate nonlinearity generates a dynamical thermalization process which drives the system to the quantum Gibbs distribution of probabilities, or average oscillation amplitudes. The effective dynamical temperature of the lattice varies from large positive to large negative values depending on the energy of the initially excited modes. This quantum Gibbs distribution is drastically different from the usually expected energy equipartition over linear modes corresponding to a regime of classical thermalization. Possible experimental observations of this dynamical thermalization are discussed for cold atoms in optical lattices, nonlinear photonic lattices and optical fiber arrays.
We study numerically time evolution in classical lattices with weak or moderate nonlinearity which leads to interactions between linear modes. Our results show that in a certain strength range a moderate nonlinearity generates a dynamical thermalization process which drives the system to the quantum Gibbs distribution of probabilities, or average oscillation amplitudes. The effective dynamical temperature of the lattice varies from large positive to large negative values depending on the energy of the initially excited modes. This quantum Gibbs distribution is drastically different from the usually expected energy equipartition over linear modes corresponding to a regime of classical thermalization. Possible experimental observations of this dynamical thermalization are discussed for cold atoms in optical lattices, nonlinear photonic lattices and optical fiber arrays.