卡开发发 发表于 2022-2-8 11:56
巡游与否看最后算出来局域性如何，强行去让电荷局域在某个原子附近是很主观的做法。

明白了，感谢老师指导！

Clifford 发表于 2022-2-8 11:33
感谢卡开发发老师！我是在如下文献中看到这个做法，作者似乎是通过调整体系总电子数实现对不同电荷态空位 ...

巡游与否看最后算出来局域性如何，强行去让电荷局域在某个原子附近是很主观的做法。

在qe中的操作应该仅需要设置体系总（per cell）电荷。对于bulk体系，qe自身没办法进行的修正，因此有两种方式可以进行尝试：
（1）超胞法。体系介电常数稍大一些，库伦作用屏蔽比较强，带电空位的过剩电荷比较局域化，此时足够的超晶格就可以使得各项随超晶格尺寸增加收敛及校正项足够小，或者也可以适当外推超胞尺寸来逼近Ncell->∞的情况。
（2）修正项。如果上述外推不容易实现的情况，可以试试文献中提到CoFFEE这个程序，不过我没有使用过这个程序不确定会存在什么坑，例如版本匹配或其他问题。如果存在问题你可以尝试查找是否有其他程序可以使用，如果实在不行我们也可以考虑怎么修改CoFFEE这个程序。