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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP pDOS相加应不应该等于总DOS
51/1返回列表
查看: 19711  |  回复: 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的回帖

luzihen

金虫 (正式写手)
引用回帖:
21楼: Originally posted by beyond7765 at 2016-11-26 10:05:12
大神们知道怎样把正交的paw basis 转为正交的原子轨道吗？需要用原子轨道描述的状态做材料性质的分析。

不会.你搜一下lobster...可能可以...
赞一下(1人) 回复此楼
25楼2016-11-28 20:41:01
已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域 +67 AI8RaGaPaSCR 2026-08-07 70/3500 2026-08-14 07:15 by XFfZ89kVC80Q
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +12 医学老男孩 2026-08-13 30/1500 2026-08-14 06:02 by zyjwdhcs
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +3 启萌科技 2026-08-12 11/550 2026-08-13 22:57 by 启萌科技
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] filecode +14 documentary 2026-08-10 16/800 2026-08-13 16:39 by 润生71
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读？ +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
信息提示
关闭
请填处理意见
关闭