| 查看: 1100 | 回复: 2 | ||
| 【奖励】 本帖被评价2次,作者zdhlover增加金币 1.5 个 | ||
| 当前主题已经存档。 | ||
[资源]
【转帖】如何用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 |
2楼2009-06-27 10:00:33
简单回复
aylayl083楼
2009-11-03 19:01
回复












回复此楼