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 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.

Abstract: In recent years, the ability of cold atmospheric pressure plasmas (CAPS) to selectively induce cell death in cancer cells has been widely established. This selectivity has been assigned to the reactive oxygen and nitrogen species (RONS) created in CAPs. To provide new insights in the search for an explanation
for the observed selectivity, we calculate the transfer free energy of multiple ROS across membranes containing a varying amount of cholesterol. The cholesterol fraction is investigated as a selectivity parameter because membranes of cancer cells are known to contain lower fractions of cholesterol compared to healthy cells. We find that cholesterol has a significant effect on the permeation of
reactive species across a membrane. Indeed, depending on the specific reactive species, an increasing cholesterol fraction can lead to (i) an increase of the transfer free energy barrier height and width, (ii) the formation of a local free energy minimum in the center of the membrane and (iii) the creation of extra free energy barriers due to the bulky sterol rings. In the context of plasma oncology, these observations suggest that the increased ingress of RONS in cancer cells can be explained by the decreased cholesterol fraction of their cell membrane.

Abstract: In recent years, the potential use of hydrogen as a clean energy source has gained considerable attention. Especially H2 formation by Ni-catalyzed reforming of methane at elevated temperatures is an attractive process. However, a more fundamental knowledge at the atomic level is needed for a full comprehension of the reactions at the catalyst surface. In this contribution, we therefore investigate the H2 formation after CHx impacts on a Ni(1 1 1) surface in the temperature range 4001600 K, by means of reactive molecular dynamics (MD) simulations using the ReaxFF potential. While some H2 formation is already observed at the lower temperatures, substantial H2 formation is only obtained at elevated temperatures of 1400 K and above. At 1600 K, the H2 molecules are even the most frequently formed species. In direct correlation with the increasing dehydrogenation at elevated temperatures, an increased surface-to-subsurface C-diffusivity is observed as well. This study highlights the major importance of the temperature on the H2 formation.

Abstract: In the last decades, non-thermal plasma has been extensively investigated as a relevant tool for various biomedical applications, ranging from tissue decontamination to regeneration and from skin treatment to tumor therapies. This high versatility is due to the different kinds and amount of reactive oxygen and nitrogen species that can be generated during a plasma treatment and put in contact with the biological target. Some recent studies report that solutions of biopolymers with the ability to generate hydrogels, when treated with plasma, can enhance the generation of reactive species and influence their stability, resulting thus in the ideal media for indirect treatments of biological targets. The direct effects of the plasma treatment on the structure of biopolymers in water solution, as well as the chemical mechanisms responsible for the enhanced generation of RONS, are not yet fully understood. In this study, we aim at filling this gap by investigating, on the one hand, the nature and extent of the modifications induced by plasma treatment in alginate solutions, and, on the other hand, at using this information to explain the mechanisms responsible for the enhanced generation of reactive species as a consequence of the treatment. The approach we use is twofold: (i) investigating the effects of plasma treatment on alginate solutions, by size exclusion chromatography, rheology and scanning electron microscopy and (ii) study of a molecular model (glucuronate) sharing its chemical structure, by chromatography coupled with mass spectrometry and by molecular dynamics simulations. Our results point out the active role of the biopolymer chemistry during direct plasma treatment. Short-lived reactive species, such as OH radicals and O atoms, can modify the polymer structure, affecting its functional groups and causing partial fragmentation. Some of these chemical modifications, like the generation of organic peroxide, are likely responsible for the secondary generation of long-lived reactive species such as hydrogen peroxide and nitrite ions. This is relevant in view of using biocompatible hydrogels as vehicles for storage and delivery reactive species for targeted therapies.

Abstract: In this computational study, a gliding arc plasma reactor with a reverse-vortex flow stabilization is modelled for the first time by a fluid plasma description. The plasma reactor operates with argon gas at atmospheric pressure. The gas flow is simulated using the k-ε Reynolds-averaged Navier–Stokes turbulent model. A quasi-neutral fluid plasma model is used for computing the plasma properties. The plasma arc movement in the reactor is observed, and the results for the gas flow, electrical characteristics, plasma density, electron temperature, and gas temperature are analyzed.

Abstract: In this experimental study, a flowing dielectric barrier discharge operating at atmospheric pressure is used for the splitting of CO2 into O2 and CO. The influence of the applied frequency and plasma power on the microdischarge properties is investigated to understand their role on the CO2 conversion. Electrical measurements are carried out to explain the conversion trends and to characterize the microdischarges through their number, their lifetime,
their intensity and the induced electrical charge. Their influence on the gas and electrode temperatures is also evidenced through optical emission spectroscopy and infrared imaging. It is shown that, in our configuration, the conversion depends mostly on the charge delivered in the plasma and not on the effective plasma voltage when the applied power is modified. Similarly, at constant total current, a better conversion is observed at low frequencies, where a less filamentary discharge regime with a higher effective plasma voltage than that at a higher
frequency is obtained.

Abstract: In this letter, we investigate the effect of various levels of nitrogen impurity on the electrical performance of an atmospheric pressure dielectric barrier discharge in helium. We illustrate the different current profiles that are obtained, which exhibit one or more discharge pulses per half cycle and evaluate their performance in ionizing the discharge and dissipating the power. It is shown that flat and broad current profiles perform the best in ionizing the discharge and use the least amount of power per generated charged particle.

Abstract: In this paper we study dry reforming of methane (DRM) in a gliding arc plasmatron (GAP) in the presence of N₂and O₂. N₂is added to create a stable plasma at equal fractions of CO₂and CH₄, and because emissions from industrial plants typically contain N₂, while O₂is added to enhance the process. We test different gas mixing ratios to evaluate the conversion and energy cost. We obtain conversions between 31 and 52% for CO₂and between 55 and 99% for CH₄, with total energy costs between 3.4 and 5.0 eV per molecule, depending on the gas mixture. This is very competitive when benchmarked with the literature. In addition, we present a chemical kinetics model to obtain deeper insight in the underlying plasma chemistry. This allows determination of the major reaction pathways to convert CO₂and CH₄, in the presence of O₂and N₂, into CO and H₂. We show that N₂assists in the CO₂conversion, but part of the applied energy is also wasted in N₂excitation. Adding O₂enhances the CH₄conversion, and lowers the energy cost, while the CO₂conversion remains constant, and only slightly drops at the highest O₂fractions studied, when CH₄is fully oxidized into CO₂.

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: In this paper, a new implementation of the Metropolis Monte Carlo (MMC) algorithm is presented. When combining the MMC model with a molecular dynamics (MD) code, crystal growth by plasma-enhanced chemical vapor deposition can be simulated. As the MD part simulates impacts of growth species onto the surface on a time scale of picoseconds, the MMC algorithm simulates the slower adatom surface processes. The implementation includes a criterion for the selection of atoms that are allowed to be displaced during the simulation, and a criterion of after how many MMC cycles the simulation is stopped. We performed combined MD-MMC simulations for hydrocarbon species that are important for the growth of ultrananocrystalline diamond (UNCD) films at partially hydrogenated diamond surfaces, since this implementation is part of a study of the growth mechanisms of (ultra)nanocrystalline diamond films. Exemplary for adatom arrangements during the growth of UNCD, the adatom surface behavior of C and C2H2 at diamond (111)1 × 1, C and C4H2 at diamond (111)1 × 1 and C3 at diamond (100)2 × 1 has been investigated. For all cases, the diamond crystal structure is pursued under the influence of MMC simulation. Additional longer time-scale MD simulations put forward very similar structures, verifying the MMC algorithm. Nevertheless, the MMC simulation time is typically one order of magnitude shorter than the MD simulation time.

Abstract: In this paper, an O2 inductively coupled plasma used for plasma enhanced atomic layer deposition of Al2O3 thin films is investigated by means of modeling. This work intends to provide more information about basic plasma properties such as species densities and species fluxes to the substrate as a function of power and pressure, which might be hard to measure experimentally. For this purpose, a hybrid model developed by Kushner et al is applied to calculate the plasma characteristics in the reactor volume for different chamber pressures ranging from 1 to 10 mTorr and different coil powers ranging from 50 to 500 W. Density profiles of the various oxygen containing plasma species are reported as well as fluxes to the substrate under various operating conditions. Furthermore, different orientations of the substrate, which can be placed vertically or horizontally in the reactor, are taken into account. In addition, special attention is paid to the recombination process of atomic oxygen on the different reactor walls under the stated operating conditions. From this work it can be concluded that the plasma properties change significantly in different locations of the reactor. The plasma density near the cylindrical coil is high, while it is almost negligible in the neighborhood of the substrate. Ion and excited species fluxes to the substrate are found to be very low and negligible. Finally, the orientation of the substrate has a minor effect on the flux of O2, while it has a significant effect on the flux of O. In the horizontal configuration, the flux of atomic oxygen can be up to one order of magnitude lower than in the vertical configuration.

Abstract: In this paper, an overview is given of modeling activities going on in our research group, for describing the plasma chemistry and plasmasurface interactions in reactive plasmas. The plasma chemistry is calculated by a fluid approach or by hybrid Monte Carlo (MC)fluid modeling. An example of both is illustrated in the first part of the paper. The example of fluid modeling is given for a dielectric barrier discharge (DBD) in CH4/O2, to describe the partial oxidation of CH4 into value-added chemicals. The example of hybrid MCfluid modeling concerns an inductively coupled plasma (ICP) etch reactor in Ar/Cl2/O2, including also the description of the etch process. The second part of the paper deals with the treatment of plasmasurface interactions on the atomic level, with molecular dynamics (MD) simulations or a combination of MD and MC simulations.

Abstract: In this paper, simulations are performed to gain a better insight into the properties of a Cl2/Ar plasma, with and without O2, during plasma etching of Si. Both plasma and surface properties are calculated in a self-consistent manner. Special attention is paid to the behavior of etch products coming from the wafer or the walls, and how the chamber walls can affect the plasma and the resulting etch process. Two modeling cases are considered. In the first case, the reactor walls are defined as clean (Al2O3), whereas in the second case a SiO2 coating is introduced on the reactor walls before the etching process, so that oxygen will be sputtered from the walls and introduced into the plasma. For this reason, a detailed reaction set is presented for a Cl2/O2/Ar plasma containing etched species, as well as an extensive reaction set for surface processes, including physical and chemical sputtering, chemical etching and deposition processes. Density and flux profiles of various species are presented for a better understanding of the bulk plasma during the etching process. Detailed information is also given on the composition of the surfaces at various locations of the reactor, on the etch products in the plasma and on the surface loss probabilities of the plasma species at the walls, with different compositions. It is found that in the clean chamber, walls are mostly chlorinated (Al2Cl3), with a thin layer of etch products residing on the wall. In the coated chamber, an oxy-chloride layer is grown on the walls for a few nanometers during the etching process. The Cl atom wall loss probability is found to decrease significantly in the coated chamber, hence increasing the etch rate. SiCl2, SiCl4 and SiCl3 are found to be the main etch products in the plasma, with the fraction of SiCl2 being always slightly higher. The simulation results compare well with experimental data available from the literature.

Abstract: In this paper, the effects of prolonged bias-enhanced nucleation (prolonged BEN) on the growth mechanisms of diamond are investigated by molecular dynamics (MD) and combined MD-Metropolis Monte Carlo (MD-MMC) simulations. First, cumulative impacts of CxHy+ and Hx+ on an a-C:H/nanodiamond composite were simulated; second, nonconsecutive impacts of the dominant ions were simulated in order to understand the observed phenomena in more detail. As stated in the existing literature, the growth of diamond structures during prolonged BEN is a process that takes place below the surface of the growing film. The investigation of the penetration behavior of CxHy+ and Hx+ species shows that the carbon-containing ions remain trapped within this amorphous phase where they dominate mechanisms like precipitation of sp3 carbon clusters. The H+ ions, however, penetrate into the crystalline phase at high bias voltages (>100 V), destroying the perfect diamond structure. The experimentally measured reduction of grain sizes at high bias voltage, reported in the literature, might thus be related to penetrating H+ ions. Furthermore, the CxHy+ ions are found to be the most efficient sputtering agents, preventing the build up of defective material.