24小时热门版块排行榜    

查看: 1335  |  回复: 4
本帖产生 2 个 1ST强帖 ,点击这里进行查看

dingfengbo

新虫 (小有名气)

[求助] encut测试是否需要修改要提交的作业脚本

我想做encut测试这是我的encut脚本:
#!/bin/sh
rm WAVECAR
for a in 400 450 500 550 600 650 700 750
do
cat >INCAR < SYSTEM=CGa
LPLANE=.TRUE.
NPAR=8
LREAL = .TRUE.
PREC    = Accurate
EDIFF   =1E-6
EDIFFG  = -0.03      
IALGO   = 48         
NELMIN  = 4         
ISYM    = 0         
ISPIN = 1
NBANDS = 240
LCHARG  = .TRUE.
LWAVE   = .TRUE.     
LVTOT   = .FALSE.
NBLOCK  = 1         
NSW     = 10
IBRION  = 2
POTIM=0.3
ISMEAR = 0   
SIGMA = 0.05
ENCUT=$a
ISTART=0;ICHARG=2
!
echo "ENCUT=$a eV"
vasp
E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'`
echo $a $E >>comment
done
测试的时候我想将之提交,下面是我的作业脚本:
#!/bin/sh
#PBS -l nodes=2:ppn=8
#PBS -u XX
#PBS -N YY

JOBINFOR=$PBS_JOBID
MASTERNODE=`hostname`
SCRATCHDIR=$PBS_JOBID
NCPU=`wc -l < $PBS_NODEFILE`

SERVER=$PBS_O_HOST
WORKDIR=$PBS_O_WORKDIR
MKDIR=/bin/mkdir
RSH=/usr/bin/rsh
CP=/bin/cp
cat $PBS_NODEFILE > $WORKDIR/mpihosts
LAUNCH="mpirun --hostfile  $WORKDIR/mpihosts -np $NCPU"
PROGRAMEXEC="/opt/vasp/vasp"

${RSH} ${MASTERNODE}
cd ${WORKDIR}
${MKDIR} ${SCRATCHDIR}
cd ${SCRATCHDIR}
${CP} ${WORKDIR}/* .
${LAUNCH}  ${PROGRAMEXEC} > out 2>&1
请问我怎样修改作业脚本才能进行encut测试的计算呢?谢谢
回复此楼

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

Photochemistry

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)

【答案】应助回帖


youzhizhe(金币+1, 1ST强帖+1): 谢谢指导。 2011-04-23 08:13:11
如果知识测试小体系的话,建议不要提交,直接单机运行。如果是大体系的话,估计一个ecut就可以了,如果全部是金属,400eV足够;如果有O,N可能要700eV。其他的都可以采用赝势中默认值。
如果一定要提交,可以考虑把这句
PROGRAMEXEC="/opt/vasp/vasp"
改为
PROGRAMEXEC='脚本位置”
2楼2011-04-22 22:58:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

haowenping

木虫 (著名写手)

引用回帖:
Originally posted by dingfengbo at 2011-04-22 16:37:10:
我想做encut测试这是我的encut脚本:
#!/bin/sh
rm WAVECAR
for a in 400 450 500 550 600 650 700 750
do
cat >INCAR <<!
SYSTEM=CGa
LPLANE=.TRUE.
NPAR=8
LREAL = .TRUE.
PREC    = Accura ...

你好,请问你有没有实现在pbs下优化你的那个截断能脚本?我现在也被这个问题困惑,如果你会了的话麻烦指教下我,谢谢了!
人有多大胆,地有多大产!
3楼2011-05-07 19:07:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingfengbo

新虫 (小有名气)

引用回帖:
Originally posted by haowenping at 2011-05-07 19:07:03:
你好,请问你有没有实现在pbs下优化你的那个截断能脚本?我现在也被这个问题困惑,如果你会了的话麻烦指教下我,谢谢了!

你好!按照2楼说的就可以实现了,不知道你想问什么?
4楼2011-05-08 12:52:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiao72379

木虫 (正式写手)

宅男甲

【答案】应助回帖


bingmou(金币+1, 1ST强帖+1): 感谢交流! 2011-05-10 16:52:58
引用回帖:
Originally posted by bingmou at 2011-04-22 22:58:54:
如果知识测试小体系的话,建议不要提交,直接单机运行。如果是大体系的话,估计一个ecut就可以了,如果全部是金属,400eV足够;如果有O,N可能要700eV。其他的都可以采用赝势中默认值。
如果一定要提交,可以考虑 ...

个人意见:PROGRAMEXEC="/opt/vasp/vasp" 改为 PROGRAMEXEC='脚本位置”是不行的,最后一行其实就是一个运行并行的vasp程序的指令。如果改成脚本位置,那运行的程序也就是vasp都不能调用。只是楼主注意下变量PROGRAMEXEC="/opt/vasp/vasp"设置在你的机器上是正确的就行。还有就是我觉得这个脚本写得罗里罗嗦, 我整合个简单点的 看着清楚的 供楼主参考下:
#!/bin/bash
#PBS -N 5-zmos2nrbig
#PBS -j oe
#PBS -q vasp
#PBS -l walltime=9999:00:00
#PBS -l nodes=2:ppn=8

rm WAVECAR
for a in 400 450 500 550 600 650 700 750
do
cat >INCAR < SYSTEM=CGa
LPLANE=.TRUE.
NPAR=8
LREAL = .TRUE.
PREC    = Accurate
EDIFF   =1E-6
EDIFFG  = -0.03      
IALGO   = 48         
NELMIN  = 4         
ISYM    = 0         
ISPIN = 1
NBANDS = 240
LCHARG  = .TRUE.
LWAVE   = .TRUE.     
LVTOT   = .FALSE.
NBLOCK  = 1         
NSW     = 10
IBRION  = 2
POTIM=0.3
ISMEAR = 0   
SIGMA = 0.05
ENCUT=$a
ISTART=0;ICHARG=2
!
echo "ENCUT=$a eV"
mpirun XXX(此处是mpi的设置参数,根据自己的机器设置) ./vasp > ./out(此处将vasp程序 vasp必须的文件和脚本放同一目录,要嫌麻烦就在脚本里加几个变量把路径写进去)
E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'`
echo $a $E >>comment
done
天行健君子以自强不息
5楼2011-05-10 11:04:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dingfengbo 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856调剂,是学校就去 +7 sllhht 2026-03-19 8/400 2026-03-20 09:30 by 每天只摆一小会
[考研] 265求调剂 +6 梁梁校校 2026-03-17 6/300 2026-03-20 08:59 by ZHANG0tao
[考研] 307求调剂 +9 冷笙123 2026-03-17 9/450 2026-03-19 22:44 by 学员8dgXkO
[考研] 294求调剂材料与化工专硕 +14 陌の森林 2026-03-18 14/700 2026-03-19 22:38 by 学员8dgXkO
[考研] 0703化学调剂 ,六级已过,有科研经历 +12 曦熙兮 2026-03-15 12/600 2026-03-19 19:42 by maocaozhuxi
[考博] 东华理工大学化材专业26届硕士博士申请 +8 zlingli 2026-03-13 8/400 2026-03-19 16:32 by 轻松不少随
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 10/500 2026-03-19 15:08 by 阳阳哇塞
[考研] 324分 085600材料化工求调剂 +3 llllkkkhh 2026-03-18 3/150 2026-03-19 14:22 by houyaoxu
[考研] 一志愿中海洋材料工程专硕330分求调剂 +7 小材化本科 2026-03-18 7/350 2026-03-19 10:46 by Linda Hu
[教师之家] 焦虑 +9 水冰月月野兔 2026-03-13 13/650 2026-03-19 09:50 by otani
[考研] 一志愿华中科技大学,080502,354分求调剂 +4 守候夕阳CF 2026-03-18 4/200 2026-03-18 22:16 by li123456789.
[考研] 085601专硕,总分342求调剂,地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 070300化学319求调剂 +6 锦鲤0909 2026-03-17 6/300 2026-03-18 13:22 by Iveryant
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考研] 278求调剂 +5 烟火先于春 2026-03-17 5/250 2026-03-18 08:43 by 星空星月
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
信息提示
请填处理意见