Journal article
Journal of Chemical Physics, 2026
APA
Click to copy
Neuman, V., Wesołowski, P., Bojarski, K., Dewan, D., Schäffler, M., Smardz, P., & Wales, D. J. (2026). Visualizing the energy landscape for a molecular dynamics trajectory. Journal of Chemical Physics.
Chicago/Turabian
Click to copy
Neuman, V., P. Wesołowski, K. Bojarski, Diksha Dewan, Moritz Schäffler, Pamela Smardz, and D. J. Wales. “Visualizing the Energy Landscape for a Molecular Dynamics Trajectory.” Journal of Chemical Physics (2026).
MLA
Click to copy
Neuman, V., et al. “Visualizing the Energy Landscape for a Molecular Dynamics Trajectory.” Journal of Chemical Physics, 2026.
BibTeX Click to copy
@article{v2026a,
title = {Visualizing the energy landscape for a molecular dynamics trajectory.},
year = {2026},
journal = {Journal of Chemical Physics},
author = {Neuman, V. and Wesołowski, P. and Bojarski, K. and Dewan, Diksha and Schäffler, Moritz and Smardz, Pamela and Wales, D. J.}
}
We introduce an open-source program that converts molecular dynamics trajectories into disconnectivity graphs, providing a concise and interpretable visualization of the energy landscape that has been traversed. Our approach applies Savitzky-Golay smoothing to per-frame thermodynamic traces (potential energy in NVE/NVT or enthalpy in NPT ensembles) to identify local extrema as proxies for minima and transition states, and generates the necessary files for disconnectivity graph construction. This workflow requires no additional geometry optimization. The method is ensemble-agnostic and compatible with both all-atom and coarse-grained simulations. For some representative biochemical systems, it processes 104-105 frames in seconds on a standard laptop and produces an approximate representation of the underlying landscape topology. The resulting graphs capture the structures that are visited and pathways between them for a selected energy and time resolution, offering an interpretable structural summary of conformational hierarchies with minimal postprocessing. Because extrema are detected directly from the trajectory, the graphs reflect the organization of the explored region of the landscape on the molecular dynamics timescale. Our approach basically substitutes local minima and maxima from the smoothed time series as proxies for the true stationary points of the underlying landscape.