24小时热门版块排行榜    

查看: 8891  |  回复: 32
【奖励】 本帖被评价28次,作者beefly增加金币 26.2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

[资源] 【beefly个人文集】用GaussView或Molden画NBO轨道

画NBO轨道,很多人会想到NBOview程序,以及后起之秀ChemCraft。或许有人还用过Molekel 4.3,不过这个不方便,还有Bug。

其实,从Gaussian 03.B开始,增加了一个新功能:把NBO存入checkpoint文件。因此,用GaussView或Molden也可以画NBO。

第一步,把NBO存入checkpoint文件。以计算HeCuF为例,输入如下。

%chk=nbo.chk
#p b3lyp/6-31g Integral(Grid=UltraFine) pop=(nboread,SAVENBOS)

title

0 1
CU    0.0   0.0  -0.3035456403
F     0.0   0.0   1.4288489797
HE    0.0   0.0  -1.9629084803

$nbo $end

在这个输入中,需要指定pop=(nboread,SAVENBOS)。SAVENBOS的作用是把NBO存入checkpoint文件(由%chk命令定义的nbo.chk),替换原来的MO。nboread用于执行NBO分析,它要求一个额外的$nbo输入行,你可以在这里加入一些感兴趣的关键词,如bndidx,plot,等。这里用了NBO的默认设置,因此也可以把nboread换成nbo,并去掉$nbo输入行。

做完这个计算后,就可以用任何版本的GaussView打开nbo.chk文件,画NBO了。细节略。如果不放心,还可以和不加SAVENBOS选项的MO图,以及用ChemCraft程序画的NBO作对比。ChemCraft用法见:

http://muchong.com/bbs/viewthread.php?tid=1585328

到这里就算完了。但是,如果你和我一样,都属于买不起GaussView、ChemCraft、NBOview的贫下中农,或者,想得到更漂亮的矢量图,那么还需要进行第二步。

第二步,把前面保存在checkpoint文件的NBO读出来,以标准格式写入Gaussian的输出文件。输入如下:

%Chk=nbo.chk
#P b3lyp/6-31g Geom=Check IOP(6/7=3) gfinput Guess=(Read,Only)

title

0  1

这里的gfinput和IOP是MOLDEN程序必须的。Geom和Guess分别用于从checkpoint获取分子结构和轨道,因此不需要定义原子坐标。还可以加上ChkBasis,可以略去基组的定义。如果第一步是自定义基组,这会很有用。

这个计算完成后,用MOLDEN打开Gaussian输出文件,就可以画NBO了。除了用MOLDEN,还可以用GabEdit,都是免费的。

参考资料:
http://www.ccl.net/chemistry/res ... .002-dir/index.html
http://www.chem.ucsb.edu/~kalju/chem226/public/task2B.html
http://www.sciencemag.org/cgi/data/322/5906/1360/DC1/1

更正说明:
pop=(nbo,SAVENBOS)没有问题。
GaussView可以用任何版本。

[ Last edited by yjcmwgk on 2010-6-16 at 21:36 ]
回复此楼

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

高斯计算总结 高斯及相关软件的使用 量化辅导 量化软件学习
Gaussian HOMO与LUMO 量子化学+实验研究 高斯日常学习总结
计算 自己感觉好的东西 量化经验 宝贝
精华 计算菜鸟 NBO分析 我的收藏

» 猜你喜欢

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

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

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

SCU清微

新虫 (初入文坛)


引用回帖:
16楼: Originally posted by sjzxbe at 2010-08-13 17:50:54
很好用的方法,与Chemcraft比对看,我觉得gaussview作图蛮漂亮的。谢谢beefly.

gaussianview打开chk文件后就能直接作NBO图了么?
32楼2017-03-01 21:45:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 33 个回答

樱花

木虫 (正式写手)


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

谢谢,学习了。
7楼2010-01-14 19:26:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

niuniu

铁虫 (初入文坛)


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

谢谢  太感谢了马上试试看啊
10楼2010-01-19 23:46:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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

谢谢!
很实用。
14楼2010-01-23 15:59:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 和导师撕破脸了 +18 圆满123456 2024-11-13 22/1100 2024-11-16 09:49 by zyhou
[基金申请] 求助!请问海外博后依托内地单位申请青基的优劣? 100+3 2280999712 2024-11-15 9/450 2024-11-15 23:08 by dxcharlary
[论文投稿] chemical science和advanced science哪个好一点 +6 yly150 2024-11-12 6/300 2024-11-15 21:06 by hcj9909053
[考博] 中南大学 粉末冶金国家重点实验室 陈超教授课题组拟招收 1~2位博士研究生通知 +7 中南大学-金属材 2024-11-14 18/900 2024-11-15 20:22 by 中南大学-金属材
[论文投稿] 职称论文 +10 xiaoqi_66 2024-11-09 13/650 2024-11-15 16:45 by sugarmei
[论文投稿] 核心期刊那个阶段最难 +4 公共茅斯 2024-11-11 4/200 2024-11-15 16:44 by sugarmei
[论文投稿] 小论文影响大论文查重吗? +3 阿道夫深度 2024-11-14 4/200 2024-11-15 16:19 by sugarmei
[基金申请] 博后面上太难了 +14 不安定因素 2024-11-14 17/850 2024-11-15 16:05 by lupbao
[教师之家] 大学老师 +11 考研一路顺风 2024-11-13 14/700 2024-11-15 15:24 by LNP@mRNA
[基金申请] 博后资助名单出来了 +12 Shxjjxjkx 2024-11-14 14/700 2024-11-15 14:40 by fanfenggui
[考博] 联系的博导说有时间聊一聊 我回复了自己可以的时间后 博导就没再回复 +4 evan0489 2024-11-15 4/200 2024-11-15 13:10 by hanweiguo
[硕博家园] reject后小感 +6 sioc-sunj 2024-11-14 8/400 2024-11-15 11:50 by 畅21
[硕博家园] 课题和专业一点都不相关怎么办 +7 imt4n9 2024-11-11 10/500 2024-11-15 01:55 by 成长虫
[硕博家园] 研究生的生活该是什么样 +4 lqy0719 2024-11-14 4/200 2024-11-14 16:45 by 阿荣喝酒
[有机交流] 同一个反应回流情况不同 20+3 1853846 2024-11-12 3/150 2024-11-14 15:42 by 太阳谷
[教师之家] 处在人生职业的分水岭 +4 otani 2024-11-13 4/200 2024-11-14 14:17 by mddzwo
[论文投稿] 要不要撤稿另投 10+4 wangzhesd 2024-11-09 8/400 2024-11-14 09:09 by LIU_V
[论文投稿] 求JACS和ACIE的最新文章模板 10+4 Jane060390 2024-11-10 4/200 2024-11-13 09:13 by 北京莱茵润色
[论文投稿] 爱思唯尔投稿系统里的通讯作者可以和文章里的通讯作者标注不同吗 +7 Omnissiah 2024-11-10 7/350 2024-11-12 14:07 by holypower
[基金申请] 要求论文发表日期在项目执行期内,论文发表日期是在线日期还是见刊日期 +6 Jiangnanyu1 2024-11-09 7/350 2024-11-10 17:49 by jurkat.1640
信息提示
请填处理意见