求助：关于物理吸附中对mof结构优化以及磁性的问题

ljn

我研究的物质是含有fe的mof(fe3O16C122H80),然后使用vasp对结构进行优化。模拟了吸附量后与实验值具有较大的差异。
1.我发现在计算结束后的电荷文件里，有几个连在有机苯环上的氢原子跑了出来，不知道究竟是怎么造成的。（我认为是造成与实验值吸附量差异的主要原因） 2.后来在查阅的过程中，发现有老师说结构优化后想获得更好的电子结构，可把ISMEAR=-5。然后在计算资源有限的情况下，我将k点设置为2 2 4。发现模拟得到的吸附量还是与实验值有较大差异。
3.请问含有fe的mof在计算的过程中，是否该考虑磁性。然后又将ISPIN=2,MAGMOM=3*5 16*0 122*0 80*0,发现结果并没有明显改善。我的INCAR文件如下

不知道应该怎么改进，尤其是电荷计算的时候，有哪些参数设置需要改进，烦请各位老师批评指教。（本人学习时间不久，询问问题不合适的地方还请见谅，谢谢大家！）