A Caltech Library Service

Acceleration of diffusive molecular dynamics simulations through mean field approximation and subcycling time integration

Sun, X. and Ariza, M. P. and Ortiz, M. and Wang, K. G. (2017) Acceleration of diffusive molecular dynamics simulations through mean field approximation and subcycling time integration. Journal of Computational Physics, 350 . pp. 470-492. ISSN 0021-9991.

Full text is not posted in this repository. Consult Related URLs below.

Use this Persistent URL to link to this item:


Diffusive Molecular Dynamics (DMD) is a class of recently developed computational models for the simulation of long-term diffusive mass transport at atomistic length scales. Compared to previous atomistic models, e.g., transition state theory based accelerated molecular dynamics, DMD allows the use of larger time-step sizes, but has a higher computational complexity at each time-step due to the need to solve a nonlinear optimization problem at every time-step. This paper presents two numerical methods to accelerate DMD based simulations. First, we show that when a many-body potential function, e.g., embedded atom method (EAM), is employed, the cost of DMD is dominated by the computation of the mean of the potential function and its derivatives, which are high-dimensional random variables. To reduce the cost, we explore both first- and second-order mean field approximations. Specifically, we show that the first-order approximation, which uses a point estimate to calculate the mean, can reduce the cost by two to three orders of magnitude, but may introduce relatively large error in the solution. We show that adding an approximate second-order correction term can significantly reduce the error without much increase in computational cost. Second, we show that DMD can be significantly accelerated through subcycling time integration, as the cost of integrating the empirical diffusion equation is much lower than that of the optimization solver. To assess the DMD model and the numerical approximation methods, we present two groups of numerical experiments that simulate the diffusion of hydrogen in palladium nanoparticles. In particular, we show that the computational framework is capable of capturing the propagation of an atomically sharp phase boundary over a time window of more than 30 seconds. The effects of the proposed numerical methods on solution accuracy and computation time are also assessed quantitatively.

Item Type:Article
Related URLs:
URLURL TypeDescription
Ariza, M. P.0000-0003-0266-0216
Ortiz, M.0000-0001-5877-4824
Additional Information:© 2017 Elsevier Inc. Received 8 June 2017, Revised 10 August 2017, Accepted 31 August 2017, Available online 11 September 2017. The authors gratefully acknowledge the support of the Institute for Critical Technologies and Applied Sciences (ICTAS) at Virginia Tech through a Junior Faculty Collaboration (JFC) project, the Ministerio de Economía y Competitividad of Spain under grant number DPI2015-66534-R, and the U.S. Army Research Laboratory (ARL) through the Materials in Extreme Dynamic Environments (MEDE) Collaborative Research Alliance (CRA) under Award Number W911NF-11-R-0001.
Funding AgencyGrant Number
Ministerio de Economía y Competitividad (MINECO)DPI2015-66534-R
Army Research Office (ARO)W911NF-11-R-0001
Subject Keywords:Diffusive molecular dynamics; Mean field approximation; Subcycling time integration; Long-term processes; Hydrogen diffusion; Phase transformation
Record Number:CaltechAUTHORS:20170912-141252063
Persistent URL:
Official Citation:X. Sun, M.P. Ariza, M. Ortiz, K.G. Wang, Acceleration of diffusive molecular dynamics simulations through mean field approximation and subcycling time integration, Journal of Computational Physics, Volume 350, 1 December 2017, Pages 470-492, ISSN 0021-9991, (
Usage Policy:No commercial reproduction, distribution, display or performance rights in this work are provided.
ID Code:81373
Deposited By: Tony Diaz
Deposited On:12 Sep 2017 21:28
Last Modified:09 Mar 2020 13:18

Repository Staff Only: item control page