24小时热门版块排行榜    

查看: 3557  |  回复: 17

dfjk-123

新虫 (小有名气)

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

各位:
我在两个不同的LINUX集群下用vasp计算同一个体系,输入文件完全相同。两个集群安装的都是vasp5.2.2版本,结构优化完成后能量居然相差了近10ev,几何结构两个基本相同。大家说说这是怎么回事?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

贺仪

铁杆木虫 (著名写手)

【答案】应助回帖


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的回帖

souledge

专家顾问 (著名写手)


WDD880227: 金币+1, 谢谢提示 2012-11-28 14:31:59
如果计算过程没有警告问题,总能(别管你取哪个)有差别都无所谓,只要生成能之类都没变化就可了。
思想重于技巧,内涵重于表象
9楼2012-11-27 19:35:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

longwen36

铁杆木虫 (正式写手)

【答案】应助回帖

引用回帖:
2楼: Originally posted by dfjk-123 at 2012-11-24 19:13:05
刚检查了两个地方生成的OUTCAR文件,发现NBANDS数是不一样的,这可能是导致能量差的原因。但是怎么会差这么大呢?我的INCAR中并没有给定NBANDS数量啊。同样版本的vasp,给出的能带数应该一样才对啊。这是为什么呢?

貌似也我也有这样的问题,
不同的机器NBANDS数不一样,所以算出来的结果也不一样。

INCAR里面没有给NBANDS,它好像会根据你算的CPU数来设定。貌似与并行的CPU数有关。
所以最好在INCAR里面最好固定NBANDS数,就好像不会有这样的问题。

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

16楼2012-12-04 16:53:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

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的回帖

贺仪

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★
franch: 金币+2, 谢谢回帖交流, 2012-11-26 17:25:56
引用回帖:
4楼: Originally posted by dfjk-123 at 2012-11-25 12:08:56
我是取的energy without entropy后面的能量。我又测试了下,将两个测试的NBANDS数设置为一样,计算出的结果基本上差不多了。我认为还是NBANDS的原因。但是,INCAR文件没有设置NBANDS数时,同样的VASP版本,在不同机 ...

你的两个电脑的设置等不一样计算出来的都会有差别,影响很大的就是你编译时以及计算时所调用的数据库。对于绝对能量的比较最好是能够在同一台机械上计算。相对能量到无所谓。
嘟啦啦嘟啦啦嘟
5楼2012-11-25 20:37:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dfjk-123

新虫 (小有名气)

引用回帖:
5楼: Originally posted by 贺仪 at 2012-11-25 20:37:33
你的两个电脑的设置等不一样计算出来的都会有差别,影响很大的就是你编译时以及计算时所调用的数据库。对于绝对能量的比较最好是能够在同一台机械上计算。相对能量到无所谓。...

好的,谢谢。还有个问题想请教下,我在OUTCAR中发现了下面的警告:
|                                                                             |
|           W    W    AA    RRRRR   N    N  II  N    N   GGGG   !!!           |
|           W    W   A  A   R    R  NN   N  II  NN   N  G    G  !!!           |
|           W    W  A    A  R    R  N N  N  II  N N  N  G       !!!           |
|           W WW W  AAAAAA  RRRRR   N  N N  II  N  N N  G  GGG   !            |
|           WW  WW  A    A  R   R   N   NN  II  N   NN  G    G                |
|           W    W  A    A  R    R  N    N  II  N    N   GGGG   !!!           |
|                                                                             |
|      Your highest band is occupied at some k-points! Unless you are         |
|      performing a calculation for an insulator or semiconductor, without    |
|      unoccupied bands, you have included TOO FEW BANDS!! Please increase    |
|      the parameter NBANDS in file 'INCAR' to ensure that the highest band   |
|      is unoccupied at all k-points. It is always recommended to             |
|      include a few unoccupied bands to accelerate the convergence of        |
|      molecular dynamics runs (even for insulators or semiconductors).       |
|      Because the presence of unoccupied bands improves wavefunction         |
|      prediction, and helps to suppress 'band-crossings.'                    |
|      Following all k-points will be listed (with the Fermi weights of       |
|      the highest band given in paranthesis) ... :                           |
|                                                                             |
|                          1       ( 0.12000)                                 |
|                          2       ( 0.05967)                                 |
|                          4       ( 1.07820)                                 |
|                          5       ( 0.22570)                                 |
|                          8       ( 0.19163)                                 |
|                          9       ( 1.07400)                                 |
|                         10       ( 1.19262)                                 |
|                         11       ( 0.11518)                                 |
|                         12       ( 0.13723)                                 |
|                         13       ( 0.37150)                                 |
|                                                                             |
|      The total occupancy of band no.   448 is   0.72092 electrons ...       |
|                                                                             |
这个需要增大NBANDS到多少啊?按照VASP的说明手册,NBANDS = NELECT/2 + NIONS/2 (non-spinpolarized) 这个值是足够的啊,为什么会有这样的警告?
6楼2012-11-27 15:45:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dfjk-123

新虫 (小有名气)

引用回帖:
5楼: Originally posted by 贺仪 at 2012-11-25 20:37:33
你的两个电脑的设置等不一样计算出来的都会有差别,影响很大的就是你编译时以及计算时所调用的数据库。对于绝对能量的比较最好是能够在同一台机械上计算。相对能量到无所谓。...

我后来试了NBANDS数不变,加大了K点,能量上升了。减少K点的话,warning就没有了。
另外一个测试表明:增加了NBANDS数,K点数不变,能量下降了。在此基础上又增加了K点数,能量在此基础上又有下降。搞不清楚K点和NBANDS的关系到底是怎么一回事啊?
多谢回复!
7楼2012-11-27 15:58:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dfjk-123

新虫 (小有名气)

送鲜花一朵
引用回帖:
8楼: Originally posted by 贺仪 at 2012-11-27 19:14:36
‘我后来试了NBANDS数不变,加大了K点,能量上升了。减少K点的话,warning就没有了。‘
这明显是K点不够,warming没有了是因为在某些K点的最高带上已经没有占据了,但这个没有是错误的提示。

对于K点的够不够你 ...

NELECT和NIONS都在OUTCAR文件中可以查到,我没有自己去数。另外,K点收敛性测试,我在1*1单胞的时候做的。变成3*3超胞后,我就将K点除以了3,有必要在超胞时还做K点测试吗? 谢谢回复!
10楼2012-11-28 08:39:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dfjk-123 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 329求调剂 +5 想上学吖吖 2026-03-19 5/250 2026-03-19 23:56 by 23Postgrad
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 一志愿苏州大学材料求调剂,总分315(英一) +3 sbdksD 2026-03-19 3/150 2026-03-19 23:21 by fmesaito
[考研] 本人考085602 化学工程 专硕 +17 不知道叫什么! 2026-03-15 19/950 2026-03-19 15:06 by 尽舜尧1
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 286求调剂 +6 lemonzzn 2026-03-16 10/500 2026-03-19 14:31 by lemonzzn
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 332求调剂 +3 ydfyh 2026-03-17 3/150 2026-03-19 10:14 by 功夫疯狂
[考研] 304求调剂 +12 小熊joy 2026-03-14 13/650 2026-03-18 12:34 by Linda Hu
[考研] 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 轻松不少随
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 考研化学学硕调剂,一志愿985 +4 张vvvv 2026-03-15 6/300 2026-03-17 17:15 by ruiyingmiao
[考研] 材料与化工专硕调剂 +5 heming3743 2026-03-16 5/250 2026-03-17 14:03 by 勇敢太监王公公
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 材料与化工 323 英一+数二+物化,一志愿:哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
信息提示
请填处理意见