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: Germanium combined with high-κ dielectrics is investigated for the next generations of CMOS devices. Therefore, we study reaction mechanisms for Al2O3 atomic layer deposition on sulfur passivated Ge using calculations based on density functional theory and total reflection X-ray fluorescence (TXRF). TXRF indicates 6 S/nm2 and 4 Al/nm2 after the first TMA/H2O reaction cycle, and growth inhibition from the second reaction cycle on. Calculations are performed on molecular clusters representing −GeSH surface sites. The calculations confirm that the TMA reaction does not affect the S content. On fully SH-terminated Ge, TMA favorably reacts with up to three −GeSH sites, resulting in a near tetrahedral Al coordination. Electron deficient structures with a GeS site shared between two Al atoms are proposed. The impact of the cluster size on the structures and reaction energetics is systematically investigated.

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: A molecular dynamics (MD) study of MgxAlyOz thin films grown by magnetron sputtering is presented using an ionic model and comparing two potential sets with formal and partial charges. The applicability of the model and the reliability of the potential sets for the simulation of thin film growth are discussed. The formal charge potential set was found to reproduce the thin film structure in close agreement with the structure of the experimentally grown thin films. Graphical abstract A molecular dynamics study of growth of MgxAlyOz thin films is presented using an ionic model and comparing two potential sets with formal and partial charges. The simulation results with the formal charge potential set showed a transition in the film from a crystalline to an amorphous structure, when the Mg metal content decreases below 50% in very close agreement with the structure of the experimentally deposited films.

Abstract: The reaction behavior of species that may affect the growth of ultrananocrystal line and nanocrystalline diamond ((U)NCD) films is investigated by means of molecular dynamics simulations. Impacts of CHx (x = 0 – 4), C2Hx (x=0-6), C3Hx (x=0-2), C4Hx (x = 0 – 2), H, and H-2 on clean and hydrogenated diamond (100)2 x 1 and (111) 1 x 1 surfaces at two different substrate temperatures are simulated. We find that the different bonding structures of the two surfaces cause different temperature effects on the sticking efficiency. These results predict a temperature-dependent ratio of diamond (100) and (111) growth. Furthermore, predictions of which are the most important hydrocarbon species for (U)NCD growth are made.

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: 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.

Abstract: RuCp2 (ruthenocene) and RuCpPy (cyclopentadienyl pyrrolyl ruthenium) complexes are used in ruthenium (Ru) atomic layer deposition (ALD) but exhibit a markedly different reactivity with respect to the substrate and co-reactant. In search of an explanation, we report here the results of a comparative study of the heterolytic and homolytic dissociation enthalpy of these two ruthenium complexes, making use of either density functional theory (DFT) or multiconfigurational perturbation theory (CASPT2). While both methods predict distinctly different absolute dissociation enthalpies, they agree on the relative values between both molecules. A reduced heterolytic dissociation enthalpy is obtained for RuCpPy compared to RuCp2, although the difference obtained from CASPT2 (19.9 kcal/mol) is slightly larger than the one obtained with any of the DFT functionals (around 17 kcal/mol). Both methods also agree on the more pronounced stability of the Cp- ligand in RuCpPy than in RuCp2 (by around 9 kcal/mol with DFT and by 6 kcal/mol with CASPT2).

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.