| 查看: 5497 | 回复: 19 | |||||
| 【奖励】 本帖被评价8次，作者nkleof增加金币 7.5 个 | |||||
[资源]
【分享】处理DOSCAR文件的脚本和程序
|
|||||
|
根据JMS NRL的split_dos脚本稍做修改，使得能够处理从s到d所有18个轨道（包括自旋上下）的DOSCAR，生成每个原子的轨道态密度，然后自己写了一个C程序，可以把同类原子的p轨道和d的t2g，eg轨道分别求和生成新的轨道态密度文件。适用于对含过渡金属元素体系的电子态密度进行处理，可以很方便的得到p轨道和t2g，eg轨道的总态密度。 使用方法： 仅限linux平台 1、把三个脚本程序解压后放到一个/bin目录并把路径添加到环境变量里（或者把三个文件都放到要处理的数据所在目录）。 2、键入split_dos.ksh运行分解DOSCAR的脚本，得到体系总的和每个原子的分波态密度。 3、键入sumdos *1 *2 *3执行求p以及t2g，eg轨道态密度和，其中*1, *2, *3为三个参数，依次为要处理的元素符号（如Fe），起始原子序号，最终原子序号（比如Fe原子在POSCAR里面是第5到13号原子，则*2和*3分别为5 13）。 http://d.namipan.com/d/bd4ee2bb9 ... 5c2334a0697fb150000 注意！！！！：sumdos的使用仅限于t2g和eg分别简并得很好的情况，即它们内部不再分裂，因为程序只是简单的把dxy, dyz, dzx和dz2, dx2分别加起来！如果有较强退简并情况，如晶体形变造成t2g和eg再分裂，请勿用此程序！ [ Last edited by nkleof on 2010-5-24 at 11:14 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 纳米技术与能源及模拟 | VASP | 各种各样的软件 |
» 猜你喜欢
一志愿武理材料305分求调剂
已经有5人回复
-
材料专硕英一数二306
已经有3人回复
-
085600材料与化工调剂 324分
已经有7人回复
-
311求调剂
已经有3人回复
-
26调剂/材料/英一数二/总分289/已过A区线
已经有7人回复
-
295求调剂
已经有5人回复
-
354求调剂
已经有5人回复
-
材料专业求调剂
已经有3人回复
-
求材料调剂
已经有7人回复
-
一志愿天大材料与化工（085600）总分338
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教Windows系统下perl脚本的ARGV应用问题，谢谢
已经有5人回复
- band dos 小程序的一些疑问？？？
已经有19人回复
- 【活动】Show 脚本，赢奖励
已经有23人回复
- 请教ORCA的并行及集群上的脚本文件
已经有14人回复
- 【分享】态密度处理程序
已经有17人回复
- 【求助】求将vasp的输出文件转化成xyz坐标的脚本
已经有11人回复
- 【求助】VASP计算电子态密度
已经有34人回复
- 【原创】将POSCAR转化为cif文件的一个脚本
已经有145人回复
2楼2010-05-24 11:42:09
3楼2010-05-25 09:09:23
4楼2010-05-25 09:13:25
5楼2010-05-25 10:05:40
7楼2010-05-25 22:31:35
8楼2010-05-26 07:10:32
9楼2010-05-26 08:09:31
10楼2010-05-26 08:28:46
11楼2010-05-26 10:34:34
★ ★
zzy870720z(金币+2):感谢解答 2010-05-26 17:59:34
zzy870720z(金币+2):感谢解答 2010-05-26 17:59:34
|
我看过他的回帖，实际上就是，假如把dxy, dyz, dzx, dx2, dz2看作5个线性无关项，那么当vasp所采用的坐标系不和所定义的x, y, z重合时，输出的5个轨道可以看做是这5个轨道的线性组合。我不确定vasp是不是会采用一套很特别的坐标系，但是我通常建模型都是用内坐标也就是direct模式，所以得到的结果应该是xyz就是晶体场坐标。 我写的那个sumdos程序只是对特定的几列进行简单加和，所以只适用于特定情况，否则不保证加出来的t2g和eg正确。我在尖晶石结构中使用过，是适用于内坐标建模的尖晶石结构的。 因此提醒注意，如果使用那个程序，应在使用前看看10和11，12和13，16和17这三组列数据（分别对应自旋向上和向下）是不是对应简并的t2g三个轨道态密度，14和15，18和19这两组列数据是不是对应eg轨道态密度。 [ Last edited by nkleof on 2010-5-26 at 17:14 ] |
12楼2010-05-26 14:55:55
13楼2010-05-26 22:15:33
14楼2010-05-27 10:31:22
15楼2010-07-30 15:12:08
16楼2010-09-02 17:03:12
17楼2010-09-05 20:08:11
18楼2010-09-06 10:17:32
19楼2010-10-09 16:14:06
20楼2018-07-10 21:38:29
简单回复
2010-05-25 14:07
回复
赞。。。。。
5