24小时热门版块排行榜    

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

neweroica

木虫 (著名写手)


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

nice
3楼2010-01-06 21:35:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 beefly 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 论文返修状态变成了awaiting AE recommendation 10+4 猪小耍 2024-11-13 15/750 2024-11-16 09:24 by guxue
[基金申请] 求助!请问海外博后依托内地单位申请青基的优劣? 80+3 2280999712 2024-11-15 9/450 2024-11-15 23:08 by dxcharlary
[教师之家] 这种人痛苦吗 +12 2671 2024-11-12 14/700 2024-11-15 22:15 by 仇恨毁灭财富
[教师之家] 评正教授需要两个国家级项目,有人用子课题糊弄,结果在评审前资格公示时被举报拿下了 +22 瞬息宇宙 2024-11-12 31/1550 2024-11-15 22:02 by 潇湘之迷
[论文投稿] 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 中南大学-金属材
[论文投稿] angew三个小修,返修过去已经20多天了没动静,有朋友遇到这种情况吗? +5 迟迟未到场 2024-11-14 5/250 2024-11-15 20:16 by 秦时明月zy
[论文投稿] 小论文影响大论文查重吗? +3 阿道夫深度 2024-11-14 4/200 2024-11-15 16:19 by sugarmei
[基金申请] 博后资助名单出来了 +12 Shxjjxjkx 2024-11-14 14/700 2024-11-15 14:40 by fanfenggui
[有机交流] 二甲胺的使用 20+5 太阳谷 2024-11-14 13/650 2024-11-15 13:11 by xiaomei1031
[硕博家园] reject后小感 +6 sioc-sunj 2024-11-14 8/400 2024-11-15 11:50 by 畅21
[基金申请] 求助 +4 Enenenene 2024-11-15 4/200 2024-11-15 09:36 by 榨菜拌青椒
[论文投稿] OE返修遇到expired +3 隔壁老王来了 2024-11-14 7/350 2024-11-14 20:46 by 隔壁老王来了
[论文投稿] 投稿系统中的通讯作者和文章中的通讯作者不一样,文章目前被录用了? +4 babybabygo 2024-11-12 5/250 2024-11-14 19:13 by 走了002
[有机交流] 同一个反应回流情况不同 20+3 1853846 2024-11-12 3/150 2024-11-14 15:42 by 太阳谷
[硕博家园] 大龄已婚想读博如何 +15 米娅阳 2024-11-11 18/900 2024-11-14 14:18 by xiaomi0401
[论文投稿] 求JACS和ACIE的最新文章模板 10+4 Jane060390 2024-11-10 4/200 2024-11-13 09:13 by 北京莱茵润色
[基金申请] 浙江省青基怎么看A还是B? +8 ToallGoof 2024-11-11 14/700 2024-11-13 09:11 by Praybest
[论文投稿] 投稿求助 40+3 liukun11 2024-11-12 4/200 2024-11-12 19:13 by liukun11
[基金申请] 要求论文发表日期在项目执行期内,论文发表日期是在线日期还是见刊日期 +6 Jiangnanyu1 2024-11-09 7/350 2024-11-10 17:49 by jurkat.1640
信息提示
请填处理意见