| 查看: 19725 | 回复: 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后的代码验证其规律
已经有46人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有4人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有16人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有11人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有4人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有8人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有37人回复
-
好奇怪的filecode
已经有7人回复
-
filecode
已经有9人回复
-
请教兼职经验
已经有3人回复
» 抢金币啦！回帖就可以得到:
-
2026年财务管理、经济与信息技术国际会议 (FMEIT 2026)
+1/903
-
上海男征女 条件优秀可加微信联系
+1/437
-
真诚找女友
+1/164
-
胆酸锂 （Lithium cholate）
+1/87
-
东南大学杰青/海优团队招聘催化方向博士后
+1/86
-
中国科学院空天院精密干涉测量团队 精密机械结构工程师招聘启事
+1/84
-
国防科技大学空天科学学院招聘启事陶瓷基复合材料方向科研助理招聘
+1/83
-
广西大学 化学化工学院 新入职教师招生 化学、化学工程与技术、材料化工专业
+1/54
-
【坐标北京，诚征男友，结婚对象】
+1/50
-
请问谁有2025最新影响因子Excel表？金币求助！
+1/27
-
希望这篇论文能中
+2/24
-
坐标郑州，征女友
+2/16
-
西北工业大学 生科院 秦绪军课题组诚招博士后、博士/硕士研究生【长期有效】
+1/14
-
论文祈福
+1/11
-
关于Ghaffar-Parkins催化剂催化氰基水解为酰胺的问题
+1/6
-
北理工集成电路杰青团队 | 诚招科助理
+1/4
-
论文万能框架
+1/4
-
中山大学计算数学方向招2027届考核制博士生
+1/3
-
国家级人才团队诚聘生物信息相关专业的科研助理
+1/3
-
国家纳米科学中心王振兴课题组招收微/纳流控方向博士后研究人员
+1/1
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
回复
