Abstract: A hybrid Monte Carlofluid modeling network is developed for an argon-hydrogen mixture, to predict the effect of small amounts of hydrogen added to a dc argon glow discharge. The species considered in the model include the Ar gas atoms, electrons, Ar+ ions and fast Ar atoms, ArH+, H+, H+2 and H+3 ions, and H atoms and H2 molecules, as well as Ar metastable atoms, sputtered Cu atoms, and the corresponding Cu+ ions. Sixty-five reactions between these species are incorporated in the model. The effect of hydrogen on various calculation results is investigated, such as the species densities, the relative role of different production and loss processes for the various species, the cathode sputtering rate and contributions by different bombarding species, and the dissociation degree of H2 and the ionization degree of Ar and Cu. The calculation results are presented and discussed for 1% H2 addition, and comparison is also made with a pure argon discharge and with only 0.1% H2 addition.

Abstract: The electron distribution function (EDF) in an electron cyclotron resonance (ECR) discharge is far from Maxwellian. The self-consistent simulation of ECR discharges requires a calculation of the EDF on every magnetic line for various ion density profiles. The straightforward self-consistent simulation of ECR discharges using the Monte Carlo technique for the EDF calculation is very computer time expensive, since the electron and ion time scales are very different. An electron Boltzmann kinetic equation averaged over the fast electron bouncing and pitch-angle scattering was derived in order to develop an effective and operative tool for the fast modeling (FM) of low-pressure ECR discharges. An analytical solution for the EDF in a loss cone was derived. To check the validity of the FM, one-dimensional (in coordinate) and two-dimensional (in velocity) Monte Carlo simulation codes were developed. The validity of the fast modeling method is proved by comparison with the Monte Carlo simulations. The complete system of equations for FM is presented and ready for use in a comprehensive study of ECR discharges. The variations of plasma density and of wall and sheath potentials are analyzed by solving a self-consistent set of equations for the EDF.

Abstract: The global relaxation process after pulsed laser induced photodetachment in a rf electronegative SIH4 discharge is studied by a self-consistent kinetic one-dimensional particle-in-cell-Monte Carlo model. Our results reveal a comprehensive physical picture of the relaxation process, including the main plasma variables, after a perturbation up to the full recovery of the steady state. A strong influence of the photodetachment on the discharge is found, which results from an increase of the electron density, leading to a weaker bulk field, and hence to a drop in the high energy tail of the electron energy distribution function (EEDF), a reduction of the reaction rates of electron impact attachment and ionization, and a subsequent decrease of the positive and negative ion densities. All the plasma quantities related to electrons recover synchronously. The recovery time of the ion densities is about 1-2 orders of magnitude longer than that of the electrons due to different recovery mechanisms. The modeled behavior of all the charged particles agrees very well with experimental results from the literature. In addition, our work clarifies some unclear processes assumed in the literature, such as the relaxation of the EEDF, the evolution of the electric field, and the recovery of negative ions.