Search for: All records

SPPARKS is an open-source parallel simulation code for developing and running various kinds of on-lattice Monte Carlo models at the atomic or meso scales. It can be used to study the properties of solid-state materials as well as model their dynamic evolution during processing. The modular nature of the code allows new models and diagnostic computations to be added without modification to its core functionality, including its parallel algorithms. A variety of models for microstructural evolution (grain growth), solid-state diffusion, thin film deposition, and additive manufacturing (AM) processes are included in the code. SPPARKS can also be used to implement grid-based algorithms such as phase field or cellular automata models, to run either in tandem with a Monte Carlo method or independently. For very large systems such as AM applications, the Stitch I/O library is included, which enables only a small portion of a huge system to be resident in memory. In this paper we describe SPPARKS and its parallel algorithms and performance, explain how new Monte Carlo models can be added, and highlight a variety of applications which have been developed within the code.

We demonstrate using theoretical, computational, and experimental studies a morphological instability, in which a polycrystalline nanorod breaks up at grain boundaries (GBs) into an array of isolated domains. Our theoretical model is used to establish a neutral stability surface demarcating stable and unstable perturbations. It is shown that GBs play a destabilizing role in which the critical wavelength for the instability decreases with the increase in the GB energy. We carry out phase field simulations, which reveal accelerated pinch-off kinetics with the increase in the GB energy and predict temporal evolution of interfacial profiles in quantitative agreement with experimental observations.

The interaction of alloying elements with grain boundaries (GBs) influences many phenomena, such as microstructural evolution and transport. While GB solute segregation has been the subject of active research in recent years, most studies focus on ground-state GB structures, i.e., lowest energy GBs. The impact of GB metastability on solute segregation remains poorly understood. Herein, we leverage atomistic simulations to generate metastable structures for a series of [001] and [110] symmetric tilt GBs in a model Al–Mg system and quantify Mg segregation to individual sites within these boundaries. Our results show large variations in the atomic Voronoi volume due to GB metastability, which are found to influence the segregation energy. The atomistic data are then used to train a Gaussian Process machine learning model, which provides a probabilistic description of the GB segregation energy in terms of the local atomic environment. In broad terms, our approach extends existing GB segregation models by accounting for variability due to GB metastability, where the segregation energy is treated as a distribution rather than a single-valued quantity.