skip to main content


Title: Modeling for Structural Engineering and Synthesis of Two Dimensional WSe2 Using a Newly Developed ReaxFF Reactive Force Field
Atomistic simulation techniques have become an indispensable tool to acquire a fundamental understanding of growth and structural characteristics of two-dimensional (2D) materials of interest, thereby accelerating experimental research in the same field. A new ReaxFF reactive force field presented here is the first comprehensive empirical potential that is explicitly designed to capture the most prominent features of 2D WSe2 solid-phase chemistry, such as defect formation as a function of local geometry and chalcogen chemical potential, vacancy migration and phase transition, thus enabling cost-effective and reliable characterization of 2D WSe2 at large length scales and time scales much longer than what is accessible by first-principles theory. This potential, validated using extensive first-principles energetics data on both periodic and nonperiodic systems and experimental measurements, can accurately describe the mechanochemical coupling between monolayer deformations and vacancy energetics, providing valuable atomistic insights into the morphological evolution of a monolayer in different environments in terms of loading conditions and various concentrations and distributions of defects. Since understanding how growth is affected by the local chemical environment is vital to fabricating efficient and functional atomically thin 2D WSe2, the new ReaxFF description enables investigations of edge-controlled growth of single crystals of 2D WSe2 using reactive environments closely matching experimental conditions at a low computational cost.  more » « less
Award ID(s):
1654107
NSF-PAR ID:
10231403
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Journal of physical chemistry
Volume:
124
Issue:
51
ISSN:
1932-7447
Page Range / eLocation ID:
28285–28297
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The successful discovery and isolation of graphene in 2004, and the subsequent synthesis of layered semiconductors and heterostructures beyond graphene have led to the exploding field of two-dimensional (2D) materials that explore their growth, new atomic-scale physics, and potential device applications. This review aims to provide an overview of theoretical, computational, and machine learning methods and tools at multiple length and time scales, and discuss how they can be utilized to assist/guide the design and synthesis of 2D materials beyond graphene. We focus on three methods at different length and time scales as follows: (i) nanoscale atomistic simulations including density functional theory (DFT) calculations and molecular dynamics simulations employing empirical and reactive interatomic potentials; (ii) mesoscale methods such as phase-field method; and (iii) macroscale continuum approaches by coupling thermal and chemical transport equations. We discuss how machine learning can be combined with computation and experiments to understand the correlations between structures and properties of 2D materials, and to guide the discovery of new 2D materials. We will also provide an outlook for the applications of computational approaches to 2D materials synthesis and growth in general.

     
    more » « less
  2. Abstract

    Reproducible wafer-scale growth of two-dimensional (2D) materials using the Chemical Vapor Deposition (CVD) process with precise control over their properties is challenging due to a lack of understanding of the growth mechanisms spanning over several length scales and sensitivity of the synthesis to subtle changes in growth conditions. A multiscale computational framework coupling Computational Fluid Dynamics (CFD), Phase-Field (PF), and reactive Molecular Dynamics (MD) was developed – called the CPM model – and experimentally verified. Correlation between theoretical predictions and thorough experimental measurements for a Metal-Organic CVD (MOCVD)-grown WSe2model material revealed the full power of this computational approach. Large-area uniform 2D materials are synthesized via MOCVD, guided by computational analyses. The developed computational framework provides the foundation for guiding the synthesis of wafer-scale 2D materials with precise control over the coverage, morphology, and properties, a critical capability for fabricating electronic, optoelectronic, and quantum computing devices.

     
    more » « less
  3. null (Ed.)
    Defects have a profound impact on the electronic and physical properties of crystals. For two-dimensional (2D) materials, many intrinsic point defects have been reported, but much remains to be understood about their origin. Using scanning transmission electron microscopy imaging, this study discovers various linear arrays of W-vacancy defects that are explained in the context of the crystal growth of coalesced, monolayer WS2. Atomistic-scale simulations show that vacancy arrays can result from steric hindrance of bulky gas-phase precursors at narrowly separated growth edges, and that increasing the edge separation leads to various intact and defective growth modes, which are driven by competition between the catalytic effects of the sapphire substrate and neighboring growth edge. Therefore, we hypothesize that the arrays result from combined growth modes, which directly result from film coalescence. The connections drawn here will guide future synthetic and processing strategies to harness the engineering potential of defects in 2D monolayers. 
    more » « less
  4. null (Ed.)
    Molecular insights into graphene-catalyst surface interactions can provide useful information for the efficient design of copper current collectors with graphitic anode interfaces. As graphene bending can affect the local electron density, it should reflect its local reactivity as well. Using ReaxFF reactive molecular simulations, we have investigated the possible bending of graphene in vacuum and near copper surfaces. We describe the energy cost for graphene bending and the binding energy with hydrogen and copper with two different ReaxFF parameter sets, demonstrating the relevance of using the more recently developed ReaxFF parameter sets for graphene properties. Moreover, the draping angle at copper step edges obtained from our atomistic simulations is in good agreement with the draping angle determined from experimental measurements, thus validating the ReaxFF results. 
    more » « less
  5. Type II germanium clathrates have recently been investigated for potential applications as anodes in batteries due to their cage-like structures that can accommodate electrochemical insertion of guest ions. To synthesize type II Ge clathrates (Ge136), several experimental routes use thermal or electrochemical desodiation of the Zintl phase compound Na4Ge4. However, the mechanism by which Na atoms are removed from the precursor to form clathrates is not well understood. Herein, we use first-principles density functional theory and nudged elastic band calculations to understand the reaction mechanism and formation energies of the products typically observed in the synthesis, namely, NaδGe136 (0 < δ < 24) type II clathrates and hexagonal phase Na1–xGe3+z. Specifically, we confirm the energetic feasibility of Na vacancy formation in Na4Ge4 and find that the barrier for Na vacancy migration is only 0.37 eV. This relatively low energy barrier is consistent with the ease with which Na4Ge4 can be desodiated to form the products. We also discuss the energetics, sodium migration pathways, and potential electrochemical performance of Ge136 as anode material for Na-ion batteries. Overall, this study highlights how first-principles calculations can be used to understand the synthesis mechanism and desodiation processes in clathrate materials and will help guide researchers in the design and evaluation of new open framework compounds as viable materials for energy storage applications. 
    more » « less