Abstract: Using reactive molecular dynamics simulations by means of the ReaxFF potential, we studied the growth mechanism of ultrathin silica (SiO2) layers during hyperthermal oxidation at room temperature. Oxidation of Si(100){2 × 1} surfaces by both atomic and molecular oxygen was investigated in the energy range 15 eV. The oxidation mechanism, which differs from thermal oxidation, is discussed. In the case of oxidation by molecular O2, silica is quickly formed and the thickness of the formed layers remains limited compared to oxidation by atomic oxygen. The Si/SiO2 interfaces are analyzed in terms of partial charges and angle distributions. The obtained structures of the ultrathin SiO2 films are amorphous, including some intrinsic defects. This study is important for the fabrication of silica-based devices in the micro- and nanoelectronics industry, and more specifically for the fabrication of metal oxide semiconductor devices.

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: 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: Since the discovery of graphene, a broad range of two-dimensional (2D) materials has captured the attention of the scientific communities. Materials, such as hexagonal boron nitride (hBN) and the transition metal dichalcogenides (TMDs) family, have shown promising semiconducting and insulating properties that are very appealing for the semiconductor industry. Recently, the possibility of taking advantage of the properties of 2D-based heterostructures has been investigated for low-power nanoelectronic applications. In this work, we aim at evaluating the relation between the nature of the materials used in such heterostructures and the amplitude of the layer-to-layer charge transfer induced by an external electric field, as is typically present in nanoelectronic gated devices. A broad range of combinations of TMDs, graphene, and hBN has been investigated using density functional theory. Our results show that the electric field induced charge transfer strongly depends on the nature of the 2D materials used in the van der Waals heterostructures and to a lesser extent on the relative orientation of the materials in the structure. Our findings contribute to the building of the fundamental understanding required to engineer electrostatically the doping of 2D materials and to establish the factors that drive the charge transfer mechanisms in electron tunneling-based devices. These are key ingredients for the development of 2D -based nanoelectronic devices.

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 and vibrational properties of 2D honeycomb structures of molybdenum disulfide (MoS2) subjected to strain have been investigated using first-principles calculations based on density functional theory. We have studied the evolution of the electronic properties of bulk and layered MoS2, going down from a few layers up to a mono-layer, and next investigated the effect of bi-axial strain on their electronic structure and vibrational frequencies. Both for tensile and compressive biaxial strains, the shrinking of the energy band-gap of MoS2 with increasing level of applied strain is observed and a transition limit of the system from semiconducting to metallic is predicted to occur for strains in the range of 8-10%. We also found a progressive downshift (upshift) of both the E-2g(1) and A(1g) Raman active modes with increasing level of applied tensile (compressive) strain. Interestingly, significant changes in the curvature of the conduction and valence band near their extrema upon the application of strain are also predicted, with correlated variations of the electron and hole effective masses. These changes present interesting possibilities for engineering the electronic properties of 2D structures of MoS2. (C) 2012 Elsevier B.V. All rights reserved.

Abstract: First-principles simulations were employed to gain atomistic insights on the working principles of amorphous HfO2 based Resistive Random Access Memory stack: the nature of the defect responsible for the switching between the High and Low Resistive States has been unambiguously identified to be the substoichiometric Hf sites (commonly called oxygen vacancy-V-O) and the kinetics of the process have been investigated through the study of O diffusion. Also the role of each material layer in the TiN/HfO2/Hf/TiN RRAM stack and the impact of the deposition techniques have been examined: metallic Hf sputtering is needed to provide an oxygen exchange layer that plays the role of defect buffer. TiN shall be a good defect barrier for O but a bad defect buffer layer. A possible scenario to explain the device degradation (switching failure) mechanism has been proposed – the relaxation of the metastable amorphous phase towards crystalline structure leads to denser, more structured cluster that can increase the defect migration barriers. (C) 2013 Elsevier B.V. All rights reserved.

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: Silicene and germanene are the silicon and germanium counterparts of graphene, respectively. Recent experimental works have reported the growth of silicene on (1 1 1)Ag surfaces with different atomic configurations, depending on the growth temperature and surface coverage. We first theoretically study the structural and electronic properties of silicene on (1 1 1) Ag surfaces, focusing on the (4 x 4) silicene/Ag structure. Due to symmetry breaking in the silicene layer (nonequivalent number of top and bottom Si atoms), the corrugated silicene layer, with the Ag substrate removed, is predicted to be semiconducting, with a computed energy bandgap of about 0.3 eV. However, the hybridization between the Si 3p orbitals and the Ag 5s orbital in the silicene/(1 1 1)Ag slab model leads to an overall metallic system, with a distribution of local electronic density of states, which is related to the slightly disordered structure of the silicene layer on the (1 1 1)Ag surface. We next study the interaction of silicene and germanene with different hexagonal non-metallic substrates, namely ZnS and ZnSe. On reconstructed (0 0 0 1)ZnS or ZnSe surfaces, which should be more energetically stable for very thin layers, silicene and germanene are found to be semiconducting. Remarkably, the nature and magnitude of their energy bandgap can be controlled by an out-of-plane electric field, an important finding for the potential use of these materials in nanoelectronic devices. (C) 2013 Elsevier B. V. All rights reserved.

Abstract: The electronic and vibrational properties of three different reconstructions of silicene on Ag(1 1 1) are calculated and compared to experimental results. The 2D epitaxial silicon layers, namely the (4 x 4), (root 13 x root 13) and (2 root 3 x 2 root 3) phases, exhibit different electronic and vibrational properties. Few peaks in the experimental Raman spectrum are identified and attributed to the vibrational modes of the silicene layers. The position and behavior of the Raman peaks with respect to the excitation energy are shown to be a fundamental tool to investigate and discern different phases of silicene on Ag( 1 1 1). (C) 2013 Elsevier B.V. All rights reserved.

Abstract: The structural and vibrational properties of two-dimensional hexagonal silicon (silicene) and germanium (germanene) are investigated by means of first-principles calculations. It is predict that the silicene (germanene) structure with a small buckling of 0.44 (0.7 ) and bond lengths of 2.28 (2.44 ) is energetically the most favorable, and it does not exhibit imaginary phonon mode. The calculated non-resonance Raman spectra of silicene is characterized by a main peak at about 575 cm(-1), namely the G-like peak. For germanene, the highest peak is at about 290 cm(-1). Extensive calculations on armchair silicene nanoribbons and armchair germanene nanoribbons are also performed, with and without hydrogenation of the edges. The studies reveal other Raman peaks mainly distributed at lower frequencies than the G-like peak which could be attributed to the defects at the edges of the ribbons, thus not present in the Raman spectra of non-defective silicene and germanene. Particularly the Raman peak corresponding to the D mode is found to be located at around 515 cm(-1) for silicene and 270 cm(-1) for germanene. The calculated G-like and the D peaks are likely the fingerprints of the Raman spectra of the low-buckled structures of silicene and germanene.

Abstract: <script type='text/javascript'>document.write(unpmarked('The resistive random-access memory (RRAM) device concept is close to enabling the development of a new generation of non-volatile memories, provided that their reliability issues are properly understood. The design of a RRAM operating with extrinsic defects based on metallic inclusions, also called conductive bridge RAM, allows the use of a large spectrum of solid electrolytes. However, when scaled to device dimensions that meet the requirements of the latest technological nodes, the discrete nature of the atomic structure of the materials impacts the device operation. Using density functional theory simulations, we evaluated the migration kinetics of Cu conducting species in amorphous and solid electrolyte materials, and established that atomic disorder leads to a large variability in terms of defect stability and kinetic barriers. This variability has a significant impact on the filament resistance and its dynamics, as evidenced during the formation step of the resistive filament. Also, the atomic configuration of the formed filament can age/relax to another metastable atomic configuration, and lead to a modulation of the resistivity of the filament. All these observations are qualitatively explained on the basis of the computed statistical distributions of the defect stability and on the kinetic barriers encountered in RRAM materials.'));

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: Resistive random access memories are emerging as a new type of memory that has the potential to combine both the speed of volatile and the retention of nonvolatile memories. It operates based on the formation/dissolution of a low-resistivity filament being constituted of either metallic ions or atomic vacancies within an insulating matrix. At present, the mechanisms and the parameters controlling the performances of the device remain unclear. In that respect, first-principles simulations provide useful insights on the atomistic mechanisms, the thermodynamic and kinetics factors that modulate the material conductivity, providing guidance into the engineering of the operation of the device. In this paper, we review the current state-of-the-art knowledge on the atomistic switching mechanisms driving the operation of copper-based conductive bridge RRAM and HfOx valence change RRAM. [GRAPHICS] Conceptual illustration of the RRAM device with the filament formation and disruption during its operation. AE/IM/CE are the active electrode/insulating matrix/counterelectrode. The blue circles represent the conducting defects. (C) 2014 WILEY-VCH Verlag GmbH Co. KGaA, Weinheim

Abstract: The amorphous oxide semiconductor Indium-Gallium-Zinc-Oxide (a-IGZO) has gained a large technological relevance as a semiconductor for thin-film transistors in active-matrix displays. Yet, major questions remain unanswered regarding the atomic origin of threshold voltage control, doping level, hysteresis, negative bias stress (NBS), and negative bias illumination stress (NBIS). We undertake a systematic study of the effects of oxygen vacancies on the properties of a-IGZO by relating experimental observations to microscopic insights gained from first-principle simulations. It is found that the amorphous nature of the semiconductor allows unusually large atomic relaxations. In some cases, oxygen vacancies are found to behave as perfect shallow donors without the formation of structural defects. Once structural defects are formed, their transition states can vary upon charge and discharge cycles. We associate this phenomenon to a possible presence of hysteresis in the transfer curve of the devices. Under NBS, the creation of oxygen vacancies becomes energetically very stable, hence thermodynamically very likely. This generation process is correlated with the occurrence of the negative bias stress instabilities observed in a-IGZO transistors. While oxygen vacancies can therefore be related to NBS and hysteresis, it appears unlikely from our results that they are direct causes of NBIS, contrary to common belief.