Abstract: Using reactive molecular dynamics simulations, we have investigated the effect of single-impact, low-energy (thermal-100 eV) bombardment of a Si(100){2 × 1} surface by atomic and molecular oxygen. Penetration probability distributions, as well as defect formation distributions, are presented as a function of the impact energy for both species. It is found that at low impact energy, defects are created chemically due to the chemisorption process in the top layers of the surface, while at high impact energy, additional defects are created by a knock-on displacement of Si. These results are of particular importance for understanding device performances of silica-based microelectronic and photovoltaic devices.

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: Despite significant progress in single walled carbon nanotube (SWCNT) production by plasma enhanced chemical vapor deposition (PECVD), the growth mechanism in this method is not clearly understood. We employ reactive molecular dynamics simulations to investigate how plasma-based deposition allows growth at low temperature. We first investigate the SWCNT growth mechanism at low and high temperatures under conditions similar to thermal CVD and PECVD. We then show how ion bombardment during the nucleation stage increases the carbon solubility in the catalyst at low temperature. Finally, we demonstrate how moderate energy ions sputter amorphous carbon allowing for SWCNT growth at 500 K. (C) 2013 Elsevier B. V. All rights reserved.

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.