Abstract: Dielectric barrier discharges (DBDs) typically operate in the filamentary regime and thus exhibit great spatial and temporal non-uniformity. In order to optimize DBDs for various applications, such as in plasma catalysis, more fundamental insight is needed. Here, we consider how the millions of microdischarges, characteristic for a DBD, influence individual gas molecules. We use a Monte Carlo approach to determine the number of microdischarges to which a single molecule would be exposed, by means of particle tracing simulations through a full-scale packed bed DBD reactor, as well as an empty DBD reactor. We find that the fraction of microdischarges to which the molecules are exposed can be approximated as the microdischarge volume over the entire reactor gas volume. The use of this concept provides good agreement between a plasma-catalytic kinetics model and experiments for plasma-catalytic NH3 synthesis. We also show that the concept of the fraction of microdischarges indicates the efficiency by which the plasma power is transferred to the gas molecules. This generalised concept is also applicable for other spatially and temporally non-uniform plasmas.

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: Ablation-controlled plasmas have been used in a range of technical applications where local thermodynamic equilibrium (LTE) is often violated near the wall due to the strong cooling effect caused by the ablation of wall materials. The thermodynamic and transport properties of ablated polytetrafluoroethylene (PTFE) vapor, which determine the flowing plasma behavior in such applications, are calculated based on a two-temperature model at atmospheric pressure. To our knowledge, no data for PTFE have been reported in the literature. The species composition and thermodynamic properties are numerically determined using the two-temperature Saha equation and the Guldberg-Waage equation according to van de Sanden et al's derivation. The transport coefficients, including viscosity, thermal conductivity and electrical conductivity, are calculated with the most recent collision interaction potentials using Devoto's electron and heavy-particle decoupling approach but expanded to the third-order approximation (second-order for viscosity) in the frame of the Chapman-Enskog method. Results are computed for different degrees of thermal non-equilibrium, i.e. the ratio of electron to heavy-particle temperatures, from 1 to 10, with electron temperature ranging from 300 to 40 000 K. Plasma transport properties in the LTE state obtained from the present work are compared with existing published results and the causes for the discrepancy analyzed. The two-temperature plasma properties calculated in the present work enable the modeling of wall ablation-controlled plasma processes.

Abstract: Packed bed dielectric barrier discharge (PB-DBD) plasma reactors are very promising for various plasma catalysis applications, but the exact mechanisms of plasma-catalyst interaction are far from understood, because the plasma discharge and catalyst/packing properties are mutually dependent. To better understand the effect of packing dielectric material on the electrical plasma properties, we study here a single bead DBD plasma reactor operating in dry air, with beads of different dielectric constant and for different applied voltages, by means of fluid modelling validated by optical imaging experiments. Our study reveals that the plasma in the single bead DBD reactor can manifest itself in two different modalities, i.e. (a) polar discharges at the bead poles in contact with the electrodes, and (b) a streamer discharge caused by surface ionization waves, which bridges the gas gap. Beads with high dielectric constant result in localised electric field enhancement and hence yield a reduction of the applied voltage required for plasma production. At low applied voltage, the discharge appears as polar discharges between the bead and the electrodes, and upon higher voltage it undergoes a transition into a bridging streamer discharge. The transition voltage to the streamer mode rises for beads with higher dielectric constant. These observations are important for plasma catalysis applications. A higher dielectric constant yields a higher electric field and thus higher average electron energy and density, giving rise to more reactive species, but it also yields a confined discharge near the contact points of packing beads, limiting the interaction area between the catalyst and the active plasma species. In addition, our model reveals that the dielectric bead behaves as a capacitor and traps charges, which can explain the significant occurrence of partial discharging in PB-DBDs and non-parallelogram shaped Lissajous plots. Hence, equivalent circuit modelling of PB-DBDs should take into account the role of packing beads in charge trapping as a capacitor.

Abstract: The application of atmospheric pressure plasmas in medicine is increasingly gaining attention in recent years, although very little is currently known about the plasma-induced processes occurring on the surface of living organisms. It is known that most bio-organisms, including bacteria, are coated by a liquid film surrounding them, and there might be many interactions between plasma species and the liquid layer before the plasma species reach the surface of the bio-organisms. Therefore, it is essential to study the behaviour of the reactive species in a liquid film, in order to determine whether these species can travel through this layer and reach the biomolecules, or whether new species are formed along the way. In this work, we investigate the interaction of reactive oxygen species (i.e. O, OH, HO2 and H2O2) with water, which is assumed as a simple model system for the liquid layer surrounding biomolecules. Our computational investigations show that OH, HO2 and H2O2 can travel deep into the liquid layer and are hence in principle able to reach the bio-organism. Furthermore, O, OH and HO2 radicals react with water molecules through hydrogen-abstraction reactions, whereas no H-abstraction reaction takes place in the case of H2O2. This study is important to gain insight into the fundamental operating mechanisms in plasma medicine, in general, and the interaction mechanisms of plasma species with a liquid film, in particular.

Abstract: Both geometrical and electrical asymmetry effects in capacitive argon discharges are investigated using a two-dimensional particle-in-cell coupled with Monte Carlo collision model. When changing the ratio of the top and bottom electrode surface areas and the phase shift between the two applied harmonics, the induced self-bias was found to develop separately. By adjusting the ratio between the high and low harmonic amplitudes, the electrical asymmetry effect at a fixed phase shift can be substantially optimized. However, the self-bias caused by the geometrical asymmetry hardly changed. Moreover, the separate control of these two asymmetry effects can also be demonstrated from their power absorption profiles. Both the axial and radial plasma density distributions can be modulated by the electrical asymmetry effect.

Abstract: A two-dimensional self-consistent fluid model combined with the full set of Maxwell equations is employed to investigate an Ar/CF4 capacitively coupled plasma, focusing on the phase-shift effect on the plasma characteristics at various frequencies and gas mixture ratios. When the discharge is sustained by a single frequency at 13.56 MHz in an Ar/CF4 mixture with a ratio of 0.9/0.1, no obvious difference is detected between the electron densities obtained in the so-called electrostatic model (with only the static electric fields taken into account) and the electromagnetic model (which includes the electromagnetic effects). However, as the frequency increases to 60 and 100 MHz, the difference becomes distinct, due to the significant influence of the electromagnetic effects. The phase-shift effect on the plasma radial uniformity has also been investigated in a dual frequency discharge, i.e. when the top driven source is switched on with a phase difference phiv ranging from 0 to π, in the frequency range 13.56100 MHz. At low concentration of CF4 (10%), Ar+ ions are the major positive ions in the entire range of frequencies. When the frequency is low, i.e. 13.56 MHz, the Ar+ density exhibits an off-axis peak at phiv = 0 due to the edge effect, and a better uniformity caused by the phase-shift modulation is obtained at phiv = π. At 60 MHz, the Ar+ density varies from edge-peaked at phiv = 0 to uniform (i.e. at phiv = 0.53π), and finally at phiv = π, a broad maximum is observed at the centre due to the standing-wave effect. As the frequency increases to 100 MHz, the best radial uniformity is reached at 0.25π, and the maximum moves again towards the radial wall in the reverse-phase case (phiv = π) due to the dominant skin effect. When the frequency is fixed at 100 MHz, the phase-shift control shows a different behaviour at a high concentration of CF4. For instance, the ${\rm CF}_3

Abstract: A two-dimensional fluid model, including the full set of Maxwell equations, has been developed and applied to investigate the effect of a phase shift between two power sources on the radial uniformity of several plasma characteristics in a hydrogen capacitively coupled plasma. This study was carried out at various frequencies in the range 13.56200 MHz. When the frequency is low, at 13.56 MHz, the plasma density is characterized by an off-axis peak when both power sources are in-phase (phgr = 0), and the best radial uniformity is obtained at phgr = π. This trend can be explained because the radial nonuniformity caused by the electrostatic edge effect can be effectively suppressed by the phase-shift effect at a phase difference equal to π. When the frequency rises to 60 MHz, the plasma density profiles shift smoothly from edge-peaked over uniform to centre-peaked as the phase difference increases, due to the pronounced standing-wave effect, and the best radial uniformity is reached at phgr = 0.3π. At a frequency of 100 MHz, a similar behaviour is observed, except that the maximum of the plasma density moves again towards the radial edge at the reverse-phase case (phgr = π), because of the dominant skin effect. When the frequency is 200 MHz, the bulk plasma density increases significantly with increasing phase-shift values, and a better uniformity is obtained at phgr = 0.4π. This is because the density in the centre increases faster than at the radial edge as the phase difference rises, due to the increasing power deposition Pz in the centre and the decreasing power density Pr at the radial edge. As the phase difference increases to π, the maximum near the radial edge becomes obvious again. This is because the skin effect has a predominant influence on the plasma density under this condition, resulting in a high density at the radial edge. Moreover, the axial ion flux increases monotonically with phase difference, and exhibits similar profiles to the plasma density. The calculation results illustrate that the radial uniformity of the various plasma characteristics is strongly dependent on the applied frequency and the phase shift between both power sources, which is important to realize, for controlling the uniformity of the plasma etch and deposition processes.