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.
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.
We study the relaxation dynamics of interacting one-dimensional fermions with band curvature after a weak quench in the interaction parameter at zero temperature. Our model lies within the class of interacting Luttinger liquids, where the harmonic Luttinger theory is extended by a weak-integrability-breaking phonon scattering term. In order to solve for the nonequilibrium time evolution, we use quantum kinetic equations exploiting the resonant but subleading character of the phonon interaction term. The interplay between phonon scattering and the quadratic Luttinger theory leads to the emergence of three distinct spatiotemporal regimes for the fermionic real-space correlation function. It features the crossover from a prequench to a prethermal state, finally evolving towards a thermal state on increasing length and time scales. The characteristic algebraically decaying real-space correlations in the prethermalized regime become modulated by an amplitude that is decaying in time according to a stretched exponential as an effect of the interactions. The asymptotic thermalization dynamics is governed by energy transport over large distances from the thermalized to the nonthermalized regions via macroscopic, dynamical slow modes. This is revealed in an algebraic decay of the system's effective temperature. The numerical value of the associated exponent agrees with the dynamical critical exponent of the Kardar-Parisi-Zhang universality class. We also discuss a criterion for the applicability of this theory away from the integrable limit of noninteracting fermions.
Predicting the exact many-body quantum dynamics of polarons in materials with strong carrier-phonon interactions presents a fundamental challenge, often necessitating one to adopt approximations that sacrifice the ability to predict the transition from nonequilibrium relaxation to thermodynamic equilibrium. Here, we exploit a recent breakthrough that generalizes the concept of memory beyond its conventional temporal meaning to also encompass space. Specifically, we leverage our finding that the dynamics of observables in systems with local couplings satisfy Green's functions with kernels that are local in time and space. This enables us to employ the dynamics of small lattices over short times to predict the dynamics of thermodynamically large lattices over arbitrarily long timescales while circumventing the deleterious impacts of finite-size effects. We thus interrogate the exact nonequilibrium formation and migration of polarons in one- (1D) and two-dimensional (2D) systems, revealing that their motion approaches diffusive transport only asymptotically in time and system size. We also compare transport in 1D and 2D systems to investigate the effect of dimension in polaron migration physics, illustrating how energy variations can cause localization-a phenomenon observable via current microscopy experiments.
We consider a 3D homogeneous superfluid at low temperature
T with 2 types of excitations, gapless phonons with a linear dispersion relation at low wavenumber, and gapped quasiparticles with a quadratic dispersion relation around extrema. We calculate the scattering amplitude of a phonon on a quasiparticle to leading order in
T for all subsonic quasiparticle velocities, with a
S-matrix formalism between exact asymptotic states dressed by virtual phonons. We then characterize the erratic motion of the quasiparticle in the superfluid due to its unceasing collisions with thermal phonons through mean force
F(k), longitudinal and transverse
k-dependent momentum diffusion coefficients, and spatial diffusion coefficient. At the minimum location
k0 of the dispersion relation, where the velocity vanishes,
F(k) varies linearly with velocity with an isotropic friction coefficient; if
k0=0, the momentum diffusion is also isotropic and
F(k0)=0; if
k0>0, it is not, and
F(k0) is nonzero but subleading with respect to friction by one order in
T. The velocity time correlation function, whose integral is the spatial diffusion coefficient, decays with the mean velocity damping rate if
k0=0; if
k0>0, it has a second exponential component, with an amplitude and a damping rate lower by a factor
∝T (it is the velocity direction thermalization rate). We also characterize force and momentum diffusion close to the stability domain sonic edge. Our general expressions are expected to be exact to leading order in
T. We illustrate them in the BCS approximation, for a fermionic quasiparticle (an unpaired fermion) in a superfluid of spin 1/2 fermions, realisable with cold atoms in flat bottom traps. We also refute the statement of Lerch, Bartosch and Kopietz (2008), that there would be no fermionic quasiparticle in such a superfluid.
We study the dynamics and stability in a strongly interacting resonantly driven two-band model. Using exact numerical simulations, we find a stable regime at large driving frequencies where the time evolution is governed by a local Floquet Hamiltonian that is approximately conserved out to very long times. For slow driving, on the other hand, the system becomes unstable and heats up to infinite temperature. While thermalization is relatively fast in these two regimes (but to different ``temperatures''), in the crossover between them we find slow nonthermalizing time evolution: temporal fluctuations become strong and temporal correlations long lived. Microscopically, we trace back the origin of this nonthermalizing time evolution to the properties of rare Floquet many-body resonances, whose proliferation at lower driving frequency removes the approximate energy conservation, and thus produces thermalization to infinite temperature.