| 查看: 18926 | 回复: 27 | |||||||
| 本帖产生 1 个 1ST强帖 ，点击这里进行查看 | |||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||||
[交流]
VASP pDOS相加应不应该等于总DOS
|
|||||||
|
今天有人讨论vasp算的pDOS加起来为什么不等于总DOS,做了点测试查了点文献,记录一下想法.如果有错请一定指正我. 首先给结论,不应该相等. 先看DOS定义[1]  这里用的eigenstate是Kohn-Sham eigenstate.是正交归一的. 再看pDOS定义[1]  这里|i><i|是投影算符.如果有|i>是正交归一的,那么有  就有  即各个轨道pDOS的和等于总DOS. 所以看pDOS到底等不等于总DOS就要看投影pDOS的轨道是不是正交归一的. vasp的说明书[2]并没有明确给出vasp的计算的算法,但是看计算结果,一般能得到的是投影在s,p,d,f等轨道的pDOS,也就是atomic orbital的pDOS.然而由于原子轨道并非正交归一,所得的和并不一定保证是总DOS. vasp说明书里面还提到了Wagner Seitz radius,也没有给出具体算法,GAPW网站给出了一个Wigner-Seitz LDOS的定义,如下  可见是画一个球,球里面的原子轨道pDOS都算,之外的都不算.基础依然是原子轨道,所以加起来依然不能保证等于总DOS.这么一看， 球不能把整个空间填满,所以也不应该相等.另外一篇文章给了一段话也说了.[3]  随手画了一个体系的总DOS和pDOS和，确实不相等.  上面说的有问题请指正，有兴趣的拿自己的体系测试一下看看相等不。 [1]https://wiki.fysik.dtu.dk/gpaw/documentation/pdos/pdos.html#density-of-states [2]https://cms.mpi.univie.ac.at/vasp/vasp/LORBIT.html [3]Maintz, Stefan, et al. "Analytic projection from plane‐\wave and PAW wavefunctions and application to chemical‐\bonding analysis in solids." Journal of computational chemistry 34.29 (2013): 2557-2567. |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理 | 第一性原理计算-tg | 表界面计算 | vasp计算方法收集 |
» 本帖已获得的红花（最新10朵）
weiql08» 猜你喜欢
论文终于录用啦！满足毕业条件了
已经有16人回复
-
求个博导看看
已经有19人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
» 抢金币啦！回帖就可以得到:
-
上海交通大学医学院王戈林课题组招聘博士后和助理研究员
+1/272
-
限广州，征女友
+2/96
-
江苏科技大学能源材料化学课题组张俊豪教授招收博士研究生1-2名
+1/71
-
山东征女友，坐标济南
+1/63
-
2026博士申请——有机化学\计算化学\药物化学方向
+1/44
-
中国地质大学（北京）王琳课题组招收2026年硕转博/申请-考核博士研究生-1月8日截止
+1/34
-
北京林业大学木质素高值化利用创新团队招收2026年入学博士生
+1/31
-
南科大薛亚辉课题组诚聘离子输运、低维器件、原子力显微镜等方向“快响行动”博士生
+1/28
-
【AI、水文方向】香港科技大学（广州）研究助理招聘
+1/26
-
福建师范大学柔性电子学院招收2026年博士（储能材料与柔性电子器件）
+2/26
-
SCI，计算机相关可以写
+1/24
-
王志博教授课题组招收硕士研究生(本招收信息长期有效)
+2/20
-
大叔征婚
+1/11
-
吉林大学材料物理本科生求问调剂信息
+1/11
-
【博士后/科研助理招聘-北京理工大学-集成电路与电子学院-国家杰青团队】
+1/7
-
乙酸乙酯如何除丙烯酸甲酯
+1/6
-
中国科学院苏州纳米所院士团队博士后岗位招聘
+1/4
-
废旧塑料热解油采购
+1/4
-
中国矿业大学黄赳课题组联合中国科学院南京土壤研究所朱晓芳研究员诚聘博士后
+1/2
-
上海理工大学 生物医学工程专业 招收2026年全日制博士生一名
+1/1
21楼2016-11-26 10:05:12
★
luzihen(金币+1): 谢谢参与
luzihen(金币+1): 谢谢参与
|
本帖内容被屏蔽 |
2楼2016-11-25 20:57:14
6楼2016-11-25 21:10:45
10楼2016-11-25 21:28:48
简单回复
吸血人狼城9楼
2016-11-25 21:19
回复
