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: 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: A method has been developed for damage-free cesium (Cs) encapsulation within single-walled carbon nanotubes (SWNTs) with fine position selectivity. Precise energy tuning of Cs-ion irradiation revealed that there is a clear energy window (2060 eV) for the efficient encapsulation of Cs through the hexagonal network of SWNT sidewalls without causing significant damage. This minimum energy threshold of Cs-ion encapsulation (∼20 eV) matches well with the value obtained by ab initio simulation (∼22 eV). Furthermore, position-selective Cs encapsulation was carried out, resulting in the successful formation of pn-junction SWNT thin films with excellent environmental stability.

Abstract: Methane reforming by plasma catalysis is a complex process that is far from understood. It requires a multidisciplinary approach which ideally takes into account all effects from the plasma on the catalyst, and vice versa. In this contribution, we focus on the interactions of CHx (x = {1,2,3}) radicals that are created in the plasma with several nickel catalyst surfaces. To this end, we perform reactive molecular dynamics simulations, based on the ReaxFF potential, in a wide temperature range of 4001600 K. First, we focus on the H2 formation as a function of temperature and surface structure. We observe that substantial H2 formation is obtained at 1400 K and above, while the role of the surface structure seems limited. Indeed, in the initial stage, the type of nickel surface influences the CH bond breaking efficiency of adsorbed radicals; however, the continuous carbon diffusion into the surface gradually diminishes the surface crystallinity and therefore reduces the effect of surface structure on the H2 formation probability. Furthermore, we have also investigated to what extent the species adsorbed on the catalyst surface can participate in surface reactions more in general, for the various surface structures and as a function of temperature. These results are part of the ongoing research on the methane reforming by plasma catalysis, a highly interesting yet complex alternative to conventional reforming processes.

Abstract: In this paper, we quantitatively investigate with atom probe tomography, the effect of temperature on the interfacial transition layer suboxide species due to the thermal oxidation of silicon. The chemistry at the interface was measured with atomic scale resolution, and the changes in chemistry and intermixing at the interface were identified on a nanometer scale. We find an increase of suboxide (SiOx) concentration relative to SiO2 and increased oxygen ingress with elevated temperatures. Our experimental findings are in agreement with reactive force field molecular dynamics simulations. This work demonstrates the direct comparison between atom probe derived chemical profiles and atomistic-scale simulations for transitional interfacial layer of suboxides as a function of temperature.

Abstract: Carbon nanotubes (CNTs) are nowadays routinely grown in a thermal CVD setup. State-of-the-art plasma-enhanced CVD (PECVD) growth, however, offers advantages over thermal CVD. A lower growth temperature and the growth of aligned freestanding single-walled CNTs (SWNTs) makes the technique very attractive. The atomic scale growth mechanisms of PECVD CNT growth, however, remain currently entirely unexplored. In this contribution, we employed molecular dynamics simulations to focus on the effect of applying an electric field on the SWNT growth process, as one of the effects coming into play in PECVD. Using sufficiently strong fields results in (a) alignment of the growing SWNTs, (b) a better ordering of the carbon network, and (c) a higher growth rate relative to thermal growth rate. We suggest that these effects are due to the small charge transfer occurring in the Ni/C system. These simulations constitute the first study of PECVD growth of SWNTs on the atomic level.

Abstract: This paper reports on the incorporation of three commercial fluorescent dyes, i.e., rhodamine 6G, fluorescein, and fluorescent brightener 184, in plasma coatings, by utilizing a dielectric barrier discharge (DBD) reactor, and the subsequent monitoring of the coatings homogeneity based on the emitted fluorescent light. The plasma coatings are qualitatively characterized with fluorescence microscopy, UVvis spectroscopy and profilometry for the determination of the coating thickness. The emitted fluorescent light of the coating correlates to the amount of dye per area, and deviations of these factors can hence be observed by monitoring the intensity of this light. This allows monitoring the homogeneity of the plasma coatings in a fast and simple way, without making major adjustments to the process.

Abstract: Using reactive molecular dynamics simulations, we have investigated the effect of single-impact, low-energy (thermal-100 eV) bombardment of a Si(100){2 × 1} surface by atomic and molecular oxygen. Penetration probability distributions, as well as defect formation distributions, are presented as a function of the impact energy for both species. It is found that at low impact energy, defects are created chemically due to the chemisorption process in the top layers of the surface, while at high impact energy, additional defects are created by a knock-on displacement of Si. These results are of particular importance for understanding device performances of silica-based microelectronic and photovoltaic devices.

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: The interaction of thermal and hyperthermal Ni ions with gas-phase C60 fullerene was investigated at two temperatures with classical molecular dynamics simulations using a recently developed interatomic many-body potential. The interaction between Ni and C60 is characterized in terms of the NiC60 binding sites, complex formation, and the collision and temperature induced deformation of the C60 cage structure. The simulations show how ion implantation theoretically allows the synthesis of both endohedral Ni@C60 and exohedral NiC60 metallofullerene complexes.

Abstract: Uniform acceptance force biased Monte Carlo (UFMC) simulations have previously been shown to be a powerful tool to simulate atomic scale processes, enabling one to follow the dynamical path during the simulation. In this contribution, we present a simple proof to demonstrate that this uniform acceptance still complies with the condition of detailed balance, on the condition that the characteristic parameter lambda = 1/2 and that the maximum allowed step size is chosen to be sufficiently small. Furthermore, the relation to Metropolis Monte Carlo (MMC) is also established, and it is shown that UFMC reduces to MMC by choosing the characteristic parameter lambda = 0 [Rao, M. et al. Mol. Phys. 1979, 37, 1773]. Finally, a simple example compares the UFMC and MMC methods.

Abstract: The behavior of hydrocarbon species at step edges of diamond terraces is investigated by means of combined molecular dynamics−Metropolis Monte Carlo simulations. The results show that the formation of ballas-like diamond films (like UNCD) and well-faceted diamond films (like NCD) can be related to the gas phase concentrations of CxHy in a new manner: Species that have high concentrations above the growing UNCD films suppress the extension of step edges through defect formation. The species that are present above the growing NCD film, however, enhance the extension of diamond terraces, which is believed to result in well-faceted diamond films. Furthermore, it is shown that, during UNCD growth, CxHy species with x ≥ 2 play an important role, in contrast to the currently adopted CVD diamond growth mechanism. Finally, the probabilities for the extension of the diamond (100) terrace are much higher than those for the diamond (111) terrace, which is in full agreement with the experimental observation that diamond (100) facets are more favored than diamond (111) facets during CVD diamond growth.