Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Universal Time Scale for Thermalization in Two-dimensional Systems

Zhen Wang, Weicheng Fu, Yong Zhang, Hong Wang Zhao

Published May 7, 2020• 1 citation
OpenAlex

Abstract

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.

Related papers

The Fermi-Pasta-Ulam-Tsingou problem after 70 years: toward universal laws of thermalization in lattice systems in the thermodynamic limit

Weicheng Fu, Zhen Wang, Wei Lin +4 more

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 $T_{\mathrm{eq}}\propto 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 $T_{\mathrm{eq}}\propto 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.

arXiv:2603.23347•Mar 24, 2026•0 citations

A general multi-wave quasi-resonance theory for lattice energy diffusion

Lin, Wei, Weicheng Fu, Wang, Zhen +2 more

In this letter, a multi-wave quasi-resonance framework is established to analyze energy diffusion in classical lattices, uncovering that it is fundamentally determined by the characteristics of eigenmodes. Namely, based on the presence and the absence of extended modes, lattices fall into two universality classes with qualitatively different thermalization behavior. In particular, we find that while the one with extended modes can be thermalized under arbitrarily weak perturbations in the thermodynamic limit, the other class can be thermalized only when perturbations exceed a certain threshold, revealing for the first time the possibility that a lattice cannot be thermalized, violating the hypothesis of statistical mechanics. Our study addresses conclusively the renowned Fermi-Pasta-Ulam-Tsingou problem for large systems under weak perturbations, underscoring the pivotal roles of both extended and localized modes in facilitating energy diffusion and thermalization processes.

arXiv:2404.15147•Apr 23, 2024•0 citations

Effects of weak disorder on the thermalization of Fermi–Pasta–Ulam–Tsingou model

Lulu Sun, Zhenjun Zhang, Peiqing Tong

Abstract We study the effects of two kinds of weak disorders on the thermalization of the Fermi–Pasta–Ulam–Tsingou model by extensive numerical simulations. The disorders are introduced to the mass of atom or coefficient of the quadratic term of potential energy. The initial energy is distributed equally among some lowest frequency modes. We find that the energy transports to high-frequency modes with time t and eventually approaches energy equipartition faster with either weak disorder than that in the homogeneous case. That means weak disorders accelerate the process of thermalization. We further study the effects of two kinds of disorders on the scaling law of equipartition time T eq. We find that T eq satisfies the following scaling law: T eq ∼ (ɛ)−a (|α|)−b for different disorder strengths in the thermodynamic limit. It is found that the exponent a ≈ 1.0 while b depends on the strength of disorder, which are different from b = 2a in the homogeneous case.

New Journal of Physics•May 28, 2020•5 citations

Wave-Turbulence Origin of the Instability of Anderson Localization against Many-Body Interactions

Zhen Wang, Weicheng Fu, Yong Zhang +1 more

Whether Anderson localization is robust against many-body interactions and, closely related, whether a disordered many-body system can be thermalized are long outstanding issues. In this Letter, we address these issues with the wave-turbulence theory. We show that, in general, the thermalization time in one-dimensional disordered lattice systems is inversely proportional to the squared interaction strength in the thermodynamic limit. It leads to the conclusion that such systems can always be thermalized by arbitrarily weak many-body interactions and thus the localized states are unstable.

arXiv:1903.09502•Physical Review Letters•May 7, 2020•40 citations

Nonintegrability and thermalization of one-dimensional diatomic lattices

Weicheng Fu, Yong Zhang, Hong Wang Zhao

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.

arXiv:1907.00308•Physical review. E•Nov 4, 2019•18 citations

Dynamic crossover towards energy equipartition in the Fermi-Pasta-Ulam-Tsingou β model with long-range interactions

Jian Wang, Ai-chen Li

Energy equipartition can be established in short-range systems after the dynamic process of thermalization. However, energy distribution between different degrees of freedom in systems with long-range interactions is unclear. We study the dynamics of energy relaxation in the Fermi-Pasta-Ulam-Tsingou β model with long-range quartic interactions, which decay as 1/d^{δ} with d being the lattice distance. The dynamic crossover of a mode-energy distribution from localized to equipartitioned with the increase of the power δ is observed. A transition of mode-energy distribution is identified around the value of δ=1, which usually serves as the distinction between strong and weak long-range couplings. We elucidate that the varying frequency overlapping of the mode-energy power spectrum is responsible for this dynamic crossover. Through further calculation of the spectral entropy, the minimum duration of quasistationary states, τ_{QSS}, is found at δ=2, which may provide possible dynamic explanations for the peculiar behavior of heat transport in long-range lattice chains. In addition, the double scaling in τ_{QSS} as a function of energy density is also observed in our long-range lattices. Our results not only contribute to understanding the dynamics of energy relaxation in long-range systems, but also shed light on the longstanding problem of thermalization and low-dimensional heat transport in short-range systems.

Physical review. E•Jul 25, 2022•6 citations
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−2T_{\mathrm{eq}}\propto g^{-2}Teq​∝g−2, where ggg 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 ggg decreases, and thermalization becomes controlled by higher-order processes. Numerical studies reveal successive regimes Teq∝g−γT_{\mathrm{eq}}\propto g^{-γ}Teq​∝g−γ with γ=2,4,6γ=2,4,6γ=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.
In this letter, a multi-wave quasi-resonance framework is established to analyze energy diffusion in classical lattices, uncovering that it is fundamentally determined by the characteristics of eigenmodes. Namely, based on the presence and the absence of extended modes, lattices fall into two universality classes with qualitatively different thermalization behavior. In particular, we find that while the one with extended modes can be thermalized under arbitrarily weak perturbations in the thermodynamic limit, the other class can be thermalized only when perturbations exceed a certain threshold, revealing for the first time the possibility that a lattice cannot be thermalized, violating the hypothesis of statistical mechanics. Our study addresses conclusively the renowned Fermi-Pasta-Ulam-Tsingou problem for large systems under weak perturbations, underscoring the pivotal roles of both extended and localized modes in facilitating energy diffusion and thermalization processes.
Abstract We study the effects of two kinds of weak disorders on the thermalization of the Fermi–Pasta–Ulam–Tsingou model by extensive numerical simulations. The disorders are introduced to the mass of atom or coefficient of the quadratic term of potential energy. The initial energy is distributed equally among some lowest frequency modes. We find that the energy transports to high-frequency modes with time t and eventually approaches energy equipartition faster with either weak disorder than that in the homogeneous case. That means weak disorders accelerate the process of thermalization. We further study the effects of two kinds of disorders on the scaling law of equipartition time T eq. We find that T eq satisfies the following scaling law: T eq ∼ (ɛ)−a (|α|)−b for different disorder strengths in the thermodynamic limit. It is found that the exponent a ≈ 1.0 while b depends on the strength of disorder, which are different from b = 2a in the homogeneous case.
Whether Anderson localization is robust against many-body interactions and, closely related, whether a disordered many-body system can be thermalized are long outstanding issues. In this Letter, we address these issues with the wave-turbulence theory. We show that, in general, the thermalization time in one-dimensional disordered lattice systems is inversely proportional to the squared interaction strength in the thermodynamic limit. It leads to the conclusion that such systems can always be thermalized by arbitrarily weak many-body interactions and thus the localized states are unstable.
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.
Energy equipartition can be established in short-range systems after the dynamic process of thermalization. However, energy distribution between different degrees of freedom in systems with long-range interactions is unclear. We study the dynamics of energy relaxation in the Fermi-Pasta-Ulam-Tsingou β model with long-range quartic interactions, which decay as 1/d^{δ} with d being the lattice distance. The dynamic crossover of a mode-energy distribution from localized to equipartitioned with the increase of the power δ is observed. A transition of mode-energy distribution is identified around the value of δ=1, which usually serves as the distinction between strong and weak long-range couplings. We elucidate that the varying frequency overlapping of the mode-energy power spectrum is responsible for this dynamic crossover. Through further calculation of the spectral entropy, the minimum duration of quasistationary states, τ_{QSS}, is found at δ=2, which may provide possible dynamic explanations for the peculiar behavior of heat transport in long-range lattice chains. In addition, the double scaling in τ_{QSS} as a function of energy density is also observed in our long-range lattices. Our results not only contribute to understanding the dynamics of energy relaxation in long-range systems, but also shed light on the longstanding problem of thermalization and low-dimensional heat transport in short-range systems.