Kinetic Monte Carlo simulations of Pd deposition and island growth on MgO(100).
Abstract
The research described in this product was performed in part in the Environmental Molecular Sciences Laboratory, a national scientific user facility sponsored by the Department of Energy's Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory. The deposition and ripening of Pd atoms on the MgO(10 0) surface are modeled using kinetic Monte Carlo simulations. The density of Pd islands is obtained by simulating the deposition of 0.1 ML in 3 min. Two sets of kinetic parameters are tested and compared with experiment over a 200–800 K temperature range. One model is based upon parameters obtained by fitting rate equations to experimental data and assuming the Pd monomer is the only diffusing species. The other is based upon transition rates obtained from density functional theory calculations which show that small Pd clusters are also mobile. In both models, oxygen vacancy defects on the MgO surface provide strong traps for Pd monomers and serve as nucleation sites for islands. Kinetic Monte Carlo simulations show that both models reproduce the experimentally observed island density versus temperature, despite large differences in the energetics and different diffusion mechanisms. The low temperature Pd island formation at defects is attributed to fastmore »
- Authors:
- Publication Date:
- Research Org.:
- Pacific Northwest National Laboratory (PNNL), Richland, WA (US), Environmental Molecular Sciences Laboratory (EMSL)
- Sponsoring Org.:
- USDOE
- OSTI Identifier:
- 1012315
- DOE Contract Number:
- AC05-76RL01830
- Resource Type:
- Journal Article
- Resource Relation:
- Journal Name: Surface Science, 601(14):3133-3142; Journal Volume: 601; Journal Issue: 14
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 36 MATERIALS SCIENCE; PALLADIUM; DEPOSITION; SUBSTRATES; MAGNESIUM OXIDES; MONTE CARLO METHOD; KINETIC EQUATIONS; DENSITY FUNCTIONAL METHOD; BINDING ENERGY; VACANCIES; Environmental Molecular Sciences Laboratory
Citation Formats
Xu, Lijun, Campbell, Charles T., Jonsson, Hannes, and Henkelman, Graeme A. Kinetic Monte Carlo simulations of Pd deposition and island growth on MgO(100).. United States: N. p., 2007.
Web. doi:10.1016/j.susc.2007.05.027.
Xu, Lijun, Campbell, Charles T., Jonsson, Hannes, & Henkelman, Graeme A. Kinetic Monte Carlo simulations of Pd deposition and island growth on MgO(100).. United States. doi:10.1016/j.susc.2007.05.027.
Xu, Lijun, Campbell, Charles T., Jonsson, Hannes, and Henkelman, Graeme A. Thu .
"Kinetic Monte Carlo simulations of Pd deposition and island growth on MgO(100).". United States.
doi:10.1016/j.susc.2007.05.027.
@article{osti_1012315,
title = {Kinetic Monte Carlo simulations of Pd deposition and island growth on MgO(100).},
author = {Xu, Lijun and Campbell, Charles T. and Jonsson, Hannes and Henkelman, Graeme A.},
abstractNote = {The research described in this product was performed in part in the Environmental Molecular Sciences Laboratory, a national scientific user facility sponsored by the Department of Energy's Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory. The deposition and ripening of Pd atoms on the MgO(10 0) surface are modeled using kinetic Monte Carlo simulations. The density of Pd islands is obtained by simulating the deposition of 0.1 ML in 3 min. Two sets of kinetic parameters are tested and compared with experiment over a 200–800 K temperature range. One model is based upon parameters obtained by fitting rate equations to experimental data and assuming the Pd monomer is the only diffusing species. The other is based upon transition rates obtained from density functional theory calculations which show that small Pd clusters are also mobile. In both models, oxygen vacancy defects on the MgO surface provide strong traps for Pd monomers and serve as nucleation sites for islands. Kinetic Monte Carlo simulations show that both models reproduce the experimentally observed island density versus temperature, despite large differences in the energetics and different diffusion mechanisms. The low temperature Pd island formation at defects is attributed to fast monomer diffusion to defects in the rate-equation-based model, whereas in the DFT-based model, small clusters form already on terraces and diffuse to defects. In the DFT-based model, the strong dimer and trimer binding energies at charged oxygen vacancy defects prevent island ripening below the experimentally observed onset temperature of 600 K.},
doi = {10.1016/j.susc.2007.05.027},
journal = {Surface Science, 601(14):3133-3142},
number = 14,
volume = 601,
place = {United States},
year = {Thu May 24 00:00:00 EDT 2007},
month = {Thu May 24 00:00:00 EDT 2007}
}
-
The deposition and ripening of Pd atoms on the MgO(10 0) surface are modeled using kinetic Monte Carlo simulations. The density of Pd islands is obtained by simulating the deposition of 0.1 ML in 3 min. Two sets of kinetic parameters are tested and compared with experiment over a 200–800 K temperature range. One model is based upon parameters obtained by fitting rate equations to experimental data and assuming the Pd monomer is the only diffusing species. The other is based upon transition rates obtained from density functional theory calculations which show that small Pd clusters are also mobile.
-
Hydrogenation of Acetylene-Ethylene Mixtures over Pd and Pd-Ag Alloys: First-Principles Based Kinetic Monte Carlo Simulations
The kinetics for the selective hydrogenation of acetylene-ethylene mixtures over model Pd(111) and bimetallic Pd-Ag alloy surfaces were examined using first principles based kinetic Monte Carlo (KMC) simulations to elucidate the effects of alloying as well as process conditions (temperature and hydrogen partial pressure). The mechanisms that control the selective and unselective routes which included hydrogenation, dehydrogenation and C-C bond breaking pathways were analyzed using first-principle density functional theory (DFT) calculations. The results were used to construct an intrinsic kinetic database that was used in a variable time step kinetic Monte Carlo simulation to follow the kinetics and the molecularmore » -
Non-Arrhenius temperature dependence of the island density of one-dimensional Al chains on Si(100): A kinetic Monte Carlo study
Classical nucleation theory predicts that the evolution of mean island density with temperature during growth in one-dimensional systems obeys the Arrhenius relation. In this study, kinetic Monte Carlo simulations of a suitable atomistic lattice-gas model were performed to investigate the experimentally observed non-Arrhenius scaling behavior of island density in the case of one-dimensional Al islands grown on Si(100). Previously, it was proposed that adatom desorption resulted in a transition temperature signaling the departure from classical predictions. Here, the authors demonstrate that desorption above the transition temperature is not possible. Instead, the authors posit that the existence of a transition temperaturemore » -
Monte Carlo simulation of the nucleation and growth of binary-alloy particles of Au, Ag, and Pd on NaCl(100) substrates
The initial phase of thin-film formation from two components that are simultaneously deposited on a substrate was simulated by Monte Carlo calculations of adatom processes, e.g., surface diffusion, reevaporation, nucleation, and cluster growth by collisions and capture. Experiments had revealed that the condensation of Pd is nearly complete on NaCl at 300 {degree}C, whereas Au and Ag exhibit very low initial condensation coefficients. Therefore, the composition of the nuclei differs strongly from that of the vapor beam. Monte Carlo calculations of the composition at the early stage of Au-Pd and of Ag-Pd depositions were fitted to experimental data from earliermore » -
Dissolution and Growth of (10(1) over-bar4) Calcite in Flowing Water: Estimation of Back Reaction Rates via Kinetic Monte Carlo Simulations
Although calcite is an important mineral for many processes, there have been relatively few simulations of it's growth and dissolution behavior. Such simulations are complicated by the multitude of defect types and by the asymmetry of the crystal. The present work combined a kinetic Monte Carlo (KMC) technique with the Kossel crystal (100) simple cubic concept and the Blasius boundary layer model to simulate the simultaneous growth and dissolution of the (1014)calcite cleavage surface in flowing water. The objective was to determine the activation energies of the back reaction (growth) from those of the forward reaction (dissolution) by obtaining agreementmore »