Abstract: Thermal-catalytic gas processing is integral to many current industrial processes. Ever-increasing demands on conversion and energy efficiencies are a strong driving force for the development of alternative approaches. Similarly, synthesis of several functional materials (such as nanowires and nanotubes) demands special processing conditions. Plasma catalysis provides such an alternative, where the catalytic process is complemented by the use of plasmas that activate the source gas. This combination is often observed to result in a synergy between plasma and catalyst. This Review introduces the current state-of-the-art in plasma catalysis, including numerous examples where plasma catalysis has demonstrated its benefits or shows future potential, including CO2 conversion, hydrocarbon reforming, synthesis of nanomaterials, ammonia production, and abatement of toxic waste gases. The underlying mechanisms governing these applications, as resulting from the interaction between the plasma and the catalyst, render the process highly complex, and little is known about the factors leading to the often-observed synergy. This Review critically examines the catalytic mechanisms relevant to each specific application.

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: 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: In this work, the role of vibrationally excited HBr (HBr(vib)) is computationally investigated for a HBr/He inductively coupled plasma applied for Si etching. It is found that at least 50% of all dissociations of HBr occur through HBr(vib). This additional dissociation pathway through HBr(vib) makes the plasma significantly more atomic. It also results in a slightly higher electron temperature (i.e. about 0.2 eV higher compared to simulation results where HBr(vib) is not included), as well as a higher gas temperature (i.e. about 50 K higher than without including HBr(vib)), due to the enhanced Franck–Condon heating through HBr(vib) dissociation,
at the conditions investigated. Most importantly, the calculated etch rate with HBr(vib) included in the model is a factor 3 higher than in the case without HBr(vib), due to the higher fluxes of etching species (i.e. H and Br), while the chemical composition of the wafer surface shows no significant difference. Our calculations clearly show the importance of including HBr(vib) for accurate modeling of HBr-containing plasmas.

Abstract: Although CO2 conversion by plasma technology is gaining increasing interest, the
underlying mechanisms for an energy-efficient process are still far from understood. In this work, a reduced non-equilibrium CO2 plasma chemistry set, based on level lumping of the vibrational levels, is proposed and the reliability of this level-lumping method is tested by a self-consistent zero-dimensional code. A severe reduction of the number of equations to be solved is achieved, which is crucial to be able to model non-equilibrium CO2 plasmas by 2-dimensional models. Typical conditions of pressure and power used in a microwave plasma for CO2 conversion are investigated. Several different sets, using different numbers of lumped groups, are considered. The lumped models with 1, 2 or 3 groups are able to reproduce the gas temperature, electron density and electron temperature profiles, as calculated by the full model treating all individual excited levels, in the entire pressure range investigated. Furthermore, a 3-groups model is also able to reproduce the shape of the vibrational distribution function (VDF) and gives the most reliable prediction of the CO2 conversion. A strong influence of the vibrational excitation on the plasma characteristics is observed. Finally, the limitations of the lumped-levels method are discussed.

Abstract: Experiments show that the etch rate of Si with SF6 inductively coupled plasma (ICP) is significantly influenced by the absolute gas flow rate in the range of 50–600 sccm, with a maximum at around 200 sccm. Therefore, we numerically investigate the effects of the gas flow rate on the bulk plasma properties and on the etch rate, to obtain more insight in the underlying reasons of this effect. A hybrid Monte Carlo—fluid model is applied to simulate an SF6 ICP. It is found that the etch rate is influenced by two simultaneous effects: (i) the residence time of the gas and (ii) the temperature profile of the plasma in the ICP volume, resulting indeed in a maximum etch rate at 200 sccm.

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: CO 2 conversion by a gliding arc plasma is gaining increasing interest, but the underlying mechanisms for an energy-efficient process are still far from understood. Indeed, the chemical complexity of the non-equilibrium plasma poses a challenge for plasma modeling due to the huge computational load. In this paper, a one-dimensional (1D) gliding arc model is developed in a cylindrical frame, with a detailed non-equilibrium CO 2 plasma chemistry set, including the CO 2 vibrational kinetics up to the dissociation limit. The model solves a set of time- dependent continuity equations based on the chemical reactions, as well as the electron energy balance equation, and it assumes quasi-neutrality in the plasma. The loss of plasma species and heat due to convection by the transverse gas flow is accounted for by using a characteristic frequency of convective cooling, which depends on the gliding arc radius, the relative velocity of the gas flow with respect to the arc and on the arc elongation rate. The calculated values for plasma density and plasma temperature within this work are comparable with experimental data on gliding arc plasma reactors in the literature. Our calculation results indicate that excitation to the vibrational levels promotes efficient dissociation in the gliding arc, and this is consistent with experimental investigations of the gliding arc based CO 2 conversion in the literature. Additionally, the dissociation of CO 2 through collisions with O atoms has the largest contribution to CO 2 splitting under the conditions studied. In addition to the above results, we also demonstrate that lumping the CO 2 vibrational states can bring a significant reduction of the computational load. The latter opens up the way for 2D or 3D models with an accurate description of the CO 2 vibrational kinetics.

Abstract: Plasma catalysis is gaining increasing interest for various environmental applications, but the crucial question is whether plasma can be created inside catalyst pores and under which conditions. In practice, various catalytic support materials are used, with various dielectric constants. We investigate here the influence of the dielectric constant on the plasma properties inside catalyst pores and in the sheath in front of the pores, for various pore sizes. The calculations are performed by a two-dimensional fluid model for an atmospheric pressure dielectric barrier discharge in helium. The electron impact ionization rate, electron temperature, electron and ion density, as well as the potential distribution and surface charge density, are analyzed for a better understanding of the discharge behavior inside catalyst pores. The results indicate that, in a 100 μm pore, the electron impact ionization in the pore, which is characteristic for the plasma generation inside the pore, is greatly enhanced for dielectric constants below 300. Smaller pore sizes only yield enhanced ionization for smaller dielectric constants, i.e., up to εr = 200, 150, and 50 for pore sizes of 50, 30, and 10 μm. Thus, the most common catalyst supports, i.e., Al2O3 and SiO2, which have dielectric constants around εr = 8−11 and 4.2, respectively, should allow more easily that microdischarges can be formed inside catalyst pores, even for smaller pore sizes. On the other hand, ferroelectric materials with dielectric constants above 300 never seem to yield plasma enhancement inside catalyst pores, not even for 100 μm pore sizes. Furthermore, it is clear that the dielectric constant of the material has a large effect on the extent of plasma enhancement inside the catalyst pores, especially in the range between εr = 4 and εr = 200. The obtained results are explained in detail based on the surface charge density at the pore walls,
and the potential distribution and electron temperature inside and above the pores. The results obtained with this model are
important for plasma catalysis, as the production plasma species in catalyst pores might affect the catalyst properties, and thus
improve the applications of plasma catalysis.

Abstract: CO2 conversion into value-added chemicals is gaining increasing interest in recent years, and a gliding arc plasma has great potential for this purpose, because of its high energy efficiency. In this study, a chemical reaction kinetics model is presented to study the CO2 splitting in a gliding arc discharge. The calculated
conversion and energy efficiency are in good agreement with experimental data in a range of different operating conditions. Therefore, this reaction kinetics model can be used to elucidate the dominant chemical reactions contributing to CO2 destruction and formation. Based on this reaction pathway analysis, the restricting factors for CO2 conversion are figured out, i.e., the reverse reactions and the small treated gas fraction. This allows us to propose some solutions in order to improve the CO2 conversion, such as decreasing the gas temperature, by using a high frequency discharge, or increasing the power
density, by using a micro-scale gliding arc reactor, or by removing the reverse reactions, which could be realized in practice by adding possible scavengers for O atoms, such as CH4. Finally, we compare our results with other types of plasmas in terms of conversion and energy efficiency, and the results illustrate that gliding arc discharges are indeed quite promising for CO2 conversion, certainly when keeping in mind the possible solutions for further performance improvement.

Abstract: tWe investigated the plasma-assisted catalytic reactions for the production of value-added chemicalsfrom Ni-catalyzed plasma dry reforming of methane by means of density functional theory (DFT). Weinspected many activation barriers, from the early stage of adsorption of the major chemical fragmentsderived fromCH4andCO2molecules up to the formation of value-added chemicals at the surface, focusingon the formation of methanol, as well as the hydrogenation of C1and C2hydrocarbon fragments. Theactivation barrier calculations show that the presence of surface-bound H atoms and in some cases alsoremaining chemical fragments at the surface facilitates the formation of products. This implies that thehydrogenation of a chemical fragment on the hydrogenated crystalline surface is energetically favouredcompared to the simple hydrogenation of the chemical fragment at the bare Ni(111) surface. Indeed, thepresence of hydrogen modifies the electronic structure of the surface and the course of the reactions.We therefore conclude that surface-bound H atoms, and to some extent also the remaining chemicalfragments at the crystalline surface, induce the following effects: they facilitate associative desorption ofmethanol and ethane by increasing the rate of H-transfer to the adsorbed fragments while they impedehydrogenation of ethylene to ethane, thus promoting again the desorption of ethylene. Overall, they thusfacilitate the catalytic conversion of the formed fragments from CH4and CO2, into value-added chemicals.Finally, we believe that the retention of methane fragments, especially CH3, in the presence of surface-boundHatoms (as observed here for Ni) can be regarded as an identifier for the proper choice of a catalystfor the production of value-added chemicals.

Abstract: Here, we present routes to increase CO2 conversion into CO using an atmospheric pressure dielectric-barrier discharge. The change in conversion as a function of simple plasma parameters, such as power, flow rate, but also frequency, on-and-off power pulse, thickness and the chemical nature of the dielectric, wall and gas temperature, are described. By means of an in-depth electrical characterization of the discharge (effective plasma voltage, dielectric voltage, plasma current, number and lifetime of the microdischarges), combined with infrared analysis of the walls of the reactor, optical emission spectroscopy for the gas temperature, and mass spectrometry for the CO2 conversion, we propose a global interpretation of the effect of all the experimental parameters on the conversion and efficiency of the reaction.

Abstract: Cold atmospheric pressure plasmas in and in contact with liquids represent a growing field of research for various applications. Understanding the interactions between the plasma generated species and the liquid is crucial. In this work we perform molecular dynamics (MD) simulations based on a quantum mechanical method, i.e. density-functional based tight-binding (DFTB), to examine the interactions of OH radicals and O atoms in bulk water. Our calculations reveal that the transport of OH radicals through water is not only governed by diffusion, but also by an equilibrium reaction of H-abstraction with water molecules. Furthermore, when two OH radicals encounter each other, they either form a stable cluster, or react, resulting in the formation of a new water molecule and an O atom. In addition,
the O atoms form either oxywater (when in singlet configuration) or they remain stable in solution (when in triplet configuration), stressing the important role that O atoms can play in aqueous solution, and in contact with biomolecules. Our observations are in line with both experimental and ab initio results from the literature.

Abstract: Selective etching allows for obtaining carbon nanotubes with a specific chirality. While plasma-assisted etching has already been used to separate metallic tubes from their semiconducting counterparts, little is known about the nanoscale mechanisms of the etching process. We combine (reactive) molecular dynamics (MD) and force-bias Monte Carlo (tfMC) simulations to study H-etching of CNTs. In particular, during the hydrogenation and subsequent etching of both the carbon cap and the tube, they sequentially transform to different carbon nanostructures, including carbon nanosheet, nanowall, and polyyne chains, before they are completely removed from the surface of a substrate-bound Ni-nanocluster.We also found that onset of the etching process is different in the cases of the cap and the tube, although the overall etching scenario is similar in both cases. The entire hydrogenation/etching process for both cases is analysed in detail, comparing with available theoretical and experimental evidences.

Abstract: One of the most challenging and recurring problems when modeling plasmas is the lack of data on the key atomic and molecular reactions that drive plasma processes. Even when there are data for some reactions, complete and validated datasets of chemistries are rarely available. This hinders research on plasma processes and curbs development of industrial applications. The QDB project aims to address this problem by providing a platform for provision, exchange, and validation of chemistry datasets. A new data model developed for QDB is presented. QDB collates published data on both electron scattering and heavy-particle reactions. These data are formed into reaction sets, which are then validated against experimental data where possible. This process produces both complete chemistry sets and identifies key reactions that are currently unreported in the literature. Gaps in the datasets can be filled using established theoretical methods. Initial validated chemistry sets for SF6/CF4/O2 and SF6/CF4/N2/H2 are presented as examples.