Multidimensional Langevin Modeling of Nonoverdamped Dynamics

Based on a given time series, data-driven Langevin modeling aims to construct a low-dimensional dynamical model of the underlying system. When dealing with physical data as provided by, e.g., all-atom molecular dynamics simulations, effects due to small damping may be important to correctly describe the statistics (e.g., the energy landscape) and the dynamics (e.g., transition times). To include these effects in a dynamical model, an algorithm that propagates a second-order Langevin scheme is derived, which facilitates the treatment of multidimensional data. Adopting extensive molecular dynamics simulations of a peptide helix, a five-dimensional model is constructed that successfully forecasts the complex structural dynamics of the system. Neglect of small damping effects, on the other hand, is shown to lead to significant errors and inconsistencies.