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

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

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

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

CONTROL OF UNIFORMITY IN CAPACITIVELY COUPLED PLASMAS CONSIDERING EDGE EFFECTS*

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

EFFECT OF REACTOR GEOMETRY ON ION ENERGY DISTRIBUTIONS FOR PULSED PLASMA DOPING (P 2 LAD)*

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

DOE WEB SEMINAR,

65 th GEC, October 22-26, 2012

SPUTTER-WIND HEATING IN IONIZED METAL PVD+

INVESTIGATION of Si and SiO 2 ETCH MECHANISMS USING an INTEGRATED SURFACE KINETICS MODEL

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

MWP MODELING AND SIMULATION OF ELECTROMAGNETIC EFFECTS IN CAPACITIVE DISCHARGES

arxiv: v1 [physics.plasm-ph] 10 Nov 2014

SIMULATIONS OF ECR PROCESSING SYSTEMS SUSTAINED BY AZIMUTHAL MICROWAVE TE(0,n) MODES*

SIMULATION OF POROUS LOW-k DIELECTRIC SEALING BY COMBINED He AND NH 3 PLASMA TREATMENT *

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

Physique des plasmas radiofréquence Pascal Chabert

SCALING OF HOLLOW CATHODE MAGNETRONS FOR METAL DEPOSITION a)

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

Feature-level Compensation & Control

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

Model for noncollisional heating in inductively coupled plasma processing sources

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

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

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

Plasma atomic layer etching using conventional plasma equipment

MODELING AND SIMULATION OF LOW TEMPERATURE PLASMA DISCHARGES

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

Plasma Processing in the Microelectronics Industry. Bert Ellingboe Plasma Research Laboratory

INTRODUCTION TO THE HYBRID PLASMA EQUIPMENT MODEL

Effect of sputter heating in ionized metal physical vapor deposition reactors

Plasma Modeling with COMSOL Multiphysics

Plasma Technology September 15, 2005 A UC Discovery Project

Electron Current Extraction and Interaction of RF mdbd Arrays

Dual-RadioFrequency Capacitively-Coupled Plasma Reactors. Tomás Oliveira Fartaria nº58595

Angular anisotropy of electron energy distributions in inductively coupled plasmas

The Plasma Simulation System Brochure.

Control of Ion Energy Distributions on Plasma Electrodes

Wave propagation and power deposition in magnetically enhanced inductively coupled and helicon plasma sources

The plasma simulation system Brochure.

Plasma abatement of perfluorocompounds in inductively coupled plasma reactors

Argon metastable densities in radio frequency Ar, Ar/O 2 and Ar/CF 4 electrical discharges

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

Plasma Chemistry Study in an Inductively Coupled Dielectric Etcher

Exploration COMSOL in Modeling RLSA TM CVD Processes

Plasmas rf haute densité Pascal Chabert LPTP, Ecole Polytechnique

Plasma Technology. FLCC Workshop & Review September 13, 2006 FLCC

DEPOSITION AND COMPOSITION OF POLYMER FILMS IN FLUOROCARBON PLASMAS*

NONLINEAR ELECTROMAGNETICS MODEL OF AN ASYMMETRICALLY DRIVEN CAPACITIVE DISCHARGE

I. INTRODUCTION J. Vac. Sci. Technol. A 16 4, Jul/Aug /98/16 4 /2454/9/$ American Vacuum Society 2454

Feature Profile Evolution during Shallow Trench Isolation (STI) Etch in Chlorine-based Plasmas

Copyright 1996, by the author(s). All rights reserved.

Multi-fluid Simulation Models for Inductively Coupled Plasma Sources

Plasma Processing of Large Curved Surfaces for SRF Cavity Modification

Investigation of the DSMC Approach for Ion/neutral Species in Modeling Low Pressure Plasma Reactor

E SC 412 Nanotechnology: Materials, Infrastructure, and Safety Wook Jun Nam

Numerical Simulation: Effects of Gas Flow and Rf Current Direction on Plasma Uniformity in an ICP Dry Etcher

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

Simulation of a two-dimensional sheath over a flat insulator conductor interface on a radio-frequency biased electrode in a high-density plasma

Lee Chen, Merritt Funk, and Radha Sundararajan Tokyo Electron America, Austin, Texas 78741

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

A3D Hybrid Model of ahelicon Source +

Consequences of asymmetric pumping in low pressure plasma processing reactors: A three-dimensional modeling study

Spatial distributions of power and ion densities in RF excited remote plasma reactors

Combinatorial RF Magnetron Sputtering for Rapid Materials Discovery: Methodology and Applications

Plasma Chemistry and Kinetics in Low Pressure Discharges

PRINCIPLES OF PLASMA DISCHARGES AND MATERIALS PROCESSING

PIC-MCC simulations for complex plasmas

Application of Rarefied Flow & Plasma Simulation Software

The Plasma Simulation System Brochure.

A Working Electron Impact Cross Section Set for CHF 3. Mark J. Kushner a) and Da Zhang b) University of Illinois 1406 W. Green St Urbana, IL 61801

Extrel Application Note

Shapes of agglomerates in plasma etching reactors

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

PIC/MCC Simulation of Radio Frequency Hollow Cathode Discharge in Nitrogen

Electron Temperature Modification in Gas Discharge Plasma

The Gaseous Electronic Conference GEC reference cell as a benchmark for understanding microelectronics processing plasmas*

Influence of driving frequency on the metastable atoms and electron energy distribution function in a capacitively coupled argon discharge

Effects of fast atoms and energy-dependent secondary electron emission yields in PIC/ MCC simulations of capacitively coupled plasmas

MONTE CARLO SIMULATION OF RADIATION TRAPPING IN ELECTRODELESS LAMPS: A STUDY OF COLLISIONAL BROADENERS*

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

Spatially resolved mass spectrometric sampling of inductively coupled plasmas using a movable sampling orifice

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

ETCHING Chapter 10. Mask. Photoresist

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

MICRODISCHARGES AS SOURCES OF PHOTONS, RADICALS AND THRUST*

Simulation of a two-dimensional sheath over a flat wall with an insulatorõconductor interface exposed to a high density plasma

Electromagnetics in COMSOL Multiphysics is extended by add-on Modules

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

The Role of Secondary Electrons in Low Pressure RF Glow Discharge

Plasma Etching: Atomic Scale Surface Fidelity and 2D Materials

A time dependent propagator method for long mean free path transport of neutral particles in plasma processing reactors

4 Modeling of a capacitive RF discharge

Striations in electronegative capacitively coupled radio-frequency plasmas: effects of the pressure, voltage, and electrode gap

FINAL REPORT. DOE Grant DE-FG03-87ER13727

NARROW GAP ELECTRONEGATIVE CAPACITIVE DISCHARGES AND STOCHASTIC HEATING

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

PARTICLE MODELLING OF PARALLEL PLATE RADIO FREQUENCY DISCHARGE PLASMAS IN HYDROGEN: INFLUENCE OF PRESSURE

M o d e l l i n g o f D u a l F r e q u e n c y C a p a c i t i v e l y. C o u p l e d P l a s m a D e v i c e s. P a u l C. Boyle

Transcription:

EFFECT OF PRESSURE AND ELECTRODE SEPARATION ON PLASMA UNIFORMITY IN DUAL FREQUENCY CAPACITIVELY COUPLED PLASMA TOOLS * Yang Yang a) and Mark J. Kushner b) a) Department of Electrical and Computer Engineering Iowa State University, Ames, IA 50011, USA yangying@eecs.umich.edu b) Department of Electrical Engineering and Computer Science, Ann Arbor, MI 48109, USA mjkush@umich.edu http://uigelz.eecs.umich.edu June 2009 * Work supported by Semiconductor Research Corp., Applied Materials and Tokyo Electron Ltd. YY_MJK_ICOPS2009_01

AGENDA Optimization of multiple frequency plasma etching reactors Description of the model Scaling with: Pressure Electrode separation Concluding remarks YY_MJK_ICOPS2009_02

MULTI-FREQUENCY PLASMA ETCHING REACTORS State of the art plasma etching reactors use multiple frequencies to create the plasma and accelerate ions into the wafer. Voltage finds its way into the plasma propagating around electrodes (not through them). YY_MJK_ICOPS2009_03 Ref: S. Rauf, AMAT

WAVE EFFECTS CHALLENGE SCALING As wafer size and frequencies increase - and wavelength decreases, electrostatic applied voltage takes on wavelike effects. Plasma shortened wavelength: λ = λ ( ) 1/ 2 0 1+ s = min(half plasma thickness, skin depth), s = sheath thickness YY_MJK_ICOPS2009_04 Lieberman, et al PSST 11 (2002) A. Perret, APhL 83 (2003) http://mrsec.wisc.edu

AN EXAMPLE: ADJUSTABLE GAP CONTROL Adjusting the gap (electrode separation) of capacitively coupled plasmas (CCPs) enables customization of the radical fluxes. Enables different processes, such as mask opening and trench etching, to be separately optimized. V. Vahedi, M. Srinivasan, A. Bailey, Solid State Technology, 51, November, 2008. YY_MJK_ICOPS2009_05

COUPLED EFFECTS IN HIGH FREQUENCY CCPs Electromagnetic wave effects impact processing uniformity in high frequency CCPs. When coupled with changing gap and pressure, controlling the plasma uniformity could be more difficult. Results from a computational investigation of impacts of pressure and gap on plasma uniformity in dual frequency CCPs (DF-CCPs) will be discussed. YY_MJK_ICOPS2009_06

HYBRID PLASMA EQUIPMENT MODEL (HPEM) Electron Energy Transport Module Te,S,μ E, N Fluid Kinetics Module Fluid equations (continuity, momentum, energy) Maxwell Equations Plasma Chemistry Monte Carlo Module YY_MJK_ICOPS2009_07 Electron Energy Transport Module: Electron Monte Carlo Simulation provides EEDs of bulk electrons Separate MCS used for secondary, sheath accelerated electrons Fluid Kinetics Module: Heavy particle and electron continuity, momentum, energy Maxwell s Equations Plasma Chemistry Monte Carlo Module: IEADs onto wafer

METHODOLOGY OF THE MAXWELL SOLVER Full-wave Maxwell solvers are challenging due to coupling between electromagnetic (EM) and sheath forming electrostatic (ES) fields. EM fields are generated by rf sources and plasma currents ES fields originate from charges. We separately solve for EM and ES fields and sum the fields for plasma transport. E v = v E EM Φ Boundary conditions (BCs): EM field: Determined by rf sources. ES field: Determined by blocking capacitor (DC bias) or applied DC voltages. ES YY_MJK_ICOPS2009_08

REACTOR GEOMETRY 2D, cylindrically symmetric. Base conditions Ar/CF 4 =90/10, 400 sccm High frequency (HF) upper electrode: 150 MHz, 300 W Low frequency (LF) lower electrode: 10 MHz, 300 W Specify power, adjust voltage. Main species in Ar/CF 4 mixture Ar, Ar*, Ar + CF 4, CF 3, CF 2, CF, C 2 F 4, C 2 F 6, F, F 2 CF 3+, CF 2+, CF +, F + e, CF 3-, F - YY_MJK_ICOPS2009_09

Ar PLASMA IN SINGLE FREQUENCY CCP [e] 10 mtorr, Max = 9.6 x 10 9 cm -3 50 mtorr, Max = 4.3 x 10 10 cm -3 With increasing Ar pressure, electron density transitions from center high to edge high. Agrees with experimental trend, albeit in a different geometry. DF-CCP at higher frequency, with electronegative gas trends? V. N. Volynets, et al., J. Vac. Sci. Technol. A 26, 406, 2008. YY_MJK_ICOPS2009_10 80 mtorr, Max = 1.5 x 10 11 cm -3 Ar, 100 MHz/750 W from upper electrode.

EM EFFECTS: FIELD IN SHEATHS HF = 50 MHz, Max = 410 V/cm Ar/CF 4 =90/10, 50 mtorr, 400 sccm HF: 300 W, LF: 10 MHz/300 W HF = 150 MHz, Max = 355 V/cm LF = 10 MHz, Max = 750 V/cm Low frequency electrostatic edge effect. High Frequency Constructive interference of waves in center of reactor. YY_MJK_ICOPS2009_11

SCALING WITH PRESSURE IN DF-CCP 10 mtorr, Max = 2.5 x 10 10 cm -3 50 mtorr, Max = 1.1 x 10 11 cm -3 75 mtorr, Max = 1.2 x 10 11 cm -3 With increasing pressure: Concurrent increase in [e]. Shift in maximum of [e] towards the HF electrode and the center of the reactor. The shift is a result of Shorter energy relaxation distance. Combination of finite wavelength and skin effect. 150 mtorr, Max = 1.6 x 10 11 cm -3 Ar/CF 4 =90/10 400 sccm HF: 150 MHz/300 W LF: 10 MHz/300 W YY_MJK_ICOPS2009_12

ELECTRON ENERGY DISTRIBUTIONS (EEDs) 10 mtorr 150 mtorr d = distance to the upper electrode. EEDF as a function of height: 10 mtorr no change in bulk plasma with tail lifted in sheath. 150 mtorr Tails of EEDs lift as HF electrode is approached. Produce different spatial distribution of ionization sources. Ar/CF 4 =90/10 400 sccm YY_MJK_ICOPS2009_13 HF: 150 MHz/300 W LF: 10 MHz/300 W

ELECTRON IMPACT IONIZATION SOURCE (S e ) Axial Direction Radial Direction (In the HF Sheath) With increasing pressure: Axial direction: Energy relaxation distance decreases and so sheath heating is dissipated close to electrode transition to net attachment. Radial direction: As energy relaxation distance decreases, S e mirrors the constructively interfered HF field - more center peaked. Ar/CF 4 =90/10 400 sccm YY_MJK_ICOPS2009_14 HF: 150 MHz/300 W LF: 10 MHz/300 W

ION FLUX INCIDENT ON WAFER Flux of Ar + Flux of CF 3 + With increasing pressure, ionization source increases but moves further from wafer.. Ar + flux is depleted by charge exchange reactions while diffusing to wafer and is maximum at 25-50 mtorr. Ar/CF 4 =90/10 400 sccm YY_MJK_ICOPS2009_15 HF: 150 MHz/300 W LF: 10 MHz/300 W

TOTAL ION IEADs INCIDENT ON WAFER: Ar/CF 4 = 90/10 Center Edge 10 mtorr Center Edge 150 mtorr Center Edge IEADs are separately collected over center and edge of wafer. Bimodal to single peak transition with increasing pressure. 10 mtorr: uniform 50 mtorr: larger radial variation. Ar/CF 4 =90/10, 400 sccm HF: 150 MHz LF: 10 MHz/300 W YY_MJK_ICOPS2009_16

SCALING WITH PRESSURE: Ar/CF 4 =80/20 10 mtorr, Max = 2.5 x 10 10 cm -3 50 mtorr, Max = 4.8 x 10 10 cm -3 100 mtorr, Max = 4.5 x 10 10 cm -3 With increasing pressure: [e] decreases from 50 to 150 mtorr owing to increasing attachment losses. Maximum of [e] still shifts towards the HF electrode and the reactor center a less dramatic shift than Ar/CF 4 =90/10. Electrostatic component remains dominant due to lower conductivity. 150 mtorr, Max = 4.2 x 10 10 cm -3 Ar/CF 4 =80/20 400 sccm HF: 150 MHz/300 W LF: 10 MHz/300 W YY_MJK_ICOPS2009_17

ION FLUX INCIDENT ON WAFER: Ar/CF 4 =80/20 Flux of Ar + Flux of CF 3 + Compared with Ar/CF 4 = 90/10 More rapid depletion of Ar + flux by charge exchange. CF 3+ flux also maximizes at intermediate pressure consequence of more confined plasma. Ar/CF 4 =80/20 400 sccm YY_MJK_ICOPS2009_18 HF: 150 MHz/300 W LF: 10 MHz/300 W

TOTAL ION IEADs INCIDENT ON WAFER: Ar/CF 4 = 80/20 Center Edge 10 mtorr Center Edge 150 mtorr Center Edge At Ar/CF 4 =80/20 plasma is peaked near HF electrode edge, and largely uniform over the surface of wafer. Improved uniformity of IEADs at all pressures. Ar/CF 4 =80/20, 400 sccm HF: 150 MHz LF: 10 MHz/300 W YY_MJK_ICOPS2009_19

SCALING WITH GAP: Ar/CF 4 =90/10 Gap = 1.5 cm, Max = 3.4 x 10 10 cm -3 2.5 cm, Max = 1.1 x 10 11 cm -3 3.5 cm, Max = 1.5 x 10 11 cm -3 5.5 cm, Max = 1.6 x 10 11 cm -3 With increasing gap: [e] increases as diffusion length increases and loss decreases. Edge peaked [e] at gap = 1.5 cm, due to electrostatic edge effect. Maximum of [e] shifts towards the HF electrode. For gap > 2.5 cm, radial [e] profile is not sensitive to gap. Electrode spacing exceeds energy relaxation length and power deposition mechanism does not change. Ar/CF 4 =90/10 50 mtorr, 400 sccm HF: 150 MHz/300 W LF: 10 MHz/300 W YY_MJK_ICOPS2009_20

EEDs vs GAP Gap = 1.5 cm Gap = 5.5 cm Ar/CF 4 =90/10 50 mtorr, 400 sccm HF: 150 MHz/300 W LF: 10 MHz/300 W YY_MJK_ICOPS2009_21 2.5 cm: Little change across bulk plasma; tail in LF sheath lifted owing to HF wave penetration. 5.5 cm: Systematic tail enhancement towards the HF electrode larger separation between HF and LF waves, system functions more linearly.

ION FLUX INCIDENT ON WAFER Flux of Ar + Flux of CF 3 + 1.5 cm: edge peaked flux due to electrostatic edge effect. 2.5-5.5 cm: middle peaked flux due to electrostatic and wave coupling. 6.5 cm: center peaked flux ( with a middle peaked [e] ): edge effect reduced at larger gap. Ar/CF 4 =90/10 50 mtorr, 400 sccm YY_MJK_ICOPS2009_22 HF: 150 MHz/300 W LF: 10 MHz/300 W

TOTAL ION IEADs INCIDENT ON WAFER vs GAP Center 1.5 cm Center Edge Edge 5.5 cm Center Edge Narrow gap has large center-to-edge nonuniformity due to change in sheath width. Narrower sheath near edge produces broaded IEAD. Large gap enables more diffusive and uniform sheath properties and so more uniform IEADs. Ar/CF 4 =90/10, 50 mtorr, 400 sccm HF: 150 MHz/300W LF: 10 MHz/300 W YY_MJK_ICOPS2009_23

CONCLUDING REMARKS For DF-CCPs sustained in Ar/CF 4 =90/10 mixture with HF = 150 MHz: With increasing pressure, maximum of ionization source (S e ) shifts towards the HF electrode as energy relaxation distance decreases. S e mirrors EM field, which is center peaked from constructive interference and [e] profile transitions from edge high to center high. Increasing fraction of CF 4 to 20% results in more uniform ion fluxes and IEADs incident on wafer. Effects of gap size in Ar/CF 4 =90/10 mixture: Between 2.5 and 6.5 cm, [e] profile is not sensitive to gap size since larger than energy relaxation distance. Small gaps have more edge-to-center non-uniformity in IEADs due to strong edge effects. YY_MJK_ICOPS2009_24