24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1054  |  回复: 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的回帖
回帖置顶 ( 共有2个 )

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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿中南大学 0855 机械 286 求调剂 +10 不会吃肉 2026-04-12 10/500 2026-04-12 22:51 by 零零二
[考研] 材料相关专业344求调剂双非工科学校或课题组 +9 hualkop 2026-04-12 9/450 2026-04-12 22:41 by kyle12138
[考研] 272分材料子求调剂 +40 Loy0361 2026-04-10 52/2600 2026-04-12 19:23 by laoshidan
[考研] 一志愿哈工大 085600 277 12材科基求调剂 5+5 chenny174 2026-04-10 31/1550 2026-04-12 19:23 by qingfeng258
[考研] 297求调剂 +17 ORCHID1 2026-04-10 18/900 2026-04-12 09:59 by zhouwenxian
[考研] 298求调剂 +5 残荷新柳 2026-04-07 5/250 2026-04-11 11:02 by 紫曦紫棋
[考研] 农学0904 312求调剂 +6 Say Never 2026-04-10 6/300 2026-04-11 10:33 by wwj2530616
[考研] 087100初试311求调剂 +4 任雅琴 2026-04-09 4/200 2026-04-11 10:33 by zhq0425
[考研] 本科211 工科085400 280分求调剂 可跨专业 +11 LZH(等待调剂中 2026-04-10 11/550 2026-04-11 08:39 by zhq0425
[考研] 调剂 +19 2261744733 2026-04-08 19/950 2026-04-09 19:11 by vgtyfty
[考研] 材料工程调剂 +12 小刘同学吖吖 2026-04-06 13/650 2026-04-09 17:07 by luoyongfeng
[考研] 368化学求调剂 +13 wwwwabcde 2026-04-07 14/700 2026-04-09 14:47 by heaven_jay
[考研] 一志愿西南大学生物学学硕344 求生物学相关调剂/生物与医药 +7 超人不会飞@ 2026-04-08 7/350 2026-04-09 09:35 by gong120082
[考研] 机械专硕273请求调剂 +6 庚申壬申 2026-04-07 6/300 2026-04-08 22:41 by bljnqdcc
[考研] 考研求调剂 +4 雯??? 2026-04-08 4/200 2026-04-08 21:44 by 土木硕士招生
[考研] 0703化学调剂 348分 +14 唉我超真没招了 2026-04-06 15/750 2026-04-08 19:16 by 我减肥1
[考研] 求调剂,现在还能填的 +3 上岸小莹加油 2026-04-08 3/150 2026-04-08 14:30 by zhq0425
[考研] 338求调剂 +8 wxygxsaaaaa 2026-04-06 8/400 2026-04-08 06:58 by 无际的草原
[考研] 专硕085403,291分,有两篇专利,一国一奖 +3 哈吉咪哈吉咪 2026-04-07 3/150 2026-04-07 18:21 by 蓝云思雨
[考研] 328求调剂 +4 ghhh88888 2026-04-06 5/250 2026-04-07 14:45 by ghhh88888
信息提示
请填处理意见