24小时热门版块排行榜    

查看: 3561  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工 求调剂 +17 enenenhui 2026-03-13 18/900 2026-03-20 08:55 by xingguangj
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 311求调剂 +4 冬十三 2026-03-18 4/200 2026-03-18 21:47 by 尽舜尧1
[考研] 085600材料与化工 +5 安全上岸! 2026-03-16 5/250 2026-03-18 15:33 by cmz0325
[考研] 085601专硕,总分342求调剂,地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 070300化学319求调剂 +6 锦鲤0909 2026-03-17 6/300 2026-03-18 13:22 by Iveryant
[考研] 299求调剂 +5 △小透明* 2026-03-17 5/250 2026-03-18 11:49 by 尽舜尧1
[考研] 268求调剂 +6 简单点0 2026-03-17 6/300 2026-03-18 09:04 by 无际的草原
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 332求调剂 +6 Zz版 2026-03-13 6/300 2026-03-17 17:03 by ruiyingmiao
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 复试调剂 +3 呼呼?~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
信息提示
请填处理意见