No Arabic abstract
We present a stochastic simulation method designed to study at an atomic resolution the growth kinetics of compounds characterized by the sp3-type bonding symmetry. Formalization and implementation details are discussed for the particular case of the 3C-SiC material. A key feature of our numerical tool is the ability to simulate the evolution of both point-like and extended defects, whereas atom kinetics depend critically on process-related parameters. In particular, the simulations can describe the surface state of the crystal and the generation/evolution of defects as a function of the initial substrate condition and the calibration of the simulation parameters. We demonstrate that quantitative predictions of the microstructural evolution of the studied systems can be readily compared with the structural characterization of actual processed samples.
We demonstrate how first-principles calculations using density-functional theory (DFT) can be applied to gain insight into the molecular processes that rule the physics of materials processing. Specifically, we study the molecular beam epitaxy (MBE) of arsenic compound semiconductors. For homoepitaxy of GaAs on GaAs(001), a growth model is presented that builds on results of DFT calculations for molecular processes on the beta2-reconstructed GaAs(001) surface, including adsorption, desorption, surface diffusion and nucleation. Kinetic Monte Carlo simulations on the basis of the calculated energetics enable us to model MBE growth of GaAs from beams of Ga and As_2 in atomistic detail. The simulations show that island nucleation is controlled by the reaction of As_2 molecules with Ga adatoms on the surface. The analysis reveals that the scaling laws of standard nucleation theory for the island density as a function of growth temperature are not applicable to GaAs epitaxy. We also discuss heteroepitaxy of InAs on GaAs(001), and report first-principles DFT calculations for In diffusion on the strained GaAs substrate. In particular we address the effect of heteroepitaxial strain on the growth kinetics of coherently strained InAs islands. The strain field around an island is found to cause a slowing-down of material transport from the substrate towards the island and thus helps to achieve more homogeneous island sizes.
Puckered honeycomb Sb monolayer, the structural analog of black phosphorene, has been recently successfully grown by means of molecular beam epitaxy. However, little is known to date about the growth mechanism for such puckered honeycomb monolayer. In this study, by using scanning tunneling microscopy in combination with first-principles density functional theory calculations, we unveil that the puckered honeycomb Sb monolayer takes a kinetics-limited two-step growth mode. As the coverage of Sb increases, the Sb atoms firstly form the distorted hexagonal lattice as the half layer, and then the distorted hexagonal half-layer transforms into the puckered honeycomb lattice as the full layer. These results provide the atomic-scale insight in understanding the growth mechanism of puckered honeycomb monolayer, and can be instructive to the direct growth of other monolayers with the same structure.
In this work we report on kinetic Monte-Carlo calculations of resistive switching and the underlying growth dynamics of filaments in an electrochemical metallization device consisting of an Ag/TiO2/Pt sandwich-like thin film system. The developed model is not limited to i) fast time scale dynamics and ii) only one growth and dissolution cycle of metallic filaments. In particular, we present results from the simulation of consecutive cycles. We find that the numerical results are in excellent agreement with experimentally obtained data. Additionally, we observe an unexpected filament growth mode which is in contradiction to the widely acknowledged picture of filament growth, but consistent with recent experimental findings.
Twelve two-dimensional group-IV monochalcogenide monolayers (SiS, SiSe, SiTe, GeS, GeSe, GeTe, SnS, SnSe, SnTe, PbS, PbSe, and PbTe) with a buckled honeycomb atomistic structure--belonging to symmetry group P3m1--and an out-of-plane intrinsic electric polarization are shown to be metastable by three independendent methods. First, we uncover a coordination-preserving structural transformation from the low-buckled honeycomb structure onto the lower-energy Pnm2$_1$ (or Pmmn for PbS, PbSe, and PbTe) phase to estimate {em energy barriers} $E_B$ that must be overcome during such structural transformation. Using the curvature of the local minima and $E_B$ as inputs to Kramers escape formula, large escape times are found, implying the structural metastability of the buckled honeycomb phase (nevertheless, and with the exception of PbS and PbSe, these phases display escape times ranging from 700 years to multiple times the age of the universe, and can be considered stable for practical purposes only in that relative sense). The second demonstration is provided by phonon dispersion relations that include the effect of long-range Coulomb forces and display no negative vibrational modes. The third and final demonstration of structural metastability is furnished by room-temperature {em ab initio} molecular dynamics for selected compounds. The magnitude of the electronic band gap evolves with chemical composition. Different from other binary two-dimensional compounds such as transition metal dichalcogenide monolayers and hexagonal boron nitride monolayers which only develop an in-plane piezoelectric response, the twelve group-IV monochalcogenide monolayers with a buckled honeycomb structure also display out-of-plane piezoelectric properties.
Understanding the microscopic mechanism of chemical vapor deposition (CVD) growth of two-dimensional molybdenum disulfide (2D MoS2) is a fundamental issue towards the function-oriented controlled growth. In this work, we report results on revealing the growth kinetics of 2D MoS2 via capturing the nucleation seed, evolution morphology, edge structure and terminations at the atomic scale during CVD growth using the transmission electron microscopy (TEM) and scanning transmission electron microscopy (STEM) studies. The direct growth of few- and mono-layer MoS2 onto graphene based TEM grids allow us to perform the subsequent TEM characterization without any solution-based transfer. Two forms of seeding centers are observed during characterizations: (i) Mo-oxysulfide (MoOxS2-y) nanoparticles either in multi-shelled fullerene-like structures or in compact nanocrystals for the growth of fewer-layer MoS2; (ii) Mo-S atomic clusters in case of monolayer MoS2. In particular, for the monolayer case, at the early stage growth, the morphology appears in irregular polygon shape comprised with two primary edge terminations: S-Mo Klein edge and Mo zigzag edge, approximately in equal numbers, while as the growth proceeds, the morphology further evolves into near-triangle shape in which Mo zigzag edge predominates. Results from density-functional theory calculations are also consistent with the inferred growth kinetics, and thus supportive to the growth mechanism we proposed. In general, the growth mechanisms found here should also be applicable in other 2D materials, such as MoSe2, WS2 and WSe2 etc.