Diffusive molecular dynamics is a novel model for materials incorporating atomistic resolution and reaching diffusive time scales. The main ideas of diffusive molecular dynamics are to first minimize an approximate variational Gaussian free energy of the system with respect to the mean atomic coordinates (averaging over many vibrational periods), and to then perform a diffusive step where atoms and vacancies (or two species in a binary alloy) flow on a diffusive time scale via a master equation. We present a mathematical framework for studying this algorithm based on relative entropy, also known as the Kullback-Leibler divergence. This adds flexibility in how the algorithm is implemented and interpreted. We then compare our formulation, relying on relative entropy and absolute continuity of measures, to existing formulations and find agreement.
Bibliographical noteFunding Information:
The first and second authors were supported by U.S. Department of Energy award DE-SC0012733.
© 2016 Society for Industrial and Applied Mathematics.
- Molecular dynamics
- Relative entropy