The Role of Secondary Electrons in Low Pressure RF Glow Discharge
|
|
- Bethanie McCarthy
- 6 years ago
- Views:
Transcription
1 WDS'05 Proceedings of Contributed Papers, Part II, , ISBN MATFYZPRESS The Role of Secondary Electrons in Low Pressure RF Glow Discharge O. Brzobohatý and D. Trunec Department of Physical Electronics, Faculty of Science, Masaryk University, Kotlářská 2, Brno, Czech Republic Abstract. We concerned on the role of secondary electrons in low pressure rf glow discharge. The Particle in Cell/Monte Carlo computer simulations (PIC/MC) was applied to simulate the discharge. The influence of secondary electron emission yield (SEY) on the plasma density was studied. Furthermore the relative density and the electron energy probability function of secondary particles, i.e., secondary electrons and the products of ionizing collisions of secondary electrons was determined. We found that the secondary particles play very important role in this type of discharge. The relative density of secondary electrons created on electrodes is 2-10%, 5-15% for the value of SEY = 0.1, 0.2 respectively. But due to the high electric intensity in electrode sheaths these electrons have relatively high energy and therefore they can ionize neutral atoms. As consequence of this the relative density of the secondary electrons (secondaries) and electrons created in ionizing collisions of the secondaries is more than 40% and the relative density of ions created in ionizing collisions of the secondaries is more than 20%, 40% for value of SEY = 0.1, 0.2, respectively. Introduction Low pressure radio frequency (rf) discharges are often used for sputtering, etching and plasma enhanced chemical vapor deposition (PECVD). It is well-known that the secondary electron emission on the electrodes affects the plasma density [Ruas, et al., 2002, Lieberman and Lichtenber, 1994]. In the present paper the role of the secondary electrons and the influence of the value of the secondary emission yield (SEY) on the plasma density are studied via computer simulation. Our simulation model was based on Particle In Cell/Monte Carlo (PIC/MC) method. The PIC/MC method was used extensively for the study of plasma and the method is well documented in [Birdsall et al.,1991] and [Vahedi and Surendra,1991]. Modified simulation program PDP1 [Verboncoeur et al., 1993] was used. Model description We studied low pressure rf capacitively coupled discharge. One dimensional electrostatic model was applied, i.e., the electric field was only in one direction perpendicular to electrodes. Our model system was composed of two planar parallel electrodes, the diameter of the electrode was chosen 10 cm, the distance between the electrodes was chosen 3 cm. One of the electrode was powered by rf voltage with amplitude U 0 = 150 V and frequency f = MHz, the second one was grounded (zero voltage). Argon (pressure 3 20 Pa) was used as working gas. For the background neutral gas we considered Maxwellian velocity distribution with temperature 300 K. Density of neutral gas was uniform and depended only on the gas pressure. So we have plasma which is composed from electrons and positive ions, they move in the electric field and collide with the neutral atoms of the working gas. Their movement can be written by the Lorentz equation. Only the binary collisions were considered, as can be seen in the Table 1. Table 1. Collision processes used in the simulation. The energy necessary to the ionization is about 15.7 ev. electrons + Ar elastic collision e + Ar e + Ar excitation e + Ar e + Ar (11.55 ev) ionization e + Ar 2e + Ar +, ε i = 15.7 ev Ar + + Ar scattering Ar + + Ar Ar + + Ar charge transfer Ar + + Ar Ar + Ar + 306
2 Interactions of charged particles with the electrodes were considered, e.g., absorption, secondary emission. The secondary emission was included through the secondary emission yield (SEY). Generally the SEY depends on the energy and impact angle of impacting electron or ion. In our case we studied mainly the role of secondary electrons in the discharge. The dependence of the SEY on impact energy and angle was not considered in our calculations. We considered the same value of the SEY for both electrons and ions. Our computer model is based on two methods, i.e., Particle In Cell method and Monte Carlo method, so a particle computer simulation was applied. The PIC method was used for calculation of motion of charged particles and the MC method was used for simulation random processes, e.g., collisions of charged particles. Both used methods are well-known in the field of computer simulation. In next two sections we briefly describe both methods. Particle In Cell method If we want to calculate the motion of N charged particles in electromagnetic field (in our case only electric field) then we have to solve self-consistently equations of motion with Poisson equation m i r i = q i E = q i φ, i = 1,...,N, φ = ρ(r 1,...,r N )/ε 0. Both equations were solved by finite-difference method. The equations of motion were solved by well-known leap-frog method : v n+1/2 i m i v n 1/2 i t = q i E n, r n+1 i ri n t = v n+1/2 i, i = 1,..., N, where v n+1/2 i is velocity, r n is position of ith particle, n denotes time step (t n+1 = t n + t). The Poisson equation has form: φ j 1 2φ j + φ j+1 ( x) 2 = ρ j ε 0, j = 1,...,M 1 where φ j is potential and ρ j is density on the jth grid point of mathematical mesh. Boundary conditions φ 0 = U 0 cos(2πft), φ M = 0 were used. We need to know charge density ρ j at the grid points (used in our finite-difference schema) for the solution of Poisson equation. This charge density was obtained from the positions of charged particles (r 1,...,r N ). The method of determination of electric charge density from the positions of particles gives the name of the whole method Particle In Cell. There are many ways how to do it. We used so called Cloud In Cell Particle In Cell method that linearly assigns the charge of the particles to the grid points. Note that for the spatial and the time step exist the conditions of precision. The spatial step is then given by the condition x l rfsh < 1, where l rfsh is diameter of plasma sheath. The spatial step x is about m. The time step is given by the condition ω pl t < 0.2, where ω pl is the plasma frequency. This condition is in our case fullfiled for t < s. Note that the plasma frequency depends on the plasma density and its value was then about 10 9 s 1 for our plasma density ( m 3. More detailed description of the PIC method can be found in [Birdsall et al., 1991]. The computation cycle was composed from four steps, i.e., 1. (r1 n,...,rn n ) ρ(rn 1,...,rN n ), 2. ρ(r1 n,...,rn N ) φn j E i n, 3. collissions of charged particles using Monte Carlo method and 4. E n i (r1 n+1,...,r n+1 N ), tn+1 = t n + t, jump to the step 1. Monte Carlo method The MC method is method which make use of (pseudo)random numbers with a given probability function. The MC method is very often used for simulation of random processes, e.g., collisions of charged particles with neutrals, secondary emission on the electrodes and generation of initial velocities. The MC method used in our model for simulation of collisions is based on the calculation of probability of collision in time (t, t+ t). Time element t is used in finite difference schema of equations of motion. In this method the probability of collision in (t, t + t) is given by P( t) = 1 exp( ν t). (1) 307
3 The collision frequency of ith type of collision is given generally by ν i (r, ε) = nσ(ε) i v r (ε), where n is the density of background neutral gas, σ is collision cross section and v r is the relative velocity of colliding particles, ε is the energy of the colliding particle. The null collision method was used and therefore the total collision frequency ν t = i ν i is constant. From this probability a number of particles that will collide in time (t, t + t) is determined by following equation N coll = P( t)n tot, where N tot is total number of particles electrons or ions. This method is more detailed discussed in [Nanbu et al., 2000]. This method determines when a collision occurs but says nothing about its type. The type of collision is given by i c i=1 ν i 1 ν t r < i c i=1 ν i ν t, (2) where i c is the type of collision that will occur and r is the random number uniformly distributed in (0, 1). The knowledge of differential cross section and total cross section is very important for the MC method because these quantities determine the probability (2) of the given type of collision and the solid angle for scattering. The collision data for argon can be found in [Lieberman and Lichtenberg, 1994] or [Vahedi and Surendra, 1994]. The SEY gives the probability of secondary emission after impact of electrons, ions, respectively. The condition for realization of the secondary emission is then given by r SEY, where r is again the random number uniformly distributed in (0, 1). Results of the simulations In our simulations we monitored the secondary electrons and products from the ionizing collisions of secondary electrons with neutrals (their trajectories and collisions) in accordance with model assumptions in the some way as other plasma charged particles. The influence of the value of the SEY on the plasma density was studied. The plasma density and the electron energy probability function were measured in the middle of the system. We chose two values of the SEY, i.e., 0.1 and 0.2 and we studied the influence of change of the value of the SEY on the plasma density. We studied three groups of electrons. First studied group was composed of secondary electrons created on the electrode secondaries. Second studied group was composed of secondaries and electrons created in ionizing collisions of secondaries SECI and third group of electrons was all electrons in system. Note that the SECI included electrons created in the ionizing collisions of SECI. It is known that in the sheaths (electrode regions) is the highest electric field in the system and therefore the secondary electrons and electrons created in sheaths gain a lot of energy. The energy that secondaries gained in the plasma sheath is proportional to the plasma potential (in our case about 70 V) and to phase of the rf field, in which the secondaries were created on the electrode. We calculated the electron energy probability function (EEPF) for secondaries to show the relatively high temperature of secondaries, further we calculated the EEPF for SECI and finally we calculated the EEPF for all electrons in our system. The EEPFs are shown in Fig. 1 (pressure 4 Pa). As can be seen the secondaries have Maxwellian-like distribution with relatively high temperature about 3 ev. Note that the EEPF is defined by f(e) = F(E)/ E, where F(E) is electron energy distribution function (EEDF) and E is( energy ) of electrons. The EEDF of the Maxwellian distribution has this form F(E) = 2 1 π (kt) E exp E 3/2 kt. So the EEPF of Maxwellian distribution is represented by a line in the semi-log plot. The temperature of electrons is then simply given by tangent of this line. The SECI have bi-maxwellian-like distribution as well as all electrons in our system, the lower temperature is about 0.5 ev and the higher is about 3 ev. This bi-maxwellian distribution was measured by Godyak et al., [1990] too, here the lower temperature was 0.34 ev and the higher temperature was 3.1 ev (the distance between the electrodes was 2 cm, the linear diameter was 14.2 cm). The high-temperature group of electrons is created in sheaths and the low-temperature group of electrons is created in ionizing collisions. Furthermore this bi-maxwellian effect is forced by Ramsauer effect in argon. The low-temperature group of electrons has mean energy close to the energy of the Ramsauer minimum (in the cross section of the elastic collision of electrons with neutrals) therefore they move almost collisionlessly. Due to their low energy they cannot overcome ambipolar potential barrier in the plasma body, so they cannot reach the oscillating plasma sheath interface, where stochastic heating take place, thence they cannot gain a lot of energy. Number of low-temperature electrons is 80 % for low pressure and the number decreases with pressure. On the other hand the high-temperature group of electrons effectively interacts with argon atoms (in 308
4 secondaries SECI all electrons f(e) [ev -3/2 m -3 ] E [ev] Figure 1. The electron energy probability functions of studied groups of electrons secondaries SECI all electrons kt H kt [ev] kt L Figure 2. The pressure dependence of temperatures of studied groups of electrons in the system. collisions) and boundaries (electrodes). They compensates their energy losses through stochastic heating on the oscillating plasma sheath boundaries. Unlike the low-temperature electrons the high-temperature electrons easily overcome the ambipolar potential barrier and collide more frequently with the axial plasma boundaries [Godyak et al., 1990]. The bi-maxwellian-like distribution is observed only at low pressures (up to 12 Pa). The Ohmic heating is more important for higher pressures and therefore the bi-maxwellian like distribution came to Mawellian like distribution, i.e., the low-temperature group of electrons are thermalized by Ohmic heating. The pressure dependence of temperatures of all groups of electrons in the system is shown in the Fig. 2. It can be seen that both temperatures merged with growing pressure. The influence of the value of the SEY on the plasma density was determined. The pressure dependence of the plasma density, calculated in the middle of the system, is shown in Fig. 3. If we increase the value of the SEY from 0.1 to 0.2 the calculated plasma density increases of about 30%, see Fig. 4. The relatively high errorbars are mainly given by number of particles in the computer cell in mathematic mesh. The experimental determination of the influence of the material of the electrodes (different values of the SEY) can be found in [Ruas et al., 2002]. But the values of the SEYs of electrodes material were 309
5 n [10 15 m -3 ] SEY = 0.1 SEY = 0.2 Figure 3. The influence of the value of SEY on the plasma density. 40 n rel n rel [%] Figure 4. The influence of the SEC on the plasma density in dependence on pressure. On the ordinate there is plotted the relative increase of plasma density when increasing the SEC from 0.1 to 0.2. The relatively increase of plasma density n rel = n 0.2 /n 0.1 1, where n 0.1, n 0.2 is plasma density in the middle of the system for the SEY = 0.1, 0.2, respectively. not determined. So the comparison with experiment could not be done. The density (in the middle of the system) of secondaries was relatively low, i.e., n secele /n total = 3 10 %, 5 16, for the value of the SEY = 0.1, 0.2 respectively. Newertheless they play, due to their relatively high temperature, very important role. They have enough energy to ionize neutral atoms and therefore the relative concentration of the SECI is more than 40%, as you can see in Fig. 5. Ions created in ionizing collisions of secondary electrons comprise more than 20%, 40% for SEY =0.1, 0.2 respectively, as can be seen in Fig. 6. So the change of the value of the SEY from 0.1 to 0.2 caused that the relative density (n sec ion /n tot ) of ions created in ionizing collisions of secondaries increase by about 20 %. This is caused by increase of secondaries by about 2 6 %. 310
6 35 30 SEY = 0.1 SEY = n sec ele /n [%] SECI secondaries Figure 5. The relative density (n rel = n sec ele /n, where n rel is relative density in the middle of the system, n sec ele is density of secondaries or SECI and n is total plasma density) versus pressure n sec ion /n [%] SEY = 0.1 SEY = 0.2 Figure 6. The relative concentration of ions created in the ionizing collisions of the SECI versus pressure (n sec. ion /n, where n sec. ion is density of ions created in the ionizing collisions of secondaries, n is total plasma density). Conclusion In the present contribution the role of secondary electrons was studied via computer simulation. The computer model was based on the Particle In Cell and Monte Carlo method. We shown that the secondary electrons are very important due to their relatively high temperature. More than 40% of all electrons in the system was created on electrodes and in ionizing collisions of secondaries. More than 25% of all ions are ions created in ionizing collisions of secondaries. The EEPFs for different groups of electrons, i.e., for secondaries, for secondaries and electrons created in ionizing collisions of secondaries and for all electrons in the system, were calculated. We shown that the secondaries have Maxwellian-like EEPF. The secondaries and their ionizing products have bi-maxwellian-like distribution (for pressure less than 10 Pa) as all electrons in the system. These interesting results could be important in plasma processing like plasma deposition or sputtering. 311
7 Acknowledgments. This research has been supported by the Czech Science Foundation under project No. 202/03/0827 and grant FRVŠ No. 2910/2005. References Birdsall, et al.: Plasma Physics Via Computer Simulation. Bristol: Adam Hilger,1991. Godyak V A and Piejak R B, Phys. Rev. Lett., 65, 996, Nanbu, K: IEEE Trans. on Plasma Science, 28, 971, Lieberman M. A. and Lichtenberg A. J. Principles of plasma Discharges and Material Processing, John Wiley & Sons, New York, Ruas R., Verdonck P., Mansano R. D. and Braithwaite N. S. J. Czech. J. Phys., 52, D465, Vahedi, V, Surendra, M., Computer Physics Communication, 87, 179, Verboncoeur, et al.,j. Comp. Physics, 104, 321,
PIC-MCC simulations for complex plasmas
GRADUATE SUMMER INSTITUTE "Complex Plasmas August 4, 008 PIC-MCC simulations for complex plasmas Irina Schweigert Institute of Theoretical and Applied Mechanics, SB RAS, Novosibirsk Outline GRADUATE SUMMER
More information4 Modeling of a capacitive RF discharge
4 Modeling of a capacitive discharge 4.1 PIC MCC model for capacitive discharge Capacitive radio frequency () discharges are very popular, both in laboratory research for the production of low-temperature
More informationPIC-MCC/Fluid Hybrid Model for Low Pressure Capacitively Coupled O 2 Plasma
PIC-MCC/Fluid Hybrid Model for Low Pressure Capacitively Coupled O 2 Plasma Kallol Bera a, Shahid Rauf a and Ken Collins a a Applied Materials, Inc. 974 E. Arques Ave., M/S 81517, Sunnyvale, CA 9485, USA
More informationPIC-MCC/Fluid Hybrid Model for Low Pressure Capacitively Coupled O 2 Plasma
PIC-MCC/Fluid Hybrid Model for Low Pressure Capacitively Coupled O 2 Plasma Kallol Bera a, Shahid Rauf a and Ken Collins a a Applied Materials, Inc. 974 E. Arques Ave., M/S 81517, Sunnyvale, CA 9485, USA
More informationMonte Carlo Collisions in Particle in Cell simulations
Monte Carlo Collisions in Particle in Cell simulations Konstantin Matyash, Ralf Schneider HGF-Junior research group COMAS : Study of effects on materials in contact with plasma, either with fusion or low-temperature
More informationElectron Density and Ion Flux in Diffusion Chamber of Low Pressure RF Helicon Reactor
WDS'06 Proceedings of Contributed Papers, Part II, 150 155, 2006. ISBN 80-86732-85-1 MATFYZPRESS Electron Density and Ion Flux in Diffusion Chamber of Low Pressure RF Helicon Reactor R. Šmíd Masaryk University,
More informationNARROW GAP ELECTRONEGATIVE CAPACITIVE DISCHARGES AND STOCHASTIC HEATING
NARRW GAP ELECTRNEGATIVE CAPACITIVE DISCHARGES AND STCHASTIC HEATING M.A. Lieberman, E. Kawamura, and A.J. Lichtenberg Department of Electrical Engineering and Computer Sciences University of California
More informationInfluence of driving frequency on the metastable atoms and electron energy distribution function in a capacitively coupled argon discharge
Influence of driving frequency on the metastable atoms and electron energy distribution function in a capacitively coupled argon discharge S. Sharma Institute for Plasma Research, Gandhinagar -382428,
More informationDOE WEB SEMINAR,
DOE WEB SEMINAR, 2013.03.29 Electron energy distribution function of the plasma in the presence of both capacitive field and inductive field : from electron heating to plasma processing control 1 mm PR
More informationFINAL REPORT. DOE Grant DE-FG03-87ER13727
FINAL REPORT DOE Grant DE-FG03-87ER13727 Dynamics of Electronegative Plasmas for Materials Processing Allan J. Lichtenberg and Michael A. Lieberman Department of Electrical Engineering and Computer Sciences
More informationPIC/MCC Simulation of Radio Frequency Hollow Cathode Discharge in Nitrogen
PIC/MCC Simulation of Radio Frequency Hollow Cathode Discharge in Nitrogen HAN Qing ( ), WANG Jing ( ), ZHANG Lianzhu ( ) College of Physics Science and Information Engineering, Hebei Normal University,
More informationSolution of Time-dependent Boltzmann Equation
WDS'5 Proceedings of Contributed Papers, Part III, 613 619, 5. ISBN 8-8673-59- MATFYZPRESS Solution of Time-dependent Boltzmann Equation Z. Bonaventura, D. Trunec, and D. Nečas Masaryk University, Faculty
More informationP. Diomede, D. J. Economou and V. M. Donnelly Plasma Processing Laboratory, University of Houston
P. Diomede, D. J. Economou and V. M. Donnelly Plasma Processing Laboratory, University of Houston 1 Outline Introduction PIC-MCC simulation of tailored bias on boundary electrode Semi-analytic model Comparison
More informationarxiv: v1 [physics.plasm-ph] 10 Nov 2014
arxiv:1411.2464v1 [physics.plasm-ph] 10 Nov 2014 Effects of fast atoms and energy-dependent secondary electron emission yields in PIC/MCC simulations of capacitively coupled plasmas A. Derzsi 1, I. Korolov
More informationTwo-Dimensional Particle-in-Cell Simulation of a Micro RF Ion Thruster
Two-Dimensional Particle-in-Cell Simulation of a Micro RF Ion Thruster IEPC--7 Presented at the nd International Electric Propulsion Conference, Wiesbaden Germany September 5, Yoshinori Takao, Koji Eriguchi,
More informationHong Young Chang Department of Physics, Korea Advanced Institute of Science and Technology (KAIST), Republic of Korea
Hong Young Chang Department of Physics, Korea Advanced Institute of Science and Technology (KAIST), Republic of Korea Index 1. Introduction 2. Some plasma sources 3. Related issues 4. Summary -2 Why is
More informationElectron Transport Behavior in a Mirror Magnetic Field and a Non-uniform Electric Field
Commun. Theor. Phys. (Beijing, China) 35 (2001) pp. 207 212 c International Academic Publishers Vol. 35, No. 2, February 15, 2001 Electron Transport Behavior in a Mirror Magnetic Field and a Non-uniform
More informationSolution of time-dependent Boltzmann equation for electrons in non-thermal plasma
Solution of time-dependent Boltzmann equation for electrons in non-thermal plasma Z. Bonaventura, D. Trunec Department of Physical Electronics Faculty of Science Masaryk University Kotlářská 2, Brno, CZ-61137,
More informationCollisionless electron heating by capacitive radio-frequency plasma sheaths 2 and Lieberman[2, 3, 4], where the electrons moving towards the sheath ar
Collisionless electron heating by capacitive radio-frequency plasma sheaths G. Gozadinosyx, D. Vendery, M.M. Turnery and M.A. Liebermanz yplasma Research Laboratory, School of Physical Sciences Dublin
More informationPRINCIPLES OF PLASMA DISCHARGES AND MATERIALS PROCESSING
PRINCIPLES OF PLASMA DISCHARGES AND MATERIALS PROCESSING Second Edition MICHAEL A. LIEBERMAN ALLAN J, LICHTENBERG WILEY- INTERSCIENCE A JOHN WILEY & SONS, INC PUBLICATION CONTENTS PREFACE xrrii PREFACE
More informationEffects of fast atoms and energy-dependent secondary electron emission yields in PIC/ MCC simulations of capacitively coupled plasmas
(14pp) Plasma Sources Science and Technology doi:10.1088/0963-0252/24/3/034002 Effects of fast atoms and energy-dependent secondary electron emission yields in PIC/ MCC simulations of capacitively coupled
More informationScaling laws verification for capacitive rf-discharge Ar plasma using particle-in-cell simulations
Scaling laws verification for capacitive rf-discharge Ar plasma using particle-in-cell simulations T H. Chunga) and H. S. Yoon Department of Physics, Dong-A University, Pusan 64-714, Korea J. K. Lee Department
More informationP. Diomede, D. J. Economou and V. M. Donnelly Plasma Processing Laboratory, University of Houston
P. Diomede, D. J. Economou and V. M. Donnelly Plasma Processing Laboratory, University of Houston Acknowledgements: DoE Plasma Science Center, NSF Presented at the 57 th AVS Conference, Albuquerque, NM
More informationContents Motivation Particle In Cell Method Projects Plasma and Ion Beam Simulations
PIC Method for Numerical Simulation Ninad Joshi NNP Group 1 Contents Motivation Particle In Cell Method Projects Plasma and Ion Beam Simulations Motivation 3 Particle simulation Ion beams and Plasmas Accelerators
More informationSimulation of a two-dimensional sheath over a flat insulator conductor interface on a radio-frequency biased electrode in a high-density plasma
JOURNAL OF APPLIED PHYSICS VOLUME 95, NUMBER 7 1 APRIL 2004 Simulation of a two-dimensional sheath over a flat insulator conductor interface on a radio-frequency biased electrode in a high-density plasma
More informationThe Computational Simulation of the Positive Ion Propagation to Uneven Substrates
WDS' Proceedings of Contributed Papers, Part II, 5 9,. ISBN 978-8-778-85-9 MATFYZPRESS The Computational Simulation of the Positive Ion Propagation to Uneven Substrates V. Hrubý and R. Hrach Charles University,
More informationPARTICLE MODELLING OF PARALLEL PLATE RADIO FREQUENCY DISCHARGE PLASMAS IN HYDROGEN: INFLUENCE OF PRESSURE
Journal of Optoelectronics and Advanced Materials Vol. 7, No. 5, October 2005, p. 2363-2369 PARTICLE MODELLING OF PARALLEL PLATE RADIO FREQUENCY DISCHARGE PLASMAS IN HYDROGEN: INFLUENCE OF PRESSURE P.
More informationOne dimensional hybrid Maxwell-Boltzmann model of shearth evolution
Technical collection One dimensional hybrid Maxwell-Boltzmann model of shearth evolution 27 - Conferences publications P. Sarrailh L. Garrigues G. J. M. Hagelaar J. P. Boeuf G. Sandolache S. Rowe B. Jusselin
More informationFrequency variation under constant power conditions in hydrogen radio frequency discharges
JOURNAL OF APPLIED PHYSICS VOLUME 89, NUMBER 3 1 FEBRUARY 2001 Frequency variation under constant power conditions in hydrogen radio frequency discharges E. Amanatides and D. Mataras a) Plasma Technology
More informationMWP MODELING AND SIMULATION OF ELECTROMAGNETIC EFFECTS IN CAPACITIVE DISCHARGES
MWP 1.9 MODELING AND SIMULATION OF ELECTROMAGNETIC EFFECTS IN CAPACITIVE DISCHARGES Insook Lee, D.B. Graves, and M.A. Lieberman University of California Berkeley, CA 9472 LiebermanGEC7 1 STANDING WAVES
More informationarxiv: v1 [physics.plasm-ph] 16 Apr 2018
Consistent simulation of capacitive radio-frequency discharges and external matching networks Frederik Schmidt Institute of Theoretical Electrical Engineering, arxiv:1804.05638v1 [physics.plasm-ph] 16
More informationSimulation of a two-dimensional sheath over a flat wall with an insulatorõconductor interface exposed to a high density plasma
JOURNAL OF APPLIED PHYSICS VOLUME 94, NUMBER 5 1 SEPTEMBER 2003 Simulation of a two-dimensional sheath over a flat wall with an insulatorõconductor interface exposed to a high density plasma Doosik Kim
More informationADVENTURES IN TWO-DIMENSIONAL PARTICLE-IN-CELL SIMULATIONS OF ELECTRONEGATIVE DISCHARGES
ADVENTURES IN TWO-DIMENSIONAL PARTICLE-IN-CELL SIMULATIONS OF ELECTRONEGATIVE DISCHARGES PART 1: DOUBLE LAYERS IN A TWO REGION DISCHARGE E. Kawamura, A.J. Lichtenberg, M.A. Lieberman and J.P. Verboncoeur
More informationDepartment of Aerospace Engineering and Engineering Mechanics, The University of Texas at Austin, Austin, Texas 78712, USA
1 MAGNETIZED DIRECT CURRENT MICRODISCHARGE, I: EFFECT OF THE GAS PRESSURE Dmitry Levko and Laxminarayan L. Raja Department of Aerospace Engineering and Engineering Mechanics, The University of Texas at
More information65 th GEC, October 22-26, 2012
65 th GEC, October 22-26, 2012 2D Fluid/Analytical Simulation of Multi-Frequency Capacitively-Coupled Plasma Reactors (CCPs) E. Kawamura, M.A. Lieberman, D.B. Graves and A.J. Lichtenberg A fast 2D hybrid
More informationPlasma Chemistry and Kinetics in Low Pressure Discharges
Plasma Chemistry and Kinetics in Low Pressure Discharges Jón Tómas Guðmundsson Science Institute, University of Iceland, Iceland tumi@hi.is 12o. Encontro Brasileiro de Física de Plasmas Brasilia, Brazil
More informationContents: 1) IEC and Helicon 2) What is HIIPER? 3) Analysis of Helicon 4) Coupling of the Helicon and the IEC 5) Conclusions 6) Acknowledgments
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 IEC:! IEC at UIUC modified into a space thruster.! IEC has
More informationStudy of DC Cylindrical Magnetron by Langmuir Probe
WDS'2 Proceedings of Contributed Papers, Part II, 76 8, 22. ISBN 978-737825 MATFYZPRESS Study of DC Cylindrical Magnetron by Langmuir Probe A. Kolpaková, P. Kudrna, and M. Tichý Charles University Prague,
More informationDual-RadioFrequency Capacitively-Coupled Plasma Reactors. Tomás Oliveira Fartaria nº58595
Dual-RadioFrequency Capacitively-Coupled Plasma Reactors Tomás Oliveira Fartaria nº58595 Index Capacitive Reactors Dual Frequency Capacitively-Coupled reactors o Apparatus for improved etching uniformity
More informationFloating probe for electron temperature and ion density measurement applicable to processing plasmas
JOURNAL OF APPLIED PHYSICS 101, 033305 2007 Floating probe for electron temperature and ion density measurement applicable to processing plasmas Min-Hyong Lee, Sung-Ho Jang, and Chin-Wook Chung a Department
More informationElectron Series Resonant Discharges: Part II: Simulations of Initiation
Electron Series Resonant Discharges: Part II: Simulations of Initiation K. J. Bowers 1 W. D. Qiu 2 C. K. Birdsall 3 Plasma Theory and Simulation Group Electrical Engineering and Computer Science Department
More informationPlasma parameter evolution in a periodically pulsed ICP
Plasma parameter evolution in a periodically pulsed ICP V. Godyak and B. Alexandrovich OSRAM SYLVANIA, 71 Cherry Hill Drive, Beverly, MA 01915, USA The electron energy probability function (EEPF) has been
More informationEquilibrium model for two low-pressure electronegative plasmas connected by a double layer
PHYSICS OF PLASMAS 13, 093504 2006 Equilibrium model for two low-pressure electronegative plasmas connected by a double layer P. Chabert, a N. Plihon, C. S. Corr, and J.-L. Raimbault Laboratoire de Physique
More informationElectron Energy Distributions in a Radiofrequency Plasma. Expanded by Permanent Magnets
J. Plasma Fusion Res. SERIES, Vol. 9 (21) Electron Energy Distributions in a Radiofrequency Plasma Expanded by Permanent Magnets Tomoyo SASAKI, Kazunori TAKAHASHI, and Tamiya FUJIWARA Department of Electrical
More informationNARROW GAP ELECTRONEGATIVE CAPACITIVE DISCHARGES AND STOCHASTIC HEATING
NARRW GAP ELECTRNEGATIVE CAPACITIVE DISCHARGES AND STCHASTIC HEATING M.A. Lieberman Deartment of Electrical Engineering and Comuter Sciences University of California Berkeley, CA 9472 Collaborators: E.
More informationElectron heating in capacitively coupled radio frequency discharges
Electron heating in capacitively coupled radio frequency discharges Dissertation zur Erlangung des Grades eines Doktors der Naturwissenschaften in der Fakultät für Physik und Astronomie der Ruhr-Universität
More informationKinetic interpretation of resonance phenomena in low pressure capacitively coupled radio frequency plasmas
PHYSICS OF PLASMAS 23, 063514 (2016) Kinetic interpretation of resonance phenomena in low pressure capacitively coupled radio frequency plasmas Sebastian Wilczek, 1 Jan Trieschmann, 1 Denis Eremin, 1 Ralf
More informationAngular anisotropy of electron energy distributions in inductively coupled plasmas
JOURNAL OF APPLIED PHYSICS VOLUME 94, NUMBER 9 1 NOVEMBER 2003 Angular anisotropy of electron energy distributions in inductively coupled plasmas Alex V. Vasenkov a) and Mark J. Kushner b) Department of
More informationStriations in electronegative capacitively coupled radio-frequency plasmas: effects of the pressure, voltage, and electrode gap
Striations in electronegative capacitively coupled radio-frequency plasmas: effects of the pressure, voltage, and electrode gap Yong-Xin Liu 1, Ihor Korolov 2, Edmund Schüngel 3, You-Nian Wang 1, Zoltán
More informationSimulation of Low Pressure Plasma Processing Reactors: Kinetics of Electrons and Neutrals
Simulation of Low Pressure Plasma Processing Reactors: Kinetics of Electrons and Neutrals R. R. Arslanbekov and V. I. Kolobov CFD Research Corporation, Huntsville, AL, USA Abstract. In this paper, we illustrate
More informationFeature-level Compensation & Control
Feature-level Compensation & Control 2 Plasma Eray Aydil, UCSB, Mike Lieberman, UCB and David Graves UCB Workshop November 19, 2003 Berkeley, CA 3 Feature Profile Evolution Simulation Eray S. Aydil University
More informationEffect of Gas Flow Rate and Gas Composition in Ar/CH 4 Inductively Coupled Plasmas
COMSOL CONFERENCE BOSTON 2011 Effect of Gas Flow Rate and Gas Composition in Ar/CH 4 Inductively Coupled Plasmas Keisoku Engineering System Co., Ltd., JAPAN Dr. Lizhu Tong October 14, 2011 1 Contents 1.
More informationMULTIPACTOR ON A DIELECTRIC SURFACE WITH LONGITUDINAL RF ELECTRIC FIELD ACTION
Progress In Electromagnetics Research Letters, Vol. 24, 177 185, 211 MULTIPACTOR ON A DIELECTRIC SURFACE WITH LONGITUDINAL RF ELECTRIC FIELD ACTION F. Zhu *, Z. Zhang, J. Luo, and S. Dai Key Laboratory
More informationLecture 6 Plasmas. Chapters 10 &16 Wolf and Tauber. ECE611 / CHE611 Electronic Materials Processing Fall John Labram 1/68
Lecture 6 Plasmas Chapters 10 &16 Wolf and Tauber 1/68 Announcements Homework: Homework will be returned to you on Thursday (12 th October). Solutions will be also posted online on Thursday (12 th October)
More informationElectron Temperature Modification in Gas Discharge Plasma
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
More informationTwo-dimensional Fluid Simulation of an RF Capacitively Coupled Ar/H 2 Discharge
Two-dimensional Fluid Simulation of an RF Capacitively Coupled Ar/H 2 Discharge Lizhu Tong Keisoku Engineering System Co., Ltd., Japan September 18, 2014 Keisoku Engineering System Co., Ltd., 1-9-5 Uchikanda,
More informationPlasma Technology September 15, 2005 A UC Discovery Project
1 Feature-level Compensation & Control Plasma Technology September 15, 2005 A UC Discovery Project 9/15/05 - Plasma Technology 2 Plasma Technology Professors Jane P. Chang (UCLA), Michael A. Lieberman,
More informationAssessment of fluctuation-induced and wall-induced anomalous electron transport in HET
Assessment of fluctuation-induced and wall-induced anomalous electron transport in HET IEPC-2015-418 Presented at Joint Conference of 30th International Symposium on Space Technology and Science 34th International
More informationLee Chen, Merritt Funk, and Radha Sundararajan Tokyo Electron America, Austin, Texas 78741
Measurement of electron temperatures and electron energy distribution functions in dual frequency capacitively coupled CF 4 /O 2 plasmas using trace rare gases optical emission spectroscopy Zhiying Chen,
More informationParticle-In-Cell Simulations of a Current-Free Double Layer
Particle-In-Cell Simulations of a Current-Free Double Layer S. D. Baalrud 1, T. Lafleur, C. Charles and R. W. Boswell American Physical Society Division of Plasma Physics Meeting November 10, 2010 1Present
More informationInplane Microwave Plasma
Inplane Microwave Plasma Introduction Wave heated discharges may be very simple where a plane wave is guided into a reactor using a waveguide or very complicated as in the case with ECR (electron cyclotron
More informationEFFECT OF PRESSURE AND ELECTRODE SEPARATION ON PLASMA UNIFORMITY IN DUAL FREQUENCY CAPACITIVELY COUPLED PLASMA TOOLS *
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
More informationApplication of Rarefied Flow & Plasma Simulation Software
2016/5/18 Application of Rarefied Flow & Plasma Simulation Software Yokohama City in Japan Profile of Wave Front Co., Ltd. Name : Wave Front Co., Ltd. Incorporation : March 1990 Head Office : Yokohama
More informationChapter VI: Cold plasma generation
Introduction This photo shows the electrical discharge inside a highpressure mercury vapor lamp (Philips HO 50) just after ignition (Hg + Ar) Chapter VI: Cold plasma generation Anode Positive column Cathode
More informationMeasurement of electron energy distribution function in an argon/copper plasma for ionized physical vapor deposition
Measurement of electron energy distribution function in an argon/copper plasma for ionized physical vapor deposition Z. C. Lu, J. E. Foster, T. G. Snodgrass, J. H. Booske, and A. E. Wendt a) Engineering
More informationCONSEQUENCES OF RADIATION TRAPPING ON ELECTRON ENERGY DISTRIBUTIONS IN LOW PRESSURE INDUCTIVELY COUPLED Hg/Ar DISCHARGES*
CONSEQUENCES OF RADIATION TRAPPING ON ELECTRON ENERGY DISTRIBUTIONS IN LOW PRESSURE INDUCTIVELY COUPLED Hg/Ar DISCHARGES* Kapil Rajaraman**, Alex Vasenkov*** and Mark J. Kushner*** **Department of Physics
More informationScattering in Cold- Cathode Discharges
Simulating Electron Scattering in Cold- Cathode Discharges Alexander Khrabrov, Igor Kaganovich*, Vladimir I. Demidov**, George Petrov*** *Princeton Plasma Physics Laboratory ** Wright-Patterson Air Force
More informationDriving frequency effects on the mode transition in capacitively coupled argon discharges
Driving frequency effects on the mode transition in capacitively coupled argon discharges Liu Xiang-Mei( ), Song Yuan-Hong( ), and Wang You-Nian( ) School of Physics and Optoelectronic Technology, Dalian
More informationPlasma Technology. FLCC Workshop & Review September 13, 2006 FLCC
1 Plasma Technology Professors Jane P. Chang (UCLA), Michael A. Lieberman, David B. Graves (UCB) and Allan J. Lichtenberg, John P. Verboncoeur, Alan Wu, Emi Kawamura, Chengche Hsu, Joe Vegh, Insook Lee
More informationTHE INFLUENCE OF EXTERNAL MAGNETIC FIELD ON THE RADIATION EMITTED BY NEGATIVE GLOW OF A DC GLOW DISCHARGE
PLASMA PHYSICS THE INFLUENCE OF EXTERNAL MAGNETIC FIELD ON THE RADIATION EMITTED BY NEGATIVE GLOW OF A DC GLOW DISCHARGE M. TOMA, I. A. RUSU, D. O. DOROHOI Plasma Physics Department, A. I. Cuza University,
More informationBeams and magnetized plasmas
Beams and magnetized plasmas 1 Jean-Pierre BOEUF LAboratoire PLAsma et Conversion d Energie LAPLACE/ CNRS, Université Paul SABATIER, TOULOUSE Beams and magnetized plasmas 2 Outline Ion acceleration and
More informationTwo-dimensional simulation of a miniaturized inductively coupled plasma reactor
JOURNAL OF APPLIED PHYSICS VOLUME 95, NUMBER 5 1 MARCH 2004 Two-dimensional simulation of a miniaturized inductively coupled plasma reactor Sang Ki Nam and Demetre J. Economou a) Plasma Processing Laboratory,
More informationNumerical Study of Power Deposition, Transport and Acceleration Phenomena in Helicon Plasma Thrusters
Numerical Study of Power Deposition, Transport and Acceleration Phenomena in Helicon Plasma Thrusters M. Magarotto, M. Manente, P. de Calro, F. Trezzolani, D. Pavarin, and D. Melazzi CISAS, Padova, Italy
More informationCharacteristics and classification of plasmas
Characteristics and classification of plasmas PlasTEP trainings course and Summer school 2011 Warsaw/Szczecin Indrek Jõgi, University of Tartu Partfinanced by the European Union (European Regional Development
More informationSimulation of the Interaction Between Two Counterflowing Rarefied Jets
Simulation of the Interaction Between Two Counterflowing Rarefied Jets Cyril Galitzine and Iain D. Boyd Department of Aerospace Engineering, University of Michigan, Ann Arbor, MI 48109 Abstract. A preliminary
More informationThe low-field density peak in helicon discharges
PHYSICS OF PLASMAS VOLUME 10, NUMBER 6 JUNE 2003 Francis F. Chen a) Electrical Engineering Department, University of California, Los Angeles, Los Angeles, California 90095-1597 Received 10 December 2002;
More informationModeling nonthermal plasmas generated in glow discharges*
Pure Appl. Chem., Vol. 71, No. 10, pp. 1837±1844, 1999. Printed in Great Britain. q 1999 IUPAC Modeling nonthermal plasmas generated in glow discharges* I. Revel, Ph. Belenguer, J. P. Boeuf and L. C. Pitchford²
More informationMONTE CARLO SIMULATION OF RADIATION TRAPPING IN ELECTRODELESS LAMPS: A STUDY OF COLLISIONAL BROADENERS*
MONTE CARLO SIMULATION OF RADIATION TRAPPING IN ELECTRODELESS LAMPS: A STUDY OF COLLISIONAL BROADENERS* Kapil Rajaraman** and Mark J. Kushner*** **Department of Physics ***Department of Electrical and
More informationModeling and Simulation of Plasma Based Applications in the Microwave and RF Frequency Range
Modeling and Simulation of Plasma Based Applications in the Microwave and RF Frequency Range Dr.-Ing. Frank H. Scharf CST of America What is a plasma? What is a plasma? Often referred to as The fourth
More informationM 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
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 S ep tem b er 6, 2004 Modelling of Dual frequency Capacitively Coupled
More informationChapter 7 Plasma Basic
Chapter 7 Plasma Basic Hong Xiao, Ph. D. hxiao89@hotmail.com www2.austin.cc.tx.us/hongxiao/book.htm Hong Xiao, Ph. D. www2.austin.cc.tx.us/hongxiao/book.htm 1 Objectives List at least three IC processes
More informationKinetic properties of particle-in-cell simulations compromised by Monte Carlo collisions. Abstract
Kinetic properties of particle-in-cell simulations compromised by Monte Carlo collisions M. M. Turner School of Physical Sciences and National Centre for Plasma Science and Technology, Dublin City University,
More informationMODELING OF AN ECR SOURCE FOR MATERIALS PROCESSING USING A TWO DIMENSIONAL HYBRID PLASMA EQUIPMENT MODEL. Ron L. Kinder and Mark J.
TECHCON 98 Las Vegas, Nevada September 9-11, 1998 MODELING OF AN ECR SOURCE FOR MATERIALS PROCESSING USING A TWO DIMENSIONAL HYBRID PLASMA EQUIPMENT MODEL Ron L. Kinder and Mark J. Kushner Department of
More informationSelf-consistent Simulation of the Coupling Between Plasma and Neutral Gas in µn-rit
Self-consistent Simulation of the Coupling Between Plasma and Neutral Gas in µn-rit IEPC-2011-323 Presented at the 32 nd International Electric Propulsion Conference, Wiesbaden, Germany R. Henrich, D.
More informationA MINI-COURSE ON THE PRINCIPLES OF LOW-PRESSURE DISCHARGES AND MATERIALS PROCESSING
A MINI-COURSE ON THE PRINCIPLES OF LOW-PRESSURE DISCHARGES AND MATERIALS PROCESSING Michael A. Lieberman Department of Electrical Engineering and Computer Science, CA 94720 LiebermanMinicourse07 1 OUTLINE
More informationarxiv: v2 [physics.plasm-ph] 31 May 2013
arxiv:128.6519v2 [physics.plasm-ph] 31 May 213 Ionization by bulk heating of electrons in capacitive radio frequency atmospheric pressure microplasmas T Hemke 1, D Eremin 1, T Mussenbrock 1, A Derzsi 2,
More informationCopyright 1996, by the author(s). All rights reserved.
Copyright 1996, by the author(s). All rights reserved. Permission to make digital or hard copies of all or part of this work for personal or classroom use is granted without fee provided that copies are
More informationCOMSOL Multiphysics Simulations of the Electric Field and Gas Flow in a Microwave Axial Injection Torch
COMSOL Multiphysics Simulations of the Electric Field and Gas Flow in a Microwave Axial Injection Torch A. Obrusník, P. Synek, L. Zajíčková Faculty of Science/CEITEC, Masaryk University, Kotlářská 2, 611
More informationPhysique des plasmas radiofréquence Pascal Chabert
Physique des plasmas radiofréquence Pascal Chabert LPP, Ecole Polytechnique pascal.chabert@lpp.polytechnique.fr Planning trois cours : Lundi 30 Janvier: Rappels de physique des plasmas froids Lundi 6 Février:
More informationEffect of Applied Electric Field and Pressure on the Electron Avalanche Growth
Effect of Applied Electric Field and Pressure on the Electron Avalanche Growth L. ZEGHICHI (), L. MOKHNACHE (2), and M. DJEBABRA (3) () Department of Physics, Ouargla University, P.O Box.5, OUARGLA 3,
More informationLANGMUIR PROBE MEASUREMENTS IN 13.56MHz DISCHARGES
LANGMUIR PROBE MEASUREMENTS IN 13.56MHz DISCHARGES A thesis for the degree of PhD Presented to DUBLIN CITY UNIVERSITY by JOHN V. SCANLAN School of Physical Sciences DUBLIN CITY UNIVERSITY Research supervisor
More informationEffects of Plasma Chamber Pressure on the Etching of Micro Structures in SiO 2 With the Charging Effects
IEEE TRANSACTIONS ON PLASMA SCIENCE, VOL. 31, NO. 4, AUGUST 2003 703 Effects of Plasma Chamber Pressure on the Etching of Micro Structures in SiO 2 With the Charging Effects Hye Sang Park, Sung Jin Kim,
More informationEffect of He and Ar Addition on N 2 Glow Discharge Characteristics and Plasma Diagnostics
Arab Journal of Nuclear Science and Applications, 6(1), (116-15) 13 Effect of He and Ar Addition on N Glow Discharge Characteristics and Plasma Diagnostics M. M. Mansour*, N. M. El-Sayed, O. F. Farag and
More informationCurrent sheath formation in the plasma focus
Plasma Science and Applications (ICPSA 2013) International Journal of Modern Physics: Conference Series Vol. 32 (2014) 1460321 (8 pages) The Author DOI: 10.1142/S2010194514603214 Current sheath formation
More informationElectron temperature is the temperature that describes, through Maxwell's law, the kinetic energy distribution of the free electrons.
10.3.1.1 Excitation and radiation of spectra 10.3.1.1.1 Plasmas A plasma of the type occurring in spectrochemical radiation sources may be described as a gas which is at least partly ionized and contains
More informationCopyright 1994, by the author(s). All rights reserved.
Copyright 1994, by the author(s). All rights reserved. Permission to make digital or hard copies of all or part of this work for personal or classroom use is granted without fee provided that copies are
More information16.55 Ionized Gases Problem Set #5
16.55 Ionized Gases Problem Set #5 Problem 1: A probe in a non-maxellian plasma The theory of cold probes, as explained in class, applies to plasmas with Maxwellian electron and ion distributions. However,
More informationNonlinear effects in inductively coupled plasmas a
PHYSICS OF PLASMAS VOLUME 10, NUMBER 5 MAY 003 Nonlinear effects in inductively coupled plasmas a A. I. Smolyakov b) Department of Physics and Engineering Physics, University of Saskatchewan, Saskatoon,
More informationOverview of Accelerated Simulation Methods for Plasma Kinetics
Overview of Accelerated Simulation Methods for Plasma Kinetics R.E. Caflisch 1 In collaboration with: J.L. Cambier 2, B.I. Cohen 3, A.M. Dimits 3, L.F. Ricketson 1,4, M.S. Rosin 1,5, B. Yann 1 1 UCLA Math
More informationLow Temperature Plasma Technology Laboratory
Low Temperature Plasma Technology Laboratory Equilibrium theory for plasma discharges of finite length Francis F. Chen and Davide Curreli LTP-6 June, Electrical Engineering Department Los Angeles, California
More information