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: 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: In this paper, an overview is given of modeling investigations carried out in our research group for a better understanding of plasmas used for medical, environmental and nano applications. The focus is both on modeling the plasma chemistry and the plasma-surface interactions. The plasma chemistry provides the densities and fluxes of the important plasma species. This information can be used as input when modeling the plasma-surface interactions. The combination of plasma simulations and plasma – surface interaction simulations provides a more comprehensive understanding of the underlying processes for these applications.

Abstract: In this paper we study two cold atmospheric pressure plasma jets, operating in Ar + 2% air, with a different electrode geometry but with the same power dissipated in the plasma. The density profiles of the biomedically active NO and O species throughout the plasma jet, previously obtained by laser diagnostics, are calculated by means of a zero-dimensional semi-empirical reaction kinetics model. A good agreement between the calculated and measured data is demonstrated. Furthermore, the most probable spatial power distribution in an RF driven plasma jet is obtained for the first time by comparing measured and calculated species density profiles. This was possible due to the strong effect of the power distribution on the NO and O density profiles. In addition the dominant reaction pathways for both the NO and the O species are identified. The model allows us to obtain key information on the reactive species production inside the jet, which is difficult to access by laser diagnostics in a coaxial geometry. Finally, we demonstrate that water impurities in the order of 100 ppm in the gas feed can have a significant effect on the spatial distribution of the NO and O density.

Abstract: The physical processes in an Ar/O2 magnetron discharge used for the reactive sputter deposition of TiOx thin films were simulated with a 2d3v particle-in-cell/Monte Carlo collisions (PIC/MCC) model. The plasma species taken into account are electrons, Ar+ ions, fast Arf atoms, metastable Arm* atoms, Ti+ ions, Ti atoms, O+ ions, O2+ ions, O− ions and O atoms. This model accounts for plasmatarget interactions, such as secondary electron emission and target sputtering, and the effects of target poisoning. Furthermore, the deposition process is described by an analytical surface model. The influence of the O2/Ar gas ratio on the plasma potential and on the species densities and fluxes is investigated. Among others, it is shown that a higher O2 pressure causes the region of positive plasma potential and the O− density to be more spread, and the latter to decrease. On the other hand, the deposition rates of Ti and O are not much affected by the O2/Ar proportion. Indeed, the predicted stoichiometry of the deposited TiOx film approaches x=2 for nearly all the investigated O2/Ar proportions.

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: 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: 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: Plasma catalysis holds great promise for environmental applications, provided that the process viability can be maximized in terms of energy efficiency and product selectivity. This requires a fundamental understanding of the various processes taking place and especially the mutual interactions between plasma and catalyst. In this review, we therefore first examine the various effects of the plasma on the catalyst and of the catalyst on the plasma that have been described in the literature. Most of these studies are purely experimental. The urgently needed fundamental understanding of the mechanisms underpinning plasma catalysis, however, may also be obtained through modelling and simulation. Therefore, we also provide here an overview of the modelling efforts that have been developed already, on both the atomistic and the macroscale, and we identify the data that can be obtained with these models to illustrate how modelling and simulation may contribute to this field. Last but not least, we also identify future modelling opportunities to obtain a more complete understanding of the various underlying plasma catalytic effects, which is needed to provide a comprehensive picture of plasma catalysis.

Abstract: In recent years there is increased interest in plasma-based CO2 conversion. Several plasma setups are being investigated for this purpose, but the most commonly used ones are a dielectric barrier discharge (DBD), a microwave (MW) plasma and a gliding arc (GA) reactor. In this proceedings paper, we will show results from our experiments in a (packed bed) DBD reactor and in a vortex-flow GA reactor, as well as from our model calculations for the detailed plasma chemistry in a DBD, MW and GA, for pure CO2 as well as mixtures of CO2 with N-2, CH4 and H2O.