We present the experimental evidence of the collisionless electron bounce resonance heating (BRH) in low-pressure dual-frequency capacitively coupled plasmas. In capacitively coupled plasmas at low pressures when the discharge frequency and gap satisfy a certain resonant condition, the high energy beamlike electrons can be generated by fast sheath expansion, and heated by the two sheaths coherently, thus the BRH occurs. By using a combined measurement of a floating double probe and optical emission spectroscopy, we demonstrate the effect of BRH on plasma properties, such as plasma density and light emission, especially in dual-frequency discharges.

We developed an implicit Particle-in-cell/Monte Carlo model in two-dimensional and axisymmetric geometry for the simulations of the radio-frequency discharges, by introducing several numerical schemes which include variable weights, multigrid field solver, etc. Compared to the standard explicit models, we found that the computational efficiency is significantly increased and the accuracy is still kept. Numerical schemes are discussed and benchmark results are shown. The code can be used to simulate practical reactors.

The electron heating mode transitions in capacitively coupled CF 4 discharges were studied by synergistically using two diagnostic methods in combination with Particle-in-Cell/Monte Carlo collision (PIC/MCC) simulations. Based on the method of phase resolved optical emission spectroscopy of trace rare gas, the spatiotemporal evolutions of energetic electrons were presented. The time-average electron density at the discharge center was measured by using a hairpin probe. All the experimental results were compared with those obtained from PIC/MCC simulations. Two different electron heating modes were observed depending on the discharge conditions: (1) the α mode (or electropositive mode), in which the electron heating maximum occurs near the sheath boundary, dominated by the sheath electric field during its expansion phase, (2) the drift-ambipolar (DA) mode (or electronegative mode), in which the electron heating maxima occur inside the entire bulk plasma and near the collapsing sheath edge, dominated by the drift field inside the bulk and the ambipolar fields near the collapsing sheath edge, respectively. The transitions between the two modes were presented when changing the rf power, working pressure and driving frequency.By increasing the power, the heating mode experiences a transition from DA to α mode. This is ascribed to the fact that at high powers, the sheath heating is enhanced, leading to a drastic decrease in the electronegativity, and consequently the DA electric field is significantly reduced. By increasing the pressure, a heating mode transition from a pure α mode, then a combination of α and DA modes, finally into a DA mode is induced. We found that the mode transition is much more sensitive to the change of working pressure than that of rf power. When increasing the pressure, there is an evident enhancement in the electron attachment, which can generate the negative ions and deplete the electrons, resulting in a higher electronegativity as well as a higher DA field, and therefore the excitation and ionization in the bulk are enhanced. The driving frequency is found to significantly affect the electronegativity, i.e. as the driving frequency increases, the discharge becomes more electropositive, and the sheath heating (α mode) dominates. Furthermore, we conclude that as the driving frequency is increased, the pressure, at which the mode transition occurs, is increased, while the power, at which the mode transition occurs, is decreased.

Recently a so-called electrical asymmetry effect (EAE), which could achieve high-degree separate control of ion flux and energy in dual-frequency capacitively coupled radio-frequency (CCRF) discharges, was discovered theoretically by Heil et al. [J. Phys. D: Appl. Phys. 41, 165202 (2008)] and was confirmed by experiments and theory/numerical simulations later on for electropositive argon discharges. In this work simulations based on particle-in-cell/Monte Carlo collision are performed to study the EAE on electronegative oxygen plasmas in geometrically symmetric CCRF discharges. Dual frequency discharges operating at 13.56 and 27.12 MHz are simulated for different pressures and the results are compared with those of electropositive argon discharges at the same conditions. It is found that in general the EAE on oxygen discharges has similar behavior as on argon discharge: The self-bias voltage η increases monotonically and almost linearly with the increase in the phase angle θ between the two driving voltages in the range 0<θ<90°, and the maximum ion energy varies by a factor of 3 by adjusting θ. However, the ion flux varies with θ by ±12% for low pressure and by ±15% for higher pressure, due primarily to an enhanced plasma series resonance, which then leads to dramatic changes in plasma density, power absorption and consequently the electronegativity. This may place a limitation for achieving separate control of ion energy and flux for electronegative plasma via the EAE.