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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » vasp 的mbj 泛函不收敛
111/2返回列表12下一页
查看: 4229  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

小木虫于

金虫 (正式写手)

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

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

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

VASP

» 猜你喜欢
1楼 2016-06-26 22:20:30
已阅   回复此楼   关注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手册)来提高收敛性。
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新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)
赞一下(1人) 回复此楼
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的回帖
相关版块跳转 我要订阅楼主 小木虫于 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭