24小时热门版块排行榜    

Znn3bq.jpeg
查看: 23700  |  回复: 188
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

neweroica

木虫 (著名写手)


[资源] 【neweroica个人文集】原创程序 vaspmo (v0.0) ――可视化VASP的分子轨道



注意:

程序在转换坐标时会有错误!这个bug已经在即将发布的新版本v0.1里更正。

最新消息:

明后天将推出vaspmo v0.1版,届时大家可以又任何轨道可视化软件来绘制VASP计算的轨道。

我先把Chemcraft和GaussView可视化VASP轨道的效果图贴上来(见附件中新添的两张图片。)




终于弄成功了!!!

现在我们可以用Molekel来绘制VASP计算的分子轨道了!

这个想法是我不久前产生的(详见  http://muchong.com/bbs/viewthread.php?tid=1787325 ):

VASP可以把波函数分解成中心在各个原子的原子轨道(在INCAR中使用LORBIT =2或LORBIT =12),于是我们就能得到近似的分子轨道。轨道系数保存在PROCAR文件中。于是我就有了一个想法:将PROCAR文件转化成Molekel能读的文件格式,这样我们就能可视化VASP投影出的分子轨道了!

现在,我终于写好了整个程序,取名为vaspmo。目前适用于周期表中原子序数从1到83的元素(La系元素除外)。可将任何K点的所有能带对应的晶体轨道全部转换出来。

在Linux下编译使用如下命令:
gcc -o vaspmo -lm vaspmo.c

在Windows下可以使用任何C或C++编译器编译(因为我用标准C编写的)。
当然,附件中也附带了我在Windows下编译好来的vaspmo.exe可执行文件。

使用方法,在VASP计算输出的目录下(必须有PROCAR和CONTCAR文件)执行Ñ
./vaspmo
然后按照程序提示输入第几个K点即可。

程序输出示例如下:

vaspmo v0.0 --  Written by ...

Use "vaspmo -h" for help and more information.

Information obtained from CONTCAR:
    There are 54 atoms, 4 ion types:
        6 C atoms
        4 H atoms
        4 S atoms
        40 Cu atoms
    Lattice vectors (Cartesian, in Angstrom):
            5.147560            0.000000            0.000000
            0.000000           12.868895            0.000000
            0.000000            0.000000           30.000000

Cartesian coordinates have been written into VASPMO.g03.

Basis sets information has been written into VASPMO.g03.

Information obtained from PROCAR:
    There are 5 k-points and 320 bands.

Please choose k-point # (1-5):1

Reading phase factors from PROCAR...
################################################################
Done.


程序结束后就生成VASPMO.g03文件,用Molekel打开就能可视化分子轨道。

具体注意事项见README.txt或源文件的开头部分。

如果你在发表文章或做报告时用到了我的这个vaspmo程序,敬请注明。

有任何建议或bug,请email我(地址在程序里)。


附图是我用Molekel画的VASP计算得到的分子轨道:
1. 苯分子的HOMO和LUMO
2. 吸附在金属表面的分子与金属的一条成键轨道。

[ Last edited by yjcmwgk on 2010-6-16 at 21:41 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

紫檀色

新虫 (小有名气)


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

感谢啦!!!!!
142楼2013-03-13 22:52:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 189 个回答

rfvuhb

铜虫 (小有名气)


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

原创,了不起!!!
感谢!!
4楼2010-01-26 08:39:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cationly

木虫 (正式写手)


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

★★★★★ 五星级
5楼2010-01-26 09:22:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雁儿霏霏

木虫 (正式写手)


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

楼主果然说到做到,严重支持
6楼2010-01-26 10:38:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
allawwp10楼
2010-01-26 19:11   回复  
 支持,支持,支持
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料299专硕求调剂 +17 +21 2026-04-09 17/850 2026-04-15 17:20 by zs92450
[考研] 药学305求调剂 +7 玛卡巴卡boom 2026-04-11 7/350 2026-04-15 13:21 by 西北望—风沙
[考研] 药学求调剂 +11 喽哈加油 2026-04-14 13/650 2026-04-14 21:14 by qingfeng258
[考研] 考研求调剂 +6 ban班小七 2026-04-11 6/300 2026-04-14 14:06 by 哆啦A梦只是个梦
[考研] 0856专硕求调剂 希望是a区院校 +24 好好休息好不好 2026-04-09 27/1350 2026-04-13 22:22 by pies112
[考研] 2本,初试303,0860求调剂 +6 floriea 2026-04-12 8/400 2026-04-12 18:13 by zhouxiaoyu
[考研] 291求调剂 +11 关忆北. 2026-04-09 12/600 2026-04-12 10:32 by 逆水乘风
[考研] 一志愿华中农微生物,288分,三年实验经历 +11 代fish 2026-04-09 11/550 2026-04-12 10:21 by Hayaay
[考研] 一志愿西北工业大学289 085602 +33 yang婷 2026-04-10 34/1700 2026-04-12 08:11 by Art1977
[找工作] 山东高校教师考核超级无底线,员工过不下去啦 +4 qut2026 2026-04-09 9/450 2026-04-12 00:54 by qut2026
[考研] 一志愿厦大0856,306求调剂 +15 Bblinging 2026-04-11 15/750 2026-04-11 22:53 by 314126402
[考研] 调剂 +6 青灯不负 2026-04-09 6/300 2026-04-11 20:35 by dongdian1
[考研] 求调剂,一志愿大连理工大学354分 +5 雨声余生 2026-04-11 6/300 2026-04-11 16:12 by 雨声余生
[考研] 调剂 +5 文道星台 2026-04-11 5/250 2026-04-11 15:01 by 凯凯要变帅
[考研] 调剂 +19 小张ZA 2026-04-10 20/1000 2026-04-10 22:08 by 猪会飞
[考研] 289 分105500药学专硕求调剂(找B区学校) +6 白云123456789 2026-04-09 8/400 2026-04-10 21:13 by zhouxiaoyu
[考研] 263能源动力专硕求调剂 +4 加大号饭盒袋 2026-04-10 4/200 2026-04-10 20:52 by gong120082
[考研] 一志愿京区985,085401电子信息,本科电子信息 +3 阳光开朗的男孩 2026-04-10 3/150 2026-04-10 16:29 by sophia_93
[考研] 085800 能源动力求调剂 +6 阿biu啊啊啊啊啊 2026-04-10 6/300 2026-04-10 15:03 by hemengdong
[考研] 348求调剂 +3 candyyyi 2026-04-09 3/150 2026-04-09 17:20 by 段伟艳
信息提示
请填处理意见