请问如何用vasp计算Gd的单个原子能量?INCAR里面如何设置? [ Last edited by zeoliters on 2009-6-13 at 17:00 ] 返回小木虫查看更多
楼上回答的是对的。 结合能=晶体中每个原子的能量-孤立的单个原子能量。 我说的晶体中的单个原子能量,就是指晶体能量,但由于元胞设置问题,故而强调是平均到每个原子的能量。 你列出的公式也是对的。即晶体总能-孤立的单个原子能量×晶体中原子个数。显然这两个公式结果相差n倍,n=晶体中原子个数。 但一般所说的结合能的量纲是 能量/每原子。 DFT计算,有意义的不是能量绝对值而是相对的能量值。你算出来的晶体能量和孤立原子能量的值都不重要,关键是二者之差是否符合实验值。 [ Last edited by 老虎大王 on 2009-1-20 at 12:02 ],
老虎, 我计算的结果不怎么理想,我把我计算的结合能始终是负值。
如果我的理解没有错误,Ec的值应该就是一个负值,我们一般说的是Ec的绝对值。孤立的单个原子的能量应该比这个原子在晶体中的能量要高些的吧(即绝对值要小些,since它们都是负数)?
不管正负,那只是谁减谁的问题,只要绝对值对,应该就可以的吧。但我说过稀土元素是很难对付的。我估计要用全电子方法来计算。
请问计算单个原子能量时INCAR中如何设置?
谁能讲讲VASP自带的USPP,PAW,PBE 等赝势有何异同,在具体应用时以什么为选择依据? [ Last edited by ustbmars on 2009-2-18 at 16:11 ]
USPP是超软赝势, PAW投影缀加赝势 PBE 是交换关联相互作用GGA的一种,不是赝势。 另外,还有一种超软赝势是模守恒赝势。 我个人觉得paw赝势比较好, 但是往往不是很全, uspp也不错。 模守恒赝势稍显逊色
请问大家说的求体系的(例如楼上说的La7Gd1Mo8O36)的结合能或是形成能的话,需要先求单个原子的结合能而不仅是单个原子的能量能量吗?
楼上回答的是对的。
结合能=晶体中每个原子的能量-孤立的单个原子能量。
我说的晶体中的单个原子能量,就是指晶体能量,但由于元胞设置问题,故而强调是平均到每个原子的能量。
你列出的公式也是对的。即晶体总能-孤立的单个原子能量×晶体中原子个数。显然这两个公式结果相差n倍,n=晶体中原子个数。 但一般所说的结合能的量纲是 能量/每原子。
DFT计算,有意义的不是能量绝对值而是相对的能量值。你算出来的晶体能量和孤立原子能量的值都不重要,关键是二者之差是否符合实验值。
[ Last edited by 老虎大王 on 2009-1-20 at 12:02 ],
老虎, 我计算的结果不怎么理想,我把我计算的结合能始终是负值。
如果我的理解没有错误,Ec的值应该就是一个负值,我们一般说的是Ec的绝对值。孤立的单个原子的能量应该比这个原子在晶体中的能量要高些的吧(即绝对值要小些,since它们都是负数)?
不管正负,那只是谁减谁的问题,只要绝对值对,应该就可以的吧。但我说过稀土元素是很难对付的。我估计要用全电子方法来计算。
请问计算单个原子能量时INCAR中如何设置?
谁能讲讲VASP自带的USPP,PAW,PBE 等赝势有何异同,在具体应用时以什么为选择依据?
[ Last edited by ustbmars on 2009-2-18 at 16:11 ]
USPP是超软赝势,
PAW投影缀加赝势
PBE 是交换关联相互作用GGA的一种,不是赝势。
另外,还有一种超软赝势是模守恒赝势。
我个人觉得paw赝势比较好,
但是往往不是很全,
uspp也不错。
模守恒赝势稍显逊色
请问大家说的求体系的(例如楼上说的La7Gd1Mo8O36)的结合能或是形成能的话,需要先求单个原子的结合能而不仅是单个原子的能量能量吗?