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.
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.
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.
Abstract We investigate the equilibration process of the strongly coupled quartic Fermi–Pasta–Ulam–Tsingou model by adding Langevin baths to the ends of the chain. The time evolution of the system is investigated by means of extensive numerical simulations and shown to match the results expected from equilibrium statistical mechanics in the time-asymptotic limit. Upon increasing the nonlinear coupling, the thermalization of the energy spectrum displays an increasing asymmetry in favor of small-scale, high-frequency modes, which relax significantly faster than the large-scale, low-frequency ones. The global equilibration time is found to scale linearly with system size and shown to exhibit a power-law decay with the strength of the nonlinearity and temperature. Nonlinear interaction adds to energy distribution among modes, thus speeding up the thermalization process.
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.
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.
Previously we numerically showed that thermalization can occur in an oscillating Bose-Einstein condensate (BEC) with a disordered harmonic trap when healing length ξ of the condensate is shorter than the correlation length σD of the disorder [see, for example, the experiment reported in Phys. Rev. A \textbf{82}, 033603 (2010)]. In this work, we investigate the weakly interacting or Anderson localization regime ξ>σD and show that the oscillating BEC can also exhibit a relaxation process from nonequilibrium to equilibrium. In such an isolated quantum system, energy and particle number are conserved and the irreversible evolution towards thermodynamic equilibrium is induced by the disorder. The thermodynamic equilibrium is evidenced by the maximized entropy S[nk] in which the waveaction spectrum nk follows the Rayleigh-Jeans distribution. Besides, unlike a monotonic irreversible process of thermalization to equilibrium, the Fermi-Pasta-Ulam-Tsingou recurrence arises in this system, manifested by the oscillation of the non-equilibrium entropy.