24小时热门版块排行榜    

Znn3bq.jpeg
查看: 4391  |  回复: 10

小木虫于

金虫 (正式写手)

[求助] vasp 的mbj 泛函不收敛 已有1人参与

大家好,我最近用vasp的mbj方法算体系的能带,怎么都不收敛。如果ialgo设置成normal,能量增加了又减少,而且总能一直都是正的。如果ialgo设置成all,总能看着向收敛,但是一直增加,而且总能也一直是正的。大家谁做过mbj的计算,该怎么处理呢?多谢!
回复此楼

» 收录本帖的淘帖专辑推荐

VASP

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

emilyoyang

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
小木虫于: 金币+20, ★★★★★最佳答案 2016-06-27 23:40:25
你先确定一下 你的POTCAR对不对,要用新的pseudopotential,里面kinetic energy
老的不行。

我的经验是  MBJ的收敛性跟元素有关,有的体系很容易收敛,有的体系非常难。难收敛的体系 尝试一下mixing参数。还不收敛,只能放弃了。

某些体系,用vasp 算MBJ 不收敛是很正常的。即使用wien2k算  也需要一定的技巧。Fabien Tran(MBJ算法的developer)跟我说过  他后来又专门针对这个难收敛 搞个几个mixer (详见wien2k手册)来提高收敛性。

» 本帖已获得的红花(最新10朵)

3楼2016-06-27 09:09:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emilyoyang

木虫 (正式写手)

引用回帖:
4楼: Originally posted by 小木虫于 at 2016-06-27 14:04:46
多谢 我想如果POTCAR中没有kinetic energy,应该报错了吧。而且体系总能一直为正值,虽说mbj的总能是没用的,但是一直为正,还是有些犯嘀咕。我用ALGO=All的时候,看着像收敛的感觉,但是能量是一直增加的,只是增 ...

MBJ一开始算的时候  可以把  NELM 设置的大点  e.g.: 2000

MBJ band gap没有说算出来一定跟实验很接近。它算出来的band gap 跟c(CMBJ refers to vasp manual) 有关。可以调节一下c 来获得跟实验compare的值(有点HSE的味道,哈哈)。建议你看看这篇文章PHYSICAL REVIEW B 83, 195134 (2011)
8楼2016-06-28 01:55:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

蛋蛋小童鞋

版主 (著名写手)

vasp也可以加mbj吗,不知道,只知道wien2k可以

发自小木虫Android客户端
愿不负最苦的时光
2楼2016-06-27 00:17:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫于

金虫 (正式写手)

引用回帖:
3楼: Originally posted by emilyoyang at 2016-06-27 09:09:01
你先确定一下 你的POTCAR对不对,要用新的pseudopotential,里面kinetic energy
老的不行。

我的经验是  MBJ的收敛性跟元素有关,有的体系很容易收敛,有的体系非常难。难收敛的体系 尝试一下mixing参数。还不收 ...

多谢 我想如果POTCAR中没有kinetic energy,应该报错了吧。而且体系总能一直为正值,虽说mbj的总能是没用的,但是一直为正,还是有些犯嘀咕。我用ALGO=All的时候,看着像收敛的感觉,但是能量是一直增加的,只是增加的逐渐减小,从10^6 逐渐到 10^0。 最后超出NELM重算了。我的体系中包含Yb Al Sb,但是Yb用的没有4f电子的赝势。
4楼2016-06-27 14:04:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jx在路上

金虫 (著名写手)

请问,MBJ算能带和其它方法有啥不同。

发自小木虫Android客户端
5楼2016-06-27 23:29:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

蛋蛋小童鞋

版主 (著名写手)

引用回帖:
5楼: Originally posted by jx在路上 at 2016-06-27 23:29:37
请问,MBJ算能带和其它方法有啥不同。

好像是因为比较准确

发自小木虫Android客户端
愿不负最苦的时光
6楼2016-06-27 23:38:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫于

金虫 (正式写手)

送红花一朵
引用回帖:
3楼: Originally posted by emilyoyang at 2016-06-27 09:09:01
你先确定一下 你的POTCAR对不对,要用新的pseudopotential,里面kinetic energy
老的不行。

我的经验是  MBJ的收敛性跟元素有关,有的体系很容易收敛,有的体系非常难。难收敛的体系 尝试一下mixing参数。还不收 ...

确实是赝势的问题,我换了最新的赝势,收敛容易多了. 但是结果还是对不上实验.我在从其他角度入手吧.
7楼2016-06-27 23:40:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Physicser

木虫 (小有名气)

精度设置低一点。
MBJ的目的是对半导体的能隙估值要精确点,计算量比HSE小一点
9楼2016-07-01 09:36:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

周廷伟周廷伟

新虫 (初入文坛)

引用回帖:
6楼: Originally posted by 蛋蛋小童鞋 at 2016-06-27 23:38:19
好像是因为比较准确
...

mbj计算AB异质结能带的问题!现在我想对A用mbj,就是普通的PBE计算。在一个INCAR如何实现?
Iamme!
10楼2018-09-25 17:21:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小木虫于 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿211,化学310分,本科重点双非,求调剂 +23 努力奋斗112 2026-04-08 23/1150 2026-04-10 23:29 by 314126402
[考研] 求调剂288 +6 ioodiiij 2026-04-10 8/400 2026-04-10 21:07 by zhouxiaoyu
[考研] 求调剂 +16 MAX怅惘 2026-04-09 18/900 2026-04-10 20:24 by 小鬼喔~~
[考研] 0858求调剂 5+5 Gky09300550, 2026-04-10 8/400 2026-04-10 19:13 by chemisry
[考研] 本9 一志愿西工大085601 324求调剂 +5 wysyjs25 2026-04-10 5/250 2026-04-10 16:57 by luoyongfeng
[考研] 一志愿西交机械专硕求调剂 +8 求上岸的小王 2026-04-10 8/400 2026-04-10 15:09 by hemengdong
[考研] 085402通信工程调剂,有4项学科竞赛国奖(电赛国二),硕士研究生调剂自荐信。 +4 m永o不v言o弃m 2026-04-09 4/200 2026-04-10 11:30 by asy1wn
[考研] 347材料专硕求调剂 +19 zj8215216 2026-04-06 19/950 2026-04-10 09:36 by 690616278
[考研] 293求调剂 +5 勇远库爱314 2026-04-08 5/250 2026-04-10 08:46 by vgtyfty
[考研] 085500求调剂材料 +8 易11122 2026-04-09 8/400 2026-04-09 23:15 by parmtree
[考研] 269求调剂 +7 跪求收留。 2026-04-04 7/350 2026-04-09 19:06 by 探123
[考研] 305求调剂 +4 77Qi 2026-04-07 4/200 2026-04-09 17:27 by wp06
[考研] 070300化学279求调剂 +17 哈哈哈^_^ 2026-04-08 18/900 2026-04-09 10:49 by 三七七想上岸
[考研] 325 调剂 +6 QQ小虾 2026-04-07 6/300 2026-04-07 15:17 by Ccclqqq
[考研] 求调剂 +4 电气小神童 2026-04-04 6/300 2026-04-07 00:14 by guanxin1001
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +11 大火山小火山 2026-04-05 11/550 2026-04-06 22:55 by yunlongyang
[考研] 一志愿武汉理工大学080200机械工程308分,求调剂 +4 终不似从前 2026-04-05 4/200 2026-04-06 11:46 by 考研学校招点人
[考研] 283求调剂 +5 baiiyu 2026-04-05 6/300 2026-04-05 20:35 by 啵啵啵0119
[考研] 调剂 +3 李广火 2026-04-05 3/150 2026-04-05 18:57 by 蓝云思雨
[考研] 考研调剂生寻找导师 +3 顾瞻考研啊 2026-04-05 3/150 2026-04-05 18:18 by 啵啵啵0119
信息提示
请填处理意见