Abstract: The chemical complexity of non-equilibrium plasmas poses a challenge for plasma modeling because of the computational load. This paper presents a dimension reduction method for such chemically complex plasmas based on principal component analysis (PCA). PCA is used to identify a low-dimensional manifold in chemical state space that is described by a small number of parameters: the principal components. Reduction is obtained since continuity equations only need to be solved for these principal components and not for all the species. Application of the presented method to a CO2 plasma model including state-to-state vibrational kinetics of CO2 and CO demonstrates the potential of the PCA method for dimension reduction. A manifold described by only two principal components is able to predict the CO2 to CO conversion at varying ionization degrees very accurately.

Abstract: An inductively coupled plasma connected to a mass spectrometer interface (sampling cone) is computationally investigated. Typical plasma characteristics, such as gas flow velocity, plasma temperature and electron density, are calculated in two dimensions (cylindrical symmetry) and compared with and without a mass spectrometer sampling interface. The results obtained from our model compare favorably with experimental data reported in the literature. A dramatic increase in the plasma velocity is reported in the region close to the interface. Furthermore, a cooled metal interface lowers the plasma temperature and electron density on the axial channel very close to the sampling cone but the corresponding values in the off axial regions are increased. Therefore, the effect of the interface strongly depends on the measurement position. It is shown that even a small shift from the actual position of the sampler leads to a considerable change of the results.

Abstract: This paper demonstrates that the CO2 conversion in a dielectric barrier discharge rises drastically upon addition of Ar or He, and the effect is more pronounced for Ar than for He. The effective CO2 conversion, on the other hand, drops upon addition of Ar or He, which is logical due to the lower CO2 content in the gas mixture, and the same is true for the energy efficiency, because a considerable fraction of the energy is then consumed into ionization/excitation of Ar or He atoms. The higher absolute CO2 conversion upon addition of Ar or He can be explained by studying in detail the Lissajous plots and the current profiles. The breakdown voltage is lower in the CO2/Ar and CO2/He mixtures, and the discharge gap is more filled with plasma, which enhances the possibility for CO2 conversion. The rates of electron impact excitationdissociation of CO2, estimated from the electron densities and mean electron energies, are indeed higher in the CO2/Ar and (to a lower extent) in the CO2/He mixtures, compared to the pure CO2 plasma. Moreover, charge transfer between Ar+ or Ar2+ ions and CO2, followed by electron-ion dissociative recombination of the CO2+ ions, might also contribute to, or even be dominant for the CO2 dissociation. All these effects can explain the higher CO2 conversion, especially upon addition of Ar, but also upon addition of He.

Abstract: The He-Ar-Cu+ IR laser operates in a hollow-cathode discharge, typically in a mixture of helium with a few-% Ar. The population inversion of the Cu+ ion levels, responsible for laser action, is attributed to asymmetric charge transfer between He+ ions and sputtered Cu atoms. The Ar gas is added to promote sputtering of the Cu cathode. In this paper, a hybrid modeling network consisting of several different models for the various plasma species present in a He-Ar-Cu hollow-cathode discharge is applied to investigate the effect of Ar concentration in the gas mixture on the discharge behavior, and to find the optimum He/Ar gas ratio for laser operation. It is found that the densities of electrons, Ar+ ions, Ar-m* metastable atoms, sputtered Cu atoms and Cu+ ions increase upon the addition of more Ar gas, whereas the densities of He+ ions, He-2(+) ions and He-m* metastable atoms drop considerably. The product of the calculated Cu atom and He+ ion densities, which determines the production rate of the upper laser levels, and hence probably also the laser output power, is found to reach a maximum around 1-5% Ar addition. This calculation result is compared to experimental measurements, and reasonable agreement has been reached.

Abstract: The electron behaviour in an Ar/CF4 inductively coupled plasma is investigated by a Langmuir probe and a hybrid model. The simulated and measured results include electron density, temperature and electron energy distribution function for different values of Ar/CF4 ratio, coil power and gas pressure. The hybrid plasma equipment model simulations show qualitative agreement with experiment. The effect of F2 electron attachment on the electron behaviour is explored by comparing two sets of data based on different F atom boundary conditions. It is demonstrated that electron attachment at F2 molecules is responsible for the depletion of low-energy electrons, causing a density decrease as well as a temperature increase when CF4 is added to an Ar plasma.

Abstract: A 2d3v particle-in-cell/Monte Carlo collisions model was applied to study the influence of the magnetic field strength on the cathode sheath region of a direct current (dc) magnetron discharge. When applying a magnetic field of 520-730 G, the cathode sheath width decreases with magnetic field strength, whereas, if a stronger magnetic field is applied (i. e. from 730 to 2600 G), the sheath width increases. This is explained by studying the structure of the sheath in different magnetic field strengths in terms of the electron and ion densities. The consequences of sheath structure on the sputter deposition process are also investigated. It is found that the magnetic field strength can control the erosion profile and the sputter rate.

Abstract: A hybrid model is developed for describing the effects of oxygen addition to argon glow discharges. The species taken into account in the model include Ar atoms in the ground state and the metastable level, O2 gas molecules in the ground state and two metastable levels, O atoms in the ground state and one metastable level, O3 molecules, Ar+, O+, O2+ and O− ions, as well as the electrons. The hybrid model consists of a Monte Carlo model for electrons and fluid models for the other plasma species. In total, 87 different reactions between the various plasma species are taken into account. Calculation results include the species densities and the importance of their production and loss processes, as well as the dissociation degree of oxygen. The effect of different O2 additions on these calculation results, as well as on the sputtering rates, is discussed.

Abstract: In traditional capacitively coupled plasmas, the discharge can be described by an electrostatic model, in which the Poisson equation is employed to determine the electrostatic electric field. However, current plasma reactors are much larger and driven at a much higher frequency. If the excitation wavelength k in the plasma becomes comparable to the electrode radius, and the plasma skin depth d becomes comparable to the electrode spacing, the electromagnetic (EM) effects will become significant and compromise the plasma uniformity. In this regime, capacitive discharges have to be described by an EM model, i.e., the full set of Maxwells equations should be solved to address the EM effects. This paper gives an overview of the theory, simulation and experiments that have recently been carried out to understand these effects, which cause major uniformity problems in plasma processing for microelectronics and flat panel display industries. Furthermore, some methods for improving the plasma uniformity are also described and compared.

Abstract: The purpose of this Brief Report is to point out the mistake in a formula for anisotropic electron scattering, previously published in Phys. Rev. A 41, 1112 (1990), which is widely used in Monte Carlo models of gas discharges. Anisotropic electron scattering is investigated based on the screened Coulomb potential between electrons and neutral atoms. The approach is also applied for electron scattering by nonpolar neutral molecules. Differential cross sections for electron scattering by Ar, N2, and CH4 are constructed on the basis of momentum and integrated cross sections. The formula derived in this paper is useful for Monte Carlo simulations of gas discharges.

Abstract: The electron bounce resonance heating (BRH) in dual-frequency capacitively coupled plasmas operated in oxygen is studied by different experimental methods and a particle-in-cell/Monte Carlo collision (PIC/MCC) simulation, and compared with the electropositive argon discharge. In comparison with argon, the experimental results show that in an oxygen discharge the resonance peaks in positive-ion density and light intensity tend to occur at larger electrode gaps. Moreover, at electrode gaps L > 2.5 cm, the positive-ion (and electron) density and the light emission drop monotonically in the oxygen discharge upon increasing L, whereas they rise (after an initial drop) in the argon case. At resonance gap the electronegativity reaches its maximum due to the BRH. All these experimental observations are explained by PIC/MCC simulations, which show that in the oxygen discharge the bulk electric field becomes quite strong and is out of phase with the sheath field. Therefore, it retards the resonance electrons when traversing the bulk, resulting in a suppressed BRH. Both experiment and simulation results show that this effect becomes more pronounced at lower high-frequency power, when the discharge mode changes from electropositive to electronegative. In a pure oxygen discharge, the BRH is suppressed with increasing pressure and almost diminishes at 12 Pa. Finally, the driving frequency significantly affects the BRH, because it determines the phase relation between bulk electric field and sheath electric field.