24小时热门版块排行榜    

查看: 27825  |  回复: 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的回帖

mei1988925

金虫 (小有名气)


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

楼主很强势!学习了!
72楼2011-04-02 08:13:43
已阅   回复此楼   关注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   回复  
 支持,支持,支持
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 FChKX8meu1ds 2026-11-09 7/350 2026-11-11 13:51 by cdqGD6ftHmbl
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 cqv1zKnzd7hr 2026-11-09 8/400 2026-11-11 13:39 by cdqGD6ftHmbl
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 cqv1zKnzd7hr 2026-11-09 7/350 2026-11-11 13:33 by cdqGD6ftHmbl
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 C6CRn66RYYXS 2026-11-10 3/150 2026-11-11 11:41 by gRND4NVAZ8VG
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 FChKX8meu1ds 2026-11-09 5/250 2026-11-11 10:46 by 1Nb0IVlmV6v0
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 FChKX8meu1ds 2026-11-09 5/250 2026-11-11 10:43 by 1Nb0IVlmV6v0
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 cqv1zKnzd7hr 2026-11-09 6/300 2026-11-11 10:32 by 1Nb0IVlmV6v0
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 zSbvVX9JFl9J 2026-11-09 5/250 2026-11-11 10:10 by 1Nb0IVlmV6v0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 G2rhSvqvhNzl 2026-11-10 4/200 2026-11-11 08:34 by 1Nb0IVlmV6v0
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 07:06 by Z7qfCsdIiD1T
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 01:10 by Z7qfCsdIiD1T
[考博] 申请环境,微生物方向博士 +3 浅笑09 2026-11-08 3/150 2026-11-10 12:44 by 可淡不可忘
[考博] 27博士申请 +3 白白白ly 2026-11-09 3/150 2026-11-09 21:39 by szq0819
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 5/250 2026-11-09 16:53 by 张丹枫1991
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 5sJFDKp2kwuI 2026-11-06 14/700 2026-11-09 11:52 by GwrQVVR6Oenx
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
信息提示
请填处理意见