| 查看: 1402 | 回复: 1 | ||
[求助]
求助：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等的生成热相差很大 希望有人提出看法，看看是否出现问题 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？
已经有20人回复
-
今天务委会开完了，明天出结果吗
已经有22人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有9人回复
-
明天应该可查了！？
已经有6人回复
-
有没有大神帮我看看基金代码
已经有10人回复
-
放榜前的不淡定
已经有19人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- QE计算在某一压强下的晶体参数要比相同条件下的VASP小的很多，是怎么回事呢？
已经有9人回复
- 求助：VASP单点能计算没有找到能量最低点
已经有6人回复
- 相同的vasp版本，不同机器上计算同一个体系，能量相差近10ev
已经有17人回复
- 关于HF方法和DFT方法中轨道能量的差别
已经有11人回复
- 在vasp计算里，有没有好的只优化原子位置和晶胞体积的方法？
已经有8人回复
- 【活动or资源】vasp晶体结构优化
已经有142人回复
- 为什么有限元法求出结构的固有频率通常比精确值大？
已经有9人回复
- 【求助】vasp 算出的0k下的能量如何扩展到其他温度压力
已经有12人回复
- 【求助】VASP里的几个能量怎么区别呢？？？？？？？？
已经有23人回复
- 【求助】vasp静态计算能量值错误
已经有21人回复
- 【求助】怎样用vasp正确的做出晶体场中的t2g-eg轨道PDOS?
已经有23人回复
- 【求助】计算得到能量后，参考点是什么呢？
已经有6人回复
- 【求助】VASP怎样精确计算出费米能级处的态密度值？
已经有22人回复
- 【求助】算VASP弹性常数时的能量-应变出错
已经有18人回复
- 【求助】vasp 离子驰豫与电子自旋的关系
已经有11人回复
- 【求助】能否发一份用vasp计算超晶胞某个表面的输入文件给我
已经有7人回复
- 【求助】Castep中结构优化后晶胞参数差别较大 vasp高手也帮我看一下
已经有22人回复
- 【转帖】用VASP计算氢分子的能量
已经有8人回复
eqianlicao
铜虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 475.6
- 红花: 1
- 帖子: 441
- 在线: 116.2小时
- 虫号: 1305350
- 注册: 2011-05-24
- 性别: GG
- 专业: 金属结构材料
2楼2019-07-22 09:54:57
