24小时热门版块排行榜    

查看: 4302  |  回复: 10

贺仪

铁杆木虫 (著名写手)

[交流] 【求助】vasp怎么计算PDOS? 已有1人参与

vasp怎么计算PDOS?INCAR怎么设置呢?结果输出在哪个文件里?谢谢咯
回复此楼
嘟啦啦嘟啦啦嘟
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

sunyang1988

至尊木虫 (著名写手)

★ ★
贺仪(金币+3):说道要点了! 2010-10-12 10:53:49
zzy870720z(金币+2):谢谢详细解答 2010-10-12 14:10:18
不知道你的图是怎么画的,请说清楚一点
    -25.406  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -25.300  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -25.193  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -25.086  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -24.979  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -24.872  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -24.765  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -24.658  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -24.552  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -24.445  0.1097E-10  0.1100E-10  0.1173E-11  0.1176E-11
    -24.338  0.1766E-06  0.1770E-06  0.1888E-07  0.1892E-07
    -24.231  0.3092E-03  0.3097E-03  0.3306E-04  0.3311E-04
    -24.124  0.6140E-01  0.6146E-01  0.6594E-02  0.6601E-02
第一列是能量,第二列是对应该能量自旋向上的态密度第三列是自旋向下的,第四五列分别是自旋向上与自旋向下态密度的积分,你好像画了一个态密度随能量分布还有一个态密度积分随能量分布
这是DOSCAR中的第一部分,是总的态密度。如果你设置了LORBIT=10-12,后面还有每个原子周围的Local DOS,用个小程序把他们分割,然后可以作图
心之所向无惧无悔
6楼2010-10-12 10:49:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

sunyang1988

至尊木虫 (著名写手)

你说的PDOS是什么意思,局域在每个原子周围的Local DOS,还是某个能量范围内的DOS,或者是某几个能带的DOS
心之所向无惧无悔
2楼2010-10-12 09:21:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

引用回帖:
Originally posted by sunyang1988 at 2010-10-12 09:21:56:
你说的PDOS是什么意思,局域在每个原子周围的Local DOS,还是某个能量范围内的DOS,或者是某几个能带的DOS

这个东东:projected density of states
嘟啦啦嘟啦啦嘟
3楼2010-10-12 10:04:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)

★ ★
zzy870720z(金币+2):谢谢指教 2010-10-12 14:09:53
引用回帖:
Originally posted by 贺仪 at 2010-10-12 10:04:34:



这个东东:projected density of states

就是局域在每个原子周围的态密度吧
在INCAR中设置LORBIT=10-12
10每个原子局域态密度,按照量子数l分解,即每个原子周围的spd轨道投影态密度
11同上,在10的基础上给出lm分解,即s,pxpypz这样的波函数投影
12不怎么懂,貌似是11基础上加上相位信息,好像用的不多
最后给出的态密度都在DOSCAR里,用个小程序分解就可以画图了
心之所向无惧无悔
4楼2010-10-12 10:12:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

引用回帖:
Originally posted by sunyang1988 at 2010-10-12 10:12:43:

就是局域在每个原子周围的态密度吧
在INCAR中设置LORBIT=10-12
10每个原子局域态密度,按照量子数l分解,即每个原子周围的spd轨道投影态密度
11同上,在10的基础上给出lm分解,即s,pxpypz这样的波函数投影
...

我是设置了 但是得到
  50  50   0   0
  0.4190046E+02  0.1229999E-08  0.1229999E-08  0.1598999E-08  0.5000000E-15
  0.000000000000000E+000
  CAR
graphite                                
      6.65092865    -25.40643856  301     -1.76092371      1.00000000
    -25.406  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -25.300  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -25.193  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -25.086  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -24.979  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -24.872  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -24.765  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -24.658  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -24.552  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -24.445  0.1097E-10  0.1100E-10  0.1173E-11  0.1176E-11
    -24.338  0.1766E-06  0.1770E-06  0.1888E-07  0.1892E-07
    -24.231  0.3092E-03  0.3097E-03  0.3306E-04  0.3311E-04
    -24.124  0.6140E-01  0.6146E-01  0.6594E-02  0.6601E-02


后面几行意思不是很明白画出的图
嘟啦啦嘟啦啦嘟
5楼2010-10-12 10:22:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

引用回帖:
Originally posted by sunyang1988 at 2010-10-12 10:49:27:
不知道你的图是怎么画的,请说清楚一点
    -25.406  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -25.300  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -25.193  0.0000E+00  0.0000E+00   ...

非常感谢的你的回答,问题是我的doscar里没有后面的东西了……


不好意思是我的lorit这个词写错了……   搞不懂怎么没提示错误呢?

还有就是 得到这些后 怎么看体系的磁性呢?

[ Last edited by 贺仪 on 2010-10-12 at 10:59 ]
嘟啦啦嘟啦啦嘟
7楼2010-10-12 10:53:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)


ym23(金币+1):欢迎常来交流 2010-10-12 15:59:33
输错了就被系统无视了,无效语句不会有反应的
心之所向无惧无悔
8楼2010-10-12 10:59:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

引用回帖:
Originally posted by sunyang1988 at 2010-10-12 10:59:01:
输错了就被系统无视了,无效语句不会有反应的

还有就是 得到这些后 怎么看体系的磁性呢?
嘟啦啦嘟啦啦嘟
9楼2010-10-12 10:59:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)

贺仪(金币+2): 2010-10-14 22:15:16
引用回帖:
Originally posted by 贺仪 at 2010-10-12 10:59:56:



还有就是 得到这些后 怎么看体系的磁性呢?

在你的另一个帖子中回答了,看OSZICAR
心之所向无惧无悔
10楼2010-10-12 11:08:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 贺仪 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +5 process2012 2026-08-23 6/300 2026-08-25 08:44 by 白大多亮
[基金申请] 只有每年这种时候来逛逛小木虫 +25 yaoyewhu2008 2026-08-20 27/1350 2026-08-25 08:01 by Equinoxhua
[基金申请] 明天应该可查了!? +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +7 scalable 2026-08-24 7/350 2026-08-24 23:01 by anata1209
[基金申请] 没有任何消息-是不是就凉了 +7 图啦图啦 2026-08-24 8/400 2026-08-24 22:00 by maomao_da
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 19/950 2026-08-24 16:52 by iaeyuan
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见