Modeling the vertical growth of van der Waals stacked 2D materials using the diffuse domain method

TL;DR

This paper develops a multiscale model using the diffuse domain method to predict and control the vertical growth of 2D materials, providing insights into growth conditions and layer morphology.

Contribution

It introduces a novel phase-field/diffuse domain approach to model 2D material growth, extending previous models to arbitrarily-shaped layers and enabling better control of multilayer fabrication.

Findings

Vertical growth conditions follow an analytic thermodynamic criterion.

Layer boundaries can develop significant curvature during growth.

The model offers a mechanistic framework for optimizing 2D material growth.

Abstract

Vertically-stacked monolayers of graphene and other atomically-thin 2D materials have attracted considerable research interest because of their potential in fabricating materials with specifically-designed properties. Chemical vapor deposition has proved to be an efficient and scalable fabrication method. However, a lack of mechanistic understanding has hampered efforts to control the fabrication process beyond empirical trial-and-error approaches. In this paper, we develop a general multiscale Burton-Cabrera-Frank (BCF) type model of the vertical growth of 2D materials to predict the necessary growth conditions for vertical versus in-plane (monolayer) growth of arbitrarily-shaped layers. This extends previous work where we developed such a model assuming the layers were fully-faceted (Ye et al., ACS Nano, 11, 12780-12788, 2017). To solve the model numerically, we reformulate the system using the phase-field/diffuse domain method that enables the equations to be solved in a fixed regular domain. We use a second-order accurate, adaptive finite-difference/nonlinear multigrid algorithm to discretize and solve the discrete system. We investigate the effect of parameters, including the van der Waals interaction energies between the layers, the kinetic attachment rates, the edge-energies and the deposition flux, on layer growth and morphologies. While the conditions that favor vertical growth generally follow an analytic thermodynamic criterion we derived for circular layers, the layer boundaries may develop significant curvature during growth, consistent with experimental observations. Our approach provides a mechanistic framework for controlling and optimizing the growth multilayered 2D materials.