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 numerical parameter study has been performed for a cylindrical atmospheric pressure dielectric barrier discharge (DBD) in helium with nitrogen impurities using a two-dimensional time-dependent fluid model. The calculated electric currents and gap voltages as a function of time for a given applied potential are presented, as well as the number densities of the various plasma species. This study shows that for the geometry under consideration the applied voltage parameters have a large impact on the electric current profiles and that the discharge current is always determined by the electron and ion conduction currents while the displacement current is nearly negligible. A relative broadening of the current profiles (compared with the duration of the half cycle of the applied voltage) with an increase in the applied frequency is obtained. Nearly sinusoidal current wave forms, usually typical for radio frequency DBDs, are observed while still operating at the frequencies of tens of kilohertz. For the setup under investigation, the Townsend mode of the DBD is observed in the entire range of applied voltage amplitudes and frequencies. It is shown that the average power density dissipated in the discharge increases with rising applied voltage and frequency. An increase in applied voltage frequency leads to an increase in the electron density and a decrease in electron energy, while increasing the voltage amplitude has the opposite effect.

Abstract: A two-dimensional fluid model, including the full set of Maxwell equations, has been developed and applied to investigate the effect of a phase shift between two power sources on the radial uniformity of several plasma characteristics in a hydrogen capacitively coupled plasma. This study was carried out at various frequencies in the range 13.56200 MHz. When the frequency is low, at 13.56 MHz, the plasma density is characterized by an off-axis peak when both power sources are in-phase (phgr = 0), and the best radial uniformity is obtained at phgr = π. This trend can be explained because the radial nonuniformity caused by the electrostatic edge effect can be effectively suppressed by the phase-shift effect at a phase difference equal to π. When the frequency rises to 60 MHz, the plasma density profiles shift smoothly from edge-peaked over uniform to centre-peaked as the phase difference increases, due to the pronounced standing-wave effect, and the best radial uniformity is reached at phgr = 0.3π. At a frequency of 100 MHz, a similar behaviour is observed, except that the maximum of the plasma density moves again towards the radial edge at the reverse-phase case (phgr = π), because of the dominant skin effect. When the frequency is 200 MHz, the bulk plasma density increases significantly with increasing phase-shift values, and a better uniformity is obtained at phgr = 0.4π. This is because the density in the centre increases faster than at the radial edge as the phase difference rises, due to the increasing power deposition Pz in the centre and the decreasing power density Pr at the radial edge. As the phase difference increases to π, the maximum near the radial edge becomes obvious again. This is because the skin effect has a predominant influence on the plasma density under this condition, resulting in a high density at the radial edge. Moreover, the axial ion flux increases monotonically with phase difference, and exhibits similar profiles to the plasma density. The calculation results illustrate that the radial uniformity of the various plasma characteristics is strongly dependent on the applied frequency and the phase shift between both power sources, which is important to realize, for controlling the uniformity of the plasma etch and deposition processes.

Abstract: Plasma medicine has been claimed to provide a novel route to heal wounds and regenerate skin, although very little is currently known about the elementary processes taking place. We carried out a series of ReaxFF-based reactive molecular dynamics simulations to investigate the interaction of O and OH radicals with lipids, more specifically with α-linolenic acid as a model for the free fatty acids present in the upper skin layer. Our calculations predict that the O and OH radicals most typically abstract a H atom from the fatty acids, which can lead to the formation of a conjugated double bond, but also to the incorporation of alcohol or aldehyde groups, thereby increasing the hydrophilic character of the fatty acids and changing the general lipid composition of the skin. Within the limitations of the investigated model, no formation of possibly toxic products was observed.

Abstract: A parameter study is carried out for an inductively coupled plasma used for the synthesis of carbon nanotubes or carbon nanofibres (CNTs/CNFs), by means of the Hybrid Plasma Equipment Model. The influence of processing parameters including gas ratio for four different gas mixtures typically used for CNT/CNF growth (i.e. CH4/H2, CH4/NH3, C2H2/H2 and C2H2/NH3), inductively coupled plasma (ICP) power (501000 W), operating pressure (10 mTorr1 Torr), bias power (01000 W) and temperature of the substrate (01000 °C) on the plasma chemistry is investigated and the optimized conditions for CNT/CNF growth are analysed. Summarized, our calculations suggest that a lower fraction of hydrocarbon gases (CH4 or C2H2, i.e. below 20%) and hence a higher fraction of etchant gases (H2 or NH3) in the gas mixture result in more 'clean' conditions for controlled CNT/CNF growth. The same applies to a higher ICP power, a moderate ICP gas pressure above 100 mTorr (at least for single-walled carbon nanotubes), a high bias power (for aligned CNTs) and an intermediate substrate temperature.

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: In this paper, the initial mechanisms of nanoparticle formation and growth in radiofrequency acetylene (C2H2) plasmas are investigated by means of a comprehensive self-consistent one-dimensional (1D) fluid model. This model is an extension of the 1D fluid model, developed earlier by De Bleecker et al. Based on the comparison of our previous results with available experimental data for acetylene plasmas in the literature, some new mechanisms for negative ion formation and growth are proposed. Possible routes are considered for the formation of larger (linear and branched) hydrocarbons C2nH2 (n = 3, 4, 5), which contribute to the generation of C2nH− anions (n = 3, 4, 5) due to dissociative electron attachment. Moreover, the vinylidene anion (H2CC−) and higher anions (n = 24) are found to be important plasma species.

Abstract: A hybrid model, called the hybrid plasma equipment model, was used to study inductively coupled SF6 plasmas used for Si etching applications. The plasma properties such as number densities of electrons, positive and negative ions, and neutrals are calculated under typical etching conditions. The electron kinetics is analysed by means of the electron energy probability function. The plasma chemistry taking place in pure SF6 and in an Ar/SF6 mixture is also discussed, and finally the effect of the argon fraction on the plasma properties is investigated.

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: 'Bohm diffusion' causes the electrons to diffuse perpendicularly to the magnetic field lines. However, its origin is not yet completely understood: low and high frequency electric field fluctuations are both named to cause Bohm diffusion. The importance of including this process in a Monte Carlo (MC) model is demonstrated by comparing calculated ionization rates with particle-in-cell/Monte Carlo collisions (PIC/MCC) simulations. A good agreement is found with a Bohm diffusion parameter of 0.05, which corresponds well to experiments. Since the PIC/MCC method accounts for fast electric field fluctuations, we conclude that Bohm diffusion is caused by fast electric field phenomena.

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