ADVENTURES IN TWO-DIMENSIONAL PARTICLE-IN-CELL SIMULATIONS OF ELECTRONEGATIVE DISCHARGES
|
|
- Hollie Alexander
- 5 years ago
- Views:
Transcription
1 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 University of California Berkeley, CA 9472 Download this talk: lieber LiebEcolePoly9 1
2 OUTLINE Introduction Experiments 2D particlin-cell (PIC) simulations with rescaled oxygen cross sections Theory 1D collisionless model of double layer (DL) Global (volumaveraged) model of upstream and downstream particle and energy balances Comparisons Reasonable agreement but also some differences Slow and fast waves The DL has timvarying structure LiebEcolePoly9 2
3 y z 1 cm 1 cm Upstream Chamber 2 E z = E z cos ωt DL INTRODUCTION Downstream Chamber cm E z = 3. cm 3.5 cm Oxygen Plasma x x (m) Timaveraged Potential (V) (2D PIC 6 mtorr O 2 ) DL Upstream edge 2 Charge density ρ DL + Downstream edge 1 Potential V = Vs Upstream edge 2 DL Downstream edge 1 Why does a DL form at low pressures? The particle loss rate is greater upstream than downstream due to the smaller upstream radius A higher ionization rate (and T e ) is needed upstream than downstream A DL both insulates the low downstream T e from the high upstream T e, and it accelerates electrons upstream to increase the ionization rate there V= LiebEcolePoly9 3
4 PIC SIMULATION METHOD Self-consistent results from first principles with no assumptions about electron and ion velocity distributions Upstream heating at MHz to the plane of the simulation RF field adjusted to keep number of upstream electrons constant Stability, speed and accuracy require small plasma reactors with low densities and large Debye lengths: n e m 3, λ D.8 mm Due to low densities, rescaled oxygen cross sections were used: Positivnegative ion recombination 2 Dissociative attachment 5 A typical simulation takes 1 2 weeks The pressure range explored is.5 24 mtorr A DL was observed for 1 24 mtorr LiebEcolePoly9 4
5 PIC RESULTS FOR 6 mtorr O 2 DISCHARGE Axial Potential (V) Radially averaged T (V) Th 2 15 DL 4 3 Tc 1 2 O- 5 4e+15 3e+15 2e O- O + 2 x (m) Axial Density (m-3) n /ne ~ O x (m) K Axial Collision Rates (m 3 /s) iz high Kiz 1e x (m) 516 K att low Kiz x (m) LiebEcolePoly9 5
6 RESULTS FOR 6 mtorr O 2 DISCHARGE (CONT D) Charge Density Components in DL (V/m 2 ) DL charge Transverse ion loss in DL Total ρ/ε Ex/ x Ey/ y Ex/ x + Ey/ y x (m) Up edge x-component EEDF (a.u.) 1e+18 1e+17 1e+16 1e Up region Up edge Down region 3.95 EEPF (a.u.) Down edge e ε (ev) Down-edge x-component EEDF (a.u.).1.1 up-directed Accelerated component Tc Th down-directed.1.1 up-directed Tc down-directed Th x-energy (V) x-energy (V) LiebEcolePoly9 6
7 DL MODEL Charge density ρ and potential V within the DL are found by solving Poisson s equation Potential V = Vs Upstream edge 2 Six types of particles contribute to ρ: thermal positive ions, negative ions, and electrons positive ions accelerated downstream negative ions and electrons accelerated upstream The particle motions are 1D and collisionless The boundary conditions are that ρ and dρ/dv vanish at the DL edges DL Downstream edge 1 An additional condition is that the sum of positive and negative charge in the double layer vanishes; equivalently, the total force acting on the double layer vanishes A final condition that upstream and downstream-directed electron fluxes nearly balance determines the equilibrium value of V s /T h. V= LiebEcolePoly9 7
8 PARTICLE AND ENERGY BALANCES Upstream region 2 DL Downstream region 1 We use a 2D rectangular geometry To determine the equilibrium quantities we use global particle balance upstream global particle balance downstream global energy balance downstream Upstream energy balance (which determines the upstream electron density) depends on the input power; we use the PIC simulation value for the comparisons LiebEcolePoly9 8
9 MODEL (TOP) and PIC (BOTTOM) EEDF S (6 mtorr O 2 ) Up edge x-component EEDF up-directed Accelerated component Tc Th down-directed Down edge x-component EEDF up-directed Tc/Th=.668 Tc down-directed Th x-energy (V) Up edge x-component EEDF (a.u.) x-energy (V) Down-edge x-component EEDF (a.u.).1.1 up-directed Accelerated component Tc Th down-directed.1.1 up-directed Tc down-directed Th x-energy (V) x-energy (V) LiebEcolePoly9 9
10 MODEL (WHITE) AND PIC (BLACK) FOR 6 mtorr O 2 Te (V) Up elec bal (#/s) 1 2e e e e+15 Th Thdn Tc Tcdn Iz Att Gam- Gam+ Wall 8e+17 7e+17 6e+17 5e+17 4e+17 3e+17 2e+17 1e+17 Down energy bal (ev/s) 3.5e+15 3e e+15 2e e+15 1e+15 5e+14 Down elec balance (#/s) No transverse ion loss in DL Edn Eup Idn Nup Coll Wall Iz Att Gam- Gam+ Wall LiebEcolePoly9 1
11 MODEL (DASH) AND PIC (CIRCLES) VS PRESSURE 1 Vs/Th 1 1 p (mtorr) Th (V) 1 1 p (mtorr) Vs (V) p (mtorr) ne2/ne LiebEcolePoly9 p (mtorr) 11
12 SLOW AND FAST WAVES: INSTABILITIES (6 mtorr O 2 ) 1 Ex(t) at x=1 cm (dn DL edge) time interval = 59 microsecs 1 FFT of Ex(t) at x = 1 cm (dn DL edge) 8.475e4 Hz e5 Hz time (secs) 1 2e+5 4e+5 6e+5 8e+5 1e+6 frequency (Hz) At 2 12 mtorr, the DL coexists with an unstable slow wave that originates downstream and propagates upstream as it grows The wave frequency is 5 1 khz with a wavelength of order 1 cm The wave produces 2% oscillations in the double layer potential and.5 cm oscillations in the DL position We believe the wave is driven by counter-streaming flows of positive and negative ions (Tuszewski and Gary, 23) There is also a fast wave that increases spatially along with the slow wave at the higher pressures, with f 1 MHz and λ 1 cm LiebEcolePoly9 12
13 MOVIE SHOWING SLOW AND FAST WAVES Red solid line: fast waves averaged over.1475 µs intervals Blue solid line: slow waves averaged over 1.18 µs intervals Slow Fast DL LiebEcolePoly9 13
14 ADVENTURES IN TWO-DIMENSIONAL PARTICLE-IN-CELL SIMULATIONS OF ELECTRONEGATIVE DISCHARGES PART 2: TRANSPORT IN A MAGNETIZED DISCHARGE E. Kawamura, M.A. Lieberman and A.J. Lichtenberg University of California Berkeley, CA 9472 LiebEcolePoly9 14
15 OUTLINE Introduction Experiments with uniform electron temperature 2D PIC simulations with rescaled oxygen cross sections Model and comparisons New model (differs from Gary et al, to appear in J. Phys. D) Experiments with hot core and cold edge 2D PIC simulations with iodinlike cross sections Two-region model and comparisons LiebEcolePoly9 15
16 INTRODUCTION Rectangular grounded chamber with insulating walls Spatially varying axial rf current density J x y L=2 cm Jx 2R=3 cm Jx(x,y,t) B = 6 G Plasma (2 mtorr oxygen) Jx Dielectric κ=4.25 cm x For an unmagnetized discharge, negative ions are confined by the (positive) plasma potential For a magnetized discharge with uniform axial B : What electron temperature profiles can be obtained? What is the negative ion flux to the transverse walls? What are the electron fluxes to the axial and transverse walls? What are the transverse diffusion coefficients for positive ions, negative ions, and electrons? LiebEcolePoly9 16
17 PIC SIMULATION METHOD MHz axial RF current density J x Adjust amplitude and spatial variation to obtain T e const J x = 1 sin 4πx πy sin L 2R [A/m2 ] Transverse heating (E y s or E z s) did not work Larger amplitudes or uniform J x s large T e s at the transverse and/or axial edges B = 6 G (sometimes 3 or 12 G) Stability, speed, accuracy, and T e const require low densities and large Debye lengths: n e m 3, λ D.8 cm (!) Due to low density, rescaled oxygen cross sections were used: Positivnegative ion recombination 2, 1, etc A typical simulation takes 1 2 weeks LiebEcolePoly9 17
18 PIC RESULTS FOR 1 K rec O 2 DISCHARGE Radial Potential (V) Radial T (ev) xO- 1xO Radial Density (m-3) Axial Density (m-3) 1e+14 O2+ 1e+14 O2+ 8e+13 O- 8e+13 O- 6e+13 6e+13 4e+13 2e+13 4e+13 2e x (m) LiebEcolePoly9 18
19 PIC RESULTS FOR OTHER O 2 DISCHARGES 2e e+14 Radial Density (m-3) (a) B x =6 G, L=19.5 cm, K recfac=2 Mag. ions O2+ 2e e+14 Radial Density (m-3) (b) B x =12 G, L=19.5 cm, K recfac=2 Unmag. ions O- O2+ 1e+14 1e+14 O- 5e+13 5e e e+14 Radial Density (m-3) (c) B x =6 G, L= 9.5 cm, K recfac =2 Mag. ions O2+ O- 2e e+14 Radial Density (m-3) (d) B x =6 G, L= 9.5 cm, K recfac =2 Mag. ions 1e+14 1e+14 O2+ O- 5e+13 5e LiebEcolePoly9 19
20 CONCLUSIONS FROM PIC SIMULATIONS Cosine transverse electron and negative ion density profiles constant electronegativity α = n n /n e n n n n n e n e Leray et al model Transverse direction What PIC shows Transverse direction Transverse electron/negative ion flux Γ e /Γ n D e /α D n Leray et al model confirmed for axial loss frequency ν L 2D ( i 1 + T ) 1/2 e RL T i Transverse diffusion coefficients D e, D i, and D n are classical LiebEcolePoly9 2
21 1D MODEL Electron balance D e d 2 n e dy 2 = (ν iz ν att )n e ν L (n e + n n ) Negative ion balance D n d 2 n n dy 2 = ν attn e K rec n n (n e + n n ) D e and D n are the transverse ambipolar diffusion coefficients, obtained from Γ i = Γ e + Γ n in terms of the species mobilities and diffusion coefficients Cosine density profiles with n n = α n e (from PIC) Boundary conditions n e = n n = at transverse wall good agreement with PIC results LiebEcolePoly9 21
22 SCALINGS FOR 1D MODEL Transverse wall negative ion current = 2π L R D nα n e Transverse wall electron current = 2π L R D e n e End wall particle current = 2D i (1 + α )(1 + T e T i ) 1/2 n e Negative ion balance: 4R π ν att = π R D nα + K rec n e α (1 + α ) Range of solutions: < α < 4 π 2 R 2 D n LiebEcolePoly9 22
23 TWO-REGION IODINE SIMULATIONS Goal was to obtain a high T e core with a cold T e periphery Rectangular grounded chamber with insulating walls MHz J x peaked at transverse midplane NOTE Jx! y 2R=3 cm Jx(x,y,t) Jx L=2 cm B =6 G Plasma (2 mtorr iodine) Rescaled oxygen cross sections were used: Positivnegative ion recombination 2 Attachment 5 (w/wo 4.2 V threshold energy) Iodine masses were used Dielectric κ=4.25 cm x LiebEcolePoly9 23
24 PIC RESULTS FOR IODINE DISCHARGE Radial Density (m ) -3 Radial T (ev) mt 2xKrec, 5xKatt(No Thr.) Dielectrics, 3x2 cm, B= 6G x=5 to 15 cm 5xI mt Iodine 2xKrec, 5xKatt(No Thr.) Discharge 1e+15 1e+14 1e+13 1e+12 Dielectrics, 3x2cm, B= 6 G 1e11α I+ I- Radial Density (m ) -3 Radial T (ev) mt Iodine 2xKrec, 5xKatt(4.2V Thr.) Dielectrics, 3x2 cm, B= 6G x=5 to 15 cm 5xI- 5xI- 5xI mt Iodine 2xKrec, 5xKatt(4.2V Thr.) Discharge 1e+15 1e+14 1e+13 1e+12 Dielectrics, 3x2cm, B= 6 G 1e11α I+ I- 1e No attachment threshold 1e V attachment threshold LiebEcolePoly9 24
25 TWO-REGION 1D MODEL Electron balance D e d 2 n e dy 2 = (ν iz ν att )n e ν L (n e + n n ) (ν iz in the core, ν iz = in the periphery) Negative ion balance D n d 2 n n dy 2 = ν attn e (Neglect recombination compared to negative ion transverse loss) Match densities and fluxes at the corperiphery transition Boundary conditions n e = n n = at transverse wall This linear system can be solved exactly LiebEcolePoly9 25
26 COMPARISON OF MODEL WITH PIC 1 15 unmagions with no thr, Ri=.5 cm 1 15 unmagions with thr, Ri=.5 cm n n, n i n n, n i E11α n e n e 1E11α Radial Density (m ) -3 1e+15 1e+14 1e+13 1e mt Iodine 2xKrec, 5xKatt(No Thr.) Discharge Dielectrics, 3x2cm, B= 6 G I+ I- 1e11α Radial Density (m ) -3 1e+15 1e+14 1e+13 1e mt Iodine 2xKrec, 5xKatt(4.2V Thr.) Discharge Dielectrics, 3x2cm, B= 6 G I+ I- 1e11α 1e No attachment threshold 1e V attachment threshold LiebEcolePoly9 26
27 CONCLUSIONS 2D PIC simulations can be powerful tools to study the physics of electronegative discharges Experimental parameters can be easily varied; various physics can be turned on and off These simulations can provide diagnostics that would be very difficult to do in laboratory experiments The models, verified by PIC, can give important scalings that can be used to design new PIC and experimental configurations Download this talk: lieber LiebEcolePoly9 27
MODELING AND SIMULATION OF LOW TEMPERATURE PLASMA DISCHARGES
MODELING AND SIMULATION OF LOW TEMPERATURE PLASMA DISCHARGES Michael A. Lieberman University of California, Berkeley lieber@eecs.berkeley.edu DOE Center on Annual Meeting May 2015 Download this talk: http://www.eecs.berkeley.edu/~lieber
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 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 informationMODELING PLASMA PROCESSING DISCHARGES
MODELING PROCESSING DISCHARGES M.A. Lieberman Department of Electrical Engineering and Computer Sciences University of California Berkeley, CA 94720 Collaborators: E. Kawamura, D.B. Graves, and A.J. Lichtenberg,
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 informationNONLINEAR ELECTROMAGNETICS MODEL OF AN ASYMMETRICALLY DRIVEN CAPACITIVE DISCHARGE
NONLINEAR ELECTROMAGNETICS MODEL OF AN ASYMMETRICALLY DRIVEN CAPACITIVE DISCHARGE M.A. Lieberman Department of Electrical Engineering and Computer Sciences University of California Berkeley, CA 94720 Collaborators:
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 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 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 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 informationA global (volume averaged) model of a chlorine discharge
A global (volume averaged) model of a chlorine discharge Eyþór Gísli Þorsteinsson and Jón Tómas Guðmundsson Science Institute, University of Iceland, Iceland Department of Electrical and Computer Engineering,
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 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 informationPlasmas rf haute densité Pascal Chabert LPTP, Ecole Polytechnique
Plasmas rf haute densité Pascal Chabert LPTP, Ecole Polytechnique chabert@lptp.polytechnique.fr Pascal Chabert, 2006, All rights reserved Programme Introduction Généralité sur les plasmas Plasmas Capacitifs
More informationDiffusion during Plasma Formation
Chapter 6 Diffusion during Plasma Formation Interesting processes occur in the plasma formation stage of the Basil discharge. This early stage has particular interest because the highest plasma densities
More informationLow Temperature Plasma Technology Laboratory
Low Temperature Plasma Technology Laboratory CENTRAL PEAKING OF MAGNETIZED GAS DISCHARGES Francis F. Chen and Davide Curreli LTP-1210 Oct. 2012 Electrical Engineering Department Los Angeles, California
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 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 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 informationThe Role of Secondary Electrons in Low Pressure RF Glow Discharge
WDS'05 Proceedings of Contributed Papers, Part II, 306 312, 2005. ISBN 80-86732-59-2 MATFYZPRESS The Role of Secondary Electrons in Low Pressure RF Glow Discharge O. Brzobohatý and D. Trunec Department
More informationNonlinear Diffusion in Magnetized Discharges. Francis F. Chen. Electrical Engineering Department
Nonlinear Diffusion in Magnetized Discharges Francis F. Chen Electrical Engineering Department PPG-1579 January, 1998 Revised April, 1998 Nonlinear Diffusion in Magnetized Discharges Francis F. Chen Electrical
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 informationCharacterization of an Oxygen Plasma by Using a Langmuir Probe in an Inductively Coupled Plasma
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
More information2D Hybrid Fluid-Analytical Model of Inductive/Capacitive Plasma Discharges
63 rd GEC & 7 th ICRP, 2010 2D Hybrid Fluid-Analytical Model of Inductive/Capacitive Plasma Discharges E. Kawamura, M.A. Lieberman, and D.B. Graves University of California, Berkeley, CA 94720 This work
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 informationThe ideal Maxwellian plasma
The ideal Maxwellian plasma Dr. L. Conde Departamento de Física Aplicada. E.T.S. Ingenieros Aeronáuticos Universidad Politécnica de Madrid Plasmas are,... The plasma state of matter may be defined as a
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 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 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 informationA theory for formation of a low pressure, current-free double layer
INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS D: APPLIED PHYSICS J. Phys. D: Appl. Phys. 39 (6) 394 334 doi:.88/-377/39/5/ A theory for formation of a low pressure, current-free double layer M A Lieberman,
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 informationPlasma Astrophysics Chapter 1: Basic Concepts of Plasma. Yosuke Mizuno Institute of Astronomy National Tsing-Hua University
Plasma Astrophysics Chapter 1: Basic Concepts of Plasma Yosuke Mizuno Institute of Astronomy National Tsing-Hua University What is a Plasma? A plasma is a quasi-neutral gas consisting of positive and negative
More informationLow-field helicon discharges
Plasma Phys. Control. Fusion 39 (1997) A411 A420. Printed in the UK PII: S0741-3335(97)80958-X Low-field helicon discharges F F Chen, X Jiang, J D Evans, G Tynan and D Arnush University of California,
More informationEffect of a dielectric layer on plasma uniformity in high frequency electronegative capacitive discharges
Effect of a dielectric layer on plasma uniformity in high frequency electronegative capacitive discharges Emi KawamuraDe-Qi WenMichael A. Lieberman and Allan J. Lichtenberg Citation: Journal of Vacuum
More informationTwo-dimensional Particle-In-Cell model of the extraction region of the PEGASES ion-ion plasma source
Two-dimensional Particle-In-Cell model of the extraction region of the PEGASES ion-ion plasma source IEPC-2013-249 Presented at the 33rdInternational Electric Propulsion Conference, The George Washington
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 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 informationExperimental investigation of double layers in expanding plasmas
PHYSICS OF PLASMAS 14, 013506 2007 Experimental investigation of double layers in expanding plasmas N. Plihon, a P. Chabert, and C. S. Corr Laboratoire de Physique et Technologie des Plasmas, Ecole Polytechnique,
More informationPIC-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 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 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 informationExperiments and theory of an upstream ionization instability excited by an accelerated electron beam through a current-free double layer
PHYSICS OF PLASMAS 13, 122101 2006 Experiments and theory of an upstream ionization instability excited by an accelerated electron beam through a current-free double layer A. Aanesland, a M. A. Lieberman,
More information2D OOPIC Simulations of the Helicon Double Layer
2D OOPIC Simulations of the Helicon Double Layer IEPC-2007-146 Presented at the 30 th International Electric Propulsion Conference, Florence, Italy I. Musso * Center for Studies and Activities for Space,
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 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 informationAn Improved Global Model for Electronegative Discharge and Ignition Conditions for Peripheral Plasma Connected to a Capacitive Discharge
An Improved Global Model for Electronegative Discharge and Ignition Conditions for Peripheral Plasma Connected to a Capacitive Discharge Sungjin Kim Electrical Engineering and Computer Sciences University
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 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 informationIncreased Upstream Ionization Due to Spontaneous Formation of a Double Layer in an Expanding Plasma
Increased Upstream Ionization Due to Spontaneous Formation of a Double Layer in an Expanding Plasma Earl E. Scime* November, 2009 APS Division of Plasma Physics Meeting * with Costel Biloiu, Ioana Biloiu,
More informationExperimental Investigations of Magnetic Reconnection. J Egedal. MIT, PSFC, Cambridge, MA
Experimental Investigations of Magnetic Reconnection J Egedal MIT, PSFC, Cambridge, MA Coronal Mass Ejections Movie from NASA s Solar Dynamics Observatory (SDO) Space Weather The Solar Wind affects the
More informationLow Temperature Plasma Technology Laboratory
Low Temperature Plasma Technology Laboratory CHARACTERIZATION AND MODELING OF THE MØRI SOURCE Francis F. Chen and Donald Arnush FINAL REPORT, UC MICRO PROJECT 98-2 JOINTLY FUNDED BY APPLIED MATERIALS LTP-2
More informationCharacterization of Ion Cyclotron Resonance Acceleration for Electric Propulsion with Interferometry
JOURNAL OF PROPULSION AND POWER Vol. 27, No. 2, March April 2011 Characterization of Ion Cyclotron Resonance Acceleration for Electric Propulsion with Interferometry Christopher Davis ElectroDynamic Applications,
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 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 informationSupplementary Information
1 Supplementary Information 3 Supplementary Figures 4 5 6 7 8 9 10 11 Supplementary Figure 1. Absorbing material placed between two dielectric media The incident electromagnetic wave propagates in stratified
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 informationExperiments with a Supported Dipole
Experiments with a Supported Dipole Reporting Measurements of the Interchange Instability Excited by Electron Pressure and Centrifugal Force Introduction Ben Levitt and Dmitry Maslovsky Collisionless Terrella
More informationDesign of an Experiment for Studying Ion Acceleration by Beating Waves
Design of an Experiment for Studying Ion Acceleration by Beating Waves R. Spektor and E.Y. Choueiri Electric Propulsion and Plasma Dynamics Laboratory (EPPDyL) Mechanical and Aerospace Engineering Department
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 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 informationTURBULENT TRANSPORT THEORY
ASDEX Upgrade Max-Planck-Institut für Plasmaphysik TURBULENT TRANSPORT THEORY C. Angioni GYRO, J. Candy and R.E. Waltz, GA The problem of Transport Transport is the physics subject which studies the physical
More informationModélisation de sources plasma froid magnétisé
Modélisation de sources plasma froid magnétisé Gerjan Hagelaar Groupe de Recherche Energétique, Plasma & Hors Equilibre (GREPHE) Laboratoire Plasma et Conversion d Énergie (LAPLACE) Université Paul Sabatier,
More informationPropagating double layers in electronegative plasmas
PHYSICS OF PLASMAS 14, 053508 2007 Propagating double layers in electronegative plasmas A. Meige a Laboratoire Plasma et Conversion d Energie (LAPLACE), Université Paul Sabatier, 118 route de Narbonne,
More informationMeasurement of EEDF and Distribution of Primary Electron in a Bucket Ion Source
Fourth IAEA Technical Meeting on "Negative Ion Based NBIs" May 9-11 2005, Padova - Italy Measurement of EEDF and Distribution of Primary Electron in a Bucket Ion Source S. Miyamoto, F. Kanayama, T. Minami,
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 informationDr. A A Mamun Professor of Physics Jahangirnagar University Dhaka, Bangladesh
SOLITARY AND SHOCK WAVES IN DUSTY PLASMAS Dr. A A Mamun Professor of Physics Jahangirnagar University Dhaka, Bangladesh OUTLINE Introduction Static Dust: DIA Waves DIA Solitary Waves DIA Shock Waves Mobile
More informationMODELING THE ELECTROMAGNETIC FIELD
MODELING THE ELECTROMAGNETIC FIELD IN ANISOTROPIC INHOMOGENEOUS MAGNETIZED PLASMA OF ECR ION SOURCES G. Torrisi (1), D. Mascali (1), A. Galatà (2), G. Castro (1), L. Celona (1), L. Neri (1), G. Sorbello
More informationDynamic Phenomena in Complex Plasmas
The University of Sydney Dynamic Phenomena in Complex Plasmas N.F. Cramer, S.V. Vladimirov, A.A. Samarian and B.W. James School of Physics, University of Sydney, Australia Dusty Plasmas at the University
More informationOne-dimensional particle-in-cell simulation of a current-free double layer in an expanding plasma
One-dimensional particle-in-cell simulation of a current-free double layer in an expanding plasma Albert Meige, Rod W. Boswell, Christine Charles, and Miles M. Turner Citation: Physics of Plasmas (1994-present)
More informationAxial symmetric open magnetic traps with depressed transversal losses of plasmas
Axial symmetric open magnetic traps with depressed transversal losses of plasmas A. Sidorov, S. Golubev, I. Izotov, S. Razin, V. Skalyga and A. Vodopyanov Institute of Applied Physics, RAS, 603950 Nizhny
More informationRECENT PROGRESS ON THE PHYSICS
RECENT PROGRESS ON THE PHYSICS OF CAPACITIVE DISCHARGES M.A. Lieberman Department of Electrical Engineering and Computer Sciences University of California Berkeley, CA 94720 Download this talk: http://www.eecs.berkeley.edu/
More informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
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 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 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 informationPenning Traps. Contents. Plasma Physics Penning Traps AJW August 16, Introduction. Clasical picture. Radiation Damping.
Penning Traps Contents Introduction Clasical picture Radiation Damping Number density B and E fields used to increase time that an electron remains within a discharge: Penning, 936. Can now trap a particle
More informationVolume Production of D - Negative Ions in Low-Pressure D 2 Plasmas - Negative Ion Densities versus Plasma Parameters -
Volume Production of D - Negative Ions in Low-Pressure D 2 Plasmas - Negative Ion Densities versus Plasma Parameters - Osamu Fukumasa and Shigefumi Mori Department of Electrical and Electronic Engineering,
More informationAdditional Heating Experiments of FRC Plasma
Additional Heating Experiments of FRC Plasma S. Okada, T. Asai, F. Kodera, K. Kitano, T. Suzuki, K. Yamanaka, T. Kanki, M. Inomoto, S. Yoshimura, M. Okubo, S. Sugimoto, S. Ohi, S. Goto, Plasma Physics
More informationModelling of low-temperature plasmas: kinetic and transport mechanisms. L.L. Alves
Modelling of low-temperature plasmas: kinetic and transport mechanisms L.L. Alves llalves@tecnico.ulisboa.pt Instituto de Plasmas e Fusão Nuclear Instituto Superior Técnico, Universidade de Lisboa Lisboa,
More information26. Non-linear effects in plasma
Phys780: Plasma Physics Lecture 26. Non-linear effects. Collisionless shocks.. 1 26. Non-linear effects in plasma Collisionless shocks ([1], p.405-421, [6], p.237-245, 249-254; [4], p.429-440) Collisionless
More informationNon-Linear Plasma Wave Decay to Longer Wavelength
Non-Linear Plasma Wave Decay to Longer Wavelength F. Anderegg 1, a), M. Affolter 1, A. Ashourvan 1, D.H.E. Dubin 1, F. Valentini 2 and C.F. Driscoll 1 1 University of California San Diego Physics Department
More informationModeling neutral-plasma interactions in scrape-off layer (SOLT) simulations*
Modeling neutral-plasma interactions in scrape-off layer (SOLT) simulations* D. A. Russell and J. R. Myra Research Corporation Boulder CO USA Presented at the US Transport Task Force Workshop Williamsburg
More informationComsol Multiphysics 在低溫電漿模擬之應用 ~My Little Journal with Simulation
Comsol Multiphysics 在低溫電漿模擬之應用 ~My Little Journal with Simulation 多重物理 CAE 分析軟體 COMSOL Multiphysics CONFERENCE 用戶研討會 Nov 9, 2012 台灣大學化工系徐振哲 Start of Plasma Plasma: ionized gas with equal amount of positive
More informationDiffusion equation, flux, diffusion coefficient scaling. Diffusion in fully ionized plasma vs. weakly ionized plasma. n => Coulomb collision frequency
Last Time Diffusion in plasma: the main reason why we need to control it (i.e. using magnetic field) Diffusion equation, flux, diffusion coefficient scaling o o t nx,t Dn D2 nx,t o D ~ L 2 T Diffusion
More informationGradient drift instability in Hall plasma devices
Gradient drift instability in Hall plasma devices IEPC-3-75 Presented at the 33 rd International Electric Propulsion Conference, The George Washington University, Washington, D.C., USA October 6, 3 W.Frias
More informationStability of strong waves and its implications for pulsar wind shocks
Stability of strong waves and its implications for pulsar wind shocks Iwona Mochol in collaboration with John Kirk Max-Planck-Institut für Kernphysik Heidelberg, Germany Madrid, 22/05/2013 Pulsar winds
More informationTurbulence and Transport The Secrets of Magnetic Confinement
Turbulence and Transport The Secrets of Magnetic Confinement Presented by Martin Greenwald MIT Plasma Science & Fusion Center IAP January 2005 FUSION REACTIONS POWER THE STARS AND PRODUCE THE ELEMENTS
More informationSpatial distributions of power and ion densities in RF excited remote plasma reactors
Plasma Sources Sci. Technol. 5 (1996) 499 509. Printed in the UK Spatial distributions of power and ion densities in RF excited remote plasma reactors Irène Pérès and Mark J Kushner University of Illinois,
More informationThe classical model of a Hall thruster is based on electron conduction across magnetic field lines being
Self-Consistent Calculation of Electron Transport in a Cylindrical Hall Thruster Lubos Brieda, Michael Keidar The George Washington University, Washington, D.C. 20052 Yevgeny Raitses and Nathaniel J. Fisch
More informationIs plasma important? Influence molecule formation?
Is plasma important? Influence molecule formation? Plasma Structure (space & time) Influence? Daan Schram Eindhoven University of Technology d.c.schram@tue.nl http://www.tue.nl/en/employee/ep/e/d/ep-uid/19780797/?no_cache=1&chash=e23e831cf0c6bebeac6023f04dd3c4b6
More informationCluster fusion in a high magnetic field
Santa Fe July 28, 2009 Cluster fusion in a high magnetic field Roger Bengtson, Boris Breizman Institute for Fusion Studies, Fusion Research Center The University of Texas at Austin In collaboration with:
More informationFundamentals of Plasma Physics Transport in weakly ionized plasmas
Fundamentals of Plasma Physics Transport in weakly ionized plasmas APPLAuSE Instituto Superior Técnico Instituto de Plasmas e Fusão Nuclear Luís L Alves (based on Vasco Guerra s original slides) 1 As perguntas
More informationPhysics and Modelling of a Negative Ion Source Prototype for the ITER Neutral Beam Injection
1 ITR/P1-37 Physics and Modelling of a Negative Ion Source Prototype for the ITER Neutral Beam Injection J.P. Boeuf a, G. Fubiani a, G. Hagelaar a, N. Kohen a, L. Pitchford a, P. Sarrailh a, and A. Simonin
More informationCharacteristics of Positive Ions in the Sheath Region of Magnetized Collisional Electronegative Discharges
Plasma Science and Technology, Vol.6, No.6, Jun. 204 Characteristics of Positive Ions in the Sheath Region of Magnetized Collisional Electronegative Discharges M. M. HATAMI, A. R. NIKNAM 2 Physics Department
More informationState-to-State Kinetics of Molecular and Atomic Hydrogen Plasmas
State-to-State Kinetics of Molecular and Atomic Hydrogen Plasmas MARIO CAPITELLI Department of Chemistry, University of Bari (Italy) CNR Institute of Inorganic Methodologies and Plasmas Bari (Italy) MOLECULAR
More informationUltra-Cold Plasma: Ion Motion
Ultra-Cold Plasma: Ion Motion F. Robicheaux Physics Department, Auburn University Collaborator: James D. Hanson This work supported by the DOE. Discussion w/ experimentalists: Rolston, Roberts, Killian,
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 informationPositive ion flux from a low-pressure electronegative discharge
Plasma Sources Sci. Technol. 8 (1999) 457 462. Printed in the UK PII: S0963-0252(99)03994-8 Positive ion flux from a low-pressure electronegative discharge T E Sheridan, P Chabert and R W Boswell Space
More information1. INTRODUCTION 2. EXPERIMENTAL SET-UP CHARACTERIZATION OF A TUBULAR PLASMA REACTOR WITH EXTERNAL ANNULAR ELECTRODES
Romanian Reports in Physics, Vol. 57, No. 3, P. 390-395, 2005 CHARACTERIZATION OF A TUBULAR PLASMA REACTOR WITH EXTERNAL ANNULAR ELECTRODES C. PETCU, B. MITU, G. DINESCU National Institute for Lasers,
More informationEffect of parallel velocity shear on the excitation of electrostatic ion cyclotron waves
4 February 2002 Physics Letters A 293 (2002) 260 265 www.elsevier.com/locate/pla Effect of parallel velocity shear on the excitation of electrostatic ion cyclotron waves E.P. Agrimson, N. D Angelo, R.L.
More information