24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

专家经验: +224
1ST强帖: 2
应助: 47 (小学生)
金币: 3488.8
散金: 14
红花: 223
帖子: 2477
在线: 1253.7小时
虫号: 1369423
注册: 2011-08-16
性别: GG
专业: 金属材料的磨损与磨蚀
管辖: 第一性原理

[求助] 关于（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却解决不了，实在不知到为什么。希望各位前辈能给予帮助。

回复此楼

» 猜你喜欢

垃圾破二本职称评审标准已经有18人回复
职称评审没过，求安慰已经有53人回复
毕业后当辅导员了，天天各种学生超烦已经有5人回复
26申博自荐已经有3人回复
A期刊撤稿已经有4人回复
EST投稿状态问题已经有7人回复

不一定挂在论坛，计算问题问题欢迎留言。

1楼 2014-04-15 15:44:28

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

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

专家经验: +224
1ST强帖: 2
应助: 47 (小学生)
金币: 3488.8
散金: 14
红花: 223
帖子: 2477
在线: 1253.7小时
虫号: 1369423
注册: 2011-08-16
性别: GG
专业: 金属材料的磨损与磨蚀
管辖: 第一性原理

引用回帖:

17楼: Originally posted by kjsun at 2014-04-16 17:18:57
这种情况对于磁性体系很常见。
磁性体系由于存在多种简并或者非简并态，导致电子有N多种选择，所以会造成收敛困难问题。简单的理解，就是他有唯一选择是（非磁性体系），冲着目标去就行了。如果它有1万种选择，它就 ...

正解，从scf迭代的观测中可以看出来能量值在几个值来回飘，确实是这样的，磁性体系的通病。
我先做个测试。

谢谢。