This is a Stillinger-Weber (SW) parameterization for molybdenum disulfide. The parameters were calibrated to atomic forces from a DFT calculation of a molybdenum disulfide monolayer at 750K, as explained in Wen et al. (2017). Parameters $q_{IJ}$ are set to be zero for all pair-wise interactions and $\gamma$ to be the same for all three-body interactions. Additionally, to explore the full parameter space, $p_{IJ}$ are allowed to take any positive real value. This also allows the two-body interaction term to be a continuous function in $p_{IJ}$, which is needed in the information geometry analysis. The relation between $\sigma_{IJ}$ and the equilibrium lattice constants of the system is also eliminated. Further, we don't require $d\phi_2/dr|_{r=d}=0$ at the equilibrium bond length $d$, i.e., the atoms are not required to be in a pair-wise equilibrium state, which removes the constraint on $B_{IJ}$. Potential fitting was done to minimize the weighted least-squares loss function. We set the tolerance for the force components corresponding to each atom to be 10\% of the force magnitude acting on the same atom, a reasonable choice of fractional tolerance. We use these values as the non-uniform error bars, i.e., inverse weights, in the loss function.