Abstract: Superimposition of dual frequencies (DFs) is one of the methods used for controlling plasma distribution in an inductively coupled plasma (ICP) source. The effects of a superimposed DF on the argon plasma characteristics have been investigated using a two-dimensional self-consistent fluid model. When both currents are fixed at 6A, the plasma density drops with decrease in one of the source frequencies due to less efficient heating and the plasma uniformity improves significantly. Moreover, for ICP operated with superimposed DFs (i.e., 4.52MHz/13.56MHz and 2.26MHz/13.56MHz), the current source exhibits the same period as the low frequency (LF) component, and the plasma density is higher than that obtained at a single frequency (i.e., 4.52 and 2.26MHz) with the same total current of 12A. However, at superimposed current frequencies of 6.78MHz/13.56MHz, the plasma density is lower than that obtained at a single frequency of 6.78MHz due to the weaker negative azimuthal electric field between two positive maxima during one period of 6.78MHz. When the superimposed DF ICP operates at 2.26 and 13.56MHz, the rapid oscillations of the induced electric field become weaker during one period of 2.26MHz as the current ratio of 2.26MHz/13.56MHz rises from 24A/7 A to 30A/1 A, and the plasma density drops with the current ratio due to weakened electron heating. The uniformity of plasma increases due to sufficient diffusion under the low-density condition.

Abstract: A systematic investigation of voltage effect on the plasma parameters, especially the species densities and chemical reaction mechanisms, in the capacitive N-2 discharges is performed by employing a two-dimensional self-consistent fluid model. The validity of the numerical model is first demonstrated by the qualitative agreement of the calculated and experimental results. Then, the densities, production mechanisms, and loss mechanisms of species from simulation are examined at various voltages. It is found that all the species densities increase monotonically with the voltage, whereas their spatial profiles at lower voltages are quite different from those at higher voltages. The electrons and Nthorn 2 ions are mainly generated by the electron impact ionization of N-2 gas, while the Nthorn ions, whose density is one or two orders of magnitude lower, are mostly formed by the ionization of N atoms. The electron impact dissociation of N-2 gas dominates the generation of N atoms, which are mostly destroyed for the Nthorn ion production. As for the excited N-2 levels, the level conversion processes play a very important role in their production and depletion mechanisms, except for the electron impact excitation of the ground state N-2 molecules. Published under license by AIP Publishing.

Abstract: Self-organized striated structure has been observed experimentally and numerically in CF4 plasmas in radio-frequency capacitively coupled plasmas recently (Liu et al 2016 Phys. Rev. Lett. 116 255002). In this work, the striated structure is investigated in a capacitively coupled oxygen discharge with the introduction of the effect from the secondary electron emission, based on a particle-in-cell/Monte Carlo collision model. As we know, the transport of positive and negative ions plays a key role in the formation of striations in electronegative gases, for which, the electronegativity needs to be large enough. As the secondary electron emission increases, electrons in the sheaths gradually contribute more ionization to the discharge. Meanwhile, the increase of the electron density, especially in the plasma bulk, leads to an increased electrical conductivity and a reduced bulk electric field, which would shield the ions' mobility. These changes result in enlarged striation gaps. And then, with more emitted electrons, obvious disruption of the striations is observed accompanied with a transition of electron heating mode. Due to the weakened field, the impact ionization in the plasma bulk is attenuated, compared with the enhanced ionization caused by secondary electrons. This would lead to the electron heating mode transition from striated (STR) mode to gamma-mode. Besides, our investigation further reveals that gamma-mode is more likely to dominate the discharge under high gas pressures or driving voltages.

Abstract: Plasma catalysis is gaining increasing interest for various environmental applications, but the crucial question is whether plasma can be created inside catalyst pores and under which conditions. In practice, various catalytic support materials are used, with various dielectric constants. We investigate here the influence of the dielectric constant on the plasma properties inside catalyst pores and in the sheath in front of the pores, for various pore sizes. The calculations are performed by a two-dimensional fluid model for an atmospheric pressure dielectric barrier discharge in helium. The electron impact ionization rate, electron temperature, electron and ion density, as well as the potential distribution and surface charge density, are analyzed for a better understanding of the discharge behavior inside catalyst pores. The results indicate that, in a 100 μm pore, the electron impact ionization in the pore, which is characteristic for the plasma generation inside the pore, is greatly enhanced for dielectric constants below 300. Smaller pore sizes only yield enhanced ionization for smaller dielectric constants, i.e., up to εr = 200, 150, and 50 for pore sizes of 50, 30, and 10 μm. Thus, the most common catalyst supports, i.e., Al2O3 and SiO2, which have dielectric constants around εr = 8−11 and 4.2, respectively, should allow more easily that microdischarges can be formed inside catalyst pores, even for smaller pore sizes. On the other hand, ferroelectric materials with dielectric constants above 300 never seem to yield plasma enhancement inside catalyst pores, not even for 100 μm pore sizes. Furthermore, it is clear that the dielectric constant of the material has a large effect on the extent of plasma enhancement inside the catalyst pores, especially in the range between εr = 4 and εr = 200. The obtained results are explained in detail based on the surface charge density at the pore walls,
and the potential distribution and electron temperature inside and above the pores. The results obtained with this model are
important for plasma catalysis, as the production plasma species in catalyst pores might affect the catalyst properties, and thus
improve the applications of plasma catalysis.

Abstract: Computer simulations are performed for an argon inductively coupled plasma (ICP) with a capacitive radio-frequency bias power, to investigate the bias effect on the discharge mode transition and on the plasma characteristics at various ICP currents, bias voltages, and bias frequencies. When the bias frequency is fixed at 13.56 MHz and the ICP current is low, e.g., 6A, the spatiotemporal averaged plasma density increases monotonically with bias voltage, and the bias effect is already prominent at a bias voltage of 90 V. The maximum of the ionization rate moves toward the bottom electrode, which indicates clearly the discharge mode transition in inductive/capacitive discharges. At higher ICP currents, i.e., 11 and 13 A, the plasma density decreases first and then increases with bias voltage, due to the competing mechanisms between the ion acceleration power dissipation and the capacitive power deposition. At 11 A, the bias effect is still important, but it is noticeable only at higher bias voltages. At 13 A, the ionization rate is characterized by a maximum at the reactor center near the dielectric window at all selected bias voltages, which indicates that the ICP power, instead of the bias power, plays a dominant role under this condition, and no mode transition is observed. Indeed, the ratio of the bias power to the total power is lower than 0.4 over a wide range of bias voltages, i.e., 0300V. Besides the effect of ICP current, also the effect of various bias frequencies is investigated. It is found that the modulation of the bias power to the spatiotemporal distributions of the ionization rate at 2MHz is strikingly different from the behavior observed at higher bias frequencies. Furthermore, the minimum of the plasma density appears at different bias voltages, i.e., 120V at 2MHz and 90V at 27.12 MHz.

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: A two-dimensional self-consistent fluid model combined with the full set of Maxwell equations is employed to investigate an Ar/CF4 capacitively coupled plasma, focusing on the phase-shift effect on the plasma characteristics at various frequencies and gas mixture ratios. When the discharge is sustained by a single frequency at 13.56 MHz in an Ar/CF4 mixture with a ratio of 0.9/0.1, no obvious difference is detected between the electron densities obtained in the so-called electrostatic model (with only the static electric fields taken into account) and the electromagnetic model (which includes the electromagnetic effects). However, as the frequency increases to 60 and 100 MHz, the difference becomes distinct, due to the significant influence of the electromagnetic effects. The phase-shift effect on the plasma radial uniformity has also been investigated in a dual frequency discharge, i.e. when the top driven source is switched on with a phase difference phiv ranging from 0 to π, in the frequency range 13.56100 MHz. At low concentration of CF4 (10%), Ar+ ions are the major positive ions in the entire range of frequencies. When the frequency is low, i.e. 13.56 MHz, the Ar+ density exhibits an off-axis peak at phiv = 0 due to the edge effect, and a better uniformity caused by the phase-shift modulation is obtained at phiv = π. At 60 MHz, the Ar+ density varies from edge-peaked at phiv = 0 to uniform (i.e. at phiv = 0.53π), and finally at phiv = π, a broad maximum is observed at the centre due to the standing-wave effect. As the frequency increases to 100 MHz, the best radial uniformity is reached at 0.25π, and the maximum moves again towards the radial wall in the reverse-phase case (phiv = π) due to the dominant skin effect. When the frequency is fixed at 100 MHz, the phase-shift control shows a different behaviour at a high concentration of CF4. For instance, the ${\rm CF}_3

Abstract: A two-dimensional fluid model, including the full set of Maxwell equations, has been developed and applied to investigate the effect of a phase shift between two power sources on the radial uniformity of several plasma characteristics in a hydrogen capacitively coupled plasma. This study was carried out at various frequencies in the range 13.56200 MHz. When the frequency is low, at 13.56 MHz, the plasma density is characterized by an off-axis peak when both power sources are in-phase (phgr = 0), and the best radial uniformity is obtained at phgr = π. This trend can be explained because the radial nonuniformity caused by the electrostatic edge effect can be effectively suppressed by the phase-shift effect at a phase difference equal to π. When the frequency rises to 60 MHz, the plasma density profiles shift smoothly from edge-peaked over uniform to centre-peaked as the phase difference increases, due to the pronounced standing-wave effect, and the best radial uniformity is reached at phgr = 0.3π. At a frequency of 100 MHz, a similar behaviour is observed, except that the maximum of the plasma density moves again towards the radial edge at the reverse-phase case (phgr = π), because of the dominant skin effect. When the frequency is 200 MHz, the bulk plasma density increases significantly with increasing phase-shift values, and a better uniformity is obtained at phgr = 0.4π. This is because the density in the centre increases faster than at the radial edge as the phase difference rises, due to the increasing power deposition Pz in the centre and the decreasing power density Pr at the radial edge. As the phase difference increases to π, the maximum near the radial edge becomes obvious again. This is because the skin effect has a predominant influence on the plasma density under this condition, resulting in a high density at the radial edge. Moreover, the axial ion flux increases monotonically with phase difference, and exhibits similar profiles to the plasma density. The calculation results illustrate that the radial uniformity of the various plasma characteristics is strongly dependent on the applied frequency and the phase shift between both power sources, which is important to realize, for controlling the uniformity of the plasma etch and deposition processes.