Modeling Chemical Exfoliation of Non-van der Waals Chromium Sulfides by Machine Learning Interatomic Potentials and Monte Carlo Simulations

TL;DR

This paper introduces a neural network potential framework to model and predict the exfoliation process of non-van der Waals chromium sulfides, revealing structural transformations and strain effects that facilitate ultrathin nanosheet formation.

Contribution

The study develops a transferable neural network potential for non-vdW materials, enabling accurate modeling of atomic migrations and phase transformations during exfoliation, surpassing traditional methods.

Findings

NNP outperforms cluster expansion in accuracy and transferability

Identifies structural transition at Cr₀.₅S with Cr migration to vdW gaps

Strain engineering can promote vdW gap formation in non-vdW CrS₂

Abstract

The chemical exfoliation of non-van der Waals (vdW) materials to ultrathin nanosheets remains a formidable challenge. This difficulty arises from the strong preference of these materials to engage in three-dimensional chemical bonding, resulting in uncontrolled atomic migration into the vdW gaps during cation deintercalation from the bulk structure, ultimately leading to unpredictable structural disorder. We propose a generic framework using neural network potentials (NNPs) to accurately model the widespread nonstoichiometric environments resulting from disordered atomic migrations during exfoliation of non-vdW materials. We apply our framework to investigate the crystal structures and phase transformations occurring during the exfoliation of non-vdW nonstoichiometric Cr$_{(1-x)}$S systems, a compelling material category with substantial potential for two-dimensional (2D) magnetic applications. The efficacy of the NNP outperforms the conventional cluster expansion, exhibiting superior accuracy and transferability to unexplored crystal structures and compositions. By employing the NNP in simulated annealing optimizations, we predict low-energy Cr$_{(1-x)}$S structures anticipated to result from experimental synthesis. A notable structural transition is discerned at the Cr$_{0.5}$S composition, with half of the Cr atoms preferentially migrating to vdW gaps. This aligns with experimental observations in the chemical exfoliation of 2D CrS$_2$, and emphasizes the vital role of excess Cr atoms beyond the Cr/S = $1/2$ composition ratio in stabilizing vdW gaps. Additionally, we utilize the NNP in a vacancy diffusion Monte Carlo simulation to illustrate the impact of lateral compressive strains in catalyzing the formation of vdW gaps within non-vdW CrS$_2$ slabs. This provides a direct pathway for more facile exfoliation of ultrathin nanosheets from non-vdW materials through strain engineering.