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: High density polyethylene surfaces were exposed to the atmospheric post-discharge of a radiofrequency plasma torch supplied in helium and oxygen. Dynamic water contact angle measurements were performed to evaluate changes in surface hydrophilicity and angle resolved x-ray photoelectron spectroscopy was carried out to identify the functional groups responsible for wettability changes and to study their subsurface depth profiles, up to 9 nm in depth. The reactions leading to the formation of CO, C = O and OC = O groups were simulated by molecular dynamics. These simulations demonstrate that impinging oxygen atoms do not react immediately upon impact but rather remain at or close to the surface before eventually reacting. The simulations also explain the release of gaseous species in the ambient environment as well as the ejection of low molecular weight oxidized materials from the surface.

Abstract: Plasma medicine is a rapidly evolving multidisciplinary field at the intersection of chemistry, biochemistry, physics, biology, medicine and bioengineering. It holds great potential in medical, health care, dentistry, surgical, food treatment and other applications. This multidisciplinary nature and variety of possible applications come along with an inherent and intrinsic complexity. Advancing plasma medicine to the stage that it becomes an everyday tool in its respective fields requires a fundamental understanding of the basic processes, which is lacking so far. However, some major advances have already been made through detailed experiments over the last 15 years. Complementary, computer simulations may provide insight that is difficultif not impossibleto obtain through experiments. In this review, we aim to provide an overview of the various simulations that have been carried out in the context of plasma medicine so far, or that are relevant for plasma medicine. We focus our attention mostly on atomistic simulations dealing with plasmabiomolecule interactions. We also provide a perspective and tentative list of opportunities for future modelling studies that are likely to further advance the field.

Abstract: A one-dimensional, implicit particle-in-cell Monte Carlo collision model is used to simulate the plasma kinetic properties at a steady state in a parallel-plate direct current argon glow microdischarge under various operating conditions, such as driving voltage (301000 V) and gap size (101000 µm) at atmospheric pressure. First, a comparison between rf and dc modes is shown for the same pressure, driving voltage and gap spacing. Furthermore, the effect of gap size scaling (in the range of 101000 µm) on the breakdown voltage, peak electron density and peak electron current density at the breakdown voltage is examined. The breakdown voltage is lower than 150 V in all gaps considered. The microdischarge is found to have a neutral bulk plasma region and a cathode sheath region with size varying with the applied voltage and the discharge gap. In our calculations, the electron and ion densities are of the order of 10181023 m−3, which is in the glow discharge limit, as the ionization degree is lower than 1% . The electron energy distribution function shows a two-energy group distribution at a gap of 10 µm and a three-energy group distribution at larger gaps such as 200 µm and 1000 µm, emphasizing the importance of the gap spacing in dc microdischarges.

Abstract: The discharge properties of the plasma bulk flow near the surface of heated compound-materials strongly affects the kinetic layer parameters modeled and manifested in the Knudsen layer. This paper extends the widely used two-layer kinetic ablation model to the ablation controlled non-equilibrium discharge due to the fact that the local thermodynamic equilibrium (LTE) approximation is often violated as a result of the interaction between the plasma and solid walls. Modifications to the governing set of equations, to account for this effect, are derived and presented by assuming that the temperature of the electrons deviates from that of the heavy particles. The ablation characteristics of one typical material, polytetrafluoroethylene (PTFE) are calculated with this improved model. The internal degrees of freedom as well as the average particle mass and specific heat ratio of the polyatomic vapor, which strongly depends on the temperature, pressure and plasma non-equilibrium degree and plays a crucial role in the accurate determination of the ablation behavior by this model, are also taken into account. Our assessment showed the significance of including such modifications related to the non-equilibrium effect in the study of vaporization of heated compound materials in ablation controlled arcs. Additionally, a two-temperature magneto-hydrodynamic (MHD) model accounting for the thermal non-equilibrium occurring near the wall surface is developed and applied into an ablation-dominated discharge for an electro-thermal chemical launch device. Special attention is paid to the interaction between the non-equilibrium plasma and the solid propellant surface. Both the mass exchange process caused by the wall ablation and plasma species deposition as well as the associated momentum and energy exchange processes are taken into account. A detailed comparison of the results of the non-equilibrium model with those of an equilibrium model is presented. The non-equilibrium results show a non-equilibrium region near the plasma-wall interaction region and this indicates the need for the consideration of the influence of the possible departure from LTE in the plasma bulk on the determination of ablation rate.

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: A hybrid model, called the hybrid plasma equipment model, was used to study Ar/Cl(2) inductively coupled plasmas used for the etching of Si. The effects of substrate bias, source power and gas pressure on the plasma characteristics and on the fluxes and energies of plasma species bombarding the substrate were observed. A comparison with experimentally measured etch rates was made to investigate how the etch process is influenced and which plasma species mainly account for the etch process. First, the general plasma characteristics are investigated at the following operating conditions: 10% Ar 90% Cl(2) gas mixture, 5mTorr total gas pressure, 100 sccm gas flow rate, 250W source power, -200V dc bias at the substrate electrode and an operating frequency of 13.56MHz applied to the coil and to the substrate electrode. Subsequently, the pressure is varied from 5 to 80mTorr, the substrate bias from -100 to -300V and the source power from 250 to 1000W. Increasing the total gas pressure results in a decrease of the etch rate and a less anisotropic flux to the substrate due to more collisions of the ions in the sheath. Increasing the substrate bias has an effect on the energy of the ions bombarding the substrate and to a lesser extent on the magnitude of the ion flux. When source power is increased, it was found that, not the energy, but the magnitude of the ion flux is increased. The etch rate was more influenced by a variation of the substrate bias than by a variation of the source power, at these operating conditions. These results suggest that the etch process is mainly affected by the energy of the ions bombarding the substrate and the magnitude of the ion flux, and to a lesser extent by the magnitude of the radical flux.

Abstract: A 2d3v particle-in-cell/Monte Carlo collisions model was applied to study the influence of the magnetic field strength on the cathode sheath region of a direct current (dc) magnetron discharge. When applying a magnetic field of 520-730 G, the cathode sheath width decreases with magnetic field strength, whereas, if a stronger magnetic field is applied (i. e. from 730 to 2600 G), the sheath width increases. This is explained by studying the structure of the sheath in different magnetic field strengths in terms of the electron and ion densities. The consequences of sheath structure on the sputter deposition process are also investigated. It is found that the magnetic field strength can control the erosion profile and the sputter rate.

Abstract: Using a molecular dynamics model the crystallinity of MgxAlyOz thin films with a variation in the stoichiometry of the thin film is studied at operating conditions similar to the experimental operating conditions of a dual magnetron sputter deposition system. The films are deposited on a crystalline or amorphous substrate. The Mg metal content in the film ranged from 100% (i.e. MgO film) to 0% (i.e. Al2O3 film). The radial distribution function and density of the films are calculated. The results are compared with x-ray diffraction and transmission electron microscopy analyses of experimentally deposited thin films by the dual magnetron reactive sputtering process. Both simulation and experimental results show that the structure of the MgAlO film varies from crystalline to amorphous when the Mg concentration decreases. It seems that the crystalline MgAlO films have a MgO structure with Al atoms in between.

Abstract: The atmospheric pressure dielectric barrier discharge in helium is a pulsed discharge in nature. If during the electrical current pulse a glow discharge is reached, then this pulse will last only a few microseconds in operating periods of sinusoidal voltage with lengths of about 10 to 100 µs. In this paper we demonstrate that right before a glow discharge is reached, the discharge very closely resembles the commonly assumed Townsend discharge structure, but actually contains some significant differing features and hence should not be considered as a Townsend discharge. In order to clarify this, we present calculation results of high time and space resolution of the pulse formation. The results indicate that indeed a maximum of ionization is formed at the anode, but that the level of ionization remains high and that the electric field at that time is significantly disturbed. Our results also show where this intermediate structure comes from.

Abstract: In this paper, some modelling approaches to describe direct current (dc) magnetron discharges developed in our research groups will be presented, including an analytical model, Monte Carlo simulations for the electrons and for the sputtered atoms, a hybrid Monte Carlo-fluid model and particle-in-cell-Monte Carlo collision simulations. The strengths and limitations of the various modelling approaches will be explained, and some characteristic simulation results will be illustrated. Furthermore, some other simulation methods related to the magnetron device will be briefly explained, more specifically for calculating the magnetic field distribution inside the discharge, and for describing the (reactive) sputtering.

Abstract: A hybrid model, called the hybrid plasma equipment model (HPEM), was used to study an inductively coupled plasma in gas mixtures of H2 or NH3 with CH4 or C2H2 used for the synthesis of carbon nanotubes or carbon nanofibres (CNTs/CNFs). The plasma properties are discussed for different gas mixtures at low and moderate pressures, and the growth precursors for CNTs/CNFs are analysed. It is found that C2H2, C2H4 and C2H6 are the predominant molecules in CH4 containing plasmas besides the feedstock gas, and serve as carbon sources for CNT/CNF formation. On the other hand, long-chain hydrocarbons are observed in C2H2-containing plasmas. Furthermore, the background gases CH4 and C2H2 show a different decomposition rate with H2 or NH3 addition at moderate pressures.

Abstract: As electron energy-loss spectroscopy (EELS) and x-ray absorption spectroscopy (XAS) probe the same transitions from coreshell states to unoccupied states above the Fermi energy, it should always be possible to apply the two techniques to the same physical phenomena, such as magnetic dichroism, and obtain the same information. Indeed, the similarity in the expression of the electron and x-ray cross-sections had been already exploited to prove the equivalence of x-ray magnetic linear dichroism and anisotropy in EELS, by noting that the polarization vector of a photon plays the same role as the momentum transfer in electron scattering. Recently, the same was proven true for x-ray magnetic circular dichroism (XMCD) by establishing a new TEM technique called EMCD (electron energy-loss magnetic chiral dichroism) (Schattschneider P et al 2006 Nature 441 4868), which makes use of special electron scattering conditions to force the absorption of a circularly polarized virtual photon. The intrinsic advantage of EMCD over XMCD is the high spatial resolution of electron microscopes, which are readily available. Among the particular obstacles in EMCD that do not exist for synchrotron radiation, is the notoriously low signal and the very particular scattering conditions necessary to observe a chiral dichroic signal. In spite of that, impressive progress has been made in recent years. The signal strength could be considerably increased, and some innovations such as using a convergent beam have been introduced. EMCD has evolved into several techniques, which make full use of the versatility of the TEM and energy filtering, spectroscopy or STEM conditions (Rubino S 2007 Magnetic circular dichroism in the transmission electron microscope PhD Thesis Vienna University of Technology, Vienna, Austria).