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 review paper, we present several examples of reactive molecular dynamics simulations, which contribute to a better understanding of the underlying mechanisms in plasma medicine on the atomic scale. This includes the interaction of important reactive oxygen plasma species with the outer cell wall of both gram-positive and gram-negative bacteria, and with lipids present in human skin. Moreover, as most biomolecules are surrounded by a liquid biofilm, the behavior of these plasma species in a liquid (water) layer is presented as well. Finally, a perspective for future atomic scale modeling studies is given, in the field of plasma medicine in general, and for cancer treatment in particular.

Abstract: In this paper, we show that the apparent delocalization of the conduction band reported from first-principles simulations for the high-mobility amorphous oxide semiconductor InGaZnO4 (a-IGZO) is an artifact induced by the periodic conditions imposed to the model. Given a sufficiently large unit-cell dimension (over 40 angstrom), the conduction band becomes localized. Such a model size is up to four times the size of commonly used models for the study of a-IGZO. This finding challenges the analyses done so far on the nature of the defects and on the interpretation of numerous electrical measurements. In particular, we re-interpret the meaning of the computed effective mass reported so far in literature. Our finding also applies to materials such as SiZnSnO, ZnSnO, InZnSnO, In2O3 or InAlZnO4 whose models have been reported to display a fully delocalized conduction band in the amorphous phase.

Abstract: In this paper, we present the results of combined molecular dynamics−Metropolis Monte Carlo (MD-MMC) simulations of hydrocarbon species at flat diamond (100)2 × 1 and (111)1 × 1 surfaces. The investigated species are considered to be the most important growth species for (ultra)nanocrystalline diamond ((U)NCD) growth. When applying the MMC algorithm to stuck species at monoradical sites, bonding changes are only seen for CH2. The sequence of the bond breaking and formation as put forward by the MMC simulations mimics the insertion of CH2 into a surface dimer as proposed in the standard growth model of diamond. For hydrocarbon species attached to two adjacent radical (biradical) sites, the MMC simulations give rise to significant changes in the bonding structure. For UNCD, the combinations of C3 and C3H2, and C3 and C4H2 (at diamond (100)2 × 1) and C and C2H2 (at diamond (111)1 × 1) are the most successful in nucleating new crystal layers. For NCD, the following combinations pursue the diamond structure the best: C2H2 and C3H2 (at diamond (100)2 × 1) and CH2 and C2H2 (at diamond (111)1 × 1). The different behaviors of the hydrocarbon species at the two diamond surfaces are related to the different sterical hindrances at the diamond surfaces.

Abstract: In this paper, two types of melanoma and glioblastoma cancer cell lines are treated with cold atmospheric plasma to assess the effect of several parameters on the cell viability. The cell viability decreases with treatment duration and time until analysis in all cell lines with varying sensitivity. The majority of dead cells stains both AnnexinV (AnnV) and propidium iodide, indicating that the plasma-treated non-viable cells are mostly late apoptotic or necrotic. Genetic mutations might be involved in the response to plasma. Comparing the effects of two gas mixtures, as well as indirect plasma-activated medium versus direct treatment, gives different results per cell line. In conclusion, this study confirms the potential of plasma for cancer therapy and emphasizes the influence of experimental parameters on therapeutic outcome.

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

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, 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, 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, a self-consistent fluid model is developed focusing on the plasma parameters in capacitively coupled 20% N 2-80% Ar discharges. Measurements of ion density are performed with the help of a floating double probe, and the emission intensities from Ar(4p) and N 2 ( B ) transitions are detected by an optical emission spectroscopy to estimate their relative densities. The consistency between the numerical and experimental results confirms the reliability of the simulation. Then the plasma characteristics, specifically the reaction mechanisms of active species, are analyzed under various voltages. The increasing voltage leads to a monotonous increase in species density, whereas a less homogeneous radial distribution is observed at a higher voltage. Due to the high concentration of Ar gas, Ar + becomes the main ion, followed by the N 2 +</mml:msubsup> ion. Besides the electron impact ionization of neutrals, the charge transfer processes of Ar +/ N 2 and N 2 +</mml:msubsup>/Ar are found to have an impact on the ionic species. The results indicate that adopting the lower charge transfer reaction rate coefficients weakens the Ar + ion density and yields a higher N 2 +</mml:msubsup> ion density. However, the effect on the species spatial distributions and other species densities is limited. As for the excited-state species, the electron impact excitation of background gases remains overwhelming in the formation of Ar(4p), N 2 ( B ), and N 2 ( a ' ), whereas the <mml:msub> N 2 ( A ) molecules are mainly formed by the decay of <mml:msub> N 2 ( B ). In addition, the dissociation of <mml:msub> N 2 collided by excited-state Ar atoms dominates the N generation, which are mostly depleted to produce N + ions.

Abstract: In this paper we used a one dimensional fluid model, for the simulations of a Dielectric Barrier Discharge at atmospheric pressure. From the current and voltage profiles and the density profiles, we notice that two different regimes can be obtained in a uniform DBD. Furthermore a two dimensional flud model was developed and we describe how the gasflow can be included in such a model.

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: In this paper we employ first-principles calculations to investigate the effect of substitutional Si doping in the amorphous calcium-phosphate (a-HAP) structure on the work of adhesion, integral charge transfer, charge density difference and theoretical tensile strengths between an a-HAP coating and amorphous titanium dioxide (a-TiO2) substrate systemically. Our calculations demonstrate that substitution of a P atom by a Si atom in a-HAP (a-Si-HAP) with the creation of OH-vacancies as charge compensation results in a significant increase of the bonding strength of the coating to the substrate. The work of adhesion of the optimized Si-doped interfaces reaches a value of up to -2.52 J m(-2), which is significantly higher than for the stoichiometric a-HAP/a-TiO2. Charge density difference analysis indicates that the dominant interactions at the interface have significant covalent character, and in particular two Ti-O and three Ca-O bonds are formed for a-Si-HAP/a-TiO2 and one Ti-O and three Ca-O bonds for a-HAP/a-TiO2. From the stress-strain curve, the Young's modulus of a-Si-HAP/a-TiO2 is calculated to be about 25% higher than that of the a-HAP/a-TiO2, and the yielding stress is about 2 times greater than that of the undoped model. Our calculations therefore demonstrate that the presence of Si in the a-HAP structure strongly alters not only the bioactivity and resorption rates, but also the mechanical properties of the a-HAP/a-TiO2 interface. The results presented here provide an important theoretical insight into the nature of the chemical bonding at the a-HAP/a-TiO2 interface, and are particularly significant for the practical medical applications of HAP-based biomaterials.

Abstract: In this paper the various elementary plasma—surface interaction processes occurring in plasma catalysis are critically evaluated. Specifically, plasma catalysis at atmospheric pressure is considered. The importance of the various processes is analyzed for the most common plasma catalysis sources, viz. the dielectric barrier discharge and the gliding arc. The role and importance of surface chemical reactions (including adsorption, surface-mediated association and dissociation reactions, and desorption), plasma-induced surface modification, photocatalyst activation, heating, charging, surface discharge formation and electric field enhancement are discussed in the context of plasma catalysis. Numerous examples are provided to demonstrate the importance of the various processes.