24小时热门版块排行榜

返回列表

小小鑫000

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 6.6
帖子: 267
在线: 111.8小时
虫号: 2648268
注册: 2013-09-11
性别: MM
专业: 理论和计算化学

[求助] 关于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 ]

回复此楼

好好研究，天天向上！

1楼2013-11-23 21:48:46

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

jpchou

专家顾问 (著名写手)

专家经验: +20
1ST强帖: 5
应助: 231 (大学生)
金币: 2342
散金: 4367
红花: 103
帖子: 1061
在线: 757.3小时
虫号: 585915
注册: 2008-08-08
性别: GG
专业: 凝聚态物性 II ：电子结构
管辖: 第一性原理

420 eV 為什麼不能再大上去?

赞一下(1人)

回复此楼

2楼2013-12-03 17:21:20

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

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左右的啊～

赞一下(1人)

回复此楼

思想重于技巧，内涵重于表象

3楼2013-12-03 23:54:48

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