24小时热门版块排行榜

返回列表

【奖励】本帖被评价34次，作者redsnowolf增加金币 26.2 个

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

redsnowolf

银虫 (小有名气)

应助: 7 (幼儿园)
金币: 416.5
帖子: 213
在线: 689.8小时
虫号: 1332218

[资源] 从PROCAR中提取态密度的shell代码

压缩包中getTdos.sh、getPdos.sh和getLdos.sh分别用来计算体系的total dos、partial dos和单个原子（可指定）的

态密度，后两个需要POSCAR文件。

输入：依次输入所要计算的能量范围、费米能级、采样点数、高斯展宽，getLdos.sh还需要指定原子（默认每个原子都

计算）。

本人用的是vasp 5.2，INCAR中设置了LSORBIT = .TRUE. 和 LORBIT = 11，所得文件格式如下

PROCAR file：
PROCAR lm decomposed
# of k-points: 2925       # of bands:  28       # of ions: 2

k-point 1 : 0.00000000 0.00000000 0.00000000    weight = 0.00006400

band 1 # energy  -13.01545745 # occ.  1.00000000

ion    s    py    pz    px dxy dyz dz2 dxz dx2 tot
  1  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000
  2  0.000  0.000  0.000  0.000  0.182  0.210  0.209  0.193  0.182  0.976
tot  0.000  0.000  0.000  0.000  0.182  0.210  0.209  0.193  0.183  0.976
  1  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000
  2  0.000  0.000  0.000  0.000  0.008  0.210  0.154 -0.016 -0.008  0.348
tot  0.000  0.000  0.000  0.000  0.008  0.210  0.154 -0.016 -0.008  0.348
  1  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000
  2  0.000  0.000  0.000  0.000  0.017  0.009  0.128  0.192 -0.017  0.329
tot  0.000  0.000  0.000  0.000  0.017  0.009  0.128  0.192 -0.017  0.329
  1  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000
  2  0.000  0.000  0.000  0.000  0.181 -0.006 -0.058  0.002  0.182  0.301
tot  0.000  0.000  0.000  0.000  0.181 -0.006 -0.058  0.002  0.182  0.301

band 2 # energy  -13.01542280 # occ.  1.00000000

ion    s    py    pz    px dxy dyz dz2 dxz dx2 tot
  1  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000
  2  0.000  0.000  0.000  0.000  0.208  0.180  0.182  0.197  0.209  0.976
tot  0.000  0.000  0.000  0.000  0.208  0.180  0.182  0.197  0.209  0.976
  1  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000
  2  0.000  0.000  0.000  0.000 -0.053  0.172  0.097  0.000  0.053  0.269
tot  0.000  0.000  0.000  0.000 -0.053  0.172  0.097  0.000  0.053  0.269
  1  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000
  2  0.000  0.000  0.000  0.000 -0.198 -0.008 -0.029  0.048  0.199  0.012
tot  0.000  0.000  0.000  0.000 -0.198 -0.008 -0.029  0.048  0.199  0.012
  1  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000
  2  0.000  0.000  0.000  0.000  0.034 -0.053  0.152 -0.191  0.035 -0.023
tot  0.000  0.000  0.000  0.000  0.034 -0.053  0.152 -0.191  0.035 -0.023

band 3 # energy  -13.01541343 # occ.  1.00000000
......

如果是几个PROCAR合并到一起的，需要把k-points那儿改成总的k点数。不知道为什么PROCAR中每个band下每个原子有4组值（对应代码里的flag）。。。

POSCAR file：
POSCAR第6行需要是每种原子的数目，有些第6行是原子种类名称第7行才是数目，这样把原子名称那行删了就行了。

图中是DOSCAR中的和从PROCAR中计算的total dos比较，高斯展宽gauss=0.1 eV。

PS. 之前找了个getdos.f的Fortran程序，结果悲剧的发现完全不能用，不得已自己写了，本人Linux盲，现学的shell，希望有高手可以改得更效高些～
dos.jpg

[ Last edited by redsnowolf on 2012-9-8 at 18:53 ]

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : get_dos_from_procar_file.rar

2012-09-08 18:47:46, 4.81 K

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

第一性原理相关文档

@个人收集@

分子动力学

» 猜你喜欢

药学硕士，第一、第二作者已发表6 篇 SCI，药理方向及相关方向2026年/2027年博士申请已经有6人回复
一篇MDPI论文改变了学习工作和生活已经有5人回复
26年博士申请自荐-电催化已经有3人回复
中国地质大学（北京）博士招生补录，数理学院材料科学与工程专业和材料与化工专业已经有6人回复
收到国自然专家邀请后几年才会有本子送过来评已经有4人回复
考博已经有5人回复
26年申博自荐-计算机视觉已经有4人回复
药化及相关博士的申请已经有3人回复

1楼 2012-09-08 18:50:15

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

redsnowolf

银虫 (小有名气)

应助: 7 (幼儿园)
金币: 416.5
帖子: 213
在线: 689.8小时
虫号: 1332218

引用回帖:

3楼: Originally posted by chenweiguang at 2012-09-09 09:32:36
这相当于自己又重新做了一遍DOS
不解的是DOSCAR里面不是已经有了吗？而且vasp提供了多种方法
我看到有些人使用PROCAR重新得到DOS，这样做的意义是什么呢

算dos需要较高的k点密度，有时候体系大了计算资源又有限，不能一下子算出来，得把k点分开算，然后再通过procar得到态密度，我就遇到这种情况了。。。