应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP pDOS相加应不应该等于总DOS
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 19103  |  回复: 27
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 1ST强帖 ，点击这里进行查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

luzihen

金虫 (正式写手)

[交流] VASP pDOS相加应不应该等于总DOS

今天有人讨论vasp算的pDOS加起来为什么不等于总DOS,做了点测试查了点文献,记录一下想法.如果有错请一定指正我.

首先给结论,不应该相等.

先看DOS定义[1]
VASP pDOS相加应不应该等于总DOS
这里用的eigenstate是Kohn-Sham eigenstate.是正交归一的.

再看pDOS定义[1]
VASP pDOS相加应不应该等于总DOS-1
这里|i><i|是投影算符.如果有|i>是正交归一的,那么有
VASP pDOS相加应不应该等于总DOS-2
就有
VASP pDOS相加应不应该等于总DOS-3
即各个轨道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的定义,如下
VASP pDOS相加应不应该等于总DOS-4
可见是画一个球,球里面的原子轨道pDOS都算,之外的都不算.基础依然是原子轨道,所以加起来依然不能保证等于总DOS.这么一看，
球不能把整个空间填满,所以也不应该相等.另外一篇文章给了一段话也说了.[3]
VASP pDOS相加应不应该等于总DOS-5

随手画了一个体系的总DOS和pDOS和，确实不相等.
VASP pDOS相加应不应该等于总DOS-6

上面说的有问题请指正，有兴趣的拿自己的体系测试一下看看相等不。


[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 obaica 竹雨_朝颜 Sahira

» 猜你喜欢

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2016-11-25 20:55:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

obaica

禁言 (著名写手)
送红花一朵
本帖内容被屏蔽
18楼2016-11-25 23:11:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 28 个回答

obaica

禁言 (著名写手)
★
luzihen(金币+1): 谢谢参与
本帖内容被屏蔽
2楼2016-11-25 20:57:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weiql08

银虫 (小有名气)
★
luzihen(金币+1): 谢谢参与
送红花一朵
LDOS跟PDOS还不是一回事，昨天看到一个提法，对于金属，提到了赝带隙，用LDOS分析。不知对不对？

[ 发自手机版 http://muchong.com/3g ]
赞一下(1人) 回复此楼
6楼2016-11-25 21:10:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luzihen

金虫 (正式写手)
引用回帖:
6楼: Originally posted by weiql08 at 2016-11-25 21:10:45
LDOS跟PDOS还不是一回事，昨天看到一个提法，对于金属，提到了赝带隙，用LDOS分析。不知对不对？

LDOS是不一样的,但是看vasp里面说的意思好像是这么做的..
赞一下 回复此楼
10楼2016-11-25 21:28:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
吸血人狼城9楼
2016-11-25 21:19   回复  
luzihen(金币+1): 谢谢参与
发自小木虫Android客户端
查看全部 28 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 08工科求调剂286 +4 tgs_001 2026-03-28 4/200 2026-03-30 21:03 by dophin1985
[考研] 277跪求调剂 +8 1915668 2026-03-27 12/600 2026-03-30 21:01 by dophin1985
[考研] 一志愿北化085600材料专硕275|有文章专利｜求调剂 +8 Micky11223 2026-03-25 8/400 2026-03-30 20:55 by dophin1985
[考研] 297求调剂 +3 田洪有 2026-03-30 3/150 2026-03-30 20:43 by wxiongid
[考研] 各位老师好，我的一志愿为北京科技大学085601材料专硕 +10 Koxui 2026-03-28 10/500 2026-03-30 19:33 by 源_2020
[考研] 085602 化学工程专硕 340分求调剂 +4 qianbai11 2026-03-29 4/200 2026-03-30 11:34 by 唐沐儿
[考研] 318求调剂 +7 陈晨79 2026-03-30 7/350 2026-03-30 10:49 by 探123
[考研] 337求调剂 +6 《树》 2026-03-29 6/300 2026-03-30 10:15 by herarysara
[考研] 349求调剂 +6 李木子啊哈哈 2026-03-25 6/300 2026-03-29 12:47 by 无际的草原
[考研] 279求调剂 +4 蝶舞轻绕 2026-03-29 4/200 2026-03-29 09:45 by laoshidan
[考研] 一志愿厦门大学化学学硕307求调剂 +10 y7czhao 2026-03-26 10/500 2026-03-28 14:23 by 唐沐儿
[考研] 304求调剂 +6 曼殊2266 2026-03-27 6/300 2026-03-28 14:10 by 唐沐儿
[考研] 286求调剂 +12 PolarBear11 2026-03-26 12/600 2026-03-28 12:14 by zllcz
[考研] 265求调剂11408 +3 刘小鹿lu 2026-03-27 3/150 2026-03-27 20:53 by nihaoar
[考研] 化学调剂 +4 爱吃番茄的旭 2026-03-24 5/250 2026-03-27 17:50 by kiokin
[考研] 一志愿 西北大学 总分282 英语一62 求调剂 +7 18419759900 2026-03-25 8/400 2026-03-27 16:38 by 18419759900
[考研] 08开头275求调剂 +4 拉谁不重要 2026-03-26 4/200 2026-03-27 14:12 by Delta2012
[考研] 315调剂 +4 0860求调剂 2026-03-26 5/250 2026-03-27 11:23 by wangjy2002
[考研] 314求调剂 +3 溪云珂 2026-03-26 3/150 2026-03-27 11:20 by sanrepian
[考研] 材料专硕331求调剂 +4 鲜当牛 2026-03-24 4/200 2026-03-24 15:58 by JourneyLucky
信息提示
关闭
请填处理意见
关闭