版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (6452)
>导师招生 (2520)
>考研 (2385)
>虫友互识 (516)
>文献求助 (453)
>考博 (348)
>硕博家园 (265)
>招聘信息布告栏 (133)
>论文投稿 (116)
>休闲灌水 (116)
>博后之家 (88)
>基金申请 (84)
>公派出国 (50)
>绿色求助(高悬赏) (43)
>找工作 (42)
>教师之家 (41)

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

本帖产生 1 个 1ST强帖，点击这里进行查看

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

luzihen

金虫 (正式写手)

1ST强帖: 1
应助: 47 (小学生)
金币: 833
帖子: 605
在线: 329.9小时
虫号: 2243418

[交流] 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.

回复此楼

» 收录本帖的淘帖专辑推荐

» 本帖已获得的红花（最新10朵）

weiql08

obaica

竹雨_朝颜

Sahira

» 猜你喜欢

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

查看全部散金贴

1楼 2016-11-25 20:55:42

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 luzihen 的主题更新