24小时热门版块排行榜    

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

tubaxu

金虫 (小有名气)


引用回帖:
Originally posted by neweroica at 2010-04-05 22:59:56:


INCAR中使用LORBIT =2或LORBIT =12

多谢楼主回复,我在INCAR里面随便一个位置(不知应该不应该),填了一行LORBIT=2,运行后,仍然没有出现PROCAR文件。
   另外运行./vaspmo文件,错误提示为:
Error:found 0 atom for ion type 1 in CONTCAR(Line 7).

不知什么原因,是不是我的INCAR文件不够规范呢?
49楼2010-04-06 12:50:03
已阅   回复此楼   关注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   回复  
 支持,支持,支持
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 333求调剂 +6 wfh030413@ 2026-03-23 6/300 2026-03-26 22:45 by 学员8dgXkO
[考研] 一志愿武汉理工,总分321,英一数二,求老师收留。 +3 nnnnnnn5 2026-03-25 3/150 2026-03-26 22:38 by 不吃魚的貓
[考研] 324求调剂 +3 hanamiko 2026-03-26 3/150 2026-03-26 21:00 by sanrepian
[考研] 08开头275求调剂 +3 拉谁不重要 2026-03-26 3/150 2026-03-26 20:22 by barlinike
[考研] 材料调剂 5+4 想要一壶桃花水 2026-03-25 10/500 2026-03-26 19:56 by 不吃魚的貓
[考研] 312求调剂 +8 上岸吧ZJY 2026-03-22 12/600 2026-03-26 18:49 by muchong357
[考研] 274求调剂 +14 顾九笙要谦虚 2026-03-24 20/1000 2026-03-26 18:15 by 乐呵呵的追梦人
[考研] 一志愿 西北大学 总分282 英语一62 求调剂 +7 18419759900 2026-03-25 7/350 2026-03-26 16:07 by 不吃魚的貓
[考研] 22 350 本科985求调剂,求老登收留 +4 李轶男003 2026-03-20 4/200 2026-03-26 16:05 by 哇啦啦啦xtj
[考研] 07化学303求调剂 +5 睿08 2026-03-25 5/250 2026-03-25 22:46 by 418490947
[考研] 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
[考研] 一志愿北化315 求调剂 +3 akrrain 2026-03-24 3/150 2026-03-24 19:35 by 了了了了。。
[考研] 085404电子信息284分求调剂 +4 13659058978 2026-03-24 4/200 2026-03-24 12:15 by syl20081243
[考研] 341求调剂(一志愿湖南大学070300) +5 番茄头--- 2026-03-22 6/300 2026-03-23 23:45 by Txy@872106
[考研] 336求调剂 +4 收到VS 2026-03-20 4/200 2026-03-23 19:02 by macy2011
[考研] 275求调剂 +6 shansx 2026-03-22 8/400 2026-03-22 15:27 by barlinike
[考研] 求调剂院校信息 +6 CX 330 2026-03-21 6/300 2026-03-22 15:25 by 无懈可击111
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
信息提示
请填处理意见