The ABINIT Users Mailing List ( CLOSED )

[myid520@sohu.com]

我现在正在使用ABINIT做计算，使用时遇到了一个问题，或者说是一类问题--算金属表面能使遇到的问题，用ABINIT所给的教程验证，这个问题也同样存在。以教程第四课为例，问题是这样的：教程第三节得到的包含单个原子的原胞的总能量为：etotal24=-2.0979739819E+00
 Ha（下面用etot24表示），第四节得到包含两个原子的超晶胞的总能量为：-4.1962972596E+00 
ha,这几乎正好是包含单个原子的原胞的总能量的2倍（相差小于0.0005ha），然而，当我以t44.in为基础，把rprim改为 rprim  0.5 
-0.5  0.0  ，0.5  0.5  0.0， 0.0  0.0  1.5原子坐标改为xcart  3*0.0 ，0.0  2*3.75279545 
， 2*0.0  7.5055909  来计算包含三个原子的超晶胞的总能时，它的总能量（6.224191304 
ha）和etot24的三倍相差0.0697ha,相差达1.12%，继续改变rprim和xcart,增加超晶胞包含的原子的个数到4和5，总能量增为8.394872163和10.422060491，这和etot24的4倍和5倍分别相差0.003ha,0.0678ha和0.035%，0.65%，继续增加下去，包含奇数个原子的超晶格的总能量一直相差0.06ha左右，这么大的差别我认为是不应该有的，至于为什么会产生这么大的差别，是我得原因，还是程序的原因，我百思不得其解，特向您请教，望不吝赐教