24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 23542  |  回复: 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的回帖

valenhou001

至尊木虫 (职业作家)


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

小小的提示:vasp在输出PROCAR,将晶体波函数(WAVECAR里面)投影到原子轨道上时,并不是利用Standard basis: LANL2MB 来投影的,而是利用POTCAR里面给出原子球区域之类的原子波函数。不同原子,这个半径是不一样的。vaspmo v0.1在组合PROCAR里面的数据时,可能对所有的原子都采用同样的半径。这一点可能要注意。
53楼2010-05-22 09:28:20
已阅   回复此楼   关注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   回复  
 支持,支持,支持
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿太原理工安全工程300分,求调剂 +3 0857求调剂. 2026-03-24 3/150 2026-03-27 01:18 by 求知若渴lz
[考研] 349求调剂 +5 杰斯塔里斯 2026-03-21 5/250 2026-03-27 00:31 by wxiongid
[考研] 336材料求调剂 +7 陈滢莹 2026-03-26 9/450 2026-03-27 00:20 by wxiongid
[考研] 286求调剂 +5 PolarBear11 2026-03-26 5/250 2026-03-26 23:59 by 运气yunqi
[考研] 305求调剂 +4 哇卢卡库 2026-03-26 4/200 2026-03-26 20:58 by sanrepian
[考研] 321求调剂 +6 Ymlll 2026-03-24 6/300 2026-03-26 20:50 by 不吃魚的貓
[考研] 08开头275求调剂 +3 拉谁不重要 2026-03-26 3/150 2026-03-26 20:22 by barlinike
[考研] 271求调剂 +6 生如夏花… 2026-03-22 6/300 2026-03-26 16:48 by 张凯十八号
[考研] 考研调剂 +8 小蜡新笔 2026-03-26 8/400 2026-03-26 16:18 by dick_runner
[考研] 22 350 本科985求调剂,求老登收留 +4 李轶男003 2026-03-20 4/200 2026-03-26 16:05 by 哇啦啦啦xtj
[考研] 调剂310 +3 温柔的晚安 2026-03-25 4/200 2026-03-25 23:16 by peike
[考研] 303求调剂 +6 蓝山月 2026-03-25 6/300 2026-03-25 22:47 by 418490947
[考研] 07化学303求调剂 +5 睿08 2026-03-25 5/250 2026-03-25 22:46 by 418490947
[考研] 材料专硕 335 分求调剂 +4 拒绝冷暴力 2026-03-25 4/200 2026-03-25 18:45 by haxia
[考研] 085600材料与化工调剂 +9 A-哆啦Z梦 2026-03-23 15/750 2026-03-25 11:18 by Ainin_
[考研] 0703化学求调剂 +6 奶油草莓. 2026-03-22 7/350 2026-03-25 10:00 by shangxh
[考研] 277分求调剂,跨调材料 +3 考研调剂lxh 2026-03-24 3/150 2026-03-24 13:52 by JourneyLucky
[考研] 求调剂院校信息 +6 CX 330 2026-03-21 6/300 2026-03-22 15:25 by 无懈可击111
[考研] 336求调剂 +5 rmc8866 2026-03-21 5/250 2026-03-21 17:24 by 学员8dgXkO
[考研] 南昌大学材料专硕311分求调剂 +6 77chaselx 2026-03-20 6/300 2026-03-21 07:24 by JourneyLucky
信息提示
请填处理意见