| 查看: 3620 | 回复: 17 | |||
dfjk-123新虫 (小有名气)
|
[求助]
相同的vasp版本，不同机器上计算同一个体系，能量相差近10ev
|
||
|
各位： 我在两个不同的LINUX集群下用vasp计算同一个体系，输入文件完全相同。两个集群安装的都是vasp5.2.2版本，结构优化完成后能量居然相差了近10ev，几何结构两个基本相同。大家说说这是怎么回事？ |
回复此楼» 猜你喜欢
材料复试求调剂
已经有8人回复
-
一志愿0703化学招61最终排名62化学求调剂
已经有17人回复
-
考研调剂
已经有5人回复
-
085600材料与化工301分求调剂院校
已经有33人回复
-
284求调剂
已经有9人回复
-
269求调剂
已经有5人回复
-
285求调剂
已经有6人回复
-
283求调剂
已经有17人回复
-
292求调剂
已经有8人回复
-
070300化学 求调剂
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 不同版本vasp 计算相同的超胞的能量 差了几十个meV 这算是合理或者是有问题？
已经有10人回复
- VASP静态自洽计算，能量升高7eV，然后局域磁矩消失，请问是什么问题？
已经有8人回复
- vasp计算单个氧原子和分子的能量
已经有20人回复
- 求助：vasp 计算固体相变，能量。
已经有3人回复
- VASP计算反铁磁和铁磁态能量差高达1eV以上，是否可信？
已经有10人回复
- 求助：vasp 计算氢原子能量是-1.16ev 而不是 -13.6ev，why ?
已经有18人回复
- 【求助】vasp静态计算能量值错误
已经有21人回复
- 【其他】vasp结合phonon计算材料不同温度下能量的问题
已经有9人回复
- 【转帖】用VASP计算氢分子的能量
已经有8人回复
【答案】应助回帖
★
WDD880227: 金币+1, 多谢指教 2012-11-28 14:31:47
WDD880227: 金币+1, 多谢指教 2012-11-28 14:31:47
|
‘我后来试了NBANDS数不变，加大了K点，能量上升了。减少K点的话，warning就没有了。‘ 这明显是K点不够，warming没有了是因为在某些K点的最高带上已经没有占据了，但这个没有是错误的提示。 对于K点的够不够你测试的时候不知道的参数就用默认的，比如说NBANDS。测试的时候增加K点的密度，直到能量达到一个稳定的值（误差在0.几个eV内，越小越好）。 总的来讲我认为你的NBANDS数量不够，NELECT/2 + NIONS/2 (non-spinpolarized)这个没错，但你的电子数是不是数对了呢？ |
» 本帖已获得的红花（最新10朵）
dfjk-123
8楼2012-11-27 19:14:36
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
9楼2012-11-27 19:35:22
longwen36
铁杆木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 6165.9
- 红花: 16
- 帖子: 344
- 在线: 208.4小时
- 虫号: 766090
- 注册: 2009-05-08
- 性别: GG
- 专业: 功能陶瓷
16楼2012-12-04 16:53:22
dfjk-123
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 29.9
- 散金: 757
- 帖子: 119
- 在线: 24.5小时
- 虫号: 665485
- 注册: 2008-12-01
- 性别: MM
- 专业: 理论和计算化学
2楼2012-11-24 19:13:05
3楼2012-11-25 10:04:22
dfjk-123
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 29.9
- 散金: 757
- 帖子: 119
- 在线: 24.5小时
- 虫号: 665485
- 注册: 2008-12-01
- 性别: MM
- 专业: 理论和计算化学
4楼2012-11-25 12:08:56
5楼2012-11-25 20:37:33
dfjk-123
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 29.9
- 散金: 757
- 帖子: 119
- 在线: 24.5小时
- 虫号: 665485
- 注册: 2008-12-01
- 性别: MM
- 专业: 理论和计算化学
|
好的，谢谢。还有个问题想请教下，我在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
dfjk-123
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 29.9
- 散金: 757
- 帖子: 119
- 在线: 24.5小时
- 虫号: 665485
- 注册: 2008-12-01
- 性别: MM
- 专业: 理论和计算化学
7楼2012-11-27 15:58:54
dfjk-123
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 29.9
- 散金: 757
- 帖子: 119
- 在线: 24.5小时
- 虫号: 665485
- 注册: 2008-12-01
- 性别: MM
- 专业: 理论和计算化学
10楼2012-11-28 08:39:37
100