高斯计算铁的PBC问题请教

ovovovo

根据大佬对于pbc计算博文http://bbs.keinsci.com/forum.php ... d=146&highlight=PBC学到了以下输入文件，文件正常结束，但是算出来的结果和同样泛函基组下单独Fe原子计算E值差了一百kcal/mol，差距太大感到疑惑，是我输入文件没有学好嘛？

zjxitcc

（1）算不了，不要用小刀砍大树。周期性计算就要用适合 周期性计算的程序，例如CP2K/QE/VASP/ABACUS。注意，我没有让你说只能学习VASP，不要误解了。学习道路一旦偏了，只能在毕不了业延期的时候（侥幸可能）算出你想要的东西。
（2）孤立Fe原子 对应的是 一个很大的单胞中间放一个Fe原子，不是你这种铁单胞，您结构搞错了。
（3）关键词严重不合理。iop(5/13=1,5/33=1)都跑出来了。

Uus/pMeC6H4-/キ

别乱用关键词，参考
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331（http://bbs.keinsci.com/thread-3460-1-1.html）
如果刚开始做计算，仔细阅读
谈谈学量子化学如何正确地入门
http://sobereva.com/355（http://bbs.keinsci.com/thread-4447-1-1.html）

你链接的帖子12楼的更新的回复都指出了做正经的周期性体系计算该用什么、如何学习，我实在不理解为何不好好读完帖子……

Lance先生

我看到了不得了的东西。“iop(5/13)=1”，这可不兴学啊。

sobereva

差值对应内聚能
北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）里专门讲了。如果你要长期做正经的第一性原理计算，十分建议通过此培训系统性学学。另外，Gaussian的PBC功能属于玩具级的，能表现出实用性的场合很少。

ovovovo

后续是用高斯对Fen(n=1-16)都进行了计算，最后取一个相对能量最低的，审稿人没有再说这个问题，不过一般高斯计算不会遇到这个问题。

Uus/pMeC6H4-/キ

ovovovo 发表于 2025-5-10 12:55
后续是用高斯对Fen(n=1-16)都进行了计算，最后取一个相对能量最低的，审稿人没有再说这个问题，不过一般高 ...

这意思是把铁晶胞换成一堆小型孤立铁原子团簇，再除以原子数来计算所谓“单原子铁”的能量？听起来还是不靠谱，孤立过渡金属团簇这种强相关体系的能量计算绝对没这么简单。自旋多重度不唯一，每个自旋多重度下稳定波函数解也可能不唯一，都不知道怎么考虑的。而且小团簇的表面效应太明显了，和晶体的体相环境差异巨大。

审稿人不质疑不代表结果真的合理，说难听点，这次被放过去不能保证以后还可以蒙混过关，用错误的方法发文章等于同时给自己和后人挖坑。