vasp之mBJ计算能带方法交流 - 第一性原理 - 热点前沿问题讨论

版块导航: 正在加载中...

客户端APP下载

论文辅导

调剂小程序

登录注册

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

24小时热门版块排行榜

返回列表

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

第一性原理

» 猜你喜欢

球磨粉体时遇到了大的问题，请指教！已经有8人回复
江汉大学解明教授课题组招博士研究生/博士后已经有3人回复

1楼 2017-02-22 17:34:04

已阅关注TA 给TA发消息送TA红花 TA的回帖

回帖置顶 ( 共有1个 )

emilyoyang

木虫 (正式写手)

应助: 232 (大学生)
金币: 4340.4
红花: 42
帖子: 667
在线: 189.9小时
虫号: 1630291
注册: 2012-02-20
性别: MM
专业: 理论和计算化学

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
obaica: 金币+10 2017-02-23 14:56:28
obaica: 回帖置顶 2017-02-23 14:57:08

上面vasp官网已经说明的很详细了。你按照上面的instructions，设置参数计算就是了。

简单来说， VASP MBJ band gap计算只需添加 METAGGA = MBJ; LASPH = True 就行了
Prof Georg Kresse 说“The MBJLDA is just an XC potential, not an XC energy functional.
Thus E XC is taken from LSDA and the forces cannot be used from such calculations.” （言外之意就是，不要用它做结构优化，没意义，算算band gap就好了）

CMBJ = 1.3 （vasp默认值）如果你读了Dr F. Tran的文章（F. Tran and P. Blaha, Phys. Rev. Lett. 102, 226401 (2009)以及其他关于改方法的测试文章，自行google之。F.Tran 是MBJ的developer，MBJ在libxc中叫TB09）你会发现该值是可以调节的。通过调节该值可以得到跟实验值一样的band gap（很扯吧，哈哈。。。。）

VASP MBJ计算band gap最大的问题是收敛性。如果你在计算时候发现很难收敛，（那恭喜你中奖了！）解决办法就是 Prof Georg Kresse 说的
For fast convergence, it is suggested to use ALGO = N IMIX = 1 ! Kerker mixing. 如过还不收敛，那没招了！
收敛性比较困难是该方法的问题，Prof Georg Kresse也不想就MBJ做过多研究（不想开发新的加速收敛方法）。记得Dr F. Tran说过他专门针对MBJ对某些体系难收敛的问题开发了几个特殊的方法来加速收敛，但不幸的是这些方法只在Wien2k中有。

总而言之，VASP中使用MBJ有点“鸡肋”的感觉，遇到问题问Prof Georg Kresse他不搭理你。所以呢如果你想玩MBJ 还是去学wien2k吧。Dr F. Tran是很热心帮你解决所有MBJ的问题的。（他自己的“孩子”嘛哈哈）
Good Luck！