应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 相同的vasp版本，不同机器上计算同一个体系，能量相差近10ev
101/1返回列表
查看: 3402  |  回复: 17
查看全部回帖 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

dfjk-123

新虫 (小有名气)

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

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-11-24 19:06:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dfjk-123

新虫 (小有名气)
刚检查了两个地方生成的OUTCAR文件，发现NBANDS数是不一样的，这可能是导致能量差的原因。但是怎么会差这么大呢？我的INCAR中并没有给定NBANDS数量啊。同样版本的vasp，给出的能带数应该一样才对啊。这是为什么呢？
赞一下 回复此楼
2楼2012-11-24 19:13:05
已阅   回复此楼   关注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的回帖

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

新虫 (小有名气)
引用回帖:
9楼: Originally posted by souledge at 2012-11-27 19:35:22
如果计算过程没有警告问题，总能（别管你取哪个）有差别都无所谓，只要生成能之类都没变化就可了。

warning给出的信息“Unless you are  performing a calculation for an insulator or semiconductor, without  unoccupied bands, you have included TOO FEW BANDS!! 我的体系是半导体。是不是意味着如果是半导体或者绝缘体的话，可以忽略这个警告？现在对于NBANDS数我还是用了VASP默认提供的数，不过都大于NELECT/2 + NIONS/2 (non-spinpolarized)了。是不是认为就是合理的了？求指导！谢谢
赞一下 回复此楼
11楼2012-11-28 08:47:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dfjk-123

新虫 (小有名气)
引用回帖:
13楼: Originally posted by 贺仪 at 2012-11-28 09:10:38
你先增加一些带试试吧 看有没有warming...

我就是看到这个waring才增加的NBANDS数量的，增加后发现能量有比较大的差异，就是前面帖子描述的现象。我测试了两个Nbands都有不同程度的能量变化，于是就想问下NBANDS的设置有没有什么标准。那个warning，如果是半导体或者导体的话，是不是就不需要理会，计算出来的相对能量是不是可信？谢谢
赞一下 回复此楼
14楼2012-11-29 14:14:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dfjk-123

新虫 (小有名气)
送鲜花一朵
引用回帖:
16楼: Originally posted by longwen36 at 2012-12-04 16:53:22
貌似也我也有这样的问题，
不同的机器NBANDS数不一样，所以算出来的结果也不一样。

INCAR里面没有给NBANDS，它好像会根据你算的CPU数来设定。貌似与并行的CPU数有关。
所以最好在INCAR里面最好固定NBANDS数， ...

嗯，有道理。谢谢你
回复此楼
17楼2012-12-05 15:22:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dfjk-123 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭