Efficient Parallel Algorithm for Statistical Ion Track Simulations in Crystalline Materials

We present an efficient parallel algorithm for statistical Molecular Dynamics simulations of ion tracks in solids. The method is based on the Rare Event Enhanced Domain following Molecular Dynamics (REED-MD) algorithm, which has been successfully applied to studies of, e.g., ion implantation into crystalline semiconductor wafers. We discuss the strategies for parallelizing the method, and we settle on a host-client type polling scheme in which a multiple of asynchronous processors are continuously fed to the host, which, in turn, distributes the resulting feed-back information to the clients. This real-time feed-back consists of, e.g., cumulative damage information or statistics updates necessary for the cloning in the rare event algorithm. We finally demonstrate the algorithm for radiation effects in a nuclear oxide fuel, and we show the balanced parallel approach with high parallel efficiency in multiple processor configurations.Comment: 17 pages, seven figures, four table

An Efficient Molecular Dynamics Scheme for the Calculation of Dopant Profiles due to Ion Implantation

We present a highly efficient molecular dynamics scheme for calculating the concentration depth profile of dopants in ion irradiated materials. The scheme incorporates several methods for reducing the computational overhead, plus a rare event algorithm that allows statistically reliable results to be obtained over a range of several orders of magnitude in the dopant concentration. We give examples of using this scheme for calculating concentration profiles of dopants in crystalline silicon. Here we can predict the experimental profile over five orders of magnitude for both channeling and non-channeling implants at energies up to 100s of keV. The scheme has advantages over binary collision approximation (BCA) simulations, in that it does not rely on a large set of empirically fitted parameters. Although our scheme has a greater computational overhead than the BCA, it is far superior in the low ion energy regime, where the BCA scheme becomes invalid.Comment: 17 pages, 21 figures, 2 tables. See: http://bifrost.lanl.gov/~reed