求助DFT+U—U值的计算

2572094680

各位大佬，我在vasp上重复了官网上NiO的U值的计算，结果是吻合的，现在想算一个新的结构体系，在这个结构中有120个C原子，吸附了4个3d金属原子，我想算一下单个原子微扰和整体微扰这两种情况的U值来进行对比，现在的疑问是单个原子微扰就是只有1个Hubbard位点，整体微扰就是多个Hubbard位点吗。单个原子微扰和整体微扰具体要怎么操作啊，我是直接分别把MAGMOM=120*0  4*4和MAGMOM=120*0  4  3*4，POTCAR 中分别有C  Fe  Fe和C Fe，这样处理有道理吗，入门小白，可能表述的不清楚，谢谢大家了

卡开发发

扰动不是通过设置magmom来实现的，在LDAUTYPE=3的情况所设置的U对应扰动势的大小。多个Hubbard位点的情况需要在某个位点扰动的时候收集其他位点的占据数，这样可以排列成响应矩阵，将求倒数的过程改为矩阵求逆即可。

2572094680

卡开发发 发表于 2023-10-13 02:44
扰动不是通过设置magmom来实现的，在LDAUTYPE=3的情况所设置的U对应扰动势的大小。多个Hubbard位点的情况需 ...

非常感谢老师的回复，我还有以下几个问题：
1、如果我的结构是120个C,4个Fe，poscar写成了这样：
C    Fe    Fe  
120     1    3
然后非自洽计算的时候INCAR是：
LDAU=T
LDAUTYPE=3
LDAUL= -1 2 -1
LDAUU= 0 0.1 0
LDAUJ= 0 0.1 0
ICHARG=11
我通过改变LDAUU和LDAUJ的值来做响应计算，是这样吗
2、另外，我之前以为算整体微扰的时候是把加的4个金属原子看成一样的，就是poscar中是：
C    Fe     
120     4
非自洽计算的INCAR中是：
LDAU=T
LDAUTYPE=3
LDAUL= -1 2
LDAUU= 0 0.1
LDAUJ= 0 0.1
ICHARG=11
我这样是不对的是吗
3、您说的排列成矩阵，具体怎么排呢，是几行几列的矩阵呀
4、如果我只想计算一下在这个结构中3d金属的U值，以您的经验，单原子微扰和整体微扰，哪个结果更可取呀

卡开发发

2572094680 发表于 2023-10-13 11:43
非常感谢老师的回复，我还有以下几个问题：
1、如果我的结构是120个C,4个Fe，poscar写成了这样：
C    ...

1、暂时建议不要响应J，那个比较复杂。
2、不行，你可能需要对POTCAR的元素名称做修改使它和其他同元素原子标记不同，虽然这种操作非常繁琐。
3、只考虑U的响应，原则上应当是四行四列。
4、原理上说应当是进行整体响应较为合理，实际情况是，（1）复杂环境下响应过程太过于复杂（2）即便响应也可能得不到合理结果。