算一个氢气分子的能量总是报错是为何？

calvincui

我在一个17*17*17的真空盒子中放入一个H2分子，结构优化。结果总是报错Call to ZHEGV failed. Returncode = 11 1 12
以下是我的INCAR POSCAR KPOINTS

INCAR:

ISTART=0
ICHARG=2
ISPIN=2
#MAGMOM=2*0
PREC=Accurate
LREAL=.FALSE.

ALGO=N
NELM=100
EDIFF=1E-5
ENCUT=250
NCORE=12

IBRION=2
NSW=100
ISIF=2
IVDW=11
EDIFFG=-0.03

AMIX=0.05
#BMIX=0.0001

ISMEAR=0
SIGMA=0.1


POSCAR:

H                                      
1.0
        17.000000000         0.0000000000         0.0000000000
        0.0000000000         17.000000000         0.0000000000
        0.0000000000         0.0000000000         17.000000000
   H
   2
Cartesian
     0.000000000         0.000000000         0.000000000
     0.000000000         0.000000000         0.760000000


KPOINTS:

KPT-Resolved Value to Generate K-Mesh: 0.030
0
Gamma
   1   1   1
0.0  0.0  0.0

calvincui

收敛了，得到的结果是E0= -.66922947E+01 请问各位大佬这个结果靠谱吗？

wypkdhd

算气体分子格子不要是正方形。另，能量毫无意义，相对能量才有意义。比如你从nist看到氢气氧气生成水的自由能是多少，你把这三个分子算一下，能量矫正后相减，如果差不多，说明你算的没问题。

乘风万里

氢分子计算或优化时一旦打开ISPIN=2，需要设置NUPDOWN=0，要不就直接关闭自旋极化（ISPIN=1）。
详细参考VASP官方回复
以及researchgate

calvincui

乘风万里 发表于 2022-1-22 14:14
氢分子计算或优化时一旦打开ISPIN=2，需要设置NUPDOWN=0，要不就直接关闭自旋极化（ISPIN=1）。
详细参考V ...

感谢大佬的分享！

calvincui

wypkdhd 发表于 2022-1-19 19:45
算气体分子格子不要是正方形。另，能量毫无意义，相对能量才有意义。比如你从nist看到氢气氧气生成水的自由 ...

感谢！

Vergelius

calvincui 发表于 2022-1-19 00:00
**** 作者被禁止或删除 内容自动屏蔽 ****

我算出来-6.75832177eV，差不多吧

Vergelius

wypkdhd 发表于 2022-1-19 19:45
算气体分子格子不要是正方形。另，能量毫无意义，相对能量才有意义。比如你从nist看到氢气氧气生成水的自由 ...

请问为何不能用正方形的格子呢？应该用什么样的格子呢？

乐平

Vergelius 发表于 2022-10-17 01:17
请问为何不能用正方形的格子呢？应该用什么样的格子呢？

不用立方体的格子是为了避免由于对称性引起的能量不准确。一般不对称的盒子得到的能量会更低一点。

可以参考大师兄网页里的例子
https://www.bigbrosci.com/2017/10/30/ex08/