Abstract: Understanding the nature and effect of the multitude of plasma–surface interactions in plasma catalysis is a crucial requirement for further process development and improvement. A particularly intriguing and rather unique property of a plasma-catalytic setup is the ability of the plasma to modify the electronic structure, and hence chemical properties, of the catalyst through charging, i.e. the absorption of excess electrons. In this work, we develop a quantum chemical model based on density functional theory to study excess negative surface charges in a heterogeneous catalyst exposed to a plasma. This method is specifically applied to investigate plasma-catalytic CO2 activation on supported M/Al2O3 (M=Ti, Ni, Cu) single atom catalysts. We find that (1) the presence of a negative surface charge dramatically improves the reductive power of the catalyst, strongly promoting the splitting of CO2 to CO and oxygen, and (2) the relative activity of the investigated transition metals is also changed upon charging, suggesting that controlled surface charging is a powerful additional parameter to tune catalyst activity and selectivity. These results strongly point to plasma-induced surface charging of the catalyst as an important factor contributing to the plasma-catalyst synergistic effects frequently reported for plasma catalysis.

Abstract: Plasma catalysis is gaining increasing interest, but the underlying mechanisms are far from understood. Different catalyst materials will have different chemical effects, but in addition, they might also have different dielectric constants, which will affect surface charging, and thus the plasma behavior. In this work, we demonstrate that surface charging plays an important role in the streamer propagation and discharge enhancement inside catalyst pores, and in the plasma distribution along the dielectric surface, and this role greatly depends on the dielectric constant of the material. For εr50, surface charging causes the plasma to spread along the dielectric surface and inside the pores, leading to deeper plasma streamer penetration, while for εr>50 or for metallic coatings, the discharge is more localized, due to very weak surface charging. In addition, at εr=50, the significant surface charge density near the pore entrance causes a large potential drop at the sharp pore edges, which induces a strong electric field and results in most pronounced plasma enhancement near the pore entrance.

Abstract: In this work we study energy transport in a gliding arc discharge with two diverging flat
electrodes in argon gas at atmospheric pressure. The discharge is ignited at the shortest electrode
gap and it is pushed downstream by a forced gas flow. The current values considered are
relatively low and therefore a non-equilibrium plasma is produced. We consider two cases, i.e.
with high and low discharge current (28 mA and 2.8mA), and a constant gas flow of 10 lmin −1 ,
with a significant turbulent component to the velocity. The study presents an analysis of the
various energy transport mechanisms responsible for the redistribution of Joule heating to the
plasma species and the moving background gas. The objective of this work is to provide a
general understanding of the role of the different energy transport mechanisms in arc formation
and sustainment, which can be used to improve existing or new discharge designs. The work is
based on a three-dimensional numerical model, combining a fluid plasma model, the shear stress
transport Reynolds averaged Navier–Stokes turbulent gas flow model, and a model for gas
thermal balance. The obtained results show that at higher current the discharge is constricted
within a thin plasma column several hundred kelvin above room temperature, while in the low-
current discharge the combination of intense convective cooling and low Joule heating prevents
discharge contraction and the plasma column evolves to a static non-moving diffusive plasma,
continuously cooled by the flowing gas. As a result, the energy transport in the two cases is
determined by different mechanisms. At higher current and a constricted plasma column, the
plasma column is cooled mainly by turbulent transport, while at low current and an unconstricted
plasma, the major cooling mechanism is energy transport due to non-turbulent gas convection. In
general, the study also demonstrates the importance of turbulent energy transport in
redistributing the Joule heating in the arc and its significant role in arc cooling and the formation
of the gas temperature profile. In general, the turbulent energy transport lowers the average gas
temperature in the arc, thus allowing additional control of thermal non-equilibrium in the
discharge.

Abstract: Plasma catalysis is gaining increasing interest for various environmental applications. Catalytic
material can be inserted in different shapes in the plasma, e.g., as pellets, (coated) beads, but also
as honeycomb monolith and 3DFD structures, also called ‘structured catalysts’, which have high
mass and heat transfer properties. In this work, we examine the streamer discharge propagation
and the interaction between plasma and catalysts, inside the channels of such structured catalysts,
by means of a two-dimensional particle-in-cell/Monte Carlo collision model. Our results reveal
that plasma streamers behave differently in various structured catalysts. In case of a honeycomb
structure, the streamers are limited to only one channel, with low or high plasma density when
the channels are parallel or perpendicular to the electrodes, respectively. In contrast, in case of a
3DFD structure, the streamers can distribute to different channels, causing discharge
enhancement due to surface charging on the dielectric walls of the structured catalyst, and
especially giving rise to a broader plasma distribution. The latter should be beneficial for plasma
catalysis applications, as it allows a larger catalyst surface area to be exposed to the plasma.

Abstract: We present a zero-dimensional plasma kinetics model, including both surface and gas phase kinetics, to determine the role of vibrationally excited states in plasma-catalytic ammonia synthesis. We defined a new method to systematically capture the conditions of dielectric barrier discharges (DBDs), including those found in packed bed DBDs. We included the spatial and temporal nature of such discharges by special consideration of the number of micro-discharges in the model. We introduce a parameter that assigns only a part of the plasma power to the microdischarges, to scale the model conditions from filamentary to uniform plasma. Because of the spatial and temporal behaviour of the micro-discharges, not all micro-discharges occurring in the plasma reactor during a certain gas residence time are affecting the molecules. The fraction of power considered in the model ranges from 0.005 %, for filamentary plasma, to 100 %, for uniform plasma. If vibrational excitation is included in the plasma chemistry, these different conditions, however, yield an ammonia density that is only varying within one order of magnitude. At only 0.05 % of the power put into the uniform plasma component, a model neglecting vibrational excitation clearly does not result in adequate amounts of ammonia. Thus, our new model, which accounts for the concept in which not all the power is deposited by the micro-discharges, but some part may also be distributed in between them, suggests that vibrational kinetic processes are really important in (packed bed) DBDs. Indeed, vibrational excitation takes place in both the uniform plasma between the micro-discharges and in the strong micro-discharges, and is responsible for an increased N2 dissociation rate. This is shown here for plasma-catalytic ammonia synthesis, but might also be valid for other gas conversion applications.

Abstract: Plasma has been studied for several years to convert CO2 into value-added products. If CO2 could be converted in the presence of H2O as a cheap H-source for making syngas and oxygenates, it would mimic natural photosynthesis. However, CO2/H2O plasmas have not yet been extensively studied, not by experiments, and certainly not computationally. Therefore, we present here a kinetic modelling study to obtain a greater understanding of the vibrational kinetics of a CO2/H2O microwave plasma. For this purpose, we first created an electron impact cross section set for H2O, using a swarm-derived method. We added the new cross section set and CO2/H2O-related chemistry to a pure CO2 model. While it was expected that H2O addition mainly causes quenching of the CO2 asymmetric mode vibrational levels due to the additional CO2/H2O vibrational-translational relaxation, our model shows that the modifications in the vibrational kinetics are mainly induced by the strong electron dissociative attachment to H2O molecules, causing a reduction in electron density, and the corresponding changes in the input of energy into the CO2 vibrational levels by electron impact processes.

Abstract: In recent years there has been growing interest in the use of low-temperature atmospheric pressure plasmas for biomedical applications. Currently, however, there is very little fundamental knowledge regarding the relevant interaction mechanisms of plasma species with living cells. In this paper, we investigate the interaction of important plasma species, such as O3, O2 and O atoms, with bacterial peptidoglycan (or murein) by means of reactive molecular dynamics simulations. Specifically, we use the peptidoglycan structure to model the gram-positive bacterium Staphylococcus aureus murein. Peptidoglycan is the outer protective barrier in bacteria and can therefore interact directly with plasma species. Our results demonstrate that among the species mentioned above, O3 molecules and especially O atoms can break important bonds of the peptidoglycan structure (i.e. CO, CN and CC bonds), which subsequently leads to the destruction of the bacterial cell wall. This study is important for gaining a fundamental insight into the chemical damaging mechanisms of the bacterial peptidoglycan structure on the atomic scale.

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: Cold atmospheric pressure plasmas have proven to provide an alternative treatment of cancer by targeting tumorous cells while leaving their healthy counterparts unharmed. However, the underlying mechanisms of the plasma–cell interactions are not yet fully understood. Reactive oxygen species, and in particular hydroxyl radicals (OH), are known to play a crucial role in plasma driven apoptosis of
malignant cells. In this paper we investigate the interaction of OH radicals, as well as H2O2 molecules and HO2 radicals, with DNA by means of reactive molecular dynamics simulations using the ReaxFF force field. Our results provide atomic-scale insight into the dynamics of oxidative stress on DNA caused by the OH radicals, while H2O2 molecules appear not reactive within the considered timescale. Among the observed processes are the formation of 8-OH-adduct radicals, forming the first stages towards the formation of 8-oxoGua and 8-oxoAde, H-abstraction reactions of the amines, and the partial opening of loose DNA ends in aqueous solution.

Abstract: The plasma behavior in a parallel-plate dielectric barrier discharge (DBD) is simulated by a two-dimensional particle-in-cell/Monte Carlo collision model, comparing for the first time an unpacked (empty) DBD with a packed bed DBD, i.e., a DBD filled with dielectric spheres in the gas gap. The calculations are performed in air, at atmospheric pressure. The discharge is powered by a pulse with a voltage amplitude of −20 kV. When comparing the packed and unpacked DBD reactors with the same dielectric barriers, it is clear that the presence of the dielectric packing leads to a transition in discharge behavior from a combination of negative streamers and unlimited surface streamers on the bottom dielectric surface to a combination of predominant positive streamers and limited surface discharges on the dielectric surfaces of the beads and plates. Furthermore, in the packed bed DBD, the electric field is locally enhanced inside the dielectric material, near the contact points between the beads and the plates, and therefore also in the plasma between the packing beads and between a bead and the dielectric wall, leading to values of $4\times {10}

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 packed bed dielectric barrier discharge (DBD) is widely used for plasma catalysis applications, but the exact plasma characteristics in between the packing beads are far from understood. Therefore, we study here these plasma characteristics by means of fluid modelling and experimental observations using ICCD imaging, for packing materials with different dielectric constants. Our study reveals that a packed bed DBD reactor in dry air at atmospheric pressure may show three types of discharges, i.e. positive restrikes, filamentary microdischarges, which can also be localized between two packing beads, and surface discharges (so-called surface ionization
waves). Restrikes between the dielectric surfaces result in the formation of filamentary microdischarges, while surface charging creates electric field components parallel to the dielectric surfaces, leading to the formation of surface discharges. A transition in discharge mode occurs from surface discharges to local filamentary discharges between the packing beads when the dielectric constant of the packing rises from 5 to 1000. This may have implications for the efficiency of plasma catalytic gas treatment, because the catalyst activation may be limited by constraining the discharge to the contact points of the beads. The production of reactive species occurs most in the positive restrikes, the surface discharges and the local microdischarges in between the beads, and is less significant in the longer filamentary microdischarges. The faster streamer propagation and discharge development with higher dielectric constant of the packing beads leads to a faster production of reactive species. This study is of great interest for plasma catalysis, where packing beads with different dielectric constants are often used as supports for the catalytic materials. It allows us to better understand how different packing materials can influence the performance of packed bed plasma reactors for environmental applications.

Abstract: We propose a novel Power to Synthesis Gas (P2SG) approach, composed of two high-efficiency and renewable electricity-driven units, i.e., plasma catalytic reforming (PCR) and water electrolysis (WE), to produce high quality syngas from CH4, CO2 and H2O. As WE technology is already commercial, we mainly focus on the PCR unit, consisting of gliding arc plasma and Ni-based catalyst, for oxidative dry reforming of methane. An energy efficiency of 78.9% and energy cost of 1.0 kWh/Nm3 at a CH4 conversion of 99% and a CO2 conversion of 79% are obtained. Considering an energy efficiency of 80% for WE, the P2SG system yields an overall energy efficiency of 79.3% and energy cost of 1.8 kWh/Nm3. High-quality syngas is produced without the need for posttreatment units, featuring the ideal stoichiometric number of 2, with concentration of 94.6 vol%, and a desired CO2 fraction of 1.9 vol% for methanol synthesis. The PCR unit has the advantage of fast response to adapting to fluctuation of renewable electricity, avoiding local hot spots in the catalyst bed and coking, in contrast to conventional catalytic processes. Moreover, pure O2 from the WE unit is directly utilized by the PCR unit for oxidative dry reforming of methane, and thus, no air separation unit, like in conventional processes, is required. This work demonstrates the viability of the P2SG approach for large-scale energy storage of renewable electricity via electricity-to-fuel conversion.