| 查看: 2082 | 回复: 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 ] |
回复此楼» 猜你喜欢
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有5人回复
-
博士读完未来一定会好吗
已经有29人回复
-
博士申请都是内定的吗？
已经有5人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有12人回复
-
投稿精细化工
已经有4人回复
-
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
huangyc
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 2242.1
- 散金: 12
- 红花: 2
- 帖子: 411
- 在线: 312.5小时
- 虫号: 624241
- 注册: 2008-10-12
- 专业: 理论和计算化学
15楼2009-02-19 16:16:22
★ ★ ★
lei0736(金币+2,VIP+0):谢谢持续关注
leongxf(金币+1,VIP+0):多谢
lei0736(金币+2,VIP+0):谢谢持续关注
leongxf(金币+1,VIP+0):多谢
|
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. 出现这个错误是因为你用的是超晶胞，ISMEAR=-5 是不能用的 很奇怪啊，为什么你的不收敛呢，我也做过Ni的，好像很容易的 你用的是什么势呢， 把超晶胞取大点吧，降低对称性，如 10X15X20 试试 |
6楼2008-06-29 20:03:43
beefly
专家顾问 (职业作家)
地沟油冶炼专家
-
专家经验: +457 - 应助: 86 (初中生)
- 金币: 12741.4
- 散金: 27307
- 红花: 300
- 帖子: 4663
- 在线: 925.3小时
- 虫号: 408372
- 注册: 2007-06-21
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 计算模拟
8楼2008-06-30 10:28:18
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
2楼2008-06-27 20:28:20
3楼2008-06-27 20:31:03
sw-2384
木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 3614.3
- 散金: 37
- 红花: 4
- 帖子: 386
- 在线: 325.9小时
- 虫号: 263155
- 注册: 2006-07-01
- 性别: MM
- 专业: 高分子科学
4楼2008-06-29 08:52:47
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
天地沙鸥
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 4592.2
- 红花: 1
- 帖子: 148
- 在线: 100.5小时
- 虫号: 555358
- 注册: 2008-05-09
- 性别: GG
- 专业: 高分子合成化学
7楼2008-06-30 08:17:54
9楼2008-07-01 17:54:40
10楼2008-07-02 10:05:00