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: The compositional influence of Cr and Y on the microstructure of Mg―Cr―O, and Mg―Y―O films synthesized by reactive magnetron sputtering has been investigated by transmission electron microscopy, X-ray diffraction and molecular dynamics simulations. A decrease in crystallinity is observed in these films as the M (Cr or Y) content is increased. It is found that M forms a solid solution with MgO for metal ratios up to ~ 70% and ~ 50% for Cr and Y respectively. Above ~ 70% Cr metal ratio the Mg―Cr―O films are found to be completely amorphous. The Mg―Y―O films are composed of Mg(Y)O and Y2O3 nano crystallites, up to ~ 50% Y metal ratio. Above this ratio, only Y2O3 nano crystallites are found. The preferential < 111> MgO grain alignment is strongly affected by the increase in M content. For M metal ratios up to ~ 50%, there is a selective promotion of the < 100> MgO grain alignments and a decline in the < 111> grain alignments.

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: A model for an atmospheric pressure inductively coupled plasma (ICP) is developed which allows rather easy extension to a variable number of species and ionisation degrees. This encompasses an easy calculation of transport parameters for mixtures, ionisation and heat capacity. The ICP is modeled in an axisymmetric geometry, taking into account the gas streaming into a flowing ambient gas. A mixture of argon and helium is applied in the injector gas stream as it is often done in laser ablation ICP spectrometry. The results show a strong influence of the added helium on the center of the ICP, which is important for chemical analysis. The length of the central channel is significantly increased and the temperature inside is significantly higher than in the case of pure argon. This means that higher gas volume flow rates can be applied by addition of helium compared to the use of pure argon. This has the advantage that the gas velocity in the transport system towards the ICP can be increased, which allows shorter washout-times. Consequently, shorter measurement times can be achieved, e.g. for spatial mapping analyses in laser ablation ICP spectrometry. Furthermore, the higher temperature and the longer effective plasma length will increase the maximum size of droplets or particles injected into the ICP that are completely evaporated at the detection site. Thus, we expect an increase of the analytical performance of the ICP by helium addition to the injector gas.

Abstract: A hybrid model, called the hybrid plasma equipment model, was used to study inductively coupled SF6 plasmas used for Si etching applications. The plasma properties such as number densities of electrons, positive and negative ions, and neutrals are calculated under typical etching conditions. The electron kinetics is analysed by means of the electron energy probability function. The plasma chemistry taking place in pure SF6 and in an Ar/SF6 mixture is also discussed, and finally the effect of the argon fraction on the plasma properties is investigated.

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.

Abstract: This contribution presents the effects of a mild O2 plasma treatment on the structural, optical, and electrical properties of single-layer (SLG) and bilayer graphene (BLG). Unexpectedly, we observe only photoluminescence in the SLG parts of a graphene flake composed of regions of various thickness upon O2 plasma treatment, whereas the BLG and few-layer graphene (FLG) parts remain optically unchanged. Confirmed with X-ray photoelectron spectroscopy (XPS) that O2 plasma induces epoxide and hydroxyl-like groups in graphene, density functional theory (DFT) calculations are carried out on representative epoxidized and hydroxylated SLG and BLG models to predict density of states (DOS) and band structures. Sufficiently oxidized SLG shows a bandgap and thus loss of semimetallic behavior, while oxidized BLG maintains its semimetallic behavior even at high oxygen density in agreement with the results of the photoluminescence spectroscopy (PL) experiments. DFT calculations confirm that the Fermi velocity in epoxidized BLG is remarkably comparable with that of pristine SLG, pointing to a similarity of electronic band structure. The similarity is also experimentally demonstrated by the electrical characterization of a plasma-treated BLG-FET. As expected from the electronegative oxygen adatoms in the graphene, epoxidized BLG presents conductive features typical of hole doping. Moreover, the electrical characteristics suggest band structures closely related to that of epoxidized graphene while deviating from that of hydroxylated graphene. Finally, upon O2 plasma treatment of BLG, the four-component 2D peak around 2700 cm1 in the Raman spectrum evolves into a single Lorentzian line, very like the 2D peak of pristine SLG. Summarizing, the data in this contribution recommend that a controlled O2 plasma treatment, which is compatible with CMOS process flow in contrast to wet chemical oxidation methods, provides an efficient and valuable technique to exploit the transport properties of the bottom layer of BLG.

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 vibrational level populations of the electronic ground state of the nitrogen molecule have been calculated for typical glow discharge conditions in argonnitrogen mixtures with nitrogen concentrations between 0.1 and 1%. Stationary solutions of the master equations of the vibrational levels have been obtained using numerical methods. The main mechanisms responsible for the population and depopulation of the vibrational levels, and for the overall shape of the vibrational distribution function are pointed out. It has been found that vibrationvibration collisions play only a minor role and therefore the population of the vibrational levels is basically determined by the electron temperature.

Abstract: Understanding the fundamental dynamics in carbon nanotube (CNT) catalysts is of primary importance to understand CNT nucleation. This Letter reports on calculated bond switching (BS) rates in pure and carbon containing nickel nanoclusters. The rates are analyzed in terms of their temperature dependent spatial distribution and the mobility of the cluster atoms. The BS mechanism is found to change from vibrational to diffusional at around 900 K, with a corresponding strong increase in activation energy. Furthermore, the BS activation energy is observed to decrease as the carbon content in the cluster increases, resulting in an effective liquification of the cluster.

Abstract: Metal-catalyzed growth mechanisms of carbon nanotubes (CNTs) were studied by hybrid molecular dynamics−Monte Carlo simulations using a recently developed ReaxFF reactive force field. Using this novel approach, including relaxation effects, a CNT with definable chirality is obtained, and a step-by-step atomistic description of the nucleation process is presented. Both root and tip growth mechanisms are observed. The importance of the relaxation of the network is highlighted by the observed healing of defects.

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: The influence of the composition on the crystallographic properties of deposited Mg(M)O (with M=Al, Cr, Ti, Y, and Zr) films is studied. For a flexible control of the composition, dual reactive magnetron sputtering was used as deposition technique. Two different approaches to predict the composition are discussed. The first is an experimental way based on the simple relationship between the deposition rate and the target-substrate distance. The second is a route using a Monte Carlo based particle trajectory code. Both methods require a minimal experimental input and enable the user to quickly predict the composition of complex thin films. Good control and flexibility allow us to study the compositional effects on the growth of Mg(M)O films. Pure MgO thin films were grown with a (111) preferential out-of-plane orientation. When adding M to MgO, two trends were noticed. The first trend is a change in the MgO lattice parameters compared to pure MgO. The second tendency is a decrease in the crystallinity of the MgO phase. The experimentally determined crystallographic properties are shown to be in correspondence with the predicted properties from molecular dynamics simulations.

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: The electronic properties of a set of intrinsic and extrinsic point defects in gamma-Al2O3 are investigated using quasiparticle calculations within the G(0)W(0) approximation. We find that the electronic signature of atomic vacancies lie deep in the band gap, close to the top of the valence band edge. The introduction of C, Si, and N impurities induces defective levels that are located close to the conduction band edge and near the middle of the band gap of the oxide. The comparison with electrical measurements reveals that the energy levels of some of these defects match with the electronic fingerprint of the defects reported in gamma-Al2O3 based nonvolatile memories. (C) 2010 American Institute of Physics. [doi:10.1063/1.3507385]