Abstract: Controlling and steering the growth of single walled carbon nanotubes is often believed to require controlling of the nucleation stage. Yet, little is known about the microscopic mechanisms governing the nucleation from hydrocarbon molecules. Specifically, we address here the dehydrogenation of hydrocarbon molecules and the formation of all-carbon graphitic islands on metallic nanoclusters from hydrocarbon molecules under conditions typical for carbon nanotube growth. Employing reactive molecular dynamics simulations, we demonstrate for the first time that the formation of a graphitic network occurs through the intermediate formation of vertically oriented, not fully dehydrogenated graphitic islands. Upon dehydrogenation of these vertical graphenes, the islands curve over the surface, thereby forming a carbon network covering the nanoparticle. The results indicate that controlling the extent of dehydrogenation offers an additional parameter to control the nucleation of carbon nanotubes.

Abstract: Ion irradiation of carbon nanotubes very often leads to defect formation. However, we have recently shown that Ar ion irradiation in a limited energy window of 1025 eV may enhance the initial cap nucleation process, when the carbon network is in contact with the metal nanocatalyst. Here, we employ reactive molecular dynamics simulations to demonstrate that ion irradiation in a higher energy window of 1035 eV may also heal network defects after the nucleation stage through a non-metal-mediated mechanism, when the carbon network is no longer in contact with the metal nanocatalyst. The results demonstrate the possibility of beneficially utilizing ions in e.g. plasma-enhanced chemical vapour deposition of carbon nanotubes.

Abstract: Using reactive molecular dynamics simulations, the melting behavior of nickelcarbon nanoclusters is examined. The phase diagrams of icosahedral and Wulff polyhedron clusters are determined using both the Lindemann index and the potential energy. Formulae are derived for calculating the equilibrium constants and the solid and liquid fractions during a phase transition, allowing more rational determination of the melting temperature with respect to the arbitrary Lindemann value. These results give more insight into the properties of nickelcarbon nanoclusters in general and can specifically be very useful for a better understanding of the synthesis of carbon nanotubes using the catalytic chemical vapor deposition method.

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 reactive ReaxFF force field is developed for studying the influence of intrinsic point defects on the chemistry with TiO2 condensed phases. The force field parameters are optimized to ab initio data for the equations of state, relative phase stabilities for titanium and titanium dioxide, potential energy differences for (TiO2)n-clusters (n = 116). Also data for intrinsic point defects in anatase were added. These data contain formation energies for interstitial titanium and oxygen vacancies, diffusion barriers of the oxygen vacancies and molecular oxygen adsorption on a reduced anatase (101) surface. Employing the resulting force field, we study the influence of concentration of oxygen vacancies and expansion or compression of an anatase surface on the diffusion of the oxygen vacancies. Also the barrier for oxygen diffusion in the subsurface region is evaluated using this force field. This diffusion barrier of 27.7 kcal/mol indicates that the lateral redistribution of oxygen vacancies on the surface and in the subsurface will be dominated by their diffusion in the subsurface, since both this barrier as well as the barriers for diffusion from the surface to the subsurface and vice versa (17.07 kcal/mol and 21.91 kcal/mol, respectively, as calculated with DFT), are significantly lower than for diffusion on the surface (61.12 kcal/mol as calculated with DFT).

Abstract: Two-dimensional group-IV lattices silicene and germanene are known to share many of graphene's remarkable mechanical and electronic properties. Due to the out-of-plane buckling of the former materials, there are more means of electronic funtionalization, e.g. by applying uniaxial strain or an out-of-plane electric field. We consider monolayer hexagonal Sn (stanene) as an ideal candidate to feasibly implement and exploit graphene physics for nanoelectronic applications: with increased out-of-plane buckling and sizable spin-orbit coupling it lends itself to improved Dirac cone engineering. We investigate the ballistic charge transport regime of armchair Sn nanoribbons, classified according to the ribbon width W = {3m – 1, 3m, 3m + 1} with integer m. We study transport through (non-magnetic) armchair ribbons using a combination of density functional theory and non-equilibrium Green's functions. Sn ribbons have earlier current onsets and carry currents 20% larger than C/Si/Ge-nanoribbons as the contact resistance of these ribbons is found to be comparable. ((c) 2014 WILEY-VCH Verlag GmbH &Co. KGaA, Weinheim)

Abstract: We present a zero-dimensional kinetic model of CO2 splitting in non-equilibrium plasmas. The model includes a description of the CO2 vibrational kinetics (25 vibrational levels up to the dissociation limit of the molecule), taking into account state-specific VT and VV relaxation reactions and the effect of vibrational excitation on other chemical reactions. The model is applied to study the reaction kinetics of CO2 splitting in an atmospheric-pressure dielectric barrier discharge (DBD) and in a moderate-pressure microwave discharge. The model results are in qualitative agreement with published experimental works. We show that the CO2 conversion and its energy efficiency are very different in these two types of discharges, which reflects the important dissociation mechanisms involved. In the microwave discharge, excitation of the vibrational levels promotes efficient dissociation when the specific energy input is higher than a critical value (2.0 eV/molecule under the conditions examined). The calculated energy efficiency of the process has a maximum of 23%. In the DBD, vibrationally excited levels do not contribute significantly to the dissociation of CO2 and the calculated energy efficiency of the process is much lower (5%).

Abstract: In metallized film capacitors, there exists an air gap of about 0.2 μm between the films, with a pressure ranging generally from 130 atm. Because of the created potential difference between the two films, a microdischarge is formed in this gap. In this paper, we use an implicit particle-in-cell Monte Carlo collision simulation method to study the discharge properties in this direct-current microdischarge with 0.2 μm gap in a range of different voltages and pressures. The discharge process is significantly different from a conventional high pressure discharge. Indeed, the high electric field due to the small gap sustains the discharge by field emission. At low applied voltage (~15 V), only the electrons are generated by field emission, while both electrons and ions are generated as a stable glow discharge at medium applied voltage (~50 V). At still higher applied voltage (~100 V), the number of electrons and ions rapidly multiplies, the electric field reverses, and the discharge changes from a glow to an arc regime.

Abstract: Cyrogenic etching of silicon is envisaged to enable better control over plasma processing in the microelectronics industry, albeit little is known about the fundamental differences compared to the room temperature process. We here present molecular dynamics simulations carried out to obtain sticking probabilities, thermal desorption rates, surface diffusion speeds, and sputter yields of F, F2, Si, SiF, SiF2, SiF3, SiF4, and the corresponding ions on Si(100) and on SiF13 surfaces, both at cryogenic and near room temperature. The different surface behavior during conventional etching and cryoetching is discussed. F2 is found to be relatively reactive compared to other species like SiF03. Thermal desorption occurs at a significantly lower rate under cryogenic conditions, which results in an accumulation of physisorbed species. Moreover, ion incorporation is often observed for ions with energies of 30400 eV, which results in a relatively low net sputter yield. The obtained results suggest that the actual etching of Si, under both cryogenic and near room temperature conditions, is based on the complete conversion of the Si surface to physisorbed SiF4, followed by subsequent sputtering of these molecules, instead of direct sputtering of the SiF03 surface.

Abstract: A two-dimensional hybrid Monte Carlofluid model is applied to study HBr/He inductively coupled plasmas used for etching of Si. Complete sets of gas-phase and surface reactions are presented and the effects of the gas mixing ratio on the plasma characteristics and on the etch rates are discussed. A comparison with experimentally measured etch rates is made to validate the modelling results. The etch rate in the HBr plasma is found to be quite low under the investigated conditions compared to typical etch rates of Si with F- or Cl-containing gases. This allows for a higher control and fine-tuning of the etch rate when creating ultra-small features. Our calculations predict a higher electron temperature at higher He fraction, because the electrons do not lose their energy so efficiently in vibrational and rotational excitations. As a consequence, electron impact ionization and dissociation become more important, yielding higher densities of ions, electrons and H atoms. This results in more pronounced sputtering of the surface. Nevertheless, the overall etch rate decreases upon increasing He fraction, suggesting that chemical etching is still the determining factor for the overall etch rate.

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: Plasma technology is gaining increasing interest for the splitting of CO2 into CO and O2. We have performed experiments to study this process in a dielectric barrier discharge (DBD) plasma with a wide range of parameters. The frequency and dielectric material did not affect the CO2 conversion and energy efficiency, but the discharge gap can have a considerable effect. The specific energy input has the most important effect on the CO2 conversion and energy efficiency. We have also presented a plasma chemistry model for CO2 splitting, which shows reasonable agreement with the experimental conversion and energy efficiency. This model is used to elucidate the critical reactions that are mostly responsible for the CO2 conversion. Finally, we have compared our results with other CO2 splitting techniques and we identified the limitations as well as the benefits and future possibilities in terms of modifications of DBD plasmas for greenhouse gas conversion in general.

Abstract: Reactive molecular dynamics simulations are performed to study the interaction of reactive oxygen species, such as OH, HO2 and H2O2, with the endotoxic biomolecule lipid A of the gram-negative bacterium Escherichia coli. It is found that the aforementioned plasma species can destroy the lipid A, which consequently results in reducing its toxic activity. All bond dissociation events are initiated by hydrogen-abstraction reactions. However, the mechanisms behind these dissociations are dependent on the impinging plasma species, i.e. a clear difference is observed in the mechanisms upon impact of HO2 radicals and H2O2 molecules on one hand and OH radicals on the other hand. Our simulation results are in good agreement with experimental observations.

Abstract: The adsorption of C and CHx radicals on anatase (001) was studied using DFT within the generalized gradient approximation using the Perde-Burke-Ernzerhof (PBE) functional. We have studied the influence of oxygen vacancies in and at the surface on the adsorption properties of the radicals. For the oxygen vacancies in anatase (001), the most stable vacancy is located at the surface. For this vacancy, the maximal adsorption strength of C and CH decreases compared to the adsorption on the stoichiometric surface, but it increases for CH2 and CH3. If an oxygen vacancy is present in the first subsurface layer, the maximal adsorption strength increases for C, CH, CH2, and CH3. When the vacancy is present in the next subsurface layer, we find that only the CH3 adsorption is enhanced, while the maximal adsorption energies for the other radical species decrease. Not only does the precise location of the oxygen vacancy determine the maximal adsorption interaction, it also influences the adsorption strengths of the radicals at different surface configurations. This determines the probability of finding a certain adsorption configuration at the surface, which in turn influences the possible surface reactions. We find that C preferentially adsorbs far away from the oxygen vacancy, while CH2 and CH3 adsorb preferentially at the oxygen vacancy site. A fraction of CH partially adsorbs at the oxygen vacancy, and another fraction adsorbs further away from the vacancy.

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.