Stress Evolution during Film Growth: CTC! Stress Evolution during Film Growth: CTC!

Size: px
Start display at page:

Download "Stress Evolution during Film Growth: CTC! Stress Evolution during Film Growth: CTC!"

Transcription

1 When Vapor Deforms Metal: Thermodynamics of Deposition Flux Dependent Intrinsic Film Stress M.J. Rost Kamerlingh Onnes Laboratory, Leiden University (NL) Stress Evolution during Film Growth: CTC! stress x thickness Deposition OFF Shull & Spaepen, J. Appl. Phys. 80 (996 ) 643 Stress Evolution during Film Growth: CTC! Capillary Forces Shull & Spaepen, J. Appl. Phys. 80 (996 ) 643 nucleation & 3D island growth (Volmer-Weber Growth)

2 Stress Evolution during Film Growth: CTC! Boundary Zipping Shull & Spaepen, J. Appl. Phys. 80 (996 ) 643 coalescence & grain boundary formation Stress Evolution during Film Growth: CTC! Shull & Spaepen, J. Appl. Phys. 80 (996 ) 643 film growth Further Surprise: Reversible Stress Jumps deposition interrupts Shull & Spaepen, J. Appl. Phys. 80 (996 ) 643 film growth

3 Deposition Interrupt: Tensile Relaxation Deposition OFF deposition interrupts Shull & Spaepen, J. Appl. Phys. 80 (996 ) 643 deposition interrupt Deposition Resumption: Jump to Initial Stress State deposition interrupts Shull & Spaepen, J. Appl. Phys. 80 (996 ) 643 deposition resumption Why is this happening? pre-coalescence surface tension continuation combined with ongoing grain growth 3

4 Boundary Adatom Insertion Model: Increase of ace Chemical Potential Diffusion of Atoms into Boundaries (GBs) Deposition OFF ON Development of Compressive Stress in s Boundary Adatom Insertion Model: Decrease of ace Chemical Potential (atoms incorporate at steps and kinks) Deposition OFF Atoms Diffuse out of Boundaries Relaxation of Compressive Stress Can this be true??? Can this be True: Values? ~ 50 MPa ~ 0. ML/s Aim: testing the GB insertion model thermodynamically deposition interrupts Shull & Spaepen, J. Appl. Phys. 80 (996 ) 643 deposition resumption 4

5 Thermodynamic Equilibrium GB GB GB Thermodynamic Equilibrium GB GB GB deposition interrupts Chemical Potential of the ace O,,,... O,,,... surface adatoms step adatoms kinks steps surface adatoms step adatoms kinks steps 5

6 Chemical Potential of the ace O,,,... O,,,... surface adatoms step adatoms kinks steps surface adatoms step adatoms kinks steps s N TV, N N N μ due to Adatoms Fs Uadatom Uinter S T q: Coverage Adatom Lattice Model Ad. Latt. Ead ktln Ignoring Interactions between Adatoms D Adatom Gas Model DAG Ead kt ln Aad h thermal de Broglie wavelength mkt μ due to Adatoms Difference in absolute value of entropy Same slope for q < 0.ML Important: interest in Dμ For q < 0.ML: adadtom kt ln Eq. Ad Ead kt. Latt. ln Ead kt ln Aad DAG 6

7 adadtom kt ln Eq. distribution of during growth? Island Nucleation Kinetic Monte Carlo Simulation of Film Growth: with courtesy of Vladimir Kaganer, Paul-Drude-Institut für Festkörperelektronik (see also Step Flow Growth 7

8 Step Flow Growth Kinetic Monte Carlo Simulation of Film Growth: with courtesy of Vladimir Kaganer, Paul-Drude-Institut für Festkörperelektronik (see also Adatom Density on a Terrace E att E form E diff E ES energy landscape θ n Fˆ adatom density 0 V s ν d ν e n N Adatom Density on a Terrace n n n ˆ d V s en F t n n diffusion step movement evaporation deposition flux Simplifications: room temperature: e 0 for typical deposition situations: Vs 0 θ n adatom density 0 E att E form V s Fˆ ν d N E diff ν e E ES energy landscape n 8

9 Definitions: Result: exp E / kt eq form s s0 exp EES / kt Adatom Density on a Terrace FN ˆ ( an)( sn) Fn ˆ n eq d asn a s d Cu(): 0 Hz E ev diff s0 5 EES 0.4eV Eatt 0eV E 0.74eV form exp d 0 Ediff kt aexp E / kt att Eq. experimentally determined θ n adatom density 0 E att E form V s Fˆ ν d N E diff ν e E ES energy landscape n adadtom kt ln FN ˆ ( an)( sn) Fn ˆ n eq d asn a s d Eq. Eq. chemical potential of a grain under stress? Wt grain E 0 grain Chemical Potential of Interior stress-strain relation Fgrain W t lim N N t0 TV,, Nt 0 t t ( ) 0 A dl W E lim lim t t grain t0 N t0 t L0 E E grain lim tlim t t0 L t0 0 L0 t L 0 σ Gr σ Gr+ dσ σ Gr+ dσ Eq. 3 grain 9

10 Chemical Potential of Interior Eq. 3 grain L 0 σ Gr σ Gr+ dσ t σ Gr+ dσ R. Sandström, J. Hallgren, J. Nuc. Mater. 4 (0) 5 adadtom kt ln FN ˆ ( an)( sn) Fn ˆ n eq d asn a s d Eq. Eq. Eq. 3 grain Combining the Equations kt ln ˆ ˆ FN( an )( sn ) Fn n eq asn a s d kt FN ˆ ( an )( sn ) Fn ˆ comp ln asn a s d eq d eq which terrace position n? d 0

11 ace Equilibrium Structure Phys.Rev.Lett. 9 (003) 060 from 5 0 C to 75 0 C in 394 min. [film coarsening] ace Equilibrium Structure Phys.Rev.Lett. 9 (003) 060 surface in equilibrium grains with convex shape

12 Ehrlich-Schwoebel barrier / Step Flow Zeno Effect Phys.Rev.Lett. 99 (007) 660 the last terrace position counts kt FN ˆ ( an )( sn ) Fn ˆ comp ln asn a s d eq d eq n N calculate for different widths N θ n adatom density 0 E att E form V s Fˆ ν d E diff ν e E ES energy landscape n N Boundary

13 the result the result downward funneling M.Giesen, G. Schulze Icking-Konert, H.Ibach, Phys.Rev.Lett. 8 (999) 30 3

14 downward funneling DE ES =0 for N 6 M.Giesen & H.Ibach,.Sci. 464 (000) L697 funneling in the vicinity of grain boundaries lowering the barrier for atoms to diffuse towards/into the GB, decreases the stress N kt ln eq DE ES =0.4eV Boundary DE ES =0eV the final result best fit: w eff =5 only last 3 terraces show funneling 4

15 Roughening instead of Step Flow Growth Kinetic Monte Carlo Simulation of Film Growth: with courtesy of Vladimir Kaganer, Paul-Drude-Institut für Festkörperelektronik (see also Roughening instead of Step Flow Growth Summary A.Saedi & M.J. Rost, Nature Communication 06 huge stress entropic effect ~ 50 MPa ~ 0. ML/s 5

Lecture 30: Kinetics of Epitaxial Growth: Surface Diffusion and

Lecture 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 information

Kinetic Monte Carlo: from transition probabilities to transition rates

Kinetic Monte Carlo: from transition probabilities to transition rates 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

More information

Surface Physics Surface Diffusion. Assistant: Dr. Enrico Gnecco NCCR Nanoscale Science

Surface 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 information

Phase-field modeling of step dynamics. University of California, Irvine, CA Caesar Research Center, Friedensplatz 16, 53111, Bonn, Germany.

Phase-field modeling of step dynamics. University of California, Irvine, CA Caesar Research Center, Friedensplatz 16, 53111, Bonn, Germany. Phase-field modeling of step dynamics ABSTRACT J.S. Lowengrub 1, Zhengzheng Hu 1, S.M. Wise 1, J.S. Kim 1, A. Voigt 2 1 Dept. Mathematics, University of California, Irvine, CA 92697 2 The Crystal Growth

More information

Kinetics. Rate of change in response to thermodynamic forces

Kinetics. 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 information

DIFFUSION AND GROWTH ON STRAINED SURFACES

DIFFUSION 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 information

Lecture 10 Thin Film Growth

Lecture 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 information

Lecture 1 Modeling and simulation for the growth of thin films

Lecture 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 information

Multiscale Modeling of Epitaxial Growth Processes: Level Sets and Atomistic Models

Multiscale 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 information

QUASI-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 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 information

3.320 Lecture 23 (5/3/05)

3.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 information

Diffusion in Dilute Alloys

Diffusion in Dilute Alloys Chapter 3 Diffusion in Dilute Alloys Our discussion of the atomistic mechanisms for diffusion has been confined to the situations where the diffusing species is chemically identical to the host atom. We

More information

INFLUENCE 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) 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 information

Mound slope and shape selection during unstable multilayer growth: Analysis of step-dynamics models including downward funneling

Mound slope and shape selection during unstable multilayer growth: Analysis of step-dynamics models including downward funneling Physics and Astronomy Publications Physics and Astronomy 3-2006 Mound slope and shape selection during unstable multilayer growth: Analysis of step-dynamics models including downward funneling Maozhi Li

More information

Pattern formation by step edge barriers: The growth of spirals and wedding cakes

Pattern formation by step edge barriers: The growth of spirals and wedding cakes Pattern formation by step edge barriers: The growth of spirals and wedding cakes Joachim Krug Institut für Theoretische Physik, Universität zu Köln MRS Fall Meeting, Boston, 11/26/2007 Kinetic growth modes

More information

Reconstruction and intermixing in thin Ge layers on Si 001

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 information

DYNAMICS OF THE EARLY STAGES IN METAL-ON-INSULATOR THIN FILM DEPOSITION

DYNAMICS OF THE EARLY STAGES IN METAL-ON-INSULATOR THIN FILM DEPOSITION Linköping Studies in Science and Technology Licentiate Thesis No. 1687 DYNAMICS OF THE EARLY STAGES IN METAL-ON-INSULATOR THIN FILM DEPOSITION Bo Lü Nanoscale Engineering Division Department of Physics,

More information

Step Bunching in Epitaxial Growth with Elasticity Effects

Step Bunching in Epitaxial Growth with Elasticity Effects Step Bunching in Epitaxial Growth with Elasticity Effects Tao Luo Department of Mathematics The Hong Kong University of Science and Technology joint work with Yang Xiang, Aaron Yip 05 Jan 2017 Tao Luo

More information

Effects of substrate rotation in oblique-incidence metal(100) epitaxial growth

Effects of substrate rotation in oblique-incidence metal(100) epitaxial growth W Recently, using a simplified ballistic deposition model 6 which includes the effects of downward funneling (DF) 8 and rapid diffusion on () facets, two of us have shown that many of the qualitative and

More information

Island 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 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 information

Supplementary Materials

Supplementary Materials Supplementary Materials Wedding Cake Growth Mechanism in One-Dimensional and Two-Dimensional Nanostructure Evolution Xin Yin, 1 Jian Shi, 1,2 Xiaobin Niu, 3,4 Hanchen Huang, 3 Xudong Wang 1 1. Department

More information

Christian Ratsch, UCLA

Christian 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 information

Adatom density kinetic Monte Carlo: A hybrid approach to perform epitaxial growth simulations

Adatom 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 information

Phase-field modeling of nanoscale island dynamics

Phase-field modeling of nanoscale island dynamics Title of Publication Edited by TMS (The Minerals, Metals & Materials Society), Year Phase-field modeling of nanoscale island dynamics Zhengzheng Hu, Shuwang Li, John Lowengrub, Steven Wise 2, Axel Voigt

More information

Monte Carlo Simulation of Pulsed Laser Deposition

Monte 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 information

Kinetic Monte Carlo Simulation of Epitaxial Thin Film Growth: Formation of Submonolayer Islands and Multilayer Mounds

Kinetic Monte Carlo Simulation of Epitaxial Thin Film Growth: Formation of Submonolayer Islands and Multilayer Mounds Chemistry Conference Papers, Posters and Presentations Chemistry 8-2007 Kinetic Monte Carlo Simulation of Epitaxial Thin Film Growth: Formation of Submonolayer Islands and Multilayer Mounds James W. Evans

More information

Morphological evolution of single-crystal ultrathin solid films

Morphological 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 information

Development and ordering of mounds during metal(100) homoepitaxy

Development and ordering of mounds during metal(100) homoepitaxy Development and ordering of mounds during metal(100) homoepitaxy K. Caspersen, A. R. Layson, C. R. Stoldt, V. Fournée, P. A. Thiel, J. W. Evans To cite this version: K. Caspersen, A. R. Layson, C. R. Stoldt,

More information

Short-range attraction, surface currents, and mound formation in metal 111 epitaxial growth

Short-range attraction, surface currents, and mound formation in metal 111 epitaxial growth Short-range attraction, surface currents, and mound formation in metal 111 epitaxial growth Jianguo Yu and Jacques G. Amar Department of Physics & Astronomy, University of Toledo, Toledo, Ohio 43606, USA

More information

From Initial to Late Stages of Epitaxial Thin Film Growth: STM Analysis and Atomistic or CoarseGrained Modeling

From Initial to Late Stages of Epitaxial Thin Film Growth: STM Analysis and Atomistic or CoarseGrained Modeling Chemistry Conference Papers, Posters and Presentations Chemistry 8-2010 From Initial to Late Stages of Epitaxial Thin Film Growth: STM Analysis and Atomistic or CoarseGrained Modeling James W. Evans Iowa

More information

Phase transitions for particles in R 3

Phase transitions for particles in R 3 Phase transitions for particles in R 3 Sabine Jansen LMU Munich Konstanz, 29 May 208 Overview. Introduction to statistical mechanics: Partition functions and statistical ensembles 2. Phase transitions:

More information

Universality of step bunching behavior in systems with non-conserved dynamics

Universality of step bunching behavior in systems with non-conserved dynamics Universality of step bunching behavior in systems with non-conserved dynamics Joachim Krug Institute for Theoretical Physics, University of Cologne Electromigration-induced step bunching on Si(111) [Yang,

More information

Heinonen, J.; Koponen, I.; Merikoski, J.; Ala-Nissilä, Tapio Island diffusion on metal fcc(100) surfaces

Heinonen, 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 information

Nucleation theory and the early stages of thin film growth

Nucleation 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 information

Plan for Lectures #4, 5, & 6. Theme Of Lectures: Nano-Fabrication

Plan 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 information

Monte Carlo simulation of thin-film growth on a surface with a triangular lattice

Monte 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 information

Two simple lattice models of the equilibrium shape and the surface morphology of supported 3D crystallites

Two 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 information

Module 16. Diffusion in solids II. Lecture 16. Diffusion in solids II

Module 16. Diffusion in solids II. Lecture 16. Diffusion in solids II Module 16 Diffusion in solids II Lecture 16 Diffusion in solids II 1 NPTEL Phase II : IIT Kharagpur : Prof. R. N. Ghosh, Dept of Metallurgical and Materials Engineering Keywords: Micro mechanisms of diffusion,

More information

Chapter 5. On-line resource

Chapter 5. On-line resource Chapter 5 The water-air heterogeneous system On-line resource on-line analytical system that portrays the thermodynamic properties of water vapor and many other gases http://webbook.nist.gov/chemistry/fluid/

More information

Why 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 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 information

Level set approach to reversible epitaxial growth

Level set approach to reversible epitaxial growth PHYSICAL REVIEW E, VOLUME 64, 061602 Level set approach to reversible epitaxial growth M. Petersen, 1,2 C. Ratsch, 2 R. E. Caflisch, 2 and A. Zangwill 1 1 School of Physics, Georgia Institute of Technology,

More information

Atomic Transport & Phase Transformations Lecture III-2

Atomic Transport & Phase Transformations Lecture III-2 Atomic Transport & Phase Transformations Lecture III-2 PD Dr. Nikolay Zotov zotov@imw.uni-stuttgart.de Atomic Transport & Phase Transformations Part III Lectures Solid State Reactions Short Description

More information

submitted to Physical Review E Numerical Tests of Nucleation Theories for the Ising Models Seunghwa Ryu 1 and Wei Cai 2

submitted to Physical Review E Numerical Tests of Nucleation Theories for the Ising Models Seunghwa Ryu 1 and Wei Cai 2 submitted to Physical Review E Numerical Tests of Nucleation Theories for the Ising Models Seunghwa Ryu 1 and Wei Cai 2 1 Department of Physics, Stanford University, Stanford, California 94305 2 Department

More information

Formation and coarsening of Ag(110) bilayer islands on NiAl(110): STM analysis and atomistic lattice-gas modeling

Formation 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 information

Lecture 6 Surface Diffusion Driven by Surface Energy

Lecture 6 Surface Diffusion Driven by Surface Energy Lecture 6 Surface Diffusion Driven by Surface Energy Examples Flattening a surface. Spherodizing. Rayleigh instability. Grain boundary grooving. Sintering Self-assembled quantum dots Atomic Flux and Surface

More information

Modellierung molekularer Prozesse beim Wachstum organischer Schichten

Modellierung molekularer Prozesse beim Wachstum organischer Schichten Modellierung molekularer Prozesse beim Wachstum organischer Schichten Slide 1 Motivation OLED para Sexiphenyl (6P) (C36H26) OFET Pentacene (5A) (C22H14) Slide 2 Outline Methods and Materials Cohesive,

More information

Diffusion in multicomponent solids. Anton Van der Ven Department of Materials Science and Engineering University of Michigan Ann Arbor, MI

Diffusion in multicomponent solids. Anton Van der Ven Department of Materials Science and Engineering University of Michigan Ann Arbor, MI Diffusion in multicomponent solids nton Van der Ven Department of Materials Science and Engineering University of Michigan nn rbor, MI Coarse graining time Diffusion in a crystal Two levels of time coarse

More information

Kinetic Monte Carlo simulation of semiconductor quantum dot growth

Kinetic 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 information

115.3 Assignment #9 Solutions

115.3 Assignment #9 Solutions 115. Assignment #9 Solutions-1 115. Assignment #9 Solutions 8.1-12 Solve the differential equation d dx = 2(1 ), where 0 = 2 for x 0 = 0. d 1 = 2dx d 1 = 2dx ln 1 =2x + C Find C b inserting the Initial

More information

A Reduced Order Model for Epitaxial Growth

A Reduced Order Model for Epitaxial Growth Contemporary Mathematics A Reduced Order Model for Epitaxial Growth Russel E. Caflisch and David G. Meyer This paper is dedicated to Stan Osher on the occasion of his 6th birthday. Abstract. We formulate

More information

Kinetic 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 information

Oxide growth model. Known as the Deal-Grove or linear-parabolic model

Oxide growth model. Known as the Deal-Grove or linear-parabolic model Oxide growth model Known as the Deal-Grove or linear-parabolic model Important elements of the model: Gas molecules (oxygen or water) are incident on the surface of the wafer. Molecules diffuse through

More information

LECTURE 11: Monte Carlo Methods III

LECTURE 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 information

Supplementary Figure 1 Transient absorption (TA) spectrum pumped at 400 nm in the FAPbI3 sample with different excitation intensities and initial

Supplementary Figure 1 Transient absorption (TA) spectrum pumped at 400 nm in the FAPbI3 sample with different excitation intensities and initial Supplementary Figure 1 Transient absorption (TA) spectrum pumped at 400 nm in the FAPbI3 sample with different excitation intensities and initial carrier concentrations: (a) N0 = 4.84 10 18 cm -3 ; (c)

More information

Thermodynamic evolution of phase explosion during high-power nanosecond laser ablation

Thermodynamic evolution of phase explosion during high-power nanosecond laser ablation Thermodynamic evolution of phase explosion during high-power nanosecond laser ablation Quanming Lu* School of Earth and Space Sciences, University of Science and Technology of China, Hefei, 230026, China

More information

Kinetic Monte Carlo Simulation of Two-dimensional Semiconductor Quantum Dots Growth

Kinetic 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 information

Uniform-acceptance force-biased Monte Carlo: A cheap way to boost MD

Uniform-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 information

Molecular Dynamics Simulations of Fusion Materials: Challenges and Opportunities (Recent Developments)

Molecular 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 information

Ripening of monolayer vacancy pits on metal surfaces: Pathways, energetics, and size-scaling for Ag(111) versus Ag(100)

Ripening 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 information

Lecture 1: Atomic Diffusion

Lecture 1: Atomic Diffusion Part IB Materials Science & Metallurgy H. K. D. H. Bhadeshia Course A, Metals and Alloys Lecture 1: Atomic Diffusion Mass transport in a gas or liquid generally involves the flow of fluid (e.g. convection

More information

Kinetic 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 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 information

Monte Carlo Simulation of SrTiO 3 Thin Film Growth

Monte Carlo Simulation of SrTiO 3 Thin Film Growth Monte Carlo Simulation of SrTiO 3 Thin Film Growth Alexander Zakjevskii, Carolyn Worley Advisor: Dr. Anter El Azab Department of Scientific Computing Florida State University Abstract Monte Carlo simulations

More information

.O. Demokritov niversität Münster, Germany

.O. Demokritov niversität Münster, Germany Quantum Thermodynamics of Magnons.O. Demokritov niversität Münster, Germany Magnon Frequency Population BEC-condensates http://www.uni-muenster.de/physik/ap/demokritov/ k z k y Group of NonLinea Magnetic

More information

Vacuum Technology and film growth. Diffusion Resistor

Vacuum Technology and film growth. Diffusion Resistor Vacuum Technology and film growth Poly Gate pmos Polycrystaline Silicon Source Gate p-channel Metal-Oxide-Semiconductor (MOSFET) Drain polysilicon n-si ion-implanted Diffusion Resistor Poly Si Resistor

More information

First passage time Markov chain analysis of rare events for kinetic Monte Carlo: double kink nucleation during dislocation glide

First 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 information

Structure analysis: Electron diffraction LEED TEM RHEED

Structure analysis: Electron diffraction LEED TEM RHEED Structure analysis: Electron diffraction LEED: Low Energy Electron Diffraction SPA-LEED: Spot Profile Analysis Low Energy Electron diffraction RHEED: Reflection High Energy Electron Diffraction TEM: Transmission

More information

Rubber elasticity. Marc R. Roussel Department of Chemistry and Biochemistry University of Lethbridge. February 21, 2009

Rubber elasticity. Marc R. Roussel Department of Chemistry and Biochemistry University of Lethbridge. February 21, 2009 Rubber elasticity Marc R. Roussel Department of Chemistry and Biochemistry University of Lethbridge February 21, 2009 A rubber is a material that can undergo large deformations e.g. stretching to five

More information

In#uence of microstructure of substrate surface on early stage of thin "lm growth

In#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 information

Surface stress and relaxation in metals

Surface stress and relaxation in metals J. Phys.: Condens. Matter 12 (2000) 5541 5550. Printed in the UK PII: S0953-8984(00)11386-4 Surface stress and relaxation in metals P M Marcus, Xianghong Qian and Wolfgang Hübner IBM Research Center, Yorktown

More information

Diffusion in Potential Fields: Time-Dependent Capture on Radial and Rectangular Substrates. J. A. Venables 1,2 and P. Yang 1

Diffusion in Potential Fields: Time-Dependent Capture on Radial and Rectangular Substrates. J. A. Venables 1,2 and P. Yang 1 Diffusion in Potential Fields: Time-Dependent Capture on Radial and Rectangular Substrates J. A. Venables 1,2 and P. Yang 1 1 Department of Physics & Astronomy, Arizona State University, Tempe AZ 85287-1504,

More information

Introduction Critical state models Pinning regimes Kinds of pinning sites HTS Results on YBCO Conclusions. Flux pinning.

Introduction Critical state models Pinning regimes Kinds of pinning sites HTS Results on YBCO Conclusions. Flux pinning. Department of Physics and Astronomy 14.6.2011 Contents Introduction Critical state models Pinning regimes Kinds of pinning sites HTS Results on YBCO Type II superconductors and vortices Type I ξ < λ S/N

More information

arxiv: v1 [cond-mat.stat-mech] 23 Aug 2011

arxiv: v1 [cond-mat.stat-mech] 23 Aug 2011 arxiv:1108.4655v1 [cond-mat.stat-mech] 23 Aug 2011 Kinetic modelling of epitaxial film growth with upand downward step barriers F.F. Leal 1,2, T. J. Oliveira 1 and S. C. Ferreira 1 1 Departamento de Física

More information

Optimizing Graphene Morphology on SiC(0001)

Optimizing Graphene Morphology on SiC(0001) Optimizing Graphene Morphology on SiC(0001) James B. Hannon Rudolf M. Tromp Graphene sheets Graphene sheets can be formed into 0D,1D, 2D, and 3D structures Chemically inert Intrinsically high carrier mobility

More information

Off-Lattice KMC Simulation of Quantum Dot Formation

Off-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 information

7. Kinetics controlled by fluctuations: Kramers theory of activated processes

7. Kinetics controlled by fluctuations: Kramers theory of activated processes 7. Kinetics controlled by fluctuations: Kramers theory of activated processes Macroscopic kinetic processes (time dependent concentrations) Elementary kinetic process Reaction mechanism Unimolecular processes

More information

Equilibria in Materials

Equilibria in Materials 009 fall Advanced Physical Metallurgy Phase Equilibria in Materials 10. 13. 009 Eun Soo Park Office: 33-316 Telephone: 880-71 Email: espark@snu.ac.kr Office hours: by an appointment 1 Contents for previous

More information

Nanowires and nanorods

Nanowires and nanorods Nanowires and nanorods One-dimensional structures have been called in different ways: nanowires, nanorod, fibers of fibrils, whiskers, etc. These structures have a nanometer size in one of the dimensions,

More information

arxiv:cond-mat/ v2 28 Jul 1999

arxiv: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 information

Available online at Physics Procedia 15 (2011) Stability and Rupture of Alloyed Atomic Terraces on Epitaxial Interfaces

Available 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 information

"Enhanced Layer Coverage of Thin Films by Oblique Angle Deposition"

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 information

On the correlation between the self-organized island pattern and substrate elastic anisotropy

On 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 information

STANDARD SAMPLE. Reduced section " Diameter. Diameter. 2" Gauge length. Radius

STANDARD SAMPLE. Reduced section  Diameter. Diameter. 2 Gauge length. Radius MATERIAL PROPERTIES TENSILE MEASUREMENT F l l 0 A 0 F STANDARD SAMPLE Reduced section 2 " 1 4 0.505" Diameter 3 4 " Diameter 2" Gauge length 3 8 " Radius TYPICAL APPARATUS Load cell Extensometer Specimen

More information

Three-dimensional simulations of self-assembly of hut shaped Si-Ge quantum dots

Three-dimensional simulations of self-assembly of hut shaped Si-Ge quantum dots University of Massachusetts Amherst From the SelectedWorks of Ashwin Ramasubramaniam 2004 Three-dimensional simulations of self-assembly of hut shaped Si-Ge quantum dots Ashwin Ramasubramaniam, University

More information

Modeling and computational studies of surface shape dynamics of ultrathin single-crystal films with anisotropic surface energy, part I

Modeling and computational studies of surface shape dynamics of ultrathin single-crystal films with anisotropic surface energy, part I Modeling and computational studies of surface shape dynamics of ultrathin single-crystal films with anisotropic surface energy, part I Mikhail Khenner Department of Mathematics, Western Kentucky University

More information

PHYSICAL SELF-ASSEMBLY AND NANO-PATTERNING*

PHYSICAL 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 information

Electrical conduction in solids

Electrical conduction in solids Equations of motion Electrical conduction in solids Electrical conduction is the movement of electrically charged particles through a conductor or semiconductor, which constitutes an electric current.

More information

J 12 J 23 J 34. Driving forces in the nano-magnetism world. Intra-atomic exchange, electron correlation effects: Inter-atomic exchange: MAGNETIC ORDER

J 12 J 23 J 34. Driving forces in the nano-magnetism world. Intra-atomic exchange, electron correlation effects: Inter-atomic exchange: MAGNETIC ORDER Driving forces in the nano-magnetism world Intra-atomic exchange, electron correlation effects: LOCAL (ATOMIC) MAGNETIC MOMENTS m d or f electrons Inter-atomic exchange: MAGNETIC ORDER H exc J S S i j

More information

Thermally Activated Processes

Thermally Activated Processes General Description of Activated Process: 1 Diffusion: 2 Diffusion: Temperature Plays a significant role in diffusion Temperature is not the driving force. Remember: DRIVING FORCE GRADIENT of a FIELD VARIALE

More information

Nucleation, Growth, and Relaxation of Thin Films: Metal(100) Homoepitaxial Systems

Nucleation, Growth, and Relaxation of Thin Films: Metal(100) Homoepitaxial Systems Chemistry Publications Chemistry 2000 Nucleation, Growth, and Relaxation of Thin Films: Metal(100) Homoepitaxial Systems Patricia A. Thiel Iowa State University, thiel@ameslab.gov James W. Evans Iowa State

More information

Molecular Dynamics and Accelerated Molecular Dynamics

Molecular 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 information

Surface stress-induced self-assembly

Surface stress-induced self-assembly Surface stress-induced self-assembly Zhigang Suo Work with Wei Lu Solid surface Liquid surface Stretched solid surface Energy per surface atom depends on stretch Stretched liquid surface Energy per surface

More information

Heteroepitaxial growth of Co/ Cu 001 : An accelerated molecular dynamics simulation study

Heteroepitaxial growth of Co/ Cu 001 : An accelerated molecular dynamics simulation study PHYSICAL REVIEW B 72, 035415 2005 Heteroepitaxial growth of Co/ Cu 001 : An accelerated molecular dynamics simulation study Radu A. Miron Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6,

More information

of Long-Term Transport

of 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 information

Growth morphology evolution in real time and real space

Growth morphology evolution in real time and real space Growth morphology evolution in real time and real space Joachim Krug Institut für Theoretische Physik, Universität zu Köln Multilayer growth modes and interlayer transport The Ehrlich-Schwoebel effect

More information

Ostwald ripening of two-dimensional homoepitaxial islands

Ostwald ripening of two-dimensional homoepitaxial islands PHYSICAL REVIEW B 72, 035437 2005 Ostwald ripening of two-dimensional homoepitaxial islands Frank Haußer* and Axel Voigt Crystal Growth Group, Research Center caesar, Ludwig-Erhard-Allee 2, 53175 Bonn,

More information

Supporting Online Material (1)

Supporting Online Material (1) Supporting Online Material The density functional theory (DFT) calculations were carried out using the dacapo code (http://www.fysik.dtu.dk/campos), and the RPBE (1) generalized gradient correction (GGA)

More information

Magnetic domain theory in dynamics

Magnetic domain theory in dynamics Chapter 3 Magnetic domain theory in dynamics Microscale magnetization reversal dynamics is one of the hot issues, because of a great demand for fast response and high density data storage devices, for

More information

Lecture 6: Ideal gas ensembles

Lecture 6: Ideal gas ensembles Introduction Lecture 6: Ideal gas ensembles A simple, instructive and practical application of the equilibrium ensemble formalisms of the previous lecture concerns an ideal gas. Such a physical system

More information

N = N A Pb A Pb. = ln N Q v kt. = kt ln v N

N = N A Pb A Pb. = ln N Q v kt. = kt ln v N 5. Calculate the energy for vacancy formation in silver, given that the equilibrium number of vacancies at 800 C (1073 K) is 3.6 10 3 m 3. The atomic weight and density (at 800 C) for silver are, respectively,

More information

Set 3: Thermal Physics

Set 3: Thermal Physics Set 3: Thermal Physics Equilibrium Thermal physics describes the equilibrium distribution of particles for a medium at temperature T Expect that the typical energy of a particle by equipartition is E kt,

More information