Abstract: A combined Monte Carlo (MC)/analytical surface model is developed to study the plasma processes occurring during the reactive sputter deposition of TiOx thin films. This model describes the important plasma species with a MC approach (i.e. electrons, Ar+ ions, {\rm O}_2

Abstract: The conversion of methane to value-added chemicals and fuels is considered to be one of the challenges of the 21st century. In this paper we study, by means of fluid modeling, the conversion of methane to higher hydrocarbons or oxygenates by partial oxidation with CO2 or O2 in a dielectric barrier discharge. Sixty-nine different plasma species (electrons, ions, molecules, radicals) are included in the model, as well as a comprehensive set of chemical reactions. The calculation results presented in this paper include the conversion of the reactants and the yields of the reaction products as a function of residence time in the reactor, for different gas mixing ratios. Syngas (i.e. H2 + CO) and higher hydrocarbons (C2Hx) are typically found to be important reaction products.

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: In this paper, we investigate the etch process of so-called low-k organic material by means of a N2/H2 capacitively coupled plasma, as applied in the micro-electronics industry for the manufacturing of computer chips. In recent years, such an organic material has emerged as a possible alternative for replacing bulk SiO2 as a dielectric material in the back-end-of-line, because of the smaller parasitic capacity between adjacent conducting lines, and thus a faster propagation of the electrical signals throughout the chip. Numerical simulations with a hybrid plasma model, using an extensive plasma and surface chemistry set, as well as experiments are performed, focusing on the plasma properties as well as the actual etching process, to obtain a better insight into the underlying mechanisms. Furthermore, the effects of gas pressure, applied power and gas composition are investigated to try to optimize the etch process. In general, the plasma density reaches a maximum near the wafer edge due to the so-called 'edge effect'. As a result, the etch rate is not uniform but will also reach its maximum near the wafer edge. The pressure seems not to have a big effect. A higher power increases the etch rate, but the uniformity becomes (slightly) worse. The gas mixing ratio has no significant effect on the etch process, except when a pure H2 or N2 plasma is used, illustrating the synergistic effects of a N2/H2 plasma. In fact, our calculations reveal that the N2/H2 plasma entails an ion-enhanced etch process. The simulation results are in reasonable agreement with the experimental values. The microscopic etch profile shows the desired anisotropic shape under all conditions under study.

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: We present a zero-dimensional kinetic model of CO2 splitting in non-equilibrium plasmas. The model includes a description of the CO2 vibrational kinetics (25 vibrational levels up to the dissociation limit of the molecule), taking into account state-specific VT and VV relaxation reactions and the effect of vibrational excitation on other chemical reactions. The model is applied to study the reaction kinetics of CO2 splitting in an atmospheric-pressure dielectric barrier discharge (DBD) and in a moderate-pressure microwave discharge. The model results are in qualitative agreement with published experimental works. We show that the CO2 conversion and its energy efficiency are very different in these two types of discharges, which reflects the important dissociation mechanisms involved. In the microwave discharge, excitation of the vibrational levels promotes efficient dissociation when the specific energy input is higher than a critical value (2.0 eV/molecule under the conditions examined). The calculated energy efficiency of the process has a maximum of 23%. In the DBD, vibrationally excited levels do not contribute significantly to the dissociation of CO2 and the calculated energy efficiency of the process is much lower (5%).

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: In this work we have analyzed the properties of a gliding dc discharge in argon at atmospheric pressure. Despite the usual designation of these discharges as ‘gliding arc discharges’, it was found previously that they operate in two different regimes—glow and arc. Here we analyze the differences in both regimes by means of two dimensional fluid modeling. In order to address different aspects of the discharge operation, we use two models—Cartesian and axisymmetric in a cylindrical coordinate system. The obtained results show that the two types of discharges produce a similar plasma column for a similar discharge current. However, the different mechanisms of plasma channel attachment to the cathode could produce certain differences in the plasma parameters (i.e. arc elongation), and this can affect gas treatments applications.

Abstract: In this tutorial review paper, we illustrate how computer modeling can contribute to a better insight in inductively coupled plasma-mass spectrometry (ICP-MS). We start with a brief overview on previous efforts, studying the fundamentals of the ICP and ICP-MS, with main focus on previous modeling activities. Subsequently, we explain in detail the model that we developed in previous years, and we show typical calculation results, illustrating the plasma characteristics, gas flow patterns and the sample transport, evaporation and ionization. We also present the effect of various experimental parameters, such as operating conditions, geometrical aspects and sample characteristics, to illustrate how modeling can help to elucidate the optimal conditions for improved analytical performance.