Junlei Zhao

Formation Mechanism of Fe Nanocubes by Magnetron Sputtering Inert Gas Condensation

In this work, we study the formation mechanisms of iron nanoparticles (Fe NPs) grown by magnetron sputtering inert gas condensation and emphasize the decisive kinetics effects that give rise specifically to cubic morphologies. Our experimental results, as well as computer simulations carried out by two different methods, indicate that the cubic shape of Fe NPs is explained by basic differences in the kinetic growth modes of {100} and {110} surfaces rather than surface formation energetics. Both our experimental and theoretical investigations show that the final shape is defined by the combination of the condensation temperature and the rate of atomic deposition onto the growing nanocluster. We, thus, construct a comprehensive deposition rate-temperature diagram of Fe NP shapes and develop an analytical model that predicts the temporal evolution of these properties. Combining the shape diagram and the analytical model, morphological control of Fe NPs during formation is feasible; as such, our method proposes a roadmap for experimentalists to engineer NPs of desired shapes for targeted applications.

Migration barriers for surface diffusion on a rigid lattice: Challenges and solutions

Abstract Atomistic rigid lattice Kinetic Monte Carlo is an efficient method for simulating nano-objects and surfaces at timescales much longer than those accessible by molecular dynamics. A laborious part of constructing any Kinetic Monte Carlo model is, however, to calculate all migration barriers that are needed to give the probabilities for any atom jump event to occur in the simulations. One of the common methods of barrier calculations is Nudged Elastic Band. The number of barriers needed to fully describe simulated systems is typically between hundreds of thousands and millions. Calculations of such a large number of barriers of various processes is far from trivial. In this paper, we will discuss the challenges arising during barriers calculations on a surface and present a systematic and reliable tethering force approach to construct a rigid lattice barrier parameterization of face-centred and body-centred cubic metal lattices. We have produced several different barrier sets for Cu and for Fe that can be used for KMC simulations of processes on arbitrarily rough surfaces. The sets are published as Data in Brief articles and available for the use.

Data sets of migration barriers for atomistic Kinetic Monte Carlo simulations of Cu self-diffusion via first nearest neighbour atomic jumps

Atomistic rigid lattice Kinetic Monte Carlo (KMC) is an efficient method for simulating nano-objects and surfaces at timescales much longer than those accessible by molecular dynamics. A laborious and non-trivial part of constructing any KMC model is, however, to calculate all migration barriers that are needed to give the probabilities for any atom jump event to occur in the simulations. We have calculated three data sets of migration barriers for Cu self-diffusion with two different methods. The data sets were specifically calculated for rigid lattice KMC simulations of copper self-diffusion on arbitrarily rough surfaces, but can be used for KMC simulations of bulk diffusion as well.

Data sets of migration barriers for atomistic Kinetic Monte Carlo simulations of Fe self-diffusion

Atomistic rigid lattice Kinetic Monte Carlo (KMC) is an efficient method for simulating nano-objects and surfaces at timescales much longer than those accessible by molecular dynamics. A laborious and non-trivial part of constructing any KMC model is, however, to calculate all migration barriers that are needed to give the probabilities for any atom jump event to occur in the simulations. We calculated three data sets of migration barriers for Fe self-diffusion: barriers of first nearest neighbour jumps, second nearest neighbours hop-on jumps on the Fe {100} surface and a set of barriers of the diagonal exchange processes for various cases of the local atomic environments within the 2nn coordination shell.