Abstract: The continuous miniaturization of nanodevices, such as transistors, solar cells, and optical fibers, requires the controlled synthesis of (ultra)thin gate oxides (<10 nm), including Si gate-oxide (SiO2) with high quality at the atomic scale. Traditional thermal growth of SiO2 on planar Si surfaces, however, does not allow one to obtain such ultrathin oxide due to either the high oxygen diffusivity at high temperature or the very low sticking ability of incident oxygen at low temperature. Two recent techniques, both operative at low (room) temperature, have been put forward to overcome these obstacles: (i) hyperthermal oxidation of planar Si surfaces and (ii) thermal or plasma-assisted oxidation of nonplanar Si surfaces, including Si nanowires (SiNWs). These nanooxidation processes are, however, often difficult to study experimentally, due to the key intermediate processes taking place on the nanosecond time scale.
In this Account, these Si nano-oxidation techniques are discussed from a computational point of view and compared to both hyperthermal and thermal oxidation experiments, as well as to well-known models of thermal oxidation, including the Deal−Grove, Cabrera−Mott, and Kao models and several alternative mechanisms. In our studies, we use reactive molecular dynamics (MD) and hybrid MD/Monte Carlo simulation techniques, applying the Reax force field. The incident energy of oxygen species is chosen in the range of 1−5 eV in hyperthermal oxidation of planar Si surfaces in order to prevent energy-induced damage. It turns out that hyperthermal growth allows for two growth modes, where the ultrathin oxide thickness depends on either (1) only the kinetic energy of the incident oxygen species at a growth temperature below Ttrans = 600 K, or (2) both the incident energy and the growth temperature at a growth temperature above Ttrans. These modes are specific to such ultrathin oxides, and are not observed in traditional thermal oxidation, nor theoretically considered by already existing models. In the case of thermal or plasma-assisted oxidation of small Si nanowires, on the other hand, the thickness of the ultrathin oxide is a function of the growth temperature and the nanowire diameter. Below Ttrans, which varies with the nanowire diameter, partially oxidized SiNW are formed, whereas complete oxidation to a SiO2 nanowire occurs only above Ttrans. In both nano-oxidation processes at lower temperature (T < Ttrans), final sandwich c-Si|SiOx|a-SiO2 structures are obtained due to a competition between overcoming the energy barrier to penetrate into Si subsurface layers and the compressive stress (∼2−3 GPa) at the Si crystal/oxide interface. The overall atomic-simulation results strongly indicate that the thickness of the intermediate SiOx (x < 2) region is very limited (∼0.5 nm) and constant irrespective of oxidation parameters. Thus, control over the ultrathin SiO2 thickness with good quality is indeed possible by accurately tuning the oxidant energy, oxidation temperature and surface curvature.
In general, we discuss and put in perspective these two oxidation mechanisms for obtaining controllable ultrathin gate-oxide films, offering a new route toward the fabrication of nanodevices via selective nano-oxidation.

Abstract: Microdischarges are an important type of plasma discharge that possess several unique characteristics, such as the presence of a stable glow discharge, high plasma density and intense excimer radiation, leading to several potential applications. The intense and controllable gas heating within the extremely small dimensions of microdischarges has been exploited in microthruster technologies by incorporating a micro-nozzle to generate the thrust. This kind of microthruster has a significantly improved specific impulse performance compared to conventional cold gas thrusters, and can meet the requirements arising from the emerging development and application of micro-spacecraft. In this paper, we performed a self-consistent 2D particle-in-cell simulation, with a Monte Carlo collision model, of a microdischarge operating in a prototype micro-plasma thruster with a hollow cylinder geometry and a divergent micro-nozzle. The model takes into account the thermionic electron emission including the Schottky effect, the secondary electron emission due to cathode bombardment by the plasma ions, several different collision processes, and a non-uniform argon background gas density in the cathode-anode gap. Results in the high-pressure (several hundreds of Torr), high-current (mA) operating regime showing the behavior of the plasma density, potential distribution, and energy flux towards the hollow cathode and anode are presented and discussed. In addition, the results of simulations showing the effect of different argon gas pressures, cathode material work function and discharge voltage on the operation of the microdischarge thruster are presented. Our calculated properties are compared with experimental data under similar conditions and qualitative and quantitative agreements are reached.

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: 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: In recent years, the interest in treating cancer cells with plasma treated media (PTM) and plasma treated water (PTW) has increased tremendously. However, the actions of PTM and PTW are still not entirely understood. For instance, it is not clear whether the action of PTM is due to a modification in proteins/amino acids after plasma treatment of the media, or due to reactive oxygen and nitrogen species (RONS) generated from the plasma, or a combination of both effects. To differentiate between the actions of RONS and modified proteins/amino acids on the treatment of cancer cells, we compared the effects of PTM and PTW on two different pancreatic ductal adenocarcinomas (MiaPaca-2, BxPc3) and pancreatic stellate cells
(PSCs) (hPSC128-SV). PSCs closely interact with cancer cells to create a tumor-promoting environment that stimulates local tumor progression and metastasis. We treated culture media and deionized water with a cold atmospheric plasma (CAP) jet, and subsequently applied this PTM/PTW at various ratios to the pancreatic cancer and PSC cell lines. We evaluated cell death, intracellular ROS concentrations and the mRNA expression profiles of four oxidative stress-related genes, i.e. Mitogen-activated protein kinase 7 (MAPK7), B-cell lymphoma 2 (BCL2), Checkpoint kinase 1 (CHEK1) and DNA damage-inducible transcript 3, also known as C/EBP homologous protein (CHOP). Our findings demonstrate that PTM and PTW have a similar efficacy to kill pancreatic cancer cells, while PTW is slightly more effective in killing PSCs, as compared to PTM. Furthermore, we observed an enhancement of the intracellular ROS concentrations in both pancreatic cancer cells and PSCs. Thus, it is likely that under our experimental conditions, the anti-cancer activity of PTM can be attributed more to the RONS present in the treated liquid, than to the modification of proteins/amino acids in the media. Furthermore, the fact that the chemo-resistant PSCs were killed by PTM/PTW may offer possibilities for new anti-cancer therapies for pancreatic cancer cells, including PSCs.

Abstract: A gliding arc plasmatron (GAP), which is very promising for purification and gas conversion,
is characterized in nitrogen using optical emission spectroscopy and high-speed photography,
because the cross sections of electron impact excitation of N 2 are well known. The gas
temperature (of about 5500 K), the electron density (up to 1.5 × 10 15 cm −3 ) and the reduced
electric field (of about 37 Td) are determined using an absolutely calibrated intensified charge-
coupled device (ICCD) camera, equipped with an in-house made optical arrangement for
simultaneous two-wavelength diagnostics, adapted to the transient behavior of a GA channel
in turbulent gas flow. The intensities of nitrogen molecular emission bands, N 2 (C–B,0–0) as
well as N +
2 (B–X,0–0), are measured simultaneously. The electron density and the reduced
electric field are determined at a spatial resolution of 30 µm, using numerical simulation and
measured emission intensities, applying the Abel inversion of the ICCD images. The temporal
behavior of the GA plasma channel and the formation of plasma plumes are studied using a
high-speed camera. Based on the determined plasma parameters, we suggest that the plasma
plume formation is due to the magnetization of electrons in the plasma channel of the GAP by
an axial magnetic field in the plasma vortex.

Abstract: Recent advances in biomedical research in cancer immunotherapy have identified the use of an oxidative stress-based approach to treat cancers, which works by inducing immunogenic cell death (ICD) in cancer cells. Since the anti-cancer effects of non-thermal plasma (NTP) are largely attributed to the reactive oxygen and nitrogen species that are delivered to and generated inside the target cancer cells, it is reasonable to postulate that NTP would be an effective modality for ICD induction. NTP treatment of tumors has been shown to destroy cancer cells rapidly and, under specific treatment regimens, this leads to systemic tumorspecific immunity. The translational benefit of NTP for treatment of cancer relies on its ability to enhance the interactions between NTP-exposed minor cells and local immune cells which initiates subsequent protective immune responses. This review discusses results from recent investigations of NTP application to induce ICD in cancer cells. With further optimization of clinical devices and treatment protocols, NTP can become an essential part of the therapeutic armament against cancer.

Abstract: Plasma catalysis is gaining increasing interest for various gas conversion applications, such as CO2 conversion into value-added chemicals and fuels, CH4 activation into hydrogen, higher hydrocarbons or oxygenates, and NH3 synthesis. Other applications are already more established, such as for air pollution control, e.g. volatile organic compound remediation, particulate matter and NOx removal. In addition, plasma is also very promising for catalyst synthesis and treatment. Plasma catalysis clearly has benefits over ‘conventional’ catalysis, as outlined in the Introduction. However, a better insight into the underlying physical and chemical processes is crucial. This can be obtained by experiments applying diagnostics, studying both the chemical processes at the catalyst surface and the physicochemical mechanisms of plasma-catalyst interactions, as well as by computer modeling. The key challenge is to design cost-effective, highly active and stable catalysts tailored to the plasma environment. Therefore, insight from thermal catalysis as well as electro- and photocatalysis is crucial. All these aspects are covered in this Roadmap paper, written by specialists in their field, presenting the state-of-the-art, the current and future challenges, as well as the advances in science and technology needed to meet these challenges.

Abstract: Dielectric barrier discharges (DBDs) typically operate in the filamentary regime and thus exhibit great spatial and temporal non-uniformity. In order to optimize DBDs for various applications, such as in plasma catalysis, more fundamental insight is needed. Here, we consider how the millions of microdischarges, characteristic for a DBD, influence individual gas molecules. We use a Monte Carlo approach to determine the number of microdischarges to which a single molecule would be exposed, by means of particle tracing simulations through a full-scale packed bed DBD reactor, as well as an empty DBD reactor. We find that the fraction of microdischarges to which the molecules are exposed can be approximated as the microdischarge volume over the entire reactor gas volume. The use of this concept provides good agreement between a plasma-catalytic kinetics model and experiments for plasma-catalytic NH3 synthesis. We also show that the concept of the fraction of microdischarges indicates the efficiency by which the plasma power is transferred to the gas molecules. This generalised concept is also applicable for other spatially and temporally non-uniform plasmas.

Abstract: Packed bed dielectric barrier discharge (PB-DBD) plasma reactors are very promising for various plasma catalysis applications, but the exact mechanisms of plasma-catalyst interaction are far from understood, because the plasma discharge and catalyst/packing properties are mutually dependent. To better understand the effect of packing dielectric material on the electrical plasma properties, we study here a single bead DBD plasma reactor operating in dry air, with beads of different dielectric constant and for different applied voltages, by means of fluid modelling validated by optical imaging experiments. Our study reveals that the plasma in the single bead DBD reactor can manifest itself in two different modalities, i.e. (a) polar discharges at the bead poles in contact with the electrodes, and (b) a streamer discharge caused by surface ionization waves, which bridges the gas gap. Beads with high dielectric constant result in localised electric field enhancement and hence yield a reduction of the applied voltage required for plasma production. At low applied voltage, the discharge appears as polar discharges between the bead and the electrodes, and upon higher voltage it undergoes a transition into a bridging streamer discharge. The transition voltage to the streamer mode rises for beads with higher dielectric constant. These observations are important for plasma catalysis applications. A higher dielectric constant yields a higher electric field and thus higher average electron energy and density, giving rise to more reactive species, but it also yields a confined discharge near the contact points of packing beads, limiting the interaction area between the catalyst and the active plasma species. In addition, our model reveals that the dielectric bead behaves as a capacitor and traps charges, which can explain the significant occurrence of partial discharging in PB-DBDs and non-parallelogram shaped Lissajous plots. Hence, equivalent circuit modelling of PB-DBDs should take into account the role of packing beads in charge trapping as a capacitor.