12 research outputs found
Parallel kinetic Monte Carlo simulation framework incorporating accurate models of adsorbate lateral interactions
Ab initio kinetic Monte Carlo (KMC) simulations have been successfully applied for over two decades to elucidate the underlying physico-chemical phenomena on the surfaces of heterogeneous catalysts. These simulations necessitate detailed knowledge of the kinetics of elementary reactions constituting the reaction mechanism, and the energetics of the species participating in the chemistry. The information about the energetics is encoded in the formation energies of gas and surface-bound species, and the lateral interactions between adsorbates on the catalytic surface, which can be modeled at different levels of detail. The majority of previous works accounted for only pairwise-additive first nearest-neighbor interactions. More recently, cluster-expansion Hamiltonians incorporating long-range interactions and many-body terms have been used for detailed estimations of catalytic rate [C. Wu, D. J. Schmidt, C. Wolverton, and W. F. Schneider, J. Catal. 286, 88 (2012)]. In view of the increasing interest in accurate predictions of catalytic performance, there is a need for general-purpose KMC approaches incorporating detailed cluster expansion models for the adlayer energetics. We have addressed this need by building on the previously introduced graph-theoretical KMC framework, and we have developed Zacros, a FORTRAN2003 KMC package for simulating catalytic chemistries. To tackle the high computational cost in the presence of long-range interactions we introduce parallelization with OpenMP. We further benchmark our framework by simulating a KMC analogue of the NO oxidation system established by Schneider and co-workers [J. Catal. 286, 88 (2012)]. We show that taking into account only first nearest-neighbor interactions may lead to large errors in the prediction of the catalytic rate, whereas for accurate estimates thereof, one needs to include long-range terms in the cluster expansion
A Caching Scheme to Accelerate Kinetic Monte Carlo Simulations of Catalytic Reactions
Kinetic Monte Carlo (KMC) simulations have been instrumental in advancing our fundamental understanding of heterogeneously catalyzed reactions, with particular emphasis on structure sensitivity, ensemble effects, and the interplay between adlayer structure and adsorbate-adsorbate lateral interactions in shaping the observed kinetics. Yet, the computational cost of KMC remains high, thereby motivating the development of acceleration schemes that would improve the simulation effciency. We present an exact such scheme, which implements a caching algorithm along with shared-memory parallelization to improve the computational performance of simulations incorporating long-range adsorbate-adsorbate lateral interactions. This scheme is based on caching information about the energetic interaction patterns associated with the products of each possible lattice process (adsorption, desorption, reaction etc). Thus, every time a reaction occurs ("ongoing reaction") it enables fast updates of the rate constants of "affected reactions", i.e. possible reactions in the region of influence of the "ongoing reaction". Benchmarks on KMC simulations of NOx oxidation/reduction, yield acceleration factors of up to 20Ă— when comparing single-thread runs without caching to runs on 16 threads with caching, for simulations with a cluster expansion Hamiltonian that incorporates up to 8th nearest-neighbor interactions
Recommended from our members
Band-structure, optical properties, and defect physics of the photovoltaic semiconductor SnS
SnS is a potential earth-abundant photovoltaic (PV) material. Employing both theory and experiment to assess the PV relevant properties of SnS, we clarify on whether SnS has an indirect or direct band gap and what is the minority carrier effective mass as a function of the film orientation. SnS has a 1.07 eV indirect band gap with an effective absorption onset located 0.4 eV higher. The effective mass of minority carrier ranges from 0.5 m(0) perpendicular to the van der Waals layers to 0.2 m0 into the van der Waals layers. The positive characteristics of SnS feature a desirable p-type carrier concentration due to the easy formation of acceptor-like intrinsic Sn vacancy defects. Potentially detrimental deep levels due to Sn(S) antisite or S vacancy defects can be suppressed by suitable adjustment of the growth condition towards S-rich. (C) 2012 American Institute of Physics. [doi:10.1063/1.3675880]Keywords: Solar cell, Fil
Chemical Descriptors of Yttria-Stabilized Zirconia at Low Defect Concentration: An <i>ab Initio</i> Study
Yttria-stabilized zirconia (YSZ)
is an important oxide ion conductor with applications in solid oxide
fuel cells (SOFCs) and oxygen sensing devices. Doping the cubic phase
of zirconia (c-ZrO<sub>2</sub>) with yttria (Y<sub>2</sub>O<sub>3</sub>) is isoelectronic, as two Zr<sup>4+</sup> ions are replaced by two
Y<sup>3+</sup> ions, plus a charge compensating oxygen vacancy (O<sub>vac</sub>). Typical doping concentrations include 3, 8, 10, and 12
mol %. For these concentrations, and all below 40 mol %, no phase
with long-range order has been observed in either X-ray or neutron
diffraction experiments. The prediction of local defect structure
and the interaction between defects is therefore of great interest.
This has not been possible to date as the number of possible defect
topologies is very large and to perform reliable total energy calculations
for all of them would be prohibitively expensive. Previous theoretical
studies have only considered a selection of representative structures.
In this study, a comprehensive search for low-energy defect structures
using a combined classical modeling and density functional theory
approach is used to identify the low-energy isolated defect structures
at the dilute limit, 3.2 mol %. Through analysis of energetics computed
using the best available Born–Mayer–Huggins empirical
potential model, a point charge model, DFT, and a local strain energy
estimated in the harmonic approximation, the main chemical and physical
descriptors that correlate to the low-energy DFT structures are discussed.
It is found that the empirical potential model reproduces a general
trend of increasing DFT energetics across a series of locally strain
relaxed structures but is unreliable both in predicting some incorrect
low-energy structures and in finding some metastable structures to
be unstable. A better predictor of low-energy defect structures is
found to be the total electrostatic energy of a simple point charge
model calculated at the unrelaxed geometries of the defects. In addition,
the strain relaxation energy is estimated effectively in the harmonic
approximation to the imaginary phonon modes of undoped c-ZrO<sub>2</sub> but is found to be unimportant in determining the low-energy defect
structures. These results allow us to propose a set of easily computed
descriptors that can be used to identify the low-energy YSZ defect
structures, negating the combinatorial complexity and number of defect
structures that need to be considered