24小时热门版块排行榜

返回列表

baoqinfu

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 160.6
散金: 457
红花: 1
帖子: 98
在线: 40.3小时
虫号: 1551666
注册: 2011-12-26
专业: 极端条件下使用的金属材料

[交流] VASP扫描计算H2分子的能量

我感觉我这样计算好像很有问题，请问大家是怎么回事？
我是通过固定两个H原子之间的距离，计算这样体系（盒子是30x30x30）的总能量，不知道为什么计算的结果感觉不对，在0.74埃附近应该是H2分子最稳定的,但能量却是大的，这到底是怎么回事？
下面是INCAR (基本上与vasp手册上静态能计算推荐的设置一致):
SYSTEM = H2 dimer in a box
ISTART = 0
ICHARG = 2
INIWAV = 1
ISMEAR = 0
SIGMA = 0.05
ISPIN = 2
IBRION = -1
ENCUT = 350
NELM = 50
EDIFF = 10E-4
NELMIN = 8
# POTIM = 0.01
-----------------------------------------------------------------------------
下面是KPOINTS
Gamma-point only
0
Monkhorst Pack
1    1    1
0    0    0
-----------------------------------------------------------------------------
脚本大概如下：
for i in 0.5 0.6 0.65 0.68 0.69 0.7 0.71 0.72 0.73 0.74 0.75 0.76 0.77 0.78 0.79 0.8 0.82 0.85 0.9 1.0 1.2 1.4 1.7 2.0 2.3 2.6 2.9 3.1 3.3 3.5 3.7 3.9 4.1 4.3 4.5 4.7 4.9 5.1 5.5
do
cat >POSCAR <<!
H2
1.0
30.0 0.0    0.0
0.0 30.0 0.0
0.0 0.0    30.0
H
2
Cartesian
0.0  0.0  0.0
0.0  0.0  $i
!
rm C* D* E* IB* O* PC* R* v* W* X*
echo "a= $i"
./run_vasp
date >> H2spin.log

E=`awk '/F=/ {print $0}' OSZICAR`
echo $i $E  >>SUMMARY
done
cat SUMMARY
------------------------------------------------------------------------------------------------
结果是
0.5
0.6
0.65 1 F= 0.71932739E+01 E0= 0.71932739E+01 d E =-.989847E-15 mag= -0.3223
0.68 1 F= 0.46978484E+01 E0= 0.46978484E+01 d E =-.281632E-15 mag= 2.0089
0.69 1 F= 0.42175507E+01 E0= 0.42175507E+01 d E =-.381031E-11 mag= 2.0115
0.7 1 F= 0.40547571E+01 E0= 0.40547571E+01 d E =-.135693E-10 mag= 2.0179
0.71 1 F= 0.38644594E+01 E0= 0.38644594E+01 d E =-.167107E-11 mag= 1.9952
0.72 1 F= 0.37176841E+01 E0= 0.37176841E+01 d E =-.140140E-10 mag= 1.9950
0.73 1 F= 0.35638797E+01 E0= 0.35638797E+01 d E =-.210662E-14 mag= 1.9952
0.74 1 F= 0.34831807E+01 E0= 0.34831807E+01 d E =-.472721E-15 mag= 2.0096
0.75 1 F= 0.33228758E+01 E0= 0.33228758E+01 d E =-.336994E-14 mag= 1.9957
0.76 1 F= 0.37295166E+01 E0= 0.37295166E+01 d E =-.615765E-11 mag= 1.9932
0.77 1 F= 0.33996558E+01 E0= 0.33996558E+01 d E =-.141141E-13 mag= 1.9990
0.78 1 F= 0.31511279E+01 E0= 0.31511279E+01 d E =-.121456E-13 mag= 1.9930
0.79 1 F= 0.27768658E+01 E0= 0.27768658E+01 d E =-.252151E-11 mag= 2.0029
0.8 1 F= 0.26422858E+01 E0= 0.26422858E+01 d E =-.490012E-15 mag= 1.9967
0.82 1 F= 0.24182024E+01 E0= 0.24182024E+01 d E =-.120547E-14 mag= 1.9974
0.85 1 F= 0.20995430E+01 E0= 0.20995430E+01 d E =-.852705E-15 mag= 2.0000
0.9 1 F= 0.16261543E+01 E0= 0.16261543E+01 d E =-.173925E-11 mag= 2.0000
1.0 1 F= 0.83133423E+00 E0= 0.83133423E+00 d E =-.975193E-11 mag= 2.0000
1.2 1 F= -.30706131E+00 E0= -.30706131E+00 d E =-.460604E-14 mag= 2.0000
1.4 1 F= -.10368067E+01 E0= -.10368067E+01 d E =-.928676E-11 mag= 2.0000
1.7 1 F= -.16641998E+01 E0= -.16641998E+01 d E =-.364711E-11 mag= 2.0000
2.0 1 F= -.19729514E+01 E0= -.19729514E+01 d E =-.934003E-11 mag= 2.0000
2.3 1 F= -.21190846E+01 E0= -.21190846E+01 d E =-.482783E-12 mag= 2.0000
2.6 1 F= -.21867386E+01 E0= -.21867386E+01 d E =-.293903E-11 mag= 2.0000
2.9 1 F= -.22070472E+01 E0= -.22070472E+01 d E =-.484678E-12 mag= 2.0000
3.1 1 F= -.22171237E+01 E0= -.22171237E+01 d E =-.546002E-12 mag= 2.0000
3.3 1 F= -.22229638E+01 E0= -.22229638E+01 d E =-.186034E-11 mag= 2.0000
3.5 1 F= -.22356005E+01 E0= -.22356005E+01 d E =-.345572E-11 mag= 2.0000
3.7 1 F= -.22370247E+01 E0= -.22370247E+01 d E =-.160644E-11 mag= 2.0000
3.9 1 F= -.22374602E+01 E0= -.22374602E+01 d E =-.131266E-10 mag= 2.0000
4.1 1 F= -.22374297E+01 E0= -.22374297E+01 d E =-.469901E-11 mag= 2.0000
4.3 1 F= -.22370057E+01 E0= -.22370057E+01 d E =-.909261E-11 mag= 2.0000
4.5 1 F= -.22365721E+01 E0= -.22365721E+01 d E =-.196378E-12 mag= 2.0000
4.7 1 F= -.22360999E+01 E0= -.22360999E+01 d E =-.376436E-11 mag= 2.0000
4.9 1 F= -.22356867E+01 E0= -.22356867E+01 d E =-.492664E-12 mag= 2.0000
5.1 1 F= -.22352891E+01 E0= -.22352891E+01 d E =-.922265E-13 mag= 2.0000
5.5 1 F= -.22345708E+01 E0= -.22345708E+01 d E =-.412528E-11 mag= 2.0000
0.5 1 F= 0.33525596E+01 E0= 0.33525596E+01 d E =-.139746E-10 mag= 1.1364
0.6
0.65 1 F= 0.17727474E+02 E0= 0.17727474E+02 d E =-.150644E-11 mag= 0.4014
0.68 1 F= 0.47059381E+01 E0= 0.47059381E+01 d E =-.338158E-12 mag= 2.0069
0.69 1 F= 0.42175507E+01 E0= 0.42175507E+01 d E =-.381031E-11 mag= 2.0115
0.7 1 F= 0.40550365E+01 E0= 0.40558384E+01 d E =-.160378E-02 mag= 2.0161
0.71 1 F= 0.38643956E+01 E0= 0.38643956E+01 d E =-.324766E-11 mag= 1.9948
0.72 1 F= 0.37176841E+01 E0= 0.37176841E+01 d E =-.140140E-10 mag= 1.9950
0.73 1 F= 0.35638797E+01 E0= 0.35638797E+01 d E =-.210662E-14 mag= 1.9952
0.74 1 F= 0.35105921E+01 E0= 0.35105921E+01 d E =-.127873E-12 mag= 2.0060
0.75 1 F= 0.33360920E+01 E0= 0.33360920E+01 d E =-.104135E-14 mag= 2.0070
0.76 1 F= 0.39488110E+01 E0= 0.39488110E+01 d E =-.182635E-14 mag= 1.9975
0.77 1 F= 0.32917067E+01 E0= 0.32917067E+01 d E =-.698864E-14 mag= 1.9986
0.78 1 F= 0.29468898E+01 E0= 0.29468898E+01 d E =-.208986E-15 mag= 1.9913
0.79 1 F= 0.27881217E+01 E0= 0.27881217E+01 d E =-.332189E-14 mag= 1.9989
0.8 1 F= 0.26421814E+01 E0= 0.26421814E+01 d E =-.884492E-14 mag= 1.9971
0.82 1 F= 0.24196240E+01 E0= 0.24196240E+01 d E =-.109690E-13 mag= 1.9978
0.85 1 F= 0.21172659E+01 E0= 0.21172659E+01 d E =-.303546E-15 mag= 2.0000
0.9 1 F= 0.16438610E+01 E0= 0.16438610E+01 d E =-.395586E-12 mag= 2.0000
----------------------------------------------------------------------------------------------------------------
这个在0.74附近明显不对，不知道是哪儿出了问题，请各位讨论。

回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

1楼 2017-09-23 10:30:42

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

相关版块跳转我要订阅楼主 baoqinfu 的主题更新

返回列表