24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 23607  |  回复: 188
【奖励】 本帖被评价146次,作者neweroica增加金币 125.3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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 ]
回复此楼

» 收录本帖的淘帖专辑推荐

纳米技术与能源及模拟 精华网帖收集 仿真建模与计算 研究生资料集
VASP and MS 材料计算模拟实用技巧 VASP 第一性原理相关文档
先进材料与分子模拟 第一性原理和电化学 VASP使用技巧 第一性原理
VASP算例 first principle 计算科学经验总结 VASP计算
计算化学-辅助工具 initio 量化经验 第一性
VASP学习 VASP模拟 VASP相关

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

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

yulingaode

木虫 (初入文坛)


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

还没学会,细心研究
23楼2010-01-27 21:13:24
已阅   回复此楼   关注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   回复  
 支持,支持,支持
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +6 白QF 2026-03-21 7/350 2026-03-28 00:42 by 白QF
[考研] 0703化学 +6 妮妮ninicgb 2026-03-27 6/300 2026-03-27 17:18 by shangxh
[考博] 26申博 +3 加油冲啊! 2026-03-26 3/150 2026-03-27 15:38 by cls512
[考研] 安徽大学专硕生物与医药专业(086000)324分,英语已过四六级,六级521,求调剂 +4 美味可乐鸡翅 2026-03-26 4/200 2026-03-27 15:27 by 星空星月
[考研] 308求调剂 +6 墨墨漠 2026-03-27 6/300 2026-03-27 15:20 by caszguilin
[考研] 269专硕求调剂 +10 金恩贝 2026-03-21 10/500 2026-03-27 15:10 by caszguilin
[考研] 333求调剂 +3 question挽风 2026-03-23 3/150 2026-03-27 11:29 by 不吃魚的貓
[考研] 材料求调剂 +5 .m.. 2026-03-25 5/250 2026-03-27 11:08 by 不吃魚的貓
[考研] 349求调剂 +4 李木子啊哈哈 2026-03-25 4/200 2026-03-26 22:49 by fmesaito
[考研] 324求调剂 +4 wysyjs25 2026-03-21 4/200 2026-03-26 20:38 by fmesaito
[考研] 289求调剂 +17 硕星赴 2026-03-23 17/850 2026-03-26 16:18 by 不吃魚的貓
[考研] 一志愿 南京邮电大学 288分 材料考研 求调剂 +3 jl0720 2026-03-26 3/150 2026-03-26 13:39 by zzll406
[考研] 环境专硕324分求调剂推荐 +5 轩小宁—— 2026-03-26 5/250 2026-03-26 12:05 by i_cooler
[考研] 284求调剂 +15 Zhao anqi 2026-03-22 15/750 2026-03-25 12:51 by wht0531
[考研] 318求调剂 +5 plum李子 2026-03-21 8/400 2026-03-25 09:26 by aa331100
[考研] 一志愿武理085500机械专业总分300求调剂 +3 an10101 2026-03-24 7/350 2026-03-25 00:00 by 山鬼0-
[考研] 求调剂 +6 研研,接电话 2026-03-24 7/350 2026-03-24 17:01 by barlinike
[考研] 材料考研调剂生 +3 黄粱一梦千年 2026-03-24 3/150 2026-03-24 17:00 by barlinike
[考研] 305分求调剂(食品工程) +5 Sxy112 2026-03-21 7/350 2026-03-24 12:27 by 544594351
[考研] 求调剂 +3 .m.. 2026-03-21 4/200 2026-03-21 16:25 by barlinike
信息提示
请填处理意见