| 查看: 2572 | 回复: 15 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
leongxf银虫 (初入文坛)
|
[交流]
[求助]vasp计算单个Ni原子能量不收敛问题
|
||
|
以下是计算使用的POSCAR，和INCAR POSCAR： Ni 1 10.0 0.0 0.0 0.0 10.0 0.0 0.0 0.0 10.0 ！ 周期取10A 1 Cartesian 5.00 5.00 5.00 ！ 原子在原胞中心 INCAR： System = fcc Ni ISTART = 0 ; ICHARG=2 PREC= n ISMEAR = 0; SIGMA = 0.01; ISPIN=1 IBRION=-1 NSW=0 NGX=70 NGY=70 NGZ=70 K点取 Automatic mesh 0 Monkhorst Pack 1 1 1 0 0 0 Iteration 60次，没办法收敛。请问大家应该怎样修改文件，新手请教。 [ Last edited by wuchenwf on 2009-6-22 at 22:11 ] |
回复此楼» 猜你喜欢
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有5人回复
-
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有5人回复
-
有机合成以后会不会被AI改变？做科研的虫友怎么看
已经有8人回复
-
求合成方法
已经有7人回复
-
现代”学阀”该如何界定
已经有12人回复
-
我的奶奶
已经有3人回复
leongxf
银虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 464.5
- 散金: 1
- 帖子: 30
- 在线: 6.2小时
- 虫号: 490352
- 注册: 2008-01-03
- 专业: 金属材料的合金相、相变及
采取各位意见，还是不能收敛
|
使用 建议1 ISPIN=2 AMIX=0.2 BMIX=0.0001 NELM=100 结果还是不能收敛 使用建议2 ISMEAR=-5 计算没有进行，显示下列错误：可否帮我分析下出错原因 VERY BAD NEWS! internal error in subroutine IBZKPT: Tetrahedron method fails for NKPT<4. NKPT = 1 FORTRAN STOP MPI Id 0 is using GM port 2, board 0 (MAC 0060dd4925fc). Received data from all 1 MPI processes. Sending mapping to MPI Id 0. Data sent to all processes. Reap remote processes: rsh comp65 -n kill -9 13191 2>/dev/null All remote MPI processes have exited. 使用建议3 K点选取10 10 10，但还是不能收敛 我以前自己加大K点，但一直不能收敛，所以还请大人帮我分析下该怎么计算。 |
5楼2008-06-29 17:59:00
leongxf
银虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 464.5
- 散金: 1
- 帖子: 30
- 在线: 6.2小时
- 虫号: 490352
- 注册: 2008-01-03
- 专业: 金属材料的合金相、相变及
计算结果
|
采取各位的建议，模拟了很多次，终于收敛了，也找到了不收敛的原因，分析一下，请各位多指教。 1. 刚开始的KEPOINTS，取1 1 1，是因为看到论坛上有贴的计算H原子能的例子，我test了这个列子，发现 1 1 1是收敛的，故Ni原子的也取一样的，后来不收敛，我把k点增加到10 10 10，还是不收敛，所以我分析可能是INCAR的参数设置问题。 根据楼上建议，分别尝试以下的INCAR（以下所有k点都为10 10 10 ） 1. 结果不收敛（使用PW91势） ISTART = 0 ; ICHARG=2 PREC= n ISMEAR = -5; SIGMA = 0.2; ISPIN=1 IBRION=-1 NSW=0 2。收敛 ISTART = 0 ; ICHARG=2 PREC= n ISMEAR = 0; SIGMA = 0.2; ISPIN=2 VOSKOWN=1 IBRION=-1 NSW=0 3。收敛（使用PBE势）***********因为有人说用此势比较好，故尝试 ISTART = 0 ; ICHARG=2 PREC= n ISMEAR = 0; SIGMA = 0.2; ISPIN=2 IBRION=-1 NSW=0 4。收敛（使用PW91势） ISTART = 0 ; ICHARG=2 PREC= n ISMEAR = 0; SIGMA = 0.2; ISPIN=1 IBRION=-1 NSW=0 综上所述，以及我另外做的其他尝试，排除了poscar的影响，以及势函数的影响，发现SIGMA = 0.2是收敛的主要因素。当然我也在以上条件3，4 INCAR中其他条件未变的情况下，仅改变SIGMA = 0.01，发现不收敛了。所以我想 SIGMA = 0.2是收敛的主要因素 2，3，4计算中，4没有计算磁矩，2，3考虑自旋，但计算结果发现 E2=-0.31eV E3=-0.54ev 2，3区别在于势函数不一样，所以采用不同的势函数计算了体材料Ni的单原子能分别为（我拟和了NI的晶格参数，所以结果都为取平衡晶格参数的值） PW91 -5.416 PBE -5.454 这样计算的NI的结合能是否为 PW91-E2=-5.106 PBE-E3=-4.914 试验上的数据参考《潜入原子方法理论及其在材料科学中的应用》为4.46ev，所以是不是我的计算还是有问题？还是我的结合能计算方法有问题？ |
12楼2008-07-03 11:49:37