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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】怎么根据PROCAR作DOS和PDOS图
51/1返回列表
查看: 2486  |  回复: 21
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

huangyc

木虫 (正式写手)

[交流] 【求助】怎么根据PROCAR作DOS和PDOS图

大家好，计算DOS时得到PROCAR，准备学习画DOS和PDOS，共126个K点，每个K点100 bands，用lev00处理时提示缺少PROCAR或PROCAR错误，用p4vasp也打不开，请帮我看看，谢谢
PROCAR new format
# of k-points:  126         # of bands: 100         # of ions:  24

k-point    1 :    0.04166667 0.04166667 0.04166667     weight = 0.00462963

band   1 # energy  -29.57997136 # occ.  2.00000000

ion      s      p      d      f    tot
  1  0.233  0.000  0.000  0.000  0.233
  2  0.233  0.000  0.000  0.000  0.233
  3  0.233  0.000  0.000  0.000  0.233
  4  0.233  0.000  0.000  0.000  0.233
  5  0.000  0.000  0.000  0.000  0.000
  6  0.000  0.000  0.000  0.000  0.000
  7  0.000  0.000  0.000  0.000  0.000
  8  0.000  0.000  0.000  0.000  0.000
  9  0.000  0.000  0.000  0.000  0.000
10  0.000  0.000  0.000  0.000  0.000
11  0.000  0.000  0.000  0.000  0.000
12  0.000  0.000  0.000  0.000  0.000
13  0.000  0.000  0.000  0.000  0.000
14  0.000  0.000  0.000  0.000  0.000
15  0.000  0.000  0.000  0.000  0.000
16  0.000  0.000  0.000  0.000  0.000
17  0.000  0.000  0.000  0.000  0.000
18  0.000  0.000  0.000  0.000  0.000
19  0.000  0.000  0.000  0.000  0.000
20  0.000  0.000  0.000  0.000  0.000
21  0.000  0.000  0.000  0.000  0.000
22  0.000  0.000  0.000  0.000  0.000
23  0.000  0.000  0.000  0.000  0.000
24  0.000  0.000  0.000  0.000  0.000
tot  0.937  0.001  0.000  0.000  0.938
这是其中一部分，下面还有band2-100，接着就是k-points2-126，很长的文件
不知道大家都怎么处理DOS？学习VASP时间不长，还希望虫友们指点

[ Last edited by SHY31 on 2008-12-4 at 20:41 ]
回复此楼

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

VASP

» 猜你喜欢
1楼2008-11-27 22:17:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hym5701

木虫 (小有名气)
★ ★
zdhlover(金币+2,VIP+0):非常感谢，欢迎常来
引用回帖:
Originally posted by huangyc at 2008-12-7 19:38:
我用
>ifort -o dos-procar dos-procar.f90编译的。楼上的两个版本我都编译不通啊

奇怪了，你为什么要将f改成f90呢，改了当然就错了下面是我编译的信息，没有任何错误信息你看看吧
[root@localhost ~]# ls
anaconda-ks.cfg  Desktop  dos-procar.f  install.log  install.log.syslog  xcrys_tmp
[root@localhost ~]# ifort -o dos-procar dos-procar.f
[root@localhost ~]# ls
anaconda-ks.cfg  Desktop  dos-procar  dos-procar.f  install.log  install.log.syslog  xcrys_tmp
[root@localhost ~]# ./dos-procar
Spin polarized calculation? (no=1,yes=2):
1
  Enter fermi energy:
.......................................
.......

[ Last edited by hym5701 on 2008-12-7 at 21:35 ]
赞一下(10人) 回复此楼
10楼2008-12-07 21:34:54
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 22 个回答

intro987

铁杆木虫 (小有名气)
★ ★
jghe(金币+2,VIP+0):感谢你参与讨论！
我使用DOSCAR处理的，应该和PROCAR是一起出来的结果，里面很多列数据，第一列是能量，后面对应s, p_y, p_z, p_x, d_xy, d_yz, d_z2, d_xz, d_x2-y2轨道（VASP出来的结果应该顺序都是一样的）。数据上下顺序对应POSCAR里面的原子顺序
lev00是处理电荷密度的，p4vasp我只用来打开过POSCAR，没试过PROCAR，可能不支持这种类型吧~
赞一下(10人) 回复此楼
2楼2008-12-04 18:43:48
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hym5701

木虫 (小有名气)
★ ★
jghe(金币+2,VIP+0):谢谢你提供的资源！
引用回帖:
Originally posted by huangyc at 2008-11-27 22:17:
大家好，计算DOS时得到PROCAR，准备学习画DOS和PDOS，共126个K点，每个K点100 bands，用lev00处理时提示缺少PROCAR或PROCAR错误，用p4vasp也打不开，请帮我看看，谢谢
PROCAR new format
# of k-points:  126   ...

网上不是用vasp的一些tools吗，里面就有一个分析PROCAR的文件，文件名是dos-procar.f，你现将这个文件编译一下，然后就运行这个文件，他会提示你要干什么，运行完后，就可以得到totdos。dat 和ldos。dat，接着用origin等画图就可以了。
赞一下(10人) 回复此楼
3楼2008-12-04 20:22:30
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangyc

木虫 (正式写手)
楼上的这个东东编译不起来啊，有语法错误
哪位再指点，谢谢
赞一下 回复此楼
4楼2008-12-05 11:15:55
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 22 个回答
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭