Abstract: Recently, core shell silicon nanowires (Si-NWs) have been envisaged to be used for field-effect transistors and photovoltaic applications. In spite of the constant downsizing of such devices, the formation of ultrasmall diameter core shell Si-NWs currently remains entirely unexplored. We report here on the modeling of the formation of such core shell Si-NWs using a dry thermal oxidation of 2 nm diameter (100) Si nanowires at 300 and 1273 K, by means of reactive molecular dynamics simulations using the ReaxFF potential. Two different oxidation mechanisms are discussed, namely a self-limiting process that occurs at low temperature (300 K), resulting in a Si core I ultrathin SiO2 silica shell nanowire, and a complete oxidation process that takes place at a higher temperature (1273 K), resulting in the formation of an ultrathin SiO2 silica nanowire. The oxidation kinetics of both cases and the resulting structures are analyzed in detail. Our results demonstrate that precise control over the Si-core radius of such NWs and the SiOx (x <= 2.0) oxide shell is possible by controlling the growth temperature used during the oxidation process.

Abstract: Atomic scale simulations have been shown to be a powerful tool for elucidating the growth mechanisms of carbon nanotubes. The growth picture is however not entirely clear yet due to the gap between current simulations and real experiments. We here simulate for the first time the nucleation and subsequent growth of single-wall carbon nanotubes (SWNTs) from oxygen-containing hydrocarbon feedstocks using the hybrid Molecular Dynamics/Monte Carlo technique. The underlying nucleation mechanisms of Ni-catalysed SWNT growth are discussed in detail. Specifically, we find that as a function of the feedstock, different carbon fractions may emerge as the main growth species, due to a competition between the feedstock decomposition, its rehydroxylation and its contribution to etching of the growing SWNT. This study provides a further understanding of the feedstock effects in SWNT growth in comparison with available experimental evidence as well as with<italic>ab initio</italic>and other simulation data, thereby reducing the simulation–experiment gap.

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: Magnetic Pole Enhanced-Inductively Coupled Plasmas (MaPE-ICPs) in analogy to the conventional ICPs exhibit two modes of operation, depending on the power coupling mechanism, i.e., a low power mode with dominant capacitive coupling (E-mode) and a high power mode with dominant inductive coupling (H-mode). A comparative study of the electron temperature measured by a Langmuir probe (T-e(LP)) and the electron excitation temperature (T-exc(OES)) determined by Optical Emission Spectroscopy (OES) is reported in the two distinct modes of a MaPE-ICP operated in argon. The dependence of T-e(LP), T-exc(OES) and their ratio (T-e(LP)/T-exc(OES)) on applied power (5-50 W) and gas pressure (15-60 mTorr) is explored, and the validity of T-exc(OES) as an alternative diagnostic to T-e(LP) is tested in the two modes of MaPE-ICP. The OES based non-invasive measurement of the plasma parameters such as electron temperature is very useful for plasma processing applications in which probe measurements are limited. (C) 2013 Elsevier B. V. All rights reserved.

Abstract: The application of nonthermal atmospheric pressure plasma is emerging as an alternative and efficient technique for the inactivation of bacterial biofilms. In this study, reactive molecular dynamics simulations were used to examine the reaction mechanisms of hydroxyl radicals, as key reactive oxygen plasma species in biological systems, with several organic molecules (i.e., alkane, alcohol, carboxylic acid, and amine), as prototypical components of biomolecules in the biofilm. Our results demonstrate that organic molecules containing hydroxyl and carboxyl groups may act as trapping agents for the OH radicals. Moreover, the impact of OH radicals on N-acetyl-glucosamine, as constituent component of staphylococcus epidermidis biofilms, was investigated. The results show how impacts of OH radicals lead to hydrogen abstraction and subsequent molecular damage. This study thus provides new data on the reaction mechanisms of plasma species, and particularly the OH radicals, with fundamental components of bacterial biofilms.

Abstract: Cold atmospheric pressure plasmas are gaining increased interest in the medical sector and clinical trials to treat skin diseases are underway. Plasmas are capable of producing several reactive oxygen and nitrogen species (RONS). However, there are open questions how plasma-generated RONS interact on a molecular level in a biological environment, e.g. cells or cell components. The redox pair glutathione (GSH) and glutathione disulphide (GSSG) forms the most important redox buffer in organisms responsible for detoxification of intracellular reactive species. We apply Raman spectroscopy, mass spectrometry, and molecular dynamics simulations to identify the time-dependent chemical modifications on GSH and GSSG that are caused by dielectric barrier discharge under ambient conditions. We find GSSG, S-oxidised glutathione species, and S-nitrosoglutathione as oxidation products with the latter two being the final products, while glutathione sulphenic acid, glutathione sulphinic acid, and GSSG are rather reaction intermediates. Experiments using stabilized pH conditions revealed the same main oxidation products as were found in unbuffered solution, indicating that the dominant oxidative or nitrosative reactions are not influenced by acidic pH. For more complex systems these results indicate that too long treatment times can cause difficult-to-handle modifications to the cellular redox buffer which can impair proper cellular function.

Abstract: The effect of gas heating in laboratory sputter magnetrons is investigated by means of numerical modeling. The model is two-dimensional in the coordinate space and three-dimensional in the velocity space based on the particle-in-cellMonte Carlo collisions technique. It is expanded in a way that allows the inclusion of the neutral plasma particles (fast gas atoms and sputtered atoms), which makes it possible to calculate the gas temperature and its influence on the discharge behavior in a completely self-consistent way. The results of the model are compared to experimental measurements and to other existing simulation results. The results show that gas heating is pressure dependent (rising with the increase in the gas pressure) and should be taken into consideration at pressures above 10 mTorr.

Abstract: In this study we report on a 2D fluid model of a gliding arc discharge in argon. Despite the 3D nature of the discharge, 2D models are found to be capable of providing very useful information about the operation of the discharge. We employ two modelsan axisymmetric and a Cartesian one. We show that for the considered experiment and the conditions of a low current arc (around 30 mA) in argon, there is no significant heating of the cathode surface and the discharge is sustained by field electron emission from the cathode accompanied by the formation of a cathode spot. The obtained discharge power and voltage are relatively sensitive to the surface properties and particularly to the surface roughness, causing effectively an amplification of the normal electric field. The arc body and anode region are not influenced by this and depend mainly on the current value. The gliding of the arc is modelled by means of a 2D Cartesian model. The arcelectrode contact points are analysed and the gliding mechanism along the electrode surface is discussed. Following experimental observations, the cathode spot is simulated as jumping from one point to another. A complete arc cycle is modelled from initial ignition to arc decay. The results show that there is no interaction between the successive gliding arcs.

Abstract: In this work we study energy transport in a gliding arc discharge with two diverging flat
electrodes in argon gas at atmospheric pressure. The discharge is ignited at the shortest electrode
gap and it is pushed downstream by a forced gas flow. The current values considered are
relatively low and therefore a non-equilibrium plasma is produced. We consider two cases, i.e.
with high and low discharge current (28 mA and 2.8mA), and a constant gas flow of 10 lmin −1 ,
with a significant turbulent component to the velocity. The study presents an analysis of the
various energy transport mechanisms responsible for the redistribution of Joule heating to the
plasma species and the moving background gas. The objective of this work is to provide a
general understanding of the role of the different energy transport mechanisms in arc formation
and sustainment, which can be used to improve existing or new discharge designs. The work is
based on a three-dimensional numerical model, combining a fluid plasma model, the shear stress
transport Reynolds averaged Navier–Stokes turbulent gas flow model, and a model for gas
thermal balance. The obtained results show that at higher current the discharge is constricted
within a thin plasma column several hundred kelvin above room temperature, while in the low-
current discharge the combination of intense convective cooling and low Joule heating prevents
discharge contraction and the plasma column evolves to a static non-moving diffusive plasma,
continuously cooled by the flowing gas. As a result, the energy transport in the two cases is
determined by different mechanisms. At higher current and a constricted plasma column, the
plasma column is cooled mainly by turbulent transport, while at low current and an unconstricted
plasma, the major cooling mechanism is energy transport due to non-turbulent gas convection. In
general, the study also demonstrates the importance of turbulent energy transport in
redistributing the Joule heating in the arc and its significant role in arc cooling and the formation
of the gas temperature profile. In general, the turbulent energy transport lowers the average gas
temperature in the arc, thus allowing additional control of thermal non-equilibrium in the
discharge.

Abstract: The modelling of a gliding arc discharge (GAD) is studied by means of the quasineutral (QN) plasma modelling approach. The model is first evaluated for reliability and proper description of a gliding arc discharge at atmospheric pressure, by comparing with a more elaborate non-quasineutral (NQN) plasma model in two different geometries – a 2D axisymmetric and a Cartesian geometry. The NQN model is considered as a reference, since it provides a continuous self-consistent plasma description, including the near electrode regions. In general, the results of the QN model agree very well with those obtained from the NQN model. The small differences between both models are attributed to the approximations in the derivation of the QN model. The use of the QN model provides a substantial reduction of the computation time compared to the NQN model, which is crucial for the development of more complex models in three dimensions or with complicated chemistries. The latter is illustrated for (i) a reverse vortex flow(RVF) GAD in argon, and (ii) a GAD in CO2. The RVF discharge is modelled in three dimensions and the effect of the turbulent heat transport on the plasma and gas characteristics is
discussed. The GAD model in CO2 is in a 1D geometry with axial symmetry and provides results for the time evolution of the electron, gas and vibrational temperature of CO2, as well as for the molar fractions of the different species.

Abstract: The lowest-energy configurations of Cn(n≤55) clusters are obtained using the energy-minimization technique with the conjugate gradient method where a modified Brenner potential is invoked to describe the carbon and hydrocarbon interaction. We found that the ground-state configuration consists of a single ring for small number of C atoms and multiring structures are found with increasing n, which can be in planar, bowl-like or caplike form. Contrary to previous predictions, the binding energy Eb does not show even-odd oscillations and only small jumps are found in the Eb(n) curve as a consequence of specific types of edges or equivalently the number of secondary atoms. We found that hydrogenation of the edge atoms may change the ground-state configuration of the nanocluster. In both cases we determined the magic clusters. Special attention is paid to trigonal and hexagonal shaped carbon clusters and to clusters having a graphenelike configuration. Trigonal clusters are never the ground state while hexagonal-shaped clusters are only the ground state when they have zigzag edges.

Abstract: Natural deep eutectic solvents (NADES) show great promise as media for enzymatic reactions in areas where (bio)compatibility with natural or medicinal products is a must. While in theory they can be tailored to the intended reaction to ensure optimized yields, the knowledge to date is predominantly empirical, with some mechanistic reports providing a fragmented view at best. Therefore, it is not easy to explain experimental observations, let alone make predictions. The aim of this study was to develop a structured, holistic understanding of the effects of NADES media on enzymatic reactions, distinguishing between effects on solubility, solvation, viscosity, inhibition and denaturation. Experimental and computational chemistry methods were combined to separately study the interactions between enzyme, substrate, and NADES as reaction media. The initial enzyme activity and the final conversion of vinyl laurate transesterification by immobilized Candida antarctica lipase were studied experimentally. The direct effect of NADES on the same enzyme was modeled by molecular dynamics simulation. The effect of solubility was studied by both experimental and computational methods. To predict the solubility and viscosity of NADES, data-driven models were developed by combining group contribution and machine learning methods, based on the accumulated experimental knowledge on NADES found in the literature. Finally, the composed relationships and prediction models were applied to the practical example of deacetylation of mannosylerythritol lipids (MELs). The experimental findings show that the chosen NADES system has a significant effect on both the apparent initial activity and the final conversion. However, in the simulations, the enzyme retains its original structure; moreover, NADES has an additional stabilizing effect on the enzyme. In addition, changes in the molar ratio of the compounds in NADES do not show a significant effect on the stability of the enzyme. These results indicate that the main effect of NADES on the reaction is mainly related to the substrate-solvent interactions (solvation energy) and the viscosity of the system. On the other hand, the experimental results only confirmed the significance of solvation, viscosity did not show a clear correlation with the studied reaction parameters. The machine learning models built for solubility and viscosity gave quantitative predictions of these properties. The accumulated knowledge was used to optimize the yield in the deacetylation reaction of MELs. The combination of these methods provides fundamental knowledge about the effect of NADES on biocatalysis, but the results are also applicable to other uses of NADES.

Abstract: Purpose Natural deep eutectic solvents (NADES) represent a green alternative to conventional organic solvents as reaction medium, offering more benign properties. To efficiently design NADES for biocatalysis, a better understanding of their effect on these reactions is needed. We hypothesize that this effect can be described by separately considering (1) the solvent interactions with the substrates, (2) the solvent viscosities and (3) the enzyme stability in NADES. Methods We investigated the effect of substrate solvation and viscosity on the reaction rate; and the stability of the enzyme in NADES. To this end, we monitored the conversion over time of the transesterification of vinyl laurate with 1- butanol by the lipase enzyme Candida antarctica B in NADES of different compounds and molar ratios. Results The initial reaction rate is higher in most NADES ( varying between 1.14 and 15.07 mu mol min(-1) mg(-1)) than in the reference n-hexane (4.0 mu mol min(-1) mg(-1))), but no clear relationship between viscosity and initial reaction rate was found. The increased reaction rate is most likely related to the solvation of the substrate due to a change in the activation energy of the reaction or a change in the conformation of the substrate. The enzyme retained part of its activity after the first 2 h of reaction (on average 20 % of the substrate reacted in the 2-24 h period). Enzyme incubation in ethylene glycol-based NADES resulted in a reduced reaction rate ( 15.07 vs. 3.34 mu mol min(-1) mg(-1)), but this may also be due to slow dissolution of the substrate. Conclusions The effect of viscosity seems to be marginal next to the effect of solvation and possible enzyme-NADES interaction. The enzyme retains some of its activity during the 24-hour measurements, but the enzyme incubation experiments did not yield accurate, comparable values. [GRAPHICS] .