Προσομοίωση ηλεκτρομαγνητικής συμπεριφοράς σε αντιδραστήρες αερίων χαμηλής πίεσης και ασθενούς ιονισμού

Abstract

Inductively Coupled Plasma Sources (ICP’s) are capable of producing high density-low pressure plasmas in a variety of applications for the semiconductor and material processing industry. In the mTorr range, ICP discharges exhibit an extraordinary behaviour concerning the electromagnetic field propagation and wave-particle interaction: Anomalous skin effect and resonant wave-particle interaction within a superimposed static magnetic field consist two of the most typical phenomena. The efficient comprehension and mathematical description of such a complex gas discharge in order to fast and accurately simulate ICP sources, is still a challenging task. Within this context, the thesis focuses on evaluating the existing scientific knowledge in plasma computational modeling in order to develop not only rapidly converging but reliable ICP simulations: The implementation methodology consists on formulating an hypothesis (model) and repetitively inquiring its accuracy by checking the simulation results against existing experimental and/or other simulation data. The continuation of the model re-formulation process depends on the accuracy of the simulation results. Initally a simulation of a Magnetic Pole Enhanced (MaPE)-ICP plasma source was developed, under the assumptions of plasma quasineutrality and ambipolar diffusion. The simulation results were checked against the experimental data of previous workers for Argon discharges and the ability of the model to provide an elementary quantitative description of low pressure ICP sources was scrutinized. The validity of the fluid model was enhanced with the incorporation of a time effective numerical solution of the Boltzmann transport equation for electrons. Simulation results of the hybrid model were compared to the previous fluid simulation results and existing experimental data, for a 30 mTorr Argon discharge in the MaPE–ICP reactor. The major improvements of the qualitative agreement in regard to the effect of parameters with particular interest are discussed. Moreover, a fluid model of ECWR (Electron Cyclotron Wave Resonance) discharges, based on predefined boundary conditions for the electromagnetic field, was developed: The simulation results for a 15 mTorr Argon plasma within a slab (1-dimensional model) were checked against the particle in cell/Monte Carlo (PIC/MC) simulation results that can be found in the literature and also compared to the analytical theory and experimental data. In addition, the model was further developed to simulate realistic geometries as a cylindrical ECWR reactor (2-D) and the data were also compared to both simulation results and experimental data of other researchers. Finally, the model was extended in order to simulate plasma propagation from a typical ECWR plasma source to a diffusion region. The simulation results for an Argon plasma generated from a cylindrical ECWR source in a processing chamber at 5 mTorr were presented in order to verify the feasibility of model application in ECWR plasma processes.