Authors: Manish Prasad, Talid Sinno
Publication Date: August 13, 2004
Journal: Journal of Chemical Physics
Citation: M. Prasad and T. Sinno, Feature-Activated Molecular Dynamics: An Efficient Approach for Atomistic Simulation of Solid-State Aggregation Phenomena, Journal of Chemical Physics, 121 (2004) 8699-8710. Selected to be included in the Virtual Journal of Biological Physics Research, 8, (2004).
Abstract: An efficient approach is presented for performing efficient molecular dynamics simulations of solute aggregation in crystalline solids. The method dynamically divides the total simulation space into “active” regions centered about each minority species, in which regular molecular dynamics is performed. The number, size, and shape of these regions is updated periodically based on the distribution of solute atoms within the overall simulation cell. The remainder of the system is essentially static except for periodic rescaling of the entire simulation cell in order to balance the pressure between the isolated molecular dynamics regions. The method is shown to be accurate and robust for the Environment-Dependant Interatomic Potential (EDIP) for silicon and an Embedded Atom Method potential (EAM) for copper. Several tests are performed beginning with the diffusion of a single vacancy all the way to large-scale simulations of vacancy clustering. In both material systems, the predicted evolutions agree closely with the results of standard molecular dynamics simulations. Computationally, the method is demonstrated to scale almost linearly with the concentration of solute atoms, but is essentially independent of the total system size. This scaling behavior allows for the full dynamical simulation of aggregation under conditions that are more experimentally realizable than would be possible with standard molecular dynamics.