  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.
