Abstract: Molecular dynamics simulations were used to investigate the size-dependent melting mechanism of nickel nanoclusters of various sizes. The melting process was monitored by the caloric curve, the overall cluster Lindemann index, and the atomic Lindemann index. Size-dependent melting temperatures were determined, and the correct linear dependence on inverse diameter was recovered. We found that the melting mechanism gradually changes from dynamic coexistence melting to surface melting with increasing cluster size. These findings are of importance in better understanding carbon nanotube growth by catalytic chemical vapor deposition as the phase state of the catalyst nanoparticle codetermines the growth mechanism.

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: In this paper, the splitting of CO2 in a pulsed plasma system, such as a dielectric barrier discharge (DBD), is evaluated from a chemical point of view by means of numerical modeling. For this purpose, a chemical reaction set of CO2 in an atmospheric pressure plasma is developed, including the vibrational states of CO2, O2, and CO. The simulated pulses are matched to the conditions of a filament (or microdischarge) and repeated with intervals of 1 μs. The influence of vibrationally excited CO2 as well as other neutral species, ions, and electrons on the CO2 splitting is discussed. Our calculations predict that the electrons have the largest contribution to the CO2 splitting at the conditions under study, by electron impact dissociation. The contribution of vibrationally excited CO2 levels in the splitting of CO2 is found be 6.4%, when only considering one microdischarge pulse and its afterglow, but it can be much higher for consecutive discharge pulses, as is typical for a filamentary DBD, when the interpulse time is short enough and accumulation effects in the vibrationally excited CO2 densities can occur.

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: The destruction of ethylene in a dielectric barrier discharge plasma is investigated by the combination of kinetic modeling and experiments, as a case study for plasma-based gas purification. The influence of the specific energy deposition on the removal efficiency and the selectivity toward CO and CO2 is studied for different concentrations of ethylene. The model allows the identication of the destruction pathway in dry and humid air. The latter is found to be mainly initiated by metastable N2 molecules, but the further destruction steps are dominated by O atoms and OH radicals. Upon increasing air humidity, the removal efficiency drops by ±15% (from 85% to 70%), but the selectivity toward CO and CO2 stays more or less constant at 60% and 22%, respectively. Beside CO and CO2, we also identified acetylene, formaldehyde, and water as byproducts of the destruction process, with concentrations of 1606 ppm, 15033 ppm, and 185 ppm in humid air (with 20% RH), respectively. Finally, we investigated the byproducts generated by the humid air discharge itself, which are the greenhouse gases O3, N2O, and the toxic gas NO2.

Abstract: It is known that film structure may change dramatically with the extent of surface diffusion during the film growth process. In the present work, surface diffusion, induced thermally or activated by energetic impacts, is investigated theoretically under conditions appropriate for magnetron sputter-deposition of MgAlO thin films with varying stoichiometry. The distribution of surface diffusion energy barriers available to the system was determined for each stoichiometry, which allowed assessing in a qualitative way how much surface diffusion will take place on the time scale available between deposition events. The activation energy barriers increase with the Al concentration in the film, and therefore, the surface diffusion rates in the time frame of typical deposition rates drop, which can explain the decrease in crystallinity in the film structure and the transition to amorphous structure. The deposition process and the immediate surface diffusion enhanced by the energetic adatoms are simulated by means of a molecular dynamics model. The longer-time thermal surface diffusion and the energy landscape are studied by the temperature accelerated dynamics method, applied in an approximate way. The surface diffusion enhanced by the energetic impacts appears to be very important for the film structure in the low-temperature deposition regime.

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: A two-dimensional model is presented that describes the behavior of argon metastable atoms, copper atoms, and copper ions in an argon direct. current glow discharge, in the standard cell of the VG9000 glow discharge mass spectrometer for analyzing flat samples. The model is combined with a previously developed model for the electrons, argon ions, and atoms in the same cell to obtain an overall picture of the glow discharge, The results of the present model comprise the number densities of the described plasma species, the relative contributions of different production and loss processes for the argon metastable atoms, the thermalization profile of the sputtered copper atoms, the relative importance of the different ionization mechanisms for the copper atoms, the ionization degree of copper, the copper ion-to-argon ion density ratio, and the relative roles of copper ions, argon ions, and atoms in the sputtering process. All these quantities are calculated for a range of voltages and pressures, Moreover, since the sticking coefficient of copper atoms on solid surfaces is not well-known in the literature, the influence of this parameter on the results is briefly discussed.

Abstract: We have investigated the plasma processes in an atmospheric pressure glow discharge (APGD) in He used for analytical spectrometry by means of fluid and Monte Carlo (MC) simulations. Typical results include the potential and electric field distributions in the plasma, the density profiles of the various plasma species throughout the discharge, the mean electron energy, as well as the rates of the various collision processes in the plasma, and the relative importance of the different production and loss rates for the various species. The similarities and differences with low-pressure glow discharges are discussed. The main differences are a very small cathode dark space region and a large positive column as well as the dominant role of molecular ions. Some characteristic features of the APGD, such as the occurrence of the different spatial zones in the discharge, are illustrated, with links to experimental observations.

Abstract: An inductively coupled plasma (ICP) is analyzed by means of experiments and numerical simulation. Important plasma properties are analyzed, namely, the effective temperature inside the central channel and the mean flow velocity inside the plasma. Furthermore, the effect of torches with different injector diameters is studied by the model. The temperature inside the central channel is determined from the end-on collected line-to-background ratio in dependence of the injector gas flow rates. Within the limits of 3% deviation, the results of the simulation and the experiments are in good agreement in the range of flow rates relevant for the analysis of relatively large droplets, i.e., 50 μm. The deviation increases for higher gas flow rates but stays below 6% for all flow rates studied. The velocity of the gas inside the coil region was determined by side-on analyte emission measurements with single monodisperse droplet introduction and by the analysis of the injector gas path lines in the simulation. In the downstream region significantly higher velocities were found than in the upstream region in both the simulation and the experiment. The quantitative values show good agreement in the downstream region. In the upstream region, deviations were found in the absolute values which can be attributed to the flow conditions in that region and because the methods used for velocity determination are not fully consistent. Eddy structures are found in the simulated flow lines. These affect strongly the way taken by the path lines of the injector gas and they can explain the very long analytical signals found in the experiments at low flow rates. Simulations were performed for different injector diameters in order to find conditions where good analyte transport and optimum signals can be expected. The results clearly show the existence of a transition flow rate which marks the lower limit for effective analyte transport conditions through the plasma. A rule-of-thumb equation was extracted from the results from which the transition flow rate can be estimated for different injector diameters and different injector gas compositions.