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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » HSE06做自洽计算时提示增大NBANDS，怎么做？
71/1返回列表
查看: 3822  |  回复: 6
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

zyqh

铁杆木虫 (著名写手)

[求助] 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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
雨后方能见彩虹
1楼2014-11-29 19:53:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +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
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyqh

铁杆木虫 (著名写手)
引用回帖:
2楼: Originally posted by valenhou001 at 2014-11-30 08:37:10
默认的NBANDS并不是简单地按这个公式NION/2 + ELECTRON/2来计算的，它在OUTCAR里输出的值还取决其他计算参数的设置。
比如是否非共线磁性参数(LNONCOLLINEAR)、NPAR、ISPIN的值

你可以简单的读读main.F里有关NB ...

万分荣幸得到候老师的回答，感谢候老师的回答。候老师，我有两个问题请教您：
1,我后来重新算了，NBANDS 设置为默认值的2.5倍，默认144，设置成了300（根据网上经验帖子），发现计算机运行很慢了。我是在长沙超算上算的。10个节点，共120个核做的HSE06自洽计算，24个小时只跑了两步，我想请问候老师；这正常吗？
2，我看别人HSE06自洽的INCAR中，设置了力的收敛精度和能量的收敛精度，自洽计算不是重已有的结构出发直接算电荷和波函数吗？为何需要设置精度。
3，候老师，我想确认下HSE06计算能带步骤 DFT自洽 + HSE自洽+ 能带计算 对吗？
谢谢您！
赞一下 回复此楼
雨后方能见彩虹
3楼2014-11-30 13:12:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
引用回帖:
3楼: Originally posted by zyqh at 2014-11-30 14:12:36
万分荣幸得到候老师的回答，感谢候老师的回答。候老师，我有两个问题请教您：
1,我后来重新算了，NBANDS 设置为默认值的2.5倍，默认144，设置成了300（根据网上经验帖子），发现计算机运行很慢了。我是在长沙超算 ...

300设置的似乎太多了，会显著得增加计算量。除非你想得到导带很大能量范围内的信息，否则的话，建议减少它。
“10个节点，共120个核做的HSE06自洽计算，24个小时只跑了两步，我想请问候老师；这正常吗？”
这种问题，没有办法回答。尽量不要问这种问题。

“自洽计算不是重已有的结构出发直接算电荷和波函数吗？为何需要设置精度。”
这说明你根本没有理解“结构优化”和“自洽计算”的差别。
自洽计算---》 对某个结构进行电子波函数的自洽计算，即根据初始电荷密度（这个可以是已自洽过的，或未自洽的，或原子电荷密度），构造Kohn-Sham哈密顿量，求解本征波函数和本征值，根据本征波函数得到新的电荷密度，再重新构造Kohn-Sham哈密顿量，重复前两步直到所指定的物理量（总能或电密度等或波函数）的变化小于所设置的精度。看看有关Kohn-Sham方程自洽求解的简单介绍就了解。这个可以在很多资料中找到。

做自洽计算因此当然要设置收敛精度了。

结构优化是将对某个结构的原子位置或晶格基矢进行调整使得原子结构中原子所受的力小于某个值，晶格的应力小于某个值，总能变化小于某个值。在结构优化中不断调整结构，对其中的每个结构都会有电子波函数的自洽计算后，再计算原子所受的力的步骤。

“我想确认下HSE06计算能带步骤 DFT自洽 + HSE自洽+ 能带计算 对吗？”
这种确认从手册里就可以找到，手册上也明确建议了计算步骤，你顺着上面的就对了。
赞一下(1人) 回复此楼
4楼2014-11-30 13:42:03
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
HSE06的计算当然会很慢，耐心！
赞一下 回复此楼
5楼2014-11-30 13:53:19
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyqh

铁杆木虫 (著名写手)
引用回帖:
4楼: Originally posted by valenhou001 at 2014-11-30 13:42:03
300设置的似乎太多了，会显著得增加计算量。除非你想得到导带很大能量范围内的信息，否则的话，建议减少它。
“10个节点，共120个核做的HSE06自洽计算，24个小时只跑了两步，我想请问候老师；这正常吗？”
这种问 ...

谢谢候老师了，很受教。您的入门手册也给予了我这个初级菜鸟很多的影响。谢谢您
赞一下 回复此楼
雨后方能见彩虹
6楼2014-11-30 15:52:44
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

2120160150

新虫 (初入文坛)
引用回帖:
6楼: Originally posted by zyqh at 2014-11-30 15:52:44
谢谢候老师了，很受教。您的入门手册也给予了我这个初级菜鸟很多的影响。谢谢您...

你好，我刚刚入门，能分享一下您的手册吗？
赞一下 回复此楼
7楼2018-05-09 11:42:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zyqh 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭