要计算该能量主要为了计算体系的结合能，请高手详细告诉一下，我对这方面不太懂。谢谢

引用回帖:

Originally posted by wenxuan5788 at 2009-1-16 12:29:
请问如何用vasp计算Gd的单个原子能量？INCAR里面如何设置？

楼主说的不是很清楚，可以这么回答你，你要是想获得Gd组成某种结构时，单个原子的能量的话，那么你就将Gd按照某种结构进行建模，比如HCP，BCC， FCC等，然后计算总能量除以原子数目即可。
如果你想计算一个孤立原子的能量的话，那么你就要建立一个很大的包，只放一个原子于其中，比如10angstrom的大包。然后计算能量，建议看一下VASP的user guide 上面有。

我计算的体系是La7Gd1Mo8O36。我想计算一下体系的结合能，我看资料上说需要计算单个原子的能量。我计算Gd单个原子能量的步骤是：先建立一个大的单包：15X15X15。只放一个Gd原子，然后计算能量。我计算结果是-0.210eV，我觉得不太对，请问Gd，La的单个原子能量一般是多少？

你气死我了。我这很多天都在进行Gd的运算。用了好几个软件，都算不出来。可能是势函数不好。你用VASP居然一下子算出来了。我早就感到VASP可能有戏，但是我没有这个软件啊。

我来告诉你我的看法：
1.你光建立一个大单胞是不行的。你要建几个，比如15*15*15，20*20*20，25*25*25等，分别进行计算，看看能量是不是收敛到一个值，也就是说，看看元胞大到什么程度的时候，不影响计算结果。
2.元胞尺寸定下来以后，再考察K-point，以及能量截断的影响。办法就是多换几个参数去试，直到能量收敛（一般两次计算能量之差不超过1mhar），也就是这些参数不影响能量值为止。因为只算一个原子，应该很快的。
3.最后你应该得到一个你确信的值。

你问Ga、La单原子的能量是多少，这是没法回答你的。因为大家如果用不同的势或者不同的算法去算，结果绝对值可能有很大不同。你只能以如下方式检验你的结果：

A。计算Gd单个原子能量。
B。对Gd晶体进行结构优化（只优化元胞尺寸），并计算Gd晶体中单个原子能量。
C。求得Gd的结合能（cohesive energy)。
如果你的结合能算出来是对的，那么你的单原子能量就是可以使用的。

Gd的结合能Ecoh，文献值一般在4.14～4.18eV，
Gd的结构如下：
空间群：P6 3/mmC (194#)，HCP，Mg型，a=3.615ang，c=5.777ang。

以前也有人计算过如下结果：a=6.873Bohr，c=10.897Bohr， Ecoh＝4.344eV，给出的实验值则是a=10.902Bohr，c=6.858Bohr。
这篇文章是PRB， 50（1994)，1363。
还有一篇文章在PRB，44（1991），7451。你自己查吧。这两篇文章都是用LSDA算的。

对La和其它原子也是同样的处理和检验。对O，可以通过计算O2分子的原子化能来考察计算的有效性。

我以上说的是比较保险的办法，步步检验。当然这样搞的话比较慢一些，而且需要人很勤奋。呵呵

不过根据我的经验，稀土元素在DFT中是很难处理的。所以以上只供你参考。

如果你算出了好的结果，尤其是如果用GGA算出了好的结果，别忘了跟我联系一下，或者发个帖子讲一下你用的势和Gd的计算结果。谢谢。

[ Last edited by 老虎大王 on 2009-1-17 at 20:18 ]，

补充：
我以上说的是，在你没有你的体系（La7Gd1Mo8O36）的结合能实验值的情况下，你可以用这个方法来确立你对计算结果的信心。
如果你有La7Gd1Mo8O36的结合能实验数据，那就简单了，只需一股脑计算La、Gd、Mo、O的单原子能量，以及La7Gd1Mo8O36的晶体总能量就可以了，一下子算出结合能，跟实验值一比就行了。

请问INCAR里面该如何设置？我按照候博士讲义里的设置，计算出来的结果，但是感觉不太理想。“对Gd晶体进行结构优化（只优化元胞尺寸），并计算Gd晶体中单个原子能量。”你说的这句话我不太明白，晶体结合能还需要计算Gd晶体中单个原子能量吗？我看文献上说结合能的计算方法是：
晶体总能——（单个原子能量*原子个数）。

[ Last edited by wenxuan5788 on 2009-1-19 at 08:45 ]

引用回帖:

Originally posted by wenxuan5788 at 2009-1-19 08:25:
请问INCAR里面该如何设置？我按照候博士讲义里的设置，计算出来的结果，但是感觉不太理想。“对Gd晶体进行结构优化（只优化元胞尺寸），并计算Gd晶体中单个原子能量。”你说的这句话我不太明白，晶体结合能还需要 ...

我来帮着回答一下吧，对于Gd晶格进行结构优化是指的用HCP的结构进行Gd的基态总能量的优化计算，这个计算出来的是Gd在HCP结构下的基态总能量，然后要减去Gd的孤立原子能量，才能得到原子的结合能。
不清楚，再问吧