24小时热门版块排行榜    

查看: 1539  |  回复: 9

lisu555

铁杆木虫 (著名写手)

[求助] 侯老师 K点测试问题,求帮解决!

问题如下:
[root@localhost Run_k-Fe-PAWPBE]# bash run_k.sh
rm: 无法删除"WAVECAR": 没有那个文件或目录
k mesh = 5 x 5 x 5
Error reading item 'IMAGES' from file INCAR.
grep: OUTCAR: 没有那个文件或目录
grep: OUTCAR: 没有那个文件或目录
k mesh = 7 x 7 x 7
Error reading item 'IMAGES' from file INCAR.
grep: OUTCAR: 没有那个文件或目录
grep: OUTCAR: 没有那个文件或目录
k mesh = 9 x 9 x 9
Error reading item 'IMAGES' from file INCAR.
grep: OUTCAR: 没有那个文件或目录
grep: OUTCAR: 没有那个文件或目录
k mesh = 11 x 11 x 11
Error reading item 'IMAGES' from file INCAR.
grep: OUTCAR: 没有那个文件或目录
grep: OUTCAR: 没有那个文件或目录
k mesh = 13 x 13 x 13
Error reading item 'IMAGES' from file INCAR.
grep: OUTCAR: 没有那个文件或目录
grep: OUTCAR: 没有那个文件或目录
k mesh = 15 x 15 x 15
Error reading item 'IMAGES' from file INCAR.




run_k如下:
#!/bin/sh
rm WAVECAR
for i in 5  7  9  11  13  15
do
cat > KPOINTS <<!
Automatic generation
0
Monhkorst-Pack
$i $i $i
0.0 0.0 0.0
!
echo "k mesh = $i x $i x $i" ;  vasp5.2
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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lisu555

铁杆木虫 (著名写手)

其它文件INCAR POTCAR 和poscar都是和教程上一样的!
2楼2016-08-29 22:02:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

从心而觅

新虫 (小有名气)

引用回帖:
2楼: Originally posted by lisu555 at 2016-08-29 22:02:34
其它文件INCAR POTCAR 和poscar都是和教程上一样的!

请问你是怎么提交脚本的?
一切有为法,如梦幻泡影!如露亦如电,应作如是观!
3楼2016-08-30 08:01:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lisu555

铁杆木虫 (著名写手)

引用回帖:
3楼: Originally posted by 从心而觅 at 2016-08-30 08:01:28
请问你是怎么提交脚本的?...

我用的是CentOS系统:dos2unix run_k→bash run_k即可。
4楼2016-08-30 09:32:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

afan0012

木虫 (小有名气)

你贴出来的脚本里面,并没有进行VASP的计算呀

发自小木虫Android客户端
5楼2016-08-31 23:52:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

afan0012

木虫 (小有名气)

哦哦,看到了,你用vasp5.2来运行VASP。运行脚本后,产生输出文件了吗?

发自小木虫Android客户端
6楼2016-08-31 23:56:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

afan0012

木虫 (小有名气)

我猜测可能是你执行vasp5.2时,没有指定哪个节点,多少核,所以任务并没有执行(如果是并行);如果你是串行的,建议不用脚本,写好KPOINTS文件执行一下vasp5.2看能不能正常运行

发自小木虫Android客户端
7楼2016-08-31 23:59:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lisu555

铁杆木虫 (著名写手)

引用回帖:
7楼: Originally posted by afan0012 at 2016-08-31 23:59:02
我猜测可能是你执行vasp5.2时,没有指定哪个节点,多少核,所以任务并没有执行(如果是并行);如果你是串行的,建议不用脚本,写好KPOINTS文件执行一下vasp5.2看能不能正常运行
...

单机串行 ,用脚本方便啊,输出文件没有能量!
8楼2016-09-01 14:01:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

afan0012

木虫 (小有名气)

引用回帖:
8楼: Originally posted by lisu555 at 2016-09-01 14:01:32
单机串行 ,用脚本方便啊,输出文件没有能量!...

我的意思是,不要用脚本改变KPOINTS,自己写好K点文件,执行vasp5.2,看是否有输出文件

发自小木虫Android客户端
9楼2016-09-01 22:47:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

afan0012

木虫 (小有名气)

脚本执行到grep ** OUTCAR那里出错的,很可能没产生OUTCAR

发自小木虫Android客户端
10楼2016-09-01 22:49:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lisu555 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 EZ67DG6CvcKL 2026-08-31 4/200 2026-09-01 06:06 by o3EMxHUAQDtS
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 EZ67DG6CvcKL 2026-08-31 5/250 2026-09-01 05:43 by o3EMxHUAQDtS
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 wsYGJnURYvm6 2026-08-31 3/150 2026-09-01 01:52 by o3EMxHUAQDtS
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见