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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP中脚本程序问题
北京石油化工学院2026年研究生招生接收调剂公告
151/2返回列表12下一页
查看: 2617  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 复试调剂 +13 呼呼？~+123456 2026-04-05 13/650 2026-04-05 22:07 by 醉翁wl
[考研] 301求调剂 +5 细胞相关蛋白 2026-04-02 9/450 2026-04-05 20:14 by lys0704
[考研] 一志愿9材料学硕297已过六级求调剂推荐 +11 adaie 2026-04-04 12/600 2026-04-05 19:04 by 蓝云思雨
[考研] 312求调剂 +3 Say Never 2026-04-04 3/150 2026-04-05 11:19 by guoweigw
[考研] 化学357分，考研调剂 +10 .Starry. 2026-04-04 11/550 2026-04-05 10:57 by cql1109
[考研] 一志愿电子科技大学085600材料与化工 329分求调剂 +10 Naiko 2026-04-04 10/500 2026-04-05 09:40 by sam3303
[考研] 调剂 +8 熊二想上岸 2026-04-04 8/400 2026-04-05 05:27 by houyaoxu
[考研] 材料专硕322分 +11 哈哈哈吼吼吼哈 2026-04-02 11/550 2026-04-04 23:37 by 永字号
[考博] 申博 +7 IQwQl 2026-04-04 7/350 2026-04-04 23:32 by mumin1990
[考研] 368求调剂 +5 今华习 2026-04-03 7/350 2026-04-04 18:47 by imissbao
[考研] 292求调剂 +11 2022080213 2026-04-04 13/650 2026-04-04 18:38 by macy2011
[考研] 085701求调剂 +7 龚禹铭 2026-04-04 8/400 2026-04-04 13:49 by 小小树2024
[考研] 266求调剂 +8 学员97LZgn 2026-04-03 8/400 2026-04-04 09:02 by 20021109
[考研] 化学调剂求助 +6 LULONG1 2026-04-03 6/300 2026-04-03 23:13 by qzxyhcsy
[考研] 0703化学 +7 goldtt 2026-04-02 9/450 2026-04-03 09:36 by 蓝云思雨
[考研] 326求调剂 +10 崽崽仔 2026-04-02 10/500 2026-04-03 09:08 by 帕尔马拉特
[考研] 考研调剂 +3 李木子0120 2026-04-02 5/250 2026-04-02 21:45 by dongzh2009
[考研] 求调剂 302分初试 0854 +5 伶可乐 2026-04-02 5/250 2026-04-02 17:53 by 笔落锦州
[考研] 一志愿厦门大学化学工程（专硕）-数二英二406分-求调剂 +5 厦大化工 2026-04-01 5/250 2026-04-02 10:03 by jp9609
[考研] 求调剂：一志愿：南京大学 专业：0705 总分320 ，本科985，四六级已过 +3 lfy760306 2026-03-31 3/150 2026-04-01 01:57 by Creta
信息提示
关闭
请填处理意见
关闭