Averaging Molecular Dynamics simulations to study the slow-strain rate behavior of metals

TL;DR

This paper introduces the Practical Time Averaging framework to enable atomistic simulations of metals under quasistatic loading, capturing dislocation microstructures and size effects efficiently at low strain rates.

Contribution

The study develops a novel PTA method that overcomes MD timescale limitations, allowing realistic quasistatic deformation simulations of crystalline solids with full atomistic detail.

Findings

Stress-strain curves show yielding near theoretical stress for nucleation.

Size effects influence the hardness and serration prominence.

Method captures dislocation microstructure evolution on slow timescales.

Abstract

The application of molecular dynamics (MD) simulations to quasistatic loading is severely limited by the large separation between atomic vibration timescales and experimentally relevant deformation rates. In this work we employ the Practical Time Averaging (PTA) framework to overcome this limitation and enable atomistic simulations of crystalline solids under quasistatic loading conditions. PTA exploits the intrinsic separation of timescales by defining slow variables as time-averaged observables of the fast atomistic dynamics and their evolution on the slow loading timescale, thereby avoiding explicit integration of the fast dynamics. Using this approach, we simulate uniaxial deformation, in both tension and compression, of 4 to 20 nm cubic specimens of face centered cubic aluminum nanocrystals at applied strain rates approaching quasistatic conditions. We define slow variables as the averaged kinetic energy, potential energy, and normal stress in the loading direction, and track their evolution on the slow timescale. The stress-strain curves show yielding close to the theoretical stress for homogeneous nucleation, followed by successive load drops and rises caused by dislocation nucleation, motion, and exit from free surfaces. The "smaller is harder" effect is evident from the stress-strain response and from the variation of yield stress with sample size. Serrations in the response are more pronounced for smaller samples. The effects of applied strain rate and initial temperature are also studied. The method also captures the evolution of intricate dislocation microstructures on the slow timescale by tracking time-averaged atomic positions. The PTA framework enables simulations at strain rates several orders of magnitude lower than those accessible to conventional MD, demonstrating significant speedup in computational time while retaining full atomistic resolution.