The generalized Langevin equation (GLE) is widely used to model complex soft-matter systems, including biomolecular dynamics, by incorporating memory effects and colored noise into coarse-grained descriptions. However, recent results suggest that combining memory with non-linear forces, which are ubiquitous in soft matter, introduces fundamental analytical inconsistencies. Here, using a simplified model, we investigate the practical numerical consequences of these analytical results in equilibrium systems. We show that non-linear forces generate cross-correlations with the noise, modifying the fluctuation-dissipation theorem and rendering the noise position-dependent and irreversible. This implies that the commonly assumed reversible Gaussian noise in GLE simulations fails to capture essential features of the microscopic fluctuations. For weak non-linearities, these issues can be partially resolved either by using an iterative optimization of memory or by using microscopically consistent noise, which unexpectedly synchronizes GLE trajectories with the underlying microscopic dynamics. For stronger non-linearities, such as high barriers or shoulders in the external potential, however, iterative reconstruction fails and we observe desynchronization, indicating that the non-linear GLE no longer correctly reproduces the microscopic dynamics. Our results show in which situations non-linear GLEs can be accurately applied and when they fail, thus providing practical guidance for their application to coarse-grain soft-matter systems.
Jung et al. (2026) studied this question.