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: Continued miniaturization of metal-oxide-semiconductor field-effect transistors (MOSFETs) requires an ever-decreasing thickness of the gate oxide. The structure of ultrathin silicon oxide films, however, critically depends on the oxidation mechanism. Using reactive atomistic simulations, we here demonstrate how the oxidation mechanism in hyperthermal oxidation of such structures may be controlled by the oxidation temperature and the oxidant energy. Specifically, we study the interaction of hyperthermal oxygen with energies of 15 eV with thin SiOx (x ≤ 2) films with a native oxide thickness of about 10 Å. We analyze the oxygen penetration depth probability and compare with results of the hyperthermal oxidation of a bare Si(100){2 × 1} (c-Si) surface. The temperature-dependent oxidation mechanisms are discussed in detail. Our results demonstrate that, at low (i.e., room) temperature, the penetrated oxygen mostly resides in the oxide region rather than at the SiOx|c-Si interface. However, at higher temperatures, starting at around 700 K, oxygen atoms are found to penetrate and to diffuse through the oxide layer followed by reaction at the c-Si boundary. We demonstrate that hyperthermal oxidation resembles thermal oxidation, which can be described by the DealGrove model at high temperatures. Furthermore, defect creation mechanisms that occur during the oxidation process are also analyzed. This study is useful for the fabrication of ultrathin silicon oxide gate oxides for metal-oxide-semiconductor devices as it links parameters that can be straightforwardly controlled in experiment (oxygen temperature, velocity) with the silicon oxide structure.

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: We study the evolution of the structural and electronic properties of crystalline indium gallium zinc oxide (IGZO) upon amorphization by first-principles calculation. The bottom of the conduction band (BCB) is found to be constituted of a pseudo-band of molecular orbitals that resonate at the same energy on different atomic sites. They display a bonding character between the s orbitals of the metal sites and an anti-bonding character arising from the interaction between the oxygen and metal s orbitals. The energy level of the BCB shifts upon breaking of the crystal symmetry during the amorphization process, which may be attributed to the reduction of the coordination of the cationic centers. The top of the valence band (TVB) is constructed from anti-bonding oxygen p orbitals. In the amorphous state, they have random orientation, in contrast to the crystalline state. This results in the appearance of localized tail states in the forbidden gap above the TVB. Zinc is found to play a predominant role in the generation of these tail states, while gallium hinders their formation. Last, we study the dependence of the fundamental gap and effective mass of IGZO on mechanical strain. The variation of the gap under strain arises from the enhancement of the anti-bonding interaction in the BCB due to the modification of the length of the oxygen-metal bonds and/or to a variation of the cation coordination. This effect is less pronounced for the amorphous material compared to the crystalline material, making amorphous IGZO a semiconductor of choice for flexible electronics. Finally, the effective mass is found to increase upon strain, in contrast to regular materials. This counterintuitive variation is due to the reduction of the electrostatic shielding of the cationic centers by oxygen, leading to an increase of the overlaps between the metal orbitals at the origin of the delocalization of the BCB. For the range of strain typically met in flexible electronics, the induced variation in the effective mass is found to be negligible (less than 1%).

Abstract: In this work, we have studied the TMA/H(2)O (TMA Al(CH(3))(3)) atomic layer deposition (ALD) of Al(2)O(3) on hydroxyl (OH) and thiol (SH) terminated semiconductor substrates. Total reflection x-ray fluorescence reveals a complex growth-per-cycle evolution during the early ALD reaction cycles. OH and SH terminated surfaces demonstrate growth inhibition from the second reaction cycle on. Theoretical calculations, based on density functional theory, are performed on cluster models to investigate the first TMA/H(2)O reaction cycle. Based on the theoretical results, we discuss possible mechanisms for the growth inhibition from the second reaction cycle on. In addition, our calculations show that AlCH(3) groups are hydrolyzed by a H(2)O molecule adsorbed on a neighboring Al atom, independent of the type of backbonds (Si-O, Ge-O, or Ge-S) of AlCH(3). The coordination of Al remains four-fold after the first TMA/H(2)O reaction cycle. (C) 2012 American Vacuum Society. [DOI: 10.1116/1.3664090]

Abstract: The impact of interfacial chemistry occurring at dielectric/gate interface of P-MOS and N-MOS devices is reviewed through a quick literature survey. A specific emphasis is put on the way the bond polarization that occurs between a dielectric and a metal substrate impacts on the gate work function. First-principle simulations are then used to study the work function changes induced by dopant aggregation in nickel monosilicide metal gates. It is shown that the changes are a natural consequence of the variation of the interface polarization.

Abstract: Nature of the mobility-gap states in amorphous Ge-rich Ge50Se50 was found to be related to homopolar Ge bonds in the chains/clusters of Ge atoms. Threshold switching material suffers Ge-Ge bond concentration drift during material ageing, which can explain the observed reliability of the aGe(50)Se(50) selector devices. Strong Ge-N bonds were introduced to alleviate the observed instability.

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: The application of coreshell SiSiO2 nanowires as nanoelectronic devices strongly depends on their structure, which is difficult to tune precisely. In this work, we investigate the formation of the coreshell nanowires at the atomic scale, by reactive molecular dynamics simulations. The occurrence of two temperature-dependent oxidation mechanisms of ultra-small diameter Si-NWs is demonstrated. We found that control over the Si-core radius and the SiOx (x ≤ 2) oxide shell is possible by tuning the growth temperature and the initial Si-NW diameter. Two different structures were obtained, i.e., ultrathin SiO2 silica nanowires at high temperature and Si core|ultrathin SiO2 silica nanowires at low temperature. The transition temperature is found to linearly decrease with the nanowire curvature. Finally, the interfacial stress is found to be responsible for self-limiting oxidation, depending on both the initial Si-NW radius and the oxide growth temperature. These novel insights allow us to gain control over the exact morphology and structure of the wires, as is needed for their application in nanoelectronics.

Abstract: Modeling the oxidation process of silicon nanowires through reactive force field based molecular dynamics simulations suggests that the formation of Si epoxide defects occurs both at the Si/SiOx interface and at the nanowire surface, whereas for flat surfaces, this defect is experimentally observed to occur only at the interface as a result of stress. In this paper, we argue that the increasing curvature stabilizes the defect at the nanowire surface, as suggested by our density functional theory calculations. The latter can have important consequences for the opto-electronic properties of thin silicon nanowires, since the epoxide induces an electronic state within the band gap. Removing the epoxide defect by hydrogenation is expected to be possible but becomes increasingly difficult with a reduction of the diameter of the nanowires.

Abstract: The thermodynamic aspects of various 2D materials are explored using Density Functional Theory (DFT). Various metal chalcogenides (MX2, M = metal, chalcogen X = S, Se, Te) are investigated with respect to their interaction and stability under different ambient conditions met in the integration process of a transistor device. Their interaction with high- dielectrics is also addressed, in order to assess their possible integration in Complementary Metal Oxide Semiconductor (CMOS) field effect transistors. 2D materials show promise for high performance nanoelectronic devices, but the presence of defects (vacancies, grain boundaries,...) can significantly impact their electronic properties. To assess the impact of defects, their enthalpies of formation and their signature levels in the density of states have been studied. We find, consistently with literature reports, that chalcogen vacancies are the most likely source of defects. It is shown that while pristine 2D materials are in general stable whenever set in contact with different ambient atmospheres, the presence of defective sites affects the electronic properties of the 2D materials to varying degrees. We observe that all the 2D materials studied in the present work show strong reactivity towards radical oxygen plasma treatments while reactivity towards other common gas phase chemical such as O-2 and H2O and groups present at the high- surface varies significantly between species. While energy band-gaps, effective masses and contact resistivities are key criteria in selection of 2D materials for scaled CMOS and tunneling based devices, the phase and ambient stabilities might also play a very important role in the development of reliable nanoelectronic applications.

Abstract: Based on a rational classification of defects in amorphous materials, we propose a simplified model to describe intrinsic defects and hydrogen impurities in amorphous indium gallium zinc oxide (a-IGZO). The proposed approach consists of organizing defects into two categories: point defects, generating structural anomalies such as metal-metal or oxygen-oxygen bonds, and defects emerging from changes in the material stoichiometry, such as vacancies and interstitial atoms. Based on first-principles simulations, it is argued that the defects originating from the second group always act as perfect donors or perfect acceptors. This classification simplifies and rationalizes the nature of defects in amorphous phases. In a-IGZO, the most important point defects are metal-metal bonds (or small metal clusters) and peroxides (O-O single bonds). Electrons are captured by metal-metal bonds and released by the formation of peroxides. The presence of hydrogen can lead to two additional types of defects: metal-hydrogen defects, acting as acceptors, and oxygen-hydrogen defects, acting as donors. The impact of these defects is linked to different instabilities observed in a-IGZO. Specifically, the diffusion of hydrogen and oxygen is connected to positive-and negative-bias stresses, while negative-bias illumination stress originates from the formation of peroxides.

Abstract: We propose that low-strain Si/O superlattices can be constructed by connecting reconstructed Si{001} surfaces by Si-O-Si bridges. Ab initio calculations show that our models are energetically more favorable than all the models proposed so far. The part of our Si/O superlattice model is experimentally accessible just by oxidizing a Si( 001) substrate. To complete our Si/O superlattice model, we propose a three-step method. We also explore the potential of our Si/O superlattice models for new materials used in future Si electronics. We find that the location of the channel where the carriers travel can be controlled between the interfaces and the Si layers by the insertion of O atoms into the Si-Si dimers. By revealing the origins of the interface electron and hole states, we find that similar interface states should be easily achieved for Si slabs and Si substrates. Interestingly, the interface electrons and holes have small effective masses in the direction parallel to the channel and large effective masses in the direction normal to the channel, which makes the Si/O superlattices attractive to be used for channel materials. We also find that the valley splitting of Si is enhanced by the formation of the Si/O/Si interfaces, which is ideal for developing Si-based qubits. Our findings open new perspectives to design and control the electronic properties of Si.

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: Amorphous semiconductors are usually characterized by a low charge carrier mobility, essentially related to their lack of long-range order. The development of such material with higher charge carrier mobility is hence challenging. Part of the issue comes from the difficulty encountered by first-principles simulations to evaluate concepts such as the electron effective mass for disordered systems since the absence of periodicity induced by the disorder precludes the use of common concepts derived from condensed matter physics. In this paper, we propose a methodology based on first-principles simulations that partially solves this problem, by quantifying the degree of delocalization of a wave function and of the connectivity between the atomic sites within this electronic state. We validate the robustness of the proposed formalism on crystalline and molecular systems and extend the insights gained to disordered/amorphous InGaZnO4 and Si. We also explore the properties of p-type oxide semiconductor candidates recently reported to have a low effective mass in their crystalline phases [G. Hautier et al., Nat. Commun. 4, 2292 (2013)]. Although in their amorphous phase none of the candidates present a valence band with delocalization properties matching those found in the conduction band of amorphous InGaZnO4, three of the seven analyzed materials show some potential. The most promising candidate, K2Sn2O3, is expected to possess in its amorphous phase a slightly higher hole mobility than the electron mobility in amorphous silicon.