| 查看: 3969 | 回复: 6 | ||
[求助]
HSE06做自洽计算时提示增大NBANDS，怎么做？ 已有1人参与
|
|
用HSE06计算时，输出结果里出现了如下的错误，提示增加空带NBANDS。 我的INCAR里没有设置NBANDS。 用的是默认值， 我看了下OUTCAR里 NBANDS=144，根据手册NBANDS = NION/2 + ELECTRON/2 ,我计算的提醒，NIONS=48 ELECTRON =200 ，带入公式计算 NBANDS =124 ，为什么OUTCAR 里会显示144呢？ 再者，如果我要在NBANDS=144的基础上增大，我需要增大到多少较为合适呢？ 请各位大侠知道，金币奉上，感激不尽，在错误的后面我会把我的INCAR,KPOINTS 贴出来，请大家看下有啥子问题？ ----------------------------------------------------------------------------- | | | 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.' | | | ----------------------------------------------------------------------------- INCAR： ISTART = 1 ICHARG = 2 EDIFF = 0.00001 EDIFFG = -0.001 ENCUT = 400 LREAL = .FALSE. LWAVE = .TRUE. LCHARG = .TRUE. NELM = 200 NSW = 0 IBRION = -1 LMAXMIX = 4 ISMEAR = 0 SIGMA = 0.1 NSIM = 4 ISYM = 0 LHFCALC = .TRUE. HFSCREEN = 0.2 ALGO = D TIME = 0.4 ENCUTFOCK = 0 AEXX = 0.25 KPOINTS: Automatic generation 0 Monhkorst-Pack 4 4 4 0.0 0.0 0.0 |
回复此楼» 猜你喜欢
0703化学调剂
已经有12人回复
-
0703化学 305求调剂
已经有4人回复
-
0703化学调剂，求各位老师收留
已经有10人回复
-
271材料工程求调剂
已经有5人回复
-
281求调剂（0805）
已经有16人回复
-
304求调剂
已经有6人回复
-
材料工程专硕调剂
已经有6人回复
-
一志愿天大材料与化工（085600）总分338
已经有4人回复
-
085700资源与环境308求调剂
已经有3人回复
-
求材料调剂
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于（VASP）Fe100表面的SCF收敛
已经有19人回复
- GW准粒子计算会高估能隙？
已经有21人回复
- 相同的vasp版本，不同机器上计算同一个体系，能量相差近10ev
已经有17人回复
- 【已经完成】Optical properties and density functional perturbation theory 的翻译
已经有50人回复
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
zyqh: 金币+5 2014-11-30 12:52:58
感谢参与，应助指数 +1
zyqh: 金币+5 2014-11-30 12:52:58
|
默认的NBANDS并不是简单地按这个公式NION/2 + ELECTRON/2来计算的，它在OUTCAR里输出的值还取决其他计算参数的设置。 比如是否非共线磁性参数(LNONCOLLINEAR)、NPAR、ISPIN的值 你可以简单的读读main.F里有关NBANDS变量的赋值. .... NBANDS=MAX(NINT(INFO%NELECT+2)/2+MAX(T_INFO%NIONS/2,3),INT(0.6*INFO%NELECT))+NMAG IF (WDES%LNONCOLLINEAR) NBANDS = NBANDS*2 NBANDS=((NBANDS+NPAR-1)/NPAR)*NPAR ... 你可在NBANDS=144的基础在适当的增加一些（10~20），增加的太大的话，会增加计算量而已。 |
2楼2014-11-30 08:37:10
3楼2014-11-30 13:12:36
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
|
300设置的似乎太多了，会显著得增加计算量。除非你想得到导带很大能量范围内的信息，否则的话，建议减少它。 “10个节点，共120个核做的HSE06自洽计算，24个小时只跑了两步，我想请问候老师；这正常吗？” 这种问题，没有办法回答。尽量不要问这种问题。 “自洽计算不是重已有的结构出发直接算电荷和波函数吗？为何需要设置精度。” 这说明你根本没有理解“结构优化”和“自洽计算”的差别。 自洽计算---》 对某个结构进行电子波函数的自洽计算，即根据初始电荷密度（这个可以是已自洽过的，或未自洽的，或原子电荷密度），构造Kohn-Sham哈密顿量，求解本征波函数和本征值，根据本征波函数得到新的电荷密度，再重新构造Kohn-Sham哈密顿量，重复前两步直到所指定的物理量（总能或电密度等或波函数）的变化小于所设置的精度。看看有关Kohn-Sham方程自洽求解的简单介绍就了解。这个可以在很多资料中找到。 做自洽计算因此当然要设置收敛精度了。 结构优化是将对某个结构的原子位置或晶格基矢进行调整使得原子结构中原子所受的力小于某个值，晶格的应力小于某个值，总能变化小于某个值。在结构优化中不断调整结构，对其中的每个结构都会有电子波函数的自洽计算后，再计算原子所受的力的步骤。 “我想确认下HSE06计算能带步骤 DFT自洽 + HSE自洽+ 能带计算 对吗？” 这种确认从手册里就可以找到，手册上也明确建议了计算步骤，你顺着上面的就对了。 |
4楼2014-11-30 13:42:03
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
5楼2014-11-30 13:53:19
6楼2014-11-30 15:52:44
7楼2018-05-09 11:42:37
