VASP计算Ni/C(111)界面分离功/结合能求助

liamesiy

求助各位大佬，本人小白想复现一篇文献中DFT计算纯Ni/C(111)面界面结合能，文献使用MS的CASTEP模块计算，我使的是用vasp.6.4.2，但计算出结果比文献小一个数量级，文献为：1.735J/m2，我计算出来为0.1351 J/m2，下面是我的计算流程，求大佬们帮我看看哪些部分有问题
建模：使用MS- cleave surface切晶面，Build Layers功能建立界面模型(见图片)，同文献加了10埃真空层防止周期之间相互影响，固定底下两层Ni和远离界面的三层C，初始界面距离文献中并没有给出，我设置为4A
结构优化中，k点网格8*8*1，动能截止值ENCUT400，均和文献一致，具体参数已附上图片

结构优化收敛后，对Ni/C整体模型、拆开的Ni和C模型分别进行静态计算，算出各自能量，具体INCAR参数也已附上图片

zhrakapeipei

感觉界面间距给太大了，并且结构优化精度低（EDIFFG在-0.01到-0.05比较合适）