计算化学公社

标题: 用纯泛函计算固体能带结构的可行性？ [打印本页]

作者
Author: KSeGaSn 时间: 2023-9-14 23:33
标题: 用纯泛函计算固体能带结构的可行性？
各位老师大家好，我初学第一性原理的计算不久，目前想要计算一个2D COF的能带结构，晶格常数约25 x 25 x 3.5，仅在z轴方向上扩胞3倍，超胞里就有约400个原子了，因此即使在cp2k先用PBE得到了初猜再切换成HSE06也算不动。我了解到纯泛函HLE17（https://pubs.acs.org/doi/10.1021/acs.jpcc.7b01066）好像对于带隙的计算是比较好的，目前计算的体系和实测虽然有差（0.4eV）但显著好于PBE（估小了1eV），请问此时使用HLE17来计算能带结构是可以接受的吗？

作者
Author: sobereva 时间: 2023-9-15 00:57
可以接受

CP2K本来也不支持杂化泛函考虑k点，提HSE06也没意义

作者
Author: KSeGaSn 时间: 2023-9-15 15:57

sobereva 发表于 2023-9-15 00:57
可以接受

CP2K本来也不支持杂化泛函考虑k点，提HSE06也没意义

好的，谢谢sob老师。另外想问一下老师，如果想要以TDDFT继续研究固体的激发性质的话，用PBE或者HLE17会有定性的区别吗（还是因为扩成超胞用杂化泛函算不太动，有了能带结构之后对能量准确度容忍度应该较好）？

作者
Author: sobereva 时间: 2023-9-16 01:14

KSeGaSn 发表于 2023-9-15 15:57
好的，谢谢sob老师。另外想问一下老师，如果想要以TDDFT继续研究固体的激发性质的话，用PBE或者HLE17会有 ...

HLE17结合TDDFT用的情况目前还很少见

TD-杂化泛函400原子未必算不动。之前我用2*7R32 512GB机子用81核并行算250多个原子周期性体系，TDDFT算最低50个态，PBE0/DZVP-MOLOPT-SR-GTH开ADMM，也就花了三个小时而已。实在没辙再用PBE，这样才有充分理由在审稿人问及时说用PBE是迫不得已的。

作者
Author: KSeGaSn 时间: 2023-9-16 09:32

sobereva 发表于 2023-9-16 01:14
HLE17结合TDDFT用的情况目前还很少见

TD-杂化泛函400原子未必算不动。之前我用2*7R32 512GB机子用81核 ...

好的，谢谢老师！我也试一下如此计算。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)