The results of the fitted force field were compared to two previously trained ReaxFF force fields . The restraints are applied on S-S bond of dimethyl disulfide (DMDS), S-C bond of dimethyl thioether (DMTE), H-S-H angle of Hydrogen sulfide (H 2 S) and H-S-S-H torsion angle of Hydrogen disulfide (H 2 S 2 ). These potential energy graphs are shown in Fig. .