求助 周期性体系吸收谱计算

wang6312

我有一个周期性的铁基材料（FeOCl），一个晶胞比较大（82个原子），实验测得UV-vis在~13500 cm-1有一个吸收，但是用PBE泛函算100晶面的吸收谱却没有这个峰（因为我是做一个反应，推测是在100面上发生，就算这个表面的吸收谱了）。因为我们这个思路是参考Nature的一篇文献（doi:10.1038/nature19059，他们是算沸石中Fe骨架的吸收谱，整个体系13个原子，用的CASPT2算的），所以考虑可能是方法的问题导致计算不准（因为我们也用PBE重复了Nature那个体系，也没有算出他们实验中得到的峰。）。所以我想问下需要怎么调整计算方法才能算出和实验吻合的吸收峰呢？我们初步设想：（1）在晶胞中取一个小的Cluster代表我们的Fe site，10个原子左右，这样就可以用CASPT2算了；（2）保持我们的Slab模型不变，用一个比PBE精度高，但是计算量比CASPT2小的方法算（比如TDDFT？可以吗？）不知是否可行？

喵星大佬

Fe这种高自旋且很难确定自旋的体系不要用团簇来算，用CASTEP/QE/VASP之类的PW程序算周期性结构
虽然结果也不怎么可靠就是了

卡开发发

Slab算不出来可能是因为PBE低估Gap，吸收谱结果不会好看。可以考虑更好的泛函，但体系稍大一点代价还是很大。

wang6312

卡开发发 发表于 2020-8-28 10:32
Slab算不出来可能是因为PBE低估Gap，吸收谱结果不会好看。可以考虑更好的泛函，但体系稍大一点代价还是很大 ...

谢谢您的解答，那我们如果尝试用精度更高的泛函去计算的话，像我们这么大的体系估计要算多久呢？能不能粗略估计下？（几周？还是要几个月？）因为我们是主要做实验的，别人配合我们在做这个计算，我们自己对这个不太了解

biogon

纯泛函的tddft算出来的相当不靠谱
起码要用PBE0

卡开发发

wang6312 发表于 2020-8-28 14:21
谢谢您的解答，那我们如果尝试用精度更高的泛函去计算的话，像我们这么大的体系估计要算多久呢？能不能粗 ...

PW程序如果不是服务器多得不得了，至少是按周起算，甚至还可能算不动。团簇做不是不可以，一个是边界要处理，一个是太小的团簇没法复现真实的表面，自旋问题也不至于不能解决，只是繁琐些。