24小时热门版块排行榜    

查看: 9038  |  回复: 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 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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

谢谢!
很实用。
14楼2010-01-23 15:59:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 beefly 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[功能材料] 请教大神,带隙是DFT计算出来的还是实验测出来的? 10+3 大力2010 2024-12-21 6/300 2024-12-22 18:08 by Ayase_lumi
[考研] 出中科院一区文章,睋:99.839.508,董徳荃董 +5 hpgqu17 2024-12-22 5/250 2024-12-22 18:02 by qsx04152006
[论文投稿] 挑数据是造假吗 +8 电话的建设 2024-12-21 8/400 2024-12-22 18:01 by beefly
[找工作] 化学博士 +5 实名心动国民萌 2024-12-17 7/350 2024-12-22 15:33 by 大发财树
[硕博家园] 寻45-50的博士教授 +5 大发财树 2024-12-22 5/250 2024-12-22 13:31 by dh6ur@31
[教师之家] 在学校是提拔成处长难度大还是教授难度大 +9 akslis2024 2024-12-17 12/600 2024-12-22 12:32 by Quakerbird
[找工作] 柔引引进博士硕士。 +17 大发财树 2024-12-18 17/850 2024-12-22 11:28 by syl200707
[硕博家园] “认怂了”还是“认清了”:青年博士纷纷“逃离”科研的背后…… +14 柏舟0000 2024-12-16 15/750 2024-12-22 10:50 by 大发财树
[考博] 读博的欲望达到了巅峰 +15 小丸九 2024-12-17 15/750 2024-12-22 10:30 by ming9871
[硕博家园] 读研之后发现: 圈子变得越来越小了... +13 丁柯翔鸭 2024-12-19 17/850 2024-12-22 10:28 by ming9871
[论文投稿] 投稿意见求助,没弄清回答的方向 8+4 moonlig 2024-12-18 4/200 2024-12-22 08:51 by steven_198377
[论文投稿] 投稿期刊 5+4 王丽媛西西西 2024-12-17 7/350 2024-12-22 04:44 by yudaoqian88
[论文投稿] sci二区的审稿周期 22+5 分外含晴 2024-12-16 12/600 2024-12-21 23:02 by 分外含晴
[论文投稿] MDPI投稿连续拒稿是什么原因?文章质量不好吗?还是其他原因? 2+8 多听多看多学 2024-12-17 17/850 2024-12-21 18:51 by 6681981668
[硕博家园] 招聘博士 +7 大发财树 2024-12-18 8/400 2024-12-21 16:17 by yanjiaming
[考博] 华南理工大学 “新能源交叉创新团队--主动安全”课题组招收海外联合培养博士生 +4 hubble 2024-12-20 5/250 2024-12-21 16:16 by 那片叶落
[教师之家] 40多岁的教师危机感很强啊 (金币+5) +21 farmer1999 2024-12-16 28/1400 2024-12-21 16:07 by willbuilder
[硕博家园] 柔性引进硕博人才 +6 大发财树 2024-12-20 6/300 2024-12-20 15:32 by owlmac
[教师之家] 某老师在组会上总是针对我带的学生怎么办? +15 lice_1987 2024-12-17 23/1150 2024-12-20 04:50 by 胖胖的大海
[基金申请] 除了评职称,除了考评,你对社科项目是什么态度? +9 xinrenxf 2024-12-16 15/750 2024-12-19 11:31 by cwwqt
信息提示
请填处理意见