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 structural and electronic properties of silicene and germanene on metallic and non-metallic substrates are investigated theoretically, using first-principles simulations. We first study the interaction of silicene with Ag(111) surfaces, focusing on the (4x4) silicene/Ag structure. Due to symmetry breaking in the silicene layer (nonequivalent number of top and bottom Si atoms), silicene is predicted to be semiconducting, with a computed energy gap of about 0.3 eV. However, the charge transfer occurring at the silicene/Ag(111) interface leads to an overall metallic system. We next investigate the interaction of silicene and germanene with hexagonal non-metallic substrates, namely ZnS and ZnSe. On reconstructed (semiconducting) (0001) ZnS or ZnSe surfaces, silicene and germanene are found to be semiconducting. Remarkably, the nature (indirect or direct) and magnitude of their energy band gap can be controlled by an out-of-plane electric field.

Abstract: First-principle complex band structures have been computed for rutile TiO2 and tetragonal ZrO2 insulating materials that are of current technological relevance to dynamic random accessmemorymetal-insulator-metal (MIM) capacitors. From the magnitude of the complex wave vectors in different orientations, the most penetrating orientations have been identified. Tunneling effective masses m(tunnel) have been extracted, are shown to be a crucial parameter for the intrinsic leakage, and are identified to be an important parameter in further scaling of MIM capacitors.

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 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: 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: 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: 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: Transition metal oxide-based resistor random access memory (RRAM) takes advantage of oxygen-related defects in its principle of operation. Since the change in resistivity of the material is controlled by the oxygen deficiency level, it is of major importance to quantify the kinetics of the oxygen diffusion, key factor for oxide stoichiometry. Ab initio accelerated molecular dynamics techniques are employed to investigate the oxygen diffusivity in amorphous hafnia (HfOx, x = 1.97, 1.0, 0.5). The computed kinetics is in agreement with experimental measurements. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.3697690]

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: 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: 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: Germanium combined with high-κ dielectrics is investigated for the next generations of CMOS devices. Therefore, we study reaction mechanisms for Al2O3 atomic layer deposition on sulfur passivated Ge using calculations based on density functional theory and total reflection X-ray fluorescence (TXRF). TXRF indicates 6 S/nm2 and 4 Al/nm2 after the first TMA/H2O reaction cycle, and growth inhibition from the second reaction cycle on. Calculations are performed on molecular clusters representing −GeSH surface sites. The calculations confirm that the TMA reaction does not affect the S content. On fully SH-terminated Ge, TMA favorably reacts with up to three −GeSH sites, resulting in a near tetrahedral Al coordination. Electron deficient structures with a GeS site shared between two Al atoms are proposed. The impact of the cluster size on the structures and reaction energetics is systematically investigated.

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.