24小时热门版块排行榜    

查看: 8268  |  回复: 6

zheyi-1984

铜虫 (小有名气)

[求助] 高斯中的分子轨道MO,NBO计算,以及电子云密度,电子云密度差计算的一些疑问

如题,最近想通过对分子轨道的分析和电子云密度变化解释一些问题,LZ是做实验的量化小白,因为基础薄弱只好厚颜发这个伸手贴,希望大家能多多给我一些建议。

刚开始自学高斯,概念不清,所以有些表述很可能不太正确,基本想解决的问题如下:

想计算二价分子离子,如:CH2ClI的二价离子中的C-Cl键和C-I键分别断裂的两种解离过程中,随着化学键C-X键的键长变化,CH2ClI的二价分子离子轨道的变化以及电子云密度分布的变化,并以此为依据对C-Cl键和C-I键断裂的难易程度作出一些解释。

以下是我想要学习计算的东西:
1.Cl和I原子对分子离子轨道的贡献值,其随着对应的C-CL键和C-I键键长的伸长的变化。
目前通过在论坛上搜索学习,学会了计算MO和NBO的方法,但是对这两种不同的分子轨道计算方法的实际物理意义很迷糊,如若我计算的初衷是想要比较Cl原子和I原子的原子轨道各自对CH2ClI的二价分子离子轨道的贡献值,那么使用哪种方法比较合适呢?我现在的理解觉得貌似是用NBO好一些?还有一个很白丁的问题,原子对于分子homo和lomo轨道的贡献值大小,能说明什么物理意义呢?

2.Cl和I原子附近的电子云分布,及其随着对应的C-CL键和C-I键键长的伸长的变化。
貌似这点可以直接通过用高斯产生波函数,然后用multiwfn软件直接得到分子离子的整体电子云密度图,也可以自由做差值(尝试了发现挺方便快捷的)。我现在的问题是,如果我想得到homo,lomo各自对应的电子云密度图,以及差值,该如何是好呢?

3.在利用高斯以及multiwfn进行计算的过程中,我还有一些疑问一直不明:

a. 由于通过高斯获取的分子轨道信息是蕴含在chk文件中的,那么是否必须通过对分子结构的优化才能得到chk文件?如果是这样,那么我遇到的一个问题就是,我想得到不同键长下的分子离子轨道变化趋势,可是并非每一个键长都能opt得到一个对应的稳定构型呀求破或者科普。

b.利用高斯获取用于multiwfn的波函数时,可以不优化,直接只让output=wfn吗?这样做可以吗?

c.物理专业的我对于各种轨道的理解真的差不多等于白丁,找了论坛推荐的一些帖子和书,看了也不是很明白,大家能否给我科普下homo.lomo,以及他们各自的-1,-2有什么具体物理意义和应用,或者推荐一些比较好接受的书和教程,我去学习

d.NBO计算得到一堆东西,但是具体这些数据如何使用一头雾水,求如何入门呢?

一口气提了好多问题...希望路过的大神们能给些建议和帮助。LZ很穷也很小白,能给的金币数量微薄,但有一颗谦虚学习的心,请大家多多关照。
回复此楼

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

气态分子化学动力学 HOMO与LUMO

» 猜你喜欢

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

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

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

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫


【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
gkf
gkf高: 金币+5, 这么晚,讲解这么多,不容易啊!! 2013-06-26 02:20:05
zhou2009: 编辑内容 2013-06-26 06:52
zheyi-1984: 金币+20, ★★★★★最佳答案 2013-07-05 20:04:21
1 记住Multiwfn,既可以根据NBO输出的信息计算自然原子轨道(NAO)对分子轨道的贡献,也可以计算各个基函数对分子轨道的贡献并加和成原子轨道对分子轨道的贡献。一般建议用前者。这本身没有确切的物理意义,取决于你怎么讨论。

2 进入主功能6的选项26,把HOMO或LUMO外的其余轨道占据数都设为0(即只保留HOMO或LUMO轨道的占据数)。然后退回到主菜单并通过主功能5来计算电子密度,得到的就是HOMO或LUMO的电子密度。分别计算完HOMO和LUMO的密度并导出为cube文件后,再利用Multiwfn的主功能13里的功能11对两个cube文件进行相减即是差值。

3a 不需要优化,除非你研究的是分子的稳定结构。

3b 可以

3c 建议你先看看结构化学书

3d 看NBO官方主页的教程,或综述文章如WIREs Comput Mol Sci 2012, 2: 1–42 ,或相关书籍如Discovering Chemistry with Natural Bond Orbitals(Weinhold,Landis,2012)

[ Last edited by zhou2009 on 2013-6-26 at 06:52 ]
2楼2013-06-26 00:44:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zheyi-1984

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by sobereva at 2013-06-26 00:44:04
1 记住Multiwfn,既可以根据NBO输出的信息计算自然原子轨道(NAO)对分子轨道的贡献,也可以计算各个基函数对分子轨道的贡献并加和成原子轨道对分子轨道的贡献。一般建议用前者。这本身没有确切的物理意义,取决于你怎 ...

多谢发现sobereva老师一直很勤恳,大半夜的应助,非常感激。
3楼2013-06-26 08:53:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zheyi-1984

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by sobereva at 2013-06-26 00:44:04
1 记住Multiwfn,既可以根据NBO输出的信息计算自然原子轨道(NAO)对分子轨道的贡献,也可以计算各个基函数对分子轨道的贡献并加和成原子轨道对分子轨道的贡献。一般建议用前者。这本身没有确切的物理意义,取决于你怎 ...

sobereva老师你好,我尝试着用multiwfn做了一些计算:

针对我提问中的2,我按照你的建议将fchk文件用multiwfn打开后,选择了6,然后出现了很多轨道,共计64种之多,截图见附件。

multiwfn很贴心,指出了15是homo,16是lomo,我想知道的是,截图中给出的每个轨道信息那一行的各种数字是什么意思呢?

还有图中的蓝色和绿色是什么意思呢?表示不同相位吗?
高斯中的分子轨道MO,NBO计算,以及电子云密度,电子云密度差计算的一些疑问
1.jpg


高斯中的分子轨道MO,NBO计算,以及电子云密度,电子云密度差计算的一些疑问-1
2.jpg

4楼2013-06-26 15:44:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zheyi-1984

铜虫 (小有名气)

还有一个问题,MO和NBO得到的轨道分布显然不同,就我这个问题的话,选择哪一种方法比较好呢?
5楼2013-06-26 15:46:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zheyi-1984

铜虫 (小有名气)

对于图1中各轨道的参数,我猜想是au和eV为单位的轨道对应能量,还有轨道类型,我只是不明白的,能量均为0的60~64轨道是什么思路?
6楼2013-06-26 15:55:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sobereva

至尊木虫 (著名写手)

本人已永久离开小木虫


【答案】应助回帖

★ ★ ★
zhou2009: 金币+3 2013-06-26 19:12:41
引用回帖:
6楼: Originally posted by zheyi-1984 at 2013-06-26 15:55:07
对于图1中各轨道的参数,我猜想是au和eV为单位的轨道对应能量,还有轨道类型,我只是不明白的,能量均为0的60~64轨道是什么思路?

你的理解正确。那些能量为0的不用管。(那是由于高斯自动消除了线性依赖的基函数,造成实际MO数少于基函数数目。但是Multiwfn强行要求MO数等于基函数数目,所以最后一些轨道就人为地填补上了。这不会对计算结果有任何影响)

用MO还是NBO来讨论都可以,你最好先多了解NBO相关的知识,弄清楚NBO的意义,然后来决定怎么讨论。

绿色代表正相位,蓝色代表负相位
7楼2013-06-26 16:04:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zheyi-1984 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:29 by ZPa0EcMwuECS
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:58 by ZPa0EcMwuECS
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:49 by ZPa0EcMwuECS
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:47 by ZPa0EcMwuECS
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:36 by ZPa0EcMwuECS
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-30 00:59 by winsaint
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 00:52 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 12/600 2026-08-30 00:31 by ZPa0EcMwuECS
[基金申请] 中青基了要发朋友圈吗? +4 349506619 2026-08-28 4/200 2026-08-29 22:41 by alongwaytogo
[基金申请] 面上意见出来了 +8 黄鸟于飞Chao 2026-08-29 16/800 2026-08-29 20:45 by 黄鸟于飞Chao
[基金申请] 为什么到现在没收到通知? +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +9 苏东坡二世 2026-08-23 9/450 2026-08-29 14:41 by hustersqt
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 9/450 2026-08-29 09:59 by jklily
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见