Abstract: We present a 3D and 2D Cartesian quasi-neutral plasma model for a low current argon gliding arc discharge, including strong interactions between the gas flow and arc plasma column.
The 3D model is applied only for a short time of 0.2 ms due to its huge computational cost. It mainly serves to verify the reliability of the 2D model. As the results in 2D compare well with those in 3D, they can be used for a better understanding of the gliding arc basic characteristics. More specifically, we investigate the back-breakdown phenomenon induced by an artificially controlled plasma channel, and we discuss its effect on the gliding arc characteristics. The
back-breakdown phenomenon, or backward-jump motion of the arc, as observed in the experiments, results in a drop of the gas temperature, as well as in a delay of the arc velocity with respect to the gas flow velocity, allowing more gas to pass through the arc, and thus increasing the efficiency of the gliding arc for gas treatment applications.

Abstract: We perform molecular dynamics simulations to study the flip-flop motion of phosphatidylserine (PS) across the plasma membrane upon increasing oxidation degree of the membrane. Our computational results show that an increase of the oxidation degree in the lipids leads to a decrease of the free energy barrier for translocation of PS through the membrane. In other words, oxidation of the lipids facilitates PS flip-flop motion across the membrane, because in native phospholipid bilayers this is only a “rare event” due to the high energy barriers for the translocation of PS. The present study provides an atomic-scale insight into the mechanisms of the PS flip-flop upon oxidation of lipids, as produced for example by cold atmospheric plasma, in living cells.

Abstract: Dielectric Barrier Discharges operating in CO and CO2 form solid products at atmospheric pressure. The main differences between both plasmas and their deposits were analyzed, at similar energy input. GC measurements revealed a mixture of CO2, CO, and O2 in the CO2 DBD exhaust, while no O2 was found in the CO plasma. A coating of nanoparticles composed of Fe, O, and C was produced by the CO2 discharge, whereas, a microscopic dendrite-like carbon structure was formed in the CO plasma. Fe3O4 and Fe crystalline phases were found in the CO2 sample. The CO
deposition was characterized as an amorphous structure, close to polymeric CO (p-CO). Interestingly, p-CO is not formed in the CO2 plasma, in spite of the significant amounts of CO produced (up to 30% in the reactor exhaust).

Abstract: The chemical complexity of non-equilibrium plasmas poses a challenge for plasma modeling because of the computational load. This paper presents a dimension reduction method for such chemically complex plasmas based on principal component analysis (PCA). PCA is used to identify a low-dimensional manifold in chemical state space that is described by a small number of parameters: the principal components. Reduction is obtained since continuity equations only need to be solved for these principal components and not for all the species. Application of the presented method to a CO2 plasma model including state-to-state vibrational kinetics of CO2 and CO demonstrates the potential of the PCA method for dimension reduction. A manifold described by only two principal components is able to predict the CO2 to CO conversion at varying ionization degrees very accurately.

Abstract: In this work, a hybrid model is used to investigate the effect of different gas ratios on the Si etching and polymer film deposition characteristics in an Ar/CF4 inductively coupled plasma. The influence of the surface processes on the bulk plasma properties is studied, and also the spatial characteristics of important gas phase and etched species. The densities of F and CF2 decrease when the surface module is included in the simulations, due to the species consumption caused by etching and polymer deposition. The influence of the surface processes on the bulk plasma depends on the Ar/CF4 gas ratio. The deposited polymer becomes thicker at high CF4 content because of more abundant CFx radicals. As a result of the competition between the polymer thickness and the F flux, the etch rate first increases and then decreases upon increasing the CF4 content. The electron properties, more specifically the electron density profile, affect the Si etch characteristics substantially by determining the radical density and flux profiles. In fact, the radial profile of the etch rate is more uniform at low CF4 content since the electron density has a smooth distribution. At high CF4 content, the etch rate is less uniform with a minimum halfway along the wafer radius, because the electron density distribution is more localized. Therefore, our calculations predict that it is better to work at relatively high Ar/CF4 gas ratios, in order to obtain high etch rate and good profile uniformity for etch applications. This, in fact, corresponds to the typical experimental etch conditions in Ar/CF4 gas mixtures as found in the literature, where Ar is typically present at a much higher concentration than CF4.

Abstract: We report on a molecular dynamics study of small classical two-dimensional clusters with ringlike configurations. We focus on the particles motion at low temperatures before the radial and angular melting sets in. It is shown that in magic number configurations a local radial melting of subshells occur, which is related to the intershell rotation.

Abstract: In this study we quantitatively investigate for the first time the plasma characteristics of an argon gliding arc with a 3D model. The model is validated by comparison with available experimental data from literature and a reasonable agreement is obtained for the calculated gas temperature and electron density. A complete arc cycle is modeled from initial ignition to arc decay. We investigate how the plasma characteristics, i.e., the electron temperature, gas temperature,
reduced electric field, and the densities of electrons, Ar+ and Ar2+ ions and Ar(4s) excited states, vary over one complete arc cycle, including their behavior in the discharge and post-discharge. These plasma characteristics exhibit a different evolution over one arc cycle, indicating that either the active discharge stage or the post-discharge stage can be beneficial for certain applications.

Abstract: In this work we have analyzed the properties of a gliding dc discharge in argon at atmospheric pressure. Despite the usual designation of these discharges as ‘gliding arc discharges’, it was found previously that they operate in two different regimes—glow and arc. Here we analyze the differences in both regimes by means of two dimensional fluid modeling. In order to address different aspects of the discharge operation, we use two models—Cartesian and axisymmetric in a cylindrical coordinate system. The obtained results show that the two types of discharges produce a similar plasma column for a similar discharge current. However, the different mechanisms of plasma channel attachment to the cathode could produce certain differences in the plasma parameters (i.e. arc elongation), and this can affect gas treatments applications.

Abstract: A comparative study of two plasma sources (floating-electrode dielectric barrier discharge, DBD, Drexel University; atmospheric pressure argon plasma jet, kINPen, INP Greifswald) on cancer cell toxicity was performed. Cell culture protocols, cytotoxicity assays, and procedures for assessment of hydrogen peroxide (H2O2) were standardized between both labs. The inhibitory concentration 50 (IC50) and its corresponding H2O2 deposition was determined for both devices. For the DBD, IC50 and H2O2 generation were largely dependent on the total energy input but not pulsing frequency, treatment time, or total number of cells. DBD cytotoxicity could not be replicated by addition of H2O2 alone and was inhibited by larger amounts of liquid present during the treatment. Jet plasma toxicity depended on peroxide generation as well as total cell number and amount of liquid. Thus, the amount of liquid present during plasma treatment in vitro is key in attenuating short-lived species or other physical effects from plasmas. These in vitro results suggest a role of liquids in or on tissues during plasma treatment in a clinical setting. Additionally, we provide a platform for correlation between different plasma sources for a predefined cellular response.

Abstract: In the present work, reactive magnetron sputtering in DC mode was used to grow complex oxide thin films, starting from two separate pure metal targets. A series of coatings was produced with a stoichiometry of the film ranging from MgO, over MgxAlyOz to Al2O3. The surface energy, crystallinity, hardness, refractive index, and surface roughness were investigated. A relationship between all properties studied and the Mg content of the samples was found. A critical compositional region for the Mg-Al-O system where all properties exhibit a change was noticed.

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: The formation process of a microdischarge (MD) in both μm- and nm-sized catalyst pores is simulated by a two-dimensional particle-in-cell/Monte Carlo collision model. A parallel-plate dielectric barrier discharge configuration in filamentary mode is considered in ambient air. The discharge is powered by a high voltage pulse. Our calculations reveal that a streamer can penetrate into the surface features of a porous catalyst and MDs can be formed inside both μm- and nm-sized pores, yielding ionization inside the pore. For the μm-sized pores, the ionization mainly occurs inside the pore, while for the nm-sized pores the ionization is strongest near and inside the pore. Thus, enhanced discharges near and inside the mesoporous catalyst are observed. Indeed, the maximum values of the electric field, ionization rate and electron density occur near and inside the pore. The maximum electric field and electron density inside the pore first increase when the pore size rises from 4 nm to 10 nm, and then they decrease for the 100 nm pore, due to
a more pronounced surface discharge for the smaller pores. However, the ionization rate is highest for the 100 nm pore due to the largest effective ionization region.

Abstract: Low-temperature plasma chemical kinetic models are particularly important to the plasma community. These models typically require dozens of inputs, especially rate coefficients. The latter are not always precisely known and it is not surprising that the error on the rate coefficient data can propagate to the model output. In this paper, we present a model that uses N = 400 different combinations of rate coefficients based on the uncertainty attributed to each rate coefficient, giving a good estimation of the uncertainty on the model output due to the rate coefficients. We demonstrate that the uncertainty varies a lot with the conditions and the type of output. Relatively low uncertainties (about 15%) are found for electron density and temperature, while the uncertainty can reach more than an order of magnitude for the population of the vibrational levels in some cases and it can rise up to 100% for the CO2 conversion. The reactions that are mostly responsible for the largest uncertainties are identified. We show that the conditions of pressure, gas temperature and power density have a great effect on the uncertainty and on which reactions lead to this uncertainty. In all the cases tested here, while the absolute values may suffer from large uncertainties, the trends observed in previous modeling work are still valid. Finally, in accordance with the work of Turner, a number of ‘good practices’ is recommended.