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

Similar documents
Lecture 30: Kinetics of Epitaxial Growth: Surface Diffusion and

Kinetic Monte Carlo: from transition probabilities to transition rates

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

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

Kinetics. Rate of change in response to thermodynamic forces

DIFFUSION AND GROWTH ON STRAINED SURFACES

Lecture 10 Thin Film Growth

Lecture 1 Modeling and simulation for the growth of thin films

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

QUASI-EQUILIBRIUM MONTE-CARLO: OFF-LATTICE KINETIC MONTE CARLO SIMULATION OF HETEROEPITAXY WITHOUT SADDLE POINTS

3.320 Lecture 23 (5/3/05)

Diffusion in Dilute Alloys

INFLUENCE OF THE ELECTROCHEMICAL POTENTIAL ON ENERGY LANDSCAPES NEAR STEP AND ISLAND EDGES: Ag(100) and Ag(111)

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

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

Reconstruction and intermixing in thin Ge layers on Si 001

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

Step Bunching in Epitaxial Growth with Elasticity Effects

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

Island Dynamics Model for Mound Formation: The Effect of a Step-Edge Barrier

Supplementary Materials

Christian Ratsch, UCLA

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

Phase-field modeling of nanoscale island dynamics

Monte Carlo Simulation of Pulsed Laser Deposition

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

Morphological evolution of single-crystal ultrathin solid films

Development and ordering of mounds during metal(100) homoepitaxy

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

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

Phase transitions for particles in R 3

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

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

Nucleation theory and the early stages of thin film growth

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

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

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

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

Chapter 5. On-line resource

Why Heteroepitaxy? Difficult to Model Multiple Species Dislocations and point defects important Species Flux important

Level set approach to reversible epitaxial growth

Atomic Transport & Phase Transformations Lecture III-2

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

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

Lecture 6 Surface Diffusion Driven by Surface Energy

Modellierung molekularer Prozesse beim Wachstum organischer Schichten

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

Kinetic Monte Carlo simulation of semiconductor quantum dot growth

115.3 Assignment #9 Solutions

A Reduced Order Model for Epitaxial Growth

Kinetic Monte Carlo (KMC)

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

LECTURE 11: Monte Carlo Methods III

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

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

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

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

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

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

Lecture 1: Atomic Diffusion

Kinetic lattice Monte Carlo simulations of diffusion processes in Si and SiGe alloys

Monte Carlo Simulation of SrTiO 3 Thin Film Growth

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

Vacuum Technology and film growth. Diffusion Resistor

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

Structure analysis: Electron diffraction LEED TEM RHEED

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

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

Surface stress and relaxation in metals

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

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

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

Optimizing Graphene Morphology on SiC(0001)

Off-Lattice KMC Simulation of Quantum Dot Formation

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

Equilibria in Materials

Nanowires and nanorods

arxiv:cond-mat/ v2 28 Jul 1999

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

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

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

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

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

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

PHYSICAL SELF-ASSEMBLY AND NANO-PATTERNING*

Electrical conduction in solids

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

Thermally Activated Processes

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

Molecular Dynamics and Accelerated Molecular Dynamics

Surface stress-induced self-assembly

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

of Long-Term Transport

Growth morphology evolution in real time and real space

Ostwald ripening of two-dimensional homoepitaxial islands

Supporting Online Material (1)

Magnetic domain theory in dynamics

Lecture 6: Ideal gas ensembles

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

Set 3: Thermal Physics

Transcription:

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

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

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

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

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

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

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 https://www.youtube.com/watch?v=nsgrksv8yh8) Step Flow Growth 7

Step Flow Growth Kinetic Monte Carlo Simulation of Film Growth: with courtesy of Vladimir Kaganer, Paul-Drude-Institut für Festkörperelektronik (see also https://www.youtube.com/watch?v=nsgrksv8yh8) 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

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 0 0.040 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

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

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

Ehrlich-Schwoebel barrier / Step Flow 9 9 9 7 7 7 5 5 5 3 3 3 0 0 0 30 40 50 60 70 0 0 0 30 40 50 60 70 0 0 0 30 40 50 60 70 9 9 9 7 7 7 5 5 5 3 3 3 0 0 0 30 40 50 60 70 0 0 0 30 40 50 60 70 0 0 0 30 40 50 60 70 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

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

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

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 https://www.youtube.com/watch?v=nsgrksv8yh8) 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