While this method yields results for every studied system, as reported on Table , the significant mismatch between values computed with different cell types, which should not exist for the crystal, severely questions the results and illustrate once again the tendency of ReaxFF to densify systems. Additionally, the only two low-density states, the crystal and the ab initio glasses with an isotropic cell, have lower values than what was obtained with AIMD, the opposite of what yielded the finite strain difference method with ReaxFF. This seriously discourages using the finite stress difference method with ReaxFF.