Publications

Plimpton, Steven J.; Battaile, Corbett C.; Chandross, M.; Holm, Elizabeth A.; Thompson, Aidan P.; Tikare, Veena T.; Webb, Edmund B.; Zhou, Xiaowang Z.

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.

Journal of Physical Chemistry B

Chandross, M.; Webb, Edmund B.; Grest, Gary S.; Martin, Marcus G.; Thompson, Aidan P.; Roth, M.W.

We present the results of Molecular Dynamics and Monte Carlo simulations of n-butane and isobutane in silicalite. We 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. We 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. Both force fields give accurate descriptions of bulk properties. We continue with simulations in which interfaces are formed between single component gases and the zeolite. After reaching equilibrium, we examine the dynamics of exchange between the bulk gas and the zeolite. In particular, we examine the average time spent in the adsorption layer by molecules as they enter the zeolite from the gas in an attempt to probe the microscopic origins of the surface barrier. The microscopic barrier is found to be insignificant for experimental systems. Finally, we calculate the permeability of the zeolite for n-butane and isobutane as a function of pressure. Our results underestimate the experimental results by an order of magnitude, indicating a strong effect from the surface barrier in these simulations. Our simulations are performed for a number of different gas temperatures and pressures, covering a wide range of state points.