24小时热门版块排行榜    

查看: 1318  |  回复: 1

long1990

铜虫 (初入文坛)

[求助] 求助:vasp计算对晶体结构结构优化时能量可能偏高

INCAR.relax
SYSTEM = Nb8H-bcc-oct
ENCUT = 350
ISTART = 0
ICHARG = 2
ISMEAR = 1
SIGMA = 0.2
NSW = 100
IBRION = 2
ISIF = 4
POTIM = 0.2
EDIFF = 1E-6
EDIFFG = -1E-3
PREC = Accurate

KPOINTS.relax
Automatic generation
0
Monkhorst Pack
6  6  11
0  0  0

INCAR.static
SYSTEM = Nb8H-bcc-oct
ENCUT = 350
ISTART = 0
ICHARG = 2
ISMEAR = -5
PREC = Accurate

KPOINTS.static
Automatic generation
0
Monkhorst Pack
7  7  13
0  0  0


运行程序:
#PBS -N vasp
#PBS -l nodes=1:ppn=8
#PBS -j oe
#PBS -q gonghr
#PBS -l walltime=1000:00:00
cd $PBS_O_WORKDIR
NP=`cat $PBS_NODEFILE|wc -l`

for a in 3.30 3.31 3.32 3.33 3.34 3.35 3.36 3.37 3.38 3.39 3.40
do
cat > POSCAR <<!
bcc
$a
2.0 0.0 0.0
0.0 2.0 0.0
0.0 0.0 1.0
8 1
Direct
0.0 0.0 0.0
0.25 0.25 0.5
0.5 0.0 0.0
0.75 0.25 0.5
0.0 0.5 0.0
0.25 0.75 0.5
0.5 0.5 0.0
0.75 0.75 0.5
0.25 0.25 0.0
!
cp INCAR.relax INCAR
cp KPOINTS.relax KPOINTS
mpirun -np $NP -machinefile $PBS_NODEFILE --mca btl self,openib /public/home/fenmo/tahitiwinny/VASP5.2/vasp.5.2/vasp
cp CONTCAR POSCAR
cp CONTCAR CONTCAR-$a
cp OSZICAR OSZICAR-$a
cp INCAR.static INCAR
cp KPOINTS.static KPOINTS
mpirun -np $NP -machinefile $PBS_NODEFILE --mca btl self,openib /public/home/fenmo/tahitiwinny/VASP5.2/vasp.5.2/vasp
E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'`
V=`grep "volume" OUTCAR | tail -1 | awk '{printf "%12.4f \n", $5 }'`
echo $a $E >>aE
echo $V $E >>EtVo.dat
done


EtVo.dat

143.7500 -83.482626
145.0600 -83.543300
146.3800 -83.589522
147.7000 -83.621667
149.0400 -83.639831
150.3800 -83.644377
151.7300 -83.636925
153.0900 -83.616228
154.4600 -83.584004
155.8300 -83.539846
157.2200 -83.484020

感觉能量整体相对偏大,比如后来得到的生成热相对Nb4H-bcc-oct等的生成热相差很大
希望有人提出看法,看看是否出现问题
回复此楼

» 猜你喜欢

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

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

eqianlicao

铜虫 (正式写手)

CSU gonghaoran组的?
超强铝合金
2楼2019-07-22 09:54:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 long1990 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 293求调剂 +9 zjl的号 2026-03-16 14/700 2026-03-18 05:13 by qytt腾飞
[考研] 材料专硕306英一数二 +8 z1z2z3879 2026-03-16 10/500 2026-03-18 00:01 by 星空星月
[考研] 0703化学调剂 ,六级已过,有科研经历 +8 曦熙兮 2026-03-15 8/400 2026-03-17 20:31 by xilongliang
[考研] 材料与化工求调剂 +6 为学666 2026-03-16 6/300 2026-03-17 20:15 by peike
[考研] 299求调剂 +4 △小透明* 2026-03-17 4/200 2026-03-17 20:09 by peike
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 求调剂,总分315,考的生物医药,一志愿湖南师范大学。调剂到任何专业都可以 +4 小丁想进步 2026-03-11 5/250 2026-03-17 16:05 by 外星文明
[考研] 211本,11408一志愿中科院277分,曾在中科院自动化所实习 +6 Losir 2026-03-12 7/350 2026-03-17 12:09 by danranxie
[考研] 274求调剂 +5 时间点 2026-03-13 5/250 2026-03-17 07:34 by 热情沙漠
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 求调剂(材料与化工327) +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 290求调剂 +9 ADT 2026-03-11 9/450 2026-03-13 21:55 by JourneyLucky
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] 材料与化工085600调剂求老师收留 +9 jiaanl 2026-03-11 9/450 2026-03-13 20:22 by JourneyLucky
[考研] 301求调剂 +6 Liyouyumairs 2026-03-11 6/300 2026-03-13 20:11 by JourneyLucky
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
信息提示
请填处理意见