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

Similar documents
FINAL REPORT. DOE Grant DE-FG03-87ER13727

DOE WEB SEMINAR,

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

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

A global (volume averaged) model of a chlorine discharge

Equilibrium model for two low-pressure electronegative plasmas connected by a double layer

Plasma Processing of Large Curved Surfaces for SRF Cavity Modification

K. Takechi a) and M. A. Lieberman Department of Electrical Engineering and Computer Sciences, University of California, Berkeley, California 94720

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

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

DPP06 Meeting of The American Physical Society. Production of negative ion plasmas using perfluoromethylcyclohexane (C 7 F 14 )

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

Plasma properties determined with induction loop probes in a planar inductively coupled plasma source

In situ electrical characterization of dielectric thin films directly exposed to plasma vacuum-ultraviolet radiation

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

Feature-level Compensation & Control

Experimental verification of Boltzmann equilibrium for negative ions in weakly collisional electronegative plasmas

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

Plasmas rf haute densité Pascal Chabert LPTP, Ecole Polytechnique

Plasma Chemistry and Kinetics in Low Pressure Discharges

Effect of negative ions on the characteristics of plasma in a cylindrical discharge

ETCHING Chapter 10. Mask. Photoresist

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

Volume Production of D - Negative Ions in Low-Pressure D 2 Plasmas - Negative Ion Densities versus Plasma Parameters -

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

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.

Effect of wall charging on an oxygen plasma created in a helicon diffusion reactor used for silica deposition

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

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

Recombination and Decay of Plasma Produced by Washer Stacked Plasma Gun inside a Curved Vacuum Chamber

The low-field density peak in helicon discharges

Anode double layer in magnetized radio frequency inductively coupled hydrogen plasma

Contents: 1) IEC and Helicon 2) What is HIIPER? 3) Analysis of Helicon 4) Coupling of the Helicon and the IEC 5) Conclusions 6) Acknowledgments

Sheaths: More complicated than you think a

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

Helicon Plasma Thruster Experiment Controlling Cross-Field Diffusion within a Magnetic Nozzle

Trajectories of Ions inside a Faraday Cage Located in a High Density Plasma Etcher

Effect of Pressure and Discharge Voltage on Plasma Parameters in Air seeded Arc-plasma

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

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

PRINCIPLES OF PLASMA DISCHARGES AND MATERIALS PROCESSING

Physique des plasmas radiofréquence Pascal Chabert

Plasma Chemistry Study in an Inductively Coupled Dielectric Etcher

Plasma Diagnosis for Microwave ECR Plasma Enhanced Sputtering Deposition of DLC Films

The chemistry of a CCl 2 F 2 radio frequency discharge

Plasma Kinetic Study of Silicon-Dioxide Removal with Fluorocompounds in a Plasma-Enhanced Chemical Vapor Deposition Chamber

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

Nonlinear Diffusion in Magnetized Discharges. Francis F. Chen. Electrical Engineering Department

Low Temperature Plasma Technology Laboratory

Experimental investigation of double layers in expanding plasmas

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

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

Characterization of low pressure plasma-dc glow discharges (Ar, SF 6 and SF 6 /He) for Si etching

Study of DC Cylindrical Magnetron by Langmuir Probe

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

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

An Improved Global Model for Electronegative Discharge and Ignition Conditions for Peripheral Plasma Connected to a Capacitive Discharge

Chapter 7. Plasma Basics

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

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

Optical plasma emission spectroscopy of etching plasmas used in Si-based semiconductor processing

Progress in the Study of Plasma Heating, Stability, and Confinement on HANBIT Mirror Device

Optical Pumping of Rubidium

Hiden EQP Applications

Study of Electronegativity in Inductively Coupled Radio-Frequency Plasma with Langmuir Probe

Improvement of MgO Characteristics Using RF-Plasma Treatment in AC Plasma Display Panel

MODELING AND SIMULATION OF LOW TEMPERATURE PLASMA DISCHARGES

Electrical Discharges Characterization of Planar Sputtering System

A Kinetic Theory of Planar Plasma Sheaths Surrounding Electron Emitting Surfaces

Electrical characteristics and plasma diagnostics of (Ar/O 2 ) gas mixture glow Discharge

Magnetoresistance Oscillations in Thin Inclined Magnetic Films

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

Breakdown behavior in radio-frequency argon discharges

INTRODUCTION TO THE HYBRID PLASMA EQUIPMENT MODEL

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

Characteristics of Neutral Beam Generated by a Low Angle Reflection and Its Etch Characteristics by Halogen-Based Gases

Gas utilization in remote plasma cleaning and stripping applications

Chapter 7 Plasma Basic

Modeling of Ion Energy Distribution Using Time-Series Neural Network

Enhanced Chamber Management and Fault Detection in Plasma Etch Processes via SEERS(Self Excited Electron Resonance Spectroscopy)

Determination of Electron Temperatures in Rare- Gases Plasma

Electrical Breakdown in Low-Pressure Nitrogen in Parallel Electric and Magnetic Fields

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

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

Langmuir Probes as a Diagnostic to Study Plasma Parameter Dependancies, and Ion Acoustic Wave Propogation

Dynamics of Drift and Flute Modes in Linear Cylindrical ECR Plasma

Plasma Diagnostics Introduction to Langmuir Probes

LECTURE 5 SUMMARY OF KEY IDEAS

Multi-fluid Simulation Models for Inductively Coupled Plasma Sources

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

Influence of reactor wall conditions on etch processes in inductively coupled fluorocarbon plasmas

Plasma parameter evolution in a periodically pulsed ICP

Multiple steady states in a radio frequency chlorine glow discharge

Characteristics of Positive Ions in the Sheath Region of Magnetized Collisional Electronegative Discharges

STRONG DOUBLE LAYER STRUCTURE IN THERMIONIC VACUUM ARC PLASMA *

MODELING PLASMA PROCESSING DISCHARGES

Extrel Application Note

Plasma Theory for Undergraduate Education

A comparison of emissive probe techniques for electric potential measurements in a Hall thruster plasma

Transcription:

Journal of the Korean Physical Society, Vol. 38, No. 3, March 001, pp. 59 63 Characterization of an Oxygen Plasma by Using a Langmuir Probe in an Inductively Coupled Plasma Jong-Sik Kim and Gon-Ho Kim Department of Physics, Hanyang University, Ansan 45-795 Tae-Hun Chung Department of Physics, Dong-A University, Pusan 604-714 Geun-Young Yeom Department of Materials Science and Engineering, Sungkyunkwan University, Suwon 440-746 Kwang-Ho Kwon Department of Electronic Engineering, Hanseo University, Seosan 356-80 (Received 1 April 000) Oxygen plasmas were investigated by using a Langmuir probe for an inductively coupled plasma with various rf powers, 100 400 W, and operating pressures, 0.1 100 mtorr. The probe current ratio α of the positive and negative currents (I+/[I e + I ]) increased with the generation of negative ions and had its maximum value in the pressure region of 40 70 mtorr. Also, the operating pressure to achieve the maximum α shifted from the low-pressure region to the highpressure region with increasing input power because enhanced ion loss through positive-negative ion recombination. I. INTRODUCTION Oxygen plasmas have been widely used in many processes for semiconductor manufacturing. In the etching process using CF4, for example, the additive oxygen enhances the liberation of fluorine atoms in the CF4 plasma, resulting in an increased etch rate for the silicon wafer. The oxygen plasma also plays an important role in the ashing and cleaning processes. When oxygen gases are used, negative ions change the plasma potential and the behaviours of the electrons and the ions near the substrate, so it may cause dislocations on the surface during film formation [1,]. In spite of the need to understand the scaling of the oxygen plasma constituents with the control parameters, investigations have not been progressed sufficiently since the parameter space is much enlarged from that of electropositive plasmas; also, the diagnostics is very difficult in the presence of negative ions. Unlike electropositive plasma with ions and electrons, an electronegative plasma is composed of three different charges, positive ions, electrons, and negative ions, which E-mail: jongsik@newton.hanyang.ac.kr were formed from the attachment reaction of electron to the neutral atoms or molecules [3,4]. Kouznetosov et al. considered the electronegative discharge plasma to be distributed over three regions in the chamber: the electronegative bulk region, the electropositive edge region, and the sheath region. Hence, the plasma in the electronegative bulk region is composed of negative ions, positive ions, and electrons, so the plasma quasineutrality condition satisfies the charge balance n i n e + n, where n is negative ion density, n i the positive ion density, and n e the electron density. In comparison, the plasma in the electropositive edge region is composed of positive ions and electrons. Finally, the sheath region has mostly ions [3,5]. Here are the important reactions in oxygen plasmas: The ionization reaction between an oxygen molecule and an energetic electron results in a positive ion (O + e O + + e). Negative ions are produced by the dissociation reactions of an oxygen molecule and the attachment reactions between an oxygen atom and an electron (O + e O + O). The typical loss mechanism in the electronegative oxygen plasma is the recombination reaction between positive and negative ions, O + + O O + O. Note that we considered only -59-

-60- Journal of the Korean Physical Society, Vol. 38, No. 3, March 001 these three reactions in this study, so it was assumed that an increase in the number of negative ions implies a decrease in the number of electrons in the bulk plasma to satisfy the charge balance condition. In recent analytical studies by Lichtenberg et al. [6], it was shown that the negative ion production increases with operating pressure and input power. They argued that ion-flux loss to the wall and positive-negative ion recombination are the dominant ion loss mechanisms in electronegative plasmas. They also predicted that the positive-negative ion recombination enhances the ion loss with increasing number of the negative ions, so the number of ions is reduced in plasmas. As the optogalvanic [7] with spectrometer method [8] for measuring negative ions shows that the Langmuir probe technique by Amemiya [] is relatively inexpensive and less complicated. Amemiya measured of the negative ion density from the electron energy distribution function which was obtained from the second derivate of the probe data taken below the plasma potential. This method, however, needs great care to be applied in an rf plasma because the plasma potential fluctuation caused by E- mode discharge in the inductively coupled plasma is serious [9,10]. In this study, a common electrical probe was employed to study the generation mechanism of the electronegative oxygen plasma in inductively coupled plasmas (ICP). The probe technique is described in details in the next section. Here, we present measurements of electronegative oxygen plasmas generated in inductively coupled plasmas and the measurements were carried out over a wide range of operating pressures and rf input powers. The transition of the dominant positive ion loss mechanism is also presented. This paper is organized as follows: Section II describes how to measure the negative ions by using a Langmuir probe in an electronegative oxygen plasma. In Section III, the experimental setup is presented. Section IV presents experimental results and a discussions of the generation of the electronegative oxygen plasma at various powers and operating pressures. Conclusions are presented in Section V. II. TECHNIQUE OF A LANGMUIR PROBE IN ELECTRONEGATIVE OXYGEN PLASMA When Langmuir probe data is taken in a plasma with negative ions [], the I-V trace is composed of three saturated currents: The saturated positive current at a deep negative part of the probe bias of a plane Langmuir probe is I i 0.6n i Se Te m i, (1) which is mostly the ion current, and the saturated negative current at a positive part of probe bias is the electron current Te I e n e Se, () πm e and the negative ion current T I = n Se, (3) πm where e is the electronic charge, S the probe area, and m i, m e and m the masses of the ion, the electron, and the negative ion, respectively. T e and T are the temperatures of the electrons, and the negative ions, respectively. Note that if the negative ions increase, then the negative ion current increases, and, on the other hand, the electron current decreases, resulting in the negative current (sum of Eq.() + Eq.(3)) decreasing due to the negative ion velocity being very small compared to the electron velocity. Furthermore, the current ratio of the positive current to the negative current increases with the production of negative ions. We introduce new parameter α defined as the current ratio of the positive current to the negative current. Assuming that T = T + = T e /10 [5], O + (µ=3) is the major positive ion, and O (µ=16) is the major negative ion (see the details in the previous section), so m i =m [5], and the current ratio α can be expressed as α = I + 1.5n +. (4) I e + I 40n e + 0.n As just discussed, α is inversely proportional to the generation of negative ions. Because those assumptions are not always correct in real situations, it is very difficult to obtain an accurate density for the negative ions from this relation. Considerable care is required to ensure an accurate measurement from the Langmuir probe in inductively coupled plasmas. The plasma potential fluctuation due to the E-mode discharge in the ICP impedes the process of determining the plasma parameters [9,10]. The speed of data acquisition (generally less than 1 Hz) is much slower than the speed of the plasma potential oscillation (13.56 MHz in this study), so the current-voltage signals are time-averaged data, as shown in Fig. 1. Probe data taken in a real DC plasma [shown in Fig. 1(a)] show only one knee of I-V, which clearly appears in the log plot in Fig. 1(c), so the plasma potential is clearly obtained. On the other hand, the probe data taken from an ICP, shown in Fig. 1(b), has a stairway shape with two inflection points [shown in Fig. 1(d)] which correspond to fluctuations of the plasma potential. According to Hershkowitz [11], a probe I-V curve like Fig. 1(b) may be due to several causes: (1) the local ionization effect near the probe, () electron or ion beams in plasmas, and (3) rf fluctuating plasmas. To clarify the probe data, we compared the potentials from the Langmuir probe with those from the emissive probe, which is well known tool for measuring the plasma potential accurately without

Characterization of an Oxygen Plasma by Using a Langmuir Probe in Jong-Sik Kim et al. -61- Fig.. Schematic diagram of an inductively coupled plasma (ICP) source and Data Acquisition System (DAS) for Langmuir probe diagnostics. measure the saturated positive current, we fit a straight line to the current data smaller than Vp1 and then took the positive ion current at Vp. Similarity, a straight line was fitted to the current data larger than Vp, and the current value at Vp was chosen as the total negative current. The electron temperature is obtained from the line fitted to the log data below than Vp1 [as shown in Fig. 1(d)] III. EXPERIMENTAL SETUP Fig. 1. Langmuir probe data for DC and ICP plasmas: (a) example of Langmuir probe data taken in a DC plasma, (b) example of Langmuir probe data in a ICP plasma, (c) semilog plot of the DC plasma data (arrow indicates the plasma potential), and (d) semilog plot of ICP plasma data (arrows 1 and indicate the low and the high plasma potential, respectively). being influenced by local ionization or beams. Consequently, the inflection points shown in Fig. 1(b), agreed with the minimum [Vp1 in Fig. 1(d)] and maximum (Vp in Fig. 1(d)) plasma potentials, respectively. We concluded that these fluctuation may be caused by the E-mode discharge of the ICP [1]. The plasma potentials Vp1 and Vp, are chosen from the intersection of the fitted lines in the log plots, as shown in Fig. 1(d). To Characterization of the electronegative oxygen plasma was carried out in the inductively couple plasma source as shown in Fig. Up to 400 W of rf power (13.56-MHz, Henry rf power supply) was delivered to a resonant tank circuit, which consisted of a rectangular-shaped threeturn antenna (300 400 mm size and 3.4-uH inductance) and external capacitors. The forward and the reflected powers were monitored by a directional coupler (Bird Wattmeter Model 43), and the reflected power was kept to less than 10oxygen-free copper (square cross-section 10 10 mm with 5-mm diameter hole at the center for the path of the cooling water) was located on the top of the tempered glass plate (4-mm thickness). A single-sided planar Langmuir probe (Ta, 6.-mm diameter) was located 0 cm from the window and at the center of a cylindrical chamber (stainless steel, 500-mm diameter and 600-mm length). The Langmuir probe data were collected by a data acquisition system (DAS in Fig. ) composed of an oscilloscope with memory function and a personal computer. A bipolar power supply (KEPCO 00) was used to bias the probe. The plasma chamber was evacuated by a turbo-molecular pump (300-l/s pumping speed) and was maintained at a base pressure as 10 6 Torr. The ion and the convectron gauges monitored the operating pressure in the chamber, and they were used for low-pressure and high-pressure operations, respectively. The operating pressure was varied widely

-6- Journal of the Korean Physical Society, Vol. 38, No. 3, March 001 Fig. 3. Current ratio α as functions of the operating pressure and the input power in oxygen (solid line) and nitrogen (dotted line) plasmas. The dashed line shows the variation of the maximum values of the current ratio with input power. from 0.1 to 100 mtorr. Oxygen and nitrogen gases were used to study electronegative and electropositive plasmas, respectively. IV. EXPERIMENTAL RESULTS AND DISCUSSION The current ratio α of the positive current to the negative current composed of electron and negative ion currents [as shown in Eq. (4)] was obtained from the Langmuir probe directly. Note that, as mentioned in Sec. II, an increase in α implies an increase in negative-ion generation. Figure 3 shows α as functions of the operating pressure and the input power. The solid line indicates the data taken from the oxygen plasmas, and the dotted line these from the nitrogen plasmas. As shown in Fig. 3, α for the oxygen plasma is much larger and has more variation than the almost constant α, 0.0 over the whole pressure range, for the nitrogen plasma. Note that Fig. 4. Distribution of ion density as functions of the operating pressure and the input power in oxygen plasmas. The dashed line shows the variation of the maximum values of the current ratio with input power. Fig. 5. Distribution of the negative ion density as functions of the operating pressure and the input power in oxygen plasmas. Data were taken from Chung s global model[13] with experimental plasma parameters. The dashed line shows the variation of the maximum values of the current ratio with input power. the nitrogen plasma is well known to be an electropositive plasma. The ratio α of the oxygen plasma varies widely with pressure and increases with input power. The data show that α increases drastically at operating pressures of 40 70 mtorr, which is approximately 0 times larger than other values. Interestingly, as indicated by the dashed line in Fig. 3, the operating pressure to archive the maximum α shifts to higher values with increasing of input power. Figure 4 shows the distribution of the ion density as functions of the operating pressure and the input power in the oxygen plasma. The ion density increases with input power and has peak values at operating pressures of 30 50 mtorr. The operating pressure producing the maximum ion density decreases with the input power, and that is indicated by the dashed line in the figure. To consider the negative ion generation as functions of the operating pressure and the input power in the oxygen plasma, we used the measured values to apply Chung s global model [13 15]. Following Eq. (3) in Ref. 13, the negative ion density can be scaled as n O R Lu BO + n e R rec d eff n O + u BO + d eff K iz, (5) K iz and K rec are the ionization and the recombination rate constants, respectively, u BO + the Bohm velocity of positive ions, R L rate of recombination and wall flux loss, respectively. The effective plasma length d eff = 1 RL Rh L +Lh R, where R is the reactor radius, L is the reactor length, are h L and h R are defined as the ratios of the normalized axial and the radial sheath edge densities the densities over the central region of the discharge, respectively. Figure 5 shows the variation of the profile of u BO + /d eff K iz with the operating pressure and

Characterization of an Oxygen Plasma by Using a Langmuir Probe in Jong-Sik Kim et al. -63- input powers, which is proportional to the generation of negative ions. Note that this plot provides just the variation of the negative ions generation with the operating pressure and the input power. Generation of negative ions increases with the power, and maximum generation occurs at operating pressures of 30 80 mtorr. Generation of negative ions increases in the low-pressure region and decreases in the high-pressure region, resulting in the peak value that occurs at a certain pressure and increases with the input power [indicated by the dotted line in Fig. 5]. As discussed in Fig. 4, the pressure producing the maximum positive ion generation decreases with the input power; however, the pressure of maximum negative ion generation increases with the input power. It seems that this behavior is related to the ion loss mechanism in electronegative plasmas. Generation of negative ions increases with the pressure, so the major ion loss mechanism transits from the flux-to-the-wall loss (ionflux-dominated region) to positive-negative ion recombination loss (recombination-dominated region) with pressure, as predicted in the model [13 15] (see Sec. I). Consequently, with increasing operating pressure, the negative ion generation also increases, and that enhances of the ion loss through the positive-negative ion recombination. Hence, the ion current is reduced and α decreases so that the peak of current ratio α in Fig. 3 takes place at specific pressures of 40 70 mtorr to transit the ion loss mechanism. V. CONCLUSIONS The characteristics of electronegative oxygen plasma was investigated by using a Langmuir probe. The current ratio α of the ion saturation current to the negative saturation current shows peak values at certain operating pressures, typically 40 70 mtorr. In this region, the ion loss is enhanced through the positive-negative ion recombination due to an increase in the number of negative ions. Probe data also show that the generation of negative ions increases with the operating pressure and the input power. ACKNOWLEDGMENTS This work was supported by the Korean Ministry of Education through the Brain Korea (BK1) Project and Research Foundation in the program year of 1998 (Project No. 1998-003-D00100). Also, the authors wish to acknowledge the financial support of the Research Fund of Hanyang University and the 000-00 Hanbit User Development Program. REFERENCES [1] Matsuoka, Y. Hoshi and M. Naoe, J. Appl. Phys. 63, 098 (1988). [] H. Amemiya, J. J. Appl. Phys. 7, 1966 (1988). [3] I. G. Kouznetsov, A. J. Lichtenberg and M. A. Lieberman, Plasma Source Sci. Technol. 5, 66 (1996). [4] A. J. Lichtenberg, V. Vahedi, M. A. Lieberman and T. Rognlien, J. Appl. Phys. 75, 339 (1994). [5] J. T. Gudmundsson and M. A. Lieberman, Plasma Source Sci. Technol. 7, 1 (1998). [6] A. J. Lichtenberg, M. A. Lieberman, I. G. Kouznetsov and T. H. Chung, Plasma Source Sci. Technol. 9, 45 (000). [7] B. Barbieri and N. Beverini, Rev. Mod. Phys. 6, 603 (1990). [8] L. J. Overzet, L. H. Beerman and J. T. Verdeyen, J. Appl. Phys. 66, 16 (1989). [9] T. Lho, N. Hershkowitz, G-H. Kim, J. Miller and W. Steer, IEEE conf. on Plasma Science-Abstract, 5A0, p. 33 (1998). [10] T. Lho, N. Hershkowitz, G-H. Kim, W. Steer and J. Miller, Plasma Source Sci. Technol. 9, 5 (000). [11] N. Hershkowitz, How to Langmuir Probes work, in Plasma Diagnostics, Vol 1. ed. by O. Auciello, D. L. Flamm (Academic Press, New York, 1990). [1] Michale A. Lieberman and Allan J. Lichtenberg, Principle of Plasma Discharges and Materials Processing (Wiley, New York, 1994), Chap. 1. [13] T. H. Chung, H. J. Yoon and D. C. Seo, J. Appl. Phys. 86, 3536 (1999). [14] T. H. Chung, J. Korean Phys. Soc. 34, 4 (1999). [15] H. J. Yoon and T. H. Chung, J. Korean Phys. Soc. 34, 9 (1999).