| 查看: 4566 | 回复: 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 |
» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有49人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有21人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有5人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有8人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有7人回复
-
filecode
已经有10人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有11人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有8人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有37人回复
-
好奇怪的filecode
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
★★★★★ 五星级,优秀推荐
|
支持分享和原创。 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
2楼2010-08-16 10:57:56
3楼2010-08-16 11:09:21
5楼2010-08-16 12:02:14
