| 查看: 4123 | 回复: 17 | ||||||
| 【奖励】 本帖被评价12次，作者nkleof增加金币 10.8 个 | ||||||
[资源]
【分享】态密度处理程序
|
||||||
|
之前发过一个简单版本的态密度处理程序，这次更新一下。 split_dos.ksh是根据JMS NRL的脚本改的，其它程序是我自己写的。 用法：只能在linux系统下，键入各个程序名称运行就可以，如果遇到平台问题，sumdos程序无法执行，可以用 gcc -lm sum.c 编译源文件产生可执行程序。 运行split_dos.ksh会读取DOSCAR，并生成总的态密度DOS0和各个原子的各轨道态密度DOS1~DOSn。（这里注意在计算DOS时INCAR文件里的NEDOS不要设置为>9999，五位数的NEDOS值会造成错误） 然后运行sumdos可以产生按照特定原子特定轨道模式生成的态密度和（注意：适用于ISPIN=2的计算结果，即DOSCAR里面共有有20列数据）。 例： sumdos name -t 1 2 3 sumdos是程序名称，name是一个自定义字符串，会加在生成的文件名里，可以作为说明（但是长度只有8个字符），-t表示生成各个轨道的态密度和，（还有两种模式: -n 表示不加和各轨道态密度，仍旧生成从s到d各个轨道的态密度； -d 表示把p轨道，eg轨道，t2g轨道各自加和。！！仍然注意，和上次一样，这里的eg指的是第14、15和18、19列数据，t2g是剩下的3列数据，如果有各个d轨道不是按照这样排列的情形则生成的就不是eg和t2g），后面的数字表示要加和的原子标号。上面的例子就是把第1、2、3号原子的所有轨道的态密度全部加起来。程序支持采用 n~m 形式表示要处理的第n到m个标号连续的原子，例如要处理13到26号原子时，不需要一个一个写出原子标号，只需要写 13~26 就可以。 附送一个将CONTCAR文件转换为可导入Materials Studio的cif文件的程序，使用很简单，只需要在linux下，vasp任务的计算文件夹下，运行exe_Ctocif，程序会自动寻找CONTCAR和POTCAR，并生成cif文件。另外也支持在只有POSCAR文件时转换为cif，但这时需要手动输入元素名称。（！！注意：如果是vasp5版本的CONTCAR，请先删除其中的元素符号那一行。） 例如： exe_Ctocif exe_Ctocif POSCAR Si O 这两种方式都可以，前一种自动从赝势里找元素名称，后一种手动输入元素名。 http://ishare.iask.sina.com.cn/f/9316742.html [ Last edited by nkleof on 2010-8-16 at 11:10 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 材料计算模拟实用技巧 | 材料计算模拟软件 | 理论计算及软件的使用 | Materials-studio |
» 猜你喜欢
一志愿985，本科211，0817化学工程与技术319求调剂
已经有4人回复
-
311求调剂
已经有6人回复
-
0703化学调剂 290分有科研经历，论文在投
已经有4人回复
-
化学工程321分求调剂
已经有6人回复
-
本人考085602 化学工程 专硕
已经有9人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有7人回复
-
326求调剂
已经有3人回复
-
焦虑
已经有8人回复
-
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- linux shell循环的问题
已经有3人回复
- VASP计算中，分波态密度问题
已经有6人回复
- vasp计算含有磷原子结构的态密度时磷原子为什么还有d态
已经有14人回复
- VASP 计算DOS，态密度中费米能的作用
已经有3人回复
- 用什么软件把三张态密度图画在一个图中？
已经有11人回复
- 总的态密度为什么和分波态密度不相等？
已经有19人回复
- 求处理vasp能带的小程序
已经有14人回复
- ENCUT 的脚本程序怎么用
已经有15人回复
- vasp计算结果处理态密度要用到的程序
已经有18人回复
- VASP计算态密度出现问题
已经有12人回复
- 【求助】总态密度与分波态密度
已经有9人回复
- 【求助】CASTEP中分析态密度
已经有7人回复
- 【求助】求VASP计算ZnO能带和态密度的详细输入设置
已经有15人回复
- 【求助】态密度中6d贡献部分丢失
已经有3人回复
- 【求助】VASP计算电子态密度
已经有34人回复
- 【求助】关于锂原子vasp计算后的1s2s2p分波态密度处理
已经有12人回复
17楼2014-03-27 11:27:59
2楼2010-08-16 10:57:56
3楼2010-08-16 11:09:21
4楼2010-08-16 11:15:26
5楼2010-08-16 12:02:14
★★★★★ 五星级,优秀推荐
|
支持分享和原创。 split_dos源代码中的逻辑关系理清楚没？form=0/1， 及 # Total is first DOS if [ $form -eq 0 ]; then i=0 else i=0 fi 这个没理解为什么要这样写。 还有后面的 while [ $i -le 0 ]; do #echo $i $start $end if [ $i -gt 0 ]; then sed -n ''$i' p' tmp.dat2 >| DOS$i fi if [ $i -gt 0 ]; then 这个循环从i<=0开始，但是上面无论怎样，看起来i的初值是0呢。为什么这样做循环呢？奇怪，好像可以大量简化，即删掉while [ $i -le 0 ]这段循环。 你说呢？ |
6楼2010-08-16 14:38:36
7楼2010-08-16 14:48:47
8楼2010-08-16 15:49:42
9楼2010-08-16 15:59:08
14楼2011-04-15 17:10:28
15楼2011-04-15 17:16:34
16楼2011-04-15 18:47:27
简单回复
zzy870720z10楼
2010-08-16 17:24
回复
wuli811楼
2010-09-14 16:15
回复
wonderflyer12楼
2010-09-15 09:02
回复
aylayl0813楼
2010-09-15 14:45
回复
laoda7918楼
2016-04-20 11:26
回复
一般
20