青青杨柳分享 http://blog.sciencenet.cn/u/Irene

博文

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

已有 15667 次阅读 2009-4-13 15:30 |个人分类:生活点滴|系统分类:科研笔记|关键词:学者| 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 
 

 
 原文地址 http://new.quantumchemistry.net/Experience/CommonSoftwares/VASP/ConcereteExample/200512/250.html  


https://m.sciencenet.cn/blog-240700-225928.html

上一篇:VASP电荷密度差的计算
下一篇:关于VASP中光学optics.F程序的编译(ifort+mkl)

0

该博文允许实名用户评论 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-6-16 01:11

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部