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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP中脚本程序问题
151/2返回列表12下一页
查看: 2650  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

aminibuck

捐助贵宾 (正式写手)

[求助] VASP中脚本程序问题

我初学VASP，看侯老师的教程中有时候需要编写脚本程序，比如优化晶格常数有如下脚本：#!/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
1
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 '{prinf "%12.4f\n",$5}''
echo $V $E>>EtVo.dat
done

我想问的是这个是写在POSCAR文件里吗？如果不是的话就是另写一个脚本吧？我们学校的高性能计算中心利用PBS提交作业，请问如果另写一个脚本的话我的PBS文件中应该加入什么命令呢？谢谢了
回复此楼

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

第一原理资料贴

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-12-03 21:53:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yh3166

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
aminibuck: 金币+50, ★有帮助 2013-12-04 21:15:38
这种小程序是随便建个文档写入，写完保存后变成执行文件就好
最后，关于pbs上面如何提交作业的话我是这样写的，你可以参考下
#!/bin/sh

#PBS -N vasp

#PBS -d ./

#PBS -o ./

#PBS -e ./

#PBS -l nodes=1:ppn=4

#PBS -l walltime=10:00:00

#PBS -l mem=

#PBS -q serique

#PBS -V

ulimit -s unlimited
source /opt/intel/Compiler/11.1/046/bin/intel64/iccvars_intel64.sh
source /opt/intel/Compiler/11.1/046/bin/intel64/ifortvars_intel64.sh
source /opt/intel/mkl/10.2.2.025/tools/environment/mklvarsem64t.sh
source /opt/intel/impi/3.2.2/bin64/mpivars.sh
EXEC=/***/***/vasp_5.2
################################################################
rm WAVECAR
for i in 250 300 350 400 450 500 550 600
do
cat > INCAR <<!
System = Si
ENCUT = $i
ISTART = 0
ICHARG = 2
ISMEAR = 0
SIGMA = 0.1
PREC = Accurate
!
echo "ENCUT = $i eV" ;  /***/***/ vasp_5.2(这里是写vasp的路径，你自己是什么写什么)
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
赞一下 回复此楼
2楼2013-12-04 20:55:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aminibuck

捐助贵宾 (正式写手)
引用回帖:
2楼: Originally posted by yh3166 at 2013-12-04 20:55:11
这种小程序是随便建个文档写入，写完保存后变成执行文件就好
最后，关于pbs上面如何提交作业的话我是这样写的，你可以参考下
#!/bin/sh

#PBS -N vasp

#PBS -d ./

#PBS -o ./

#PBS -e ./

#PBS -l n ...

你好，非常感谢你的回答。照你的意思是直接把脚本程序写进PBS脚本里面是吧？那INCAR和KPOINTS这样的文件也还是必须的吧？比如说按照你这个PBS模板是不是依旧还需要写一个INCAR文件如下：
SYSTEM=Si
ENCUT=$i
ISTART=0
ICHARGE=2
ISMEAR=0;SIGMA=0.1
PREC=Accurate
请问是这样理解吗？
赞一下 回复此楼
3楼2013-12-04 21:15:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yh3166

金虫 (正式写手)
嗯，脚本写进pbs里，其他INCAR,KPOINTS,POTCAR都还是要有的。
赞一下 回复此楼
4楼2013-12-05 08:49:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

djh434325996

木虫 (小有名气)
引用回帖:
3楼: Originally posted by aminibuck at 2013-12-04 21:15:15
你好，非常感谢你的回答。照你的意思是直接把脚本程序写进PBS脚本里面是吧？那INCAR和KPOINTS这样的文件也还是必须的吧？比如说按照你这个PBS模板是不是依旧还需要写一个INCAR文件如下：
SYSTEM=Si
ENCUT=$i
IS ...

侯老师的cat > INCAR <<!
               。。。
               。。。
               ！             这个部分就是在写INCAR。
其他几个CAR，你要事先放在运行的目录下。
赞一下 回复此楼
5楼2013-12-05 08:57:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aminibuck

捐助贵宾 (正式写手)
引用回帖:
4楼: Originally posted by yh3166 at 2013-12-05 08:49:30
嗯，脚本写进pbs里，其他INCAR,KPOINTS,POTCAR都还是要有的。

好的，谢谢。我自己试一试。
赞一下 回复此楼
6楼2013-12-05 09:32:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aminibuck

捐助贵宾 (正式写手)
引用回帖:
5楼: Originally posted by djh434325996 at 2013-12-05 08:57:46
侯老师的cat > INCAR <<!
               。。。
               。。。
               ！             这个部分就是在写INCAR。
其他几个CAR，你要事先放在运行的目录下。...

谢谢，那PBS脚本里面是否还需要再次写一遍？
赞一下 回复此楼
7楼2013-12-05 09:33:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yh3166

金虫 (正式写手)
引用回帖:
7楼: Originally posted by aminibuck at 2013-12-05 09:33:44
谢谢，那PBS脚本里面是否还需要再次写一遍？...

不用的，你这个脚本有POSCAR吗，再说，你这个只是测试，还没到真正计算
赞一下 回复此楼
8楼2013-12-05 15:07:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aminibuck

捐助贵宾 (正式写手)
引用回帖:
8楼: Originally posted by yh3166 at 2013-12-05 15:07:53
不用的，你这个脚本有POSCAR吗，再说，你这个只是测试，还没到真正计算...

嗯，我明白了，就拿你贴出来的PBS来说它会自动生成INCAR文件，无需再自己写了。谢谢
赞一下 回复此楼
9楼2013-12-05 18:25:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

djh434325996

木虫 (小有名气)
引用回帖:
7楼: Originally posted by aminibuck at 2013-12-05 09:33:44
谢谢，那PBS脚本里面是否还需要再次写一遍？...

其他几个CAR不用写在PBS脚本里面
for i in 250 300 350 400 450 500 550 600
do
cat > INCAR <<!
System = Si
ENCUT = $i
ISTART = 0
ICHARG = 2
ISMEAR = 0
SIGMA = 0.1
PREC = Accurate
!
只是一个循环，测试不同的ECUT
赞一下 回复此楼
10楼2013-12-06 11:30:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 aminibuck 的主题更新
151/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然面上和省基金B类撒花 +18 花田半亩～白 2026-04-21 18/900 2026-04-23 11:31 by 12021227
[考研] 有没有学校收留 +3 蒋昌鹏qtj 2026-04-20 3/150 2026-04-22 20:25 by 学员JpLReM
[考研] 312求调剂 +3 山河似你温柔 2026-04-22 3/150 2026-04-22 20:17 by 学员JpLReM
[考博] 华师大读博 +3 xq83 2026-04-22 5/250 2026-04-22 10:42 by xq83
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +5 yexuqing 2026-04-19 5/250 2026-04-22 10:01 by easeheart
[论文投稿] 急需审稿人！！！ +3 陆小果画大饼 2026-04-21 3/150 2026-04-21 23:54 by jzy_123456
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 8/400 2026-04-21 16:26 by 啃面包的小书虫
[考研] 295分求调剂 +6 ?要上岸? 2026-04-17 6/300 2026-04-21 08:18 by Equinoxhua
[考研] 085600材料与化工调剂 5+3 孜孜不倦2002 2026-04-19 6/300 2026-04-20 21:25 by babero
[论文投稿] 有没有接收比较快的sci期刊呀，最好在一个月之内的，研三孩子求毕业 20+4 之护着 2026-04-16 7/350 2026-04-20 15:45 by 豆豆7758
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 湖南大学刘巧玲课题组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叶
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 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 中豫男
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[有机交流] 二苯甲酮酸类衍生物 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
信息提示
关闭
请填处理意见
关闭