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: An inductively coupled plasma, connected to the sampling cone of a mass spectrometer, is computationally investigated. The occurrence of rotational motion of the auxiliary and carrier gas flows is studied. The effects of operating parameters, i.e., applied power and gas flow rates, as well as geometrical parameters, i.e., sampler orifice diameter and injector inlet diameter, are investigated. Our calculations predict that at higher applied power the auxiliary and carrier gas flows inside the torch move more forward to the sampling cone, which is validated experimentally for the auxiliary gas flow, by means of an Elan 6000 ICP-MS. Furthermore, an increase of the gas flow rates can also modify the occurrence of rotational motion. This is especially true for the carrier gas flow rate, which has a more pronounced effect to reduce the backward motion than the flow rates of the auxiliary and cooling gas. Moreover, a larger sampler orifice (e.g., 2 mm instead of 1 mm) reduces the backward flow of the auxiliary gas path lines. Finally, according to our model, an injector inlet of 2 mm diameter causes more rotations in the carrier gas flow than an injector inlet diameter of 1.5 mm, which can be avoided again by changing the operating parameters.

Abstract: A theoretical study of ultrashort laser-induced electron emission from a negatively biased metallic cathode has been performed. Classical as well as tunneling electron emission mechanisms are considered. It was found that electron emission is governed by an interplay of processes inside as well as above the cathode. A hybrid model is proposed, where the electron distribution within the target is retrieved from Boltzmann scattering integrals, while the charge distribution above it is studied by a Particle-In-Cell simulation. The results indicate that non-equilibrium effects determine the initial emission process, whereas the space charge above the target suppresses the effectively emitted charge.

Abstract: A study of ns-laser ablation is presented, which focuses on the transient behavior of the physical processes that act in and above a copper sample. A dimensionless multiphase collisional radiative model describes the interplay between the ablation, collisional, and radiative mechanisms. Calculations are done for a 6 ns-Nd:YAG laser pulse operating at 532 nm and fluences up to 15 J/cm2. Temporal intensity profiles as well as transmissivities are in good agreement with experimental results. It is found that volumetric ablation mechanisms and photo-processes both play an essential role in the onset of ns-laser induced breakdown.

Abstract: The present study focuses on the role of mass removal mechanisms in ns-laser ablation. A copper sample is placed in argon, initially set at standard pressure and temperature. Calculations are performed for a 6 ns laser pulse with a wavelength of 532 nm and laser fluences up to 10 J/cm2. The transient behavior in and above the copper target is described by a hydrodynamic model. Transmission profiles and ablation depths are compared with experimental results and similar trends are found. Our calculations reveal an interesting self-inhibiting mechanism: volumetric mass removal in the supercritical region triggers plasma shielding and therefore stops proceeding. This self-limiting process indicates that volumetric mass removal does not necessarily result in large ablation depths.

Abstract: A zero-dimensional, semi-empirical model is used to describe the plasma chemistry in an argon plasma jet flowing into humid air, mimicking the experimental conditions of a setup from the Eindhoven University of Technology. The model provides species density profiles as a function of the position in the plasma jet device and effluent. A reaction chemistry set for an argon/humid air mixture is developed, which considers 84 different species and 1880 reactions. Additionally, we present a reduced chemistry set, useful for higher level computational models. Calculated species density profiles along the plasma jet are shown and the chemical pathways are explained in detail. It is demonstrated that chemically reactive H, N, O and OH radicals are formed in large quantities after the nozzle exit and H2, O2(1Δg), O3, H2O2, NO2, N2O, HNO2 and HNO3 are predominantly formed as 'long living' species. The simulations show that water clustering of positive ions is very important under these conditions. The influence of vibrational excitation on the calculated electron temperature is studied. Finally, the effect of varying gas temperature, flow speed, power density and air humidity on the chemistry is investigated.

Abstract: The destruction of ethylene in a dielectric barrier discharge plasma is investigated by the combination of kinetic modeling and experiments, as a case study for plasma-based gas purification. The influence of the specific energy deposition on the removal efficiency and the selectivity toward CO and CO2 is studied for different concentrations of ethylene. The model allows the identication of the destruction pathway in dry and humid air. The latter is found to be mainly initiated by metastable N2 molecules, but the further destruction steps are dominated by O atoms and OH radicals. Upon increasing air humidity, the removal efficiency drops by ±15% (from 85% to 70%), but the selectivity toward CO and CO2 stays more or less constant at 60% and 22%, respectively. Beside CO and CO2, we also identified acetylene, formaldehyde, and water as byproducts of the destruction process, with concentrations of 1606 ppm, 15033 ppm, and 185 ppm in humid air (with 20% RH), respectively. Finally, we investigated the byproducts generated by the humid air discharge itself, which are the greenhouse gases O3, N2O, and the toxic gas NO2.

Abstract: We investigate the theoretical possibility of single layer graphene formation on a nickel surface at different substrate temperatures under far-from-equilibrium high precursor flux conditions, employing state-of-the-art hybrid reactive molecular dynamics/uniform acceptance force bias Monte Carlo simulations. It is predicted that under these conditions, the formation of a single layer graphene-like film may proceed through a combined depositionsegregation mechanism on a nickel substrate, rather than by pure surface segregation as is typically observed for metals with high carbon solubility. At 900 K and above, nearly continuous graphene layers are obtained. These simulations suggest that single layer graphene deposition is theoretically possible on Ni under high flux conditions.

Abstract: In this work particle-in-cell/Monte Carlo collision simulations are performed to study the heating mechanism and plasma characteristics in direct current (dc) superposed radio-frequency (RF) capacitively coupled plasmas, operated both in single-frequency (SF) and dual-frequency (DF) regimes. An RF (60/2 MHz) source is applied on the bottom electrode to sustain the discharge, and a dc source is fixed on the top electrode. The heating mechanism appears to be very different in dc superposed SF and DF discharges. When only a single source of 60 MHz is applied, the plasma bulk region is reduced by the dc source, thus the ionization rate and hence the electron density decrease with rising dc voltage. However, when a DF source of 60 and 2 MHz is applied, the electron density can increase upon addition of a dc voltage, depending on the gap length and applied dc voltage. This is explained from the spatiotemporal ionization rates in the DF discharge. In fact, a completely different behavior is observed for the ionization rate in the two half-periods of the LF source. In the first LF half-period, the situation resembles the dc superposed SF discharge, and the reduced plasma bulk region due to the negative dc bias results in a very small effective discharge area and a low ionization rate. On the other hand, in the second half-period, the negative dc bias is to some extent counteracted by the LF voltage, and the sheath close to the dc electrode becomes particularly thin. Consequently, the amplitude of the high-frequency sheath oscillations at the top electrode is largely enhanced, while the LF sheath at the bottom electrode is in its expanding phase and can thus well confine the high-energy electrons. Therefore, the ionization rate increases considerably in this second LF half-period. Furthermore, in addition to the comparison between SF and DF discharges and the effect of gap length and dc voltage, the effect of secondary electrons is examined.

Abstract: In this paper, we investigate the etch process of so-called low-k organic material by means of a N2/H2 capacitively coupled plasma, as applied in the micro-electronics industry for the manufacturing of computer chips. In recent years, such an organic material has emerged as a possible alternative for replacing bulk SiO2 as a dielectric material in the back-end-of-line, because of the smaller parasitic capacity between adjacent conducting lines, and thus a faster propagation of the electrical signals throughout the chip. Numerical simulations with a hybrid plasma model, using an extensive plasma and surface chemistry set, as well as experiments are performed, focusing on the plasma properties as well as the actual etching process, to obtain a better insight into the underlying mechanisms. Furthermore, the effects of gas pressure, applied power and gas composition are investigated to try to optimize the etch process. In general, the plasma density reaches a maximum near the wafer edge due to the so-called 'edge effect'. As a result, the etch rate is not uniform but will also reach its maximum near the wafer edge. The pressure seems not to have a big effect. A higher power increases the etch rate, but the uniformity becomes (slightly) worse. The gas mixing ratio has no significant effect on the etch process, except when a pure H2 or N2 plasma is used, illustrating the synergistic effects of a N2/H2 plasma. In fact, our calculations reveal that the N2/H2 plasma entails an ion-enhanced etch process. The simulation results are in reasonable agreement with the experimental values. The microscopic etch profile shows the desired anisotropic shape under all conditions under study.

Abstract: Nonthermal atmospheric pressure plasmas are gaining increasing attention for biomedical applications. However, very little fundamental information on the interaction mechanisms between the plasma species and biological cells is currently available. We investigate the interaction of important plasma species, such as OH, H2O2, O, O3, as well as O2 and H2O, with bacterial peptidoglycan by means of reactive molecular dynamics simulations, aiming for a better understanding of plasma disinfection. Our results show that OH, O, O3, and H2O2 can break structurally important bonds of peptidoglycan (i.e., CO, CN, or CC bonds), which consequently leads to the destruction of the bacterial cell wall. The mechanisms behind these breakups are, however, dependent on the impinging plasma species, and this also determines the effectiveness of the cell wall destruction.

Abstract: Structural defects inevitably appear during the nucleation event that determines the structure and properties of single-walled carbon nanotubes. By combining ion bombardment experiments with atomistic simulations we reveal that ion bombardment in a suitable energy range allows these defects to be healed resulting in an enhanced nucleation of the carbon nanotube cap. The enhanced growth of the nanotube cap is explained by a nonthermal ion-induced graphene network restructuring mechanism.

Abstract: In this contribution, we present an overview of the various techniques for combining atomistic molecular dynamics with Monte Carlo simulations, mainly in the context of condensed matter systems, as well as a brief summary of the main accelerated dynamics techniques. Special attention is given to the force bias Monte Carlo technique and its combination with molecular dynamics, in view of promising recent developments, including a definable timescale. Various examples of the application of combined molecular dynamics / Monte Carlo simulations are given, in order to demonstrate the enhanced simulation efficiency with respect to either pure molecular dynamics or Monte Carlo.