24小时热门版块排行榜    

查看: 4236  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 基金啊基金 +3 longfie172 2026-08-20 3/150 2026-08-21 00:04 by wuchongjun
[基金申请] 科研孤儿太难了 +14 我4大白菜 2026-08-20 15/750 2026-08-20 22:33 by 袁向阳007
[基金申请] 让我中一个面上吧! +8 大萍1987 2026-08-20 10/500 2026-08-20 21:48 by 大萍1987
[基金申请] 言之凿凿,明天官网会通知放榜? +15 医学老男孩 2026-08-20 17/850 2026-08-20 20:35 by CaoLG
[基金申请] 时间戳又变了 +12 wuchongjun 2026-08-20 17/850 2026-08-20 20:05 by spzhao
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 只有每年这种时候来逛逛小木虫 +21 yaoyewhu2008 2026-08-20 22/1100 2026-08-20 17:03 by yaoyewhu2008
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
信息提示
请填处理意见