SPPARKS WWW Site

Pictures & Movies from SPPARKS Simulations

hafnium dioxide hafnium dioxide surface growth
abnormal abnormal grain growth
erbium hydrides diffusion in Erbium hydrides
thin film growth thin film growth
nanoporous metals diffusion in nanoporous metals
nuclear fuel sintering microstructural evolution due to sintering in a nuclear fuel
nuclear fuel gas transport fission gas transport in a nuclear fuel

Pictures on this page are shown in small size. Click on the image to view a larger version.

For movies on this page, click on the small image to trigger the animation.


Hafnium dioxide surface growth

This is work by Mahdi Shirazi (mahdi.shirazi at tyndall.ie) and Simon Elliott (simon.elliott at tyndall.ie) at the Tyndall National Institute in Cork, Ireland. Movie editing was by Kevin Berwick(kevin.berwick at dit.ie) of the Dublin Institute of Technology. The movie shows the the growth of Hafnium dioxide by Atomic Layer Deposition (ALD). Hafnium dioxide is used in microelectronics as a high-k dielectric in modern transistors.

In this study, activation energies, which affect the growth rate and the time evolution of the film, have been obtained by ab initio calculations using Density Functional Theory (DFT). Multiple proton diffusion, steric effects, structural relaxation and cooperative effects are included for the first time in an ALD growth simulation, using SPPARKS. The movie shows 40 growth cycles, simulated on the Stokes supercomputer at the Irish Centre for High-End Computing.

Image and a 36 Mb MPG movie

This paper describes the DFT methodology used to calculate kMC input parameters:

Multiple Proton Diffusion and Film Densification in Atomic Layer Deposition Modeled by Density Functional Theory, M. Shirazi and S. D. Elliott, Chemistry of Materials, 25, 878-889 (2013). (abstract)


Abnormal grain growth

This is work by Liz Holm (eaholm at sandia.gov) at Sandia and Todd Hoffmann (todd.hoffmann at asu.edu) and his advisor at Carnegie Mellon, to study the mechanisms and statistics associated with abnormal grain growth. In grain growth, when the system is annealed, grain boundaries will move and rearrange so as to increase the average grain size and decrease the grain boundary area per unit volume. In normal grain growth, this occurs relatively uniformly, but in abnormal grain growth, one grain grows at the expense of the others. Abnormal grain growth can have many causes. Here the effects of static particle dispersions on grain growth are studied. Particles should pin grain boundaries in positions that maximize their contact with the particles, so that grain growth cannot proceed. However, occasionally one grain breaks free and is able to grow at the expense of others, leading to abnormal growth.

The on-lattice SPPARKS application (app_style potts/pin) that runs this model implements a Potts model of grain growth and also allows for inert pinning particles. The pinning particles can be inserted preferentially at grain boundaries or throughtout the system. Several hundred independent runs were made on 300x300x300 lattices, each of at least a million Monte Carlo steps (sweep over all lattice sites) at finite simulation temperature, to acquire enough data for statistical analysis.

The images and movies show typical results for systems with 5 volume percent pinning particles, inserted when the average grain radius was 10 sites. In the normal case, the system becomes fully pinned. In the second case, abnormal growth eventually occurs.

Two images and two 1.5 Mb AVI movies


Diffusion in Erbium hydrides

This is work by Corbett Battaile (ccbatta at sandia.gov) and Steve Plimpton at Sandia, to model the mobility of gaseous species such as tritium and helium in metal hydrides, like erbium. Such hydrides are commonly used for neutron generation in many applications in industry, medicine, and national security, including petroleum exploration, explosives detection, toxic waste analysis, and boron neutron capture therapy. As tritium decays into helium, the helium gas clusters to form bubbles, often found to be preferentially oriented along {111} planes in the erbium lattice.

The on-lattice SPPARKS application (app_style erbium) that runs this model uses an fcc lattice with additional tetrahedral and octrahedral interstitial sites. Both the tritium and helium can diffuse through the lattice on the interstitial sites.

The image is a snapshot from a simulation on 1024 processors of a 500x500x500 unit-cell model (2 billion sites) which encodes the list of 1-, 2-, and 3-site migration mechanisms for the two mobile species. The helium (red) diffuses but is not yet aggregating due to a lack of appropriate binding and barrier energies for the various event configurations. These values (inputs to SPPARKS) are currently being calculated via density functional theory in the electronic structure code SeqQuest.


Thin film growth

This is work by Ed Webb (ebwebb at sandia.gov), Xiaowang Zhou (xzhou at sandia.gov), and Steve Plimpton (sjplimp at sandia.gov) at Sandia, to model thin film deposition and growth. For many thin film applications, a high degree of control is desired over film properties including surface/interface roughness, film texture, grain structure, and defect content. For example, controlling void content is critical to a number of applications where voids created during deposition lead to shortened device lifetimes. Interface roughness is also important to many photovoltaic, microelectronic, and optical thin film applications when multilayered films are deposited. On-lattice KMC models can be used to study many of these effects, though they cannot capture effects like dislocation formation in the depositing crystal.

The on-lattice SPPARKS application (app_style diffusion) which ran these Ni models uses an underlying triangular lattice (2d) or fcc lattice (3d) and defines two types of events (deposition and diffusion). Nearest-neighbor and/or 2nd nearest-neighbor (Schwoebel) hops can be enabled or disabled. In these simulations, all deposition was performed at normal incidence.

The first 3 images and accompanying movies are of 2d growth under 3 conditions. The first is with no Schwoebel hops (which might be suppressed experimentally by adding impurities to the system) and results in a thin film with high porosity. The 2nd and 3rd are with Schwoebel hops at low and high temperatures (250 and 600K). In the low-temperature case there are about 3 diffusive hops per deposition event. In the high-temperature case, the ratio is 72000, resulting in a much smoother film.

The next 3 images are 3d growth (on a 2d surface), again under 3 conditions. The first two are with no Schwoebel hops at low and high temperatures (300 and 400K). The third is with Schwoebel hops at 300K. Similar effects as in the 2d case are observed, though the system is harder to visualize.

Finally, a rough surface (sinusoidal hillocks) was annealed at 500K for 10,000 seconds with a diffusion-only model in a fcc model, both with and without Schwoebel events. In both case the unstable surface asperities (top of figure) turned 45 degrees to a (111) orientation to maximize the surface stablility. However, without Schwoebel hops (middle of figure), the system lacks the ability to anneal out terraces. With Schwoebel hops (bottom of figure), terraces were all annealed out and clear (111) and (100) facets developed.

Three images and three 10 Mb QuickTime movies of 2d thin film growth. Three images of 3d film growth. One image (yellow atoms) of hillock rearrangement on a 2d surface.


Diffusion in nanoporous metals

This is work by Greg Wagner (gjwagne at sandia.gov) at Sandia, to model diffusion in nanoporous metals. Such materials have a range of applications based on their very large surface-to-volume ratios. At Sandida nanoporous palladium particles are being studied for use as a hydrogen storage medium, and nanoporous gold wires are being tested for use as electrodes in super-capacitors. As these complex materials age, especially at elevated temperatures, the evolution of their nanostructures with time can greatly affect material properties and performance. KMC can model this aging over long timescales.

The on-lattice SPPARKS application (app_style diffusion) that runs this model uses an fcc lattice and energies and barrier heights calcaluted from atomistic simulations using embedded atom method (EAM) potentials. Diffusion in the KMC model is described by Kawasaki dynamics, in which each event is the motion of an individual atom. A basic atom hop is defined as movement to any vacant neighboring site. Special consideration is given to Schwoebel barriers, in which an atom diffuses up or down a ledge between atomic monolayers, which is effectively a hop to a lattice location two sites away.

The images and associated movies are a bulk diffusion model with 60x60x60 unit cells (0.864 million sites) and a 50% porosity factor. If we define the length scale of the system as the inverse of the surface-to-volume ratio, then we observe the fourth power of this quantity grows linearly with time, in agreement with theory.

Two images and two QuickTime movies (6 and 2 Mb)


Microstructural evolution due to sintering in a nuclear fuel

This is work by Veena Tikare (vtikare at sandia.gov) at Sandia and Cristina Garcia Cardona (cgarcia at sciences.sdsu.edu) at San Diego State Univ, to simulate microstructural evolution during sintering. They augmented a Potts model of grain growth with three additional physical effects: surface diffusion, vacancy diffusion, and vacancy annihilation. The implementation of the annihilation algorithm in SPPARKS was non-trivial as an internal vacancy is annihilated by collapsing all the sites along a line from it to an external surface. In parallel this involves moving all the sites along that line towards the vacancy, requiring interprocessor communication. This sintering model is able to capture all the necessary mechanism to simulate simple solid-state sintering correctly. It was demonstrated by comparing it to 3D, in-situ images taken in a high-energy synchrotron during the sintering of Cu particles. The images below are from a simulation of the sintering process for a powder of Cu particles. The first pair of images are for an initial material density of 69%. As it compacts, an intermediate density of 79% is shown. The final pair of images are for a compacted density at the end of the simulation of 83%. Both 3d pictures of the full model and 2d slices are shown. The density versus time profiles, microstructural images and grain size distributions for this system are in good agreement with experiment.

This paper gives more details of the model and simulations:

Parallel Simulation of 3D Sintering, C. G. Cardona, V. Tikare, S. J. Plimpton, Int J Computational Materials Science and Surface Engineering, 4, 37-54 (2011). (abstract)


Fission gas transport in a nuclear fuel

This is work by Veena Tikare (vtikare at sandia.gov) at Sandia to study fission gas bubble transport in Light Water Reactor (LWR) fuel. During operation, Xe and Kr gas atoms are produced as fission by-products. They are highly insoluble in the fuel matrix and precipitate to form two types of bubbles, small (~10 nm) intragranular bubbles and larger (~500 nm) intergranular bubbles. As the gas accumulates, it eventually finds it way to grain boundaries by diffusion and the intergranular bubbles grow, initially by gas accumulation and later by coalescence. Eventually coalescence of the intergranular bubbles progress to the point that they percolate, releasing the pressurized gas to the external surface. The bubbles sinter back and the process repeats itself. This process of percolation and gas release has important implications for fuel performance and longevity. To model these effects, a Potts model of grain growth was augmented to incorporate gas generation, gas atom diffusion, precipitation to form bubbles, and surface diffusion of the gas bubbles. This model was used to study the formation, growth and percolation of the intragranular gas bubbles in LWR fuels. Images from such a simulation are shown below; the first image shows only the gas bubbles in a two-grain microstructure with periodic boundary conditions. The inter- and intra-granular bubbles are clearly shown. The next 2 images show gas bubble growth at the planar grain boundary during the simulation. The final image shows the percolating bubble, in dark blue, at the end of the simulation.