| 查看: 2414 | 回复: 5 | ||||
[求助]
关于VASP中对ENCUT收敛性测试的疑问 已有1人参与
|
|
各位高手大家好： 有个问题一直困扰我很多天了，是这样的。我近期用vasp对Co进行模拟计算。首先要确定Co的ENCUT和KPOINT大小，但是文献上没告诉，所以得自己先进行收敛性测试。资料上说先要假设一个ENCUT值，然后进行KPOINT的收敛性测试，疑问出来了，那就是假设的ENCUT应该取多大？这是第一个疑问。 我当时就取了PONTCAR中的ENMAX=270，接着进行了K点的收敛性测试，最后确定为8x8x8，然后用这个KPOINT再进行ENCUT的收敛性测试，得到这个结果： 270 -27.073083 300 -27.112801 330 -27.106534 360 -27.089872 390 -27.080460 420 -27.077593 没有出现总能变化在0.001，而且420这已经是一个大的不能再大的值了，所以420后面的值更不可能收敛了，我的脚本也应该没问题呀！ #!/bin/sh rm WAVECAR for i in 15 270 300 330 360 390 420 do cat > INCAR <<! SYSTEM = Co-Diamond ENCUT = $i ISTART = 0 ; ICHARG = 2 ISMEAR = -5 PREC = Accurate ! echo "ENCUT = $i eV" ; time vasp E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'` echo $i $E >>comment done 到底哪儿出问题了,求高手解释呀！不胜感激！ [ Last edited by 小小鑫000 on 2013-11-27 at 10:32 ] |
回复此楼» 本帖@通知
» 猜你喜欢
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有3人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有21人回复
-
能否退出参与的面上项目解除限项
已经有17人回复
-
估计是周四
已经有4人回复
-
我面上完蛋了
已经有9人回复
-
2026国自然函评费到账
已经有18人回复
-
建议基金发布提前给出明确的时间点
已经有15人回复
-
范进中举一文的中心思想
已经有7人回复
-
朋友圈看到的
已经有9人回复
-
让我中一个面上吧！
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp阶段能测试求助，谢谢啊！急求。
已经有15人回复
- 对于收敛性检测
已经有24人回复
- 关于VASP计算收敛的判断
已经有25人回复
- VASP静态自洽计算收敛特慢！
已经有31人回复
- VASP NEB 不收敛问题
已经有25人回复
- vasp中，做截断能测试时，测试收敛的标准
已经有10人回复
- 很奇怪的问题，vasp计算ENCUT设大就不让计算，高手帮忙啊！
已经有26人回复
- 如何做vasp掺杂
已经有15人回复
- 关于vasp收敛和优化
已经有24人回复
- vasp优化难收敛，请大家给我看看
已经有5人回复
- VASP弛豫收敛太慢
已经有18人回复
- vasp计算中，如何选择正确的K点
已经有7人回复
- 求助：vasp结构优化不能收敛
已经有18人回复
- ENCUT测试脚本
已经有17人回复
- VASP计算时的几个细节设置问题，希望有经验的人过来指点一下！
已经有23人回复
- VASP 截断能测试 急求
已经有8人回复
- 【求助】VASP收敛问题
已经有15人回复
- 【求助成功】vasp测试K点收敛性的脚本问题
已经有20人回复
- 【原创】写了个集群下应用测试K点或截断能的脚本
已经有18人回复
- 【求助】收敛测试怎么做？？？
已经有17人回复
- 【求助】vasp不同磁结构最后的收敛问题
已经有6人回复
2楼2013-12-03 17:21:20
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
|
虽然一般可能更多是先假设K-mesh进行ECUT的收敛测试，不过反过来也没啥原则性错误。 不过问题就在于，如果要假设ECUT，那就必然要假设个“一定能够满足收敛”的ECUT，否则得到的K-mesh又能有什么意义？ 所以，这就需要把ECUT设置得比一般别人文献报告的还要略高，才比较安全。如果我是LZ，可能ECUT至少要从500eV或者更高来假设，然后对K-mesh进行测试。之后在固定合理的k-mesh的情况下，回头对ecut进行假设。 话说，为什么"420这已经是一个大得不能再大的值了"呢？很多用VASP计算的文献，给出的cut-off energy都是500eV左右的啊～ |
3楼2013-12-03 23:54:48
4楼2013-12-04 16:57:24
5楼2013-12-04 16:59:05
6楼2018-09-18 21:48:46
20