Abstract: We investigate at the atomic level the most probable phase transformations under strain, that are responsible for the ferroelectric/ antiferroelectric behavior in Hf1-xZrxO2 materials. Four different crystalline phase transformations exhibit a polar/non-polar transition: monoclinic-to-orthorhombic requires a gliding strain tensor, orthorhombic-to-orthorhombic transformation does not need strain to polarize the material, whereas tetragonal-to-cubic cell compression and tetragonal-to-orthorhombic cell elongation destabilizes the non-polar tetragonal phase, facilitating the transition towards a polar atomic configuration, therefore changing the polarization-electric field loop from antiferroelectric to ferroelectric. Oxygen vacancies can reduce drastically the polarization reversal barriers.

Abstract: Si:C and Si:C:P alloys are potential candidates for source-drain stressor applications in n-type Fin Field Effect Transistors (FinFETs). Increasing the C content to achieve high strain results in the arrangement of C atoms as third nearest neighbors (3nn) in the Si: C lattice. During thermal annealing, the presence of C atoms as 3nn may promote clustering at the interstitial sites, causing loss of stress. The concentration of C atoms as 3nn is reduced by the incorporation of a small amount of Ge atoms during the growth, whereas in-situ P doping does not influence this 3nn distribution [J Solid State Sci. Technol vol 6, p 755, 2017]. Small amounts of Ge are provided during low temperature selective epitaxial growth scheme, which are based on cyclic deposition and etching (CDE). In this work, we aim to provide physical insights into the aforementioned phenomena, to understand the behavior of 3nn C atoms and the types of defects that are formed in the annealed Si: C films. Using ab-initio simulations, the Ge-C interaction in the Si matrix is investigated and this insight is used to explain how the Ge incorporation leads to a reduced 3nn distribution of the C atoms. The interaction between C and P in the Si: C: P films is also investigated to explain why the P incorporation has not led to a reduction in the 3nn distribution. We then report on the Raman characterization of Si: C layers subjected to post epi annealing. As the penetration depth of the laser is dependent on the wavelength, Raman measurements at two different wavelengths enable us to probe the depth distribution of 3nn C atoms after applying different annealing conditions. We observed a homogeneous loss in 3nn C throughout the layer. Whereas in the kinematic modeling of high resolution X-ray diffraction spectra, a gradient in the substitutional C loss was observed close to the epitaxial layer/substrate interface. This gradient can be due to the out diffusion of C atoms into the Si substrate or to the formation of interstitial C clusters, which cannot be distinguished in HR-XRD. Deep Level Transient Spectroscopy indicated that the prominent out-diffusing species was interstitial CO complex while the interstitial C defects were also prevalent in the epi layer. (c) 2019 The Electrochemical Society.

Abstract: The scalability potential of low resistivity ternary metallic alloys (MAX) as an interconnect medium has been benchmarked against copper through first-principle simulations. We report that some carbon and nitrogen MAX phases have the potential to display a reduced sensitivity of their intrinsic resistivity to scaling, while showing improved electromigration properties.

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: 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.

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: 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: 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: 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: 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: 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: 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: 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: 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.