24小时热门版块排行榜    

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

dfjk-123

新虫 (小有名气)

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

各位:
我在两个不同的LINUX集群下用vasp计算同一个体系,输入文件完全相同。两个集群安装的都是vasp5.2.2版本,结构优化完成后能量居然相差了近10ev,几何结构两个基本相同。大家说说这是怎么回事?
回复此楼
已阅   回复此楼   关注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的回帖
查看全部 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 申博自荐 +6 Linxia林夏 2026-04-13 6/300 2026-04-20 19:26 by 想申博!
[考研] 320求调剂 +6 深郊akm 2026-04-17 6/300 2026-04-20 18:57 by fs26jie
[论文投稿] 期刊推荐 +3 材料研究生 2026-04-15 5/250 2026-04-20 16:02 by 豆豆7758
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +4 yexuqing 2026-04-19 4/200 2026-04-20 14:47 by brantleo
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 5/250 2026-04-20 10:47 by YuY66
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 通信工程求调剂!!! +7 zlb770521 2026-04-14 7/350 2026-04-19 20:56 by Equinoxhua
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂 +10 小聂爱学习 2026-04-16 12/600 2026-04-19 16:51 by 中豫男
[考研] 085404 22408 309分求调剂 +10 lzmk 2026-04-14 11/550 2026-04-19 16:42 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[考研] 300求调剂 +12 橙a777 2026-04-15 12/600 2026-04-18 23:51 by 路病情
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[考研] 一志愿华中农业071010,320求调剂 +17 困困困困坤坤 2026-04-14 19/950 2026-04-17 20:08 by 关一盏灯cd
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
[基金申请] RY:中国产出的科学垃圾论文,绝对数量和比例都世界第一 +7 zju2000 2026-04-14 18/900 2026-04-16 11:36 by 欢乐颂叶蓁
信息提示
请填处理意见