Publications

The kinetic Monte Carlo method and its variants are powerful tools for modeling materials at the mesoscale, meaning at length and time scales in between the atomic and continuum. We have completed a 3 year LDRD project with the goal of developing a parallel kinetic Monte Carlo capability and applying it to materials modeling problems of interest to Sandia. In this report we give an overview of the methods and algorithms developed, and describe our new open-source code called SPPARKS, for Stochastic Parallel PARticle Kinetic Simulator. We also highlight the development of several Monte Carlo models in SPPARKS for specific materials modeling applications, including grain growth, bubble formation, diffusion in nanoporous materials, defect formation in erbium hydrides, and surface growth and evolution.

The authors present the results of molecular dynamics simulations of n-butane and isobutane in silicalite. They begin with a comparison of the bulk adsorption and diffusion properties for two different parameterizations of the interaction potential between the hydrocarbon species, both of which have been shown to reproduce experimental gas-liquid coexistence curves. They examine diffusion as a function of the loading of the zeolite, as well as the temperature dependence of the diffusion constant at loading and for infinite dilution. They continue with simulations in which interfaces are formed between single component gases and the zeolite. After reaching equilibrium, they examine the dynamics of exchange between the bulk gas and the zeolite. Finally, they calculate the permeability of the zeolite for n-butane and isobutane as a function of pressure. Their simulations are performed for a number of different gas temperatures and pressures, covering a wide range of state points.