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: Plasma catalysis holds great promise for environmental applications, provided that the process viability can be maximized in terms of energy efficiency and product selectivity. This requires a fundamental understanding of the various processes taking place and especially the mutual interactions between plasma and catalyst. In this review, we therefore first examine the various effects of the plasma on the catalyst and of the catalyst on the plasma that have been described in the literature. Most of these studies are purely experimental. The urgently needed fundamental understanding of the mechanisms underpinning plasma catalysis, however, may also be obtained through modelling and simulation. Therefore, we also provide here an overview of the modelling efforts that have been developed already, on both the atomistic and the macroscale, and we identify the data that can be obtained with these models to illustrate how modelling and simulation may contribute to this field. Last but not least, we also identify future modelling opportunities to obtain a more complete understanding of the various underlying plasma catalytic effects, which is needed to provide a comprehensive picture of plasma catalysis.

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 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: While plasma treatment of skin diseases and wound healing has been proven highly effective, the underlying mechanisms, and more generally the effect of plasma radicals on skin tissue, are not yet completely understood. In this paper, we perform ReaxFF-based reactive molecular dynamics simulations to investigate the interaction of plasma generated OH radicals with a model system composed of free fatty acids, ceramides, and cholesterol molecules. This model system is an approximation of the upper layer of the skin (stratum corneum). All interaction mechanisms observed in our simulations are initiated by H-abstraction from one of the ceramides. This reaction, in turn, often starts a cascade of other reactions, which eventually lead to the formation of aldehydes, the dissociation of ceramides or the elimination of formaldehyde, and thus eventually to the degradation of the skin barrier function.

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

Abstract: Both geometrical and electrical asymmetry effects in capacitive argon discharges are investigated using a two-dimensional particle-in-cell coupled with Monte Carlo collision model. When changing the ratio of the top and bottom electrode surface areas and the phase shift between the two applied harmonics, the induced self-bias was found to develop separately. By adjusting the ratio between the high and low harmonic amplitudes, the electrical asymmetry effect at a fixed phase shift can be substantially optimized. However, the self-bias caused by the geometrical asymmetry hardly changed. Moreover, the separate control of these two asymmetry effects can also be demonstrated from their power absorption profiles. Both the axial and radial plasma density distributions can be modulated by the electrical asymmetry effect.

Abstract: In this work, the role of vibrationally excited HBr (HBr(vib)) is computationally investigated for a HBr/He inductively coupled plasma applied for Si etching. It is found that at least 50% of all dissociations of HBr occur through HBr(vib). This additional dissociation pathway through HBr(vib) makes the plasma significantly more atomic. It also results in a slightly higher electron temperature (i.e. about 0.2 eV higher compared to simulation results where HBr(vib) is not included), as well as a higher gas temperature (i.e. about 50 K higher than without including HBr(vib)), due to the enhanced Franck–Condon heating through HBr(vib) dissociation,
at the conditions investigated. Most importantly, the calculated etch rate with HBr(vib) included in the model is a factor 3 higher than in the case without HBr(vib), due to the higher fluxes of etching species (i.e. H and Br), while the chemical composition of the wafer surface shows no significant difference. Our calculations clearly show the importance of including HBr(vib) for accurate modeling of HBr-containing plasmas.

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: Live-cell imaging with fluorescence microscopy is a powerful tool, especially in cancer research, widely-used for capturing dynamic cellular processes over time. However, light-induced toxicity (phototoxicity) can be incurred from this method, via disruption of intracellular redox balance and an overload of reactive oxygen species (ROS). This can introduce confounding effects in an experiment, especially in the context of evaluating and screening novel therapies. Here, we aimed to unravel whether phototoxicity can impact cellular homeostasis and response to non-thermal plasma (NTP), a therapeutic strategy which specifically targets the intracellular redox balance. We demonstrate that cells incorporated with a fluorescent reporter for live-cell imaging have increased sensitivity to NTP, when exposed to ambient light or fluorescence excitation, likely through altered proliferation rates and baseline intracellular ROS levels. These changes became even more pronounced the longer the cells stayed in culture. Therefore, our results have important implications for research implementing this analysis technique and are particularly important for designing experiments and evaluating redox-based therapies like NTP.

Abstract: A two-dimensional hybrid Monte Carlofluid model is applied to study HBr/He inductively coupled plasmas used for etching of Si. Complete sets of gas-phase and surface reactions are presented and the effects of the gas mixing ratio on the plasma characteristics and on the etch rates are discussed. A comparison with experimentally measured etch rates is made to validate the modelling results. The etch rate in the HBr plasma is found to be quite low under the investigated conditions compared to typical etch rates of Si with F- or Cl-containing gases. This allows for a higher control and fine-tuning of the etch rate when creating ultra-small features. Our calculations predict a higher electron temperature at higher He fraction, because the electrons do not lose their energy so efficiently in vibrational and rotational excitations. As a consequence, electron impact ionization and dissociation become more important, yielding higher densities of ions, electrons and H atoms. This results in more pronounced sputtering of the surface. Nevertheless, the overall etch rate decreases upon increasing He fraction, suggesting that chemical etching is still the determining factor for the overall etch rate.