用纯泛函计算固体能带结构的可行性？

KSeGaSn · 发表于 Post on 2023-9-14 23:33:31

各位老师大家好，我初学第一性原理的计算不久，目前想要计算一个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来计算能带结构是可以接受的吗？

sobereva · 发表于 Post on 2023-9-15 00:57:28

可以接受

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

KSeGaSn · 发表于 Post on 2023-9-15 15:57:31

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

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

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

sobereva · 发表于 Post on 2023-9-16 01:14:52

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是迫不得已的。

KSeGaSn · 发表于 Post on 2023-9-16 09:32:45

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

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[综合讨论] 用纯泛函计算固体能带结构的可行性？