Abstract: We studied the evolution of an isolated pulsed plasma in a vortex flow stabilised microwave (MW) discharge in N2 at 25 mbar via the combination of 0D kinetics modelling, iCCD imaging and laser scattering diagnostics. Quenching of electronically excited N2 results in fast gas heating and the onset of a thermal-ionisation instability, contracting the discharge volume. The onset of a thermal-ionisation instability driven by vibrational excitation pathways is found to facilitate significantly higher N2 conversion (i.e. dissociation to atomic N2 ) compared to pre-instability conditions, emphasizing the potential utility of this dynamic in future fixation applications. The instability onset is found to be instigated by super-elastic heating of the electron energy distribution tail via vibrationally excited N2 . Radial contraction of the discharge to the skin depth is found to occur post instability, while the axial elongation is found to be temporarily contracted during the thermal instability onset. An increase in power reflection during the thermal instability onset eventually limits the destabilising effects of exothermic electronically excited N2 quenching. Translational and vibrational temperature reach a quasi-non-equilibrium after the discharge contraction, with translational temperatures reaching ∼1200 K at the pulse end, while vibrational temperatures are found in near equilibrium with the electron energy (1 eV, or ∼11 600 K). This first description of the importance of electronically excited N2 quenching in thermal instabilities gives an additional fundamental understanding of N2 plasma behaviour in pulsed MW context, and thereby brings the eventual implementation of this novel N2 fixation method one step closer.

Abstract: Plasma-enhanced chemical deposition (PECVD) of carbon nanotubes has already been shown to allow chirality control to some extent. In PECVD, however, etching may occur simultaneously with the growth, and the occurrence of intermediate processes further significantly complicates the growth process.
We here employ a computational approach with experimental support to study the plasma-based formation of Ni nanoclusters, Ni-catalyzed CNT growth and subsequent etching processes, in order to understand the underpinning nanoscale mechanisms. We find that hydrogen is the dominant factor in both the re-structuring of a Ni film and the subsequent appearance of Ni nanoclusters, as well as in the CNT nucleation and etching processes. The obtained results are compared with available theoretical and experimental studies and provide a deeper understanding of the occurring nanoscale mechanisms in plasma-assisted CNT nucleation and growth.

Abstract: Atomic scale simulations of the nucleation and growth of carbon nanotubes is essential for understanding their growth mechanism. In spite of over twenty years of simulation efforts in this area, limited progress has so far been made on addressing the role of the hydrocarbon growth precursor. Here we report on atomic scale simulations of cap nucleation of single-walled carbon nanotubes from hydrocarbon precursors. The presented mechanism emphasizes the important role of hydrogen in the nucleation process, and is discussed in relation to previously presented mechanisms. In particular, the role of hydrogen in the appearance of unstable carbon structures during in situ experimental observations as well as the initial stage of multi-walled carbon nanotube growth is discussed. The results are in good agreement with available experimental and quantum-mechanical results, and provide a basic understanding of the incubation and nucleation stages of hydrocarbon-based CNT growth at the atomic level.

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: The application of coreshell SiSiO2 nanowires as nanoelectronic devices strongly depends on their structure, which is difficult to tune precisely. In this work, we investigate the formation of the coreshell nanowires at the atomic scale, by reactive molecular dynamics simulations. The occurrence of two temperature-dependent oxidation mechanisms of ultra-small diameter Si-NWs is demonstrated. We found that control over the Si-core radius and the SiOx (x ≤ 2) oxide shell is possible by tuning the growth temperature and the initial Si-NW diameter. Two different structures were obtained, i.e., ultrathin SiO2 silica nanowires at high temperature and Si core|ultrathin SiO2 silica nanowires at low temperature. The transition temperature is found to linearly decrease with the nanowire curvature. Finally, the interfacial stress is found to be responsible for self-limiting oxidation, depending on both the initial Si-NW radius and the oxide growth temperature. These novel insights allow us to gain control over the exact morphology and structure of the wires, as is needed for their application in nanoelectronics.

Abstract: Device performance of Si|SiOx core-shell based nanowires critically depends on the exact control over the oxide thickness. Low-temperature plasma oxidation is a highly promising alternative to thermal oxidation allowing for improved control over the oxidation process, in particular for ultrasmall Si nanowires. We here elucidate the room temperature plasma oxidation mechanisms of ultrasmall Si nanowires using hybrid molecular dynamics / force-bias Monte Carlo simulations. We demonstrate how the oxidation and concurrent water formation mechanisms are a function of the oxidizing plasma species and we demonstrate how the resulting core-shell oxide thickness can be controlled through these species. A new mechanism of water formation is discussed in detail. The results provide a detailed atomic level explanation of the oxidation process of highly curved Si surfaces. These results point out a route toward plasma-based formation of ultrathin core-shell Si|SiOx nanowires at room temperature.

Abstract: This study reports on the possible effects of OH radical impact on the transmembrane domain 6 of P-glycoprotein, TM6, which plays a crucial role in drug binding in human cells. For the first time, we employ molecular dynamics (MD) simulations based on the self-consistent charge density functional tight binding (SCC-DFTB) method to elucidate the potential sites of fragmentation and mutation in this domain upon impact of OH radicals, and to obtain fundamental information about the underlying reaction mechanisms. Furthermore, we apply non-reactive MD simulations to investigate the long-term effect of this mutation, with possible implications for drug binding. Our simulations indicate that the interaction of OH radicals with TM6 might lead to the breaking of C-C and C-N peptide bonds, which eventually cause fragmentation of TM6. Moreover, according to our simulations, the OH radicals can yield mutation in the aromatic ring of phenylalanine in TM6, which in turn affects its structure. As TM6 plays an important role in the binding of a range of cytotoxic drugs with P-glycoprotein, any changes in its structure are likely to affect the response of the tumor cell in chemotherapy. This is crucial for cancer therapies based on reactive oxygen species, such as plasma treatment.

Abstract: An extensive CO2 plasmamodel is presented that is relevant for the production of ‘‘solar fuels.’’ It is based on reaction rate coefficients fromrigorously reviewed literature, and is augmented with reactionrate coefficients that are obtained fromscaling laws.The input data set,which is suitable for usage with the plasma simulation software Plasimo (https://plasimo.phys.tue.nl/), is available via the Plasimo and publisher’s websites.1 The correctness of this model implementation has been established by independent ZDPlasKin implementation (http://www.zdplaskin.
laplace.univ-tlse.fr/), to verify that the results agree. Results of these ‘‘global models’’ are presented for a DBD plasma reactor.

Abstract: Molecular dynamics simulations were used to investigate the size-dependent melting mechanism of nickel nanoclusters of various sizes. The melting process was monitored by the caloric curve, the overall cluster Lindemann index, and the atomic Lindemann index. Size-dependent melting temperatures were determined, and the correct linear dependence on inverse diameter was recovered. We found that the melting mechanism gradually changes from dynamic coexistence melting to surface melting with increasing cluster size. These findings are of importance in better understanding carbon nanotube growth by catalytic chemical vapor deposition as the phase state of the catalyst nanoparticle codetermines the growth mechanism.

Abstract: In dc magnetrons the electrons emitted from the cathode may return there due to the applied magnetic field. When that happens, they can be recaptured or reflected back into the discharge, depending on the value of the reflection coefficient (RC). A 2d3v (two-dimensional in coordinate and three-dimensional in velocity space) particle-in-cellMonte Carlo model, including an external circuit, is developed to determine the role of the electron recapture in the discharge processes. The detailed discharge structure as a function of RC for two pressures (4 and 25mtorr) is studied. The importance of electron recapture is clearly manifested, especially at low pressures. The results indicate that the discharge characteristics are dramatically changed with varying RC between 0 and 1. Thus, the electron recapture at the cathode appears to be a significant mechanism in magnetron discharges and RC a very important parameter in their correct quantitative description that should be dealt with cautiously.

Abstract: In this work we have analyzed the properties of a gliding dc discharge in argon at atmospheric pressure. Despite the usual designation of these discharges as ‘gliding arc discharges’, it was found previously that they operate in two different regimes—glow and arc. Here we analyze the differences in both regimes by means of two dimensional fluid modeling. In order to address different aspects of the discharge operation, we use two models—Cartesian and axisymmetric in a cylindrical coordinate system. The obtained results show that the two types of discharges produce a similar plasma column for a similar discharge current. However, the different mechanisms of plasma channel attachment to the cathode could produce certain differences in the plasma parameters (i.e. arc elongation), and this can affect gas treatments applications.

Abstract: CO2 dissociation and its subsequent conversion into added-value chemicals is a promising strategy for recycling CO2 gas into reusable products. One of the possible methods is direct plasma-induced dissociation. In this work we study the efficiency of CO2 dissociation in pulsed atmospheric-pressure gas discharge between two conducting electrodes by a 0-D numerical plasma model. The purpose of the study is to provide results on the optimal conditions of CO2 conversion with respect to the energy efficiency and dissociation by varying the maximum power density value and the pulse length. The power density is directly related to the discharge current and the reduced electric field in the discharge. We consider pulse lengths in the range from hundreds of nanosecond up to milliseconds. The results obtained show that the dissociation degree and energy efficiency are sensitive to the pulse length (duration) and the power density, so that a considerable improvement of the discharge performance can be achieved by fine-tuning these parameters. The study is intended to provide guidance in designing an experimental set-up and a power supply with the characteristics necessary to achieve optimal conversion.

Abstract: We use a zero-dimensional reaction kinetics model to simulate CO2 conversion in microwave discharges where the excitation of the vibrational levels plays a significant role in the dissociation kinetics. The model includes a description of the CO2 vibrational kinetics, taking into account state-specific VT and VV relaxation reactions and the effect of vibrational excitation on other chemical reactions. The model is used to simulate a general tubular microwave reactor, where a stream of CO2 flows through a plasma column generated by microwave radiation. We study the effects of the internal plasma parameters, namely the reduced electric field, electron density and the total specific energy input, on the CO2 conversion and its energy efficiency. We report the highest energy efficiency (up to 30%) for a specific energy input in the range 0.41.0 eV/molecule and a reduced electric field in the range 50100 Td and for high values of the electron density (an ionization degree greater than 10−5). The energy efficiency is mainly limited by the VT relaxation which contributes dominantly to the vibrational energy losses and also contributes significantly to the heating of the reacting gas. The model analysis provides useful insight into the potential and limitations of CO2 conversion in microwave discharges.

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.