Abstract: The interaction of silicene, the silicon counterpart of graphene, with (0001) ZnS surfaces is investigated theoretically, using first-principles simulations. The charge transfer occurring at the silicene/(0001) ZnS interface leads to the opening of an indirect energy band gap of about 0.7 eV in silicene. Remarkably, the nature (indirect or direct) and magnitude of the energy band gap of silicene can be controlled by an external electric field: the energy gap is predicted to become direct for electric fields larger than about 0.5 V angstrom(-1), and the direct energy gap decreases approximately linearly with the applied electric field. The predicted electric field tunable energy band gap of the silicene/(0001) ZnS interface is very promising for its potential use in nanoelectronic devices.

Abstract: Spectacular progress in developing advanced Si circuits with reduced size, along the track of Moore's law, has been relying on necessary developments in wet cleaning of nanopatterned Si wafers to provide contaminant free surfaces. The most efficient cleaning is achieved when complete wetting can be realized. In this work, ordered arrays of silicon nanopillars on a hitherto unexplored small scale have been used to study the wetting behavior on nanomodulated surfaces in a substantial range of surface treatments and geometrical parameters. With the use of optical reflectance measurements, the nanoscale water imbibition depths have been measured and the transition to the superhydrophobic Cassie-Baxter state has been accurately determined. For pillars of high aspect ratio (about 15), the transition occurs even when the surface is grafted with a hydrophilic functional group. We have found a striking consistent deviation between the contact angle measurements and the straightforward application of the classical wetting models. Molecular dynamics simulations show that these deviations can be attributed to the long overlooked atomic-scale surface perturbations that are introduced during the nanofabrication process. When the transition condition is approached, transient states of partial imbibition that characterize intermediate states between the Wenzel and Cassie-Baxter states are revealed in our experiments.

Abstract: Because of its excellent properties in nanotechnology applications, atomic layer deposition of ruthenium (Ru) has been the subject of numerous experimental studies. Recently, two different Ru precursors were compared for plasma-enhanced atomic layer deposition (PEALD) of Ru, and their reactivity was found to be different. Inhibition was observed for bis(ethylcyclopentadienyl)ruthenium (Ru(EtCp)(2)), while nearly linear growth behavior was observed for (methylcyclopentadienyl-pyrrolyl)ruthenium (Ru(MeCp)Py). To understand this difference in reactivity, we investigate the adsorption of RuCp, and RuCpPy (i.e., without substituents) on a TiN surface using calculations based on periodic boundary conditions density functional theory (DFT) combined with experiments based on Rutherford backscattering spectroscopy (RBS). The calculations demonstrate that the RuCpPy precursor chemisorbs on the TiN(100) surface while the RuCp2 precursor only physisorbs. We propose a reaction mechanism for the chemisorption of RuCpPy. The area density of the calculated RuCpPy surface species is compared with the experimental values from RBS. The impact of a H-plasma is also investigated. The DFT calculations and experimental results from RBS provide insight into the adsorption processes of the RuCpPy and RuCp2 precursors on the TiN(100) surface.

Abstract: The exact structure and properties of the Si vertical bar SiO2 interface are very important in microelectronics and photovoltaic devices such as metal-oxide-semiconductor field-effect transistors (MOSFETs) and solar cells. Whereas Si vertical bar SiO2 structures are traditionally produced by thermal oxidation, hyperthermal oxidation shows a number of promising advantages. However, the Si vertical bar SiO2 interface induced in hyperthermal Si oxidation has not been properly investigated yet. Therefore, in this work, the interface morphology and interfacial stresses during hyperthermal oxidation at room temperature are studied using reactive molecular dynamics simulations based on the ReaxFF potential. Interface thickness and roughness, as well as the bond length and bond angle distributions in the interface are discussed and compared with other models developed for the interfaces induced by traditional thermal oxidation. The formation of a compressive stress is observed. This compressive stress, which at the interface amounts about 2 GPa, significantly slows down the inward silica growth. This value is close to the experimental value in the Si vertical bar SiO2 interface obtained in traditional thermal oxidation.

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.