Skip to content

Computational Synthesis of MoS2 Layers by Reactive Molecular Dynamics Simulations: Initial Sulfidation of MoO3 Surfaces

Evidence and attribution

Authority of statements

Prose below (Summary, Methods, Findings) summarizes the publication identified by doi, title, and pdf_path in the front matter. For reaction sequences, numerical protocols, and figures, rely on the peer-reviewed article and Supporting Information.

Summary

Chemical vapor deposition (CVD) is the scalable route to monolayer transition metal dichalcogenides such as MoS\(_2\), but atomistic selection rules and optimal growth conditions remain poorly resolved. This letter reports first-principles-informed ReaxFF reactive molecular dynamics (RMD) of the early sulfidation of an \(\alpha\)-MoO\(_3\) layer on \(\alpha\)-Al\(_2\)O\(_3\)(0001) exposed to S\(_2\), aimed at the predeposited-oxide CVD setups used in experiment. The trajectories resolve a three-stage pathway from MoO\(_3\) and S\(_2\) toward Mo–S bonding: oxygen release and surface self-reduction, SO/SO\(_2\) chemistry with S\(_2\)-assisted reduction, and sulfidation with Mo–S bond formation. The authors relate this sequence to reported experiments and prior electronic-structure literature as a mechanistic input for rational TMD CVD.

Methods

A — Force-field training / fitting: Mo/O and Mo/S ReaxFF from Chenoweth et al. and Ostadhossein et al., reoptimized on QM data for MoO\(_3\)–S\(_2\) chemistry; S/O branch per prior CVD work—full training tables in 2017hong-venue-microsoft-word (SI PDF).

B — Molecular dynamics / sampling: ReaxFF RMD in NVT, Nosé–Hoover (25.0 fs damping), timestep 0.25–0.30 fs. Slab: α-MoO\(_3\)(001) monolayer (1152 atoms) on α-Al\(_2\)O\(_3\)(0001) (3375 atoms), cell 47.09 × 45.57 Ų, ~100 Å vacuum, S\(_2\) gas. Preequilibration: 500 K, 125 ps; cool 100 K over 25 ps. Further heated reactive stages per manuscript/figures.

C — DFT / static QM: QM used to build/refit ReaxFF (see SI); not reported as on-the-fly ab initio MD in the letter summary.

D — Review / non-simulation framing: Primary Nano Lett. application paper—not a review.

Engine: ReaxFF RMD workflow as implemented for this study (see article for code paths). System: α-MoO\(_3\)(001) monolayer (1152 atoms) on α-Al\(_2\)O\(_3\)(0001) (3375 atoms), in-plane cell 47.09 × 45.57 Ų, ~100 Å vacuum gap, with S\(_2\) gas. Boundaries / periodicity: 3D PBC is implied by the slab + vacuum supercell description; confirm fixed substrate layers in the PDF if needed. Ensemble / thermostat: NVT with Nosé–Hoover (25.0 fs damping). Timestep: 0.25–0.30 fs. Duration / staging: 500 K preequilibration for 125 ps; cool 100 K over 25 ps; additional heated reactive segments per figures (full timeline in article). Barostat / pressure: N/A — NVT gas–surface RMD without NPT control in the summarized protocol. Electric field: N/A — not used. Replica / enhanced sampling: N/A — not used.

Findings

The RMD picture is a sequential O\(_2\) evolution and MoO\(_3\) surface reduction, SO/SO\(_2\) formation with S\(_2\)-mediated reduction, then sulfidation and Mo–S bond formation—offered as the atomic-scale early-stage chemistry linking oxide precursors to MoS\(_2\) growth. The authors state consistency with experimental observations and existing theory for sulfidation under CVD-relevant conditions.

Sensitivity / design levers. Temperature ramps and S\(_2\) exposure tie directly to the staged chemistry above; quantitative coverage vs time plots are in the Letter.

Limitations / outlook (as authored). The Summary already flags the predeposited MoO\(_3\) model scope; ReaxFF accuracy limits quantitative barriers and long-time kinetics.

Corpus / PDF honesty. Definitive figure-level chemistry should be checked against pdf_path and 2017hong-venue-microsoft-word SI for the refit tables.

Limitations

  • The model emphasizes predeposited MoO\(_3\) on sapphire rather than co-evaporated MoO\(_3\) and sulfur powders; the paper discusses this distinction and still argues for transferable elementary steps.
  • ReaxFF remains an empirical reactive model: barrier heights and long-time kinetics should be checked against higher-level benchmarks where decisions hinge on quantitative rates.

Relevance to group

MoS\(_2\) CVD chemistry and ReaxFF reoptimization against QM for new oxide–chalcogen interfaces align with broader reactive MD and 2D TMD work in the corpus.

Citations and evidence anchors

Reader notes (navigation)