为什么没有人回答呢？请高手帮忙讲解一下啊 谢谢

考虑自旋的一般要修正吧！

引用回帖:

Originally posted by nyoxd at 2010-04-27 11:51:12:
考虑自旋的一般要修正吧！

谢谢，我看说明上给出的INCAR的例子写的ICHARG=1,而且要连续计算两次求差值，那第一次计算的时候ICHARG应该不用1而用0吧？

这里说的是DFT计算中可移植性的误差吧，统一体系的误差可能差不多，因此计算完之后计入这个修正就比较准确了
如果是气体的话要通过查热力学表找出相应的熵计入熵修正
具体怎么个意思我也不知道请高人指点

引用回帖:

Originally posted by sunyang1988 at 2010-08-09 20:01:58:
这里说的是DFT计算中可移植性的误差吧，统一体系的误差可能差不多，因此计算完之后计入这个修正就比较准确了
如果是气体的话要通过查热力学表找出相应的熵计入熵修正
具体怎么个意思我也不知道请高人指点

我是看到他给出的INCAR：
ISYM = 0 ! no symmetry
ISPIN = 2 ! allow for spin polarisation
VOSKOWN = 1 ! this is important, in particular for GGA
ISMEAR = 0 ! Gaussian smearing, otherwise negative occupancies
SIGMA = 0.1 ! intermid. smearing width
AMIX = 0.2 ! mixing set manually
BMIX = 0.0001
NELM = 20 ! 20 electronic steps
ICHARG = 1

这里面为什么要设置ICHARG=1呢？还有就是NELM=20会不会太小了呢？

你刚才提到的气体分子的熵修正是怎么回事呢？能否具体说一下呢，谢谢
，

引用回帖:

Originally posted by xiaojie7783 at 2010-08-10 00:31:19:


我是看到他给出的INCAR：
ISYM = 0 ! no symmetry
ISPIN = 2 ! allow for spin polarisation
VOSKOWN = 1 ! this is important, in particular for GGA
ISMEAR = 0 ! Gaussian smearing, otherwise negat ...

我认为这里是接着第一步检验赝势的好坏计算的，所以用ICHARG=1。
计算结合能时，直接优化并计算出体材的能量再减掉各原子的修正。