Abstract: In this work, a hybrid model is used to investigate the effect of different gas ratios on the Si etching and polymer film deposition characteristics in an Ar/CF4 inductively coupled plasma. The influence of the surface processes on the bulk plasma properties is studied, and also the spatial characteristics of important gas phase and etched species. The densities of F and CF2 decrease when the surface module is included in the simulations, due to the species consumption caused by etching and polymer deposition. The influence of the surface processes on the bulk plasma depends on the Ar/CF4 gas ratio. The deposited polymer becomes thicker at high CF4 content because of more abundant CFx radicals. As a result of the competition between the polymer thickness and the F flux, the etch rate first increases and then decreases upon increasing the CF4 content. The electron properties, more specifically the electron density profile, affect the Si etch characteristics substantially by determining the radical density and flux profiles. In fact, the radial profile of the etch rate is more uniform at low CF4 content since the electron density has a smooth distribution. At high CF4 content, the etch rate is less uniform with a minimum halfway along the wafer radius, because the electron density distribution is more localized. Therefore, our calculations predict that it is better to work at relatively high Ar/CF4 gas ratios, in order to obtain high etch rate and good profile uniformity for etch applications. This, in fact, corresponds to the typical experimental etch conditions in Ar/CF4 gas mixtures as found in the literature, where Ar is typically present at a much higher concentration than CF4.

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 paper we analyse the gas phase production and loss pathways for several biomedically active species, i.e. N2(A), O, O3, O2(a), N, H, HO2, OH, NO, NO2, N2O5, H2O2, HNO2 and HNO3, in an argon plasma jet flowing into an open humid air atmosphere. For this purpose, we employ a zero-dimensional reaction kinetics model to mimic the typical experimental conditions by fitting several parameters to experimentally measured values. These include ambient air diffusion, the gas temperature profile and power deposition along the jet effluent. We focus in detail on how the pathways of the biomedically active species change as a function of the position in the effluent, i.e. inside the discharge device, active plasma jet effluent and afterglow region far from the nozzle. Moreover, we demonstrate how the reaction kinetics and species production are affected by different ambient air humidities, total deposited power into the plasma and gas temperature along the jet. It is shown that the dominant pathways can drastically change as a function of the distance from the nozzle exit or experimental conditions.

Abstract: We use a zero-dimensional reaction kinetics model to simulate CO2 conversion in microwave discharges where the excitation of the vibrational levels plays a significant role in the dissociation kinetics. The model includes a description of the CO2 vibrational kinetics, taking into account state-specific VT and VV relaxation reactions and the effect of vibrational excitation on other chemical reactions. The model is used to simulate a general tubular microwave reactor, where a stream of CO2 flows through a plasma column generated by microwave radiation. We study the effects of the internal plasma parameters, namely the reduced electric field, electron density and the total specific energy input, on the CO2 conversion and its energy efficiency. We report the highest energy efficiency (up to 30%) for a specific energy input in the range 0.41.0 eV/molecule and a reduced electric field in the range 50100 Td and for high values of the electron density (an ionization degree greater than 10−5). The energy efficiency is mainly limited by the VT relaxation which contributes dominantly to the vibrational energy losses and also contributes significantly to the heating of the reacting gas. The model analysis provides useful insight into the potential and limitations of CO2 conversion in microwave discharges.

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: A packed bed dielectric barrier discharge plasma reactor is computationally studied with a fluid model. Two different complementary axisymmetric 2D geometries are used to mimic the intrinsic 3D problem. It is found that a packing enhances the electric field strength and electron temperature at the contact points of the dielectric material due to polarization of the beads by the applied potential. As a result, these contact points prove to be of direct importance to initiate the plasma. At low applied potential, the discharge stays at the contact points, and shows the properties of a Townsend discharge. When a high enough potential is applied, the plasma will be able to travel through the gaps in between the beads from wall to wall, forming a kind of glow discharge. Therefore, the inclusion of a so-called ‘channel of voids’ is indispensable in any type of packed bed modelling.

Abstract: In this experimental study, a flowing dielectric barrier discharge operating at atmospheric pressure is used for the splitting of CO2 into O2 and CO. The influence of the applied frequency and plasma power on the microdischarge properties is investigated to understand their role on the CO2 conversion. Electrical measurements are carried out to explain the conversion trends and to characterize the microdischarges through their number, their lifetime,
their intensity and the induced electrical charge. Their influence on the gas and electrode temperatures is also evidenced through optical emission spectroscopy and infrared imaging. It is shown that, in our configuration, the conversion depends mostly on the charge delivered in the plasma and not on the effective plasma voltage when the applied power is modified. Similarly, at constant total current, a better conversion is observed at low frequencies, where a less filamentary discharge regime with a higher effective plasma voltage than that at a higher
frequency is obtained.

Abstract: Dielectric barrier discharges (DBDs) are commonly used to generate cold plasmas at
atmospheric pressure. Whatever their configuration (tubular or planar), the presence of a dielectric barrier is mandatory to prevent too much charge build up in the plasma and the formation of a thermal arc. In this article, the role of the barrier thickness (2.0, 2.4 and 2.8 mm) and of the kind of dielectric material (alumina, mullite, pyrex, quartz) is investigated on the filamentary behavior in the plasma and on the CO2 conversion in a tubular flowing DBD, by means of mass spectrometry measurements correlated with electrical characterization and IR imaging. Increasing the barrier thickness decreases the capacitance, while preserving the electrical charge. As a result, the voltage over the dielectric increases and a larger number of microdischarges is generated, which enhances the CO2 conversion. Furthermore, changing the dielectric material of the barrier, while keeping the same geometry and dimensions, also affects the CO2 conversion. The highest CO2 conversion and energy efficiency are obtained for quartz and alumina, thus not following the trend of the relative permittivity. From the
electrical characterization, we clearly demonstrate that the most important parameters are the somewhat higher effective plasma voltage (yielding a somewhat higher electric field and electron energy in the plasma) for quartz, as well as the higher plasma current (and thus larger electron density) and the larger number of microdischarge filaments (mainly for alumina, but also for quartz). The latter could be correlated to the higher surface roughness for alumina and to the higher voltage over the dielectric for quartz.

Abstract: CO2 conversion into value-added products has gained significant interest over the few last years, as the greenhouse gas concentrations constantly increase due to anthropogenic activities. Here we report on experiments for CO2 conversion by means of a cold atmospheric plasma using a cylindrical flowing dielectric barrier discharge (DBD) reactor. A detailed comparison of this DBD ignited in a so-called burst mode (i.e. where an AC voltage is applied during a limited amount of time) and pure AC mode is carried out to evaluate their effect on the conversion of CO2 as well as on the energy efficiency. Decreasing the duty cycle in the burst mode from 100% (i.e. corresponding to pure AC mode) to 40% leads to a rise in the
conversion from 16–26% and to a rise in the energy efficiency from 15 to 23%. Based on a detailed electrical analysis, we show that the conversion correlates with the features of the microfilaments. Moreover, the root-mean-square voltage in the burst mode remains constant as a function of the process time for the duty cycles <70%, while a higher duty cycle or the usual pure AC mode leads to a clear voltage decay by more than 500 V, over approximately 90 s, before reaching a steady state regime. The higher plasma voltage in the burst mode yields a higher electric field. This causes the increasing the electron energy, and therefore their
involvement in the CO2 dissociation process, which is an additional explanation for the higher CO2 conversion and energy efficiency in the burst mode.

Abstract: Electrical breakdown by the application of an electric field occurs more easily in hot gases than in cold gases because of the extra electron-species interactions that occur as a result of dissociation, ionization and excitation at higher temperature. This paper discusses some overlooked physics and clarifies inaccuracies in the evaluation of the effective ionization coefficients and the critical reduced breakdown electric field of CO2 at elevated temperature, considering the influence of excited states and ion kinetics. The critical reduced breakdown electric field is obtained by balancing electron generation and loss mechanisms using the electron energy distribution function (EEDF) derived from the Boltzmann transport equation under the two-term approximation. The equilibrium compositions of the hot gas mixtures are determined based on Gibbs free energy minimization considering the ground states as well as vibrationally and electronically excited states as independent species, which follow a Boltzmann distribution with a fixed excitation temperature. The interaction cross sections between electrons and the excited species, not reported previously, are properly taken into account. Furthermore, the ion kinetics, including electron–ion recombination, associative electron detachment, charge transfer and ion conversion into stable negative ion clusters, are also considered. Our results indicate that the excited species lead to a greater population of high-energy electrons at higher gas temperature and this affects the Townsend rate coefficients (i.e. of electron impact ionization and attachment), but the critical reduced breakdown electric field strength of CO2 is only affected when also properly accounting for the ion kinetics. Indeed, the latter greatly influences the effective ionization coefficients and hence the critical reduced breakdown electric field at temperatures above 1500 K. The rapid increase of the dissociative electron attachment cross-section of molecular oxygen with rising vibrational quantum number leads to a larger electron loss rate and this enhances the critical reduced breakdown electric field strength in the temperature range where the concentration of molecular oxygen is relatively high. The results obtained in this work show reasonable agreement with experimental results from literature, and are important for the evaluation of the dielectric strength of CO2 in a highly reactive environment at elevated temperature.

Abstract: We present a 3D and 2D Cartesian quasi-neutral plasma model for a low current argon gliding arc discharge, including strong interactions between the gas flow and arc plasma column.
The 3D model is applied only for a short time of 0.2 ms due to its huge computational cost. It mainly serves to verify the reliability of the 2D model. As the results in 2D compare well with those in 3D, they can be used for a better understanding of the gliding arc basic characteristics. More specifically, we investigate the back-breakdown phenomenon induced by an artificially controlled plasma channel, and we discuss its effect on the gliding arc characteristics. The
back-breakdown phenomenon, or backward-jump motion of the arc, as observed in the experiments, results in a drop of the gas temperature, as well as in a delay of the arc velocity with respect to the gas flow velocity, allowing more gas to pass through the arc, and thus increasing the efficiency of the gliding arc for gas treatment applications.

Abstract: In this study we quantitatively investigate for the first time the plasma characteristics of an argon gliding arc with a 3D model. The model is validated by comparison with available experimental data from literature and a reasonable agreement is obtained for the calculated gas temperature and electron density. A complete arc cycle is modeled from initial ignition to arc decay. We investigate how the plasma characteristics, i.e., the electron temperature, gas temperature,
reduced electric field, and the densities of electrons, Ar+ and Ar2+ ions and Ar(4s) excited states, vary over one complete arc cycle, including their behavior in the discharge and post-discharge. These plasma characteristics exhibit a different evolution over one arc cycle, indicating that either the active discharge stage or the post-discharge stage can be beneficial for certain applications.

Abstract: The formation process of a microdischarge (MD) in both μm- and nm-sized catalyst pores is simulated by a two-dimensional particle-in-cell/Monte Carlo collision model. A parallel-plate dielectric barrier discharge configuration in filamentary mode is considered in ambient air. The discharge is powered by a high voltage pulse. Our calculations reveal that a streamer can penetrate into the surface features of a porous catalyst and MDs can be formed inside both μm- and nm-sized pores, yielding ionization inside the pore. For the μm-sized pores, the ionization mainly occurs inside the pore, while for the nm-sized pores the ionization is strongest near and inside the pore. Thus, enhanced discharges near and inside the mesoporous catalyst are observed. Indeed, the maximum values of the electric field, ionization rate and electron density occur near and inside the pore. The maximum electric field and electron density inside the pore first increase when the pore size rises from 4 nm to 10 nm, and then they decrease for the 100 nm pore, due to
a more pronounced surface discharge for the smaller pores. However, the ionization rate is highest for the 100 nm pore due to the largest effective ionization region.

Abstract: Plasma–liquid interactions represent a growing interdisciplinary area of research involving plasma science, fluid dynamics, heat and mass transfer, photolysis, multiphase chemistry and aerosol science. This review provides an assessment of the state-of-the-art of this multidisciplinary area and identifies the key research challenges. The developments in diagnostics, modeling and further extensions of cross section and reaction rate databases that are necessary to address these challenges are discussed. The review focusses on nonequilibrium plasmas.

Abstract: Low-temperature plasma chemical kinetic models are particularly important to the plasma community. These models typically require dozens of inputs, especially rate coefficients. The latter are not always precisely known and it is not surprising that the error on the rate coefficient data can propagate to the model output. In this paper, we present a model that uses N = 400 different combinations of rate coefficients based on the uncertainty attributed to each rate coefficient, giving a good estimation of the uncertainty on the model output due to the rate coefficients. We demonstrate that the uncertainty varies a lot with the conditions and the type of output. Relatively low uncertainties (about 15%) are found for electron density and temperature, while the uncertainty can reach more than an order of magnitude for the population of the vibrational levels in some cases and it can rise up to 100% for the CO2 conversion. The reactions that are mostly responsible for the largest uncertainties are identified. We show that the conditions of pressure, gas temperature and power density have a great effect on the uncertainty and on which reactions lead to this uncertainty. In all the cases tested here, while the absolute values may suffer from large uncertainties, the trends observed in previous modeling work are still valid. Finally, in accordance with the work of Turner, a number of ‘good practices’ is recommended.