24小时热门版块排行榜

返回列表

【奖励】本帖被评价2次，作者zdhlover增加金币 1.5 个

当前主题已经存档。

zdhlover

荣誉版主 (职业作家)

应助: 0 (幼儿园)
贵宾: 6.748
金币: 19449.4
帖子: 3105
在线: 119.6小时
虫号: 409529

[资源] 【转帖】如何用VASP计算单个原子的能量[权威版]

[分享] 如何用VASP计算单个原子的能量[权威版]

这是VASP作者给出的解答,所以算是权威了.原文如下:

  The energy zero
  ..
  in VASP all energies are referred to the the reference state for which
  the potential was generated!
  this is in most cases not the real groundstate of the atom
  ..
  to determined the energy of the grounstate of the atom place the atom
  in a larger non cubic box to break initial symmetry (i.e.
  11Ax10Ax9A)
  use the G point only
  INCAR:
  ISPIN = 2 ! spin polarized
  ISMEAR = 0 ; SIGMA = 0.2 ! for small sigma conv. for TM is diff.
  MAGMOM = 2 ! initial magnetic moment
  one should use the energy value energy without entropy of the OUTCAR file
  since this coverges most rapidly to the correct energy for sigma 0
  ..

Ecoh=E_total-n*Eatom

  大约翻译如下:
  在VASP中，所有的能量是相对于产生赝势的组态来说的，所以并不能将计算得到的能量当成真正的原子的基态能量。

  为了得到原子的基态能量，将原子放在一个非正方体的晶格内，比如说11Ax10Ax9A的晶格（这样做是为了消除简并）。
  KPOINTS只用G撒点
  大概的INCAR文件可以如下：
  ISPIN = 2 ! spin polarized
  ISMEAR = 0 ; SIGMA = 0.2 ! for small sigma conv. for TM is diff.
  MAGMOM = 2 ! initial magnetic moment

  计算结束后，应该用OUTCAR中不含entropy的那一项能量，因为它通常更快的收敛到sigma=0的正确值。最后计算结合能的公式应该为

Ecoh=E_total-n*Eatom

回复此楼