Kinetic Monte Carlo: from transition probabilities to transition rates
|
|
- Sharyl Dixon
- 6 years ago
- Views:
Transcription
1 Kinetic Monte Carlo: from transition probabilities to transition rates With MD we can only reproduce the dynamics of the system for 100 ns. Slow thermallyactivated processes, such as diffusion, cannot be modeled. Metropolis Monte Carlo samples configurational space and generates configurations according to the desired statistical-mechanics distribution. However, there is no time in Metropolis MC and the method cannot be used to study evolution of the system or kinetics. An alternative computational technique that can be used to study kinetics of slow processes is the kinetic Monte Carlo (kmc) method. Cartoon by Larry Gonick As compared to the Metropolis Monte Carlo method, kinetic Monte Carlo method has a different scheme for the generation of the next state. The main idea behind kmc is to use transition rates that depend on the energy barrier between the states, with time increments formulated so that they relate to the microscopic kinetics of the system. In Metropolis MC methods we decide whether to accept a move by considering the energy difference between the states. In kmc methods we use rates that depend on the energy barrier between the states.
2 Kinetic Monte Carlo: kinetics of atomic rearrangements Before starting kmc simulation we have to make a list of all possible events that can be realized during the simulation and calculate rates for each event. Input to KMC: Fast processes MD simulations Slow processes transition state theory, experiments For example, diffusion on the surface is determined by the energy barriers for the breaking the adatom-substrate bonds, E a =E saddle -E min, and the rate constant for diffusion can be calculated using a simplified transition state theory, e.g. [A. Voter, Phys. Rev. B 34, 6819 (1986)]: Energy E a Distance E a k TST n p exp kt n p is the number of possible jump directions is the harmonic frequency. The assumptions are that the diffusion of an adatom is a result of random uncorrelated hops between neighboring binding sites and that the time between hops is much longer as compared to the time for the hop. For atoms adjacent to an island/step, additional (Ehrlich-Schwoebel) barriers/rates should be specified for breaking bonds with the atoms of the island. As the system becomes more complex, the number of possible events becomes larger
3 Kinetic Monte Carlo: kinetics of atomic rearrangements When the rate constants of all processes are known, we can perform kmc simulation in the time domain. In the case of a single process, the reciprocal of the rate of the process determines the time required for the reaction to occur. This quantity can be set equal to the kmc time. In the case of a many-particle multi-process system, however, introduction of time is less straightforward and several modifications of kmc exist. For example, we can use the following scheme [A. F. Voter, Phys. Rev. B 34, (1986)]: 1. The rate of all the allowed processes can be combined to obtain the total rate and the time step is calculated as inverse of the total rate. 2. At each time step one of all the possible processes is randomly selected with probability that is the product of the time step and the rate of the individual process. Other algorithms include: Assigning to each particle in the system independent time clock (and time step) and calculating the real time step as an average over all independent time steps [e.g. Pei-Lin Cao, Phys. Rev. Lett. 73, (1994)] Choosing a single constant time step that is less than the duration of the fastest process. The processes are chosen randomly and allowed to occur with probabilities based on individual rates [Dawnkaski, Srivastava, Garrison, J. Chem. Phys. 102, (1995)]
4 Example: kmc simulation of diamond {001} (2x1):H surface under CVD growth conditions [Dawnkaski, Srivastava, Garrison, Chem. Phys. Lett. 232, 524, 1994] First, the rates of all relevant surface reaction have to be evaluated (from molecular static or dynamics simulations i 0 i exp( E a i k B ) T 0 i [s -1 ] a E i [ev] University of Virginia, MSE a 4270/6270: few Introduction more reactions to Atomistic Simulations, Leonid Zhigilei
5 Example: kmc simulation of diamond {001} (2x1):H surface under CVD growth conditions Given the rate of each individual process in a system, the probability of a process occurring within any specified time period or timestep is simply the product of the timestep and the rate of the process. A single constant timestep can be chosen so that it is less than the duration of the fastest process considered. The probabilities of all the considered processes are thus between 0 and 1. The timestep is chosen such that the acceptance probability of the fastest process is about 0.5. This gives values of t =10-5 sfor 1200 K, 10-7 s for 1500 K, and 10-8 s for 1800 K. Topographical snapshots of a growing diamond surface. The surface area is a square of side 25 A and the plots is shaded by layer number with the highest layer being the lightest in shade. [Dawnkaski, Srivastava, Garrison, Chem. Phys. Lett. 232, 524, 1994]
6 Example: Growth of fractal structures in fullerene layers MD and kinetic Monte Carlo simulations by Hui Liu (term project for MSE 6270) All possible thermally-activated events have to be considered STM images of C 60 film growing on graphite i = o exp(-e i /k B T) Ln (V j ) -5.5 MD simulations finding the -6.0 energy barriers, attempt frequencies, and probabilities of x10 University of Virginia, MSE 4270/6270: Introduction to Atomistic -3 Simulations, diffusion Leonidjump Zhigilei events 1/kT (1/meV)
7 Example: Growth of fractal structures in fullerene layers MD and kinetic Monte Carlo simulations by Hui Liu (term project for MSE 627) STM images of C 60 film growing on graphite
8 Example: Diffusion of Ge adatoms on a reconstructed Si(001) substrate MD and kinetic Monte Carlo simulations by Avinash Dongare Calculation of the potential energy surface for a Ge adatom using procedure described by Roland and Gilmer in Phys. Rev. B 46, (1992)
9 Example: Diffusion of Ge adatoms on a reconstructed Si(001) substrate MD and kinetic Monte Carlo simulations by Avinash Dongare MD simulations of Ge adatom diffusion on a Au covered Si(001) substrate Ge adatom trajectories on a Au covered Si(001) substrate Mean square displacements for Ge adatoms on a Aucovered Si(001) substrate, and logarithmic plot of the diffusion coefficient from MD simulations 2 r ( t) ~ 4Dt Ln (D [cm 2 /sec]) D D0 exp( Ed / kbt ) D-MD /k B T
10 Example: Diffusion of Ge adatoms on a reconstructed Si(001) substrate MD and kinetic Monte Carlo simulations by Avinash Dongare Surface structures after deposition of 0.07 ML of Si on a Si(001) substrate at 400 K with a deposition rate of 0.1 ML/min kmc STM, Mo et al., Phys. Rev. Lett. 66, 1998 (1991) Surface structures predicted in kmc simulations of the deposition of 0.10 ML of Ge on a Au patterned Si(001) substrate (76 nm x 76 nm) at 600 K with deposition rate of 9 ML/min.
11 Kinetic Monte Carlo: limitations One (main?) problem in kmc is that we have to specify all the barriers/rates in advance, before the simulation. But what if we have a continuous variation of the activation energies in the system? What if the activation energies are changing during the simulation? E ( r, t) k ( E, T ) k 0 exp kt Example: strain on the surface can affect the diffusion of adatoms and nucleation of islands. There could be many possible origins of strain, e. g. buried islands, mesas, dislocation patterns in heteroepitaxial systems. One can try to introduce the effect of strain on the activation energies for the diffusion of adatoms. For example, Nurminen et al., Phys. Rev. B 63, , 2000, tried several approaches. In one approach, they introduced a spatial dependence of the adatom-substrate interaction on a patterned surface, E = E S (x,y) + ne N, where E is the diffusion activation energy, E S is the contribution due to the interaction with substrate, and E N is the energy of interaction with other adatoms. In another approach, an additional hop-direction dependent diffusion barrier E D is introduced to describe the long-range interaction between adatom and domain boundaries: University of Virginia, MSE 4270/6270: Introduction E=E to Atomistic S +ne N +E Simulations, D Leonid Zhigilei
12 Kinetic Monte Carlo: example Kinetic Monte Carlo simulation of island growth on a homogeneous substrate and a substrate with nanoscale patterning (by Nurminen, Kuronen, Kaski, Helsinki University of Technology) islands on a homogeneous substrate Blue denotes the substrate and green the deposited atoms. islands on a substrate with nanoscale patterning (a checkerboard structure) The growing heteroepitaxial islands by themselves can locally modify diffusion barriers. Relaxation of nanostructures introduces local strains that constantly change the energy landscape and corresponding probabilities of Monte Carlo events. Modified approach based on the locally activated Monte Carlo techniques (Kaukonen et al. Phys. Rev. B, 61, 980, 2000) has been proposed to account for local strains.
13 Example: kmc in simulations of dislocation dynamics [Karin Lin and D. C. Chrzan, Phys. Rev. B 60, 3799 (1999)] Simulation of evolution of the collective behavior of a large number of dislocations requires velocity vs. stress law. Typically, in Dislocation Dynamics simulations an empirical law involving a damping term is assumed. Kinetic Monte Carlo can provide the needed velocity vs. stress relationship. In particular, the simulations may reflect the stochastic aspects associated with overcoming the Peierls barrier, interactions with vacancies, etc. in a natural way. Dislocation dynamics method can be used (potentially) to connect atomic scale calculations with macroscopic continuum description of plasticity. Kinetic Monte Carlo simulations can be parameterized based on atomic scale MD studies of the properties of dislocation cores, kinks, etc. The model discussed in this paper includes dislocation segment interactions in the isotropic elasticity theory limit. The Peierls potential is also included and free surface boundary conditions are used (all image forces and surface tractions are reflected in the energetics governing the dynamics). Dislocations are assumed to be composed entirely of screw and edge segments. The algorithm involves cataloging all of the possible kinetic events, and calculating the rates associated with these processes. Kinetic events include the production/annihilation of double kink pairs, as well as the University of Virginia, MSE 4270/6270: Introduction lateral motion to Atomistic of the existing Simulations, kinks. Leonid Zhigilei
14 Summary: MD, Metropolis MC and kinetic MC With MD we can only reproduce the dynamics of the system for 100 ns. Slow thermallyactivated processes, such as diffusion, cannot be modeled. An alternative computational techniques for slow processes are Monte Carlo methods. Monte Carlo method is a common name for a wide variety of stochastic techniques. These techniques are based on the use of random numbers and probability statistics to investigate problems in areas as diverse as economics, nuclear physics, and flow of traffic. There are many variations of Monte Carlo methods. In this lecture we will briefly discuss two methods that are often used in materials science - classical Metropolis Monte Carlo and kinetic Monte Carlo. Metropolis Monte Carlo generates configurations according to the desired statistical-mechanics distribution. There is no time, the method cannot be used to study evolution of the system. Equilibrium properties can be studied. Cartoon by Larry Gonick Kinetic Monte Carlo can address kinetics. The main idea behind KMC is to use transition rates that depend on the energy barrier between the states, with time increments formulated so that they relate to the microscopic kinetics of the system.
15 Summary: MD, Metropolis MC and kinetic MC Molecular Dynamics based on the solution of the equations of motion for all particles in the system. Complete information on atomic trajectories can be obtained, but time of the simulations is limited (up to nanoseconds) appropriate for fast processes (e.g. FIB local surface modification, sputtering, implantation) or for quasi-static simulations (e.g. stress distribution in nanostructures). Metropolis Monte Carlo generates random configurations with probability of each configuration defined by the desired distribution P(r N ). This is accomplished by setting up a random walk through the configurational space with specially designed choice of probabilities of going from one state to another. Equilibrium properties can be found/studied (e.g. surface reconstruction and segregation, composition variations in the surface region due to the surface or substrate induced strains, stability of nanostructures). Kinetic Monte Carlo when the rate constants of all processes are known, we can perform kmc simulation in the time domain. Time increments are defined by the rates of all processes and are formulated so that they relate to the microscopic kinetics of the system. This method should be used when kinetics rather than equilibrium thermodynamics dominates the structural and/or compositional changes in the system. Kinetic Monte Carlo vs Metropolis Monte Carlo: in MMC we decide whether to accept a move by considering the energy difference between the states, whereas in kmc methods we use rates that depend on the energy barrier between the states. The main advantages of kinetic Monte Carlo is that time is defined and only a small number of elementary reactions are considered, so the calculations are fast.
Reconstruction and intermixing in thin Ge layers on Si 001
Reconstruction and intermixing in thin Ge layers on Si 001 L. Nurminen, 1 F. Tavazza, 2 D. P. Landau, 1,2 A. Kuronen, 1 and K. Kaski 1 1 Laboratory of Computational Engineering, Helsinki University of
More information3.320 Lecture 23 (5/3/05)
3.320 Lecture 23 (5/3/05) Faster, faster,faster Bigger, Bigger, Bigger Accelerated Molecular Dynamics Kinetic Monte Carlo Inhomogeneous Spatial Coarse Graining 5/3/05 3.320 Atomistic Modeling of Materials
More informationChristian Ratsch, UCLA
Strain Dependence of Microscopic Parameters and its Effects on Ordering during Epitaxial Growth Christian Ratsch, UCLA Institute for Pure and Applied Mathematics, and Department of Mathematics Collaborators:
More informationKinetic Monte Carlo (KMC)
Kinetic Monte Carlo (KMC) Molecular Dynamics (MD): high-frequency motion dictate the time-step (e.g., vibrations). Time step is short: pico-seconds. Direct Monte Carlo (MC): stochastic (non-deterministic)
More informationMolecular Dynamics Simulations of Fusion Materials: Challenges and Opportunities (Recent Developments)
Molecular Dynamics Simulations of Fusion Materials: Challenges and Opportunities (Recent Developments) Fei Gao gaofeium@umich.edu Limitations of MD Time scales Length scales (PBC help a lot) Accuracy of
More informationKinetic Monte Carlo simulation of nucleation on patterned substrates
PHYSICAL REVIEW B, VOLUME 63, 035407 Kinetic Monte Carlo simulation of nucleation on patterned substrates L. Nurminen, A. Kuronen, and K. Kaski Helsinki University of Technology, Laboratory of Computational
More informationWhy Heteroepitaxy? Difficult to Model Multiple Species Dislocations and point defects important Species Flux important
Why Heteroepitaxy? Difficult to Model Multiple Species Dislocations and point defects important Species Flux important The Silicon Age is Ending New materials and growth required Useful in design of Devices
More informationKinetic Monte Carlo modelling of semiconductor growth
Kinetic Monte Carlo modelling of semiconductor growth Peter Kratzer Faculty of Physics, University Duisburg-Essen, Germany Time and length scales morphology Ga As 2D islands surface reconstruction Methods
More informationLECTURE 11: Monte Carlo Methods III
1 LECTURE 11: Monte Carlo Methods III December 3, 2012 In this last chapter, we discuss non-equilibrium Monte Carlo methods. We concentrate on lattice systems and discuss ways of simulating phenomena such
More informationQUASI-EQUILIBRIUM MONTE-CARLO: OFF-LATTICE KINETIC MONTE CARLO SIMULATION OF HETEROEPITAXY WITHOUT SADDLE POINTS
QUASI-EQUILIBRIUM MONTE-CARLO: OFF-LATTICE KINETIC MONTE CARLO SIMULATION OF HETEROEPITAXY WITHOUT SADDLE POINTS Henry A. Boateng University of Michigan, Ann Arbor Joint work with Tim Schulze and Peter
More informationKinetic Monte Carlo: Coarsegraining
Kinetic Monte Carlo: Coarsegraining Time and Space Peter Kratzer Faculty of Physics, University Duisburg-Essen, Germany morphology Time and length scales Ga As 2D islands surface reconstruction Methods
More informationKinetic lattice Monte Carlo simulations of diffusion processes in Si and SiGe alloys
Kinetic lattice Monte Carlo simulations of diffusion processes in Si and SiGe alloys, Scott Dunham Department of Electrical Engineering Multiscale Modeling Hierarchy Configuration energies and transition
More informationLecture 1 Modeling and simulation for the growth of thin films
Lecture 1 Modeling and simulation for the growth of thin films Russel Caflisch Mathematics Department Materials Science and Engineering Department UCLA & IPAM www.math.ucla.edu/~material 1 Outline Epitaxial
More informationSurface Physics Surface Diffusion. Assistant: Dr. Enrico Gnecco NCCR Nanoscale Science
Surface Physics 008 8. Surface Diffusion Assistant: Dr. Enrico Gnecco NCCR Nanoscale Science Random-Walk Motion Thermal motion of an adatom on an ideal crystal surface: - Thermal excitation the adatom
More informationFirst passage time Markov chain analysis of rare events for kinetic Monte Carlo: double kink nucleation during dislocation glide
INSTITUTE OF PHYSICS PUBLISHING MODELLING AND SIMULATION IN MATERIALS SCIENCE AND ENGINEERING Modelling Simul. Mater. Sci. Eng. 10 (2002) 581 596 PII: S0965-0393(02)39232-5 First passage time Markov chain
More informationG. L. KELLOGG. Sandia National Laboratories Albuquerque, NM USA. 1. Abstract
' C' 54NB 9s - q27 r, coup -C/606222-- AN ATOMC VEW OF CLUSTER DFFUSON ON METAL SURFACES G. L. KELLOGG Sandia National Laboratories Albuquerque, NM 8785-43 USA. Abstract Field ion microscope measurements
More informationKinetic Monte Carlo (KMC) Kinetic Monte Carlo (KMC)
Kinetic Monte Carlo (KMC) Molecular Dynamics (MD): high-frequency motion dictate the time-step (e.g., vibrations). Time step is short: pico-seconds. Direct Monte Carlo (MC): stochastic (non-deterministic)
More informationTwo simple lattice models of the equilibrium shape and the surface morphology of supported 3D crystallites
Bull. Nov. Comp. Center, Comp. Science, 27 (2008), 63 69 c 2008 NCC Publisher Two simple lattice models of the equilibrium shape and the surface morphology of supported 3D crystallites Michael P. Krasilnikov
More informationMultiscale Modeling of Epitaxial Growth Processes: Level Sets and Atomistic Models
Multiscale Modeling of Epitaxial Growth Processes: Level Sets and Atomistic Models Russel Caflisch 1, Mark Gyure 2, Bo Li 4, Stan Osher 1, Christian Ratsch 1,2, David Shao 1 and Dimitri Vvedensky 3 1 UCLA,
More informationMolecular Dynamics and Accelerated Molecular Dynamics
Molecular Dynamics and Accelerated Molecular Dynamics Arthur F. Voter Theoretical Division National Laboratory Lecture 3 Tutorial Lecture Series Institute for Pure and Applied Mathematics (IPAM) UCLA September
More informationStress Evolution during Film Growth: CTC! Stress Evolution during Film Growth: CTC!
When Vapor Deforms Metal: Thermodynamics of Deposition Flux Dependent Intrinsic Film Stress M.J. Rost Kamerlingh Onnes Laboratory, Leiden University (NL) Rost@physics.leidenuniv.nl www.physics.leidenuniv.nl/rost
More informationMonte Carlo simulation of thin-film growth on a surface with a triangular lattice
Vacuum 52 (1999) 435 440 Monte Carlo simulation of thin-film growth on a surface with a triangular lattice Wei Helin*, Liu Zuli, Yao Kailun Department of Physics, Huazhong University of Science and Technology,
More informationLecture 30: Kinetics of Epitaxial Growth: Surface Diffusion and
Lecture 30: Kinetics of Epitaxial Growth: Surface Diffusion and Nucleation Today s topics Understanding the basics of epitaxial techniques used for surface growth of crystalline structures (films, or layers).
More informationKinetics and Functionality of Cu-coordinated Pyridyl-porphyrin Supramolecular Self-assembly on a Au(111) Surface
Kinetics and Functionality of Cu-coordinated Pyridyl-porphyrin Supramolecular Self-assembly on a Au(111) Surface LI, Yang, MPhil candidate, Physics, HKUST Supervisor, Prof. LIN, Nian 2012-08-08 Outline
More informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
More informationFractal Structures in Fullerene Layers: Simulation of the Growth Process
J. Phys. Chem. C 2008, 112, 4687-4695 4687 Fractal Structures in Fullerene Layers: Simulation of the Growth Process Hui Liu, Zhibin Lin, Leonid V. Zhigilei, and Petra Reinke* Department of Materials Science
More informationOff-Lattice KMC Simulation of Quantum Dot Formation
Off-Lattice KMC Simulation of Quantum Dot Formation T. P. Schulze University of Tennessee Peter Smereka & Henry Boateng University of Michigan Research supported by NSF-DMS-0854920 Heteroepitaxy and Quantum
More informationDIFFUSION AND GROWTH ON STRAINED SURFACES
DIFFUSION AND GROWTH ON STRAINED SURFACES A.C. SCHINDLER AND D.E. WOLF Theoretische Physik, FB 10, Gerhard Mercator Universität Duisburg, 47048 Duisburg, Germany E-mail: A.Schindler@Uni-Duisburg.de Using
More informationKinetics. Rate of change in response to thermodynamic forces
Kinetics Rate of change in response to thermodynamic forces Deviation from local equilibrium continuous change T heat flow temperature changes µ atom flow composition changes Deviation from global equilibrium
More informationCOMPUTATIONAL INVESTIGATION OF THE EFFECT OF CLUSTER IMPACT ENERGY ON THE MICROSTRUCTURE OF FILMS GROWN BY CLUSTER DEPOSITION
COMPUTATIONAL INVESTIGATION OF THE EFFECT OF CLUSTER IMPACT ENERGY ON THE MICROSTRUCTURE OF FILMS GROWN BY CLUSTER DEPOSITION AVINASH M. DONGARE, DEREK D. HASS, AND LEONID V. ZHIGILEI Department of Materials
More informationOn the correlation between the self-organized island pattern and substrate elastic anisotropy
JOURNAL OF APPLIED PHYSICS 100, 013527 2006 On the correlation between the self-organized island pattern and substrate elastic anisotropy E. Pan a and R. Zhu Department of Civil Engineering, University
More informationDiffusion mechanisms in Cu grain boundaries
PHYSICAL REVIEW B VOLUME 62, NUMBER 6 Diffusion mechanisms in Cu grain boundaries 1 AUGUST 2000-II Mads R. So rensen Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 Yuri
More informationNucleation theory and the early stages of thin film growth
Nucleation theory and the early stages of thin film growth C. Ratsch a) Department of Mathematics, University of California, Los Angeles, California 90095-1555 J. A. Venables Department of Physics and
More informationMonte Carlo Simulation of Pulsed Laser Deposition
Monte Carlo Simulation of Pulsed Laser Deposition Pui-Man Lam @, S. J. Liu, and C.H. Woo Department of Mechanical Engineering, The Hong Kong Polytechnic University, Hong Kong ABSTRACT Using the Monte Carlo
More informationKinetic Monte Carlo Simulation of Two-dimensional Semiconductor Quantum Dots Growth
Kinetic Monte Carlo Simulation of Two-dimensional Semiconductor Quantum Dots Growth by Ernie Pan Richard Zhu Melissa Sun Peter Chung Computer Modeling and Simulation Group The University of Akron Outline
More informationHeinonen, J.; Koponen, I.; Merikoski, J.; Ala-Nissilä, Tapio Island diffusion on metal fcc(100) surfaces
Powered by TCPDF (www.tcpdf.org) This is an electronic reprint of the original article. This reprint may differ from the original in pagination and typographic detail. Heinonen, J.; Koponen, I.; Merikoski,
More informationEnergy barriers for diffusion on stepped Pt(1 1 1) surface
Vacuum 54 (1999) 113 117 Energy barriers for diffusion on stepped Pt(1 1 1) surface F. Máca *, M. Kotrla, O.S. Trushin Institute of Physics, ASCR, Praha, Czech Republic Institute of Microelectronics, RASC,
More informationffl Mineral surfaces are not static nor inert, they are dynamic and highly interactive with their environment. ffl Foundation of surface interactions
Mineral-Fluid Interactions ffl All Global-scale phenomena begin with atomic scale dissolution reactions! ffl Mineral-fluid reactions control the dissolved load ffl which also control: quality of fresh
More informationLecture 10 Thin Film Growth
Lecture 10 Thin Film Growth 1/76 Announcements Homework: Homework Number 2 is returned today, please pick it up from me at the end of the class. Solutions are online. Homework 3 will be set Thursday (2
More informationAtomistic simulation of chemical vapor deposition of (111)-oriented diamond film using a kinetic Monte Carlo method
JOURNAL OF MATERIALS SCIENCE 34 (1999)7 20 Atomistic simulation of chemical vapor deposition of (111)-oriented diamond film using a kinetic Monte Carlo method M. GRUJICIC, S. G. LAI Department of Mechanical
More informationAtomistic simulations on the mobility of di- and tri-interstitials in Si
Atomistic simulations on the mobility of di- and tri-interstitials in Si related publications (since 2001): Posselt, M., Gao, F., Zwicker, D., Atomistic study of the migration of di- and tri-interstitials
More informationKinetic Monte Carlo method for dislocation glide in silicon
Journal of Computer-Aided Materials Design, 6: 175 183, 1999. KLUWER/ESCOM 1999 Kluwer Academic Publishers. Printed in the Netherlands. Kinetic Monte Carlo method for dislocation glide in silicon WEI CAI
More informationMorphological evolution of single-crystal ultrathin solid films
Western Kentucky University From the SelectedWorks of Mikhail Khenner March 29, 2010 Morphological evolution of single-crystal ultrathin solid films Mikhail Khenner, Western Kentucky University Available
More informationMolecular Dynamics Saddle Search Adaptive Kinetic Monte Carlo
Molecular Dynamics Saddle Search Adaptive Kinetic Monte Carlo Samuel T. Chill and Graeme Henkelman Department of Chemistry and the Institute for Computational Engineering and Sciences, The University of
More informationAu-C Au-Au. g(r) r/a. Supplementary Figures
g(r) Supplementary Figures 60 50 40 30 20 10 0 Au-C Au-Au 2 4 r/a 6 8 Supplementary Figure 1 Radial bond distributions for Au-C and Au-Au bond. The zero density regime between the first two peaks in g
More informationATOMISTIC MODELING OF BORON ACTIVATION AND DIFFUSION IN STRAINED SIGE
ATOMISTIC MODELING OF BORON ACTIVATION AND DIFFUSION IN STRAINED SIGE Scott T. Dunham,, Jakyoung Song, and Chihak Ahn Dept. of Electrical Engineering, Dept. of Physics University of Washington, Box 35500,
More informationInteraction of ion beams with matter
Interaction of ion beams with matter Introduction Nuclear and electronic energy loss Radiation damage process Displacements by nuclear stopping Defects by electronic energy loss Defect-free irradiation
More informationKinetic Monte Carlo simulation of semiconductor quantum dot growth
Solid State Phenomena Online: 2007-03-15 ISSN: 1662-9779, Vols. 121-123, pp 1073-1076 doi:10.4028/www.scientific.net/ssp.121-123.1073 2007 Trans Tech Publications, Switzerland Kinetic Monte Carlo simulation
More informationDevelopment of a Kinetic Monte Carlo Code
Development of a Kinetic Monte Carlo Code Daniel Pedersen June 25, 2013 Abstract A framework for constructing kinetic monte carlo (KMC) simulations of diffusive events on a lattice was developed. This
More informationLong-Term Atomistic Simulation of Heat and Mass Transport
Long-Term Atomistic Simulation of Heat and Mass Transport Kevin G. Wang 1, Mauricio Ponga 2, Michael Ortiz 2 1 Department of Aerospace and Ocean Engineering, Virginia Polytechnic Institute and State University
More informationMultiscale Materials Modeling
Multiscale Materials Modeling Lecture 02 Capabilities of Classical Molecular Simulation These notes created by David Keffer, University of Tennessee, Knoxville, 2009. Outline Capabilities of Classical
More informationMODELING DEFECT MEDIATED DOPANT DIFFUSION IN SILICON. by Brian T. Puchala
MODELING DEFECT MEDIATED DOPANT DIFFUSION IN SILICON by Brian T. Puchala A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy (Materials Science and
More informationDepartment of Engineering Mechanics, SVL, Xi an Jiaotong University, Xi an
The statistical characteristics of static friction J. Wang, G. F. Wang*, and W. K. Yuan Department of Engineering Mechanics, SVL, Xi an Jiaotong University, Xi an 710049, China * E-mail: wanggf@mail.xjtu.edu.cn
More informationModeling the sputter deposition of thin film photovoltaics using long time scale dynamics techniques
Loughborough University Institutional Repository Modeling the sputter deposition of thin film photovoltaics using long time scale dynamics techniques This item was submitted to Loughborough University's
More informationStructural and Mechanical Properties of Nanostructures
Master s in nanoscience Nanostructural properties Mechanical properties Structural and Mechanical Properties of Nanostructures Prof. Angel Rubio Dr. Letizia Chiodo Dpto. Fisica de Materiales, Facultad
More informationWang-Landau Monte Carlo simulation. Aleš Vítek IT4I, VP3
Wang-Landau Monte Carlo simulation Aleš Vítek IT4I, VP3 PART 1 Classical Monte Carlo Methods Outline 1. Statistical thermodynamics, ensembles 2. Numerical evaluation of integrals, crude Monte Carlo (MC)
More informationScaling behavior of one-dimensional Pt chains migration on Pt 110 1=2 surface
17 April 2000 Physics Letters A 268 2000 413 417 www.elsevier.nlrlocaterphysleta Scaling behavior of one-dimensional Pt chains migration on ž / ž / Pt 110 1=2 surface GuoCe Zhuang a,b,), Wei Wang a a National
More informationMultiscale modelling of D trapping in W
CMS Multiscale modelling of D trapping in W Kalle Heinola, Tommy Ahlgren and Kai Nordlund Department of Physics and Helsinki Institute of Physics University of Helsinki, Finland Contents Background Plasma-wall
More informationCherry-Pit Structures in Binary Immiscible Alloy Under Ion Irradiation
Cherry-Pit Structures in Binary Immiscible Alloy Under Ion Irradiation Shipeng Shu, Kenneth Tussey May 8, 2011 Abstract We study an special microstructure (matrix atom riched small clusters inside the
More informationPlan for Lectures #4, 5, & 6. Theme Of Lectures: Nano-Fabrication
Plan for Lectures #4, 5, & 6 Theme Of Lectures: Nano-Fabrication Quantum Wells, SLs, Epitaxial Quantum Dots Carbon Nanotubes, Semiconductor Nanowires Self-assembly and Self-organization Two Approaches
More informationMatSci 331 Homework 4 Molecular Dynamics and Monte Carlo: Stress, heat capacity, quantum nuclear effects, and simulated annealing
MatSci 331 Homework 4 Molecular Dynamics and Monte Carlo: Stress, heat capacity, quantum nuclear effects, and simulated annealing Due Thursday Feb. 21 at 5pm in Durand 110. Evan Reed In this homework,
More information"Enhanced Layer Coverage of Thin Films by Oblique Angle Deposition"
Mater. Res. Soc. Symp. Proc. Vol. 859E 2005 Materials Research Society JJ9.5.1 "Enhanced Layer Coverage of Thin Films by Oblique Angle Deposition" * karabt@rpi.edu Tansel Karabacak *, Gwo-Ching Wang, and
More informationINFLUENCE OF THE ELECTROCHEMICAL POTENTIAL ON ENERGY LANDSCAPES NEAR STEP AND ISLAND EDGES: Ag(100) and Ag(111)
INFLUENCE OF THE ELECTROCHEMICAL POTENTIAL ON ENERGY LANDSCAPES NEAR STEP AND ISLAND EDGES: Ag(100) and Ag(111) Michael I. Haftel* + and Theodore L. Einstein + *Code 6331, Naval Research Laboratory, Washington,
More informationFundamental science and synergy of multi-species surface interactions in high-plasma--flux environments
Fundamental science and synergy of multi-species surface interactions in high-plasma--flux environments Our thanks to John Hogan Managed by UT-Battelle for the Department of Energy ReNew PMI, UCLA, March
More informationMonte Carlo Simulation of Ferroelectric Domain Structure: Electrostatic and Elastic Strain Energy Contributions
Monte Carlo Simulation of Ferroelectric Domain Structure: Electrostatic and Elastic Strain Energy Contributions B.G. Potter, Jr., B.A. Tuttle, and V. Tikare Sandia National Laboratories Albuquerque, NM
More informationDynamical Monte-Carlo Simulation of Surface Kinetics
Dynamical Monte-Carlo Simulation of Surface Kinetics V. Guerra and J. Loureiro Centro de Física dos Plasmas, Instituto Superior Técnico, 1049-001 Lisboa, Portugal Abstract. This wor presents a theoretical
More informationIsland Dynamics Model for Mound Formation: The Effect of a Step-Edge Barrier
Island Dynamics Model for Mound Formation: The Effect of a Step-Edge Barrier Joe Papac, 1 Dionisios Margetis, 2 Frederic Gibou, 3 and Christian Ratsch 4 1 Department of Mathematics, University of California,
More informationarxiv:cond-mat/ v2 28 Jul 1999
International Journal of Modern Physics C, Vol. 10, No. 0 (1999) 000 000 c World Scientific Publishing Company A SIMPLE SOLID-ON-SOLID MODEL OF EPITAXIAL FILMS GROWTH: SURFACE ROUGHNESS AND DYNAMICS arxiv:cond-mat/9905110v2
More informationMSE 4270/6270: Introduction to Atomistic Simulations Spring 2018, Tuesday and Thursday, 2:00-3:15 pm Thornton Hall D115
MSE 4270/6270: Introduction to Atomistic Simulations Spring 2018, Tuesday and Thursday, 2:00-3:15 pm Thornton Hall D115 Contact Information: Instructor: Leonid Zhigilei Office: Wilsdorf Hall 303D Office
More informationof Long-Term Transport
Atomistic Modeling and Simulation of Long-Term Transport Phenomena in Nanomaterials K.G. Wang, M. and M. Ortiz Virginia Polytechnic Institute and State University Universidad de Sevilla California Institute
More informationFormation and coarsening of Ag(110) bilayer islands on NiAl(110): STM analysis and atomistic lattice-gas modeling
Chemistry Publications Chemistry 2010 Formation and coarsening of Ag(110) bilayer islands on NiAl(110): STM analysis and atomistic lattice-gas modeling Yong Han Iowa State University, y27h@iastate.edu
More informationSTM spectroscopy (STS)
STM spectroscopy (STS) di dv 4 e ( E ev, r) ( E ) M S F T F Basic concepts of STS. With the feedback circuit open the variation of the tunneling current due to the application of a small oscillating voltage
More informationc 2011 Society for Industrial and Applied Mathematics
MULTISCALE MODEL. SIMUL. Vol. 9, No., pp. 335 354 c 20 Society for Industrial and Applied Mathematics A MULTISCALE METHOD FOR EPITAXIAL GROWTH YI SUN, RUSSEL CAFLISCH, AND BJORN ENGQUIST Abstract. In this
More informationIn#uence of microstructure of substrate surface on early stage of thin "lm growth
Vacuum 56 (2000) 185}190 In#uence of microstructure of substrate surface on early stage of thin "lm growth Helin Wei*, Zuli Liu, Kailun Yao Department of Physics, Huazhong University of Science and Technology.
More information12. Spectral diffusion
1. Spectral diffusion 1.1. Spectral diffusion, Two-Level Systems Until now, we have supposed that the optical transition frequency of each single molecule is a constant (except when we considered its variation
More informationMolecular Dynamics & Adaptive Kinetic Monte Carlo Growth Study of Hydrocarbon Flakes
Molecular Dynamics & Adaptive Kinetic Monte Carlo Growth Study of Hydrocarbon Flakes Amit R. Sharma Department of Physics Wright State University Dayton, Ohio, USA amit.sharma@wright.edu Hydrocarbon co-deposits
More informationAvailable online at Physics Procedia 15 (2011) Stability and Rupture of Alloyed Atomic Terraces on Epitaxial Interfaces
Available online at www.sciencedirect.com Physics Procedia 15 (2011) 64 70 Stability and Rupture of Alloyed Atomic Terraces on Epitaxial Interfaces Michail Michailov Rostislaw Kaischew Institute of Physical
More informationAdvanced sampling. fluids of strongly orientation-dependent interactions (e.g., dipoles, hydrogen bonds)
Advanced sampling ChE210D Today's lecture: methods for facilitating equilibration and sampling in complex, frustrated, or slow-evolving systems Difficult-to-simulate systems Practically speaking, one is
More informationQuantum Monte Carlo Benchmarks Density Functionals: Si Defects
Quantum Monte Carlo Benchmarks Density Functionals: Si Defects K P Driver, W D Parker, R G Hennig, J W Wilkins (OSU) C J Umrigar (Cornell), R Martin, E Batista, B Uberuaga (LANL), J Heyd, G Scuseria (Rice)
More informationDensity Functional Theory: from theory to Applications
Density Functional Theory: from theory to Applications Uni Mainz May 27, 2012 Large barrier-activated processes time-dependent bias potential extended Lagrangian formalism Basic idea: during the MD dynamics
More informationIntroduction To Materials Science FOR ENGINEERS, Ch. 5. Diffusion. MSE 201 Callister Chapter 5
Diffusion MSE 201 Callister Chapter 5 1 Goals: Diffusion - how do atoms move through solids? Fundamental concepts and language Diffusion mechanisms Vacancy diffusion Interstitial diffusion Impurities Diffusion
More informationRipening of monolayer vacancy pits on metal surfaces: Pathways, energetics, and size-scaling for Ag(111) versus Ag(100)
Chemistry Publications Chemistry 2007 Ripening of monolayer vacancy pits on metal surfaces: Pathways, energetics, and size-scaling for Ag(111) versus Ag(100) Mingmin Shen Iowa State University J.-M. Wen
More informationUniform-acceptance force-biased Monte Carlo: A cheap way to boost MD
Dresden Talk, March 2012 Uniform-acceptance force-biased Monte Carlo: A cheap way to boost MD Barend Thijsse Department of Materials Science and Engineering Delft University of Technology, The Netherlands
More informationDislocations in graphene
Dislocations in graphene M. Ortiz California Institute of Technology In collaboration with: M.P. Ariza, Universidad de Sevilla Symposium on Multiscale Dislocation Dynamics UCSD, La Jolla, January 16-17,
More informationImportance Sampling in Monte Carlo Simulation of Rare Transition Events
Importance Sampling in Monte Carlo Simulation of Rare Transition Events Wei Cai Lecture 1. August 1, 25 1 Motivation: time scale limit and rare events Atomistic simulations such as Molecular Dynamics (MD)
More informationBulk Structures of Crystals
Bulk Structures of Crystals 7 crystal systems can be further subdivided into 32 crystal classes... see Simon Garrett, "Introduction to Surface Analysis CEM924": http://www.cem.msu.edu/~cem924sg/lecturenotes.html
More informationPHYSICAL SELF-ASSEMBLY AND NANO-PATTERNING*
Mater. Res. Soc. Symp. Proc. Vol. 849 2005 Materials Research Society KK8.4.1 PHYSICAL SELF-ASSEMBLY AND NANO-PATTERNING* T.-M. Lu, D.-X. Ye, T. Karabacak, and G.-C. Wang, Department of Physics, Applied
More informationScanning Tunneling Microscopy. how does STM work? the quantum mechanical picture example of images how can we understand what we see?
Scanning Tunneling Microscopy how does STM work? the quantum mechanical picture example of images how can we understand what we see? Observation of adatom diffusion with a field ion microscope Scanning
More informationChapter V: Interactions of neutrons with matter
Chapter V: Interactions of neutrons with matter 1 Content of the chapter Introduction Interaction processes Interaction cross sections Moderation and neutrons path For more details see «Physique des Réacteurs
More informationAtomistic studies of thin film growth
Reprinted from: T. S. Rahman, C. Ghosh, O. Trushin, A. Kara, and A. Karim, Atomistic Studies of Thin Film Growth, Proc. SPIE Annual Meeting 2004, 5509, 1 (2004). Atomistic studies of thin film growth Talat
More informationFirst-Passage Kinetic Monte Carlo Algorithm for Complex Reaction-Diffusion Systems
First-Passage Kinetic Monte Carlo Algorithm for Complex Reaction-Diffusion Systems Aleksandar Donev 1 Lawrence Postdoctoral Fellow Lawrence Livermore National Laboratory In collaboration with: Vasily V.
More informationThermodynamics of diffusion (extracurricular material - not tested)
Thermodynamics of diffusion (etracurricular material - not tested) riving force for diffusion iffusion in ideal and real solutions Thermodynamic factor iffusion against the concentration gradient Spinodal
More informationUpper Critical Dimension for Irreversible Cluster Nucleation and Growth. Abstract
Upper Critical Dimension for Irreversible Cluster Nucleation and Growth Feng Shi, Yunsic Shim, and Jacques G. Amar Department of Physics & Astronomy University of Toledo, Toledo, OH 43606 (Dated: December
More informationThe basics of Scanning Probe Microscopy
McGill University, Montreal, Jan. 8 th 2009 The basics of Scanning Probe Microscopy Federico Rosei Canada Research Chair in Nanostructured Organic and Inorganic Materials Énergie, Matériaux et Télécommunications,
More information2 M. N. POPESCU, F. FAMILY, AND J. G. AMAR a set of deterministic, coupled reaction-diffusion equations describing the time (coverage) dependence of a
CAPTURE-NUMBERS AND ISLAND SIZE-DISTRIBUTIONS IN IRREVERSIBLE HOMOEPITAXIAL GROWTH A Rate-Equation Approach M. N. POPESCU 1,F.FAMILY Department of Physics, Emory University, Atlanta, GA 30322 AND J. G.
More informationAdatom density kinetic Monte Carlo: A hybrid approach to perform epitaxial growth simulations
PHYSICAL REVIEW B 68, 155429 2003 Adatom density kinetic Monte Carlo: A hybrid approach to perform epitaxial growth simulations L. Mandreoli* and J. Neugebauer Fritz-Haber-Institut der Max-Planck-Gesellschaft,
More informationMultiscale modelling of electrical breakdown at high electric field
CMS HIP Multiscale modelling of electrical breakdown at high electric field Flyura Djurabekova, Helga Timkó, Aarne Pohjonen, Stefan Parviainen, Leila Costelle and Kai Nordlund Helsinki Institute of Physics
More informationAtomistic aspects: the role of DFT calculations
Atomistic aspects: the role of DFT calculations Peter Kratzer Fritz-Haber-Institut der MPG, Berlin, Germany How to obtain information about the energetics at the atomic level From statistical physics to
More informationMonte Carlo. Lecture 15 4/9/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
Monte Carlo Lecture 15 4/9/18 1 Sampling with dynamics In Molecular Dynamics we simulate evolution of a system over time according to Newton s equations, conserving energy Averages (thermodynamic properties)
More information1. Robust hexagonal rings on Cu(111) Figure S1 2. Details of Monte Carlo simulations
Supporting Information for Influence of Relativistic Effects on Assembled Structures of V-Shaped Bispyridine Molecules on M(111) Surfaces where M = Cu, Ag, Au Xue Zhang, 1,ǁ Na Li, 1, Hao Wang, 1 Chenyang
More information