|
各位老师同学,我们的一个客户反应,其在某商用超算平台进行的vasp计算,所得到的能量与我们编译的vasp计算的能量,相差0.5eV。客户在超算平台提供的vasp(vasp.5.4.4)计算得到的 TOTEN=-107.503 eV。而我们编译的vasp得到结果均为 TOTEN=-108.069 eV。 注:所有输入文件(incar/kpoints/poscar/potcar)均完全一致。
按照我们的理解,正常的vasp编译结果不可能出现如此大的能量差异。因此,我们更换不同vasp版本(5.4.4/6.3.0/6.3.2/6.4.2); 使用不同的编译器(intel2017/intel2018/intel2019/oneapi2022/aocc+aocl),在不同的处理器上进行 (intel / amd)进行了;编译时均使用arch文件夹内的makefile.include.linux_intel 编译模板,未做任何更改。得到的计算结果均为-108.069 eV。
我找出了我博士期间在该超算平台账号,其中就有当时平台给我提供的vasp执行文件(超算平台利用intel2017编译的vasp.5.4.4,仅提供给我vasp_gam、vasp_ncl、vasp_std执行文件),我利用该vasp_std进行了测试,发现结果为TOTEN=-107.503 eV。这结果和我们客户的在该超算计算的结果一致。我们目前无法确认该超算平台的计算结果是否正确,但是我们非常希望找到造成这样差异的根源在哪。以及,正确的结果是-107.503 eV还是-108.069 eV。
我们开始请教专业做vasp计算的老师、同学。得到的信息是:结构优化中,能量越小,可能是越合理的,即-108.069 eV可能是更合理的结果。但是目前问题依旧停留在 “vasp在什么编译条件下或并行环境下,会造成如此的误差 ”。由于本人不是专门做第一项原理计算,对于理论及vasp软件了解不够深刻,特希望在论坛寻求答案。
我列出两种计算结果输入及OUTCAR文件,各位老师同学,可以在自己的平台下测试下该计算的能量结果(使用vasp_std计算,三代志强32核大概需要10min计算结束)。如果能找出这个问题的原因,可能能给到我们一些对于vasp计算的警示,特别是如果vasp编译可能导致计算结果巨大差异,那么我应该谨慎对待各类平台提供的vasp执行文件。
期待各位老师同学的测试及反馈
|
|