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: 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: In this work, a time modulated RF driven DBD-like atmospheric pressure plasma jet in Ar + 2%O2, operating at a time averaged power of 6.5 W is investigated. Spatially resolved ozone densities and gas temperatures are obtained by UV absorption and Rayleigh scattering, respectively. Significant gas heating in the core of the plasma up to 700 K is found and at the position of this increased gas temperature a depletion of the ozone density is found. The production and destruction reactions of O3 in the jet effluent as a function of the distance from the nozzle are obtained from a zero-dimensional chemical kinetics model in plug flow mode which considers relevant air chemistry due to air entrainment in the jet fluent. A comparison of the measurements and the models show that the depletion of O3 in the core of the plasma is mainly caused by an enhanced destruction of O3 due to a large atomic oxygen density.

Abstract: Experiments have demonstrated the potential selective anticancer capacity of cold atmospheric plasmas (CAPs), but the underlying mechanisms remain unclear. Using computer simulations, we try to shed light on the mechanism of selectivity, based on aquaporins (AQPs), i.e. transmembrane protein channels transferring external H 2 O 2 and other reactive oxygen species, created e.g. by CAPs, to the cell interior. Specifically, we perform molecular dynamics simulations for the permeation of H 2 O 2 through AQP1 (one of the members of the AQP family) and the palmitoyl-oleoyl-phosphatidylcholine (POPC) phospholipid bilayer (PLB). The free energy barrier of H 2 O 2 across AQP1 is lower than for the POPC PLB, while the permeability coefficient, calculated using the free energy and diffusion rate profiles, is two orders of magnitude higher. This indicates that the delivery of H 2 O 2 into the cell interior should be through AQP. Our study gives a better insight into the role of AQPs in the selectivity of CAPs for treating cancer cells.

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: One of the main issues in plasma chemistry modeling is that the cross sections and rate coefficients are subject to uncertainties, which yields uncertainties in the modeling results and hence hinders the predictive capabilities. In this paper, we reveal the impact of these uncertainties on the model predictions of plasma-based dry reforming in a dielectric barrier discharge. For this purpose, we performed a detailed uncertainty analysis and sensitivity study. 2000 different combinations of rate coefficients, based on the uncertainty from a log-normal distribution, are used to predict the uncertainties in the model output. The uncertainties in the electron density and electron temperature are around 11% and 8% at the maximum of the power deposition for a 70% confidence level. Still, this can have a major effect on the electron impact rates and hence on the calculated conversions of CO2 and CH4, as well as on the selectivities of CO and H2. For the CO2 and CH4 conversion, we obtain uncertainties of 24% and 33%, respectively. For the CO and H2 selectivity, the corresponding uncertainties are 28% and 14%, respectively. We also identify which reactions contribute most to the uncertainty in the model predictions. In order to improve the accuracy and reliability of plasma chemistry models, we recommend using only verified rate coefficients, and we point out the need for dedicated verification experiments.

Abstract: The size and lattice constant evolution of Pb nanoparticles (NPs) synthesized by high fluence implantation in crystalline Si have been studied with a variety of experimental techniques. Results obtained from small-angle x-ray scattering showed that the Pb NPs grow with increasing implantation fluence and annealing duration. The theory of NP growth kinetics can be applied to qualitatively explain the size evolution of the Pb NPs during the implantation and annealing processes. Moreover, the lattice constant of the Pb NPs was evaluated by conventional x-ray diffraction. The lattice dilatation was observed to decrease with increasing size of the Pb NPs. Such lattice constant tuning can be attributed to the pseudomorphism caused by the lattice mismatch between the Pb NPs and the Si matrix.

Abstract: The characteristics of spin waves in ferromagnetic waveguides with non-uniform magnetization have been investigated for situations where the shape anisotropy field of the waveguide is comparable to the external bias field. Spin-wave generation was realized by the magnetoelastic effect by applying normal and shear strain components, as well as by the Oersted field emitted by an inductive antenna. The magnetoelastic excitation field has a non-uniform profile over the width of the waveguide because of the non-uniform magnetization orientation, whereas the Oersted field remains uniform. Using micromagnetic simulations, we indicate that both types of excitation fields generate quantised width modes with both odd and even mode numbers as well as tilted phase fronts. We demonstrate that these effects originate from the average magnetization orientation with respect to the main axes of the magnetic waveguide. Furthermore, it is indicated that the excitation efficiency of the second-order mode generally surpasses that of the first-order mode due to their symmetry. The relative intensity of the excited modes can be controlled by the strain state as well as by tuning the dimensions of the excitation area. Finally, we demonstrate that the nonreciprocity of spin-wave radiation due to the chirality of an Oersted field generated by an inductive antenna is absent for magnetoelastic spin-wave excitation.

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: 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: Experiments show that the etch rate of Si with SF6 inductively coupled plasma (ICP) is significantly influenced by the absolute gas flow rate in the range of 50–600 sccm, with a maximum at around 200 sccm. Therefore, we numerically investigate the effects of the gas flow rate on the bulk plasma properties and on the etch rate, to obtain more insight in the underlying reasons of this effect. A hybrid Monte Carlo—fluid model is applied to simulate an SF6 ICP. It is found that the etch rate is influenced by two simultaneous effects: (i) the residence time of the gas and (ii) the temperature profile of the plasma in the ICP volume, resulting indeed in a maximum etch rate at 200 sccm.

Abstract: In this paper, a simulation method is described to predict the etching behaviour of Cl2/O2/Ar inductively coupled plasmas on a Si substrate, as used in shallow trench isolation for the production of electronic devices. The hybrid plasma equipment model (HPEM) developed by Kushner et al is applied to calculate the plasma characteristics in the reactor chamber and two additional Monte Carlo simulations are performed to predict the fluxes, angles and energy of the plasma species bombarding the Si substrate, as well as the resulting surface processes such as etching and deposition. The simulations are performed for a wide variety of operating conditions such as gas composition, chamber pressure, power deposition and substrate bias. It is predicted by the simulations that when the fraction of oxygen in the gas mixture is too high, the oxidation of the Si substrate is superior to the etching of Si by chlorine species, resulting in an etch rate close to zero as is also observed in the experiments.

Abstract: A modelling study is performed of high-density low-pressure inductively coupled CF4/CHF3/H2/Cl2/O2/HBr plasmas under different gas mixing ratios. A reaction set describing the complete plasma chemistry is presented and discussed. The gas fraction of each component in this mixture is varied to investigate the sensitivity of the plasma properties, like electron density, plasma potential and species densities, towards the gas mixing ratios. This research is of great interest for microelectronics applications because these gases are often combined in two (or more)-component mixtures, and mixing gases or changing the fraction of a gas can sometimes yield unwanted reaction products or unexpected changes in the overall plasma properties due to the increased chemical complexity of the system. Increasing the CF4 fraction produces more F atoms for chemical etching as expected, but also more prominently lowers the density of Cl atoms, resulting in an actual drop in the etch rate under certain conditions. Furthermore, CF4 decreases the free electron density when mixed with Cl2. However, depending on the other gas components, CF4 gas can also sometimes enhance free electron density. This is the case when HBr is added to the mixture. The addition of H2 to the gas mixture will lower the sputtering process, not only due to the lower overall positive ion density at higher H2 fractions, but also because more H+, H2 + and H3 + are present and they have very low sputter yields. In contrast, a larger Cl2 fraction results in more chemical etching but also in less physical sputtering due to a smaller abundance of positive ions. Increasing the O2 fraction in the plasma will always lower the etch rate due to more oxidation of the wafer surface and due to a lower plasma density. However, it is also observed that the density of F atoms can actually increase with rising O2 gas fraction. This is relevant to note because the exact balance between fluorination and oxidation is important for fine-tuning the overall etch rate and for control of the sidewall profile. Finally, HBr is often used as a chemical etcher, but when mixed with F- or Cl-containing gases, HBr creates the same diluting effects as Ar or He, because a
higher fraction results in less chemical etching but more (physical) sputtering.

Abstract: Non-small cell lung cancer (NSCLC) is the most common type of lung cancer, with 85% of all lung cancer reported as NSCLC. Moreover, there are no effective treatments in advanced NSCLC. This study shows for the first time that oral administration of plasma-treated water (PTW) can cure advanced NSCLC. The cold plasma in water generates a cocktail of reactive species, and oral administration of this cocktail to mice showed no toxicities even at the highest dose of PTW, after a single dose and repeated doses for 28 d in mice. In vivo studies reveal that PTW showed favorable anticancer effects on chemo-resistant lung cancer, similarly to gefitinib treatment as a reference drug in a chemo-resistant NSCLC model. The anticancer activities of PTW seem to be involved in inhibiting proliferation and angiogenesis and enhancing apoptosis in the cancer cells. Interestingly, the PTW contributes to enhanced immune response and improved cachexia in the model.

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