Thomas Danielson

First-principles investigation of helium in Y2O3

Density functional theory has been used to investigate the relative stability, migration barriers and potential energy surfaces of helium in Y2O3 in order to provide insight to the trapping mechanisms of helium at the oxide-Fe interface in nanostructured ferritic alloys. The helium interstitial sites have been located and their relative stability in Y2O3 is determined from the solution energy. The interaction of helium with the oxide has been investigated using the charge density and electron localization function. The migration barriers have been calculated and further investigated using potential energy surfaces computed on various lattice planes. Helium trapping has been attributed to very large diffusion barrier energies, especially close to oxygen atoms.

SQERTSS: Dynamic rank based throttling of transition probabilities in kinetic Monte Carlo simulations

Abstract Lattice based Kinetic Monte Carlo (KMC) simulations offer a powerful simulation technique for investigating large reaction networks while retaining spatial configuration information, unlike ordinary differential equations. However, large chemical reaction networks can contain reaction processes with rates spanning multiple orders of magnitude. This can lead to the problem of “KMC stiffness” (similar to stiffness in differential equations), where the computational expense has the potential to be overwhelmed by very short time-steps during KMC simulations, with the simulation spending an inordinate amount of KMC steps/CPU time simulating fast frivolous processes (FFPs) without progressing the system (reaction network). In order to achieve simulation times that are experimentally relevant or desired for predictions, a dynamic throttling algorithm involving separation of the processes into speed-ranks based on event frequencies has been designed and implemented with the intent of decreasing the probability of FFP events, and increasing the probability of slow process events—allowing rate limiting events to become more likely to be observed in KMC simulations. This Staggered Quasi-Equilibrium Rank-based Throttling for Steady-state (SQERTSS) algorithm is designed for use in achieving and simulating steady-state conditions in KMC simulations. As shown in this work, the SQERTSS algorithm also works for transient conditions: the correct configuration space and final state will still be achieved if the required assumptions are not violated, with the caveat that the sizes of the time-steps may be distorted during the transient period.

Generalized adsorption isotherms for molecular and dissociative adsorption of a polar molecular species on two polar surface geometries: Perovskite (100) (Pm-3m) and fluorite (111) (Fm-3m)

Lattice based kinetic Monte Carlo simulations have been used to determine a functional form for the second order adsorption isotherms on two commonly investigated crystal surfaces: the (111) fluorite surface and the (100) perovskite surface which has the same geometric symmetry as the NaCl (100) surface. The functional form is generalized to be applicable to all values of the equilibrium constant by a shift along the pressure axis. Functions have been determined for estimating the pressure at which a desired coverage would be achieved and, conversely, for estimating the coverage at a certain pressure. The generalized form has been calculated by investigating the surface adsorbate coverage across a range of thermodynamic equilibrium constants that span the range 10−26 to 1013. The equations have been shown to be general for any value of the adsorption equilibrium constant.

The F-t-Pj-RG method: An adjacent-rolling-windows based steady-state detection technique for application to kinetic Monte Carlo simulations

Abstract A window-based steady-state detection algorithm has been developed for application to kinetic Monte Carlo simulation data. The algorithm, termed F-t-Pj-RG sequentially applies an F-test, a t-test, and a projection test on adjacent windows of the data while rolling (or shifting) and growing the windows when any of the tests fail. In aggregate, the algorithm is able to (a) automatically reject the warm-up period as not being at steady-state, as well as (b) determine an appropriate window size for converged statistics when sampling the data, which is necessary for detection of steady-state, and (c) detect steady-state within a particular tolerance. The last step, the projection test, is actually an oscillating-slope projection test, and is performed on j sequential data windows (i.e., more than two adjacent windows). It requires more than simply being within the user defined tolerance: the oscillating-slope projection test includes a condition that the slope must oscillate around zero when ≥ 2, which is an additional indication of steady-state. When all three tests are passed, the F-t-Pj test is passed, indicating that the prerequisites of steady-state detection have been met and also that conditions consistent with the definition of steady-state have been realized. This algorithm is applied to a variety of data sets that correspond to the diverse type of data trends that can be produced by kinetic Monte Carlo simulations. The algorithm is shown to be robust in its ability to handle differing functional forms, and is able to detect steady-state with low computational cost. The low computational cost of this method and its robustness towards varied data trends make it suitable for on-the-fly use in kinetic Monte Carlo simulations.

Below-Room-Temperature C–H Bond Breaking on an Inexpensive Metal Oxide: Methanol to Formaldehyde on CeO2(111)

Upgrading of primary alcohols by C-H bond breaking currently requires temperatures of >200 °C. In this work, new understanding from simulation of a temperature-programmed reaction study with methanol over a CeO2(111) surface shows C-H bond breaking and the subsequent desorption of formaldehyde, even below room temperature. This is of particular interest because CeO2 is a naturally abundant and inexpensive metal oxide. We combine density functional theory and kinetic Monte Carlo methods to show that the low-temperature C-H bond breaking occurs via disproportionation of adjacent methoxy species. We further show from calculations that the same transition state with comparable activation energy exists for other primary alcohols; with ethanol, 1-propanol, and 1-butanol explicitly calculated. These findings indicate a promising class of transition states to search for in seeking low-temperature C-H bond breaking over inexpensive oxides.

Investigation of helium at a Y2Ti2O7 nanocluster embedded in a BCC Fe matrix

Nanostructured ferritic alloys (NFAs) are prime candidates for structural and first wall components of fission and fusion reactors. The main reason for this is their ability to effectively withstand high concentrations of the transmutation product helium. A high number density of oxide nanoclusters dispersed throughout a BCC Fe matrix act as trapping sites for helium and prevent its eventual delivery to high risk nucleation sites. The current study uses density functional theory to investigate the helium trapping mechanisms at the boundary between BCC iron and Y2Ti2O7, a common stoichiometry of the oxide nanoclusters in NFAs. The investigation is carried out on a structure matched oxide nanocluster that is embedded within a BCC Fe supercell. Investigation of the electronic structure and a mapping of the potential energy landscape reveals that the localized iono-covalent bonds present within the oxides create a potential energy-well within the metallically bonded BCC Fe matrix, so that trapping of helium at the oxide nanocluster is thermodynamically and kinetically favorable.