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: 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: 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: We report on the thin film resistivity of several platinum-group metals (Ru, Pd, Ir, and Pt). Platinum-group thin films show comparable or lower resistivities than Cu for film thicknesses below about 5 nm due to a weaker thickness dependence of the resistivity. Based on experimentally determined mean linear distances between grain boundaries as well as ab initio calculations of the electron mean free path, the data for Ru, Ir, and Cu were modeled within the semiclassical Mayadas-Shatzkes model [Phys. Rev. B 1, 1382 (1970)] to assess the combined contributions of surface and grain boundary scattering to the resistivity. For Ru, the modeling results indicated that surface scattering was strongly dependent on the surrounding material with nearly specular scattering at interfaces with SiO2 or air but with diffuse scattering at interfaces with TaN. The dependence of the thin film resistivity on the mean free path is also discussed within the Mayadas-Shatzkes model in consideration of the experimental findings. Published by AIP Publishing.

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: <script type='text/javascript'>document.write(unpmarked('The crystal and defect structure of SnS crystals grown using chemical vapour deposition for application in electronic devices are investigated. The structural analysis shows the presence of two distinct crystal morphologies, that is thin flakes with lateral sizes up to 50 m and nanometer scale thickness, and much thicker but smaller crystallites. Both show similar Raman response associated with SnS. The structural analysis with transmission electron microscopy shows that the flakes are single crystals of -SnS with [010] normal to the substrate. Parallel with the surface of the flakes, lamellae with varying thickness of a new SnS phase are observed. High-resolution transmission electron microscopy (TEM), scanning transmission electron microscopy (STEM), first-principles simulations (DFT) and nanobeam diffraction (NBD) techniques are employed to characterise this phase in detail. DFT results suggest that the phase is a strain stabilised \u0027 one grown epitaxially on the -SnS crystals. TEM analysis shows that the crystallites are also -SnS with generally the [010] direction orthogonal to the substrate. Contrary to the flakes the crystallites consist of two to four grains which are tilted up to 15 degrees relative to the substrate. The various grain boundary structures and twin relations are discussed. Under high-dose electron irradiation, the SnS structure is reduced and -Sn formed. It is shown that this damage only occurs for SnS in direct contact with SiO2. Lay description SnS is a p-type semiconductor, which has attracted significant interest for electronic devices due to its unique properties, low-toxicity and abundance of Sn in nature. Although in the past it has been most extensively studied as the absorber material in solar cells, it has recently garnered interest for application as a p-type two-dimensional semiconductor in nanoelectronic devices due to its anisotropic layered structure similar to the better known phosphorene. Tin sulphide can take the form of several phases and the electronic properties of the material depend strongly on its crystal structure. It is therefore crucial to study the crystal structure of the material in order to predict the electronic properties and gain insight into the growth mechanism. In this work, SnS crystals deposited using a chemical vapour deposition technique are investigated extensively for their crystal and defect structure using transmission electron microscopy (TEM) and related techniques. We find the presence of two distinct crystal morphologies, that is thin flakes with lateral sizes up to 50 m and nm scale thickness, and much thicker but smaller crystallites. The flakes are single crystals of -SnS and contain lamellae with varying thickness of a different phase which appear to be -SnS at first glance. High-resolution scanning transmission electron microscopy is used to characterise these lamellae where the annular bright field (ABF) mode better reveals the position of the sulphur columns. The sulphur columns in the lamellae are found to be shifted relative to the -SnS structure which indicates the formation of a new phase which is a distorted version of the phase which we tentatively refer to as \u0027-SnS. Simulations based on density functional theory (DFT) are used to model the interface and a similar shift of sulphur columns in the -SnS layer is observed which takes place as a result of strong interaction at the interface between the two phases resulting in strain transfer. Nanobeam electron diffraction (NBD) is used to map the lattice mismatch in the thickness of the flakes which reveals good in-plane matching and some expansion out-of-plane in the lamellae. Contrary to the flakes the crystallites are made solely of -SnS and consist of two to four grains which are tilted up to 15 degrees relative to the substrate. The various grain boundary structures and twin relations are discussed. At high electron doses, SnS is reduced to -Sn, however the damage occurs only for SnS in direct contact with SiO2.'));

Abstract: The electronic properties of two-dimensional honeycomb structures of molybdenum disulfide (MoS(2)) subjected to biaxial strain have been investigated using first-principles calculations based on density functional theory. On applying compressive or tensile bi-axial strain on bi-layer and mono-layer MoS(2), the electronic properties are predicted to change from semiconducting to metallic. These changes present very interesting possibilities for engineering the electronic properties of two-dimensional structures of MoS(2).

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: The idea of stacking multiple monolayers of different two-dimensional materials has become a global pursuit. In this work, a silicene armchair nanoribbon of width W and van der Waals-bonded to different transition-metal dichalcogenide (TMD) bilayer substrates MoX2 and WX2, where X = S, Se, Te is considered. The orbital resolved electronic structure and ballistic transport properties of these systems are simulated by employing van der Waals-corrected density functional theory and nonequilibrium Green's functions. We find that the lattice mismatch with the underlying substrate determines the electronic structure, correlated with the silicene buckling distortion and ultimately with the contact resistance of the two-terminal system. The smallest lattice mismatch, obtained with the MoTe2 substrate, results in the silicene ribbon properties coming close to those of a freestanding one. With the TMD bilayer acting as a dielectric layer, the electronic structure is tunable from a direct to an indirect semiconducting layer, and subsequently to a metallic electronic dispersion layer, with a moderate applied perpendicular electric field.

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: It is generally understood that the resistivity of metal thin films scales with film thickness mainly due to grain boundary and boundary surface scattering. Recently, several experiments and ab initio simulations have demonstrated the impact of crystal orientation on resistivity scaling. The crystal orientation cannot be captured by the commonly used resistivity scaling models and a qualitative understanding of its impact is currently lacking. In this work, we derive a resistivity scaling model that captures grain boundary and boundary surface scattering as well as the anisotropy of the band structure. The model is applied to Cu and Ru thin films, whose conduction bands are (quasi-) isotropic and anisotropic, respectively. After calibrating the anisotropy with ab initio simulations, the resistivity scaling models are compared to experimental resistivity data and a renormalization of the fitted grain boundary reflection coefficient can be identified for textured Ru.

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

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.