24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3328  |  回复: 17
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

dfjk-123

新虫 (小有名气)

[求助] 相同的vasp版本,不同机器上计算同一个体系,能量相差近10ev

各位:
我在两个不同的LINUX集群下用vasp计算同一个体系,输入文件完全相同。两个集群安装的都是vasp5.2.2版本,结构优化完成后能量居然相差了近10ev,几何结构两个基本相同。大家说说这是怎么回事?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

【答案】应助回帖


WDD880227: 金币+1, 多谢指教 2012-11-28 14:31:47
引用回帖:
7楼: Originally posted by dfjk-123 at 2012-11-27 15:58:54
我后来试了NBANDS数不变,加大了K点,能量上升了。减少K点的话,warning就没有了。
另外一个测试表明:增加了NBANDS数,K点数不变,能量下降了。在此基础上又增加了K点数,能量在此基础上又有下降。搞不清楚K点和 ...

‘我后来试了NBANDS数不变,加大了K点,能量上升了。减少K点的话,warning就没有了。‘
这明显是K点不够,warming没有了是因为在某些K点的最高带上已经没有占据了,但这个没有是错误的提示。

对于K点的够不够你测试的时候不知道的参数就用默认的,比如说NBANDS。测试的时候增加K点的密度,直到能量达到一个稳定的值(误差在0.几个eV内,越小越好)。

总的来讲我认为你的NBANDS数量不够,NELECT/2 + NIONS/2 (non-spinpolarized)这个没错,但你的电子数是不是数对了呢?

» 本帖已获得的红花(最新10朵)

嘟啦啦嘟啦啦嘟
8楼2012-11-27 19:14:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 18 个回答

dfjk-123

新虫 (小有名气)

刚检查了两个地方生成的OUTCAR文件,发现NBANDS数是不一样的,这可能是导致能量差的原因。但是怎么会差这么大呢?我的INCAR中并没有给定NBANDS数量啊。同样版本的vasp,给出的能带数应该一样才对啊。这是为什么呢?
2楼2012-11-24 19:13:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
franch: 金币+2, 谢谢回帖交流, 2012-11-26 17:25:44
你比较的是什么能量?绝对的吗?vasp属于周期性的结构,能量的零点没法定义,在计算能量中有些数值是根据你所用的数据库决定的,但是如果你取相对的能量的话不花有什么变化的
嘟啦啦嘟啦啦嘟
3楼2012-11-25 10:04:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dfjk-123

新虫 (小有名气)

引用回帖:
3楼: Originally posted by 贺仪 at 2012-11-25 10:04:22
你比较的是什么能量?绝对的吗?vasp属于周期性的结构,能量的零点没法定义,在计算能量中有些数值是根据你所用的数据库决定的,但是如果你取相对的能量的话不花有什么变化的

我是取的energy without entropy后面的能量。我又测试了下,将两个测试的NBANDS数设置为一样,计算出的结果基本上差不多了。我认为还是NBANDS的原因。但是,INCAR文件没有设置NBANDS数时,同样的VASP版本,在不同机器上却给的不一样,这是为什么?
另外,资源有限,我的作业需要使用两个机器做计算,如果得出的能量差别如此大,那算出来的相对能量还能准确吗?应该相信哪一个?
4楼2012-11-25 12:08:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见