| 查看: 19722 | 回复: 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» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有42人回复
-
filecode
已经有9人回复
-
请教兼职经验
已经有3人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有5人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有13人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有11人回复
-
读博的好处
已经有4人回复
-
filecode
已经有17人回复
-
FileCode能看出啥？
已经有32人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有15人回复
» 抢金币啦！回帖就可以得到:
-
北京信息科技大学机电工程学院硕士研究生招生线上宣讲会-8月16日
+5/280
-
东南大学杰青/海优团队招聘催化方向博士后
+1/89
-
胆酸锂 （Lithium cholate）
+1/87
-
威格科技：为化学合成打造全流程惰性气氛解决方案
+1/85
-
中国科学院空天院精密干涉测量团队 精密机械结构工程师招聘启事
+1/85
-
中国科学院空天院精密干涉测量团队 精密机械结构工程师招聘启事
+1/84
-
中国科学院空天院精密干涉测量团队招聘光学/算法/机械方向博士后/特别研究助理
+1/84
-
北京信息科技大学机电工程学院硕士研究生招生线上宣讲会——8月16日
+1/82
-
祈福！
+1/55
-
别研究filecode了，祝都有一个好结果：）
+2/50
-
物理学报录用，散币感谢合作者和审稿人
+2/34
-
韩国釜山大学储能电池与功能材料实验室诚招博后（电池/电化学/器件制备/MOF合成）
+1/13
-
浙江大学医学院附属第一医院招聘科研助理
+1/11
-
求七钼酸铵的cif文件
+1/7
-
海南大学化学院-科研助理招聘 （可读博，可与华中科技大学，天津大学联培）+博后招聘
+1/4
-
北理工集成电路杰青团队 | 诚招科助理
+1/4
-
中山大学李国强课题组招收博士生（柔性电子方向）
+1/3
-
北理工集成电路杰青团队 | 诚招科助理
+1/3
-
海南大学化学院—功能分子器件团队2027博士招生+博士后招聘
+1/3
-
国家级人才团队诚聘生物信息相关专业的科研助理
+1/3
28楼2021-04-06 09:44:37
★
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
回复
