Abstract: In this paper, an overview is given of modeling investigations carried out in our research group for a better understanding of plasmas used for medical, environmental and nano applications. The focus is both on modeling the plasma chemistry and the plasma-surface interactions. The plasma chemistry provides the densities and fluxes of the important plasma species. This information can be used as input when modeling the plasma-surface interactions. The combination of plasma simulations and plasma – surface interaction simulations provides a more comprehensive understanding of the underlying processes for these applications.

Abstract: While it is well-known that the plasma-enhanced catalytic chemical vapor deposition (PECVD) of carbon nanotubes (CNTs) offers a number of advantages over thermal CVD, the influence of the various individual contributing factors is not well understood. Especially the role of ions is unclear, since ions in plasmas are generally associated with sputtering rather than with growing a material. Even so, various studies have demonstrated the beneficial effects of ion bombardment during the growth of CNTs. This review looks at the role of the ions in plasma-enhanced CNT growth as deduced from both experimental and simulation studies. Specific attention is paid to the beneficial effects of ion bombardment. Based on the available literature, it can be concluded that ions can be either beneficial or detrimental for carbon nanotube growth, depending on the exact conditions and the control over the growth process.

Abstract: In this study, we performed density functional theory (DFT) calculations to investigate different reaction mechanisms of CO oxidation catalyzed by the Si atom embedded defective BC2N nanostructures as well as the analysis of the structural and electronic properties. The structures of all the complexes are optimized and characterized by frequency calculations at the M062X/6-31G* computational level. Also, The electronic structures and thermodynamic parameters of adsorbed CO and O-2 molecules over Si-doped BC2N nanostructures are examined in detail. Moreover, to investigate the curvature effect on the CO oxidation reaction, all the adsorption and CO oxidation reactions on a finite-sized armchair (6,6) Si-BC2NNT are also studied. Our results indicate that there can be two possible pathways for the CO oxidation with O-2 molecule: O-2(g) + CO(g) -> O-2(ads) + CO(ads) -> CO2(g) + O-(ads) and O-(ads) + CO(g) -> CO2(g). The first reaction proceeds via the Langmuir-Hinshelwood (LH) mechanism while the second goes through the Eley-Rideal (ER) mechanism. On the other hand, by increasing the tube diameter, the energy barrier increases due to the strong adsorption energy of the O-2 molecule which is related to its dissociation over the tube surface. Our calculations indicate that the two step energy barrier of the oxidation reaction over Si-BC2NNS is less than that over the Si-BC2NNT. Hence, Si-BC2NNS may serve as an efficient and highly activated substrate to CO oxidation rather than (4,4) Si-BC2NNT. (C) 2018 Elsevier B.V. All rights reserved.

Abstract: Doped Cr2O3 has been shown to be a p-type transparent conducting oxide (TCO). Its conductivity, however, is low. As for most p-type TCOs, the main problem is the high effective hole mass due to flat valence bands. We use first-principles methods to investigate whether one can increase the valence band dispersion (i.e. reduce the hole mass) by anion alloying with sulfur, while keeping the band gap large enough for transparency. The alloying concentrations considered are given by Cr(4)SxO(6-x), with x = 1-5. To be able to describe the electronic properties of these materials accurately, we first study Cr2O3, examining critically the accuracy of different density functionals and methods, including PBE, PBE+U, HSE06, as well as perturbative approaches within the GW approximation. Our results demonstrate that Cr4S2O4 has an optical band gap of 3.08 eV and an effective hole mass of 1.8 m(e). This suggests Cr4S2O4 as a new p-type TCO host candidate.

Abstract: Over the past decade, graphene oxide (GO) has emerged as a promising membrane material with superior separation performance and intriguing mechanical/chemical stability. However, its practical implementation remains very challenging primarily because of its undesirable swelling in an aqueous environment. Here, we demonstrated that dissociation of water molecules into H3O+ and OH- ions inside the interlayer gallery of a layered GO membrane can strongly affect its stability and performance. We reveal that H3O+ and OH- ions form clusters inside the GO laminates that impede the permeance of water and salt ions through the membrane. Dynamics of those clusters is sensitive to an external ac electric field, which can be used to tailor the membrane performance. The presence of H3O+ and OH- ions also leads to increased stability of the hydrogen bond (H-bond) network among the water molecules and the GO layers, which further reduces water permeance through the membrane, while crucially imparting stability to the layered GO membrane against undesirable swelling. KEYWORDS: layered graphene-oxide membrane, aqueous stability, H3O+ and OH- ions, external electric field, molecular dynamics

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: 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: Cold atmospheric pressure plasmas in and in contact with liquids represent a growing field of research for various applications. Understanding the interactions between the plasma generated species and the liquid is crucial. In this work we perform molecular dynamics (MD) simulations based on a quantum mechanical method, i.e. density-functional based tight-binding (DFTB), to examine the interactions of OH radicals and O atoms in bulk water. Our calculations reveal that the transport of OH radicals through water is not only governed by diffusion, but also by an equilibrium reaction of H-abstraction with water molecules. Furthermore, when two OH radicals encounter each other, they either form a stable cluster, or react, resulting in the formation of a new water molecule and an O atom. In addition,
the O atoms form either oxywater (when in singlet configuration) or they remain stable in solution (when in triplet configuration), stressing the important role that O atoms can play in aqueous solution, and in contact with biomolecules. Our observations are in line with both experimental and ab initio results from the literature.

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: In this paper, an overview is given of modeling activities going on in our research group, for describing the plasma chemistry and plasmasurface interactions in reactive plasmas. The plasma chemistry is calculated by a fluid approach or by hybrid Monte Carlo (MC)fluid modeling. An example of both is illustrated in the first part of the paper. The example of fluid modeling is given for a dielectric barrier discharge (DBD) in CH4/O2, to describe the partial oxidation of CH4 into value-added chemicals. The example of hybrid MCfluid modeling concerns an inductively coupled plasma (ICP) etch reactor in Ar/Cl2/O2, including also the description of the etch process. The second part of the paper deals with the treatment of plasmasurface interactions on the atomic level, with molecular dynamics (MD) simulations or a combination of MD and MC simulations.

Abstract: We present a method for area selective deposition of 2D WS2 nanoribbons with tunable thickness on a dielectric substrate. The process is based on a complete conversion of a prepatterned, H-terminated Si layer to metallic W by WF6, followed by in situ sulfidation by H2S. The reaction process, performed at 450 degrees C, yields nanoribbons with lateral dimension down to 20 nm and with random basal plane orientation. The thickness of the nanoribbons is accurately controlled by the thickness of the pre-deposited Si layer. Upon rapid thermal annealing at 900 degrees C under inert gas, the WS2 basal planes align parallel to the substrate.

Abstract: Despite significant progress in single walled carbon nanotube (SWCNT) production by plasma enhanced chemical vapor deposition (PECVD), the growth mechanism in this method is not clearly understood. We employ reactive molecular dynamics simulations to investigate how plasma-based deposition allows growth at low temperature. We first investigate the SWCNT growth mechanism at low and high temperatures under conditions similar to thermal CVD and PECVD. We then show how ion bombardment during the nucleation stage increases the carbon solubility in the catalyst at low temperature. Finally, we demonstrate how moderate energy ions sputter amorphous carbon allowing for SWCNT growth at 500 K. (C) 2013 Elsevier B. V. All rights reserved.