| 查看: 2938 | 回复: 14 | |||
[求助]
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文件中应该加入什么命令呢？谢谢了 |
回复此楼» 收录本帖的淘帖专辑推荐
 第一原理资料贴 |
» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
今天基金会出结果吗？20260819
已经有13人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
时间戳变了，能看出什么问题？
已经有16人回复
-
filecode，4个jtjc了
已经有13人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有121人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有偿求vasp指导
已经有17人回复
- VASP中如何查看正在运行程序的源文件
已经有3人回复
- VASP 运行报错
已经有11人回复
- 并行vasp运行问题
已经有10人回复
- 新手vasp安装求助
已经有11人回复
- k均值聚类程序求助
已经有7人回复
- 测试POSCAR的脚本不运行
已经有10人回复
- vasp脚本运行问题
已经有7人回复
- ENCUT 的脚本程序怎么用
已经有15人回复
- vasp运行被强制终止是什么原因引起的？
已经有7人回复
- 【求助成功】vasp测试K点收敛性的脚本问题
已经有20人回复
- 【求助】菜鸟问有关VASP作业提交的问题！
已经有21人回复
- 【求助】mpi vasp问题
已经有32人回复
- 【求助】vasp的计算时间含义
已经有12人回复
- 【求助】求将vasp的输出文件转化成xyz坐标的脚本
已经有11人回复
- 【分享】处理DOSCAR文件的脚本和程序
已经有19人回复
- 【求助】vasp运行问题
已经有7人回复
- 【求助】帮忙编个用于VASP小程序
已经有22人回复
yh3166
金虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 1207.2
- 红花: 2
- 帖子: 466
- 在线: 298.5小时
- 虫号: 1242579
- 注册: 2011-03-23
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
aminibuck: 金币+50, ★有帮助 2013-12-04 21:15:38
感谢参与，应助指数 +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
3楼2013-12-04 21:15:15
yh3166
金虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 1207.2
- 红花: 2
- 帖子: 466
- 在线: 298.5小时
- 虫号: 1242579
- 注册: 2011-03-23
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2013-12-05 08:49:30
djh434325996
木虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 4231.4
- 散金: 10
- 红花: 1
- 帖子: 268
- 在线: 116.6小时
- 虫号: 1829360
- 注册: 2012-05-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2013-12-05 08:57:46
6楼2013-12-05 09:32:44
7楼2013-12-05 09:33:44
yh3166
金虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 1207.2
- 红花: 2
- 帖子: 466
- 在线: 298.5小时
- 虫号: 1242579
- 注册: 2011-03-23
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
8楼2013-12-05 15:07:53
9楼2013-12-05 18:25:14
djh434325996
木虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 4231.4
- 散金: 10
- 红花: 1
- 帖子: 268
- 在线: 116.6小时
- 虫号: 1829360
- 注册: 2012-05-22
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
10楼2013-12-06 11:30:00
