---
title: Fermi-Pasta-Ulam-Tsingou Equation Overview
url: https://www.emergentmind.com/topics/fermi-pasta-ulam-tsingou-equation
type: topic
---

# Fermi-Pasta-Ulam-Tsingou Equation Overview

The Fermi–Pasta–Ulam–Tsingou (FPUT) equation denotes, in its primary sense, the discrete Hamiltonian equations of motion for a one-dimensional lattice of particles coupled by weakly anharmonic nearest-neighbor springs; in modern usage it also denotes the continuum and kinetic equations that arise from this lattice in distinguished asymptotic regimes. The model is central to the study of recurrence, metastability, prethermalization, resonant energy transfer, and eventual equipartition. In current treatments, the term encompasses the $\alpha$-chain, the $\beta$-chain, and the combined $\alpha+\beta$ chain, together with their KdV, mKdV, Gardner, Burgers, Toda, and wave-kinetic descriptions [2509.17753].

## 1. Discrete Hamiltonian formulation and normal-mode structure

In the periodic setting, the FPUT chain consists of $N$ identical masses with nearest-neighbor interaction and Hamiltonian
\[
H = H_0 + H_3 + H_4,
\]
\[
H_0 = \sum_{j=1}^N \left[ \frac{p_j^2}{2m} + \frac{\kappa}{2}(q_j-q_{j+1})^2 \right],\qquad
H_3 = \frac{\alpha}{3}\sum_{j=1}^N (q_j-q_{j+1})^3,\qquad
H_4 = \frac{\beta}{4}\sum_{j=1}^N (q_j-q_{j+1})^4.
\]
The corresponding Newton equation is
\[
m \ddot{q}_j = \kappa(q_{j+1}+q_{j-1}-2q_j)
+ \alpha\big[(q_{j+1}-q_j)^2-(q_j-q_{j-1})^2\big]
+ \beta\big[(q_{j+1}-q_j)^3-(q_j-q_{j-1})^3\big].
\]
This periodic formulation is the common starting point for resonance analysis and wave-kinetic theory [1810.06902].

A second standard realization uses fixed-end boundary conditions. For the dimensionless $\beta$-FPUT chain with $q_0=q_{N+1}=0$, the Hamiltonian is
\[
H=\sum_{j=1}^{N}\left[\frac{1}{2}p_j^2+\frac{1}{2}(q_{j+1}-q_j)^2+\frac{\beta}{4}(q_{j+1}-q_j)^4\right],
\]
and the equations of motion are
\[
\ddot{q}_j=(q_{j+1}-2q_j+q_{j-1})+\beta\Big[(q_{j+1}-q_j)^3-(q_j-q_{j-1})^3\Big],\qquad j=1,\dots,N.
\]
For fixed ends, the linear frequencies are
\[
\omega_k=2\sin\!\Big(\frac{k\pi}{2(N+1)}\Big),
\]
whereas for the periodic dimensionless chain they are
\[
\omega(k)=2\sin(\pi k/N).
\]
The standard linear-mode energy is
\[
E_k=\frac12\left(P_k^2+\omega_k^2Q_k^2\right),
\]
and the paradigmatic recurrence experiment excites the first mode, typically through
\[
q_j(0)=A\sin\!\Big(\frac{j\pi}{N+1}\Big),\qquad p_j(0)=0
\]
for fixed ends, or its periodic traveling-wave analogue for the ring [1908.00564].

The distinction between periodic and fixed-end chains is not merely technical. It changes the Fourier basis, the dispersion relation, and the admissible exact resonances. At the same time, the two settings display the same broad phenomenology: low-dimensional initial excitations can remain organized for unexpectedly long times, and the approach to thermalization is highly structured rather than immediate. This suggests that the “FPUT equation” is best understood as a family of closely related lattice equations rather than a single normal form.

## 2. Exact resonances, divisibility effects, and finite-$N$ energy transfer

In the weakly nonlinear regime, exact discrete resonances organize irreversible energy transfer in mode space. For the periodic FPUT chain with dispersion $\omega(k)=2\sin(\pi k/N)$, exact $M$-wave resonances satisfy modular momentum conservation together with exact frequency balance. A fundamental structural result is that $3$-wave resonances, and more generally $1\leftrightarrow(M-1)$ processes, are forbidden by strict subadditivity of the dispersion. Consequently, $4$-wave resonances must be of $2\leftrightarrow2$ type, and in the discrete finite chain they are integrable in the Birkhoff or resonant normal-form sense and do not produce spectrum-wide mixing [1810.06902].

Higher-order resonances are therefore decisive. The arithmetic classification based on cyclotomic polynomials shows that $6$-wave resonances always exist for any $N$, whereas $5$-wave resonances exist if $N$ is divisible by $3$ and $N>6$. The same analysis proves that the mode set can decompose into dynamically independent components whose number depends sensitively on the odd divisors of $N$ not divisible by $3$. This divisibility dependence explains why finite-$N$ thermalization routes vary strongly with particle number: when $3\nmid N$, quintets are absent and sextuplets dominate; when $3\mid N$, quintets enable faster component-wise mixing, while sextuplets bridge components and provide a restoring mechanism for full-scale thermalization [1810.06902].

This finite-$N$ picture connects directly to kinetic theory. For the $\beta$-FPUT chain in the kinetic limit $N\to\infty$, $\beta\to0$ with $\beta=N^{-\gamma}$ and $0<\gamma<1$, the resonant $4$-wave channel becomes the leading kinetic mechanism. A rigorous derivation for a reduced resonant evolution proves the wave kinetic equation up to the sub-kinetic time scale
\[
T=N^{-\epsilon}\min\bigl(N,N^{5\gamma/4}\bigr),
\]
after a deterministic phase renormalization removes dangerous self-interactions [2506.02948]. A later rigorous treatment of the full $\beta$-FPUT system incorporates the non-resonant terms directly into the diagrammatic expansion and reaches times up to $T_{\mathrm{kin}}^{2/3}$, with
\[
T_{\mathrm{kin}}=\frac{1}{4\pi\beta^2},
\]
thereby validating the onset of thermalization in the weakly nonlinear large-box regime [2605.19308].

A common misconception is that “the” thermalization mechanism in FPUT chains is uniform across all system sizes. The finite periodic chain is arithmetic: its exact resonant manifolds depend on the divisor structure of $N$, whereas in the thermodynamic limit Fourier space becomes dense and the kinetic description becomes appropriate. The discrete resonance theory and the kinetic theory therefore describe different but compatible asymptotic layers of the same model.

## 3. Continuum reductions and integrable asymptotics

Long-wavelength asymptotics map the discrete FPUT lattice to several integrable or asymptotically integrable continuum equations. For quasi-unidirectional waves in the $\alpha+\beta$ lattice, the first-order continuum reduction is KdV, and the second-order reduced dynamics is governed by a combination of the first two nontrivial equations in the KdV hierarchy. At third order, a combination of the first three nontrivial hierarchy equations appears only when
\[
14\alpha^3-27\alpha\beta+12\gamma=0,
\]
a condition satisfied by the Toda chain. This identifies a precise algebraic criterion for persistence of asymptotic integrability beyond leading order [2010.03520].

For the $\beta$-FPUT chain, the continuum limit is the modified Korteweg–de Vries equation
\[
u_t+\sigma\,6\,u^2u_x+u_{xxx}=0,\qquad \sigma=\operatorname{sign}(\beta).
\]
Here $\beta>0$ yields focusing mKdV with bright solitons
\[
u(x,t)=\pm\sqrt{6v}\,\operatorname{sech}\!\left(\sqrt{v}(x-vt)\right),
\]
whereas $\beta<0$ yields defocusing mKdV with kinks
\[
u(x,t)=\sqrt{3v}\,\tanh\!\left(\sqrt{v}(x-vt)\right).
\]
The sign of $\beta$ therefore changes the solitary-wave content of the continuum description and, through it, the recurrence dynamics of the discrete chain [1908.00564].

When both quadratic and cubic nearest-neighbor nonlinearities are retained, the continuum reduction is the Gardner equation
\[
4v_t+6v^2v_x+12\mu\,v v_x-v_{xxx}=0.
\]
Via Hirota bilinearization, this equation supports exact multi-soliton solutions and table-top soliton molecules generated by velocity resonance. The same framework recovers the cubic-only chain as mKdV and the quadratic-only chain as KdV through the stated Galilean and Miura interconnections [2308.16535].

A different thermodynamic scaling suppresses dispersion and produces a Burgers regime. In that setting, a near-identity Hamiltonian normal form reduces the FPUT dynamics to a pair of generalized inviscid Burgers equations for dressed left- and right-moving fields. The resulting shock time
\[
t_\mathrm{s}
=
\Big(\frac{N}{2\pi\sqrt{2}\,a\,\alpha}\Big)\,\frac{F(\mu)}{\cos\theta}
\]
predicts when the Fourier spectrum develops the universal law
\[
E_k(t_\mathrm{s})\sim k^{-8/3}.
\]
After approximately two shock times, the spectrum reaches a regime close to
\[
E_k\sim k^{-2},
\]
which persists over an extensive time window before equipartition [2208.08818]. A full Hamiltonian perturbative derivation of this Burgers reduction, including the short-time law
\[
E_k\sim t^{2k-2},
\]
has since been given for the periodic $\alpha+\beta$ chain [2407.16534].

These reductions clarify why distinct asymptotic equations coexist in the FPUT literature. KdV and mKdV describe weakly dispersive, long-wave, near-integrable dynamics; Burgers describes the zero-dispersion steepening regime; Toda provides the integrable neighbor governing adiabatic invariants and metastability. This suggests that continuum limits in FPUT theory are regime-selective rather than mutually exclusive [2509.17753].

## 4. Recurrences, metastability, chaos, and ergodization

The classical recurrence problem remains one of the most quantitative parts of FPUT theory. For the fixed-end $\beta$ chain with first-mode excitation, the first recurrence time rescales as
\[
T_r=\frac{t_r}{(N+1)^3},
\qquad
S=E\beta(N+1).
\]
For large $N$, $T_r$ depends only on $S$. In the nearly linear regime, it is linear in $S$ with positive slope on both sides of $\beta=0$, while in the highly nonlinear regime it satisfies
\[
T_r\propto |S|^{-1/2}.
\]
The prefactors differ for $\beta>0$ and $\beta<0$ because the negative-$\beta$ continuum limit involves soliton–kink interactions that accelerate rephasing. The same study reports $\beta$-dependent critical energies above which recurrences do not form within the measured times and argues that, in the thermodynamic limit with extensive energy, recurrences deteriorate and ultimately disappear [1908.00564].

Metastability can also be measured by observables tied to nearby integrable dynamics. In the FPUT-$\alpha$ chain, a Toda integral constructed from $\operatorname{tr}L^4$ acts as an adiabatic invariant and provides an operational equilibrium time. Its ergodization time is systematically longer than both the inverse maximum Lyapunov exponent and the Lyapunov saturation time. For large chains it becomes essentially system-size independent, whereas below a critical size it grows dramatically. The measured crossover is summarized by
\[
\varepsilon_c(N)\approx 11\,N^{-1.9},
\]
which the paper relates to a possible KAM-like regime [2511.08149].

A complementary viewpoint is mean-field rather than perturbative. For the quartic $\beta$ chain with fixed ends, a mean-field Hamiltonian
\[
H_H(\mathbf{q},\mathbf{p}) = K(\mathbf{p}) + U_2(\mathbf{q}) + \frac{\alpha(\varepsilon)}{N}[U_2(\mathbf{q})]^2
\]
is constructed so that
\[
h_N=\frac{H-H_H}{N}
\]
tends to zero in probability and almost surely under the Gibbs measure as $N\to\infty$. The resulting effective normal modes have renormalized frequencies
\[
\Omega_k=\omega_k\sqrt{1+\gamma(\varepsilon)},
\]
and a phenomenological Langevin closure reproduces mode-energy relaxation and spectral line shapes across energies from the quasi-integrable regime to the strongly chaotic one [2512.00929].

Chaos at low energy can itself be modeled as a random perturbation of the Toda chain. For the fixed-end $\alpha+\beta$ chain, the largest Lyapunov exponent is controlled by the interaction between a few soliton-like Toda modes and a background of radiative modes acting as an intrinsic bath. Retaining only one randomly perturbed Toda soliton-like mode reproduces the observed power law
\[
\lambda\sim \epsilon^{1.55\pm0.07}
\]
for large chains, and analogous exponents for quasi-Toda variants with higher-order matching to the Toda potential [2107.12489].

Taken together, these results resolve a second misconception: positive Lyapunov exponents do not by themselves imply rapid ergodization. Weak chaos, metastable packets, adiabatic invariants, and very slow action diffusion can coexist over parametrically long windows before true thermal equilibrium is reached.

## 5. Nonequilibrium transport, damping, and finite-temperature dynamics

In the $\beta$-FPUT chain, the microscopic lattice equations admit a mesoscopic kinetic closure in terms of a spatially nonhomogeneous wave kinetic equation for the occupation density $n(x,k,t)$:
\[
\partial_t n_k(x,t) + v_k \partial_x n_k(x,t) = \mathcal{I}_k[n](x,t),
\]
with $v_k=\cos(k/2)$ and a $4$-wave collision integral generated by the quartic nonlinearity. With thermostats imposed as incoming Rayleigh–Jeans spectra at $x=0$ and $x=L$, the stationary conductivity
\[
\mathcal{K}=\frac{JL}{\Delta T}
\]
crosses over from ballistic scaling $\mathcal{K}\propto L^1$ at small $L$ to anomalous scaling
\[
\mathcal{K}\propto L^\alpha,\qquad \alpha\approx0.4,
\]
at large $L$. The same kinetic analysis identifies a critical scale
\[
k_c(L)\propto L^{-3/10}
\]
separating low-$k$ ballistic phonons from high-$k$ diffusive phonons. In this picture, anomalous conduction arises from the coexistence of nearly noninteracting low-wavenumber modes and locally equilibrated high-wavenumber modes. The same framework yields second-sound–like ballistic peaks with speeds $\approx \pm1$ and a diffusive heat peak with effective diffusion coefficient $D\approx0.24$ in the reported setup [2203.08627].

Finite temperature changes recurrence qualitatively. For the $\alpha$-FPUT chain with a spatially sinusoidal initial temperature profile, the continuum thermoelastic description coupled to the ballistic heat equation gives
\[
T(x,t)=T_b+\Delta T\,J_0(\omega t)\sin(\lambda x)
\]
and a resonant mechanical response
\[
u(x,t)=z(t)\cos(\lambda x),\qquad
z(t)=-\frac{\beta\Delta T}{\lambda}\,\omega t\,J_1(\omega t).
\]
This “ballistic resonance” converts thermal energy into mechanical energy at early times, but at later times the mechanical energy decays monotonically. The study concludes that the well-known zero-temperature recurrence paradox is eliminated at finite temperatures in this setting [1910.12573].

Weak dissipation also changes recurrence symmetry. In an optical NLSE analog of FPUT recurrence, small linear attenuation induces separatrix crossing at multiple critical losses,
\[
\alpha_{c n} = \frac{p^2 e_+e_-}{g(\omega_0)}\frac{a^2}{n},
\]
switching the dynamics between unshifted and shifted recurrence families. This does not suppress recurrence outright, but it destroys the conservative phase-space separation that protects the undamped dynamics [2206.00912].

These nonequilibrium studies emphasize that “recurrence versus thermalization” is not a binary opposition. Boundary driving, ballistic heat transport, and weak damping all preserve recognizable FPUT structures while changing the relevant observables from return peaks to currents, symmetry classes, and transport exponents.

## 6. Generalizations: dimers, nonlocal solitary waves, and disorder

The FPUT equation extends naturally to polyatomic and heterogeneous lattices. In the dimer FPUT lattice, masses and/or spring potentials alternate with period $2$, and the traveling-wave problem can be written as
\[
\Phi_c(\phi,\omega)=c^2\omega^2 M\phi''-\Delta_-(\omega)\,\mathbf{V}'\big(\Delta_+(\omega)\phi\big)=0.
\]
For nonsymmetric dimers, the existence of small-amplitude periodic traveling waves is proved by exploiting the gradient structure and translation invariance of the traveling-wave operator. The key point is that derivative orthogonality,
\[
\langle \Phi_c(\phi,\omega),\phi'\rangle=0,
\]
eliminates one solvability condition and permits a bifurcation argument with a two-dimensional kernel [2412.17733].

Diatomic lattices also support nonlocal solitary waves. In the equal-mass limit, a monatomic solitary wave perturbs into a micropteron: a localized core together with a small periodic tail whose amplitude is
\[
a_c^\mu=\mathcal{O}_c(|\mu|),
\]
rather than beyond all algebraic orders. The construction depends on a hidden solvability condition for an advance–delay operator and on asymptotically sinusoidal Jost solutions [1909.00770]. Numerical continuation then reveals a connected landscape of micropterons, nanopterons, and solitary waves in the diatomic problem and gives evidence that selected diatomic solitary waves are stable under the full FPUT dynamics [2010.02107].

A different generalization introduces heterogeneity directly into the lattice couplings. In a non-Hamiltonian disordered $\alpha$-FPUT variant, variability is modeled by i.i.d. Gaussian factors $v_j$ with tolerance $\tau$. A two-mode multiple-scale reduction yields envelope equations whose effective coupling
\[
\widetilde{B}(\tau)\approx -0.00893 \tau^2 - 0.000084 \tau + 0.90728
\]
crosses zero at
\[
\tau_c \approx 10.0749\%.
\]
Below this threshold, energy transfer out of the first mode is strongly reduced and the chain exhibits $q$-space localization; above it, the reduced dynamics acquires an unbounded sector and the full non-Hamiltonian system can blow up in finite time. Chaos diagnostics based on the maximum Lyapunov exponent and the Smaller Alignment Index show that when localization is strong, the probability of chaos increases with system size [2212.05644].

These extensions show that the FPUT equation is not restricted to a monatomic, homogeneous, conservative chain. Alternating masses, oscillatory tails, and disorder all preserve the core competition between discreteness, nonlinearity, and dispersion, while changing the spectrum of admissible coherent structures and the mechanisms by which energy remains trapped, leaks, or delocalizes.

Source: https://www.emergentmind.com/topics/fermi-pasta-ulam-tsingou-equation