| 查看: 1973 | 回复: 10 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
参数测试能量怎么都不稳定，求救！
|
||
|
最近在做一个体系的参数测试，搞了好久了，但是一直没有得出结果，虽然可以用文献上的参数，但是想自己验证一下，掌握方法。 脚本是这样 #!/bin/sh for i in 350 400 450 500 550 600 650 700 do cat > INCAR << ! SYSTEM = LMP ENCUT = $i ISTART = 0 ; ICHARG = 2 LDAU = .TRUE LDAUTYPE = 1 LDAUL = -1 2 2 -1 -1 LDAUU = 0 5.5 5.5 0 0 LDAUJ = 0 1 1 0 0 ISPIN = 2 MAGMOM = 4*0 1*7 2*-7 1*7 4*0 16*0 ISMEAR = 0; SIGMA = 0.05 NELM = 30 LMIXMIX =4 LREAL = Auto PREC = High ! echo "ENCUT = $i eV" ; mpirun -np 16 vasp5 E=$(grep "energy without entropy" OUTCAR | tail -1| awk '{printf "%12.6f \n", $5}') echo $i $E >>comment done KPOINTS 用的 3 6 9 结果是这样： 350 -200.25397 400 -199.72309 450 -199.49069 500 -199.42421 550 -199.44941 600 -199.49795 650 -199.53697 700 -199.57734 K点也取过很大的 趋势也是这样 sigma 也试了不同的值 能量始终都到不了稳定的值 另外我这个POSCAR是用的优化好的结构，我这个体系是反铁磁结构，所以加自旋加Ｕ，但我也试过不加的，能量趋势也是这样，然后用文献的值520 来测K点 已经增大到7 15 19了 能量还是一直降低 根本看不到能量变化小与1mev的希望啊…… 各位大神出来指点一二吧，我实在没招了…… 感激不尽啊！！！  |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有10人回复
-
Ei源刊怎么投
已经有3人回复
-
面上没中，邀请各位路过的虫友分析一下分数
已经有10人回复
-
广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！
已经有6人回复
-
两块石头
已经有5人回复
-
哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生
已经有3人回复
-
找导师
已经有8人回复
-
国社科又开始会评了，不知道这次命运如何
已经有14人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有14人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
送红花一朵 
11楼2013-11-17 18:40:15
2楼2013-11-10 22:53:45
3楼2013-11-11 15:07:28
4楼2013-11-11 20:03:45