24小时热门版块排行榜    

查看: 3112  |  回复: 17
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

longchang

金虫 (小有名气)

[交流] 【求助】VASP脚本程序怎么运行?

刚学VASP,对其还不是很熟悉;现在正在试着运行脚本文件(按侯老师的例子),准备好INCAR、KPOINTS、POTCAR文件后,应该准备POSCAR文件,那脚本程序run_a0怎样运行?脚本run_a0为:
#!/bin/sh
rm WAVECAR
for i in 3.80 3.85 3.90 3.95 4.00 4.05 4.10 4.15 4.20 4.25 4.30
do
cat >POSCAR < Al-fcc
$i
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
Direct
0.0 0.0 0.0
!
echo "a =$i angstrom";time vasp
E=‘grep "TOTEN" OUTCAR |tail -1|awk ‘{printf "%12.6f \n",$5}'‘
V=‘grep "volume" OUTCAR |tail -1|awk ‘{printf "%12.4f \n",$5}'‘
echo $V  $E >>EtVo.dat
done
VASP运行需四个文件,这里难道用run_a0替代POSCAR文件?但那样的话VASP运行不了,但又不能单独运行脚本程序(出不了结果),到底该怎麽办?哪位能告诉我具体的运行步骤?谢谢!

[ Last edited by aylayl08 on 2010-4-21 at 09:01 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

longchang

金虫 (小有名气)

rm: cannot remove `WAVECAR': No such file or directory
a =3.80 angstrom
bash: vasp: command not found

real    0m0.001s
user    0m0.000s
sys     0m0.000s
bash: TOTEN: command not found
awk: cmd. line:1:
awk: cmd. line:1: ^ invalid char '�' in expression
a =3.85 angstrom
bash: vasp: command not found

real    0m0.001s
user    0m0.000s
sys     0m0.000s
bash: TOTEN: command not found
awk: cmd. line:1:
awk: cmd. line:1: ^ invalid char '�' in expression
a =3.90 angstrom
bash: vasp: command not found

real    0m0.001s
user    0m0.000s
sys     0m0.000s
bash: TOTEN: command not found
awk: cmd. line:1:
awk: cmd. line:1: ^ invalid char '�' in expression
a =3.95 angstrom
bash: vasp: command not found
real    0m0.001s
user    0m0.000s
sys     0m0.000s
bash: TOTEN: command not found
awk: cmd. line:1:
awk: cmd. line:1: ^ invalid char '�' in expression
a =4.00 angstrom
bash: vasp: command not found

real    0m0.001s
user    0m0.001s
sys     0m0.000s
bash: TOTEN: command not found
.............
...............
...............
是怎么回事?难道VASP运行脚本需要”特殊“代码?
10楼2009-05-10 21:14:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 18 个回答

emily6526

铁虫 (初入文坛)

★ ★ ★ ★
fegg7502(金币+4,VIP+0):鼓励交流! 5-9 17:10
可以自己写一个POSCAR,把$i换成一个具体的数。
Al-fcc
$i
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
Direct
0.0 0.0 0.0
2楼2009-05-09 10:54:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sars518

木虫 (著名写手)

小木虫灌水大队队长

★ ★
fegg7502(金币+2,VIP+0):鼓励交流! 5-9 17:10
其实这个脚本的作用就是 改变 缩放系数  写POSCAR
人生是一场旅行,在欣赏风景的同时,别忘了你也是别人的风景。
3楼2009-05-09 16:06:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tjpm

金虫 (正式写手)

★ ★
fegg7502(金币+2,VIP+0):thank you very much! 5-10 18:41
其实shell脚本很简单的,有C语言基础随便看看就会了,当然还需要熟悉常用的linux命令。
简单写个注释吧
====
#!/bin/sh #每种shell脚本开头需要指定起类型,比如unix下原始的就是sh,linux都标配bash了, Solaris也是用bash了。功能上现在似乎zsh比较好,目前我是把zsh做默认shell.
rm WAVECAR  
删除WAVECAR,其实这里有点不好,因为很多系统默认把rm 替换为rm -i
也就是说你要删除这个文件时会跳出一个提示问是否要删除。建议改为 rm  WAVECAR -rf
for i in 3.80 3.85 3.90 3.95 4.00 4.05 4.10 4.15 4.20 4.25 4.30
do
这个sh/bash下典型的for循环,就是让i逐次取后面的值
cat >POSCAR < Al-fcc
$i
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
Direct
0.0 0.0 0.0
!
这个就是利用管道创建文件POSCAR,两个“!”分别标识开头和结尾。
echo "a =$i angstrom";time vasp
屏幕输出 a的值,并运行vasp, time命令加在前面就是计算其运行时间。
E=‘grep "TOTEN" OUTCAR |tail -1|awk ‘{printf "%12.6f \n",$5}'‘
搜索包含"TOTEN"的所有行,利用tail -1取最后一行,然后格式化输出第5个字符段。
V=‘grep "volume" OUTCAR |tail -1|awk ‘{printf "%12.4f \n",$5}'‘
这个同上
echo $V  $E >>EtVo.dat
把V,E的值添加到EtVo.dat
done
表示结束for循环。

其实这些工具脚本都很简单的,关键还是要明白输出文件,然后就可以自己写一些适合自己用的。
4楼2009-05-09 20:20:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学调剂0703 +7 啊我我的 2026-03-11 7/350 2026-03-15 23:03 by 凌千颂111
[考研] 机械专硕调剂 +3 笨笨兔子 2026-03-12 3/150 2026-03-15 20:02 by 栗子粥?
[考博] 欢迎申博同学联系 +3 天道酬勤2026686 2026-03-10 7/350 2026-03-15 19:03 by 天道酬勤2026686
[考研] 311求调剂 +3 26研0 2026-03-15 3/150 2026-03-15 09:12 by JourneyLucky
[考研] 265求调剂 +4 威化饼07 2026-03-12 4/200 2026-03-14 17:23 by userper
[考研] 一志愿天大化工(085600)调剂总分338 +6 蔡大美女 2026-03-09 6/300 2026-03-14 02:46 by JourneyLucky
[考研] 332分材料工程调剂 +3 莓好时光海苔 2026-03-09 3/150 2026-03-14 02:03 by JourneyLucky
[考研] 一志愿天津大学,英一数二305分求调剂,四六级已过 +8 小小番的茄 2026-03-09 8/400 2026-03-14 01:53 by JourneyLucky
[考研] 265求调剂 +9 小木虫085600 2026-03-09 12/600 2026-03-14 01:11 by JourneyLucky
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
[考研] 285求调剂 +4 ytter 2026-03-12 4/200 2026-03-13 14:48 by jxchenghu
[考研] 085600材料与化工 309分请求调剂 +7 dtdxzxx 2026-03-12 8/400 2026-03-13 14:43 by jxchenghu
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-12 5/250 2026-03-13 10:56 by houyaoxu
[考研] 纺织、生物、化学、材料相关专业招生了 +4 耶耶业 2026-03-09 7/350 2026-03-12 19:05 by Equinoxhua
[考研] 0857环境调剂 +5 熠熠_11 2026-03-10 5/250 2026-03-11 10:59 by wang_dand
[考研] 求调剂材料专硕293 +6 段_(:з」∠)_ 2026-03-10 6/300 2026-03-10 18:22 by ms629
[硕博家园] 木虫好像不热闹了,是不是? +4 偏振片 2026-03-10 4/200 2026-03-10 09:51 by longwave
信息提示
请填处理意见