计算化学公社

标题: 算一个氢气分子的能量总是报错是为何？ [打印本页]

作者
Author: calvincui 时间: 2022-1-18 23:29
标题: 算一个氢气分子的能量总是报错是为何？
我在一个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

作者
Author: calvincui 时间: 2022-1-19 00:00
收敛了，得到的结果是E0= -.66922947E+01 请问各位大佬这个结果靠谱吗？

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

作者
Author: 乘风万里 时间: 2022-1-22 14:14
氢分子计算或优化时一旦打开ISPIN=2，需要设置NUPDOWN=0，要不就直接关闭自旋极化（ISPIN=1）。
详细参考VASP官方回复
以及researchgate

作者
Author: calvincui 时间: 2022-1-25 11:32

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

感谢大佬的分享！

作者
Author: calvincui 时间: 2022-1-25 11:32

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

感谢！

作者
Author: Vergelius 时间: 2022-10-17 01:13

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

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

作者
Author: Vergelius 时间: 2022-10-17 01:17

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

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

作者
Author: 乐平 时间: 2022-10-17 15:28

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

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)