Electron Temperature Modification in Gas Discharge Plasma

Similar documents
Plasma parameter evolution in a periodically pulsed ICP

DOE WEB SEMINAR,

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

Control of gas discharge plasma properties utilizing nonlocal electron energy distribution functions DOE-OFS Plasma Science Center

A Kinetic Theory of Planar Plasma Sheaths Surrounding Electron Emitting Surfaces

MODELING AND SIMULATION OF LOW TEMPERATURE PLASMA DISCHARGES

MODELING PLASMA PROCESSING DISCHARGES

Plasma-Beam Experiment with DC- RF Magnetized Plasma Discharge

Physique des plasmas radiofréquence Pascal Chabert

Study of a Micro Hollow Cathode Discharge at medium argon gas pressure

4 Modeling of a capacitive RF discharge

The distorting effect of the ion current on electron temperature measured by an electric probe

Effects of the Gas Pressure on Low Frequency Oscillations in E B Discharges

Simulation of Low Pressure Plasma Processing Reactors: Kinetics of Electrons and Neutrals

P. Diomede, D. J. Economou and V. M. Donnelly Plasma Processing Laboratory, University of Houston

MWP MODELING AND SIMULATION OF ELECTROMAGNETIC EFFECTS IN CAPACITIVE DISCHARGES

Effect of Gas Flow Rate and Gas Composition in Ar/CH 4 Inductively Coupled Plasmas

Plasma Modeling with COMSOL Multiphysics

ADVENTURES IN TWO-DIMENSIONAL PARTICLE-IN-CELL SIMULATIONS OF ELECTRONEGATIVE DISCHARGES

The role of recycling in pulsed sputtering magnetrons

Control of Ion Energy Distributions on Plasma Electrodes

EFFECT OF PRESSURE AND ELECTRODE SEPARATION ON PLASMA UNIFORMITY IN DUAL FREQUENCY CAPACITIVELY COUPLED PLASMA TOOLS *

Electron Density and Ion Flux in Diffusion Chamber of Low Pressure RF Helicon Reactor

NARROW GAP ELECTRONEGATIVE CAPACITIVE DISCHARGES AND STOCHASTIC HEATING

Effect of Noble Gas. Plasma Processing Laboratory University of Houston. Acknowledgements: DoE Plasma Science Center and NSF

Modélisation de sources plasma froid magnétisé

Electron energy distribution function measurements and plasma parameters in inductively coupled argon plasma

SCALING OF PLASMA SOURCES FOR O 2 ( 1 ) GENERATION FOR CHEMICAL OXYGEN-IODINE LASERS

CONSEQUENCES OF RADIATION TRAPPING ON ELECTRON ENERGY DISTRIBUTIONS IN LOW PRESSURE INDUCTIVELY COUPLED Hg/Ar DISCHARGES*

Experimental Investigations of Magnetic Reconnection. J Egedal. MIT, PSFC, Cambridge, MA

Plasmas rf haute densité Pascal Chabert LPTP, Ecole Polytechnique

The Role of Secondary Electrons in Low Pressure RF Glow Discharge

Global modeling of HiPIMS systems: transition from homogeneous to self organized discharges

One dimensional hybrid Maxwell-Boltzmann model of shearth evolution

Plasma diagnostics of pulsed sputtering discharge

Confinement of toroidal non-neutral plasma in Proto-RT

Two-dimensional simulation of a miniaturized inductively coupled plasma reactor

Confinement of toroidal non-neutral plasma in Proto-RT

PIC-MCC/Fluid Hybrid Model for Low Pressure Capacitively Coupled O 2 Plasma

Cluster fusion in a high magnetic field

Retarding Field Analyzer for Ion Energy Distribution Measurement Through a Radio-Frequency or Pulsed Biased Sheath

INTRODUCTION TO THE HYBRID PLASMA EQUIPMENT MODEL

Physics of High Pressure Helicon Plasma and Effect of Wavenumber Spectrum

The Q Machine. 60 cm 198 cm Oven. Plasma. 6 cm 30 cm. 50 cm. Axial. Probe. PUMP End Plate Magnet Coil. Filament Cathode. Radial. Hot Plate.

The Plasma Simulation System Brochure.

Angular anisotropy of electron energy distributions in inductively coupled plasmas

Electron Energy Distributions in a Radiofrequency Plasma. Expanded by Permanent Magnets

Large Plasma Device (LAPD)

Robert A. Meger Richard F. Fernster Martin Lampe W. M. Manheimer NOTICE

Two-dimensional Particle-In-Cell model of the extraction region of the PEGASES ion-ion plasma source

Floating probe for electron temperature and ion density measurement applicable to processing plasmas

The plasma simulation system Brochure.

LOW-PRESSURE radio-frequency (RF) inductive-coupled

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

PIC-MCC/Fluid Hybrid Model for Low Pressure Capacitively Coupled O 2 Plasma

Measurement of electron energy distribution function in an argon/copper plasma for ionized physical vapor deposition

Beams and magnetized plasmas

PIC-MCC simulations for complex plasmas

OPTIMIZATION OF PLASMA UNIFORMITY USING HOLLOW-CATHODE STRUCTURE IN RF DISCHARGES*

Nanopantography: A method for parallel writing of etched and deposited nanopatterns

Repetition: Practical Aspects

CHAPTER 8. SUMMARY AND OUTLOOK 90 Under the operational conditions used in the present work the translation temperatures can be obtained from the Dopp

Measurement of EEDF and Distribution of Primary Electron in a Bucket Ion Source

Model for noncollisional heating in inductively coupled plasma processing sources

Physics and Modelling of a Negative Ion Source Prototype for the ITER Neutral Beam Injection

65 th GEC, October 22-26, 2012

Nonlinear effects in inductively coupled plasmas a

Characteristics and classification of plasmas

ANGULAR DEPENDENCE OF ELECTRON VELOCITY DISTRIBUTIONS IN LOW-PRESSURE INDUCTIVELY COUPLED PLASMAS 1

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

Chapter VI: Cold plasma generation

Modélisation particulaire du plasma magnétron impulsionnel haute puissance

NARROW GAP ELECTRONEGATIVE CAPACITIVE DISCHARGES AND STOCHASTIC HEATING

Feature-level Compensation & Control

Formation of High-b ECH Plasma and Inward Particle Diffusion in RT-1

MODELING OF AN ECR SOURCE FOR MATERIALS PROCESSING USING A TWO DIMENSIONAL HYBRID PLASMA EQUIPMENT MODEL. Ron L. Kinder and Mark J.

KINETIC DESCRIPTION OF MAGNETIZED TECHNOLOGICAL PLASMAS

Low Temperature Plasma Technology Laboratory

The Effect of Discharge Characteristics on Dielectric Barrier Discharges According to the Relative Permittivity

Effect of electron energy distribution function on power deposition and plasma density in an inductively coupled discharge at very low pressures

Study of DC Cylindrical Magnetron by Langmuir Probe

PRINCIPLES OF PLASMA DISCHARGES AND MATERIALS PROCESSING

Magnetic Field Configuration Dependence of Plasma Production and Parallel Transport in a Linear Plasma Device NUMBER )

Comparison of the B field dependency of plasma parameters of a weakly magnetized inductive and Helicon hydrogen discharge

MODELING OF SEASONING OF REACTORS: EFFECTS OF ION ENERGY DISTRIBUTIONS TO CHAMBER WALLS*

Low-field helicon discharges

vacuum analysis plasma diagnostics surface science gas analysis

Density Collapse in Improved Confinement Mode on Tohoku University Heliac

Terry Kammash and Myoung-Jae Lee Department of Nuclear Engineering The University of Michigan Ann Arbor, MI (313)

June 5, 2012: Did you miss it? No Problem. Next chance in Venus in Front of the Sun and there is more: plasma imaging!

Derivation of dynamo current drive in a closed current volume and stable current sustainment in the HIT SI experiment

DEPOSITION OF THIN TiO 2 FILMS BY DC MAGNETRON SPUTTERING METHOD

Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas )

Multi-fluid Simulation Models for Inductively Coupled Plasma Sources

Scattering in Cold- Cathode Discharges

Two-dimensional Fluid Simulation of an RF Capacitively Coupled Ar/H 2 Discharge

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

Ion Beam Cocktail Development and ECR Ion Source Plasma Physics Experiments at JYFL

Section 5: Thin Film Deposition part 1 : sputtering and evaporation. Jaeger Chapter 6. EE143 Ali Javey

Low Temperature Plasma Technology Laboratory

Transcription:

Electron Temperature Modification in Gas Discharge Plasma Valery Godyak University of Michigan and RF Plasma Consulting egodyak@comcast.net Workshop: Control of Distribution Functions in Low Temperature Plasmas at 64th Annual Gaseous Electronics Conference November 14 18, 2011; Salt Lake City, Utah This work was supported in part by the DOE OFES (Contract No DE-SC0001939) 1

Electron Temperature in Gas Discharge (Uniform electric field, E grad(n) = 0, Maxwellian EEDF and direct ionization) plasma parameters are in equilibrium with electric field, spatial and temporal locality Ionization balance (continuity and momentum eqs.) in a steady-state, self-sustained. bounded plasma defines z, resulting in: T e = T e (pλ), independently on P d and n. Electron energy balance, P d = (3/2)VT e nξdv, results in: Re(E pl ) = const(pλ), n ~ P d Here ξ is the characteristic frequency of electron energy loss, ξ = v e /λ ε ξ = ν en 2m/M + Σ2ν*ε*/3T e + z{2ε i /3T e + (4/3) + ⅓ [1+ ln(m/2πm)]} The fundamental relations for the electric field E and the total number of electrons/ions N p follows from the electron energy balance: E 2 = 3T e mξυ eff (1+ ω 2 /υ 2 eff) and N p = <n>v = 2P d /3ξT e Thus, for given P d and pλ, T e and n should be the same for all kinds of discharges These are basically true for non-maxwellian EEDF and with non-linear processes

Non-local electron kinetics In gas discharge plasmas, electrons are not in equilibrium with molecules and ions, T e >>T i,t g, they are not in equilibrium within their own ensemble, non-maxwellian, and when λ ε = v e /ξ > L, they are not in equilibrium with a non-uniform heating electric field E. The last is domain of electron non-local kinetics where plasma parameters are not local function of the field, grad(t e ) 0 and df e (ε+ev)/dr 0 plasma density (cm -3 ) 1 10 11 6 10 10 2 10 10 Nonlocal effects in a low pressure ICP 10 11 10 10 electron temperature (ev) 8 T e 4 V p 0 0 2 4 6 8 10 axial position (cm) plasma potential (V) eepf (ev -3/2 cm -3 ) 10 9 10 8 10 7 z = 4.0, 2.0,1.0, 0.5, 0.2 cm 0 5 10 15 20 total electron energy (ev) 3

EEDF in a non uniform E-field At ω << ω p, the external electromagnetic field is localized at the plasma boundary, S, δ << L, and in plasma bulk E n -1 when E grad(n) 0 Electric field non-uniformity typically occurs when ω << ω p and electromagnetic field is localized at the plasma boundary, S, δ < L, and in plasma when E grad(n) 0 Hot electrons generated in the zone of strong electric field produce ionization in the area of week field. In the presence of some separation mechanism preventing new-born electrons mixing with hot electrons or/and to penetrate the heating zone, the new-born electrons remain cold. This results in EEDF having two electron groups (hot and cold). This cold plasma has features of a non-self-sustained discharge Thus, the local electron heating together with some separation mechanism result in plasma cooling 4

Electron cooling in negative glow of dc glow discharge Solntsev et al, 8 th ICPIG, p. 86, Vienna, 1967, have measured ultra-cold electrons (T e = 0.04-0.3 ev, 1-2 orders of magnitude lower than that in the positive column. He, 0.6-4 Torr, I d = 0.6 8 ma. Haas et al, PSST. 7, 471, 1998, have demonstrated plasma electron cooling by injecting 100 ev electron beam into CCP. 5

Heating mode transition in CCP Ar CCP at13.65 MHz, L = 2 cm Godyak and Piejak, Phys. Rev. Lett. 65, 996, 1990 6

Transition to high plasma density (γ-mode) CCP, 13.56 MHz, He 0.3 Torr T e pressure dependence Godyak et al, Phys. Rev. Lett. 68, 49, 1992 7

EEDF and plasma parameter in ICP Center of Ar ICP, 2R = 20 cm, L = 10 cm 10 13 10 12 10 11 6.78 MHz, 50 W Three-temperature structure of EEDF is due to selective electron heating at the condition of anomalous skin effect. At high plasma density EEDF trends to Maxwellian distribution eepf (ev -3/2 cm -3 ) 10 10 10 13 10 9 10 8 300 mt 10 7 100 10 ε* ε i 1 0.3 10 6 0 10 20 30 40 50 electron energy (ev) plasma density (cm -3 ) 10 12 10 11 10 10 Godyak et al, PSST, 11, 525, 2002 6.78 MHz, 50 W 8 6 4 2 0 10-1 10 0 10 1 10 2 10 3 gas pressure (mt) 10 effective electron temperature (ev) 8

Temporal non-locality: Time resolved EEDF measurement ω >> z dn/dt = 0, but ω ξ dt e /dt 0 Toroidal ICP driven with ferrite inductor at 50-250 khz, Hg-Ar, 0.3 Torr current sensor PC 10 11 50 khz, 2 A discharge vessel toroid transformer Power amplifier reference probe Langmuir probe cathode Probe station Waveform synthesizer trigger DAQ board control and analysis software Front panel (GUI) eepf (ev -3 / 2 cm -3 ) lamp voltage (abs. value V) 10 10 10 9 10 8 80 60 40 20 0 1 2 3 4 5 6 7 8 energy (ev) 1.5 50 khz 2 8 1 2 A 4 8 5 µs t = 1 µs 0.5 electron temperature (ev) Alexandrovich et al, LS10, p. 283, Toulouse, France, 2004 0 0 0 5 10 15 20 time (µs) 9

Frequency dependence of EEDF in ICP with anomalous skin effect Selective electron heating. Collisionless heating occurs at v e /δ > ω ε t ½m(δω) 2 At P d = 12 W f(mhz) ε t (ev) T eff (ev) 3.4 0.65 5.4 6.8 2.5 4.2 13.56 9.0 2.9 At high plasma density, the difference in T eff disappears, since ν ee ~ nε -3/2 Godyak and Kolobov, Phys. Rev. Lett., 81, 369, 1998 10

Temporal nonlocality: Electron temperature variation in pulse discharge and Low frequency RF Discharge EEDF in afterglow stage of ICP with internal coil Evolution T e and n in a periodically pulsed ICP 10 12 Ar, 30 mt, 50 W; off cycle Ar, CW 100 W eepf (ev -3/2 cm -3 ) 10 11 10 10 10 9 t = 2.8 µs 3.6 4.4 6.8 9.2 12.4 18.8 T = 2 µs on T = 20 µs off electrron temperature (ev) 1 0.1 3 mt 10 p Te (mt) (ev) 3.0 6.5 10 4.2 30 3.1 100 2.1 300 1.5 30 300 100 0 5 10 15 electron energy (ev) 1 10 100 1000 time (µs) Godyak and Alexandrovich, XXVII ICPIG, vol. 1, p.221, Eindhoven, The Nederland, 2005 11

EEDF Modification with discharge current constriction Ne/Hg p =1 Torr D = 5 cm, d = 1,25 and 2.5 cm Distance from orifice Godyak et al, Phys. Rev. A 38, 2044, 1988 12

T e control with negatively biased greed (Kato et al, 1994) Plasma source and diffusion zones are separated with negatively biased mesh Ikada et all, Thin Solid Films 457, 55, 2004 13

Localized ECR array reactor with trapped ECR plasma Ar/SiF 4 at 10 mtorr with microcrystalline silicon deposition. Multicusp magnetic confinement of fast electrons. An order of magnitude T e reduction in diffusion zone. Local ESR electron heating with separation between hot and cold electrons provided by magnetic filters Measurement position Bulkin et al, Ecole Polytechnique. 2010, to be published 14

Global magnetic filter with trapped ICP drive at 5 MHz 10 118 124 32 N S electron temperature (ev) 8 6 4 2 low probe, r = 0, 20 mm, 100 W B = 0 B = 0 pumping 0 1 10 100 gas pressure (mtorr) T eh Magnetic field breaks non-locality in electron kinetics, leading to plasma stratification on hot (T eh ) and cold (T ec ) electron zones, T eh >T e0 >T ec B = 0 T e0 T ec Godyak, 63 GEC, Paris, France, 2010 15

EEPF and plasma parameters along magnetic filter Ecole Polytechnique B=0 B=0 Aanesland et al, 2011, to be published 16

Conclusions In gas discharge plasma at λ ε >> Λ and large de/dr, EEDF is not in local equilibrium with E-field, plasma parameters and the field distributions are decoupled and df(ε+eφ)/dr 0 Generation of excess of high energy electrons cools down the main body of electron population leading to formation two electron groups Formation of highly non-equilibrium EEDF with two-temperature structure (T e1 << T e2 ) requires both, strong E-field localization (to produce fast electrons) and some separation mechanism preventing low energy electron heating and mixing with hot electrons. Non-equilibrium discharges with strong localization (in space and/or in time) of the heating field and with electron separation feature seems is a viable way for creation of plasma with controllable EEDF. 17