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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 相同的vasp版本，不同机器上计算同一个体系，能量相差近10ev
181/2返回列表12下一页
查看: 3322  |  回复: 17
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

dfjk-123

新虫 (小有名气)

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

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2012-11-24 19:06:11
已阅   回复此楼   关注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)这个没错，但你的电子数是不是数对了呢？
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

dfjk-123
嘟啦啦嘟啦啦嘟
8楼2012-11-27 19:14:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
★
WDD880227: 金币+1, 谢谢提示 2012-11-28 14:31:59
如果计算过程没有警告问题，总能（别管你取哪个）有差别都无所谓，只要生成能之类都没变化就可了。
赞一下(1人) 回复此楼
思想重于技巧，内涵重于表象
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数，就好像不会有这样的问题。
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

dfjk-123
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 的主题更新
181/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭