Fep1d: A Script for the Analysis of Reaction Coordinates

The dynamics of complex systems with many degrees of freedom can be analyzed by projecting it onto one or few coordinates (collective variables). The dynamics is often described then as diffusion on a free energy landscape associated with the coordinates. Fep1d is a script for the analysis of such one-dimensional coordinates. The script allows one to construct conventional and cut-based free energy profiles, to assess the optimality of a reaction coordinate, to inspect whether the dynamics projected on the coordinate is diffusive, to transform (rescale) the reaction coordinate to more convenient ones, and to compute such quantities as the mean first passage time, the transition path times, the coordinate dependent diffusion coefficient, and so forth. Here, we describe the implemented functionality together with the underlying theoretical framework.