AbstractA majority of computational mechanical analyses of nanocrystalline materials have been carried out using classical molecular dynamics (MD). Due to the fundamental reason that the MD simulations must resolve atomic level vibrations, they cannot be carried out at the timescale of the order of microseconds. Additionally, MD simulations have to be carried out at very high loading rates (∼108 s−1) rarely observed in experiments. In this investigation a modified Hybrid Monte Carlo (HMC) method that can be used to analyze time-dependent (strain rate dependent) atomistic mechanical deformation of nanocrystalline structures at higher timescales than currently possible using MD is established. In this method there is no restriction on the size of MD timestep except that it must be such that to ensure a reasonable acceptance rate between consecutive Monte-Carlo (MC) time-steps. For the purpose of comparison HMC analyses of a nanocrystalline Ni sample at a strain rate of 109 s−1 with three different timesteps, viz. 2 fs, 4fs, and 8 fs, are compared with the analyses based on MD simulations at the same strain rate and a MD timestep of 2 fs. MD simulations of nanocrystalline Ni reproduce the defect nucleation and propagation results as well as strength values reported in the literature. In addition, HMC with timestep of 8 fs correctly reproduces defect formation and stress-strain response observed in the case of MD simulations with permissible timestep of 2 fs (for the interatomic potential used 2 fs is the highest MD timestep). Simulation time analyses show that by using HMC a saving of the order of 4 can be achieved bringing the atomistic analyses closer to the continuum timescales.
Application of classical simulations for the computation of vibrational properties of free molecules