应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】VASP脚本程序怎么运行？
181/2返回列表12下一页
查看: 3179  |  回复: 17
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

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 ]
回复此楼

» 收录本帖的淘帖专辑推荐

Photochemistry 计算-vasp

» 猜你喜欢
1楼 2009-05-09 10:46:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

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循环。

其实这些工具脚本都很简单的，关键还是要明白输出文件，然后就可以自己写一些适合自己用的。
赞一下(13人) 回复此楼
4楼2009-05-09 20:20:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

405963512

铜虫 (正式写手)
★
fegg7502(金币+1,VIP+0):鼓励交流！ 5-10 18:41
我觉得他问的没这么难吧，那个脚本程序是做晶格常数测试的，看取的这些a值那个能量最低，那么这个a所对应的值作为标准的晶格常数，带入 到poscar，进行后续的计算。这是他的作用
至于运行方法，chmod +x 文件名；这个命令吧你的shell变成可执行程序，然后敲入这个文件名，就可运行
我很笨，是弄了一周才弄明白的，虽不是什么大问题，和大家分享！
赞一下(11人) 回复此楼
6楼2009-05-09 21:13:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

含草

金虫 (小有名气)
★ ★
fegg7502(金币+2,VIP+0):鼓励交流！ 5-10 18:41
运行脚本文件之后,文件里面根本就没有数据,这是什么原因?谢谢
./opt.sh: line 14: echoa=3.70 angstrom: command not found
running on    1 nodes
distr:  one band on    1 nodes,    1 groups
vasp.4.6.31 08Feb07 complex
POSCAR found :  0 types and    0 ions
p0_31941:  p4_error: interrupt SIGSEGV: 11

real    0m0.022s
user    0m0.019s
sys     0m0.003s
./opt.sh: command substitution: line 16: unexpected EOF while looking for matching `"'
./opt.sh: command substitution: line 17: syntax error: unexpected end of file
./opt.sh: line 14: echoa=3.80 angstrom: command not found
running on    1 nodes
distr:  one band on    1 nodes,    1 groups
vasp.4.6.31 08Feb07 complex
POSCAR found :  0 types and    0 ions
p0_31949:  p4_error: interrupt SIGSEGV: 11

real    0m0.022s
user    0m0.019s
sys     0m0.003s
./opt.sh: command substitution: line 16: unexpected EOF while looking for matching `"'
./opt.sh: command substitution: line 17: syntax error: unexpected end of file
./opt.sh: line 14: echoa=3.90 angstrom: command not found
running on    1 nodes
distr:  one band on    1 nodes,    1 groups
vasp.4.6.31 08Feb07 complex
POSCAR found :  0 types and    0 ions
赞一下(11人) 回复此楼
8楼2009-05-09 21:57:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sars518

木虫 (著名写手)

小木虫灌水大队队长
★
zxzj05(金币+1,VIP+0):3Q! ^_^ 期待你的更多精彩回答！ 5-12 10:11
!
echo "a =$i angstrom";time vasp


改成

!
vasp
echo "a =$i angstrom"

[ Last edited by sars518 on 2009-5-11 at 21:48 ]
赞一下(11人) 回复此楼
人生是一场旅行，在欣赏风景的同时，别忘了你也是别人的风景。
13楼2009-05-11 21:43:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

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
赞一下(10人) 回复此楼
2楼2009-05-09 10:54:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sars518

木虫 (著名写手)

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

含草

金虫 (小有名气)
学习了很多
回复此楼
5楼2009-05-09 20:30:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tjpm

金虫 (正式写手)
★ ★
fegg7502(金币+1,VIP+0):鼓励交流！ 5-10 18:41
zxzj05(金币+1,VIP+0):3Q! ^_^ 期待你的更多精彩回答！ 5-12 10:10
恩
./run_a0就可以了
引用回帖:
Originally posted by 405963512 at 2009-5-9 21:13:
我觉得他问的没这么难吧，那个脚本程序是做晶格常数测试的，看取的这些a值那个能量最低，那么这个a所对应的值作为标准的晶格常数，带入 到poscar，进行后续的计算。这是他的作用
至于运行方法，chmod +x 文件名； ...

赞一下(20人) 回复此楼
7楼2009-05-09 21:35:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tjpm

金虫 (正式写手)
★ ★ ★
fegg7502(金币+3,VIP+0):鼓励交流！ 5-10 18:41
你的脚本有问题。
1、./opt.sh: line 14: echoa=3.70 angstrom: command not found
是echo "a=$a angstrom",你少了空格。
2、/opt.sh: command substitution: line 16: unexpected EOF while looking for matching `"'
./opt.sh: command substitution: line 17: syntax error: unexpected end of file
用cat创建文件时，重定向符“<<”后 前后截断字符要相同。
赞一下(10人) 回复此楼
9楼2009-05-10 00:14:04
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 longchang 的主题更新
181/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 291求调剂 +12 关忆北. 2026-04-14 12/600 2026-04-21 15:35 by 18300969621
[考研] 一志愿中科大材料与化工，353分还有调剂学校吗 +11 否极泰来2026 2026-04-15 13/650 2026-04-20 22:31 by Equinoxhua
[论文投稿] 期刊推荐 +3 材料研究生 2026-04-15 5/250 2026-04-20 16:02 by 豆豆7758
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +4 yexuqing 2026-04-19 4/200 2026-04-20 14:47 by brantleo
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 求调剂推荐 +9 小聂爱学习 2026-04-14 9/450 2026-04-19 17:03 by 中豫男
[考研] 求调剂 +10 小聂爱学习 2026-04-16 12/600 2026-04-19 16:51 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[考研] 0854求调剂 +23 门路摸摸 2026-04-15 27/1350 2026-04-19 01:59 by 烟雨流涯
[考研] 300求调剂 +12 橙a777 2026-04-15 12/600 2026-04-18 23:51 by 路病情
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 297，工科调剂? +5 河南农业大学-能 2026-04-14 5/250 2026-04-18 15:17 by Equinoxhua
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[考研] 一志愿华中农业071010，320求调剂 +17 困困困困坤坤 2026-04-14 19/950 2026-04-17 20:08 by 关一盏灯cd
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
[考研] 一志愿沪9，生物学326求调剂 +9 刘墨墨 2026-04-15 9/450 2026-04-16 17:14 by 崔崔崔cccc
信息提示
关闭
请填处理意见
关闭