Abstract: We performed a detailed investigation of the oxidation of the technologically relevant GaAs(001)-beta 2(2x4) surface via density functional calculations. The purpose is to gain insights on the atomistic mechanisms and local bondings that underlie the degradation of the surface properties once exposed to oxygen. The study comprises the adsorption of single O atoms, through the sampling of several adsorption sites, and the subsequent formation of the O adsorbate at increasing coverage by taking into account multiple-atom adsorption. Based on the evaluation of the energetics and the structural properties of the atomistic models generated, the results here reported delineate a consistent picture of the initial stage of the surface oxidation: (i) at low coverage, in the limit of single O insertions, oxygen is incorporated on the surface forming a twofold-bridging Ga-O-As bond; (ii) at increasing coverage, as multiple O atoms are involved, this is accompanied by the formation of a threefold-coordinated bond (with two Ga and one As atoms); (iii) the latter has important implications regarding the electronic properties of the adsorbate since this O bonding may result in the formation of As dangling bonds. Moreover, a clear trend of increased energy gain for the incorporation of neighboring O atoms compared to single O insertions indicates that the formation of oxide clusters is favored over a regime of uniform oxidation. Our findings provide a detailed description of the O bonding and stress the importance of modeling the adsorption of multiple O atoms for an accurate description of the surface oxidation.

Abstract: Recently, core shell silicon nanowires (Si-NWs) have been envisaged to be used for field-effect transistors and photovoltaic applications. In spite of the constant downsizing of such devices, the formation of ultrasmall diameter core shell Si-NWs currently remains entirely unexplored. We report here on the modeling of the formation of such core shell Si-NWs using a dry thermal oxidation of 2 nm diameter (100) Si nanowires at 300 and 1273 K, by means of reactive molecular dynamics simulations using the ReaxFF potential. Two different oxidation mechanisms are discussed, namely a self-limiting process that occurs at low temperature (300 K), resulting in a Si core I ultrathin SiO2 silica shell nanowire, and a complete oxidation process that takes place at a higher temperature (1273 K), resulting in the formation of an ultrathin SiO2 silica nanowire. The oxidation kinetics of both cases and the resulting structures are analyzed in detail. Our results demonstrate that precise control over the Si-core radius of such NWs and the SiOx (x <= 2.0) oxide shell is possible by controlling the growth temperature used during the oxidation process.

Abstract: Using reactive molecular dynamics simulations based on the ReaxFF potential, we studied the growth mechanism of ultrathin silica (SiO2) layers during hyperthermal oxidation as a function of temperature in the range 100-1300 K. Oxidation of Si(100){2 x 1} surfaces by both atomic and molecular oxygen was investigated for hyperthermal impact energies in the range of 1 to 5 eV. Two different growth mechanisms are found, corresponding to a low temperature oxidation and a high temperature one. The transition temperature between these mechanisms is estimated to be about 700 K. Also, the initial step of the Si oxidation process is analyzed in detail. Where possible, we validated our results with experimental and ab initio data, and good agreement was obtained. 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: 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: 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: 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: 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: Strained SiGe heterostructure-on-insulator (0 0 1) and (1 1 0) PMOSFETs are investigated including important aspects like CV characteristics, mobility, and ON current. The simulations are based on the self-consistent solution of 6 × 6 k · p Schrödinger Equation, multi subband Boltzmann Transport Equation and Poisson Equation, and capture size quantization, strain, crystallographic orientation, and SiGe alloy effects on a solid physical basis. The simulation results are validated by comparison with different experimental data sources. The simulation results show that the strained SiGe HOI PMOSFET with (1 1 0) surface orientation has a higher gate capacitance and a much higher mobility and ON current compared to a similar device with the traditional (0 0 1) surface orientation.

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: The structural and vibrational properties of amorphous germanium oxide (a-GeO₂) are investigated using first-principles calculations based on density functional theory. We first generate an a-GeO₂ structure by first-principles molecular dynamics and analyze its structural properties. The vibrational spectra is then calculated within a density-functional approach. Both static and dynamic properties are in good agreement with experimental data. We next generate defects in our structure (oxygen vacancies with several density and charge states) and consider the most stable atomic configurations, focusing on the vibrational features of threefold coordinated O and divalent Ge centers.

Abstract: We derive an analytical model for the electrostatics and the drive current in a silicon nanowire operating in JFET mode. We show that there exists a range of nanowire radii and doping densities for which the nanowire JFET satisfies reasonable device characteristics. For thin nanowires we have developed a self-consistent quantum mechanical model to obtain the electronic structure.

Abstract: A database-driven approach combined with ab initio density functional theory (DFT) simulations is used to identify and simulate alternative ferroelectric materials beyond Hf(Zr)O-2. The database-driven screening method identifies a class of wurtzite ferroelectric materials. DFT simulations of wurtzite magnesium chalcogenides, including MgS, MgSe, and MgTe, show their potential to achieve improved ferroelectric (FE) stability, simple atomistic unit cell structure, and large FE polarization. Strain engineering can effectively modulate the FE switching barrier height for facilitating FE switching. The effect of the piezoelectric property on the FE switching barrier heights is also examined.

Abstract: This paper benchmarks various epitaxial growth schemes based on n-type group-IV materials as viable source/drain candidates for Ge nMOS devices. Si:P grown at low temperature on Ge, gives an active carrier concentration as high as 3.5 x 10(20) cm(-3) and a contact resistivity down to 7.5 x 10(-9) Omega.cm(2). However, Si:P growth is highly defective due to large lattice mismatch between Si and Ge. Within the material stacks assessed, one option for Ge nMOS source/drain stressors would be to stack Si:P, deposited at contact level, on top of a selectively grown n-SiyGe1-x-ySnx at source/drain level, in line with the concept of Si passivation of n-Ge surfaces to achieve low contact resistivities as reported in literature (Martens et al. 2011 Appl. Phys. Lett., 98, 013 504). The saturation in active carrier concentration with increasing P (or As)-doping is the major bottleneck in achieving low contact resistivities for as-grown Ge or SiyGe1-x-ySnx. We focus on understanding various dopant deactivation mechanisms in P-doped Ge and Ge1-xSnx alloys. First principles simulation results suggest that P deactivation in Ge and Ge1-xSnx can be explained both by P-clustering and donor-vacancy complexes. Positron annihilation spectroscopy analysis, suggests that dopant deactivation in P-doped Ge and Ge1-xSnx is primarily due to the formation of P-n-V and SnmPn-V clusters. (C) 2020 The Author(s). Published on behalf of The Electrochemical Society by IOP Publishing Limited.

Abstract: Density functional theory simulations are used to identify the structural factors that define the material properties of ovonic threshold switches (OTS). They show that the nature of mobility-gap trap states in amorphous Ge-rich Ge50Se50 is related to Ge-Ge bonds, whereas in Se-rich Ge30Se70 the Ge valence-alternating-pairs and Se lone-pairs dominate. To obtain a faithful description of the electronic structure and delocalization of states, it is required to combine hybrid exchange-correlation functionals with large unit-cell models. The extent of localization of electronic states depends on the applied external electric field. Hence, OTS materials undergo structural changes during electrical cycling of the device, with a decrease in the population of less exothermic Ge-Ge bonds in favor of more exothermic Ge-Se. This reduces the amount of charge traps, which translates into coordination changes, an increase in mobility-gap, and subsequently changes in the selector-device electrical parameters. The threshold voltage drift process can be explained by natural evolution of the nonpreferred Ge-Ge bonds (or “chains”/clusters thereof) in Ge-rich GexSe1-x. The effect of extrinsic doping is shown for Si and N, which introduce strong covalent bonds into the system, increase both mobility-gap and crystallization temperature, and decrease the leakage current.