24小时热门版块排行榜    

查看: 1042  |  回复: 2

GTF1902

金虫 (小有名气)

[求助] 【重新开帖】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 ]
回复此楼

» 猜你喜欢

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

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

171713294

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★
感谢参与,应助指数 +1
GTF1902: 金币+2, 有帮助 2013-05-13 09:16:41
GTF1902: 回帖置顶 2013-05-13 10:01:05
franch: 金币+2, 谢谢回帖交流,, 2013-05-13 21:19:50
楼住,给你个建议,你吧你的问题一列列弄好,别这样一大片,没人会给你看的,重新弄好 我们给你看看
坚持,不放弃
2楼2013-05-13 08:45:57
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

GTF1902

金虫 (小有名气)

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点输出了六次  是不是计算了六次啊  真的好郁闷  那就亏大了  问题出在哪里?
3楼2013-05-13 10:01:00
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 GTF1902 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 323求调剂 +6 洼小桶 2026-03-18 6/300 2026-03-23 00:29 by king123!
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +11 枫桥ZL 2026-03-18 13/650 2026-03-22 20:26 by edmund7
[考研] 北科281学硕材料求调剂 +6 tcxiaoxx 2026-03-20 6/300 2026-03-22 20:23 by edmund7
[考研] 298求调剂 +6 上岸6666@ 2026-03-20 6/300 2026-03-22 20:21 by edmund7
[考研] 284求调剂 +5 Zhao anqi 2026-03-22 5/250 2026-03-22 17:38 by barlinike
[考研] 求调剂 +6 十三加油 2026-03-21 6/300 2026-03-22 17:00 by i_cooler
[考研] 寻找调剂 +4 倔强芒? 2026-03-21 4/200 2026-03-22 16:14 by 木托莫露露
[考研] 一志愿中南化学(0703)总分337求调剂 +9 niko- 2026-03-19 10/500 2026-03-22 16:08 by ColorlessPI
[考研] 求调剂 +4 要好好无聊 2026-03-21 4/200 2026-03-21 18:57 by 学员8dgXkO
[考研] 313求调剂 +4 肆叁贰壹22 2026-03-19 4/200 2026-03-21 17:33 by ColorlessPI
[考研] 303求调剂 +5 睿08 2026-03-17 7/350 2026-03-21 03:11 by JourneyLucky
[考研] 332求调剂 +4 ydfyh 2026-03-17 4/200 2026-03-21 02:20 by JourneyLucky
[考研] 材料 336 求调剂 +3 An@. 2026-03-18 4/200 2026-03-21 01:39 by JourneyLucky
[考研] 271材料工程求调剂 +8 .6lL 2026-03-18 8/400 2026-03-21 00:58 by JourneyLucky
[考研] 350求调剂 +5 weudhdk 2026-03-19 5/250 2026-03-20 22:04 by luoyongfeng
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 295复试调剂 +8 简木ChuFront 2026-03-19 8/400 2026-03-20 20:44 by zhukairuo
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
[考研] 298-一志愿中国农业大学-求调剂 +9 手机用户 2026-03-17 9/450 2026-03-20 14:24 by 无懈可击111
信息提示
请填处理意见