Plasma modelling & reactor design

Similar documents
Influence of vibrational kinetics in a low pressure capacitively coupled hydrogen discharge

COMPARISON OF VHF, RF AND DC PLASMA EXCITATION FOR a-si:h DEPOSITION WITH HYDROGEN DILUTION

Optical properties of semiconductors. Dr. Katarzyna Skorupska

ANALYTICAL MODEL FOR THE BYPASS VALVE IN A LOOP HEAT PIPE

The UCLA/NICADD Plasma Density Transition Trapping Experiment

ONE. The Earth-atmosphere system CHAPTER

A novel helicon plasma source for negative ion beams for fusion

VI. Electrokinetics. Lecture 31: Electrokinetic Energy Conversion

A novel helicon plasma source for negative ion beams for fusion

Hydrogen in amorphous and microcrystalline silicon films prepared by hydrogen dilution

MODULE 2: DIFFUSION LECTURE NO DIFFUSION COEFFICIENT: MEASUREMENT AND PREDICTION

Vacuum. Gas kinetics, pumping, gauges

DRIFT AND HALL MOBILITY OF HOLE CARRIERS IN STRAINED SIGE FILMS GROWN ON (001) SI SUBSTRATES

Lecture 6 Plasmas. Chapters 10 &16 Wolf and Tauber. ECE611 / CHE611 Electronic Materials Processing Fall John Labram 1/68

The Effect of Oxygen Induced Degradation on Charge Carrier Dynamics in

Plasma Chamber. Fortgeschrittenes Praktikum I. Supervisors: Baran Eren, Dr. Marco Wisse, Dr. Laurent Marot. Abstract

Document Version Publisher s PDF, also known as Version of Record (includes final page, issue and volume numbers)

Optical self-energy of superconducting Pb in the terahertz region

Lecture 28: Kinetics of Oxidation of Metals: Part 1: rusting, corrosion, and

Factors Effect on the Saturation Parameter S and there Influences on the Gain Behavior of Ytterbium Doped Fiber Amplifier

Study of Current Decay Time during Disruption in JT-60U Tokamak

ETCHING Chapter 10. Mask. Photoresist

Session 12 : Monopropellant Thrusters

Modern Methods in Heterogeneous Catalysis Research: Preparation of Model Systems by Physical Methods

Chapter 9 Practical cycles

Large Storage Window in a-sinx/nc-si/a-sinx Sandwiched Structure

A global (volume averaged) model of a chlorine discharge

Modelling of non-uniform DC driven glow discharge in argon gas

Plasma deposition of thin film silicon at low substrate temperature and at high growth rate

Technology for Micro- and Nanostructures Micro- and Nanotechnology

Distribution of populations in excited states of electrodeless discharge lamp of Rb atoms

Transport Phenomena Coupled by Chemical Reactions in Methane Reforming Ducts

THE SEEBECK COEFFICIENT OF TiO 2 THIN FILMS

Plasma Deposition (Overview) Lecture 1

Lecture contents. Metals: Drude model Conductivity frequency dependence Plasma waves Difficulties of classical free electron model

Hong Young Chang Department of Physics, Korea Advanced Institute of Science and Technology (KAIST), Republic of Korea

Dust formation and charging in an Ar/SiH4 radiofrequency

MODULE 2: DIFFUSION LECTURE NO. 2

MODELING AND SIMULATION OF REFORMER AUTO- THERMAL REACTOR IN AMMONIA UNIT

k T m 8 B P m k T M T

Modelling heat transfer and fluid flow inside a pressure cooker

High performance microchanneled asymmetric

Flow Velocity Measurement Principles of Hot Film Anemometry

Setting up the Mathematical Model Review of Heat & Material Balances

Section 6: Ion Implantation. Jaeger Chapter 5

Internal Energy in terms of Properties

- A spark is passed through the Argon in the presence of the RF field of the coil to initiate the plasma

hal , version 1-7 May 2008

Introduction to Silvaco ATHENA Tool and Basic Concepts in Process Modeling Part - 1. Instructor: Dragica Vasileska

Hybrid Quantum Processors: Molecular Ensembles as Quantum Memory for Solid State Circuits

Effects of cross field diffusion in a low pressure high density oxygen/silane plasma

Chemistry 531 Spring 2009 Problem Set 6 Solutions

1. INTRODUCTION 2. EXPERIMENTAL SET-UP CHARACTERIZATION OF A TUBULAR PLASMA REACTOR WITH EXTERNAL ANNULAR ELECTRODES

Fundamental Mechanisms, Predictive Modeling, and Novel Aerospace Applications of Plasma Assisted Combustion

EE143 Fall 2016 Microfabrication Technologies. Lecture 6: Thin Film Deposition Reading: Jaeger Chapter 6

Characterization of an Oxygen Plasma by Using a Langmuir Probe in an Inductively Coupled Plasma

Homogeneous and Inhomogeneous Model for Flow and Heat Transfer in Porous Materials as High Temperature Solar Air Receivers

Chapter 2 Introductory Concepts of Wave Propagation Analysis in Structures

AT 25 C! CH10090 Thermodynamics (part 2) Enthalpy changes during reactions. Let s remember what we did in CH10089

Repetition: Practical Aspects

Chapter 1 Fundamentals

EE C245 ME C218 Introduction to MEMS Design Fall 2007

All about sparks in EDM

Plasma Chemistry Study in an Inductively Coupled Dielectric Etcher

I. INTRODUCTION. 127 J. Vac. Sci. Technol. B 15(1), Jan/Feb X/97/15(1)/127/6/$ American Vacuum Society 127

Enhancing the Efficiency of Silicon-Based Solar Cells by the

Simulation of the PIR detector active function

H loss mechanism during anneal of silicon nitride: Chemical dissociation

Chemical Vapor Deposition (CVD)

Efficiency of Microwave Heating of Weakly Loaded Polymeric Nanocomposites

EE143 Fall 2016 Microfabrication Technologies. Lecture 7: Ion Implantation Reading: Jaeger Chapter 5

High Temperature Catalysis

Analysis of High-Altitude Ionization Gauge Measurements Using the Direct Simulation Monte Carlo Method

The Vacuum Sorption Solution

CVD: General considerations.

The chemistry of a CCl 2 F 2 radio frequency discharge

Radial Basis Function Networks: Algorithms

The directivity of the forced radiation of sound from panels and openings including the shadow zone

First Name Last Name CIRCLE YOUR LECTURE BELOW: Div. 1 10:30 am Div. 2 2:30 pm Div. 3 4:30 pm Prof. Gore Prof. Udupa Prof. Chen

NARROW GAP ELECTRONEGATIVE CAPACITIVE DISCHARGES AND STOCHASTIC HEATING

2D Hybrid Fluid-Analytical Model of Inductive/Capacitive Plasma Discharges

Production and destruction of CF x radicals in radio-frequency fluorocarbon plasmas

Proportional Counters

III. Electromagnetic uniformity: finite RF wavelength in large area, VHF reactors: standing waves, and telegraph effect

PARTICLE CONTROL AT 100 nm NODE STATUS WORKSHOP: PARTICLES IN PLASMAS

Frequency variation under constant power conditions in hydrogen radio frequency discharges

Effect of a dielectric layer on plasma uniformity in high frequency electronegative capacitive discharges

3.155J/6.152J Microelectronic Processing Technology Fall Term, 2004

Lecture 8, the outline

Hotwire-assisted Atomic Layer Deposition of Pure Metals and Metal Nitrides

Spectroscopic investigation of Tm 3+ containing Lithium borate glasses

ELEMENT2 High Resolution- ICP-MS INSTRUMENT OVERVIEW

Investigation of InP etching mechanisms in a Cl 2 /H 2 inductively coupled plasma by optical emission spectroscopy

Plasma-Surface Interactions in Patterning High-k k Dielectric Materials

Structure of 11 Be studied in β-delayed neutron- and γ- decay from polarized 11 Li

Plasma based modification of thin films and nanoparticles. Johannes Berndt, GREMI,Orléans

Physique des plasmas radiofréquence Pascal Chabert

Semiconductor Nanocrystals from Nonthermal Plasmas. Rebecca J. Anthony University of Minnesota

DEPOSITION OF THIN TiO 2 FILMS BY DC MAGNETRON SPUTTERING METHOD

HiBarSens: Tunable diode laser spectroscopy for ultra barrier measurement

Transcription:

Plasma modelling & reactor design Alan Howling, Ben Strahm, Christoh Hollenstein Ecole Polytechnique Fédérale de Lausanne (EPFL), Switzerland Center for Research in Plasma Physics (CRPP) 1) The "lasma dimension" for lasma deosition of μc-si:h 2) Direct uming of a lasma reactor... and one non-intrusive lasma diagnostic for each. IWTFSSC-2 Berlin Aril 2009

1) Plasma-Enhanced Chemical Vaour Deosition of μc-si:h Φ H2 silane hydrogen Φ SiH4 black box PECVD

1) PECVD reactor setu P RF [W] f [MHz] Φ H2 silane hydrogen Φ SiH4 Plasma Box, Oerlikon rincile: Jacques Schmitt vacuum chamber RF connection gas flow showerhead RF electrode butterfly valve lasma box side view rocess ums heating grounded box glass substrate

c 1 1) Hydrogen dilution for lasma deosition of μc-si:h Φ H2 0.05 0 silane hydrogen Φ SiH4 low lasma box ΦSiH silane concentration c = 4 < 5% Φ is commonly used total

1) Reason for hydrogen dilution Need high ratio of Hto SiH x fluxes to deosit μc-si:h intuitively "Add more hydrogen than silane" Selective etching model Chemical annealing model Surface diffusion model

1)...but the otimal lasma control arameters are less clear c 1 Φ H2 f [MHz] P RF [W] 0.05 0 silane hydrogen Φ SiH4 ga Φ tot? area [mbar] lasma black box: arameter inter-deendence

1)...because the lasma is comlex. Φ H2 Φ SiH4 P RF [W] f [MHz] { } molecules H SiH { n electron diss : n } e, Te { y sec reactions } 2 4 SiH, H, cations, anions x olysilanes, dust 2 4 H SiH olysilanes, dust { transort in the lasma and the sheath} { surface reactions } SiH x H, H2

1) Consider low ressures (< 2 mbar) Φ H2 Φ SiH4 P RF [W] f [MHz] { } molecules H SiH { n electron diss : n } e, Te 2 4 SiH, H, cations, anions x H 2 SiH4 { transort in the lasma and the sheath} { } surface reactions SiH x H, H 2

1) Consider only the majority secies Φ H2 { molecules } H 2 SiH4 H 2 SiH4 Φ SiH4 All reactive secies (SiH x, H, and ions) have very low density because of their volume and surface reactions H 2 and SiH 4 have the dominant artial ressures in the lasma reactor Only H 2 and SiH 4 leave the lasma reactor The artial ressure of SiH 4 (and H 2 ) is ~the same in the uming line as in the lasma "The lasma comosition (and deosition) is determined by the artial ressures of SiH 4 and H 2 in the lasma"

1) How to measure the silane artial ressure in the lasma? "inaccessible" lasma reactor vacuum chamber load lock gate valve lasma box heated closed showerhead uming line butterfly valve rocess ums substrate loading door

1) FTIR in um line: Non-intrusive diagnostic for SiH 4 ressure Fourier transform infrared (FTIR) absortion sectroscoy via ZnSe windows. IR source 1 m single ass extended uming line IR detector

1) Silane artial ressure measurement a) Silane Q-branch absorbance; b) Integrate sectrum; c) Read off a calibrated curve Silane ressure with lasma < silane ressure without lasma, SiH 0 SiH <, 4 4 due to electron dissociation of silane ( irreversible loss). Hydrogen artial ressure increases with lasma (surface recombination, & um seed adjustment)

1) Silane inut concentration (without lasma) Φ H2 Φ SiH4 Define 0 0 SiH H + = 4 2 0 0 SiH H + = 4 2 silane inut concentration, Note 0 c 1 c = 0 SiH 4

1) Silane concentration with lasma Φ H2 Φ SiH4 + = SiH Define silane concentration with lasma, H 4 2 c l SiH = 4 + = SiH 4 2 H

1) Silane deletion due to lasma Φ H2 Φ SiH4 Define + = SiH H 4 2 + = SiH H 4 2 silane fractional deletion, Note 0 D 1 D = 0 SiH 4 4 0 SiH 4 SiH

1) Silane concentration in lasma Φ H2 Φ SiH4 Therefore + = SiH 4 2 ( ) silane concentration in lasma, c = c 1 D H l + = SiH 4 2 H 2 arameters, c & D, define the lasma comosition

1) The lasma "black box" becomes a {c,d} unit box 1 ure silane the "lasma dimension" silane inut concentration, c vanishingly weak lasma concentration infinitely strong lasma infinitely diluted silane 0 0 silane deletion fraction, D 1

1) Contours of constant lasma comosition silane inut concentration, c 1 same lasma comosition for ( ) c = c 1 D = const. l 0 0 silane deletion fraction, D 1

1) Contours of constant lasma comosition silane inut concentration, c 1 ( c = 0.05 ) D = 0.1 ( ) cl = c 1 D = 0.05 0 0 silane deletion fraction, D 1 5% silane 95% hydrogen strong inut dilution, weak deletion IN THE PLASMA ( c = 1 ) D = 0.95 5% silane IN THE 95% hydrogen PLASMA NO inut dilution, 95% deletion

1) Same lasma comosition = same film roerties Constant lasma comosition, c l, gives constant film roerties

1) Same lasma comosition = same film roerties transition material

1) Same lasma comosition = same film roerties c l = 1.2% c l = 0.5% Microcrystalline silicon can be deosited even with high silane concentration, rovided that the silane deletion fraction is sufficiently high, because then the lasma is dominated by hydrogen.

1) Two regimes of oeration silane inut concentration, c 1 REGIME 1 "conventional" : Film crystallinity governed by hydrogen dilution of silane, 'indeendent' of lasma. 0 0 silane deletion fraction, D 1 REGIME 2 "recent" : Film crystallinity governed by strong lasma deletion, even for high silane concentration.

1) Simle lasma chemistry model Based on a reduced set of gas hase reactions k1 e + SiH 4 SiH 2 + 2H + e e + + + e e + kh H2 2H + e e + + e and simlified surface reactions. SiH + 2 S Sisurf H2 H + H surf H H + vacancy H R 1 H H2 2 2 + vacancy surf +

1) Simle, analytical lasma chemistry model IN Plasma reactions & deosition OUT Φ Φ = kn n + a n SiH e SiH SiH 4 4 4 ( 2 ) + Sn + R n = k n n + an H SiH H H e H H 2 2 2 2 { radical fluxes } Sn SiH 2 Rn H { surface reaction} deosition H+H recombination A zero-dimensional model is aroriate to showerhead reactors

1) Simle, analytical lasma chemistry model Just one silane dissociation channel here: e + SiH 4 SiH 2 + 2H + e More detailed lasma chemistry only changes the numerical constants:... the general conclusions remain the same. radical flux ratio Γ Γ Rn k nh 2 k 1 = = 2 + 2 = 2 1 + 2 Sn k n k c H H H H SiH SiH SiH l 2 2 4 Hydrogen and silane are the dominant artial ressures; so we can exect radical densities to deend on them. "This shows why the lasma deosition is determined by c l, the silane concentration in the lasma."

lasma black box f [MHz] P RF [W] ga Φ tot 1) Deletion accounts for all of the lasma arameters? area [mbar] analytical model D = 1+ akn e ( 1+ c) 1 Φ -1 6 Tgas total a s = effective uming seed = 6.1 10 ; ga area ( [ ]) -1 kn e s = silane dissociation rate = function of PRF, f MHz ; ΦSiH 4 c = silane inut concentration,. Φ total { } { Φ } D if any of, ga, area, c, P, f &/or, T RF total gas

1) The "lasma dimension" for lasma deosition of μc-si:h 2) Direct uming of a lasma reactor IWTFSSC-2 Berlin Aril 2009

2) Plasma chemistry equilibration time to steady-state deletion van den Donker et al (2005) Feitknecht et al (2002) SiH normalized intensity Almost a minute to reach the lasma comosition suitable to grow μc-si:h... 0 20 40 60 time (s) time (s) μc-si:h... thick amorhous incubation layer... deteriorates solar cell erformance.

2) OES: Non-intrusive, raid diagnostic for equilibration time lasma emission um grid Otical Emission Sectrometer

2) Time-resolved otical emission sectrum from lasma ignition Emission from silane radicals, SiH*, falls as the silane becomes deleted Emission from excited molecular and atomic hydrogen rises as its artial ressure rises The lasma chemistry equilibration time is less than one second!

2) Comare oen and closed reactors (numerical) Oen laboratory reactor (axial symmetry) Closed, directly-umed lasma box (lateral view) Plasma zone (20 cm in diameter) Dead volume 57 cm Plasma zone indirect uming t=0-100 s time to steady-state > 100 s! direct uming t=0-1 s time to steady-state < 1 s! Silane back-diffusion from any intermediate dead volume increases the equilibration time.

2) Practical consequence of the uming configuration Amorhous incubation layers with oen door 409 nm 96 nm lasma box side view <10 nm su bs tra te But the door is normally closed, the amorhous incubation layer is then thinner than 10 nm.

Two Conclusions 1. Plasma comosition and deosition deend on the silane concentration in the lasma, c l =c(1-d), and not only on the silane concentration in the inut flow, c. Strong hydrogen dilution in the lasma, and μc-si:h deosition, can be obtained with high inut concentration of silane. - monitored non-intrusively by FTIR in the uming line. 2. Raid equilibration of lasma chemistry requires a closed, directly-umed showerhead reactor with a uniform lasma - avoid gas circulation between the lasma and any dead volumes. - monitored non-intrusively by OES in the uming line. Refs. B. Strahm et al, Plasma Sources, Sci. Technol. 16 80 (2007). A. A. Howling et al, Plasma Sources, Sci. Technol. 16 679 (2007). Acknowledgements: Swiss Federal Research Grant CTI 6947.1; and Dr. U. Kroll and colleagues, Oerlikon-Solar Lab SA, Neuchâtel, Switzerland.