ShangChien 发表于 2021-12-31 16:21
谢谢回答，现在GPAW.21.6.0直接用pip安装就可以基本使用，所以安装还是很方便的，效果如下：

即使想要 ...

请问你还在用gpaw吗？
gpaw的例子你都在哪找的？除了H2那个

FD和LCAO的处理都有文献，我也就简单说一下：
1、有限差分就是通过有限差分的方法构造动能和密度的梯度这些量，这样也可以将这些量表示成矩阵形式。例如d^2f/dx^2~f(x[i+1])-2f(x)+f(x[i-1])/2h，其中h是x的网格间距，这样d^2f/dx^2就可以写成类似于N对角矩阵，积分也是按照网格积分。
3、LCAO的方法其实应该用的是赝原子轨道，构造方法类似于SIESTA，是自动产生的，但极化基组采用的方案与SIESTA不同。LCAO功能比较不全，比如没办法变晶胞优化。LCAO的准确程度我没找到具体的数据，但PW和Grid的delta测试（https://molmod.ugent.be/deltacodesdft）表现应该还马马虎虎。
4、Eggbox误差，这个是引入网格后因网格不完备导致的平移对称性的破缺。之所以叫Eggbox就是因为好比你的体系是鸡蛋，鸡蛋在放在装鸡蛋的槽子当中，如果槽的尺寸很大，你把鸡蛋稍微调整一点位置，这个鸡蛋的状态和原来是不一样的；但是当槽的尺寸细小，这种差异就没有了。

GPAW我用过一阵子，感觉就是依赖多导致编译复杂，另外就是不同模式下到底哪些功能不好使可能要折腾很久，我觉得能把一种模式做完善就挺不错。