| 查看: 4664 | 回复: 42 | ||||||
| 【奖励】 本帖被评价34次，作者redsnowolf增加金币 26.2 个 | ||||||
[资源]
从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
» 收录本帖的淘帖专辑推荐
 第一性原理相关文档 | @个人收集@ | 分子动力学 |
» 猜你喜欢
今天系统多次维护，明天很可能放榜！
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
2楼2012-09-08 18:59:38
3楼2012-09-09 09:32:36
4楼2012-09-09 10:49:50
|
本帖内容被屏蔽 |
5楼2012-09-09 11:49:27
6楼2012-09-09 12:39:14
7楼2012-09-09 12:48:26
9楼2012-09-13 10:41:03
12楼2013-05-03 15:18:03
23楼2015-08-05 17:25:38
简单回复
kingsix8楼
2012-09-11 21:42
回复
五星好评 顶一下，感谢分享！
105135818710楼
2013-03-14 19:10
回复
五星好评 顶一下，感谢分享！
105135818711楼
2013-04-03 10:13
回复
顶一下，感谢分享！
欧阳平平13楼
2013-05-04 10:05
回复
顶一下，感谢分享！
2008089014楼
2014-04-09 17:55
回复
五星好评 顶一下，感谢分享！
wangdading15楼
2014-04-24 09:34
回复
五星好评 顶一下，感谢分享！
shyzhhx1316楼
2014-05-15 15:05
回复
五星好评 顶一下，感谢分享！
哈哈呵呵哈哈17楼
2014-08-09 14:18
回复
五星好评 顶一下，感谢分享！
小康康18楼
2014-10-17 11:03
回复
五星好评 顶一下，感谢分享！
pyy11119楼
2014-11-07 21:14
回复
三星好评 顶一下，感谢分享！
ridge888820楼
2014-11-22 21:17
回复
五星好评 顶一下，感谢分享！
右?J21楼
2014-11-26 14:41
回复
五星好评 顶一下，感谢分享！
dongbaojuan22楼
2014-12-21 22:49
回复
五星好评 顶一下，感谢分享！
txwang24楼
2015-08-18 22:46
回复
五星好评 顶一下，感谢分享！
cdxy25楼
2015-10-23 09:51
回复
五星好评 顶一下，感谢分享！
人生果26楼
2015-12-08 19:42
回复
五星好评 顶一下，感谢分享！
Lelims27楼
2015-12-28 16:45
回复
五星好评 顶一下，感谢分享！
lxy671abc28楼
2016-01-29 10:16
回复
三星好评 顶一下，感谢分享！
lxy671abc29楼
2016-01-29 10:16
回复
微醉的晨光30楼
2016-03-06 22:07
回复
五星好评 顶一下，感谢分享！
bandgap.BZ31楼
2016-03-14 19:19
回复
五星好评 顶一下，感谢分享！
haizhiwuyu32楼
2016-04-24 13:26
回复
五星好评 顶一下，感谢分享！
漫天飘雪33楼
2016-06-28 10:11
回复
五星好评 顶一下，感谢分享！
haizhiwuyu34楼
2016-06-29 11:00
回复
顶一下，感谢分享！
园中园35楼
2016-07-07 17:28
回复
五星好评 顶一下，感谢分享！
小龙斯坦36楼
2016-07-12 11:46
回复
五星好评 顶一下，感谢分享！
Molly_Cao37楼
2016-07-17 16:51
回复
五星好评 顶一下，感谢分享！
kings0138楼
2016-07-19 15:44
回复
五星好评 顶一下，感谢分享！
园中园39楼
2016-07-30 18:58
回复
顶一下，感谢分享！
lvxinyuabc40楼
2016-08-21 11:22
回复
五星好评 顶一下，感谢分享！
etna41楼
2016-10-01 14:51
回复
五星好评 顶一下，感谢分享！
kunijiang42楼
2017-02-16 15:34
回复
五星好评 顶一下，感谢分享！
西南flyer43楼
2017-04-03 15:39
回复
五星好评 顶一下，感谢分享！
25