求助，关于 DFT+U 和 PDOS 概念的疑问

SeriTaku

刚开始念研究生，对概念不清晰，如果有什么荒谬的见解，还请大家原谅。
在计算一个关于 Tb 元素的体系，在不考虑 DFT+U 的情况下，能带位置和 ARPES 实验有一定差距。
考虑 DFT+U 的时候出现了一些概念上的问题，Tb 的价电子里不包含 d 电子，所以一开始 DFT+U 选择了考虑 f 轨道，即 VASP 参数中选择 LDAUL=3，计算后能带形状完全没有变动。
看了一下 PDOS, 发现与实验不相符的区间， Tb 的 f 轨道几乎没有贡献，但是 d 轨道贡献很大，于是 DFT+U 选择了 d 轨道，即 VASP 参数中 LDAUL=2，凑出了和实验上比较接近的图样。
这诞生了两个疑问。
1. 我选择的赝势是 PAW_PBE Tb_3 06Sep2000 这个赝势应该是不包含 d 电子的，为什么 PDOS 中却能看到 Tb d 轨道贡献很大？这还包含了我的一个疑问，PDOS 的投影是对某个原子的某个轨道的投影，那它能说明什么事情？说明这个轨道对晶体某段区间的能量的贡献程度？
2. 这种情况下 DFT+U 选择 LDAUL=2 作为参数有合理性吗？还是说是一个完全错误的选择，做出与实验相近的结果只是凑巧而已？
还望大家不吝赐教。

SeriTaku

granvia 发表于 2022-8-19 11:33
PDOS是借助第三方软件画的吗？有没肉眼检查一下数据，排除一下软件bug？

看了一下 PROCAR 和 OUTCAR，发现关于 Tb 是没有任何 f 电子的数据的。于是把 Tb 用 Tb_3 这个赝势单独拿出来算了一下，发现 OUTCAR 的total charge 里也没有 f 电子。这让我感觉很疑惑。

granvia

PDOS是借助第三方软件画的吗？有没肉眼检查一下数据，排除一下软件bug？

SeriTaku

乐平 发表于 2022-8-19 11:04
猜测除了 Tb 元素外，体系还有过渡金属？ d 轨道的贡献是由过渡金属产生的吧？

是的，体系内有过渡金属元素，并且调整它的 U 值无法调整与实验不符的那部分能带了，画出 fatband 会发现这些能带主要由 Tb 的 d 轨道贡献，这是很困扰我的情况，以及这种情况我是否应该对 Tb 的 d 轨道加 U 计算呢？

乐平

猜测除了 Tb 元素外，体系还有过渡金属？ d 轨道的贡献是由过渡金属产生的吧？