24小时热门版块排行榜    

查看: 3054  |  回复: 17
【奖励】 本帖被评价12次,作者nkleof增加金币 10.8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

nkleof

银虫 (正式写手)


[资源] 【分享】态密度处理程序

之前发过一个简单版本的态密度处理程序,这次更新一下。

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

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

lono75

铁杆木虫 (著名写手)


★★★★★ 五星级,优秀推荐

支持原创
谢谢分享
5楼2010-08-16 12:02:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 nkleof 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 杭电、天津科技、青农和宁波工程学院如何选? +8 味道很好啊 2024-06-13 8/400 2024-06-14 21:40 by hngdxg
[基金申请] BO4的YQ答辩通知发布了吗? +4 博学笃行 2024-06-11 4/200 2024-06-14 19:44 by yifeng11
[基金申请] 2024国社科通讯评审 +6 qsd10086 2024-06-13 7/350 2024-06-14 18:02 by thesuna
[论文投稿] 投稿时忘记修改一作 +6 gll123456 2024-06-13 8/400 2024-06-14 17:07 by lihua129
[硕博家园] 关于硕博连读的一些疑问? +4 Lwenter 2024-06-14 4/200 2024-06-14 14:32 by ou0551
[考博] 申博找导师 +4 疏影横斜水清浅3 2024-06-13 6/300 2024-06-14 14:31 by zxl_1105
[有机交流] ππ堆积会发生在有机溶剂中吗 5+3 zibuyu0420 2024-06-13 4/200 2024-06-14 14:17 by 小肉干
[基金申请] 博后面上今天有bug可以看到是否资助? +19 lyfbangong 2024-06-12 30/1500 2024-06-14 12:30 by Yanchong
[基金申请] 国自然基金公布的时候基金号有吗 +8 潇洒怡惜 2024-06-13 11/550 2024-06-14 11:24 by JRfei
[基金申请] 化学会年会改成一年开一次了? +5 babu2015 2024-06-12 13/650 2024-06-14 10:46 by babu2015
[论文投稿] 返修送审后第3天,只有原来的一个审稿人接受了审稿,其他两个有没有可能跑路? 60+5 huanpo116 2024-06-08 9/450 2024-06-14 10:01 by bobvan
[教师之家] 我们学院常年位居 各学院 倒数第二。专业撤销的话,在编者有什么补偿? +12 河西夜郎 2024-06-09 13/650 2024-06-14 09:27 by dhdhgfv
[论文投稿] 文章proof要求使用机构的邮箱 5+3 不可不信缘 2024-06-11 11/550 2024-06-14 07:00 by 3001160025
[有机交流] 旋蒸请教 +6 办公室摸鱼 2024-06-09 9/450 2024-06-13 21:31 by wuxiaominggo
[基金申请] 博士后面上项目状态还是专家评审吗 10+9 Thatcheremu 2024-06-13 55/2750 2024-06-13 21:23 by 乌合麒麟
[硕博家园] 科研求助 +5 杲www 2024-06-12 6/300 2024-06-13 16:16 by 姓李名明
[基金申请] C口食品 +6 好好不过分 2024-06-09 19/950 2024-06-13 12:26 by zydemon525
[找工作] 成都产品质量检测研究院 200+3 鲸鱼663 2024-06-11 8/400 2024-06-13 08:10 by 加纳居士
[基金申请] 博后特助这周出结果吗?往年都是啥时候啊? +13 jsqy 2024-06-12 17/850 2024-06-12 19:55 by Lynn212
[论文投稿] 计算机类,请问eswa与kbs哪个相对来说容易一些啊? +3 chenrui2015 2024-06-08 4/200 2024-06-11 13:37 by p-cloud
信息提示
请填处理意见