DFT-D2优化晶格常数和计算bader电荷分析，结果与文献不符
如题！！！
在晶格常数优化时所用的截断能，精度，范德华修正的参数及KPOINTS，POSCAR，POTCAR均和文献一致，可是优化得到的结果还是和文献值有一定差距，计算值比文献值偏小，计算得到的bader电荷比文献值偏小更多，这可能会是由什么原因造成的？
晶格常数实验值： a=b=3.674 c=5.033
晶格常数文献值： a=b=3.587 c=4.906
晶格常数计算值： a=b=3.549 c=4.886
bader电荷分析文献值为0.83，计算值为0.54。
附上晶格常数中的INCAR，其中红色线内均为文献给定值，还请大家帮忙看看是什么问题造成的，是参数设置问题吗？应该怎么解决，谢谢大家！！！
2.png
返回小木虫查看更多
今日热帖
京公网安备 11010802022153号
优化得到的结果还是和文献值有一定差距，计算值比文献值偏小，计算得到的bader电荷比文献值偏小更多
>>>>你说的POSCAR也一样，你确定么？是吸附还是取代体系看bader，杂质原子的位置也完全一样？
>>>>我觉得你的晶格常数部分无影响，1.07%的误差完全没有关系的。
>>>>bader可能就是跟你原子的位置是有关的，比如你的吸附了一个原子 是否吸附到合适位置
>>>>即便bader不一样，但是我们用bader是分析得失电子的，所以看的是delta电荷。你要有俩体系做比较，如果电荷变化不大那就没有关系了。【另外，看趋势骂，如果是得失电子规律是一样的，那么0.3e也没多大区别了】
谢谢你的回复！
>>>>我计算和文献中所用的实验值都是： a=b=3.674 c=5.033；
>>>>虽然晶格常数误差不是很大，但是bader计算太离谱，只好回过头来看是不是优化的问题了；
>>>>我计算的bader是体结构的，因为差太多就没继续计算表面和吸附体系的；
>>>>bader不止差0.3，理论上每个原子是1e，文献是0.83，我计算的只有0.54；
我尝试计算下表面的bader电荷，比较下试试!
看错了，bader是差0.3，理论是1，计算得到0.54是不是差太多，还是说只要后续计算的电荷转移趋势是正确的就OK。
>>>>我计算和文献中所用的实验值都是： a=b=3.674 c=5.033；
>>>>虽然晶格常数误差不是很大，但是bader计算太离谱，只好回过头来看是不是优化的问题了；
>>>>我计算的bader是体结构的，因为差太多就没继续计算表面和吸附体系的；
>>>>bader不止差0.3，理论上每个原子是1e，文献是0.83，我计算的只有0.54；
>>>>你的意思是，你现在在算的是理想块体的（没有任何缺陷参杂的）bader？那按理说结构不会有什么大错误。我也不好解释为什么差了0.3e了。 你说的理论上每个原子都是1e这个也只是理论上。我记得我以前也看过理想块体的bader，理论上是1的也是仅仅0.几。毕竟任何化学键都是既有共价键成分也有离子键成分。这点很好理解。
所以我一直觉得bader的绝对值意义不大，用来看电荷变化推测电荷转移就ok了
，
我觉得后续趋势正确就可以。你先全算完再比较看吧 可以后续再交流
好的好的 我先算完看看趋势 有问题再向你请教
-0.01的收敛限似乎不够，恐怕还不如默认的能量收敛限
POTIM有点大了
ISIF怎么回事？要变胞吗？通常优化晶格参数可以设一系列变值，然后拟合；也可以去ISIF=7直接优化，条件是不改变晶胞对称性。
K点也会影响到优化出来的晶格参数