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: A numerical parameter study has been performed for a cylindrical atmospheric pressure dielectric barrier discharge (DBD) in helium with nitrogen impurities using a two-dimensional time-dependent fluid model. The calculated electric currents and gap voltages as a function of time for a given applied potential are presented, as well as the number densities of the various plasma species. This study shows that for the geometry under consideration the applied voltage parameters have a large impact on the electric current profiles and that the discharge current is always determined by the electron and ion conduction currents while the displacement current is nearly negligible. A relative broadening of the current profiles (compared with the duration of the half cycle of the applied voltage) with an increase in the applied frequency is obtained. Nearly sinusoidal current wave forms, usually typical for radio frequency DBDs, are observed while still operating at the frequencies of tens of kilohertz. For the setup under investigation, the Townsend mode of the DBD is observed in the entire range of applied voltage amplitudes and frequencies. It is shown that the average power density dissipated in the discharge increases with rising applied voltage and frequency. An increase in applied voltage frequency leads to an increase in the electron density and a decrease in electron energy, while increasing the voltage amplitude has the opposite effect.

Abstract: A number of recent computational material design studies based on density functional theory (DFT) calculations have put forward a new class of materials with electrically switchable chemical characteristics that can be exploited in the development of tunable gas storage and electrocatalytic applications. We find systematic flaws in almost every computational study of gas adsorption on polarized or charged surfaces, stemming from an improper and unreproducible treatment of periodicity, leading to very large errors of up to 3 eV in some cases. Two simple corrective procedures that lead to consistent results are proposed, constituting a crucial course correction to the research in the field.

Abstract: A multi-species Monte Carlo (MC) model, combined with an analytical surface model, has been developed in order to investigate the general plasma processes occurring during the sputter deposition of complex oxide films in a dual-magnetron sputter deposition system. The important plasma species, such as electrons, Ar+ ions, fast Ar atoms and sputtered metal atoms (i.e. Mg and Al atoms) are described with the so-called multi-species MC model, whereas the deposition of MgxAlyOz films is treated by an analytical surface model. Targetsubstrate distances for both magnetrons in the dual-magnetron setup are varied for the purpose of growing stoichiometric complex oxide thin films. The metal atoms are sputtered from pure metallic targets, whereas the oxygen flux is only directed toward the substrate and is high enough to obtain fully oxidized thin films but low enough to avoid target poisoning. The calculations correspond to typical experimental conditions applied to grow these complex oxide films. In this paper, some calculation results are shown, such as the densities of various plasma species, their fluxes toward the targets and substrate, the deposition rates, as well as the film stoichiometry. Moreover, some results of the combined model are compared with experimental observations. Note that this is the first complete model, which can be applied for large and complicated magnetron reactor geometries, such as dual-magnetron configurations. With this model, we are able to describe all important plasma species as well as the deposition process. It can also be used to predict film stoichiometries of complex oxide films on the substrate.

Abstract: A Monte Carlo model has been developed for investigating the electron behavior in a dual-magnetron sputter deposition system. To describe the three-dimensional (3D) geometry, different reference frames, i.e. a local and a global coordinate system, were used. In this study, the influence of both closed and mirror magnetic field configurations on the plasma properties is investigated. In the case of a closed magnetic field configuration, the calculated electron trajectories show that if an electron is emitted in (or near) the center of the cathode, where the influence of the magnetic field is low, it is able to travel from one magnetron to the other. On the other hand, when an electron is created at the race track area, it is more or less trapped in the strong magnetic field and cannot easily escape to the second magnetron region. In the case of a mirror magnetic field configuration, irrespective of where the electron is emitted from the cathode, it cannot travel from one magnetron to the other because the magnetic field lines guide the electron to the substrate. Moreover, the electron density and electron impact ionization rate have been calculated and studied in detail for both configurations.

Abstract: A modelling study is performed of high-density low-pressure inductively coupled CF4/CHF3/H2/Cl2/O2/HBr plasmas under different gas mixing ratios. A reaction set describing the complete plasma chemistry is presented and discussed. The gas fraction of each component in this mixture is varied to investigate the sensitivity of the plasma properties, like electron density, plasma potential and species densities, towards the gas mixing ratios. This research is of great interest for microelectronics applications because these gases are often combined in two (or more)-component mixtures, and mixing gases or changing the fraction of a gas can sometimes yield unwanted reaction products or unexpected changes in the overall plasma properties due to the increased chemical complexity of the system. Increasing the CF4 fraction produces more F atoms for chemical etching as expected, but also more prominently lowers the density of Cl atoms, resulting in an actual drop in the etch rate under certain conditions. Furthermore, CF4 decreases the free electron density when mixed with Cl2. However, depending on the other gas components, CF4 gas can also sometimes enhance free electron density. This is the case when HBr is added to the mixture. The addition of H2 to the gas mixture will lower the sputtering process, not only due to the lower overall positive ion density at higher H2 fractions, but also because more H+, H2 + and H3 + are present and they have very low sputter yields. In contrast, a larger Cl2 fraction results in more chemical etching but also in less physical sputtering due to a smaller abundance of positive ions. Increasing the O2 fraction in the plasma will always lower the etch rate due to more oxidation of the wafer surface and due to a lower plasma density. However, it is also observed that the density of F atoms can actually increase with rising O2 gas fraction. This is relevant to note because the exact balance between fluorination and oxidation is important for fine-tuning the overall etch rate and for control of the sidewall profile. Finally, HBr is often used as a chemical etcher, but when mixed with F- or Cl-containing gases, HBr creates the same diluting effects as Ar or He, because a
higher fraction results in less chemical etching but more (physical) sputtering.

Abstract: A method has been developed for damage-free cesium (Cs) encapsulation within single-walled carbon nanotubes (SWNTs) with fine position selectivity. Precise energy tuning of Cs-ion irradiation revealed that there is a clear energy window (2060 eV) for the efficient encapsulation of Cs through the hexagonal network of SWNT sidewalls without causing significant damage. This minimum energy threshold of Cs-ion encapsulation (∼20 eV) matches well with the value obtained by ab initio simulation (∼22 eV). Furthermore, position-selective Cs encapsulation was carried out, resulting in the successful formation of pn-junction SWNT thin films with excellent environmental stability.

Abstract: A Langmuir probe and an ICCD are employed to study the discharge mode transition in Ar inductively coupled plasma. Electron density and plasma emission intensity are measured during the E (capacitive discharge) to H (inductive discharge) mode transitions at different pressures. It is found that plasma exists with a low electron density and a weak emission intensity in the E mode, while it has a high electron density and a strong emission intensity in the H mode. Meanwhile, the plasma emission intensity spatial (2D image) profile is symmetrical in the H mode, but the 2D image is an asymmetric profile in the E mode. Moreover, the electron density and emission intensity jump up discontinuously at high pressure, but increase almost continuously at the E to H mode transition under low pressure.

Abstract: A hybrid Monte Carlofluid model is applied to simulate the wafer-temperature-dependent etching of silicon with SF6 inductively coupled plasmas (ICP). The bulk plasma within the ICP reactor volume as well as the surface reactions occurring at the wafer are self-consistently described. The calculated etch rates are validated by experiments. The calculations and experiments are performed at two different wafer temperatures, i.e. 300 and 173 K, resembling conventional etching and cryoetching, respectively. In the case of cryoetching, a physisorbed SFx layer (x = 06) is formed on the wafer, which is negligible at room temperature, because of fast thermal desorption, However, even in the case of cryoetching, this layer can easily be disintegrated by low-energy ions, so it does not affect the etch rates. In the investigated pressure range of 19 Pa, the etch rate is always slightly higher at cryogenic conditions, both in the experiments and in the model, and this could be explained in the model due to a local cooling of the gas above the wafer, making the gas denser and increasing the flux of reactive neutrals, like F and F2, towards the wafer.

Abstract: A hybrid model, called the hybrid plasma equipment model, was used to study Ar/Cl(2) inductively coupled plasmas used for the etching of Si. The effects of substrate bias, source power and gas pressure on the plasma characteristics and on the fluxes and energies of plasma species bombarding the substrate were observed. A comparison with experimentally measured etch rates was made to investigate how the etch process is influenced and which plasma species mainly account for the etch process. First, the general plasma characteristics are investigated at the following operating conditions: 10% Ar 90% Cl(2) gas mixture, 5mTorr total gas pressure, 100 sccm gas flow rate, 250W source power, -200V dc bias at the substrate electrode and an operating frequency of 13.56MHz applied to the coil and to the substrate electrode. Subsequently, the pressure is varied from 5 to 80mTorr, the substrate bias from -100 to -300V and the source power from 250 to 1000W. Increasing the total gas pressure results in a decrease of the etch rate and a less anisotropic flux to the substrate due to more collisions of the ions in the sheath. Increasing the substrate bias has an effect on the energy of the ions bombarding the substrate and to a lesser extent on the magnitude of the ion flux. When source power is increased, it was found that, not the energy, but the magnitude of the ion flux is increased. The etch rate was more influenced by a variation of the substrate bias than by a variation of the source power, at these operating conditions. These results suggest that the etch process is mainly affected by the energy of the ions bombarding the substrate and the magnitude of the ion flux, and to a lesser extent by the magnitude of the radical flux.

Abstract: A hybrid model, called the hybrid plasma equipment model (HPEM), is used to describe the plasma chemistry in an inductively coupled plasma, operating in a gas mixture of C2H2 with either H2 or NH3, as typically used for carbon nanotube (CNT) growth. Two-dimensional profiles of power density, electron temperature and density, gas temperature, and densities of some plasma species are plotted and analyzed. Besides, the fluxes of the various plasma species towards the substrate (where the CNTs can be grown), as well as the decomposition rates of the feedstock gases (C2H2, NH3 and H2), are calculated as a function of the C2H2 fraction in both gas mixtures.

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: A gliding arc plasmatron (GAP), which is very promising for purification and gas conversion,
is characterized in nitrogen using optical emission spectroscopy and high-speed photography,
because the cross sections of electron impact excitation of N 2 are well known. The gas
temperature (of about 5500 K), the electron density (up to 1.5 × 10 15 cm −3 ) and the reduced
electric field (of about 37 Td) are determined using an absolutely calibrated intensified charge-
coupled device (ICCD) camera, equipped with an in-house made optical arrangement for
simultaneous two-wavelength diagnostics, adapted to the transient behavior of a GA channel
in turbulent gas flow. The intensities of nitrogen molecular emission bands, N 2 (C–B,0–0) as
well as N +
2 (B–X,0–0), are measured simultaneously. The electron density and the reduced
electric field are determined at a spatial resolution of 30 µm, using numerical simulation and
measured emission intensities, applying the Abel inversion of the ICCD images. The temporal
behavior of the GA plasma channel and the formation of plasma plumes are studied using a
high-speed camera. Based on the determined plasma parameters, we suggest that the plasma
plume formation is due to the magnetization of electrons in the plasma channel of the GAP by
an axial magnetic field in the plasma vortex.

Abstract: A gliding arc plasma is a potential way to convert CO2 into CO and O2, due to its non-equilibrium character, but little is known about the underlying mechanisms. In this paper, a self-consistent two-dimensional (2D) gliding arc model is developed, with a detailed non-equilibrium CO2 plasma chemistry, and validated with experiments. Our calculated values of the electron number density in the plasma, the CO2 conversion and energy efficiency show reasonable agreement with the experiments, indicating that the model can provide a realistic picture of the plasma chemistry. Comparison of the results with classical thermal conversion, as well as other plasma-based technologies for CO2 conversion reported in literature, demonstrates the non-equilibrium character of the gliding arc, and indicates that the gliding arc is a promising plasma reactor for CO2 conversion. However, some process modifications should be exploited to further improve its performance. As the model provides a realistic picture of the plasma behaviour, we use it first to investigate the plasma characteristics in a whole gliding arc cycle, which is necessary to understand the underlying mechanisms. Subsequently, we perform a chemical kinetics analysis, to investigate the different pathways for CO2 loss and formation. Based on the revealed discharge properties and the underlying CO2 plasma chemistry, the model allows us to propose solutions on how to further improve the
CO2 conversion and energy efficiency by a gliding arc plasma.

Abstract: A gliding arc (GA) plasma has great potential for CO2 conversion into value-added chemicals, because of its high energy efficiency. To improve the application, a 2D/3D fluid model is needed to investigate the CO2 conversion mechanisms in the actual discharge geometry. Therefore, the complex CO2 chemical kinetics description must be reduced due to the huge computational cost associated with 2D/3D models. This paper presents a chemistry reduction method for CO2 plasmas, based on the so-called directed relation graph method. Depending on the defined threshold values, some marginal species are identified. By means of a sensitivity analysis, we can further reduce the chemistry set by removing one by one the marginal species. Based on the socalled flux-sensitivity coupling, we obtain a reduced CO2 kinetics model, consisting of 36 or 15 species (depending on whether the 21 asymmetric mode vibrational states of CO2 are explicitly included or lumped into one group), which is applied to a GA discharge. The results are compared with those predicted with the full chemistry set, and very good agreement is reached. Moreover, the range of validity of the reduced CO2 chemistry set is checked, telling us that this reduced set is suitable for low power GA discharges. Finally, the time and spatial evolution of the CO2 plasma characteristics are presented, based on a 2D model with the reduced kinetics.