We perform a thorough investigation of the first FPUT recurrence in the β-FPUT chain for both positive and negative β. We show numerically that the rescaled FPUT recurrence time Tr = tr/(N + 1) 3 depends, for large N , only on the parameter S ≡ Eβ(N + 1). Our numerics also reveal that for small |S|, Tr is linear in S with positive slope for both positive and negative β. For large |S|, Tr is proportional to |S| −1/2 for both positive and negative β but with different multiplicative constants. In the continuum limit, the β-FPUT chain approaches the modified Korteweg-de Vries (mKdV) equation, which we investigate numerically to better understand the FPUT recurrences on the lattice. In the continuum, the recurrence time closely follows the |S| −1/2 scaling and can be interpreted in terms of solitons, as in the case of the KdV equation for the α chain. The difference in the multiplicative factors between positive and negative β arises from soliton-kink interactions which exist only in the negative β case. We complement our numerical results with analytical considerations in the nearly linear regime (small |S|) and in the highly nonlinear regime (large |S|). For the former, we extend previous results using a shifted-frequency perturbation theory and find a closed form for Tr which depends only on S. In the latter regime, we show that Tr ∝ |S| −1/2 is predicted by the soliton theory in the continuum limit. We end by discussing the striking differences in the amount of energy mixing as well as the existence of the FPUT recurrences between positive and negative β and offer some remarks on the thermodynamic limit.

Classical statistical mechanics has long relied on assumptions such as the equipartition theorem to understand the behavior of the complicated systems of many particles. The successes of this approach are well known, but there are also many well-known issues with classical theories. For some of these, the introduction of quantum mechanics is necessary, e.g., the ultraviolet catastrophe. However, more recently, the validity of assumptions such as the equipartition of energy in classical systems was called into question. For instance, a detailed analysis of a simplified model for blackbody radiation was apparently able to deduce the Stefan–Boltzmann law using purely classical statistical mechanics. This novel approach involved a careful analysis of a “metastable” state which greatly delays the approach to equilibrium. In this paper, we perform a broad analysis of such a metastable state in the classical Fermi–Pasta–Ulam–Tsingou (FPUT) models. We treat both the α-FPUT and β-FPUT models, exploring both quantitative and qualitative behavior. After introducing the models, we validate our methodology by reproducing the well-known FPUT recurrences in both models and confirming earlier results on how the strength of the recurrences depends on a single system parameter. We establish that the metastable state in the FPUT models can be defined by using a single degree-of-freedom measure—the spectral entropy (η)—and show that this measure has the power to quantify the distance from equipartition. For the α-FPUT model, a comparison to the integrable Toda lattice allows us to define rather clearly the lifetime of the metastable state for the standard initial conditions. We next devise a method to measure the lifetime of the metastable state tm in the α-FPUT model that reduces the sensitivity to the exact initial conditions. Our procedure involves averaging over random initial phases in the plane of initial conditions, the P1-Q1 plane. Applying this procedure gives us a power-law scaling for tm, with the important result that the power laws for different system sizes collapse down to the same exponent as Eα2→0. We examine the energy spectrum E(k) over time in the α-FPUT model and again compare the results to those of the Toda model. This analysis tentatively supports a method for an irreversible energy dissipation process suggested by Onorato et al.: four-wave and six-wave resonances as described by the “wave turbulence” theory. We next apply a similar approach to the β-FPUT model. Here, we explore in particular the different behavior for the two different signs of β. Finally, we describe a procedure for calculating tm in the β-FPUT model, a very different task than for the α-FPUT model, because the β-FPUT model is not a truncation of an integrable nonlinear model.