Modelling silver thin film growth on zinc oxide

Ag thin film growth on ZnO substrates has been investigated theoretically using multi-timescale simulation methods. The models are based on an atomistic approach where the interactions between atoms are treated classically using a mixture of fixed and variable charge potential energy functions. After some preliminary tests it was found that existing fixed charge potential functions were unreliable for surface growth simulations. This resulted in the development of a ReaxFF variable charge potential fitted to Ag/ZnO surface interactions. Ab initio models of simple crystal structures and surface configurations were used for potential fitting and testing. The dynamic interaction of the Ag atoms with the ZnO surface was first investigated using single point depositions, via molecular dynamics, whereby the Ag impacted various points on an irreducible symmetry zone of the ZnO surface at a range of energies. This enabled the determination of the relative numbers of atoms that could penetrate, reflect or bond to the surface as a function of incident energy. The results showed that at an energy of up to 10 eV, most atoms deposited adsorbed on top of the surface layer. The second part of the dynamic interaction involved a multi-timescale technique whereby molecular dynamics (MD) was used in the initial stages followed by an adaptive kinetic Monte Carlo (AKMC) approach to model the diffusion over the surface between impacts. An impact energy of 3 eV was chosen for this investigation. Ag was grown on various ZnO surfaces including perfect polar, O-deficient and surfaces with step edges. Initial growth suggests that Ag prefers to be spread out across a perfect surface until large clusters are forced to form. After further first layer growth, subsequent Ag atoms begin to deposit on the existing Ag clusters and are unlikely to join the first layer. Ag island formation (as mentioned within the literature) can then occur via this growth mechanism. O-deficient regions of ZnO surfaces result in unfavourable Ag adsorption sites and cause cluster formation to occur away from O-vacancies. In contrast, ZnO step edges attract deposited Ag atoms and result in the migration of surface Ag atoms to under-coordinated O atoms in the step edge. Various improvements have been made to the existing methodology in which transitions are determined. A new method for determining defects within a system, by considering the coordination number of atoms, is shown to increase the number of transitions found during single ended search methods such as the relaxation and translation (RAT) algorithm. A super-basin approach based on the mean rate method is also introduced as a method of accelerating a simulation when small energy barriers dominate. This method effectively combines states connected by small energy barriers into a single large basin and calculates the mean time to escape such basin. To accelerate growth simulations further and allow larger systems to be considered, a lattice based adaptive kinetic Monte Carlo (LatAKMC) method is developed. As off-lattice AKMC and MD results suggest Ag resides in highly symmetric adsorption sites and that low energy deposition events lead to no penetrating Ag atoms or surface deformation, the on-lattice based approach is used to grow Ag on larger perfect polar ZnO surfaces. Results from the LatAKMC approach agree with off-lattice AKMC findings and predict Ag island formation. Critical island sizes of Ag on ZnO are also approximated using a mean rate approach. Single Ag atoms are placed above an existing Ag cluster and all transition states are treated as belonging to a single large super-basin . Results indicate that small Ag clusters on the perfect ZnO surface grow in the surface plane until a critical island size of around 500 atoms is reached. Once a critical island size is reached, multiple Ag ad-atoms will deposit on the island before existing Ag atoms join the cluster layer and hence islands will grow upwards. A marked difference is seen for second layer critical island sizes; second layer Ag islands are predicted to be two orders of magnitude smaller (< 7 atoms). This analysis suggests that Ag on ZnO (000 ̄1) may exhibit Stranski-Krastanov (layer plus island) growth

Critical island size for Ag thin film growth on ZnO (000-1)

Island growth of Ag on ZnO is investigated with the development of a new technique to approximate critical island sizes. Ag is shown to attach in one of three highly symmetric sites on the ZnO surface or initial monolayers of grown Ag. Due to this, a lattice based adaptive kinetic Monte Carlo (LatAKMC) method is used to investigate initial growth phases. As island formation is commonly reported in the literature, the critical island sizes of Ag islands on a perfect polar ZnO surface and a first monolayer of grown Ag on the ZnO surface are considered. A mean rate approach is used to calculate the average time for an Ag ad-atom to drop off an island and this is then compared to deposition rates on the same island. Results suggest that Ag on ZnO (0 0 0 View the MathML source1¯) will exhibit Stranski–Krastanov (layer plus island) growth

Growth of silver on zinc oxide via lattice and off-lattice adaptive kinetic Monte Carlo

The growth of Ag on ZnO was modelled using a reactive force field potential and a combination of molecular dynamics and adaptive kinetic Monte Carlo (AKMC) simulations. An adaptive lattice-based AKMC model is described as a method of extending timescales and length scales that can be simulated. Reusing previously found transitions to reduce computational time is discussed for both the lattice and offlattice AKMC approaches. With these methods, growth of over 1 monolayer’s worth of Ag is simulated corresponding to a real deposition time of up to 0.1 s. The results show that the deposited silver aggregates on the surface through mainly single atom moves with few concerted motions. Initially silver adatoms do not agglomerate and the energy barriers for silver dimers to form is larger than for them to break apart. The first layer of silver grows as a series of connected regions rather than forming well-defined centro-symmetric islands

Development of an empirical interatomic potential for the Ag–Ti system

Two interatomic potential mixing rules for the Ti–Ag system were investigated based on the embedded-atom method (EAM) elemental potentials. First principles calculations were performed using SIESTA for various configurations of the Ti–Ag system to see which model best fitted the ab initio results. The results showed that the surface energies, especially that of Ti, were not well fitted by either model and the surface binding energies differed from the ab initio calculations. As a result, the modified embedded-atom method (MEAM) was investigated. In contrast to the other models, surface energies for pure Ti calculated by MEAM were in good agreement with the experimental data and the ab initio results. The MEAM mixing rule was used to investigate Ag ad-atoms on Ti and Ti ad-atoms on Ag. The results showed good agreement with SIESTA after parameter optimisation

Modelling thin film growth in the Ag-Ti system

Simulations of thin lm growth in the Ag-Ti system are presented using molecular dynamics combined with an adaptive kinetic Monte Carlo method (AKMC) with a modi ed embedded atom potential t to ab initio data for the surface energies. For the model, atoms are assumed to deposit normally with a kinetic energy of 1-3 eV, with a typical deposition rate of around 10 monolayers per second, similar to what might be expected in a sputter deposition process. For the growth of Ti on the Ag (100) and Ag (111) surfaces, the Ti adatoms prefer to exchange with the original surface layer atoms creating a mixed Ag/Ti surface. On a silver substrate, up to four mixed layers need to be formed before a pure Ti layer is obtained. Conversely, simulations of Ag depositing onto Ti (0001) showed that in the initial phase of growth, the Ag adatoms prefer to be separated before a complete rst layer of Ag was obtained in a close-packed structure. The implementation of a super-basin method within AKMC allowed the simulation of 0.4s of Ti growth on the Ag substrates, with up to 3 new layers added

Thermal dynamics of silver clusters grown on rippled silica surfaces

Silver nanoparticles have been deposited on silicon rippled patterned templates at an angle of incidence of 70° to the surface normal. The templates are produced by oblique incidence argon ion bombardment and as the fluence increases, the periods and heights of the structures increase. Structures with periods of 20 nm, 35 nm and 45 nm have been produced. Moderate temperature vacuum annealing shows the phenomenon of cluster coalescence following the contour of the more exposed faces of the ripple for the case of 35 nm and 45 nm but not at 20 nm where the silver aggregates into larger randomly distributed clusters. In order to understand this effect, the morphological changes of silver nanoparticles deposited on an asymmetric rippled silica surface are investigated through the use of molecular dynamics simulations for different deposition angles of incidence between 0° and 70° and annealing temperatures between 500 K and 900 K. Near to normal incidence, clusters are observed to migrate over the entire surface but for deposition at 70°, a similar patterning is observed as in the experiment. The random distribution of clusters for the periodicity of 20 nm is linked to the geometry of the silica surface which has a lower ripple height than the longer wavelength structures. Calculations carried out on a surface with such a lower ripple height also demonstrate a similar effect

Reaction pathways in atomistic models of thin film growth

The atomistic processes that form the basis of thin film growth often involve complex multi-atom movements of atoms or groups of atoms on or close to the surface of a substrate. These transitions and their pathways are often difficult to predict in advance. By using an adaptive kinetic Monte Carlo (AKMC) approach many complex mechanisms can be identified so that the growth processes can be understood and ultimately controlled. Here the AKMC technique is briefly described along with some special adaptions that can speed up the simulations when, for example, the transition barriers are small. Examples are given of such complex processes that occur in different material systems especially for the growth of metals and metallic oxides

Development of a ReaxFF potential for Ag/Zn/O and application to Ag deposition on ZnO

A new empirical potential has been derived to model an Ag–Zn–O system. Additional parameters have been included into the reactive force field (ReaxFF) parameter set established for ZnO to describe the interaction between Ag and ZnO for use in molecular dynamics (MD) simulations. The reactive force field parameters have been fitted to density functional theory (DFT) calculations performed on both bulk crystal and surface structures. ReaxFF accurately reproduces the equations of state determined for silver, silver zinc alloy and silver oxide crystals via DFT. It also compares well to DFT binding energies and works of separation for Ag on a ZnO surface. The potential was then used to model single point Ag deposition on polar (000View the MathML source1¯) and non-polar (10View the MathML source1¯0) orientations of a ZnO wurtzite substrate, at different energies. Simulation results then predict that maximum Ag adsorption on a ZnO surface requires deposition energies of ≤ 10 eV

ZrH Lammps Input and Structures

Lammps input and structure files for Zr/H systems relating to modelling the effect of hydrogen on crack growth in zirconium via molecular dynamic

ZrH strain and crack velocity plots

Plots of Zr/H stress-strain relation and crack velocity of Zr/H structures under strain