计算化学公社

标题: 求助：CP2K中的DFT+U [打印本页]

作者
Author: zhukangg 时间: 2020-12-20 09:24
标题: 求助：CP2K中的DFT+U
最近在学习CP2K，计算的时候需要使用DFT+U的方法，按照之前在VASP中的参数设置了U值，使用的是MULLIKEN方法，但是显示了warning：
*** WARNING in dft_plus_u.F:1457 :: DFT+U energy contibution is negative ***
*** possibly due to unphysical Mulliken charges. Check your input, if ***
*** this warning persists or try a different method! ***

请问在CP2K中，U值该如何设置，与VASP中保持一致吗？

作者
Author: 卡开发发 时间: 2020-12-20 10:03

*** WARNING in dft_plus_u.F:1457 :: DFT+U energy contibution is negative ***
*** possibly due to unphysical Mulliken charges. Check your input, if ***
*** this warning persists or try a different method! ***

看提示来说，试试换成Lowdin的。原则上不同的基组和方法、投影态等U值可能会有些许差异。

作者
Author: zhukangg 时间: 2020-12-20 10:18

卡开发发发表于 2020-12-20 10:03
看提示来说，试试换成Lowdin的。原则上不同的基组和方法、投影态等U值可能会有些许差异。

好的，谢谢。我把method换成了Muliken chaeges，就没有这个问题了

作者
Author: Penson 时间: 2020-12-22 10:15
不同软件的U值不能等同使用，差异较大，需要参考相应软件计算的文献或者自己测试

作者
Author: 丁越 时间: 2021-7-14 11:01

Penson 发表于 2020-12-22 10:15
不同软件的U值不能等同使用，差异较大，需要参考相应软件计算的文献或者自己测试

那请教一下老师，如果用cp2k做DFT+U，这个U值该咋确定啊？（如果不可以用VASP的U值的话）

作者
Author: Penson 时间: 2021-7-15 09:20

丁越发表于 2021-7-14 11:01
那请教一下老师，如果用cp2k做DFT+U，这个U值该咋确定啊？（如果不可以用VASP的U值的话）

参考文献或自己测试不同U值结果的合理性

作者
Author: 丁越 时间: 2021-7-16 09:24

Penson 发表于 2021-7-15 09:20
参考文献或自己测试不同U值结果的合理性

老师我绘制了官网教程例子里的WO3做PDOS，在Materials project 中给的W的U值是6.2ev(VASP),那么这个U是加给W的d呢还是f呢？我把这些图都绘制了一遍，发现只要晶胞大小合理，DOS图基本是正确的，加U值有些改进甚微。

作者
Author: Penson 时间: 2021-7-16 10:24

丁越发表于 2021-7-16 09:24
老师我绘制了官网教程例子里的WO3做PDOS，在Materials project 中给的W的U值是6.2ev(VASP),那么这个U是加 ...

没看懂，到底是MS还是VASP？加在哪个轨道这个得看看文献吧，一般只有到稀土元素会考虑到f轨道的影响，之前的元素主要还是d轨道的贡献更大

作者
Author: 卡开发发 时间: 2021-7-16 13:26

丁越发表于 2021-7-16 09:24
老师我绘制了官网教程例子里的WO3做PDOS，在Materials project 中给的W的U值是6.2ev(VASP),那么这个U是加 ...

MP的结构的页面下面有这个结构获取的数据采用的输入文件，去看他给的INCAR加在哪就行了。

作者
Author: 丁越 时间: 2021-7-16 18:37

Penson 发表于 2021-7-16 10:24
没看懂，到底是MS还是VASP？加在哪个轨道这个得看看文献吧，一般只有到稀土元素会考虑到f轨道的影响，之 ...

不好意思老师，我忘说了我用的是cp2k。这个U值来自于MP（Materials project）网站，是用VASP算的。我把U加到d和f都试了一遍，发现加到f就没变化。

作者
Author: 丁越 时间: 2021-7-16 18:42

卡开发发发表于 2021-7-16 13:26
MP的结构的页面下面有这个结构获取的数据采用的输入文件，去看他给的INCAR加在哪就行了。

谢谢卡开老师，我正在学习cp2k。论坛里您关于DFT+U的解答我都看了一遍，我想请教一下对于cp2k来说，这个U该怎么获取呢（如果不能参考VASP之类算的U的话）？还有在DOS，能带，弹性模量这些性质计算的时候我直接采取杂化泛函去算，不用U去矫正，这样是否可以呢？

作者
Author: 卡开发发 时间: 2021-7-16 22:26

丁越发表于 2021-7-16 18:42
谢谢卡开老师，我正在学习cp2k。论坛里您关于DFT+U的解答我都看了一遍，我想请教一下对于cp2k来说，这个U ...

我了解的情况似乎没什么太好的办法，除非有人愿意贡献这部分的代码。DOS能带可能杂化泛函还好，但弹性模量会比较昂贵。

作者
Author: 丁越 时间: 2021-7-17 10:06

卡开发发发表于 2021-7-16 22:26
我了解的情况似乎没什么太好的办法，除非有人愿意贡献这部分的代码。DOS能带可能杂化泛函还好，但弹性模 ...

好的老师，我知道了，谢谢啦

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