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: Nanostructured carbon materials show a tremendous variety in atomic structure, morphology, properties, and applications. As all properties are ultimately determined by the structure of the material, a thorough understanding of the growth mechanisms that give rise to the particular structure is critical. On many occasions, it has been shown that plasma enhanced growth can be strongly beneficial. This review will describe the authors current understanding of plasma enhanced growth of carbon nanotubes, the prototypical example of nanostructured carbon materials, as obtained from experiments, simulations, and modeling. Specific emphasis is put on where experiments and computational approaches correspond, and where they differ. Also, the current status on simulating PECVD growth of some other carbon nanomaterials is reviewed, including amorphous carbon, graphene, and metallofullerenes. Finally, computational challenges with respect to the simulation of PECVD growth are identified.

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: 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: In this letter, CMOS-compatible Ni/HfO2/TiN resistive random access memory stacks demonstrated attractive unipolar switching properties, showing >10(3) endurance and long retention at 150 degrees C. The Ni bottom electrode (BE) improved the switching yield over the NiSiPt BE. To better understand the unipolar forming mechanism, ab initio simulation and time of flight-secondary ion mass spectroscopy were utilized. Compared to the NiSiPt BE, Ni BE gives larger Ni diffusion in the HfO2 and lower formation enthalpy of Ni2+ species during electrical forming. Both the electrical and physical results supported a Ni-injection mechanism for the filament formation. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.3695078]

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: 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: In order to minimize the dead volume in large cells for laser ablation inductively coupled plasma mass spectrometry, and improve the aerosol entrainment characteristics, the gas inlet nozzle has been set in rotation. This allowed a wider volume to be swept than with the traditional static inlet nozzle approach. Therefore, sensitivity combined with site-to-site repeatability was improved by a factor of two, together with minimization of aerosol loss within the cell and signal dispersion.