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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP pDOS相加应不应该等于总DOS
281/1返回列表
查看: 19717  |  回复: 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的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

weiql08

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

[ 发自手机版 http://muchong.com/3g ]
赞一下(1人) 回复此楼
6楼2016-11-25 21:10:45
已阅   回复此楼   关注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的回帖
普通回帖

obaica

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

青云丹鹤

至尊木虫 (文坛精英)
★
luzihen(金币+1): 谢谢参与
看不懂
回复此楼
4楼2016-11-25 21:09: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的回帖

obaica

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

beyond7765

新虫 (初入文坛)
★
luzihen(金币+1): 谢谢参与
大神们知道怎样把正交的paw basis 转为正交的原子轨道吗？需要用原子轨道描述的状态做材料性质的分析。

发自小木虫Android客户端
赞一下 回复此楼
21楼2016-11-26 10:05:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

obaica

禁言 (著名写手)
本帖内容被屏蔽
26楼2016-11-28 22:41:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

竹雨_朝颜

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
送红花一朵
顶！我在用QE计算时也碰到这个问题了！的确不相等！
赞一下 回复此楼
27楼2020-11-24 15:44:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Sahira

新虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
送红花一朵
确实是，刚刚画了一下不一样，很疑问，才过来查一下，感谢
赞一下 回复此楼
28楼2021-04-06 09:44:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
guogencai3楼
2016-11-25 21:00   回复  
luzihen(金币+1): 谢谢参与
xhmaohan5楼
2016-11-25 21:10   回复  
luzihen(金币+1): 谢谢参与
shenrenren7楼
2016-11-25 21:18   回复  
luzihen(金币+1): 谢谢参与
liubq088楼
2016-11-25 21:19   回复  
luzihen(金币+1): 谢谢参与
吸血人狼城9楼
2016-11-25 21:19   回复  
luzihen(金币+1): 谢谢参与
发自小木虫Android客户端
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): 谢谢参与
相关版块跳转 我要订阅楼主 luzihen 的主题更新
281/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +12 医学老男孩 2026-08-13 32/1600 2026-08-14 08:19 by 医学老男孩
[考博] 澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域 +71 AI8RaGaPaSCR 2026-08-07 74/3700 2026-08-14 08:15 by am8QyKH2Q7LX
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +9 Tide man 2026-08-10 10/500 2026-08-13 20:33 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 龙-樱
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +4 majunge000 2026-08-11 6/300 2026-08-13 16:50 by 四季常青藤
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[论文投稿] 职称评审，求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +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 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 关于代码变化问题，想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] fileCode有新解读？ +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
信息提示
关闭
请填处理意见
关闭