| 查看: 18943 | 回复: 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朵）
» 猜你喜欢
存款400万可以在学校里躺平吗
已经有36人回复
-
英文综述是否需要润色及查重
已经有5人回复
-
救命帖
已经有5人回复
-
限项规定
已经有5人回复
-
为什么nbs上溴 没有产物点出现呢
已经有9人回复
-
招博士
已经有3人回复
-
最失望的一年
已经有18人回复
-
求推荐英文EI期刊
已经有5人回复
» 抢金币啦！回帖就可以得到:
-
邀您投稿 Polymers 特刊
+1/487
-
东北大学-招收2026年硕士研究生2-3名（金属材料3D打印方向）
+5/125
-
保障电池制造超纯环境：威格科技电池生产全线解决方案
+1/86
-
华南师范大学胡勇军教授课题组招收2026年博士研究生
+1/84
-
招收博士生（大连理工大学，2026.09入学）
+1/65
-
中国地质大学（武汉）地质学、地质资源与地质工程、资源与环境方向招生，请尽快联系！
+1/54
-
西交利物浦大学黄彪院士招收26年全奖博士生1名（工业智能方向）
+1/31
-
2026湖南科技大学化学化工学院招新能源电池方向博士生
+1/30
-
西工大可塑性液体材料课题组招收软物质、流体力学、细胞工程等背景博士生
+1/28
-
论文投稿散金祈愿-电机学报
+2/26
-
博士后, 博士招生 美国大纽约地区 NJIT
+1/19
-
2026申博自荐 本硕双一流学科 纳米药物递送方向 一篇一区TOP 两个国家奖学金
+1/14
-
哈尔滨工业大学（深圳）-包斌招收2026年博士研究生及博士后
+1/9
-
华南理工大学宋波教授招收2026年博士生（二氧化碳有机、高分子和光电生物催化还原）
+1/8
-
【博士后/科研助理招聘-北京理工大学-集成电路与电子学院-国家杰青团队】
+1/5
-
华南理工大学宋波教授招聘材料和化学方向博士后（长期有限）
+1/4
-
0854电子信息330求调剂
+1/3
-
山东师范大学海外优青（校长团队）招聘硕士/博士/博后
+1/2
-
意大利米兰理工大学Lorenza Petrini教授招收玛丽居里MEDALLOY博士
+1/2
-
换工作
+1/2
6楼2016-11-25 21:10:45
25楼2016-11-28 20:41:01
★
luzihen(金币+1): 谢谢参与
luzihen(金币+1): 谢谢参与
|
本帖内容被屏蔽 |
2楼2016-11-25 20:57:14
4楼2016-11-25 21:09:45
10楼2016-11-25 21:28:48
送红花一朵 |
本帖内容被屏蔽 |
18楼2016-11-25 23:11:17
21楼2016-11-26 10:05:12
|
本帖内容被屏蔽 |
26楼2016-11-28 22:41:10
27楼2020-11-24 15:44:12
28楼2021-04-06 09:44:37
简单回复
2016-11-25 21:00
回复
luzihen(金币+1): 谢谢参与
xhmaohan5楼
2016-11-25 21:10
回复
luzihen(金币+1): 谢谢参与
2016-11-25 21:18
回复
luzihen(金币+1): 谢谢参与
liubq088楼
2016-11-25 21:19
回复
luzihen(金币+1): 谢谢参与
吸血人狼城9楼
2016-11-25 21:19
回复
caoquanrang11楼
2016-11-25 21:38
回复
luzihen(金币+1): 谢谢参与
lgp0112楼
2016-11-25 21:40
回复
luzihen(金币+1): 谢谢参与
顶。
dsctg13楼
2016-11-25 21:45
回复
luzihen(金币+1): 谢谢参与
psylhh14楼
2016-11-25 21:48
回复
luzihen(金币+1): 谢谢参与
tzynew15楼
2016-11-25 21:50
回复
luzihen(金币+1): 谢谢参与
tan91416楼
2016-11-25 22:04
回复
luzihen(金币+1): 谢谢参与
jixiangbao17楼
2016-11-25 22:40
回复
luzihen(金币+1): 谢谢参与
nono200919楼
2016-11-26 08:45
回复
luzihen(金币+1): 谢谢参与
·
gorgan20楼
2016-11-26 09:19
回复
luzihen(金币+1): 谢谢参与
祝福 [ 发自手机版 http://muchong.com/3g ]
gorgan22楼
2016-11-26 18:52
回复
祝福 [ 发自手机版 http://muchong.com/3g ]
youngen23楼
2016-11-27 12:47
回复
luzihen(金币+1): 谢谢参与
wtbccq50524楼
2016-11-27 12:48
回复
luzihen(金币+1): 谢谢参与