Abstract: The phase transition processes induced by ultrashort, 100 fs pulsed laser irradiation of Au, Cu, and Ni are studied by means of a combined atomistic-continuum approach. A moderately low absorbed laser fluence range, from 200 to 600 J/m2 is considered to study phase transitions by means of a local and a nonlocal order parameter. At low laser fluences, the occurrence of layer-by-layer evaporation has been observed, which suggests a direct solid to vapor transition. The calculated amount of molten material remains very limited under the conditions studied, especially for Ni. Therefore, our results show that a kinetic equation that describes a direct solid to vapor transition might be the best approach to model laser-induced phase transitions by continuum models. Furthermore, the results provide more insight into the applicability of analytical superheating theories that were implemented in continuum models and help the understanding of nonequilibrium phase transitions.

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: We present a computational study for the conversion of CH4 and CO2 into value-added chemicals, i.e., the so-called dry reforming of methane, in a dielectric barrier discharge reactor. A zero-dimensional chemical kinetics model is applied to study the plasma chemistry in a 1:1 CH4/CO2 mixture. The calculations are first performed for one microdischarge pulse and its afterglow, to study in detail the chemical pathways of the conversion. Subsequently, long time-scale simulations are carried out, corresponding to real residence times in the plasma, assuming a large number of consecutive microdischarge pulses, to mimic the conditions of the filamentary discharge regime in a dielectric barrier discharge (DBD) reactor. The conversion of CH4 and CO2 as well as the selectivity of the formed products and the energy cost and energy efficiency of the process are calculated and compared to experiments for a range of different powers and gas flows, and reasonable agreement is reached.

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: The growth mechanism and chirality formation of a single-walled carbon nanotube (SWNT) on a surface-bound nickel nanocluster are investigated by hybrid reactive molecular dynamics/force-biased Monte Carlo simulations. The validity of the interatomic potential used, the so-called ReaxFF potential, for simulating catalytic SWNT growth is demonstrated. The SWNT growth process was found to be in agreement with previous studies and observed to proceed through a number of distinct steps, viz., the dissolution of carbon in the metallic particle, the surface segregation of carbon with the formation of aggregated carbon clusters on the surface, the formation of graphitic islands that grow into SWNT caps, and finally continued growth of the SWNT. Moreover, it is clearly illustrated in the present study that during the growth process, the carbon network is continuously restructured by a metal-mediated process, thereby healing many topological defects. It is also found that a cap can nucleate and disappear again, which was not observed in previous simulations. Encapsulation of the nanoparticle is observed to be prevented by the carbon network migrating as a whole over the cluster surface. Finally, for the first time, the chirality of the growing SWNT cap is observed to change from (11,0) over (9,3) to (7,7). It is demonstrated that this change in chirality is due to the metal-mediated restructuring process.

Abstract: In this paper, the effects of prolonged bias-enhanced nucleation (prolonged BEN) on the growth mechanisms of diamond are investigated by molecular dynamics (MD) and combined MD-Metropolis Monte Carlo (MD-MMC) simulations. First, cumulative impacts of CxHy+ and Hx+ on an a-C:H/nanodiamond composite were simulated; second, nonconsecutive impacts of the dominant ions were simulated in order to understand the observed phenomena in more detail. As stated in the existing literature, the growth of diamond structures during prolonged BEN is a process that takes place below the surface of the growing film. The investigation of the penetration behavior of CxHy+ and Hx+ species shows that the carbon-containing ions remain trapped within this amorphous phase where they dominate mechanisms like precipitation of sp3 carbon clusters. The H+ ions, however, penetrate into the crystalline phase at high bias voltages (>100 V), destroying the perfect diamond structure. The experimentally measured reduction of grain sizes at high bias voltage, reported in the literature, might thus be related to penetrating H+ ions. Furthermore, the CxHy+ ions are found to be the most efficient sputtering agents, preventing the build up of defective material.

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.

Abstract: Nonthermal plasmas are currently being studied as a green alternative to the Haber-Bosch process, which is, today, the dominant industrial process allowing for the fixation of nitrogen and, as such, a fundamental component for the production of nitrogen-based industrial fertilizers. In this context, the gliding arc plasma (GAP) is considered a promising choice among nonthermal plasma options. However, its stability is still a key parameter to ensure industrial transfer of the technology. Nowadays, the conventional approach to stabilize this plasma process is to use external resistors. Although this indeed allows for an enhancement of the plasma stability, very little is reported about how it impacts the process efficiency, both in terms of NOx yield and energy cost. In this work, this question is specifically addressed by studying a DC-powered GAP utilized for nitrogen fixation into NOx at atmospheric pressure stabilized by variable external resistors. Both the performance and the stability of the plasma are reported as a function of the utilization of the resistors. The results confirm that while the use of a resistor indeed allows for a strong stabilization of the plasma without impacting the NOx yield, especially at high plasma current, it dramatically impacts the energy cost of the process, which increases from 2.82 to 7.9 MJ/mol. As an alternative approach, we demonstrate that the replacement of the resistor by an inductor is promising since it allows for decent stabilization of the plasma, while it does not affect either the energy cost of the process or the NOx yield.

Abstract: An optimized laser ablation setup, proposed for high repetition rate inductively coupled plasma mass spectrometry (ICPMS) analyses such as 2D imaging or depth profiling, is presented. For such applications, the particle washout time needs to be as short as possible to allow high laser pulse frequencies for reduced analysis time. Therefore, it is desirable to have an ablation setup that operates as a laminar flow reactor (LFR). A top-down strategy was applied that resulted in the present design. In the first step, a previously applied ablation setup was analyzed on the basis of computational fluid dynamics (CFD) results presented by D. Autrique et al. (Spectrochim. Acta, B 2008, 63, 257−270). By means of CFD simulations, the design was modified in such a way that it operated in the LFR regime. Experimental results demonstrate that the current design can indeed be regarded as an LFR. Furthermore, the operation under LFR conditions allowed some insight into the initial radial concentration distribution if the experimental ICPMS signal and analytical expressions are taken into account. Recommendations for a modified setup for more resilient spatial distributions are given. With the present setup, a washout time of 140 ms has been achieved for a 3% signal area criterion. Therefore, 7 Hz repetition rates can be applied with the present setup. Using elementary formulas of the analytical model, an upper bound for the washout times for similar setups can be predicted. The authors believe that the presented setup geometry comes close to the achievable limit for reliable short washout times.