求助-VASP计算银镍界面结合能结果为正

英雄布柳赫尔

第一次发帖，有不对的地方希望各位老师批评指正

各位老师们好，最近刚入门VASP第一性原理计算，想仿照文献计算一下Ag（晶面111）和Ni（晶面111）复合界面的结合能（文献里用的是MS CASTEP计算的），我使用的VASP，版本是6.3.0。

计算步骤是先使用MS cleave surface切晶面，然后再适当扩胞减小失配度（没有降到5%以下），用Build Layers功能建立异质结（Ag和Ni厚度都是三层，扩胞为2*2*1），中间加了20埃的真空层防止周期之间相互影响。之后正常转化为POSCAR文件放进VASP里计算，结构优化和静态自洽后从OSZICAR文件里找到体系总能量。根据界面结合能计算公式（复合后的总能量减去复合前原子之和的能量，然后再除以原子个数及接触面积），我单独对Ag晶胞和Ni晶胞进行结构优化和静态自洽，取OSZICAR里的总能量除以原子个数作为单个原子的能量，根据计算公式计算得到的界面结合能为正，说明界面不稳定明显不合理。分析觉得可能是以下几个方面有问题，不知道具体是哪里不对或者是还有其他可能的问题：

1.建模方面失配度较高（手动计算是13%左右），初始结构不稳定

2.计算参数选取不合理

3.界面结合能计算方法有问题

恳请各位老师指点迷津，万分感谢，这里放上Ag-Ni(Ag和Ni各12个原子)和Ag与Ni(均为晶胞四个原子)的OSZICAR截图以及结构优化和静态自洽的INCAR，需要我提供其他文件都可以。