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: Although plasma catalysis is gaining increasing interest for various environmental applications, the underlying mechanisms are still far from understood. For instance, it is not yet clear whether and how plasma streamers can propagate in catalyst pores, and what is the minimum pore size to make this happen. As this is crucial information to ensure good plasma-catalyst interaction, we study here the mechanism of plasma streamer propagation in a catalyst pore, by means of a twodimensional particle-in-cell/Monte Carlo collision model, for various pore diameters in the nm range to μm-range. The so-called Debye length is an important criterion for plasma penetration into catalyst pores, i.e. a plasma streamer can penetrate into pores when their diameter is larger than the Debye length. The Debye length is typically in the order of a few 100 nm up to 1 μm at the conditions under study, depending on electron density and temperature in the plasma streamer. For pores in the range of ∼50 nm, plasma can thus only penetrate to some extent and at
very short times, i.e. at the beginning of a micro-discharge, before the actual plasma streamer reaches the catalyst surface and a sheath is formed in front of the surface. We can make plasma streamers penetrate into smaller pores (down to ca. 500 nm at the conditions under study) by increasing the applied voltage, which yields a higher plasma density, and thus reduces the Debye length. Our simulations also reveal that the plasma streamers induce surface charging of the catalyst pore sidewalls, causing discharge enhancement inside the pore, depending on pore diameter and depth.

Abstract: Purified water supply for human use, agriculture and industry is the major global priority nowadays. The advanced oxidation process based on atmospheric pressure non-thermal plasma (NTP) has been used for purification of wastewater, although the underlying mechanisms of degradation of organic pollutants are still unknown. In this study we employ two needle-type atmospheric pressure non-thermal plasma jets, i.e., indirect (ID-APPJ) and direct (D-APPJ) jets operating at Ar feed gas, for the treatment of methylene blue, methyl orange and congo red dyes, for two different times (i.e., 20 min and 30 min). Specifically, we study the decolorization/degradation of all three dyes using the above mentioned plasma sources, by means of UV-Vis spectroscopy, HPLC and a density meter. We also employ mass spectroscopy to verify whether only decolorization or also degradation takes place after treatment of the dyes by the NTP jets. Additionally, we analyze the interaction of OH radicals with all three dyes using reactive molecular dynamics simulations, based on the density functional-tight binding method. This investigation represents the first report on the degradation of these three different dyes by two types of NTP setups, analyzed by various methods, and based on both experimental and computational studies.

Abstract: A fluid model is self-consistently established to investigate the harmonic effects in an inductively coupled plasma, where the electromagnetic field is solved by the finite difference time domain technique. The spatiotemporal distribution of harmonic current density, harmonic potential, and other plasma quantities, such as radio frequency power deposition, plasma density, and electron temperature, have been investigated. Distinct differences in current density have been observed when calculated with and without Lorentz force, which indicates that the nonlinear Lorentz force plays an important role in the harmonic effects, especially at low frequencies. Moreover, the even harmonics are larger than the odd harmonics both in the current density and the potential. Finally, the dependence of various plasma quantities with and without the Lorentz force on various driving frequencies is also examined. It is shown that the deposited power density decreases and the depth of penetration increases slightly because of the Lorentz force. The electron density increases distinctly while the electron temperature remains almost the same when the Lorentz force is taken into account.

Abstract: In this paper, the negative ion behavior in a C4F8 inductively coupled plasma (ICP) is investigated using a hybrid model. The model predicts a non-monotonic variation of the total negative ion density with power at low pressure (1030 mTorr), and this trend agrees well with experiments that were carried out in many fluorocarbon (fc) ICP sources, like C2F6, CHF3, and C4F8. This behavior is explained by the availability of feedstock C4F8 gas as a source of the negative ions, as well as by the presence of low energy electrons due to vibrational excitation at low power. The maximum of the negative ion density shifts to low power values upon decreasing pressure, because of the more pronounced depletion of C4F8 molecules, and at high pressure (∼50 mTorr), the anion density continuously increases with power, which is similar to fc CCP sources. Furthermore, the negative ion composition is identified in this paper. Our work demonstrates that for a clear understanding of the negative ion behavior in radio frequency C4F8 plasma sources, one needs to take into account many factors, like the attachment characteristics, the anion composition, the spatial profiles, and the reactor configuration. Finally, a detailed comparison of our simulation results with experiments is conducted.

Abstract: In this paper we study the cold atmospheric pressure plasma jet, called kinpen, operating in Ar with different admixture fractions up to 1% pure , and + . Moreover, the device is operating with a gas curtain of dry air. The absolute net production rates of the biologically active ozone () and nitrogen dioxide () species are measured in the far effluent by quantum cascade laser absorption spectroscopy in the mid-infrared. Additionally, a zero-dimensional semi-empirical reaction kinetics model is used to calculate the net production rates of these reactive molecules, which are compared to the experimental data. The latter model is applied throughout the entire plasma jet, starting already within the device itself. Very good qualitative and even quantitative agreement between the calculated and measured data is demonstrated. The numerical model thus yields very useful information about the chemical pathways of both the and the generation. It is shown that the production of these species can be manipulated by up to one order of magnitude by varying the amount of admixture or the admixture type, since this affects the electron kinetics significantly at these low concentration levels.

Abstract: Cyrogenic etching of silicon is envisaged to enable better control over plasma processing in the microelectronics industry, albeit little is known about the fundamental differences compared to the room temperature process. We here present molecular dynamics simulations carried out to obtain sticking probabilities, thermal desorption rates, surface diffusion speeds, and sputter yields of F, F2, Si, SiF, SiF2, SiF3, SiF4, and the corresponding ions on Si(100) and on SiF13 surfaces, both at cryogenic and near room temperature. The different surface behavior during conventional etching and cryoetching is discussed. F2 is found to be relatively reactive compared to other species like SiF03. Thermal desorption occurs at a significantly lower rate under cryogenic conditions, which results in an accumulation of physisorbed species. Moreover, ion incorporation is often observed for ions with energies of 30400 eV, which results in a relatively low net sputter yield. The obtained results suggest that the actual etching of Si, under both cryogenic and near room temperature conditions, is based on the complete conversion of the Si surface to physisorbed SiF4, followed by subsequent sputtering of these molecules, instead of direct sputtering of the SiF03 surface.

Abstract: The mechanism of interaction of cold nonequilibrium plasma jets with mammalian cells in physiologic liquid is reported. The major biological active species produced by an argon RF plasma jet responsible for cell viability reduction are analyzed by experimental results obtained through physical, biological, and chemical diagnostics. This is complemented with chemical kinetics modeling of the plasma source to assess the dominant reactive gas phase species. Different plasma chemistries are obtained by changing the feed gas composition of the cold argon based RF plasma jet from argon, humidified argon (0.27%), to argon/oxygen (1%) and argon/air (1%) at constant power. A minimal consensus physiologic liquid was used, providing isotonic and isohydric conditions and nutrients but is devoid of scavengers or serum constituents. While argon and humidified argon plasma led to the creation of hydrogen peroxide dominated action on the mammalian cells, argonoxygen and argonair plasma created a very different biological action and was characterized by trace amounts of hydrogen peroxide only. In particular, for the argonoxygen (1%), the authors observed a strong negative effect on mammalian cell proliferation and metabolism. This effect was distance dependent and showed a half life time of 30 min in a scavenger free physiologic buffer. Neither catalase and mannitol nor superoxide dismutase could rescue the cell proliferation rate. The strong distance dependency of the effect as well as the low water solubility rules out a major role for ozone and singlet oxygen but suggests a dominant role of atomic oxygen. Experimental results suggest that O reacts with chloride, yielding Cl2 − or ClO−. These chlorine species have a limited lifetime under physiologic conditions and therefore show a strong time dependent biological activity. The outcomes are compared with an argon MHz plasma jet (kinpen) to assess the differences between these (at least seemingly) similar plasma sources.

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: We developed a microkinetic model to reveal the effects of plasma-generated radicals, intermediates, and vibrationally excited species on the catalytic hydrogenation of CO2 to CH3OH on a Cu(111) surface. As a benchmark, we first present the mechanisms of thermal catalytic CH3OH formation. Our model predicts that the reverse water-gas shift reaction followed by CO hydrogenation, together with the formate path, mainly contribute to CH3OH formation in thermal catalysis. Adding plasma-generated radicals and intermediates results in a higher CH3OH turnover frequency (TOF) by six to seven orders of magnitude, showing the potential of plasma-catalytic CO2 hydrogenation into CH3OH, in accordance with the literature. In addition, CO2 vibrational excitation further increases the CH3OH TOF, but the effect is limited due to relatively low vibrational temperatures under typical plasma catalysis conditions. The predicted increase in CH3OH formation by plasma catalysis is mainly attributed to the increased importance of the formate path. In addition, the conversion of plasma-generated CO to HCO* and subsequent HCOO* or H2CO* formation contribute to CH3OH formation. Both pathways bypass the HCOO* formation from CO2, which is the main bottleneck in the process. Hence, our model points toward the important role of CO, but also O, OH, and H radicals, as they influence the reactions that consume CO2 and CO. In addition, our model reveals that the H pressure should not be smaller than ca. half of the O pressure in the plasma as this would cause O* poisoning, which would result in very small product TOFs. Thus, plasma conditions should be targeted with a high CO and H content as this is favorable for CH3OH formation, while the O content should be minimized.

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: 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: The electron bounce resonance heating (BRH) in dual-frequency capacitively coupled plasmas operated in oxygen is studied by different experimental methods and a particle-in-cell/Monte Carlo collision (PIC/MCC) simulation, and compared with the electropositive argon discharge. In comparison with argon, the experimental results show that in an oxygen discharge the resonance peaks in positive-ion density and light intensity tend to occur at larger electrode gaps. Moreover, at electrode gaps L > 2.5 cm, the positive-ion (and electron) density and the light emission drop monotonically in the oxygen discharge upon increasing L, whereas they rise (after an initial drop) in the argon case. At resonance gap the electronegativity reaches its maximum due to the BRH. All these experimental observations are explained by PIC/MCC simulations, which show that in the oxygen discharge the bulk electric field becomes quite strong and is out of phase with the sheath field. Therefore, it retards the resonance electrons when traversing the bulk, resulting in a suppressed BRH. Both experiment and simulation results show that this effect becomes more pronounced at lower high-frequency power, when the discharge mode changes from electropositive to electronegative. In a pure oxygen discharge, the BRH is suppressed with increasing pressure and almost diminishes at 12 Pa. Finally, the driving frequency significantly affects the BRH, because it determines the phase relation between bulk electric field and sheath electric field.