24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1372  |  回复: 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的回帖

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的回帖
查看全部 5 个回答

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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 人工智能320调剂08工类还有机会吗 +6 振—TZ 2026-04-10 6/300 2026-04-10 15:19 by zl8213662
[考研] 材料复试求调剂 +20 xhhdjdjsjks 2026-04-09 20/1000 2026-04-10 10:25 by 孙小小12457
[考研] 调剂 +10 卷卷卷心菜_ 2026-04-09 10/500 2026-04-10 09:53 by Delta2012
[考研] 070300化学学硕311分求调剂 +12 梁富贵险中求 2026-04-04 14/700 2026-04-10 08:07 by kangsm
[考研] 已调剂 +18 柴郡猫_ 2026-04-09 19/950 2026-04-09 22:10 by 柴郡猫_
[考研] 284求调剂 +7 archer.. 2026-04-09 8/400 2026-04-09 21:30 by bljnqdcc
[考研] 070300化学279求调剂 +17 哈哈哈^_^ 2026-04-08 18/900 2026-04-09 10:49 by 三七七想上岸
[考研] 0703调剂,一志愿天津大学319分 +23 haaaabcd 2026-04-05 26/1300 2026-04-08 16:19 by luoyongfeng
[考研] 一志愿南昌大学,085600,344分求调剂 +11 调剂上岸玘 2026-04-05 12/600 2026-04-08 16:17 by luoyongfeng
[考研] 专硕085403,291分,有两篇专利,一国一奖 +3 哈吉咪哈吉咪 2026-04-07 3/150 2026-04-07 18:21 by 蓝云思雨
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +11 大火山小火山 2026-04-05 11/550 2026-04-06 22:55 by yunlongyang
[考研] 工科 22408 267求推荐 +4 wanwan00 2026-04-05 5/250 2026-04-06 22:47 by chenzhimin
[考研] 297分083200求助 +9 aekx 2026-04-05 9/450 2026-04-06 20:57 by flysky1234
[考研] 377求调剂 +6 by.ovo 2026-04-05 6/300 2026-04-05 22:18 by dongzh2009
[考研] 求调剂 +7 张.1 2026-04-05 7/350 2026-04-05 20:40 by 啵啵啵0119
[考研] 求调剂 +10 Hll胡 2026-04-04 10/500 2026-04-05 20:09 by nepu_uu
[考研] 298求调剂 +7 manman511 2026-04-05 7/350 2026-04-05 10:29 by 唐沐儿
[考研] 能动调剂326专硕 +4 wan112233 2026-04-04 4/200 2026-04-04 22:47 by yu221
[考研] 考研调剂 +5 四川王涛 2026-04-04 5/250 2026-04-04 22:18 by 啵啵啵0119
[考研] 求调剂 +4 压力??大 2026-04-03 4/200 2026-04-03 21:36 by 啵啵啵0119
信息提示
请填处理意见