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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于（VASP）Fe100表面的SCF收敛
61/1返回列表
查看: 2635  |  回复: 19
查看全部回帖 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

[求助] 关于（VASP）Fe100表面的SCF收敛 已有5人参与

从我做计算开始Fe100表面就是个难题，目前这些常见的程序就没有几个能够正儿八经的把这个体系算出来。
想请教一下各位些问题，VASP采用怎么样的参数能够把Fe100计算收敛？

我自己的计算采用的是PAW-PBE的势，GGA=RP，真空层=25A，表面大小1x1x1，高度是5个原子，Monkhorst-Pack K=10x10x1 Nbands=30 PREC=High
目前采用ALGO=48、38、53都不成，其实functional换过不少，另外，对于MIX项也调整过。本人新手，不知各位前辈看看哪里的设定有错误，以及有没有解决方案。

之前本人用MS做计算，Dmol3彻底没办法（LCAO处理金属体系我也不说什么）。对于CASTEP而言PBE也是很难收敛。由于去年在做毕业论文看了一下Hammer.B的文章，指出RPBE对表面计算收敛性比较好，所以CASTEP采用RPBE+Ultrasoft pseudopotential确实能够收敛。

既然CASTEP能够处理这个问题，VASP却解决不了，实在不知到为什么。希望各位前辈能给予帮助。
回复此楼

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

VASP 第一性原理和电化学

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
不一定挂在论坛，计算问题问题欢迎留言。
1楼 2014-04-15 15:44:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hn19870519

木虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
Fe001面还算好收敛吧。。。一般电子步收敛的不好是结构没弛豫好吧。
赞一下 回复此楼
6楼2014-04-16 10:02:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hn19870519

木虫 (正式写手)

【答案】应助回帖

另外 +U和收敛不收敛没啥具体关系吧
赞一下 回复此楼
7楼2014-04-16 10:02:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hn19870519

木虫 (正式写手)

【答案】应助回帖

楼主的意思是不是设定5层，然后只做自洽？这样不收敛不奇怪。因为你的结构没有跑稳定，电子步不收敛也正常！
赞一下 回复此楼
11楼2014-04-16 15:23:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hn19870519

木虫 (正式写手)

【答案】应助回帖

初始结构不是最优的结构，开始的几步离子步不收敛很正常，我一般设置NELM=200跑，五层的Fe001面比较对称，肯定可以跑收敛的，我之前也算过类似的体系。当然就像楼主说的，最后还是要看算出来的力是不是收敛了。
赞一下 回复此楼
13楼2014-04-16 15:46:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hn19870519

木虫 (正式写手)

【答案】应助回帖

另外，我当初算的时候 用PBE就可以很好的收敛了。就用38来跑，也没有设置别的特别的参数就可以很好收敛。
赞一下 回复此楼
15楼2014-04-16 16:18:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 卡开发发 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭