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: tWe investigate microdischarge formation inside catalyst pores by a two-dimensional fluid model forvarious pore sizes in the m-range and for various applied voltages. Indeed, this is a poorly understoodphenomenon in plasma catalysis. The calculations are performed for a dielectric barrier discharge inhelium, at atmospheric pressure. The electron and ion densities, electron temperature, electric field andpotential, as well as the electron impact ionization and excitation rate and the densities of excited plasmaspecies, are examined for a better understanding of the characteristics of the plasma inside a pore. Theresults indicate that the pore size and the applied voltage are critical parameters for the formation of amicrodischarge inside a pore. At an applied voltage of 20 kV, our calculations reveal that the ionizationmainly takes place inside the pore, and the electron density shows a significant increase near and inthe pore for pore sizes larger than 200m, whereas the effect of the pore on the total ion density isevident even for 10m pores. When the pore size is fixed at 30m, the presence of the pore has nosignificant influence on the plasma properties at an applied voltage of 2 kV. Upon increasing the voltage,the ionization process is enhanced due to the strong electric field and high electron temperature, andthe ion density shows a remarkable increase near and in the pore for voltages above 10 kV. These resultsindicate that the plasma species can be formed inside pores of structured catalysts (in the m range),and they may interact with the catalyst surface, and affect the plasma catalytic process.

Abstract: We apply an implicit particle-in-cell Monte-Carlo (PIC-MC) method to study a radio-frequency argon microdischarge at steady state in the glow discharge limit, in which the microdischarge is sustained by secondary electron emission from the electrodes. The plasma density, electron energy distribution function (EEDF), and electron temperature are calculated in a wide range of operating conditions, including driving voltage, microdischarge gap, and pressure. Also, the effect of gap size scaling (in the range of 50-1000 μm) on the plasma sustaining voltage and peak electron density at atmospheric pressure is examined, which has not been explored before. In our simulations, three different EEDFs, i.e., a so-called three temperature hybrid mode, a two temperature α mode, and a two temperature γ mode distribution, are identified at different gaps and voltages. The maximum sustaining voltage to avoid a transition from the glow mode to an arc is predicted, as well as the minimum sustaining voltage for a steady glow discharge. Our calculations elucidate that secondary electrons play an essential role in sustaining the discharge, and as a result the relationship between breakdown voltage and gap spacing is far away from the Paschen law at atmospheric pressure.

Abstract: Continued miniaturization of metal-oxide-semiconductor field-effect transistors (MOSFETs) requires an ever-decreasing thickness of the gate oxide. The structure of ultrathin silicon oxide films, however, critically depends on the oxidation mechanism. Using reactive atomistic simulations, we here demonstrate how the oxidation mechanism in hyperthermal oxidation of such structures may be controlled by the oxidation temperature and the oxidant energy. Specifically, we study the interaction of hyperthermal oxygen with energies of 15 eV with thin SiOx (x ≤ 2) films with a native oxide thickness of about 10 Å. We analyze the oxygen penetration depth probability and compare with results of the hyperthermal oxidation of a bare Si(100){2 × 1} (c-Si) surface. The temperature-dependent oxidation mechanisms are discussed in detail. Our results demonstrate that, at low (i.e., room) temperature, the penetrated oxygen mostly resides in the oxide region rather than at the SiOx|c-Si interface. However, at higher temperatures, starting at around 700 K, oxygen atoms are found to penetrate and to diffuse through the oxide layer followed by reaction at the c-Si boundary. We demonstrate that hyperthermal oxidation resembles thermal oxidation, which can be described by the DealGrove model at high temperatures. Furthermore, defect creation mechanisms that occur during the oxidation process are also analyzed. This study is useful for the fabrication of ultrathin silicon oxide gate oxides for metal-oxide-semiconductor devices as it links parameters that can be straightforwardly controlled in experiment (oxygen temperature, velocity) with the silicon oxide structure.

Abstract: In this study, we investigated the diffusion of H-atoms to the subsurface and their further diffusion into the bulk of a Ni(111) crystal by means of density functional theory calculations in the context of thermal and plasma-assisted catalysis. The H-atoms at the surface can originate from the dissociative adsorption of H2 or CH4 molecules, determining the surface H-coverage. When a threshold H-coverage is passed, corresponding to 1.00 ML for the crystalline Ni(111) surface, the surface-bound H-atoms start to diffuse to the subsurface. A similar threshold coverage is observed for the interstitial H-coverage. Once the interstitial sites are filled up with a coverage above 1.00 ML of H, dissolution of interstitial H-atoms to the layer below the interstitial sites will be initiated. Hence, by applying a high pressure or inducing a reactive plasma and high temperature, increasing the H-flux to the surface, a large amount of hydrogen can diffuse in a crystalline metal like Ni and can be absorbed. The formation of metal hydride may modify the entire reaction kinetics of the system. Equivalently, the H-atoms in the bulk can easily go back to the surface and release a large amount of heat. In a plasma process, H-atoms are formed in the plasma, and therefore the energy barrier for dissociative adsorption is dismissed, thus allowing achievement of the threshold coverage without applying a high pressure as in a thermal process. As a result, depending on the crystal plane and type of metal, a large number of H-atoms can be dissolved (absorbed) in the metal catalyst, explaining the high efficiency of plasma-assisted catalytic reactions. Here, the mechanism of H-dissolution is established as a chemical identifier for the investigation of the reaction kinetics of a chemical process.

Abstract: We study the mechanisms of CO2 conversion in a nanosecond repetitively pulsed (NRP) discharge, by means of a chemical kinetics model. The calculated conversions and energy efficiencies are in reasonable agreement with experimental results over a wide range of specific energy input values, and the same applies to the evolution of gas temperature and CO2 conversion as a function of time in the afterglow, indicating that our model provides a realistic picture of the underlying mechanisms in the NRP discharge and can be used to identify its limitations and thus to suggest further improvements. Our model predicts that vibrational excitation is very important in the NRP discharge, explaining why this type of plasma yields energy-efficient CO2 conversion. A significant part of the CO2 dissociation occurs by electronic excitation from the lower vibrational levels toward repulsive electronic states, thus resulting in dissociation. However, vibration−translation (VT) relaxation (depopulating the higher vibrational levels) and CO + O recombination (CO + O + M → CO2 + M), as well as mixing of the converted gas with fresh gas entering the plasma in between the pulses, are limiting factors for the conversion and energy efficiency. Our model predicts that extra cooling, slowing down the rate of VT relaxation and of the above recombination reaction, thus enhancing the contribution of the highest vibrational levels to the overall CO2 dissociation, can further improve the performance of the NRP discharge for energy-efficient CO2 conversion.

Abstract: Cold atmospheric plasma and nanomedicine originally emerged as individual domains, but are increasingly applied in combination with each other. Most research is performed in the context of cancer treatment, with only little focus yet on the possible synergies. Many questions remain on the potential of this promising hybrid technology, particularly regarding regenerative medicine and tissue engineering. In this perspective article, we therefore start from the fundamental mechanisms in the individual technologies, in order to envision possible synergies for wound healing and tissue recovery, as well as research strategies to discover and optimize them. Among these strategies, we demonstrate how cold plasmas and nanomaterials can enhance each other’s strengths and overcome each other’s limitations. The parallels with cancer research, biotechnology and plasma surface modification further serve as inspiration for the envisioned synergies in tissue regeneration. The discovery and optimization of synergies may also be realized based on a profound understanding of the underlying redox- and field-related biological processes. Finally, we emphasize the toxicity concerns in plasma and nanomedicine, which may be partly remediated by their combination, but also partly amplified. A widespread use of standardized protocols and materials is therefore strongly recommended, to ensure both a fast and safe clinical implementation.

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.

Abstract: The initial stage of nanoparticle formation and growth in radiofrequency acetylene (C2H2) plasmas is investigated by means of a self-consistent one-dimensional fluid model. A detailed chemical kinetic scheme, containing electron impact, ion-neutral, and neutral-neutral reactions, has been developed in order to predict the underlying dust growth mechanisms and the most important dust precursors. The model considers 41 different species (neutrals, radicals, ions, and electrons) describing hydrocarbons (CnHm) containing up to 12 carbon atoms. Possible routes for particle growth are discussed. Both positive and negative ion reaction pathways are considered, as consecutive anion- and cation-molecule reactions seem to lead to a fast build up of the carbon skeleton.

Abstract: The reproducible temporal separation of ion signals generated from a single multi-element droplet, observed in previous studies, was investigated in detail in this work using an ICPTOFMS with high temporal resolution. It was shown that the signal peak intensities of individual elements temporally shift relative to each other only for droplets moving through the plasma off-axis. The magnitude of these shifts correlated with the vaporization temperatures of the analytes and depended on the radial position of the droplets as well as on the thermal properties and velocity profiles of the carrier gases of the ICP. The occurrence of the signal shifting was explained by a spatial separation of analytes already present in the vapor phase in the ICP from a yet unvaporized residue of the droplet. This separation is most likely driven by anisotropic diffusion of vaporized analytes towards the plasma axis and a radial velocity gradient. The proposed explanation is supported by modeling of the gas velocities inside the ICP and imaging of the atomic and ionic emissions produced from single droplets, whose patterns were sloping towards the center of the torch. The effects observed in these studies are important not only for the fundamental understanding of analyteplasma interactions but have also a direct impact on the signal intensities and stability.

Abstract: The electron bounce resonance heating (BRH) in dual-frequency capacitively coupled plasmas operated in oxygen is studied by different experimental methods and a particle-in-cell/Monte Carlo collision (PIC/MCC) simulation, and compared with the electropositive argon discharge. In comparison with argon, the experimental results show that in an oxygen discharge the resonance peaks in positive-ion density and light intensity tend to occur at larger electrode gaps. Moreover, at electrode gaps L > 2.5 cm, the positive-ion (and electron) density and the light emission drop monotonically in the oxygen discharge upon increasing L, whereas they rise (after an initial drop) in the argon case. At resonance gap the electronegativity reaches its maximum due to the BRH. All these experimental observations are explained by PIC/MCC simulations, which show that in the oxygen discharge the bulk electric field becomes quite strong and is out of phase with the sheath field. Therefore, it retards the resonance electrons when traversing the bulk, resulting in a suppressed BRH. Both experiment and simulation results show that this effect becomes more pronounced at lower high-frequency power, when the discharge mode changes from electropositive to electronegative. In a pure oxygen discharge, the BRH is suppressed with increasing pressure and almost diminishes at 12 Pa. Finally, the driving frequency significantly affects the BRH, because it determines the phase relation between bulk electric field and sheath electric field.

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: Carbon nanotubes (CNTs) are nowadays routinely grown in a thermal CVD setup. State-of-the-art plasma-enhanced CVD (PECVD) growth, however, offers advantages over thermal CVD. A lower growth temperature and the growth of aligned freestanding single-walled CNTs (SWNTs) makes the technique very attractive. The atomic scale growth mechanisms of PECVD CNT growth, however, remain currently entirely unexplored. In this contribution, we employed molecular dynamics simulations to focus on the effect of applying an electric field on the SWNT growth process, as one of the effects coming into play in PECVD. Using sufficiently strong fields results in (a) alignment of the growing SWNTs, (b) a better ordering of the carbon network, and (c) a higher growth rate relative to thermal growth rate. We suggest that these effects are due to the small charge transfer occurring in the Ni/C system. These simulations constitute the first study of PECVD growth of SWNTs on the atomic level.