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

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.

[1]  L. Verlet.  Computer "Experiments" on Classical Fluids. I. Thermodynamical Properties of Lennard-Jones Molecules , 1967 .

[2]  Bryan M. Wong,et al.  An embedded-atom method interatomic potential for Pd–H alloys , 2008 .

[3]  Kenichi Kato,et al.  Shape-dependent hydrogen-storage properties in Pd nanocrystals: which does hydrogen prefer, octahedron (111) or cube (100)? , 2014, Journal of the American Chemical Society.

[4]  Serge Piperno,et al.  Explicit/implicit fluid/structure staggered procedures with a structural predictor and fluid subcycling for 2D inviscid aeroelastic simulations , 1997 .

[5]  Yu-Ming Lin,et al.  An integrated purification and production of hydrogen with a palladium membrane-catalytic reactor , 1998 .

[6]  J. B. Adams,et al.  Modelling and Simulation in Materials Science and Engineering EAM potential for magnesium from quantum mechanical forces , 1996 .

[7]  J. Dionne,et al.  Reconstructing solute-induced phase transformations within individual nanocrystals. , 2016, Nature materials.

[8]  Kenichi Kato,et al.  Hydrogen storage in Pd nanocrystals covered with a metal-organic framework. , 2014, Nature materials.

[9]  D. Gudbjartsson,et al.  Variants in ELL2 influencing immunoglobulin levels associate with multiple myeloma , 2015, Nature Communications.

[10]  N. Takagi,et al.  Path and mechanism of hydrogen absorption at Pd(100) , 1998 .

[11]  Arthur F. Voter,et al.  Introduction to the Kinetic Monte Carlo Method , 2007 .

[12]  Robert W. Balluffi,et al.  Kinetics Of Materials , 2005 .

[13]  R. Miron,et al.  Heteroepitaxial growth of Co/Cu(001): An accelerated molecular dynamics simulation study , 2005 .

[14]  S. Baranton,et al.  Octahedral palladium nanoparticles as excellent hosts for electrochemically adsorbed and absorbed hydrogen , 2017, Science Advances.

[15]  P Zapol,et al.  Avalanching strain dynamics during the hydriding phase transformation in individual palladium nanoparticles , 2015, Nature Communications.

[16]  S. Plimpton,et al.  Computational limits of classical molecular dynamics simulations , 1995 .

[17]  M. Ortiz,et al.  Finite-temperature non-equilibrium quasi-continuum analysis of nanovoid growth in copper at low and high strain rates , 2015 .

[18]  Jun Chen,et al.  Magnesium nanowires: enhanced kinetics for hydrogen absorption and desorption. , 2007, Journal of the American Chemical Society.

[19]  Xiaojun Yu,et al.  A comparison of subcycling algorithms for bridging disparities in temporal scale between the fire and solid domains , 2013 .

[20]  Wenquan Lu,et al.  Silicon‐Based Nanomaterials for Lithium‐Ion Batteries: A Review , 2014 .

[21]  R. Miron,et al.  Accelerated molecular dynamics with the bond-boost method , 2003 .

[22]  David J. Srolovitz,et al.  Thermodynamics of solid and liquid embedded‐atom‐method metals: A variational study , 1991 .

[23]  Kevin G. Wang,et al.  Long-term atomistic simulation of hydrogen diffusion in metals , 2015 .

[24]  Astrid Pundt,et al.  Hydrogen in Nano‐sized Metals , 2004 .

[25]  Charles M. Lieber,et al.  Nanoelectronics from the bottom up. , 2007, Nature materials.

[26]  Kevin G. Wang,et al.  DEFORMATION-DIFFUSION COUPLED ANALYSIS OF LONG-TERM HYDROGEN DIFFUSION IN NANOFILMS , 2016 .

[27]  O. Custance,et al.  Surface diffusion of single vacancies on Ge(111)-c(2×8) studied by variable temperature scanning tunneling microscopy , 2004 .

[28]  Michael Ortiz,et al.  A variational approach to coarse-graining of equilibrium and non-equilibrium atomistic description at finite temperature , 2007 .

[29]  D. Perez,et al.  Multiscale model for microstructure evolution in multiphase materials: Application to the growth of isolated inclusions in presence of elasticity. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[30]  A. Voter,et al.  Extending the Time Scale in Atomistic Simulation of Materials Annual Re-views in Materials Research , 2002 .

[31]  Ray,et al.  Pressure-composition isotherms for palladium hydride. , 1993, Physical review. B, Condensed matter.

[32]  N. Metropolis,et al.  The Monte Carlo method. , 1949 .

[33]  J. Dionne,et al.  In situ detection of hydrogen-induced phase transitions in individual palladium nanocrystals. , 2014, Nature materials.

[34]  O. Salata,et al.  Applications of nanoparticles in biology and medicine , 2004, Journal of nanobiotechnology.

[35]  J. Nocedal.  Updating Quasi-Newton Matrices With Limited Storage , 1980 .

[36]  Alvaro L. G. A. Coutinho,et al.  On decoupled time step/subcycling and iteration strategies for multiphysics problems , 2008 .

[37]  M. Baskes,et al.  Modified embedded-atom potentials for cubic materials and impurities. , 1992, Physical review. B, Condensed matter.

[38]  Yunzhi Wang,et al.  Finding activation pathway of coupled displacive-diffusional defect processes in atomistics: Dislocation climb in fcc copper , 2012 .

[39]  M. Baskes,et al.  Embedded-atom method: Derivation and application to impurities, surfaces, and other defects in metals , 1984 .

[40]  Michael Ortiz,et al.  Atomistic long-term simulation of heat and mass transport , 2014 .

[41]  J. Dionne,et al.  Direct visualization of hydrogen absorption dynamics in individual palladium nanoparticles , 2017, Nature Communications.

[42]  Jörg Rottler,et al.  Solute-defect interactions in Al-Mg alloys from diffusive variational Gaussian calculations , 2014 .

[43]  E. Jaynes.  Information Theory and Statistical Mechanics , 1957 .

[44]  P. Plecháč,et al.  Spin-diffusions and diffusive molecular dynamics , 2017, 1702.01469.

[45]  Mitchell Luskin,et al.  A Theoretical Examination of Diffusive Molecular Dynamics , 2015, SIAM J. Appl. Math..

[46]  J. W. Kim,et al.  Three-dimensional imaging of dislocation dynamics during the hydriding phase transformation. , 2017, Nature materials.

[47]  Yunzhi Wang,et al.  Diffusive molecular dynamics and its application to nanoindentation and sintering , 2011 .

[48]  Gene H. Golub,et al.  Calculation of Gauss quadrature rules , 1967, Milestones in Matrix Computation.