Full geometry optimizations were performed within 2D periodic boundary conditions using the LAMMPS code and the Reax force field (ReaxFF) parametrized for MoS2 systems . Geometry convergence was considered to be reached at a maximum force component of 0.1 eV Å -1 and remaining pressure of the atomic system below 10 kbar.