计算化学公社

标题: CP2K使用DFT+U方法时应该如何确定U [打印本页]

作者
Author:
qma    时间: 2023-10-4 20:28
标题: CP2K使用DFT+U方法时应该如何确定U
电极材料常常涉及到各种过渡金属氧化物,为了准确计算因此要用到DFT+U。
但CP2K的参数U不像VASP可以直接用现成的,想请教各位老师在CP2K使用DFT+U方法时应该如何确定U?


作者
Author:
sobereva    时间: 2023-10-5 06:50
从文献里找适合当前体系当前研究的问题的。很多文献都给了,诸如J. Chem. Theory Comput. 2019, 15, 5987−5997提出适用于CP2K的LnPP1镧系基组和赝势时直接给了+U参数(实为Ueff)。注意CP2K用的+U和VASP、QE等平面波程序的大小明显不同,不可套用,一些CP2K用户的文章里的+U参数实际上是无脑套用其它文献里平面波程序的,因此并不合适,必须注意文中+U参数的来源。

也可以自己根据实验数据或者HSE06等杂化泛函算的质量较好的数据拟合Ueff。

CP2K做杂化泛函计算非常快,一千原子在像样的双路服务器上杂化泛函开ADMM算单点都能跑得动,除非真的算不动再说+U的事。如果是直接拿很小晶胞来算而需要考虑k点,CP2K的杂化泛函和+U都不支持k点,这种情况适合用QE等程序。





欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3