Research Article  Open Access
Karin Kleiner, Aleix ComasVives, Maryam Naderian, Jonathan E. Mueller, Donato Fantauzzi, Mostafa Mesgar, John A. Keith, Josef Anton, Timo Jacob, "Multiscale Modeling of AuIsland Ripening on Au(100)", Advances in Physical Chemistry, vol. 2011, Article ID 252591, 11 pages, 2011. https://doi.org/10.1155/2011/252591
Multiscale Modeling of AuIsland Ripening on Au(100)
Abstract
We describe a multiscale modeling hierarchy for the particular case of Auisland ripening on Au(100). Starting at the microscopic scale, density functional theory was used to investigate a limited number of selfdiffusion processes on perfect and imperfect Au(100) surfaces. The obtained structural and energetic information served as basis for optimizing a reactive forcefield (here ReaxFF), which afterwards was used to address the mesoscopic scale. Reactive force field simulations were performed to investigate more diffusion possibilities at a lower computational cost but with similar accuracy. Finally, we reached the macroscale by means of kinetic Monte Carlo (kMC) simulations. The reaction rates for the reaction process database used in the kMC simulations were generated using the reactive force field. Using this strategy, we simulated nucleation, aggregation, and fluctuation processes for monoatomic high islands on Au(100) and modeled their equilibrium shape structures. Finally, by calculating the step line tension at different temperatures, we were able to make a direct comparison with available experimental data.
1. Introduction
Surface diffusion plays a key role in dynamical processes occurring on metallic surfaces which involve physical or chemical displacements of atoms or molecules on a particular surface [1], such as adsorptiondesorption, crystal growth, coarsening, chemical reactions, wetting, spreading, or capillarity. Hence, a profound understanding of this phenomenon has implications in a great variety of fields such as electrochemistry, microelectronics, tribology, and corrosion protection [2]. Two main pathways can be distinguished for the diffusion of adatoms across a surface. The first one corresponds to the hopping of an adatom between different adsorption sites, whereas in the socalled exchange mechanism, adatoms replace atoms of the surface layer (or fall into vacancies) while surface atoms are simultaneously pushed onto adsorption sites [3]. The latter pathway was proposed for Al(001) [4] and Au(100) [5] based on theoretical studies (Müller and Ibach stated that the hopping mechanism could be also competitive on Au(100)) [6]. While for planar terraces these processes are rather well understood, [7–14] much less is known about diffusion processes on imperfect surfaces with lower coordinated adsorption sites (e.g., kink and stepsites or vacancies). However, these types of sites often determine morphological changes of the surface structure, for instance, in Ostwald ripening, metal deposition and dissolution, or island and step fluctuation processes. It should be noted that clean Au(100) shows a quasihexagonal reconstruction in the first surface layer [15], however, this reconstruction can be lifted in an electrochemical cell by applying a sufficiently positive potential. Under these conditions, chemisorbed species (e.g., anions) and electric fields are present at the metal interface. Although these parameters are not taken into account in the present work; we reported the influence of these parameters in previous experimental and theoretical studies [16–18]. Moreover, the lifting and the restoration of the reconstruction are slow processes, so that both types of surfaces can be assumed to be present over a sizable range of potentials.
In order to address diffusion phenomena at surfaces theoretically, different approaches are used, which range from experimentally based, semiempirical molecular dynamics up to first principles studies. In the present work, we describe a first principlesbased approach for investigating surface dynamics on Au(100). While improvements in experimental techniques have led to considerable progress in the investigation of surface diffusion, there is still a lack on the experimental data for the selfdiffusion on Au(100).
The basis for the present study is a recent contribution from our group where the most reasonable processes for selfdiffusion on Au(100) were studied using quantum mechanics (density functional theory, DFT). This study specifically focused on terraces and imperfect kink and stepconfigurations [19]. Since the rate of adatom migration depends mainly on its immediate environment, we concentrated on the nearest neighbor interactions only. Nevertheless, with this restriction, already more than a thousand surface configurations have to be considered [20]. Given the high computational expense of DFT calculations, with the present computing resources, one is rather limited in the number of systems that can be studied from first principles. However, we can use the binding energies and selfdiffusion processes previously studied by ab initio simulations to construct and optimize a reactive force field (here ReaxFF) that is able to reproduce this behavior. After the mandatory tests related to accuracy and transferability of the force field, ReaxFF has the advantage of being nearly as accurate as QM but with much lower computational demands [21]. Using ReaxFF to investigate all possible diffusion processes (considering only firstnearest neighbors), the obtained process rates were used with kinetic Monte Carlo simulations (kMC) for largescale simulations that allowed us to analyze island nucleation, aggregation, and fluctuation as well as the equilibrium shapes of adatom islands on Au(100). Figure 1 shows the overall computational procedure, where we started at the ab initio level of theory to obtain a reactive molecular dynamics forcefields, which was then used to provide the necessary process rates for macroscopic kinetic Monte Carlo simulations. Similar hierarchical simulation schemes have been employed by the groups of Goddard, Neurock, and van Santen, for instance [22–25].
2. First Principles Calculations
We used density functional theory to investigate different selfdiffusion mechanisms in the presence or absence of a variety of surface defects on Au(100), for example, kink and stepedges. Usually, DFT calculations in the gradient generalized approximations (such as the PBE density functional used in this work) are able to achieve a chemical accuracy of ≤0.1 eV for reactions [13]. Since the details of these studies are described in [19], here we only give a brief summary. From these first principle studies, we were able to obtain binding energies of the stable intermediates and transition states and to extract rate constants for the diffusion events. An overview of the different pathways that have been studied is given in Figure 2.
(a) Terrace
(b) Stepedge
(c) Kink
(d) Innercorner
(e) Outercorner
(f) Stepvacancy
(g) Dimer
On terraced Au(100) surfaces, we found that the lowest energy mechanism was the exchange diffusion pathway from position A to C , while hopping diffusion from A to B over a bridge site was less favorable . The last evaluated reaction pathway at terraces was the atom hopping over the top position from A to C, being the least favorable with a rather high energy barrier of 1.39 eV.
The presence of stepedges has a tremendous influence on the Auadatom diffusion, since the energy barrier is even lower for the diffusion along the stepedge (0.38 eV). As expected, diffusion perpendicular to or from the step edge has a much higher barrier of 0.84 eV (detachment) and 0.66 eV (attachment). Comparing the barriers shows that rearrangement of adatoms along the step would be faster than island growth or decay, which would lead to islands that quickly assume their equilibrium shapes before changing their size.
Afterwards, systems involving corners at stepedges were analyzed (see Figure 2). For the innercorner system, we found that diffusion along long edges is preferred over migration along short edges. This suggests that Au atoms adsorbed at terraces first reach kinksites or innercorner positions at islands by exchange diffusion, followed by a hopping diffusion to the stepedges where they rapidly diffuse along the stepedges. As far as the outercorner systems are concerned, adatoms, which are located directly next to an island corner, have essentially the same binding energy as an adatom at a stepedge. In addition, we found that detachment from the outercorner was preferred over detachment away from the stepedge.
Finally, we calculated diffusion pathways for forming Audimers, which can be considered the first step of island formation, and for extraction of an adatom out of a stepedge forming a stepvacancy. While the latter process has a relatively high energy barrier of 1.07 eV, dimer diffusion is comparable to the case of a single atom.
In addition to reaction barriers, we calculated vibrational frequencies factors for each reaction. This was accomplished by fitting a frequency to the energy profile of each reactant state to within the harmonic approximation. Using these frequencies, we calculated the rates for the different diffusion events at variable temperatures (see Table 1) by means of the transitionstate theory, as given by equation (1): where is a prefactor, is the barrier for diffusion (activation), is the Boltzmann constant, and is the temperature (all of our mechanisms are firstorder kinetic reactions, and so these terms are simply the vibrational frequency of the selfdiffusion process in units of Hz). One of the main limitations of this approach is that it does not account for anharmonicites and thermal expansion, which cause deviations at higher temperatures. Moreover, it also ignores memory effects arising from the surface excitations, which change the effective Arrhenius barrier. Nevertheless, the Arrhenius law’s validity within the limit of provides reasonable temperature ranges, for which our results are valid [26].

The obtained diffusion rate constants (Table 1) clearly illustrate the importance of the surface structure on cluster formation and island growth. The fastest rate corresponds to the filling of an innercorner by a diffusing atom (IV2 pathway), although the same system also yields the slowest rate with pathway IV3 (not considering the very unfavorable AC hopping diffusion at terraces). Concerning the terraces, the diffusion over the top position is practically negligible while the exchange diffusion rate is around 300 times larger than the bridge diffusion rate. Our calculated activation energies for the AB hopping process and for the AC exchange process on the terrace are 0.70 and 0.55 eV, respectively. These values are in reasonable agreement with the approximated experimental activation barrier of 0.5 eV for the adatom hopping on Au(100) [27]. In addition, they also agree well with previous theoretical studies. For example, Müller and Ibach reported energy barriers of 0.64 eV and 0.60 eV for the hopping and the exchange processes on an Au(100) [6].
Although ab initio studies only allow the investigation of a limited set of diffusion possibilities, the studies presented so far already provide useful insights into diffusion events on Au(100) with important implications for Ostwald ripening and island formation. However, the next step is to account for more possibilities, as described subsequently in the list of possible diffusion events section, by means of larger scale reactive forcefield calculations. Therefore, in the next section we will summarize the generation of the corresponding Auforcefield, which was then used to evaluate the rates for many more diffusion scenarios.
3. Reactive Forcefield (ReaxFF) Simulations
The next step in the multiscale approach presented in this work corresponds to the generation and application of a reactive forcefield within the ReaxFF framework [28–30]. ReaxFF is a reactive molecular dynamics method that uses a bondorderdependent potential energy formulation, similar to the Tersoff [31] or Brenner [32] potentials. These potentials are all based on Pauling’s idea of mapping bond distances onto bond orders to enable the determination of different quantum chemical states of a molecular structure [33].
The dependence of the energy contributions on the bond order also means that such energy terms implicitly contain multibody contributions. The Au forcefield employed in the present study contains the following three energy terms (2): where is the energy corresponding to interatomic bonds, is a penalty energy that corrects atomic overcoordinations, and accounts for van der Waals interactions and interatomic repulsions when atoms are too close to each other. The full expressions for these terms as well as the included parameters can be found in more detail in a previous publication from our group [21].
The formulation of the system energy is more sophisticated than those of nonreactive potentials, for example, EAM [34–37], MEAM [38], UFF [39, 40], CHARMM [41], OPLS [42], or AMBER [43, 44], and therefore computationally more demanding. Reactive forcefields represent a useful tool for overcoming limitations of these empirical potentials as it allows for the description of chemical reactions, that is, bond formation and dissociation, with almost QM accuracy. In contrast to most other forcefields, the ReaxFF potential used in the present work is completely trained against ab initio DFT calculations.
The goal of this training process is to optimize the parameters of the forcefield such that the DFT energies and structures are reproduced as accurately as possible. In this specific case, first principles DFTPBE calculations on the equations of state for several gold bulk phases, as well as binding energies and selfdiffusion processes on Au(100) surfaces as described in the previous section, were used in the training procedure. The validity and the transferability of the obtained interaction potential was successfully tested for additional diffusion processes, surface free energies of Au low index surfaces, and cohesive energies of molecular Au clusters. Details concerning these tests as well as a first application to study the morphology of larger Au nanoparticles can be found in [21]. The purpose was to obtain a large database of diffusion processes on Au(100) under a wide range of defects that one might expect to exist on real surfaces, in order to obtain the diffusion barriers which are needed to run kMC simulations. As already mentioned before, the simulation of all the diffusion processes that could occur is a task that cannot be achieved in a reasonable time using DFT or other ab initio methods.
Figure 3 shows exemplarily the performance of the ReaxFF Auforcefield after the optimization procedure. Besides simple terrace (bridge) diffusion, diffusion along step edges, around kink sites, and diffusion forming molecular dimers are compared. The comparison shows that the forcefield is indeed a useful tool for calculating all the remaining diffusion processes required to evaluate the database of diffusion rates required for the kMC approach, and thus for bridging the gap between QM and the macroscopic regime.
(a) Terrace diffusion
(b) Step diffusion
(c) Kink diffusion
(d) Dimer diffusion
4. Kinetic Monte Carlo
Kinetic Monte Carlo (kMC) is a wellknown technique for largescale simulations in various areas of research such as chemical physics or materials design [45–47]. In order to perform such kMC simulations on the dynamic behavior of twodimensional gold islands on Au(100), we used the previously described ReaxFF forcefield to generate a list of independent (Markov chain) diffusion rates for all possible diffusion scenarios (taking the firstnearest neighbors into account) [48].
Along the kMC simulations, the total rate (cumulative function) of all diffusion events of a particular system is calculated by means of (3):
In this function, all possible elementary diffusion processes are taken into account by their individual rates . Then, by using a random number, , between zero and one, we assign it to a rate out of all the processes described in the cumulative function (): where then specifies the particular process which will be executed. After each diffusion event, the simulation time is propagated by
However, for a realistic value of , we have to allow fluctuations, therefore, a second random number, , is chosen between zero and one and afterwards used in (5) in the following way:
After executing the particular diffusion process and propagating the time the iterative simulation proceeds by calculating the cumulative function () of the new structure, and so forth. As the kMC approach is based on random numbers, reasonable statistics are required to allow for reliable qualitative and quantitative conclusions.
5. List of Possible Diffusion Events
Given the structure of the Au(100) surface, the number of possible diffusion events, (considering only nextnearest neighbors and neglecting symmetry considerations) is 2^{10} = 1024, which is the number of reactions required for the rate catalog. The profile of the considered diffusion pathways is depicted in Figure 4. While the atom of interest diffuses from position 1 to 2, the nearest neighbor (NN) sites (labeled from 3 to 12) are either empty or occupied.
(a)
(b)
Although the classification described in Figure 4 only considers the nextnearest neighbors, when calculating the diffusion events with ReaxFF and afterwards the rates to be used in the following kMC simulations, a slightly modified system has been employed. For the ReaxFF studies, each system was characterized by a fourlayer slab for the underlying Au(100) surface with a lateral extension of 10 × 10 atoms. When occupying the surrounding sites, not only the positions specified by the labels 3–12 have been filled, but always the entire outgoing rows of atoms. Thus, for instance, the diffusion process indicated in Figure 4(a) has been studied with the model shown in Figure 4(b). On the basis of this structure, we used ReaxFF to map the energy profile when moving the diffusion atom (red) successively from position 1 to 2. The energy profile provided the activation energy and vibrational frequencies, which together with transition state theory (TST) provided all diffusion rates to be included in the event database for the kMC studies.
6. Input Structures
Kinetic Monte Carlo was utilized to study different systems, in which the Au(100) substrate was modeled by a square lattice of 100 × 100 positions. Periodic boundary conditions were imposed in x and ydirections (parallel to the substrate) and all atoms were assumed to be uncharged. In this work, we analyzed four different initial input configurations, their evolution are being followed at 850 K. Figure 5 shows (a) a random distribution of adatoms, (b) a twoisland structure, (c) a spherical island, and (d) a square island.
(a)
(b)
(c)
(d)
It should be noted that although during the kMC simulations the systems were always free to evolve, for simplicity, the exchange mechanism at terraces was not included in the rate catalog.
7. Nucleation and Aggregation
Subsequently, we followed the evolution of the initial input structures (see Figure 5) at 850 K. The first structure (a) starts with a randomly distributed dispersion of adatoms on the Au(100) surface. In Figure 6, we show successive snapshots of the kMC simulation of this structure. Along the structural evolution, particle nucleation is observed, as, after several hundreds of kMC steps, the randomly dispersed adatoms begin forming nucleation seeds. These nucleation seeds afterwards form small islands that successively grow during the ~4000 ns simulation period. Thus, the process indicated in Figure 6 indeed shows the first stages of island formation. As expected, we also found that the nucleation speed increases with temperature and adatom coverage. This can be explained by the fact that with randomly dispersed adatoms at the beginning, the highest rate corresponds to the free hopping of adatoms on the surface. After several simulation steps, adatoms meet each other forming slowly growing islands. This makes the stepedge diffusion possible and avoids the diffusion away from the step (the latter process has a higher energy barrier), therefore, continuously favoring the growth of the islands.
(a)
(b)
(c)
(d)
With the second system (Figure 5(b)), we followed the evolution of a pair of disparately sized islands at 850 K (see Figure 7). After several kMC steps, the smaller island starts dissolving, while the detached Au atoms quickly migrate to the larger island where they are assimilated. Hence, this simulation provides insights on how island growth proceeds.
(a)
(b)
(c)
(d)
8. Island Shape (Equilibrium and Fluctuation)
Stepedge fluctuation is the dominating process in determining the dynamic shape of islands on the Au(100) surface. Several diffusion processes contribute to stepedge fluctuation: diffusion along the stepedge, along a kinksite, or along an outercorner. As explained in the DFT section, the rate for stepedgediffusion is the highest among these processes, followed by the diffusion along the outercorner and by the diffusion along kink sites.
Figure 8 shows the evolution of an initially spherical island (Figure 5(c)) at 850 K. The snapshots show that after several iterations, the roughness of the island border increases. Qualitatively, this is in good agreement with experimental results obtained by Giesen et al. by means of STM measurements in the case of twodimensional Cu islands on Cu(100) [14, 49].
(a)
(b)
(c)
(d)
Interestingly, starting with a squareshaped island and maintaining a temperature of 850 K (Figure 5(d)) showed no significant deviations from the behavior of the spherical island, and hence the corresponding snapshots are not shown. Every migration step of an adatom out of the square island produces a vacancy. This mechanism has the lowest reaction rate among all the considered ones, and therefore this occurs very rarely. This can be explained by the degrees of freedom. As noted already by Ibach [50] stepatoms can only move in one direction and therefore the binding order on the (100) surface is three. This means that three bonds need to be broken along this process, which results in an increased energy barrier for this diffusion process.
Finally, we analyzed the equilibrium shape of the islands for a fixed number of gold adatoms but at different temperatures. To obtain the equilibrium shape of islands at a given temperature, one needs to average over many individual samples [51]. Therefore, we analyzed the island shapes during the ripening process starting with an initially random adatom distribution Figure 5(a). Islands already begin to show their equilibrium shape during growth, when the effects of island aggregation are averaged out [49, 52, 53]. We studied the temperature dependence of equilibrium fluctuations of monatomic high islands on Au(100) by analyzing the obtained island shapes at three different temperatures, 500, 750, and 1000 K. As shown in Figure 9, the equilibrium island shape becomes rougher with increasing temperature, which was also suggested by Giesen et al. [49]. This roughening was found in several STMstudies on Pb and Cu, Ag, Pt, [49, 50, 54], and also for the case of Au(100) [55], The origin of this temperature dependence of the island shape is connected with the step line tension, which becomes smaller as the temperature increases, causing the island borders to increase their roughness. The step line tension of a monoatomic step on a solid surface is defined as the work per length required for generating the step. For uncharged surfaces, this work is equal to the Helmholtz step free energy. Bombis and Ibach reported a step line tension of 170 meV per atom length based on STM images in vacuum. This contrasted with the comparative smaller value of –70 meV per atom length obtained by Giesen et al. for the unreconstructed Au(100) in an electrolyte solution (sulfuric acid) at higher potentials [56–58]. Also, previous theoretical calculations have estimated step energies on Au(100) between 65 [59] and 180 meV [60] per step atom. Hence, in order to obtain a direct comparison between the available experimental data and our calculated values, we calculated the step line tension for different temperatures ranging from 200–800 K. We calculated the step line tension according to the same procedure used by Bombis and Ibach [55] and described by Schlößer et al. [61] and hence do not describe it in detail here. In this approach, the expectation value is calculated with (7): where is the middle radius of a circular island, is the number of islands used in calculating the average value (–80), is the number of angles for which is calculated (), and is the radius at a given moment. is proportional to 1/, where depends on the direction of the step relative to the surface. In the particular case of gold clusters, is obtained by (8):
From (8), we can calculate at different temperatures. These results are summarized in Table 2. The corresponding values can also be expressed as a function of the gold atom radius since is proportional to the island radius .

The relation between the temperature and follows a quadratic dependence. Based on this quadratic dependence and the data reported in Table 2, we obtained a value of 439 meV/nm at 353 K for the step line tension . As function of the gold atom radius , we obtained a value of 127 meV. At the same temperature (353 K), the value obtained by Bombis and Ibach of 170 meV can be compared to our value [55]. The difference between both values might be due to the rate catalog (exchange processes with the surface have so far not been considered) or the difference in the quotient of the step line tensions along perpendicular steps and . This quotient provides information on how far away the islands are from reaching equilibrium conditions (at equilibrium ). The value obtained by Bombis and Ibach for this quotient [55] was 0.96 compared with our value of 0.99 [62].
Ongoing work in our group seeks to include the exchange diffusion process and other processes which could influence the kMC description of the Ostwald ripening on Au(100) in future models.
9. Conclusions
In the present work, we have described a multiscale approach for studying the model case of Auisland ripening on Au(100). Starting with the atomistic level, we performed ab initio DFT calculations on a set of selected selfdiffusion processes on perfect and imperfect Au(100) surfaces, obtaining both binding and activation energies. Afterwards, this data was used to generate a reactive ReaxFF force field for Au–Au interactions, which provided the basis for investigating many additional diffusion scenarios. Indeed, using ReaxFF, we were able to consider all possible processes involving nearest neighbors. The thus generated database of diffusion rates was finally used to study the system on the macroscopic level by means of kMC simulations.
Using kMC, we evaluated different initial surface configurations, by which we were able to investigate adatom nucleation as well as island aggregation and fluctuation processes on Au(100). Finally, by averaging over island shapes during the simulation runs at different temperatures, we studied the equilibrium shape of monoatomic high islands at different temperatures (here 500, 750, and 1000 K). In good agreement with experimental results for gold islands on Au(100), we found that increasing the temperature leads to island roughening. We explain this temperature dependence of island roughing in terms of the temperaturedependent behavior of the step line tension. Our step line tensions agree fairly well with the reported experimental values and validate the present multiscale approach, which should be easily adaptable to the study of other systems.
Acknowledgments
The authors gratefully acknowledge financial support from the Deutsche Forschungsgemeinschaft (DFG) within the Emmy Noether Program, and the EU network ELCAT (Proposal no. 2149362, 2008–2012).
References
 G. Rosenfeld, K. Morgenstern, and G. Comsa, “Diffusion and stability of large clusters on crystal surfaces,” in Surface Diffusion: Atomistic and Collective Processes, M. C. Tringides, Ed., pp. 361–377, Plenum Press, New York, NY, USA, 1997. View at: Google Scholar
 A. G. Naumovets and Z. Zhang, “Fidgety particles on surfaces: how do they jump, walk, group, and settle in virgin areas?” Surface Science, vol. 500, no. 1–3, pp. 414–436, 2002. View at: Publisher Site  Google Scholar
 D. W. Bassett and P. R. Webber, “Diffusion of single adatoms of platinum, iridium and gold on platinum surfaces,” Surface Science, vol. 70, no. 1, pp. 520–531, 1978. View at: Google Scholar
 P. J. Feibelman, “Diffusion path for an Al adatom on Al(001),” Physical Review Letters, vol. 65, no. 6, pp. 729–732, 1990. View at: Publisher Site  Google Scholar
 B. D. Yu and M. Scheffler, “Physical origin of exchange diffusion on fcc(100) metal surfaces,” Physical Review B, vol. 56, no. 24, pp. R15569–R15572, 1997. View at: Google Scholar
 J. E. Müller and H. Ibach, “Migration of point defects at charged Cu, Ag, and Au (100) surfaces,” Physical Review B, vol. 74, no. 8, article 085408, 2006. View at: Publisher Site  Google Scholar
 C. M. Chang and C. M. Wei, “Selfdiffusion of adatoms and dimers on fcc(100) surfaces,” Chinese Journal of Physics, vol. 43, no. 1 II, pp. 169–175, 2005. View at: Google Scholar
 R. Stumpf and M. Scheffler, “Ab initio calculations of energies and selfdiffusion on flat and stepped surfaces of Al and their implications on crystal growth,” Physical Review B, vol. 53, no. 8, pp. 4958–4973, 1996. View at: Google Scholar
 H. M. Polatoglou, M. Methfessel, and M. Scheffler, “Vacancyformation energies at the (111) surface and in bulk Al, Cu, Ag, and Rh,” Physical Review B, vol. 48, no. 3, pp. 1877–1883, 1993. View at: Publisher Site  Google Scholar
 G. Boisvert, L. J. Lewis, M. J. Puska, and R. M. Nieminen, “Energetics of diffusion on the (100) and (111) surfaces of Ag, Au, and Ir from first principles,” Physical Review B, vol. 52, no. 12, pp. 9078–9085, 1995. View at: Publisher Site  Google Scholar
 G. Boisvert and L. J. Lewis, “Selfdiffusion of adatoms, dimers, and vacancies on Cu(100),” Physical Review B, vol. 56, no. 12, pp. 7643–7655, 1997. View at: Google Scholar
 G. Boisvert, L. J. Lewis, and M. Scheffler, “Island morphology and adatom selfdiffusion on Pt(111),” Physical Review B, vol. 57, no. 3, pp. 1881–1889, 1998. View at: Google Scholar
 C. Lee, G. T. Barkema, M. Breeman, A. Pasquarello, and R. Car, “Diffusion mechanism of Cu adatoms on a Cu(001) surface,” Surface Science, vol. 306, no. 3, pp. L575–L578, 1994. View at: Google Scholar
 M. Giesen, “Step and island dynamics at solid/vacuum and solid/liquid interfaces,” Progress in Surface Science, vol. 68, no. 1–3, pp. 1–153, 2001. View at: Publisher Site  Google Scholar
 G. A. Somorjai and M. A. Van Hove, “Adsorbateinduced restructuring of surfaces,” Progress in Surface Science, vol. 30, no. 34, pp. 201–231, 1989. View at: Google Scholar
 D. M. Kolb, “Reconstruction phenomena at metalelectrolyte interfaces,” Progress in Surface Science, vol. 51, no. 2, pp. 109–173, 1996. View at: Publisher Site  Google Scholar
 T. Jacob, “Potentialinduced lifting of the Au(100)surface reconstruction studied with DFT,” Electrochimica Acta, vol. 52, no. 6, pp. 2229–2235, 2007. View at: Publisher Site  Google Scholar
 S. Venkatachalam, P. Kaghazchi, L. A. Kibler, D. M. Kolb, and T. Jacob, “First principles studies of the potentialinduced lifting of the Au(100) surface reconstruction,” Chemical Physics Letters, vol. 455, no. 1–3, pp. 47–51, 2008. View at: Publisher Site  Google Scholar
 K. Pötting, W. Schmickler, and T. Jacob, “Selfdiffusion on Au(100): a density functional theory study,” ChemPhysChem, vol. 11, no. 7, pp. 1395–1404, 2010. View at: Publisher Site  Google Scholar
 G. Rosenfeld, Habilitation Thesis, University of Bonn, 1998.
 J. A. Keith, D. Fantauzzi, T. Jacob, and A. C. T. Van Duin, “Reactive forcefield for simulating gold surfaces and nanoparticles,” Physical Review B, vol. 81, no. 23, Article ID 235404, 2010. View at: Publisher Site  Google Scholar
 E. W. Hansen and M. Neurock, “Firstprinciplesbased Monte Carlo simulation of ethylene hydrogenation kinetics on Pd,” Journal of Catalysis, vol. 196, no. 2, pp. 241–252, 2000. View at: Publisher Site  Google Scholar
 J. G. O. Ojwang, S. Chaudhuri, A. C. T. Van Duin et al., “Multiscale modeling of interaction of alane clusters on Al(111) surfaces: a reactive force field and infrared absorption spectroscopy approach,” Journal of Chemical Physics, vol. 132, no. 8, Article ID 084509, 2010. View at: Publisher Site  Google Scholar
 T. Verstraelen, B. M. Szyja, D. Lesthaeghe et al., “Multilevel modeling of silicatemplate interactions during initial stages of zeolite synthesis,” Topics in Catalysis, vol. 52, no. 9, pp. 1261–1271, 2009. View at: Publisher Site  Google Scholar
 W. Goddard, B. Merinov, A. Van Duin et al., “Multiparadigm multiscale simulations for fuel cell catalysts and membranes,” Molecular Simulation, vol. 32, no. 34, pp. 251–268, 2006. View at: Publisher Site  Google Scholar
 T. AlaNissila, R. Ferrando, and S. C. Ying, “Collective and single particle diffusion on surfaces,” Advances in Physics, vol. 51, no. 3, pp. 949–1078, 2002. View at: Publisher Site  Google Scholar
 S. Günther, E. Kopatzki, M. C. Bartelt, J. W. Evans, and R. J. Behm, “Anisotropy in nucleation and growth of twodimensional islands during homoepitaxy on “hex” Reconstructed Au(100),” Physical Review Letters, vol. 73, no. 4, pp. 553–556, 1994. View at: Publisher Site  Google Scholar
 A. C. T. Van Duin, S. Dasgupta, F. Lorant, and W. A. Goddard, “ReaxFF: a reactive force field for hydrocarbons,” Journal of Physical Chemistry A, vol. 105, no. 41, pp. 9396–9409, 2001. View at: Publisher Site  Google Scholar
 A. C. T. van Duin, Computational Methods in Catalysis and Materials Science, WileyVCH, Weinheim, Germany, 2009.
 K. D. Nielson, A. C. T. Van Duin, J. Oxgaard, W. Q. Deng, and W. A. Goddard, “Development of the ReaxFF reactive force field for describing transition metal catalyzed reactions, with application to the initial stages of the catalytic formation of carbon nanotubes,” Journal of Physical Chemistry A, vol. 109, no. 3, pp. 493–499, 2005. View at: Publisher Site  Google Scholar
 J. Tersoff, “Modeling solidstate chemistry: interatomic potentials for multicomponent systems,” Physical Review B, vol. 39, no. 8, pp. 5566–5568, 1989. View at: Publisher Site  Google Scholar
 D. W. Brenner, “Empirical potential for hydrocarbons for use in simulating the chemical vapor deposition of diamond films,” Physical Review B, vol. 42, no. 15, pp. 9458–9471, 1990. View at: Publisher Site  Google Scholar
 L. Pauling, “Atomic radii and interatomic distances in metals,” Journal of the American Chemical Society, vol. 69, no. 3, pp. 542–553, 1947. View at: Google Scholar
 M. S. Daw and M. I. Baskes, “Semiempirical, quantum mechanical calculation of hydrogen embrittlement in metals,” Physical Review Letters, vol. 50, no. 17, pp. 1285–1288, 1983. View at: Publisher Site  Google Scholar
 M. S. Daw and M. I. Baskes, “Embeddedatom method: derivation and application to impurities, surfaces, and other defects in metals,” Physical Review B, vol. 29, no. 12, pp. 6443–6453, 1984. View at: Publisher Site  Google Scholar
 S. M. Foiles, M. I. Baskes, and M. S. Daw, “Embeddedatommethod functions for the fcc metals Cu, Ag, Au, Ni, Pd, Pt, and their alloys,” Physical Review B, vol. 33, no. 12, pp. 7983–7991, 1986. View at: Publisher Site  Google Scholar
 S. M. Foiles, M. I. Baskes, and M. S. Daw, “Embeddedatommethod functions for the fcc metals Cu, Ag, Au, Ni, Pd, Pt, and their alloys,” Physical Review B, vol. 37, no. 17, p. 10378, 1988. View at: Publisher Site  Google Scholar
 M. I. Baskes, “Modified embeddedatom potentials for cubic materials and impurities,” Physical Review B, vol. 46, no. 5, pp. 2727–2742, 1992. View at: Publisher Site  Google Scholar
 C. J. Casewit, K. S. Colwell, and A. K. Rappé, “Application of a Universal force field to organic molecules,” Journal of the American Chemical Society, vol. 114, no. 25, pp. 10035–10046, 1992. View at: Google Scholar
 A. K. Rappé, C. J. Casewit, K. S. Colwell, W. A. Goddard, and W. M. Skiff, “UFF, a full periodic table force field for molecular mechanics and molecular dynamics simulations,” Journal of the American Chemical Society, vol. 114, no. 25, pp. 10024–10035, 1992. View at: Google Scholar
 B. R. Brooks, R. E. Bruccoleri, B. D. Olafson, D. J. States, S. Swaminathan, and M. Karplus, “CHARMM: a program for macromolecular energy, minimization, and dynamics calculations,” Journal of Computational Chemistry, vol. 4, no. 2, pp. 187–217, 1983. View at: Google Scholar
 W. Damm, A. Frontera, J. TiradoRives, and W. L. Jorgensen, “OPLS allatom force field for carbohydrates,” Journal of Computational Chemistry, vol. 18, no. 16, pp. 1955–1970, 1997. View at: Google Scholar
 W. D. Cornell, P. Cieplak, C. I. Bayly et al., “A second generation force field for the simulation of proteins, nucleic acids, and organic molecules,” Journal of the American Chemical Society, vol. 117, no. 19, pp. 5179–5197, 1995. View at: Publisher Site  Google Scholar
 W. D. Cornell, P. Cieplak, C. I. Bayly et al., “A second generation force field for the simulation of proteins, nucleic acids, and organic molecules,” Journal of the American Chemical Society, vol. 118, no. 9, p. 2309, 1996. View at: Google Scholar
 K. Binder and D. P. Landau, A Guide to Monte Carlo Simulations in Statistical Physics, Cambridge University Press, Cambridge, UK, 2000.
 A. Gross, Theoretical Surface Science: A Microscopic Perspective, Springer, Heidelberg, Germany, 2009.
 K. A. Fichthorn and W. H. Weinberg, “Theoretical foundations of dynamical Monte Carlo simulations,” The Journal of Chemical Physics, vol. 95, no. 2, pp. 1090–1096, 1991. View at: Google Scholar
 S. P. Meyn and R. L. Twedee, Markov Chains and Stochastic Stability, Springer, 1993.
 M. Giesen, C. Steimer, and H. Ibach, “What does one learn from equilibrium shapes of twodimensional islands on surfaces?” Surface Science, vol. 471, no. 1–3, pp. 80–100, 2001. View at: Publisher Site  Google Scholar
 H. Ibach, Physics of Surfaces and Interfaces, Springer, Heidelberg, Germany, 2006.
 H. P. Bonzel, “3D equilibrium crystal shapes in the new light of STM and AFM,” Physics Reports, vol. 385, no. 12, pp. 1–67, 2003. View at: Publisher Site  Google Scholar  MathSciNet
 K. Pötting, N. B. Luque, P. M. Quaino, H. Ibach, and W. Schmickler, “Island dynamics on charged silver electrodes: kinetic MonteCarlo simulations,” Electrochimica Acta, vol. 54, no. 19, pp. 4494–4500, 2009. View at: Publisher Site  Google Scholar
 M. Giesen and G. Schulze IckingKonert, “Equilibrium fluctuations and decay of step bumps on vicinal Cu (111) surfaces,” Surface Science, vol. 412413, pp. 645–656, 1998. View at: Google Scholar
 J. C. Heyraud and J. J. Métois, “Equilibrium shape and temperature; Lead on graphite,” Surface Science, vol. 128, no. 23, pp. 334–342, 1983. View at: Google Scholar
 C. Bombis and H. Ibach, “Island equilibrium shape and shape fluctuations on the reconstructed Au(100) surface,” Surface Science, vol. 564, no. 1–3, pp. 201–210, 2004. View at: Publisher Site  Google Scholar
 S. Dieluweit, H. Ibach, and M. Giesen, “Potential dependence of step and kink energies on Au(100) electrodes in sulfuric acid,” Faraday Discussions, vol. 121, pp. 27–42, 2002. View at: Google Scholar
 S. Dieluweit and M. Giesen, “STM studies on the island dynamics on Au(100) electrodes in sulfuric acid,” Journal of Physics Condensed Matter, vol. 14, no. 16, pp. 4211–4225, 2002. View at: Publisher Site  Google Scholar
 S. Baier, S. Dieluweit, and M. Giesen, “Step and island dynamics on Cu(111), Ag(111) and Au(100) electrodes in electrolyte,” Surface Science, vol. 502503, pp. 463–473, 2002. View at: Publisher Site  Google Scholar
 P. Stoltze, “Simulation of surface defects,” Journal of Physics: Condensed Matter, vol. 6, no. 45, pp. 9495–9517, 1994. View at: Publisher Site  Google Scholar
 L. Vitos, H. L. Skriver, and J. Kollár, “Formation energy for steps and kinks on cubic transition metal surfaces,” Surface Science, vol. 425, no. 2, pp. 212–223, 1999. View at: Publisher Site  Google Scholar
 D. C. Schlößer, L. K. Verheij, G. Rosenfeld, and G. Comsa, “Determination of step free energies from island shape fluctuations on metal surfaces,” Physical Review Letters, vol. 82, no. 19, pp. 3843–3846, 1999. View at: Google Scholar
 K. Kleiner, Theoretische modellierung der selbsdiffusion auf Au(100), Diploma thesis, University of Ulm, Ulm, Germany, 2010.
Copyright
Copyright © 2011 Karin Kleiner et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.