VASP+HSE06计算配体场d轨道分裂和pdos分析问题求助

上清

请教各位老师同学：
通过HSE计算得到的DOS图中轨道能级与理论不符，我算了一下钙钛矿的PDOS d轨道的五个分轨道的能级（和d带中心一样的积分算法）结果是dxy = dyz = dxz > dz2 = dx2-y2
但是八面体场不应该是 dz2 = dx2-y2 >dxy = dyz = dxz 这样的顺序吗？（图中只绘制了Ti原子的DOS，没有显示的线条是被覆盖了）



POSCAR如下：

Ca1 Ti1 O3
1.0
   3.8894709999999990    0.0000000000000000    0.0000000000000002
   0.0000000000000006    3.8894709999999990    0.0000000000000002
   0.0000000000000000    0.0000000000000000    3.8894709999999990
Ca Ti O
1 1 3
direct
   0.5000000000000000    0.5000000000000000    0.5000000000000000 Ca2+
   0.0000000000000000    0.0000000000000000    0.0000000000000000 Ti4+
   0.0000000000000000    0.5000000000000000    0.0000000000000000 O2-
   0.0000000000000000    0.0000000000000000    0.5000000000000000 O2-
   0.5000000000000000    0.0000000000000000    0.0000000000000000 O2-

上清

elpa 发表于 2025-2-26 10:06
请教一下怎么做的积分？用的什么脚本

pymatgen可以积分。我是自己写的公式所以搞错了。就是x轴指定范围内的积分就行

elpa

请教一下怎么做的积分？用的什么脚本

上清

北大-陶豫 发表于 2025-1-17 16:32
0 - 10 eV 这段确实是 dz2 = dx2-y2 >dxy = dyz = dxz ，有什么问题吗

谢谢老师回复，之前的问题已经解决了。是我之前python计算d带中心的公式写错了

北大-陶豫

0 - 10 eV 这段确实是 dz2 = dx2-y2 >dxy = dyz = dxz ，有什么问题吗