Abstract: Plasma-surface interactions are in general highly complex due to the interplay of many concurrent processes. Molecular dynamics simulations provide insight in some of these processes, subject to the accessible time and length scales, and the availability of suitable force fields. In this introductory tutorial-style review, we aim to describe the current capabilities and limitations of molecular dynamics simulations in this field, restricting ourselves to low-temperature nonthermal plasmas. Attention is paid to the simulation of the various fundamental processes occurring, including sputtering, etching, implantation, and deposition, as well as to what extent the basic plasma components can be accounted for, including ground state and excited species, electric fields, ions, photons, and electrons. A number of examples is provided, giving an bird’s eye overview of the current state of the field.

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: 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: Continued miniaturization of metal-oxide-semiconductor field-effect transistors (MOSFETs) requires an ever-decreasing thickness of the gate oxide. The structure of ultrathin silicon oxide films, however, critically depends on the oxidation mechanism. Using reactive atomistic simulations, we here demonstrate how the oxidation mechanism in hyperthermal oxidation of such structures may be controlled by the oxidation temperature and the oxidant energy. Specifically, we study the interaction of hyperthermal oxygen with energies of 15 eV with thin SiOx (x ≤ 2) films with a native oxide thickness of about 10 Å. We analyze the oxygen penetration depth probability and compare with results of the hyperthermal oxidation of a bare Si(100){2 × 1} (c-Si) surface. The temperature-dependent oxidation mechanisms are discussed in detail. Our results demonstrate that, at low (i.e., room) temperature, the penetrated oxygen mostly resides in the oxide region rather than at the SiOx|c-Si interface. However, at higher temperatures, starting at around 700 K, oxygen atoms are found to penetrate and to diffuse through the oxide layer followed by reaction at the c-Si boundary. We demonstrate that hyperthermal oxidation resembles thermal oxidation, which can be described by the DealGrove model at high temperatures. Furthermore, defect creation mechanisms that occur during the oxidation process are also analyzed. This study is useful for the fabrication of ultrathin silicon oxide gate oxides for metal-oxide-semiconductor devices as it links parameters that can be straightforwardly controlled in experiment (oxygen temperature, velocity) with the silicon oxide structure.

Abstract: The exact structure and properties of the Si vertical bar SiO2 interface are very important in microelectronics and photovoltaic devices such as metal-oxide-semiconductor field-effect transistors (MOSFETs) and solar cells. Whereas Si vertical bar SiO2 structures are traditionally produced by thermal oxidation, hyperthermal oxidation shows a number of promising advantages. However, the Si vertical bar SiO2 interface induced in hyperthermal Si oxidation has not been properly investigated yet. Therefore, in this work, the interface morphology and interfacial stresses during hyperthermal oxidation at room temperature are studied using reactive molecular dynamics simulations based on the ReaxFF potential. Interface thickness and roughness, as well as the bond length and bond angle distributions in the interface are discussed and compared with other models developed for the interfaces induced by traditional thermal oxidation. The formation of a compressive stress is observed. This compressive stress, which at the interface amounts about 2 GPa, significantly slows down the inward silica growth. This value is close to the experimental value in the Si vertical bar SiO2 interface obtained in traditional thermal oxidation.

Abstract: Uniform-acceptance force-bias Monte Carlo (fbMC) methods have been shown to be a powerful technique to access longer timescales in atomistic simulations allowing, for example, phase transitions and growth. Recently, a new fbMC method, the time-stamped force-bias Monte Carlo (tfMC) method, was derived with inclusion of an estimated effective timescale; this timescale, however, does not seem able to explain some of the successes the method. In this contribution, we therefore explicitly quantify the effective timescale tfMC is able to access for a variety of systems, namely a simple single-particle, one-dimensional model system, the Lennard-Jones liquid, an adatom on the Cu(100) surface, a silicon crystal with point defects and a highly defected graphene sheet, in order to gain new insights into the mechanisms by which tfMC operates. It is found that considerable boosts, up to three orders of magnitude compared to molecular dynamics, can be achieved for solid state systems by lowering of the apparent activation barrier of occurring processes, while not requiring any system-specific input or modifications of the method. We furthermore address the pitfalls of using the method as a replacement or complement of molecular dynamics simulations, its ability to explicitly describe correct dynamics and reaction mechanisms, and the association of timescales to MC simulations in general.

Abstract: Nanostructured carbon materials show a tremendous variety in atomic structure, morphology, properties, and applications. As all properties are ultimately determined by the structure of the material, a thorough understanding of the growth mechanisms that give rise to the particular structure is critical. On many occasions, it has been shown that plasma enhanced growth can be strongly beneficial. This review will describe the authors current understanding of plasma enhanced growth of carbon nanotubes, the prototypical example of nanostructured carbon materials, as obtained from experiments, simulations, and modeling. Specific emphasis is put on where experiments and computational approaches correspond, and where they differ. Also, the current status on simulating PECVD growth of some other carbon nanomaterials is reviewed, including amorphous carbon, graphene, and metallofullerenes. Finally, computational challenges with respect to the simulation of PECVD growth are identified.

Abstract: Plasma catalysis is an emerging new technology for the electrification and downscaling of NH3 synthesis. Increasing attention is being paid to the optimization of plasma catalysis with respect to the plasma conditions, the catalyst material, and their mutual interaction. In this work we use microkinetic models to study how the total conversion process is impacted by the combination of different plasma conditions and transition metal catalysts. We study how plasma-generated radicals and vibrationally excited N2 (present in a dielectric barrier discharge plasma) interact with the catalyst and impact the NH3 turnover frequencies (TOFs). Both filamentary and uniform plasmas are studied, based on plasma chemistry models that provided plasma phase speciation and vibrational distribution functions. The Langmuir−Hinshelwood reaction rate coefficients (i.e., adsorption reactions and subsequent reactions among adsorbates) are determined using conventional scaling relations. An additional set of Eley−Rideal reactions (i.e., direct reactions of plasma radicals with adsorbates) was added and a sensitivity analysis on the assumed reaction rate coefficients was performed. We first show the impact of different vibrational distribution functions on the catalytic dissociation of N2 and subsequent production of NH3, and we gradually include more radical reactions, to illustrate the contribution of these species and their corresponding reaction pathways. Analysis over a large range of catalysts indicates that different transition metals (metals such as Rh, Ni, Pt, and Pd) optimize the NH3TOFs depending on the population of the vibrational levels of N2. At higher concentrations of plasma-generated radicals, the NH3 TOFs become less dependent on the catalyst material, due to radical adsorptions on the more noble catalysts and Eley−Rideal reactions on the less noble catalysts.

Abstract: Low-temperature growth of carbon nanotubes (CNTs) has been claimed to provide a route towards chiral-selective growth, enabling a host of applications. In this contribution, we employ reactive molecular dynamics simulations to demonstrate how plasma-based deposition allows such low-temperature growth. We first show how ion bombardment during the growth affects the carbon dissolution and precipitation process. We then continue to demonstrate how a narrow ion energy window allows CNT growth at 500 K. Finally, we also show how CNTs in contrast cannot be grown in thermal CVD at this low temperature, but only at high temperature, in agreement with experimental data. (C) 2013 Elsevier Ltd. All rights reserved.

Abstract: The adsorption probability and reaction behavior of CHx plasma species on various nickel catalyst surfaces is investigated by means of reactive molecular dynamics (MD) simulations using the ReaxFF potential. Such catalysts are used in the reforming of hydrocarbons and in the growth of carbon nanotubes, and further insight in the underlying mechanisms of these processes is needed to increase their applicability. Single and consecutive impacts of CHx radicals (x={1,2,3}) were performed on four different Ni surfaces, at a temperature of 400 K. The adsorption probability is shown to be related to the number of free electrons, i.e. a higher number leads to more adsorptions, and the steric hindrance caused by the hydrogen atoms bonded to the impacting CHx species. Furthermore, some of the CH bonds break after adsorption, which generally leads to diffusion of the hydrogen atom over the surface. Additionally, these adsorbed H-atoms can be used in reactions to form new molecules, such as CH4 and C2Hx, although this is dependent on the precise morphology of the surface. New molecules are also formed by subtraction of H-atoms from adsorbed radicals, leading to occasional formation of H2 and C2Hx molecules.

Abstract: Industrial production of HCN from NH3 and CH4 not only uses precious Pt or Pt−Rh catalysts but also requires extremely high temperatures (∼1600 K). From an energetic, operational, and safety perspective, a drastic decrease in temperature is highly desirable. Here, we report ammonia reforming of methane for the production of HCN and H2 at 673 K by the combination of CH4/NH3 plasma and a supported Cu/silicalite-1 catalyst. 30% CH4 conversion has been achieved with 79% HCN selectivity. Catalyst characterization and plasma diagnostics reveal that the excellent reaction performance is attributed to metallic Cu active sites. In addition, we propose a possible reaction pathway, viz. E-R reactions with N, NH, NH2, and CH radicals produced in the plasma, for the production of HCN, based on density functional theory calculations. Importantly, the Cu/silicalite-1 catalyst costs less than 5% of the commercial Pt mesh catalyst.

Abstract: Nonthermal atmospheric pressure plasmas are gaining increasing attention for biomedical applications. However, very little fundamental information on the interaction mechanisms between the plasma species and biological cells is currently available. We investigate the interaction of important plasma species, such as OH, H2O2, O, O3, as well as O2 and H2O, with bacterial peptidoglycan by means of reactive molecular dynamics simulations, aiming for a better understanding of plasma disinfection. Our results show that OH, O, O3, and H2O2 can break structurally important bonds of peptidoglycan (i.e., CO, CN, or CC bonds), which consequently leads to the destruction of the bacterial cell wall. The mechanisms behind these breakups are, however, dependent on the impinging plasma species, and this also determines the effectiveness of the cell wall destruction.

Abstract: The combination of catalysis and nonthermal plasma holds promise for enabling difficult chemical conversions. The possible synergy between both depends strongly on the nature of the reactive plasma species and the catalyst material. In this paper, we show how vibrationally excited species and plasma-generated radicals interact with transition metal catalysts and how changing the catalyst material can improve the conversion rates and product selectivity. We developed a microkinetic model to investigate the impact of vibrational excitations and plasma-generated radicals on the nonoxidative coupling of methane over transition metal surfaces. We predict a significant increase in ethylene formation for vibrationally excited methane. Plasma-generated radicals have a stronger impact on the turnover frequencies with high selectivity toward ethylene on noble catalysts and mixed selectivity on non-noble catalysts. In general, we show how the optimal catalyst material depends on the desired products as well as the plasma conditions.