Abstract: In this work, a detailed structural and spectroscopic study of nanocrystalline diamond (NCD) thin films grown by a continuous bias assisted CVD growth technique is reported. This technique allows the tuning of grain size and phase purity in the deposited material. The crystalline properties of the films are characterized by SEM, TEM, EELS, and Raman spectroscopy. A clear improvement of the crystalline structure of the nanograined diamond film is observed for low negative bias voltages, while high bias voltages lead to thin films consisting of diamond grains of only ∼10 nm nanometer in size, showing remarkable similarities with so-called ultrananocrystalline diamond. These layers arecharacterized by an increasing amount of sp2-bonded carbon content of the matrix in which the diamond grains are embedded. Classical molecular dynamics simulations support the observed experimental data, giving insight in the underlying mechanism for the observed increase in deposition rate with bias voltage. Furthermore, a high atomic concentration of hydrogen has been determined in these films. Finally, Raman scattering analyses confirm that the Raman line observed at ∼1150 cm−1 cannot be attributed to trans-poly-acetylene, which continues to be reported in literature, reassigning it to a deformation mode of CHx bonds in NCD.

Abstract: A crucial step, which is still not well understood in the destruction of volatile organic compounds (VOCs) with low temperature plasmas, is the initiation of the process. Here, we present a kinetic model for the destruction of ethylene in low temperature plasmas that allows us to calculate the relative importance of all plasma species and their related reactions. Modifying the ethylene concentration and/or the SED had a major impact on the relative importance of the radicals (i.e., mainly atomic oxygen) and the metastable nitrogen (i.e., more specifically N2(equation image)) in the destruction process. Our results show that the direct destruction by electron impact reactions for ethylene can be neglected; however, we can certainly not neglect the influence of N2(equation image)).

Abstract: A Langmuir probe and an ICCD are employed to study the discharge mode transition in Ar inductively coupled plasma. Electron density and plasma emission intensity are measured during the E (capacitive discharge) to H (inductive discharge) mode transitions at different pressures. It is found that plasma exists with a low electron density and a weak emission intensity in the E mode, while it has a high electron density and a strong emission intensity in the H mode. Meanwhile, the plasma emission intensity spatial (2D image) profile is symmetrical in the H mode, but the 2D image is an asymmetric profile in the E mode. Moreover, the electron density and emission intensity jump up discontinuously at high pressure, but increase almost continuously at the E to H mode transition under low pressure.

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: A theoretical study has been conducted of ultrashort pulsed laser induced electron emission from an aluminum surface. Electron emission fluxes retrieved from the commonly employed Fowler-DuBridge theory were compared to fluxes based on a laser-induced non-equilibrium electron distribution. As a result, the two-and three-photon photoelectron emission parameters for the Fowler-DuBridge theory have been approximated. We observe that at regimes where photoemission is important, laser-induced electron emission evolves in a more smooth manner than predicted by the Fowler-DuBridge theory. The importance of the actual electron distribution decreases at higher laser fluences, whereas the contribution of thermionic emission increases. Furthermore, the influence of a space charge effect on electron emission was evaluated by a one dimensional particle-in-cell model. Depending on the fluences, the space charge reduces the electron emission by several orders of magnitude. The influence of the electron emission flux profiles on the effective electron emission was found to be negligible. However, a non-equilibrium electron velocity distribution increases the effective electron emission significantly. Our results show that it is essential to consider the non-equilibrium electron distribution as well as the space charge effect for the description of laser-induced photoemission. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4729071]

Abstract: Uniform acceptance force biased Monte Carlo (UFMC) simulations have previously been shown to be a powerful tool to simulate atomic scale processes, enabling one to follow the dynamical path during the simulation. In this contribution, we present a simple proof to demonstrate that this uniform acceptance still complies with the condition of detailed balance, on the condition that the characteristic parameter lambda = 1/2 and that the maximum allowed step size is chosen to be sufficiently small. Furthermore, the relation to Metropolis Monte Carlo (MMC) is also established, and it is shown that UFMC reduces to MMC by choosing the characteristic parameter lambda = 0 [Rao, M. et al. Mol. Phys. 1979, 37, 1773]. Finally, a simple example compares the UFMC and MMC methods.

Abstract: The possibility of application of nanomaterials is determined by our ability to control the properties of the materials, which are ultimately determined by their structure and hence their growth processes. We employ hybrid molecular dynamics / Monte Carlo (MD/MC) simulations to explore the growth of SWNTs and graphene on nickel as a catalyst, with the specific goal of unraveling the growth mechanisms. While the general observations are in agreement with the literature, we find a number of interesting phenomena to be operative which are crucial for the growth, and which are not accessible by MD simulations alone due to the associated time scale. Specifically, we observe metal mediated healing and restructuring processes to take place, reorganizing the carbon network during the initial nucleation step. In the case of carbon nanotube growth, this leads to the growth of tubes with a determinable chirality. In the case of graphene formation, we find that graphene is only formed at temperatures above 700 K. These results are of importance for understanding the growth mechanisms of these carbon nanomaterials on the fundamental level.

Abstract: Monte Carlo (MC) methods have a long-standing history as partners of molecular dynamics (MD) to simulate the evolution of materials at the atomic scale. Among these techniques, the uniform-acceptance force-bias Monte Carlo (UFMC) method [ G. Dereli Mol. Simul. 8 351 (1992)] has recently attracted attention [ M. Timonova et al. Phys. Rev. B 81 144107 (2010)] thanks to its apparent capacity of being able to simulate physical processes in a reduced number of iterations compared to classical MD methods. The origin of this efficiency remains, however, unclear. In this work we derive a UFMC method starting from basic thermodynamic principles, which leads to an intuitive and unambiguous formalism. The approach includes a statistically relevant time step per Monte Carlo iteration, showing a significant speed-up compared to MD simulations. This time-stamped force-bias Monte Carlo (tfMC) formalism is tested on both simple one-dimensional and three-dimensional systems. Both test-cases give excellent results in agreement with analytical solutions and literature reports. The inclusion of a time scale, the simplicity of the method, and the enhancement of the time step compared to classical MD methods make this method very appealing for studying the dynamics of many-particle systems.

Abstract: Transition metal oxide-based resistor random access memory (RRAM) takes advantage of oxygen-related defects in its principle of operation. Since the change in resistivity of the material is controlled by the oxygen deficiency level, it is of major importance to quantify the kinetics of the oxygen diffusion, key factor for oxide stoichiometry. Ab initio accelerated molecular dynamics techniques are employed to investigate the oxygen diffusivity in amorphous hafnia (HfOx, x = 1.97, 1.0, 0.5). The computed kinetics is in agreement with experimental measurements. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.3697690]

Abstract: A numerical model for NxOy remediation in humid air plasma produced with a dielectric barrier discharge at atmospheric pressure is presented. Special emphasis is given to NO2 and N2O reduction with the decrease of O2 content in the feedstock gas. A detailed reaction mechanism including electronic and ionic processes, as well as the contribution of radicals and excited atomic/molecular species is proposed. The temporal evolution of the densities of NO, NO2 and N2O species, and some other by-products, is analyzed, and the major pathways for the NxOy remediation are discussed for one pulse. Subsequently, simulations are presented for a multi-pulses case, where three O2 contents are tested for optimization of the remediation process. It is found that when the gas mixture O2/N2/H2O/NOx has no initial O2 content, the best NOx and N2O remediation is achieved.

Abstract: Modeling results are presented to gain a better insight in the properties of a SiH4/O2/Ar inductively coupled plasma (ICP) and how it interacts with a silicon substrate (wafer), as applied in the microelectronics industry for the fabrication of electronic devices. The SiH4/O2/Ar ICP is used for the filling of microtrenches with isolating material (SiO2), as applied in shallow trench isolation (STI). In this article, a detailed reaction set that describes the plasma chemistry of SiH4/O2/Ar discharges as well as surface processes, such as sputtering, oxidation, and deposition, is presented. Results are presented on the plasma properties during the plasma enhanced chemical vapor deposition process (PECVD) for different gas ratios, as well as on the shape of the filled trenches and the surface compositions of the deposited layers. For the operating conditions under study it is found that the most important species accounting for deposition are SiH2, SiH3O, SiH3 and SiH2O, while SiH+2, SiH+3, O+2 and Ar+ are the dominant species for sputtering of the surface. By diluting the precursor gas (SiH4) in the mixture, the deposition rate versus sputtering rate can be controlled for a desired trench filling process. From the calculation results it is clear that a high deposition rate will result in undesired void formation during the trench filling, while a small deposition rate will result in undesired trench bottom and mask damage by sputtering. By varying the SiH4/O2 ratio, the chemical composition of the deposited layer will be influenced. However, even at the highest SiH4/O2 ratio investigated (i.e., 3.2:1; low oxygen content), the bulk deposited layer consists mainly of SiO2, suggesting that low-volatile silane species deposit first and subsequently become oxidized instead of being oxidized first in the plasma before deposition. Finally, it was found that the top surface of the deposited layer contained less oxygen due to preferential sputtering of O atoms, making the top layer more Si-rich. However, this effect is negligible at a SiH4/O2 ratio of 2:1 or lower.

Abstract: Ultrathin bismuth titanate films (Bi2Ti2O7, 5-25 nm) are deposited onto SiO2/Si substrates by aqueous chemical solution deposition and their evolution during annealing is studied. The films crystallize into a preferentially oriented, pure pyrochlore phase between 500 and 700 degrees C, depending on the film thickness and the total thermal budget. Crystallization causes a strong increase of surface roughness compared to amorphous films. An increase of the interfacial layer thickness is observed after anneal at 600 degrees C, together with intermixing of bismuth with the substrate as shown by TEM-EDX. The band gap was determined to be similar to 3 eV from photoconductivity measurements and high dielectric constants between 30 and 130 were determined from capacitance voltage measurements, depending on the processing conditions. (C) 2012 Elsevier Ltd. All rights reserved.

Abstract: A two-dimensional fluid model, including the full set of Maxwell equations, has been developed and applied to investigate the effect of a phase shift between two power sources on the radial uniformity of several plasma characteristics in a hydrogen capacitively coupled plasma. This study was carried out at various frequencies in the range 13.56200 MHz. When the frequency is low, at 13.56 MHz, the plasma density is characterized by an off-axis peak when both power sources are in-phase (phgr = 0), and the best radial uniformity is obtained at phgr = π. This trend can be explained because the radial nonuniformity caused by the electrostatic edge effect can be effectively suppressed by the phase-shift effect at a phase difference equal to π. When the frequency rises to 60 MHz, the plasma density profiles shift smoothly from edge-peaked over uniform to centre-peaked as the phase difference increases, due to the pronounced standing-wave effect, and the best radial uniformity is reached at phgr = 0.3π. At a frequency of 100 MHz, a similar behaviour is observed, except that the maximum of the plasma density moves again towards the radial edge at the reverse-phase case (phgr = π), because of the dominant skin effect. When the frequency is 200 MHz, the bulk plasma density increases significantly with increasing phase-shift values, and a better uniformity is obtained at phgr = 0.4π. This is because the density in the centre increases faster than at the radial edge as the phase difference rises, due to the increasing power deposition Pz in the centre and the decreasing power density Pr at the radial edge. As the phase difference increases to π, the maximum near the radial edge becomes obvious again. This is because the skin effect has a predominant influence on the plasma density under this condition, resulting in a high density at the radial edge. Moreover, the axial ion flux increases monotonically with phase difference, and exhibits similar profiles to the plasma density. The calculation results illustrate that the radial uniformity of the various plasma characteristics is strongly dependent on the applied frequency and the phase shift between both power sources, which is important to realize, for controlling the uniformity of the plasma etch and deposition processes.

Abstract: Simulation results of electrostatics in Si cylindrical junctionless nanowire transistors with a homogenous channel are presented. Junctionless transistors including strain and arbitrary crystallographic orientations are studied. Size quantization effects are simulated by self-consistent solutions of the Poisson and Schrodinger equations. The 6 x 6 k.p method is employed for the calculation of the valence subband structure in a junctionless nanowire pFET. The influence of stress/strain and crystallographic channel orientation on to the electrostatics in terms of subband structure, charge density, and C-V curve is systematically studied. (C) 2011 Elsevier Ltd. All rights reserved.

Abstract: The electronic properties of two-dimensional honeycomb structures of molybdenum disulfide (MoS(2)) subjected to biaxial strain have been investigated using first-principles calculations based on density functional theory. On applying compressive or tensile bi-axial strain on bi-layer and mono-layer MoS(2), the electronic properties are predicted to change from semiconducting to metallic. These changes present very interesting possibilities for engineering the electronic properties of two-dimensional structures of MoS(2).