Abstract: A one-dimensional fluid model for a dielectric barrier discharge in methane, used as a chemical reactor for gas conversion, is developed. The model describes the gas phase chemistry governing the conversion process of methane to higher hydrocarbons. The spatially averaged densities of the various plasma species as a function of time are discussed. Besides, the conversion of methane and the yields of the reaction products as a function of the residence time in the reactor are shown and compared with experimental data. Higher hydrocarbons (C2Hy and C3Hy) and hydrogen gas are typically found to be important reaction products. Furthermore, the main underlying reaction pathways are determined.

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: 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.

Abstract: A fluid model is self-consistently established to investigate the harmonic effects in an inductively coupled plasma, where the electromagnetic field is solved by the finite difference time domain technique. The spatiotemporal distribution of harmonic current density, harmonic potential, and other plasma quantities, such as radio frequency power deposition, plasma density, and electron temperature, have been investigated. Distinct differences in current density have been observed when calculated with and without Lorentz force, which indicates that the nonlinear Lorentz force plays an important role in the harmonic effects, especially at low frequencies. Moreover, the even harmonics are larger than the odd harmonics both in the current density and the potential. Finally, the dependence of various plasma quantities with and without the Lorentz force on various driving frequencies is also examined. It is shown that the deposited power density decreases and the depth of penetration increases slightly because of the Lorentz force. The electron density increases distinctly while the electron temperature remains almost the same when the Lorentz force is taken into account.

Abstract: The interaction of thermal and hyperthermal Ni ions with gas-phase C60 fullerene was investigated at two temperatures with classical molecular dynamics simulations using a recently developed interatomic many-body potential. The interaction between Ni and C60 is characterized in terms of the NiC60 binding sites, complex formation, and the collision and temperature induced deformation of the C60 cage structure. The simulations show how ion implantation theoretically allows the synthesis of both endohedral Ni@C60 and exohedral NiC60 metallofullerene complexes.

Abstract: The destruction of ethylene in a dielectric barrier discharge plasma is investigated by the combination of kinetic modeling and experiments, as a case study for plasma-based gas purification. The influence of the specific energy deposition on the removal efficiency and the selectivity toward CO and CO2 is studied for different concentrations of ethylene. The model allows the identication of the destruction pathway in dry and humid air. The latter is found to be mainly initiated by metastable N2 molecules, but the further destruction steps are dominated by O atoms and OH radicals. Upon increasing air humidity, the removal efficiency drops by ±15% (from 85% to 70%), but the selectivity toward CO and CO2 stays more or less constant at 60% and 22%, respectively. Beside CO and CO2, we also identified acetylene, formaldehyde, and water as byproducts of the destruction process, with concentrations of 1606 ppm, 15033 ppm, and 185 ppm in humid air (with 20% RH), respectively. Finally, we investigated the byproducts generated by the humid air discharge itself, which are the greenhouse gases O3, N2O, and the toxic gas NO2.

Abstract: In this work particle-in-cell/Monte Carlo collision simulations are performed to study the heating mechanism and plasma characteristics in direct current (dc) superposed radio-frequency (RF) capacitively coupled plasmas, operated both in single-frequency (SF) and dual-frequency (DF) regimes. An RF (60/2 MHz) source is applied on the bottom electrode to sustain the discharge, and a dc source is fixed on the top electrode. The heating mechanism appears to be very different in dc superposed SF and DF discharges. When only a single source of 60 MHz is applied, the plasma bulk region is reduced by the dc source, thus the ionization rate and hence the electron density decrease with rising dc voltage. However, when a DF source of 60 and 2 MHz is applied, the electron density can increase upon addition of a dc voltage, depending on the gap length and applied dc voltage. This is explained from the spatiotemporal ionization rates in the DF discharge. In fact, a completely different behavior is observed for the ionization rate in the two half-periods of the LF source. In the first LF half-period, the situation resembles the dc superposed SF discharge, and the reduced plasma bulk region due to the negative dc bias results in a very small effective discharge area and a low ionization rate. On the other hand, in the second half-period, the negative dc bias is to some extent counteracted by the LF voltage, and the sheath close to the dc electrode becomes particularly thin. Consequently, the amplitude of the high-frequency sheath oscillations at the top electrode is largely enhanced, while the LF sheath at the bottom electrode is in its expanding phase and can thus well confine the high-energy electrons. Therefore, the ionization rate increases considerably in this second LF half-period. Furthermore, in addition to the comparison between SF and DF discharges and the effect of gap length and dc voltage, the effect of secondary electrons is examined.

Abstract: A computer model is developed for describing argon/nitrogen glow discharges. The species taken into account in the model include electrons, Ar atoms in the ground state and in the 4s metastable levels, N2 molecules in the ground state and in six different electronically excited levels, N atoms, Ar+ ions, N+, N2+, N3+ and N4+ ions. The fast electrons are simulated with a Monte Carlo model, whereas all other species are treated in a fluid model. 74 different chemical reactions are considered in the model. The calculation results include the densities of all the different plasma species, as well as information on their production and loss processes. The effect of different N2 additions, in the range between 0.1 and 10%, is investigated.