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.

Constitutive modeling of some 2D crystals: graphene, hexagonal BN, MoS$_2$, WSe$_2$ and NbSe$_2$

Nonlinear Elastic Behavior of Two-Dimensional Molybdenum Disulfide

The Electronic Structure and Elastic Properties of Nbsisb and Nbgesb: A First-Principles Study

A Constitutive Model for Nanomaterials Based on Spatial Secant

A multiscale computational framework for the analysis of graphene involving geometrical and material nonlinearities

Phase field crystal model for heterostructures

Bicrystallography-informed Frenkel-Kontorova model for interlayer dislocations in strained 2D heterostructures

Constitutive Matrices for 32 Typical Classes of Crystalline Solids with Couple Stress, Quadrupole, and Curvature-Based Flexoelectric Effects

Analytical Models for Predicting the Nonlinear Stress–Strain Relationships and Behaviors of Two-Dimensional Carbon Materials

Linear elastic diatomic multilattices: Three-dimensional constitutive modeling and solutions of the shift vector equation

Anisotropy of two-dimensional materials under tension: a Molecular Dynamic study across graphene, hexagonal boron nitride and molybdenum disulphide

An Atomistic-based Finite Deformation Continuum Membrane Model for Monolayer Transition Metal Dichalcogenides

Model for the on-site matrix elements of the tight-binding hamiltonian of a strained crystal: Application to silicon, germanium and their alloys

Computational modeling of the mechanics and fracture of 2D materials with defects and grain boundaries

A higher order Cauchy-Born rule based constitutive model of SWNT

Irreducible Representations for Constitutive Equations of Anisotropic Solids II: Crystal and Quasicrystal Classes D 2m+1d , D 2m+1 and C 2M+1v

The stability and the nonlinear elasticity of 2D hexagonal structures of Si and Ge from first-principles calculations

A new constitutive model for tetragonal energetic single crystals

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

Continuum modeling of twinning, amorphization, and fracture: theory and numerical simulations

Constitutive modeling of deformation behavior of high-entropy alloys with face-centered cubic crystal structure