Heat Transport of the $β$-Fermi--Pasta--Ulam--Tsingou chain in the long-wave limit

Abstract

Resolving the anomalous conductivity exponent of the symmetric $β$-Fermi--Pasta--Ulam--Tsingou chain by molecular dynamics can require very large systems because of long finite-size crossovers and thermal-contact resistance. Motivated by this computational challenge, we derive a long-wavelength continuum description and investigate whether the kinetic-theory scaling $κ\propto L^{2/5}$ becomes accessible with a moderate number of numerical degrees of freedom. The nonlinear elastic field retains the cubic stress of the microscopic interaction and exchanges heat with Langevin reservoirs. A flux-conservative spatial discretization constructs the force and energy current from the same stress, providing a consistent interior transport estimator. For $L>8$, corresponding to approximately $10^3$ mesh nodes and above at the reference resolution, the fitted exponent is $0.399\pm0.004$. Mesh refinement supports the stability of this exponent between the two finer resolutions, although the conductivity amplitude remains resolution-dependent. The main result is thus an accessible finite-size regime near $2/5$, rather than convergence of the absolute conductivity as the mesh spacing vanishes. This formulation provides a practical route to studying anomalous transport in anharmonic systems and a conservative continuum framework for investigating energy flow in nonlinear elastic media.