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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » vasp 的mbj 泛函不收敛
111/2返回列表12下一页
查看: 4323  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 情人节自我反思：在爱情中有过遗憾吗？ +5 瞬息宇宙 2026-02-15 6/300 2026-02-18 12:51 by 月下雪林
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 11:09 by lqtl9djx19
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:54 by lqtl9djx19
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:39 by lqtl9djx19
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:53 by lqtl9djx19
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:38 by lqtl9djx19
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 07:55 by lotyj5cz79
[基金申请] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:40 by lotyj5cz79
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:38 by lotyj5cz79
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:23 by lotyj5cz79
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:08 by lotyj5cz79
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 3/150 2026-02-18 06:53 by lotyj5cz79
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 00:40 by tk2gfblvuz
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 00:23 by tk2gfblvuz
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-17 23:40 by tk2gfblvuz
[基金申请] 基金正文30页指的是报告正文还是整个申请书 +3 successhe 2026-02-16 4/200 2026-02-17 20:56 by successhe
[基金申请] 今年春晚有几个节目很不错，点赞！ +5 瞬息宇宙 2026-02-16 6/300 2026-02-17 12:49 by jymy19840415
[微米和纳米] 球磨粉体时遇到了大的问题，请指教！ 10+3 6sbiam 2026-02-12 15/750 2026-02-16 15:03 by tgzxzqj
[基金申请] 过年走亲戚时感受到了所开私家车的鄙视链 +3 瞬息宇宙 2026-02-15 5/250 2026-02-16 14:23 by aspect3000
[硕博家园] 江汉大学解明教授课题组招博士研究生/博士后 +3 cleverlyy 2026-02-12 3/150 2026-02-12 21:02 by qsdf1
信息提示
关闭
请填处理意见
关闭