Critical Role of Precursor Flux in Modulating Nucleation Density in 2D Material Synthesis Revealed by Digital Twin

doi:10.21203/rs.3.rs-4957545/v1

Download PDF

Article

Critical Role of Precursor Flux in Modulating Nucleation Density in 2D Material Synthesis Revealed by Digital Twin

https://doi.org/10.21203/rs.3.rs-4957545/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted

You are reading this latest preprint version

Chemical vapor deposition (CVD) is the most widespread approach for two-dimensional (2D) material synthesis, yet control of nucleation density remains a major hurdle towards large-area growth. We find that precursor flux, a function of gas velocity and precursor concentration, is the critical parameter controlling nucleation. We observe that for a vertically aligned substrate, the presence of a cavity/slot in the substrate-supporting plate creates an enhanced growth zone for 2D-MoS₂. The effect of this confined space on nucleation density is experimentally verified by electron microscopy. To understand this intriguing observation, we developed a hyper-realistic multiphysics computational fluid dynamics model, i.e., a digital-twin of our CVD reactor, which reveals that space confinement achieves nearly-zero gas velocities. Digital twin-informed calculations indicate significantly lower metal precursor flux at the confined space during initial stages of growth, while precursor concentration is uniform across the substrate. The digital twin also makes an important prediction regarding a large time-lag between the set temperature, reactor environmental temperature, and substrate temperature, with implications for nucleation and growth. We offer a framework for designing confined spaces to control nucleation via regulating precursor flux, and for simulating reactor parameters for rapid optimization via the digital-twin model.

Physical sciences/Materials science/Nanoscale materials/Two-dimensional materials

Physical sciences/Materials science/Theory and computation

2D materials

chemical vapor deposition

Computational Fluid Dynamics (CFD)

precursor flux

nucleation density

digital-twin.

Transition metal dichalcogenides (TMDs) are exemplary two-dimensional (2D) materials, promising for electronic and optoelectronic applications^1–3, as well as for catalysis, including carbon dioxide reduction⁴ and hydrogen evolution reactions⁵. Mechanical exfoliation⁶ and chemical vapor deposition (CVD)^7–9 are the most common methods used for the synthesis of TMD monolayers (MLs). CVD can potentially provide high uniformity, repeatability, control of layer number, and scalability to large areas. However, CVD is a complex multiscale process involving heating precursors (gaseous or solid), vapor-phase transport by an inert carrier gas to the substrate, and subsequent reaction and deposition. Obtaining high-quality MLs using CVD is challenging due to reactor-to-reactor variability, and local changes in temperature, concentration, and flow fields^10,11. Even a seemingly minor variation in these experimental parameters can significantly impact the quality of the growth^12–14. However, multiphysics simulation studies that model the complex transport phenomena occurring in a CVD reactor, and use them to develop new growth strategies are lacking.

There is a critical need to connect experimental parameters, such as the set temperature and gas flow rate, to mesoscopic variables such as the local temperature, concentration, and fluid flow profiles, and finally to thermodynamic variables and experimental observations of growth. For example, temperature impacts evaporation rates of the precursors^15,16, growth rate, and shape of the synthesized 2D material^17,18. However, the temperature profile inside a CVD reactor is expected to be non-uniform, with implications on non-stoichiometric growth. Similarly, obtaining a quantitative model of the vapor pressure (concentration) profile is critical for modifying nucleation and growth kinetics. Further, the velocity profile of the carrier gas influences concentration profile at the substrate^19,20, as well as sticking coefficients²¹. Therefore, developing a digital-twin (in-silico CVD reactor) for a real CVD system will lead to an informed modification of experimental parameters and improve reproducibility, thus going beyond the currently practiced experimental trial-and-error paradigm. Previous simulation studies have mostly aimed at understanding the reaction pathways controlling growth via density functional theory^22,23 (DFT), though some recent studies focus on computational fluid dynamics (CFD) aspects^24–26. However, the effect of the complete reactor geometry in modulating gas flow dynamics, along with precursor concentration and temperature model, in a multiple temperature zone setup is still critically missing.

Importantly, the combined role of local velocity and concentration, i.e. a treatment of precursor flux is critically missing. For example, confined space geometries have been utilized in CVD processes to modulate precursor concentration and velocity profiles^27–30. Gas flow velocity and temperature field simulations have been performed for the case of single heating zone reactors³¹, and limited to horizontally positioned substrates^32,33. However, velocity, concentration, and temperature profiles have been treated separately, leading to an incomplete understanding of CVD process, including nucleation and growth.

An experimental knob to directly control nucleation density is needed for extending to large-area growth of 2D materials. Random, uninhibited nucleation results in a high nucleation density, limiting the size of obtained flakes^34,35. A higher precursor concentration at the substrate initiates a “self-seeding” nucleation mechanism, which may also result in multi-layer growth³⁶. Thus, a controlled nucleation density is the key to achieving desired ML growth and flake size³⁴.

Here, we experimentally perform oxygen-assisted three-zone atmospheric pressure CVD (AP-CVD) of ML MoS₂ and create the corresponding digital-twin: a comprehensive computational fluid dynamics (CFD) model including heat, mass, and momentum transport equations. For a vertically aligned substrate (normal to the incoming gas flow), we find that the presence of a slot (cavity) in the supporting plate creates an excellent growth zone for MoS₂, compared to no growth observed on the rest of the substrate, and on the substrate which does not have a slot. We create a COMSOL Multiphysics© finite element model by accounting for fluid flow dynamics, evaporation of chalcogen and metal oxide precursors, realistic substrate geometries, and temperature profiles. Using this hyper-realistic model, we observed no significant concentration gradients for the vertically kept substrate. However, a velocity dead-zone is created in the slot area, wherein we find direct evidence of reduced nucleation density via electron microscopy. Thus, precursor flux, which is a function of gas velocity, precursor concentration and concentration gradient, emerges as the critical parameter controlling nucleation density. Further, informed by our digital-twin, we observed a large time lag between the set temperature, reactor temperature, and substrate temperature, with consequences on growth uniformity and nucleation. Thus, using this model, we make several critical observations on nucleation density, spatial temperature profiles, and gas flow. The deep study of precursor flux, enabled by the digital-twin, will lead to the synthesis of higher quality 2D materials and faster optimization of the nucleation/growth process.

The setup used for ML MoS₂ growth (on SiO₂/Si substrates) is shown in Fig. 1a. The reactor is a quartz tube (60 mm diameter). Three heated zones provide independent control over evaporation temperature of the precursors and growth temperature at the substrate, with the zones separated by 75 mm wide insulation barriers. Inert carrier gas (Ar or N₂) transports the vaporized precursors to the substrate kept downstream. We also introduce a small amount of oxygen (O₂) into the system (see Methods for details of the growth process), as it is shown to improve the growth³⁷. Before heating, the reactor is flushed with high-purity inert gas to obtain reduced levels of residual moisture and O₂ (see Supporting Information section I and Video V1 for corresponding simulations of the tube environment). Figure 1b shows a model of the reactor, as used for the simulations with appropriate boundary conditions and sources. Chalcogen (sulfur) powder and transition metal oxide (MoO₃) powders are placed in alumina boats in the first and second heating zones respectively.

To ensure maximum precursor impingement on the growing surface, the substrate is placed vertically in the third heating zone. As we will see later, vertical placement also enables minimum spatial variation in the precursor concentration, providing a more uniform growth environment. In contrast, for a horizontally kept substrate, even for small dimensions (1 × 1 cm²), a concentration gradient is expected^32,33 (also see Supporting Information section II). This gradient is most pronounced during the nucleation stage, and critically affects the type of growth obtained. To support the substrate in a vertical configuration, we created a slot (a groove, 1 mm deep, 1–2 mm wide) in the bottom supporting alumina plate, and supported at the back by an alumina boat (see Fig. 1b, lower left).

Remarkably, presence of the slot is decisive toward obtaining ML MoS₂ growth at similar growth conditions across different growth runs and substrates. Figure 1c and 1d are optical microscopy images of 300 nm SiO₂/Si substrates after CVD growth, when the substrate is placed inside the slot, and without the slot, respectively. A Raman peak shift of ~ 20 cm^− 1, and typical optical contrast confirms the presence of ML MoS₂³⁸ (see Supporting Information section III for optical spectroscopy and SEM characterization). The presence of ML MoS₂ flakes only in the case of the slotted configuration is striking, and points toward the significant effect of the slot on modifying local growth environment.

Development of the multiphysics digital-twin. We used COMSOL Multiphysics© to develop a coupled gas-flow, precursor concentration, and temperature distribution 3D model for a realistic simulation of the CVD process. To understand the gas flow dynamics and the precursor concentration inside the slot, we carefully modeled the source geometry, substrate supports, and specific materials of various assemblies. The synthesis process can be divided into four broad parts – nitrogen flushing, reactor heating, precursor vaporization, and reaction and deposition. We developed an integrated model for the first three of the mentioned processes, simulating the actual experimental conditions as closely as possible (see Methods for mesh sizes for simulations and specific COMSOL modules used). The deposition process is modeled considering evaporation and downstream transport of powder precursors to the substrate.

In the context of CVD, transfer of physical quantities inside a physical system is a combination of convection (due to inert carrier gas flow) and diffusion^25,26 (due to concentration gradients). The governing equation for describing such a physical phenomenon is the diffusion-convection equation (DCE):

$$\:\begin{array}{c}\frac{\partial\:C\left(x,y,z,t\right)}{\partial\:t}=\:\overrightarrow{\nabla\:}.\left(D\overrightarrow{\nabla\:}C\left(x,y,z,t\right)-\overrightarrow{v}C\left(x,y,z,t\right)\right)\#\left(1\right)\end{array}$$

where C is the local fluid concentration (mol/m³), $\:\overrightarrow{v}$ is the carrier gas velocity (m/s), D is the precursor diffusivity in the carrier gas (m²/s)³⁹, and x (m) is the coordinate along the reactor. Using the DCE, we can also model precursor evaporation with the addition of a source term:

where r_v is the precursor source term, i.e., molar production rate per unit volume (mol/(m³.s)). The saturation vapor pressure (P_s) for a given species is primarily a function of temperature (T)^40,41. We used the Knudsen-Langmuir Eq. 4² to model the temperature dependent evaporation of the precursors, as detailed in Supporting Information section IV.

Understanding concentration profiles and fluid flow inside the reactor. We first examine precursor and O₂ spatiotemporal profiles at the substrate. We note that our synthesis process minimizes precursor evaporation prior to the set temperature being realized (see Supporting Information section IV for corresponding axial simulations of MoO₃ and O₂). Hence, for the purpose of modeling precursor evaporation, we allow the precursors to evaporate only after heating is completed (30 mins).

Figure 2 shows the MoO₃/S₈ ratio and O₂ profiles at the substrate after 5 minutes of the beginning of growth, with and without the slot. Elemental sulfur (S) in its vaporized form exists as S₈ (octa-sulfur) molecules⁴³. Hence, the concentrations of sulfur and oxygen depicted are those of S₈ and O₂ species respectively, and are discussed as such in the rest of the paper. Elemental concentrations may be obtained by multiplying S₈ and O₂ concentrations by a factor of 8 and 2 respectively. Importantly, the presence of the slot does not result in a spatially varying precursor profile at the substrate for any of the simulated species, and thus we have indicated the average values of precursor concentration. Nor does it affect the average surface concentration of the reacting species compared with the non-slotted configuration (MoO₃/S₈ ratio (Fig. 2a, 2c), O₂ (Fig. 2b, 2d)). A similar trend was observed throughout the growth time (see Supporting Information section IV and Videos V2, V3, and V4 for simulated profiles along the length of the reactor and at other time steps during growth).

Precursor concentration gradients are not the underlying reason for the observed drastic changes in growth between slotted and non-slotted configurations. Concentration, velocity fields, and temperature are the primary parameters affecting nucleation and growth, amongst which concentration and temperature are similar for both slotted and non-slotted configurations. Hence, we must take a careful look at the velocity profiles near the substrate.

First, we calculate the Reynolds and Mach numbers to understand the type of carrier gas flow (laminar vs. turbulent, and incompressible vs. compressible, respectively) in our reactor. The Reynolds number ($\:Re$) is a nondimensional number that quantifies the ratio of inertial to viscous forces in a flow field, calculated as

$$\:\begin{array}{c}Re=\:\frac{\rho\:uL}{\mu\:}\#\left(3\right)\end{array}$$

where $\:\rho\:$ is carrier gas density (kg/m³), $\:u$ is carrier gas speed (m/s), $\:L$ is characteristic length (m) (the diameter of the tube here), and $\:\mu\:$ is carrier gas dynamic viscosity (kg/(m.s)). A lower value of Re (typically < 10³) predicts the flow to be laminar^44,45. Further, the Mach number is defined as

$$\:\begin{array}{c}M=\:\frac{u}{c}\#\left(4\right)\end{array}$$

where $\:u$ is carrier gas flow speed (m/s) and $\:c$ is speed of sound in the inert carrier gas (m/s). Typical gas flow speeds in our CVD reactor are of the order of 10^− 2-10^− 3 m/s, and with ambient viscosity. The above equations lead to Re ~ 10 and Mach number ~ 10^− 4-10^− 5 and thus, the carrier gas flow should be laminar and incompressible. Although the calculated Mach number is low, we have employed the compressible fluid flow sub-module, because large temperature variations in the system (as revealed by our multiphysics simulations) can cause significant density changes inside the CVD reactor.

Now, we simulate the carrier gas velocity streamlines (while accounting for simulated temperature profiles), representative of the fluid flow direction at any given point. In any given volume element, the evaporated precursor atoms will be imparted momentum in the direction tangential to its nearby streamlines. Hence, streamlines serve as useful tools for visualizing and understanding the local growth environment near the substrate. For discussion, we have indicated (Fig. 3) only those set of streamlines that impinge on the surface of the substrate, and ignored the ones which fall outside or go around.

A comparison between Fig. 3a (slotted) and Fig. 3b (non-slotted) reveals that the carrier gas, upon encountering the substrate, is pushed down toward the slot, whereby its velocity decreases due to the slot geometry. Thus, a region of very low velocity, hereafter referred to as a “velocity dead-zone” (stationary zone) is created near the slot. Due to the lower flow velocity and laminar nature of fluid flow, streamlines bend inside the slot (Fig. 3d). For the non-slotted configuration (Fig. 3c), streamlines show a reversal of the flow at bottom of the substrate, leading the flow away from the surface before turning sideways to continue the forward flow. This reverse current is suppressed in Fig. 3d due to the slot geometry confining the flow inside it. This can directly affect the nucleation density at the substrate, as we will see later.

Next, we examined the velocity profiles directly at the substrate surface. Note that during heating and growth, the velocity inside the reactor is in a steady-state, and hence we plot the steady-state magnitudes of the velocity components. Figure 4a and 4d show the transverse velocity profiles over the substrate, distinctly showing an overall reduction in magnitude in the case when a slot is present. This effect is more pronounced in the region of the substrate physically inside the slot, where the velocity is almost two orders of magnitude lower than in the nearby region (physically outside the slot). Figure 4b and 4e show surface plots of the normal component of the velocity at the substrate. It can be seen in Fig. 4b that a lower velocity region exists toward the bottom of the substrate where it is in contact with the alumina plate. With the slotted configuration in Fig. 4e, this area extends out. The significant changes in the local growth environment due to modified gas flow dynamics should directly affect the delivery of precursor on the substrate.

To corroborate our understanding of the non-trivial influence of the carrier gas velocity profile and precursor concentration, we calculated the incoming precursor flux ($\:J$, at the substrate) analytically. The steady-state nucleation rate⁴⁶ goes as $\:{Q}_{0}\propto\:\:{J}^{2}$ (see Supporting Information section V for a detailed derivation). Hence, a decreased precursor flux to the substrate would result in an overall reduction in the nucleation density.

Using data obtained from the multiphysics simulations, particularly precursor concentration, surface temperature during the estimated nucleation window (roughly between 1–5 minutes of growth, see Supporting Information section IV), and the velocity components at the substrate, we obtained the precursor flux at each substrate surface element (governed by the mesh size employed for simulations). Earlier studies have indicated that TMD nucleation is governed primarily by the concentration and adsorption-desorption equilibrium of the transition metal precursor^15,47,48. With this in mind, Fig. 4c and 4f show the calculated MoO₃ flux, incoming at the substrate at 5 minutes of growth, for the non-slotted and slotted configurations, respectively (see Supporting Information section V for MoO₃ flux profiles throughout the growth, S flux profiles, and details of the flux calculations). In the presence of a slot, the MoO₃ flux reduces to 1/3rd of the value compared to the case without slot. Interestingly, the spatial variation in the precursor flux and the surface velocity follows a similar trend.

To directly measure the nucleation density, we performed scanning electron microscopy (SEM) of the substrates. Figure 4g, 4h and 4i, 4j correspond to slotted and non-slotted configurations, respectively, under similar optimized growth conditions (see Supporting Information section VI for comparison at other vertical distances from the slot). A comparison of Figs. 4g and 4i, which are taken at similar vertical positions on the substrate, well away from the slot, shows similar density, size and shape of nucleation sites. A more careful observation reveals the presence of small MoS₂ flakes with a triangular shape. Due to an increased rate of nucleation in this region, the growth of MoS₂ flakes is limited to a maximum edge length of ~ 1–2 µm. On the contrary, Figs. 4h (in the slot) and 4j (without slot), which are at the same vertical positions on the substrate, show remarkably different depositions. In Fig. 4h, the nucleation density is reduced, and MoS₂ triangular flakes with edge length $\:\ge\:\:$30 µm are observed. Interestingly, for the substrate kept without a slot (Fig. 4j), a deposition similar to Fig. 4g and 4i is observed. This is direct evidence that a slot-induced, confined space reaction environment leads to a decrease in the nucleation density. Hence, reduced precursor flux leads to reduced nucleation density in the slot. Thus, precursor flux, and not just the precursor concentration, is the decisive parameter controlling nucleation.

Predicting temperature variations in the CVD reactor. Temperature has a significant impact on nucleation and ML growth. Even in analytical calculations for the impingement rates discussed above, the rate of nucleation depends exponentially on the temperature of the substrate. Thus, a complete understanding of the dynamic heating process and local temperatures in a CVD reactor is necessary. In our COMSOL model, we accounted for the heating rate and heat diffusion between neighboring heating zones, as well as realistic thermal conductivity of materials in the reactor (alumina, quartz, and inert gas). Such a model is valuable even for simulating the post-deposition cooling process, with an effect on uncontrolled deposition.

Figure 5 and Video V5 depict the time-dependent simulated axial temperature profile along the reactor. The set temperatures of all three zones are achieved in 30 minutes (onset of the growth time in our CVD process). We observe that there is finite heat leakage via the insulation zones, and the temperatures of the areas near to the insulation zone are less than at the centers of the heating zones. Further, due to the presence of heating sources at the radial boundaries, a radial temperature gradient forms in the reactor (see Supporting Information section VII).

Remarkably, the temperature (Fig. 5c) of the substrate assembly (combination of substrate, alumina boat, and plate) is much lower than the center of the heating zone. As a critical experimental insight, we observe that by the onset of growth, the substrate is still not at the desired temperature, and doesn’t reach the set temperature even till the growth is finished (Fig. 5e and Supporting Information section VII). We note that such an observation is very challenging to be made experimentally and leads to a more comprehensive understanding of the observed growth. Through careful selection of holder materials and adjustment of the growth process, such effects may be minimized, or even tuned.

We attribute this significant time-lag to the large heat capacity of the substrate assembly, compared to the inert environment. The specific heat capacities of inert gas (N₂) and alumina are similar. However, the large difference in density leads to a larger heat capacity for alumina boats and supports, compared to that of the surrounding inert atmosphere (see Supporting Information section VII). Figure 5e shows time-dependent plots for the heat source and the average surface temperature at the substrate, with and without alumina supports. The average temperature at the substrate significantly lags behind the source due to the alumina supports, while it follows the source temperature closely in the absence of alumina (with air being replaced for alumina in the simulations) as supporting material. We note that “without alumina supports” is a hypothetical setting used in the simulation, where the substrate is unsupported. In reality, such a system can be created by using a holder with one or a few contact points on the substrate.

The large temperature differences between the substrate and environment have profound effects on the growth process, especially on the nucleation happening at the initial stage of growth. Similarly, such temperature variation could also be found in the precursor alumina boats, which implies that the actual temperature of precursor evaporation may be much lower than expected. Also, large temperature gradients might lead to deposition of unreacted material, and bulk/secondary deposition resulting in poor quality of growth. In short, low temperature differentials are expected to lead to controllable film growth, and conversely, high temperature differentials can lead to particle formation but higher growth rates. The significant temperature differentials between the substrate and local inert gas environment may also give rise to thermo-diffusion or the Ludwig-Soret effect^49,50, which can lead to the aggregation of heavier (MoO₃) and lighter (S) precursors in the colder and warmer regions of the as-formed thermal boundary layer close to the substrate^51,52.

In summary, we report the non-trivial effect of transition metal precursor flux in modulating the 2D TMD nucleation density in a confined-space microreactor (vertically placed substrate in a slot). This unique insight is enabled by the development of a hyper-realistic digital-twin CFD model, accounting for specific reactor and substrate holder geometries, evaporation and transport of precursor species, and reactor heating. Significant modifications in the carrier gas flow near the slot (velocity-dead zone), revealed by the digital-twin, impact precursor delivery during the nucleation and growth stages. Thus, the precursor flux is the critical parameter controlling nucleation, wherein the reduction of nucleation density is experimentally verified by electron microscopy. Further, we predict a significant time lag between the actual substrate temperature and the set temperature, explained by the high heat capacity of substrate holder. Thus, such a multiscale model offers a better understanding of the actual growth temperature and local temperature variations at different places in the reactor, allowing local tunability and control.

It should be noted that the digital-twin model is centered on transport phenomena (heat/mass transfer and fluid flow) and is independent of specific reaction kinetics. This allows for a generally applicable understanding of how macroscopic experimental parameters may affect microscopic variables such as the residence time and nucleation density. This digital-twin model may also be integrated with density functional theory (DFT) calculations and reactive molecular dynamics (MD) simulations to relate a set of macroscopic growth parameters to microscopic physical quantities. These simulations can be used to relate optical and structural properties (such as defects), to growth conditions and reactor configurations. This study motivates exploration of complex substrate holder geometries, and introducing additional physical features to the reactor to controllably modify the local growth environment for high-quality material synthesis. Since such studies are challenging to do in practice, the as-simulated profiles will lead to informed decision making. Overall, our combined experimental-theoretical work paves the way for informed, predictive, and reproducible growth of high-quality and large-area 2D materials via regulating precursor flux, direct control of nucleation density, and the development of a hyper-realistic digital-twin model.

Experimental details. Monolayer MoS₂ was grown in a split tube type, horizontal, three-zone furnace. Sulfur (Sigma-Aldrich, CAS no. 7704-34-9, 99.98%) and MoO₃ (Sigma-Aldrich, CAS no. 1313-27-5, 99.97%) powder precursors were used, kept in the first and second zone, respectively. The substrate was placed vertically inside a slot (groove/cavity) in an alumina plate, supported by another alumina block. Ultra-high purity inert gas (N₂) was used to flush the reactor of ambient air, moisture, and other contaminants. The exhaust is open to atmosphere and the process takes place at atmospheric pressure with constant N₂ flow.

The three zones are individually heated to 200°C, 530°C, and 750°C in 30 minutes, with 100/2 sccm N₂/O₂ flow. During growth, the flow rate is 100/2 sccm (N₂/O₂) for the first 5 minutes and 100 sccm (N₂) for 15 minutes. After 20 minutes of growth, the furnace lid was opened for rapid cooling and 200 sccm N₂ was flown to quench the growth. Surface coverage and morphology of the as-grown ML MoS₂ was observed using a scanning electron microscope (Ultra55 FE-SEM Karl Zeiss) at 5 kV accelerating voltage. Raman and photoluminescence (PL) measurements were carried out in a HORIBA LabRamHR Raman/PL setup using a 532 nm laser.

Multiphysics simulations. Simulations were performed using COMSOL Multiphysics, which utilized the CFD, Heat transfer, and Chemical Reaction Engineering modules. The submodules, Laminar Flow, Heat transfer for fluids, Transport of diluted species, and Transport of concentrated species were used to model the physics inside the CVD tube. The simulation was performed on a mesh termed "Coarser" in COMSOL with a total of 101905 elements within the volume of the tube and with an average element quality of 0.6268 in skewness (0 being the worst and 1 being the best). The mesh was fine enough to realize the slot in the model. Running each simulation on a workstation with a 12th Gen Intel i9 with 16 cores took three and a half hours.

ASSOCIATED CONTENT

The following files are available free of charge.

Supporting Information (PDF)

Section I: Inert gas flushing of CVD reactor before growth; Section II: Comparison of horizontal vs. vertical substrate placement; Section III: Characterization of as-grown monolayer MoS₂; Section IV: Details of precursor evaporation and related modeling; Section V: Calculating precursor adatom flux at the substrate for slotted and non- slotted configurations; Section VI: Scanning electron microscopy (SEM) for comparing nucleation density; Section VII: Temperature profile along the reactor and at the substrate during heating and growth.

Supporting Information Videos

Video V1: Carrier gas (N₂) concentration during heating and growth (mol/m³); Video V2: MoO₃ concentration profile during growth (mol/m³); Video V3: S₈ concentration profile during growth (mol/m³); Video V4: O₂ concentration profile during heating and growth (mol/m³); Video V5: Temperature profile during heating and growth (K).

2D, two-dimensional; TMD, transition metal dichalcogenide; ML, monolayer; CVD, chemical vapor deposition, CFD, computational fluid dynamics; MoS₂, molybdenum disulfide; DFT, density functional theory; MOCVD, metal-organic chemical vapor deposition; AP-CVD, atmospheric pressure chemical vapor deposition; SEM, scanning electron microscopy; DCE, diffusion-convection equation; RMD, reactive molecular dynamics

Author Contributions

Abhinav Sinha (ABS), MV, AS, and AGR developed the overall modeling framework. The inert gas flushing and precursor evaporation model was developed by MV, ABS, AS, and AGR. The temperature model was developed by ABS, KSK, and AS. All simulations were performed by ABS and NK. MV, KSK, and AS developed the experimental framework. MV and KSK performed chemical vapor deposition. MV and KSK carried out optical microscopy, scanning electron microscopy and optical spectroscopy characterization. MV, ABS, AS, and AGR discussed and prepared the manuscript with contributions from all authors.

Data Availability

All data is available upon reasonable request.

ACKNOWLEDGMENTS

AS acknowledges funding from DST Nanomission CONCEPT (Consortium for Collective and Engineered Phenomena in Topology) grant and Anusandhan National Research Foundation (ANRF) grant (CRG-2022-003334). The authors acknowledge the financial support received from the I-STEM program funded by the Office of the Principal Scientific Adviser to the Government of India, for providing the COMSOL license. ABS and NK acknowledge the Kotak IISc AI-ML Centre (KIAC-UG Computation Lab) for computational resources. AGR acknowledges financial support from the Department of Science and Technology (DST) and the Ministry of Electronics and Information Technology (MeitY) via the National Supercomputing Mission (NSM) under grant DST/NSM/R&D_HPC_Applications/Extension Grant/2023/16. AGR also thanks the Infosys Foundation, Bengaluru for an Infosys Young Investigator grant. KSK acknowledges the DST-INSPIRE Fellowship. The authors also acknowledge the Micro Nano Characterization Facility (MNCF), Centre for Nano Science and Engineering (CeNSE) for use of characterization facilities. The authors thank Madhavan Nampoothiri for useful discussions regarding the precursor flux calculations.

Datta, I.; Chae, S. H.; Bhatt, G. R.; Tadayon, M. A.; Li, B.; Yu, Y.; Park, C.; Park, J.; Cao, L.; Basov, D. N.; Hone, J.; Lipson, M. Low-Loss Composite Photonic Platform Based on 2D Semiconductor Monolayers. Nat. Photonics 2020, 14 (4), 256–262. https://doi.org/10.1038/s41566-020-0590-4.
Singh, A.; Jo, S. S.; Li, Y.; Wu, C.; Li, M.; Jaramillo, R. Refractive Uses of Layered and Two-Dimensional Materials for Integrated Photonics. ACS Photonics 2020, 7 (12), 3270–3285. https://doi.org/10.1021/acsphotonics.0c00915.
Pak, S.; Jang, A.-R.; Lee, J.; Hong, J.; Giraud, P.; Lee, S.; Cho, Y.; An, G.-H.; Lee, Y.-W.; Suk Shin, H.; M. Morris, S.; Cha, S.; Inn Sohn, J.; Min Kim, J. Surface Functionalization-Induced Photoresponse Characteristics of Monolayer MoS₂ for Fast Flexible Photodetectors. Nanoscale 2019, 11 (11), 4726–4734. https://doi.org/10.1039/C8NR07655C.
Zhao, Y.; Chen, Y.; Ou, P.; Song, J. Basal Plane Activation via Grain Boundaries in Monolayer MoS₂ for Carbon Dioxide Reduction. ACS Catal. 2023, 13 (19), 12941–12951. https://doi.org/10.1021/acscatal.3c03113.
Zhu, J.; Wang, Z.-C.; Dai, H.; Wang, Q.; Yang, R.; Yu, H.; Liao, M.; Zhang, J.; Chen, W.; Wei, Z.; Li, N.; Du, L.; Shi, D.; Wang, W.; Zhang, L.; Jiang, Y.; Zhang, G. Boundary Activated Hydrogen Evolution Reaction on Monolayer MoS₂. Nat. Commun. 2019, 10 (1), 1348. https://doi.org/10.1038/s41467-019-09269-9.
Ottaviano, L.; Palleschi, S.; Perrozzi, F.; D’Olimpio, G.; Priante, F.; Donarelli, M.; Benassi, P.; Nardone, M.; Gonchigsuren, M.; Gombosuren, M.; Lucia, A.; Moccia, G.; Cacioppo, O. A. Mechanical Exfoliation and Layer Number Identification of MoS₂ Revisited. 2D Mater. 2017, 4 (4), 045013. https://doi.org/10.1088/2053-1583/aa8764.
Huang, J.-K.; Pu, J.; Hsu, C.-L.; Chiu, M.-H.; Juang, Z.-Y.; Chang, Y.-H.; Chang, W.-H.; Iwasa, Y.; Takenobu, T.; Li, L.-J. Large-Area Synthesis of Highly Crystalline WSe₂ Monolayers and Device Applications. ACS Nano 2014, 8 (1), 923–930. https://doi.org/10.1021/nn405719x.
Li, T.; Guo, W.; Ma, L.; Li, W.; Yu, Z.; Han, Z.; Gao, S.; Liu, L.; Fan, D.; Wang, Z.; Yang, Y.; Lin, W.; Luo, Z.; Chen, X.; Dai, N.; Tu, X.; Pan, D.; Yao, Y.; Wang, P.; Nie, Y.; Wang, J.; Shi, Y.; Wang, X. Epitaxial Growth of Wafer-Scale Molybdenum Disulfide Semiconductor Single Crystals on Sapphire. Nat. Nanotechnol. 2021, 16 (11), 1201–1207. https://doi.org/10.1038/s41565-021-00963-8.
Pengfei Yan, Jin Wang, Guofeng Yang, Naiyan Lu, Guangyong Chu, Xiumei Zhang, Xiaowen Shen, Chemical Vapor Deposition of Monolayer MoS₂ on Sapphire, Si and GaN Substrates, Superlattices and Microstructures, Volume 120, 2018, Pages 235-240, ISSN 0749-6036, Https://Doi.Org/10.1016/j.Spmi.2018.05.049.
Tang, L.; Tan, J.; Nong, H.; Liu, B.; Cheng, H.-M. Chemical Vapor Deposition Growth of Two-Dimensional Compound Materials: Controllability, Material Quality, and Growth Mechanism. Acc. Mater. Res. 2021, 2 (1), 36–47. https://doi.org/10.1021/accountsmr.0c00063.
Senkić, A.; Bajo, J.; Supina, A.; Radatović, B.; Vujičić, N. Effects of CVD Growth Parameters on Global and Local Optical Properties of MoS₂ Monolayers. Mater. Chem. Phys. 2023, 296, 127185. https://doi.org/10.1016/j.matchemphys.2022.127185.
Cho, Y. J.; Sim, Y.; Lee, J.-H.; Hoang, N. T.; Seong, M.-J. Size and Shape Control of CVD-Grown Monolayer MoS₂. Curr. Appl. Phys. 2023, 45, 99–104. https://doi.org/10.1016/j.cap.2022.11.008.
Liu, H.; Zhu, Y.; Meng, Q.; Lu, X.; Kong, S.; Huang, Z.; Jiang, P.; Bao, X. Role of the Carrier Gas Flow Rate in Monolayer MoS₂ Growth by Modified Chemical Vapor Deposition. Nano Res. 2017, 10 (2), 643–651. https://doi.org/10.1007/s12274-016-1323-3.
Zheng, B.; Chen, Y. Controllable Growth of Monolayer MoS₂ and MoSe₂ Crystals Using Three-Temperature-Zone Furnace. IOP Conf. Ser. Mater. Sci. Eng. 2017, 274, 012085. https://doi.org/10.1088/1757-899X/274/1/012085.
Wang, W.; Zeng, X.; Wu, S.; Zeng, Y.; Hu, Y.; Ding, J.; Xu, S. Effect of Mo Concentration on Shape and Size of Monolayer MoS₂ Crystals by Chemical Vapor Deposition. J. Phys. Appl. Phys. 2017, 50 (39), 395501. https://doi.org/10.1088/1361-6463/aa81ae.
Arafat, A.; Islam, M. S.; Ferdous, N.; Islam, A. S. M. J.; Sarkar, M. M. H.; Stampfl, C.; Park, J. Atomistic Reaction Mechanism of CVD Grown MoS₂ through MoO₃ and H₂S Precursors. Sci. Rep. 2022, 12 (1), 16085. https://doi.org/10.1038/s41598-022-20531-x.
Zhu, D.; Shu, H.; Jiang, F.; Lv, D.; Asokan, V.; Omar, O.; Yuan, J.; Zhang, Z.; Jin, C. Capture the Growth Kinetics of CVD Growth of Two-Dimensional MoS₂. Npj 2D Mater. Appl. 2017, 1 (1), 1–8. https://doi.org/10.1038/s41699-017-0010-x.
Zafar, A.; Zafar, Z.; Zhao, W.; Jiang, J.; Zhang, Y.; Chen, Y.; Lu, J.; Ni, Z. Sulfur-Mastery: Precise Synthesis of 2D Transition Metal Dichalcogenides. Adv. Funct. Mater. 2019, 29 (27), 1809261. https://doi.org/10.1002/adfm.201809261.
Momeni, K.; Ji, Y.; Chen, L.-Q. Computational Synthesis of 2D Materials Grown by Chemical Vapor Deposition. J. Mater. Res. 2022, 37 (1), 114–123. https://doi.org/10.1557/s43578-021-00384-2.
Mishra, P.; Verma, N. A CFD Study on a Vertical Chemical Vapor Deposition Reactor for Growing Carbon Nanofibers. Chem. Eng. Res. Des. 2012, 90 (12), 2293–2301. https://doi.org/10.1016/j.cherd.2012.05.006.
Chakraborty, P.; Zachariah, M. R. Sticking Coefficient and Processing of Water Vapor on Organic-Coated Nanoaerosols. J. Phys. Chem. A 2008, 112 (5), 966–972. https://doi.org/10.1021/jp076442f.
Govind Rajan, A.; Warner, J. H.; Blankschtein, D.; Strano, M. S. Generalized Mechanistic Model for the Chemical Vapor Deposition of 2D Transition Metal Dichalcogenide Monolayers. ACS Nano 2016, 10 (4), 4330–4344. https://doi.org/10.1021/acsnano.5b07916.
Ghorai, S.; Govind Rajan, A. From Molecular Precursors to MoS₂ Monolayers: Nanoscale Mechanism of Organometallic Chemical Vapor Deposition. Chem. Mater. 2024, 36 (6), 2698–2710. https://doi.org/10.1021/acs.chemmater.3c02675.
Bhowmik, S.; Govind Rajan, A. Chemical Vapor Deposition of 2D Materials: A Review of Modeling, Simulation, and Machine Learning Studies. iScience 2022, 25 (3), 103832. https://doi.org/10.1016/j.isci.2022.103832.
Johari, M. H.; Sirat, M. S.; Mohamed, M. A.; Mustaffa, A. F.; Mohmad, A. R. Computational Fluid Dynamics Insights into Chemical Vapor Deposition of Homogeneous MoS₂ Film with Solid Precursors. Cryst. Res. Technol. n/a (n/a), 2300139. https://doi.org/10.1002/crat.202300139.
Momeni, K.; Ji, Y.; Nayir, N.; Sakib, N.; Zhu, H.; Paul, S.; Choudhury, T. H.; Neshani, S.; van Duin, A. C. T.; Redwing, J. M.; Chen, L.-Q. A Computational Framework for Guiding the MOCVD-Growth of Wafer-Scale 2D Materials. Npj Comput. Mater. 2022, 8 (1), 1–8. https://doi.org/10.1038/s41524-022-00936-y.
Suleman, M.; Lee, S.; Kim, M.; Nguyen, V. H.; Riaz, M.; Nasir, N.; Kumar, S.; Park, H. M.; Jung, J.; Seo, Y. NaCl-Assisted Temperature-Dependent Controllable Growth of Large-Area MoS₂ Crystals Using Confined-Space CVD. ACS Omega 2022, 7 (34), 30074–30086. https://doi.org/10.1021/acsomega.2c03108.
Tian, D.; Kang, L.; Zhang, Z.; Meng, L.; Li, X.; Yan, W. Controlled Synthesis of Continuous MoS₂ Films via Space-Confined Vapor Deposition. Chem. Phys. 2023, 571, 111923. https://doi.org/10.1016/j.chemphys.2023.111923.
Ryu, B.; Liu, Y.; Pu, H.; Wang, Y.; Shi, X.; Soo Lee, S.; Chen, J. A Facile Approach for Site-Selective and Large-Area Growth of MoS₂ through Heterogeneous Nucleation. Appl. Surf. Sci. 2023, 607, 155066. https://doi.org/10.1016/j.apsusc.2022.155066.
Duan, Z.; Chen, T.; Shi, J.; Li, J.; Song, K.; Zhang, C.; Ding, S.; Li, B.; Wang, G.; Hu, S.; He, X.; He, C.; Xu, H.; Liu, X.; Jin, C.; Zhong, J.; Hao, G. Space-Confined and Substrate-Directed Synthesis of Transition-Metal Dichalcogenide Nanostructures with Tunable Dimensionality. Sci. Bull. 2020, 65 (12), 1013–1021. https://doi.org/10.1016/j.scib.2020.03.002.
Tang, L.; Li, T.; Luo, Y.; Feng, S.; Cai, Z.; Zhang, H.; Liu, B.; Cheng, H.-M. Vertical Chemical Vapor Deposition Growth of Highly Uniform 2D Transition Metal Dichalcogenides. ACS Nano 2020, 14 (4), 4646–4653. https://doi.org/10.1021/acsnano.0c00296.
Sharifi, Y.; Achenie, L. E. K. Effect of Substrate Geometry on the Deposition Rate in Chemical Vapor Deposition. J. Cryst. Growth 2007, 304 (2), 520–525. https://doi.org/10.1016/j.jcrysgro.2007.03.046.
Huet, B.; Zhang, X.; Redwing, J. M.; Snyder, D. W.; Raskin, J.-P. Multi-Wafer Batch Synthesis of Graphene on Cu Films by Quasi-Static Flow Chemical Vapor Deposition. 2D Mater. 2019, 6 (4), 045032. https://doi.org/10.1088/2053-1583/ab33ae.
Kim, H.; Ovchinnikov, D.; Deiana, D.; Unuchek, D.; Kis, A. Suppressing Nucleation in Metal–Organic Chemical Vapor Deposition of MoS₂ Monolayers by Alkali Metal Halides. Nano Lett. 2017, 17 (8), 5056–5063. https://doi.org/10.1021/acs.nanolett.7b02311.
Chen, W.; Zhao, J.; Zhang, J.; Gu, L.; Yang, Z.; Li, X.; Yu, H.; Zhu, X.; Yang, R.; Shi, D.; Lin, X.; Guo, J.; Bai, X.; Zhang, G. Oxygen-Assisted Chemical Vapor Deposition Growth of Large Single-Crystal and High-Quality Monolayer MoS₂. J. Am. Chem. Soc. 2015, 137 (50), 15632–15635. https://doi.org/10.1021/jacs.5b10519.
Zhou, D.; Shu, H.; Hu, C.; Jiang, L.; Liang, P.; Chen, X. Unveiling the Growth Mechanism of MoS₂ with Chemical Vapor Deposition: From Two-Dimensional Planar Nucleation to Self-Seeding Nucleation. Cryst. Growth Des. 2018, 18 (2), 1012–1019. https://doi.org/10.1021/acs.cgd.7b01486.
S Kumar, K.; Dash, A. K.; Sabreen H, H.; Verma, M.; Kumar, V.; Watanabe, K.; Taniguchi, T.; Gautam, G. S.; Singh, A. Understanding Interplay of Defects, Oxygen and Strain in 2D Materials for next-Generation Optoelectronics. 2D Mater. 2024. https://doi.org/10.1088/2053-1583/ad4e44.
Mondal, M.; Dash, A. K.; Singh, A. Optical Microscope Based Universal Parameter for Identifying Layer Number in Two-Dimensional Materials. ACS Nano 2022, 16 (9), 14456–14462. https://doi.org/10.1021/acsnano.2c04833.
Hirschfelder, J. O.; Curtiss, C. F.; Bird, R. B. The Molecular Theory of Gases and Liquids; Wiley, 1954.
Gulbransen, E. A.; Andrew, K. F.; Brassart, F. A. Vapor Pressure of Molybdenum Trioxide. J. Electrochem. Soc. 1963, 110 (3), 242. https://doi.org/10.1149/1.2425719.
Peng, D.-Y.; Zhao, J. Representation of the Vapour Pressures of Sulfur. J. Chem. Thermodyn. 2001, 33 (9), 1121–1131. https://doi.org/10.1006/jcht.2001.0835.
Ohring, M. Chapter 3 - Thin-Film Evaporation Processes. In Materials Science of Thin Films (Second Edition); Ohring, M., Ed.; Academic Press: San Diego, 2002; pp 95–144. https://doi.org/10.1016/B978-012524975-1/50006-9.
Lyons, J. R. An Estimate of the Equilibrium Speciation of Sulfur Vapor over Solid Sulfur and Implications for Planetary Atmospheres. J. Sulfur Chem. 2008, 29 (3–4), 269–279. https://doi.org/10.1080/17415990802195615.
Holman, J.P. (2009) Heat Transfer. 10th Edition, McGraw-Hill, New York. - References - Scientific Research Publishing. https://scirp.org/reference/referencespapers?referenceid=2298536 (accessed 2024-06-26).
Schlichting, H.; Gersten, K. Boundary-Layer Theory; Springer Science & Business Media, 2003.
Markov, I. V. Crystal Growth for Beginners: Fundamentals of Nucleation, Crystal Growth and Epitaxy; World Scientific, 2003.
Wang, T.; Guo, J.; Zhang, Y.; Yang, W.; Feng, X.; Li, J.; Yang, P. Synthesis of High-Quality Monolayer MoS₂ via a CVD Upstream Deposition Strategy for Charge Capture and Collection. Cryst. Growth Des. 2024, 24 (7), 2755–2763. https://doi.org/10.1021/acs.cgd.3c01369.
Ding, B.; Zhu, C.; Wang, T.; Li, L.; Li, Z.; Cheng, L.; Feng, S.; Zhang, G.; Zang, Y.; Hu, J.; Li, L.; Xia, C. Large-Area Growth of Monolayer MoS₂ by Using Atmospheric-Pressure Chemical Vapor Deposition with Nucleation Controlling Process. Surf. Interface Anal. 2024, 56 (1), 22–31. https://doi.org/10.1002/sia.7264.
Rahman, M. A.; Saghir, M. Z. Thermodiffusion or Soret Effect: Historical Review. Int. J. Heat Mass Transf. 2014, 73, 693–705. https://doi.org/10.1016/j.ijheatmasstransfer.2014.02.057.
Kafoussias, N. G.; Williams, E. W. Thermal-Diffusion and Diffusion-Thermo Effects on Mixed Free-Forced Convective and Mass Transfer Boundary Layer Flow with Temperature Dependent Viscosity. Int. J. Eng. Sci. 1995, 33 (9), 1369–1384. https://doi.org/10.1016/0020-7225(94)00132-4.
YEH, H.-M. Enrichment of Heavy Water by Thermal Diffusion. Chem. Eng. Commun. 1998, 167 (1), 167–179. https://doi.org/10.1080/00986449808912699.
Ma, N.-Y. R.; Beyerlein, A. L. The 2H Thermal Diffusion Isotope Effect in Benzene and Methanol. J. Chem. Phys. 1983, 78 (11), 7010–7011. https://doi.org/10.1063/1.444651.

There is NO Competing Interest.

VideoV1.mp4
Carrier gas (N2) concentration during heating and growth (mol/m3)
VideoV2.mp4
MoO3 concentration profile during growth (mol/m3)
VideoV3.mp4
S8 concentration profile during growth (mol/m3)
VideoV4.mp4
O2 concentration profile during heating and growth (mol/m3)
VideoV5.mp4
Temperature profile during heating and growth (K)
SupportingInformationSinhaVerma.docx

Download PDF

Version 1

posted

You are reading this latest preprint version

Critical Role of Precursor Flux in Modulating Nucleation Density in 2D Material Synthesis Revealed by Digital Twin

Status:

Version 1

Abstract

Figures

Introduction

Results and Discussion

Conclusions

Methods

ASSOCIATED CONTENT

Abbreviations

Declarations

References

Additional Declarations

Supplementary Files

Status:

Version 1