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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » encut测试是否需要修改要提交的作业脚本
51/1返回列表
查看: 1382  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-04-22 16:37:10
已阅   回复此楼   关注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
赞一下(2人) 回复此楼
天行健君子以自强不息
5楼2011-05-10 11:04:26
已阅   回复此楼   关注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='脚本位置”
赞一下(1人) 回复此楼
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的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 304求调剂 +7 castLight 2026-04-16 7/350 2026-04-17 20:05 by 关一盏灯cd
[考研] 一志愿沪9，326求生物学调剂 +11 刘墨墨 2026-04-13 11/550 2026-04-17 15:35 by wanganpin
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
[考研] 279求调剂 +13 张番茄不炒蛋 2026-04-11 13/650 2026-04-17 10:38 by cuisz
[考研] 22专硕求调剂 +10 haoyun上岸 2026-04-11 12/600 2026-04-16 22:21 by 猪会飞
[考研] 化学070300 求调剂 +28 哈哈哈^_^ 2026-04-12 28/1400 2026-04-16 21:36 by 大力水手力大无�
[考研] 294求调剂 +14 淡然654321 2026-04-15 14/700 2026-04-16 21:01 by lpl364211
[考研] 求调剂推荐 +8 小聂爱学习 2026-04-14 8/400 2026-04-16 07:22 by 学员JpLReM
[考研] 085404 22408 309分求调剂 +9 lzmk 2026-04-14 10/500 2026-04-15 20:02 by 学员JpLReM
[考研] 生物学调剂 +9 纸扇zhishan 2026-04-13 9/450 2026-04-15 18:28 by AN流800
[考研] 药学305求调剂 +7 玛卡巴卡boom 2026-04-11 7/350 2026-04-15 13:21 by 西北望—风沙
[考研] 求调剂 +12 何气正 2026-04-13 13/650 2026-04-14 14:47 by zs92450
[考研] 085408光电信息工程专硕355一志愿长春光机所调剂 +6 王ymaa 2026-04-13 13/650 2026-04-14 11:33 by 王ymaa
[考研] 求调剂 +3 我爱高数高数爱� 2026-04-12 3/150 2026-04-14 01:00 by 王珺璞
[考研] 材料考研调剂 +29 云木达达 2026-04-11 31/1550 2026-04-13 13:32 by lyh鲁老师
[考研] 一志愿085802 323分求调剂 +13 drizzle_9 2026-04-12 14/700 2026-04-13 10:26 by Faiz5552
[考研] 0831一轮调剂失败求助 +10 小熊睿睿_s 2026-04-11 10/500 2026-04-12 22:43 by 长弓傲
[考研] 调剂结束 +6 floriea 2026-04-12 8/400 2026-04-12 18:13 by zhouxiaoyu
[考研] 291求调剂 +8 关忆北. 2026-04-11 8/400 2026-04-12 09:32 by 逆水乘风
[考研] 359求调剂 +5 胃痉挛累了 2026-04-11 5/250 2026-04-11 19:55 by lbsjt
信息提示
关闭
请填处理意见
关闭