On the calculation of time correlation functions by potential scaling

We present and analyze a general method to calculate time correlation functions from molecular dynamics on scaled potentials for complex systems for which simulation is affected by broken ergodicity. Depending on the value of the scaling factor, correlations can be calculated for times that can be orders of magnitude longer than those accessible to direct simulations. We show that the exact value of the time correlation functions of the original system (i.e., with unscaled potential) can be obtained, in principle, using an action-reweighting scheme based on a stochastic path-integral formalism. Two tests (involving a bistable potential model and a dipeptide bond-vector orientational relaxation) are exemplified to showcase the strengths, as well as the limitations of the approach, and a procedure for the estimation of the time-dependent standard deviation error is outlined.