| 查看: 1436 | 回复: 0 | |||
[交流]
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附近明显不对，不知道是哪儿出了问题，请各位讨论。 |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有34人回复
-
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有78人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有12人回复
-
FileCode能看出啥？
已经有30人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有15人回复
-
关于Filecode分析方法
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有10人回复
-
阿姨
已经有4人回复
-
一作与独作在应聘高校教师时区别大吗
已经有4人回复
-
静等基金结果
已经有21人回复
» 本主题相关商家推荐: (我也要在这里推广)
