| 查看: 2142 | 回复: 12 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
nkleof银虫 (正式写手)
|
[交流]
【求助】晶格优化问题 已有7人参与
|
||
|
我用vasp优化一个体系，所得到的结果8.36A远小于实验值8.51A。 所采用的赝势是PAWPBE的，过渡金属元素Fe和Mn分别试过标准势和pv势。采用的优化方法分别试过先优化原胞体积（ISIF=7），后驰豫离子位置（ISIF=2）；以及全优化（ISIF=3）。驰豫方式分别试过CG法（IBRION=2）和准牛顿法（IBRION=1）。为了消除pulay stress，每种方法都以得到的CONTCAR作为POSCAR反复优化了3次。 但是以上所有方法所得到的结果都是8.36附近，远小于实验值，也小于其它文献所给出的GGA计算值8.49。 请经验丰富的虫友帮忙想想可能是在什么地方出问题了？ |
回复此楼» 猜你喜欢
极限拉扯，如何抉择
已经有14人回复
-
申请科研项目纯靠实力能成功吗？
已经有6人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有12人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
nkleof
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2093.2
- 散金: 649
- 红花: 5
- 帖子: 813
- 在线: 165.9小时
- 虫号: 152581
- 注册: 2006-01-02
- 专业: 凝聚态物性 II ：电子结构
11楼2010-04-20 16:52:34
gfei2008
木虫 (小有名气)
- 应助: 6 (幼儿园)
- 金币: 1716.7
- 红花: 1
- 帖子: 176
- 在线: 116.2小时
- 虫号: 544222
- 注册: 2008-04-12
- 专业: 凝聚态物性 II ：电子结构
2楼2010-04-19 21:11:59
nkleof
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2093.2
- 散金: 649
- 红花: 5
- 帖子: 813
- 在线: 165.9小时
- 虫号: 152581
- 注册: 2006-01-02
- 专业: 凝聚态物性 II ：电子结构
|
我用的k网格疏密程度是k点之间距离0.09 /A，PREC=Accurate因此ENCUT是赝势内默认值的1.3倍，所以这两个优化条件我觉得精度还可以，如果有问题请经验丰富的高人再给点建议。 拟合的优化正打算做，但是我觉得vasp的自动优化功能误差应该不会太大，因此用ISIF=2或3的设置算出来的值和静态计算拟合优化出来的值应该不会有太大差异，是不是有可能我的什么设置出错？ 如下是我做优化用的INCAR PREC = Accurate LREAL = A IALGO = 48 NPAR = 8 NSIM = 8 NELMIN = 4 NSW = 200 LORBIT = 12 EDIFF = 1E-5 IBRION = 2 ISIF = 3 ISPIN = 2 MAGMOM = 8*-6 16*6 32*0 [ Last edited by nkleof on 2010-4-19 at 21:27 ] |
3楼2010-04-19 21:23:24
4楼2010-04-20 10:09:34
200