24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

GTF1902

金虫 (小有名气)

应助: 2 (幼儿园)
金币: 1045.4
散金: 20
帖子: 188
在线: 63.6小时
虫号: 1720390
注册: 2012-03-27
专业: 凝聚态物性I:结构、力学和

[求助] 【重新开帖】KPOINTS的测试脚本问题

我用vasp优化结构先测试K点时出了点问题，请大家看完下面是我的一个优化结构vasp文件 INCAR ：
System = PCCOF ISIF = 2 ISMEAR = 0 ISTART = 0 ICHARG = 2 SIGMA = 0.1 NSW = 0 PREC = Medium ENCUT = 500 GGA = 91 VOSKOWN = 1 IBRION = 2 EDIFF = 1E-4 EDIFFG = -0.03 LREAL = Auto IALGO = 48 LALANE = .TRUE NPAR = 6 LSCALU = .FALSE. NSIM = 4 RUN_K： #!/bin/sh rm -f WAVECAR for i in 1 2 3 4 5 6 7 do cat>KPOINTS<<! Automatic generation 0 Monkhorst-Pack $i $i $i 0.0 0.0 0.0 ! echo"k mesh = $i x $i x $i"; time /share/home/sci_028/sourcecode/vasp.5.2/vasp>out.log E=`grep"TOTEN"OUTCAR |tail -1| awk '{printf"%12.6f\n", $5 }'` KP=`grep"irreducible"OUTCAR |tail -1| awk '{printf"%5i\n", $2 }'` echo $i $KP $E>>comment done
还没算完但是 comment文件非常怪异：大侠请看
1 1 -569.842023 1 1 -569.842023 1 1 -569.842023 1 1 -569.842023 1 1 -569.842023 1 1 -569.842023 2 1 -600.814240 2 1 -600.814240 2 1 -600.814240 2 1 -600.814240 2 1 -600.814240 2 1 -600.814240 3 8 -599.488993 3 8 -599.488993 3 8 -599.488993 3 8 -599.488993 3 8 -599.488993 3 8 -599.488993 4 8 -599.575041 4 8 -599.575041 4 8 -599.575041 4 8 -599.575041 4 8 -599.575041 4 8 -599.575041 5 27 -599.519665 5 27 -599.519665 5 27 -599.519665 5 27 -599.519665 5 27 -599.519665 5 27 -599.519665
刚测试完第五个点
问题是为啥同样的数据要出现6次啊是不是我的测试脚本文件出问题了这是我第二遍算了第一遍测试4个K点也是同样的数据在comment文件中显示六次是不是电脑重复计算了六次啊，那就亏大了我是刚学习也不知哪里出问题了

[ Last edited by ben_ladeng on 2013-5-15 at 18:21 ]

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

VASP5.3 k-points parallelism计算, KPAR的设置，KPAR=2 已经有11人回复
VASP中如何创建 INCAR , POSCAR, KPOINTS和POTCAR文件已经有7人回复
测试POSCAR的脚本不运行已经有10人回复
关于VASP一些琐事的求助（K点，吸附优化等）已经有12人回复
ENCUT测试脚本已经有17人回复
KPOINTS网格怎么写已经有3人回复
VASP计算时的几个细节设置问题，希望有经验的人过来指点一下！已经有23人回复
【原创】CPMD3.15.1 和3.13.2在Ubuntu下的自动下载，安装，测试脚本（Gfortran）已经有12人回复
【求助成功】vasp测试K点收敛性的脚本问题已经有20人回复
【求助】VASP计算能带时KPOINTS设置出错已经有14人回复
【求助】Hse06计算时KPOINTS出错已经有30人回复
【求助成功】晶格优化的结果很诡异已经有13人回复
【求助完毕】IBZKPT,KPOINTS与DOS计算已经有4人回复

1楼 2013-05-12 23:52:08

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

GTF1902

金虫 (小有名气)

应助: 2 (幼儿园)
金币: 1045.4
散金: 20
帖子: 188
在线: 63.6小时
虫号: 1720390
注册: 2012-03-27
专业: 凝聚态物性I:结构、力学和

GTF1902: 回帖置顶 2013-05-13 10:01:55

引用回帖:

2楼: Originally posted by 171713294 at 2013-05-13 08:45:57
楼住，给你个建议，你吧你的问题一列列弄好，别这样一大片，没人会给你看的，重新弄好我们给你看看

我用vasp优化结构前先进行K测试，下面是我的incar 和K_RUN文件

System = PCCOF
ISIF = 2
ISMEAR = 0
ISTART = 0
ICHARG = 2
SIGMA = 0.1
NSW = 0
PREC = Medium
ENCUT = 500
GGA = 91
VOSKOWN = 1
IBRION = 2
EDIFF = 1E-4
EDIFFG = -0.03
LREAL = Auto
IALGO = 48
LALANE = .TRUE
NPAR = 6
LSCALU = .FALSE.
NSIM = 4

我是用先用下面的txt提交任务：
APP_NAME=scivip
OMP_NUM_THREADS=1
NP_PER_NODE=6
MY_MPI_TYPE=openmpi
MY_MPI_HOME=/share/home/sci_028/sourcecode/openmpi-1.4.2
NP=6
RUN="/share/home/sci_028/gtf/vasp/shiyansan/run_k.txt >out.log"

下面是的 run_k.txt：
#!/bin/sh
rm -f WAVECAR
for i in 1 2 3 4 5 6 7
do
cat > KPOINTS <<!
Automatic generation
0
Monkhorst-Pack
$i $i $i
0.0 0.0 0.0
!
echo "k mesh = $i x $i x $i"; time /share/home/sci_028/sourcecode/vasp.5.2/vasp >out.log
E=`grep "TOTEN" OUTCAR |tail -1| awk '{printf "%12.6f\n", $5 }'`
KP=`grep "irreducible" OUTCAR |tail -1| awk '{printf "%5i\n", $2 }'`
echo $i $KP $E >>comment
done

但是comment输出却为：

1 1 -569.842023
1 1 -569.842023
1 1 -569.842023
1 1 -569.842023
1 1 -569.842023
1 1 -569.842023
2 1 -600.814240
2 1 -600.814240
2 1 -600.814240
2 1 -600.814240
2 1 -600.814240
2 1 -600.814240
3 8 -599.488993
3 8 -599.488993
3 8 -599.488993
3 8 -599.488993
3 8 -599.488993
3 8 -599.488993
4 8 -599.575041
4 8 -599.575041
4 8 -599.575041
4 8 -599.575041
4 8 -599.575041
4 8 -599.575041
5 27 -599.519665
5 27 -599.519665
5 27 -599.519665
5 27 -599.519665
5 27 -599.519665
5 27 -599.519665
6 27 -599.540837
6 27 -599.540837
6 27 -599.540837
6 27 -599.540837
6 27 -599.540837
6 27 -599.540837
7 64 -599.557505
7 64 -599.557505
7 64 -599.557505
7 64 -599.557505
7 64 -599.557505
7 64 -599.557505
也就说没个k点输出了六次是不是计算了六次啊真的好郁闷那就亏大了问题出在哪里？