24小时热门版块排行榜    

查看: 4206  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 贺仪 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 娱乐 +3 Tide man 2026-08-03 3/150 2026-08-04 09:58 by gltch
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 17/850 2026-08-04 09:46 by archvillain
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 5/250 2026-08-04 09:27 by v1TJC4qJauAB
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 09:17 by v1TJC4qJauAB
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 cu9Nq1xK233Z 2026-08-03 10/500 2026-08-04 09:15 by v1TJC4qJauAB
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +12 winnerche 2026-07-29 16/800 2026-08-04 08:18 by Equinoxhua
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 10/500 2026-08-04 08:05 by 88817753
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 07:37 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 13/650 2026-08-04 07:30 by k7OM8YghWbkC
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 6/300 2026-08-04 06:10 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[基金申请] 微信指数没变化,科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
请填处理意见