Recent works proved a hydrodynamic limit for periodically forced atom chains with harmonic interaction and pinning, together with momentum flip [16, 17]. When energy is the only conserved quantity, one would expect similar results in the anharmonic case, as conjectured for the temperature profile and energy flux in [18]. However, outside the harmonic case, explicit computations are generally no longer possible, thus making a rigorous proof of this hydrodynamic limit difficult. Consequently, we numerically investigate the plausibility of this limit for the particular case of a chain with \(\beta \) –FPUT interactions and harmonic pinning. We present our simulation results suggesting that the PDE for the limiting temperature profile and the Green–Kubo type formula for the limiting energy current conjectured in [18] are correct. We then use this Green–Kubo type formula to investigate the relationship between the energy current and period of the forcing. This relationship is investigated in the case of significant rate of momentum flip, small rate of momentum flip and no momentum flip. We compare the relationship observed in the anharmonic case to that of the harmonic case for which explicit formulae are available [16].
