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: In this paper, we present the results of combined molecular dynamics−Metropolis Monte Carlo (MD-MMC) simulations of hydrocarbon species at flat diamond (100)2 × 1 and (111)1 × 1 surfaces. The investigated species are considered to be the most important growth species for (ultra)nanocrystalline diamond ((U)NCD) growth. When applying the MMC algorithm to stuck species at monoradical sites, bonding changes are only seen for CH2. The sequence of the bond breaking and formation as put forward by the MMC simulations mimics the insertion of CH2 into a surface dimer as proposed in the standard growth model of diamond. For hydrocarbon species attached to two adjacent radical (biradical) sites, the MMC simulations give rise to significant changes in the bonding structure. For UNCD, the combinations of C3 and C3H2, and C3 and C4H2 (at diamond (100)2 × 1) and C and C2H2 (at diamond (111)1 × 1) are the most successful in nucleating new crystal layers. For NCD, the following combinations pursue the diamond structure the best: C2H2 and C3H2 (at diamond (100)2 × 1) and CH2 and C2H2 (at diamond (111)1 × 1). The different behaviors of the hydrocarbon species at the two diamond surfaces are related to the different sterical hindrances at the diamond surfaces.

Abstract: Plasma medicine has been claimed to provide a novel route to heal wounds and regenerate skin, although very little is currently known about the elementary processes taking place. We carried out a series of ReaxFF-based reactive molecular dynamics simulations to investigate the interaction of O and OH radicals with lipids, more specifically with α-linolenic acid as a model for the free fatty acids present in the upper skin layer. Our calculations predict that the O and OH radicals most typically abstract a H atom from the fatty acids, which can lead to the formation of a conjugated double bond, but also to the incorporation of alcohol or aldehyde groups, thereby increasing the hydrophilic character of the fatty acids and changing the general lipid composition of the skin. Within the limitations of the investigated model, no formation of possibly toxic products was observed.

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: In this paper, a new implementation of the Metropolis Monte Carlo (MMC) algorithm is presented. When combining the MMC model with a molecular dynamics (MD) code, crystal growth by plasma-enhanced chemical vapor deposition can be simulated. As the MD part simulates impacts of growth species onto the surface on a time scale of picoseconds, the MMC algorithm simulates the slower adatom surface processes. The implementation includes a criterion for the selection of atoms that are allowed to be displaced during the simulation, and a criterion of after how many MMC cycles the simulation is stopped. We performed combined MD-MMC simulations for hydrocarbon species that are important for the growth of ultrananocrystalline diamond (UNCD) films at partially hydrogenated diamond surfaces, since this implementation is part of a study of the growth mechanisms of (ultra)nanocrystalline diamond films. Exemplary for adatom arrangements during the growth of UNCD, the adatom surface behavior of C and C2H2 at diamond (111)1 × 1, C and C4H2 at diamond (111)1 × 1 and C3 at diamond (100)2 × 1 has been investigated. For all cases, the diamond crystal structure is pursued under the influence of MMC simulation. Additional longer time-scale MD simulations put forward very similar structures, verifying the MMC algorithm. Nevertheless, the MMC simulation time is typically one order of magnitude shorter than the MD simulation time.

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: We present here some of our modeling efforts for PECVD growth of nanostructured carbon materials with focus on amorphous hydrogenated carbon. Experimental data from an expanding thermal plasma setup were used as input for the simulations. Attention was focused both on the film growth mechanism, as well as on the hydrocarbon reaction mechanisms during growth of the films. It is found that the reaction mechanisms and sticking coefficients are dependent on the specific surface sites, and the structural properties of the growth radicals. The film growth results are in correspondence with the experiment. Furthermore, it is found that thin a-C:H films can be densified using an additional H-flux towards the substrate.

Abstract: Magnetron sputtering is a widely used physical vapor deposition technique for deposition and formation of nanocatalyst thin films and clusters. Nevertheless, so far only few studies investigated this formation process at the fundamental level. We here review atomic scale molecular dynamics simulations aimed at elucidating the nanocatalyst growth process through magnetron sputtering. We first introduce the basic magnetron sputtering background and machinery of molecular dynamics simulations, and then describe the studies conducted in this field so far. We also present a perspective view on how the field may be developed further.

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: This paper outlines two different numerical simulation approaches, carried out by our group, used for describing hydrocarbon plasmas in their applications for either nanoparticle formation in the plasma or the growth of nanostructured thin films, such as nanocrystalline diamond (NCD). A plasma model based on the fluid approach is utilized to study the initial mechanisms giving rise to nanoparticle formation in an acetylene plasma. The growth of NCD is investigated by molecular dynamics simulations, describing the interaction of the hydrocarbon species with a substrate.