24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 7866  |  回复: 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很穷也很小白,能给的金币数量微薄,但有一颗谦虚学习的心,请大家多多关照。
回复此楼
已阅   回复此楼   关注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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 22408 312求调剂 +14 门路摸摸 2026-04-14 14/700 2026-04-14 16:29 by zhouxiaoyu
[考研] 085400电子信息类(川大控制工程)求调剂可跨专业 求老师联系 +4 626776879 2026-04-08 4/200 2026-04-14 14:49 by htcscq
[考研] 生物学调剂 +7 纸扇zhishan 2026-04-13 7/350 2026-04-14 14:21 by jyl0317
[考研] 08工学 309分求调剂 +11 Yin DY 2026-04-08 11/550 2026-04-14 12:55 by fqwang
[考研] 一志愿华南理工大学331分材料求调剂 +10 天下ww 2026-04-09 11/550 2026-04-13 23:25 by pies112
[考研] 272分材料子求调剂 +41 Loy0361 2026-04-10 53/2650 2026-04-13 14:20 by 张zhihao
[考研] 一志愿华工085600 331分 +7 天下ww 2026-04-09 7/350 2026-04-13 09:01 by lhj2009
[考研] 0831一轮调剂失败求助 +10 小熊睿睿_s 2026-04-11 10/500 2026-04-12 22:43 by 长弓傲
[教师之家] 山东双非院校考核超级无底线,领导幸灾乐祸,教师遭殃恐 +3 qut2026 2026-04-11 7/350 2026-04-12 20:24 by qut2026
[考研] 求调剂 +6 archer.. 2026-04-09 8/400 2026-04-11 10:55 by zhq0425
[考研] 农学0904 312求调剂 +6 Say Never 2026-04-10 6/300 2026-04-11 10:33 by wwj2530616
[考研] 284求调剂 +12 archer.. 2026-04-10 13/650 2026-04-11 08:44 by zhq0425
[考研] 计算机类求调剂,22408-274分 +7 上岸de小虫 2026-04-09 8/400 2026-04-10 19:56 by fxue1114
[考研] 已调剂 +18 柴郡猫_ 2026-04-09 19/950 2026-04-09 22:10 by 柴郡猫_
[考研] 278求调剂 +27 范婷娜 2026-04-07 31/1550 2026-04-09 20:49 by zhouxiaoyu
[考研] 367求调剂 +10 hffQAQ 2026-04-09 10/500 2026-04-09 18:06 by lijunpoly
[考研] 311求调剂 +6 surte 2026-04-08 13/650 2026-04-09 14:00 by surte
[考研] 334求调剂 +16 Riot2025 2026-04-08 17/850 2026-04-09 09:28 by wdyheheeh
[考研] 生物学328分求调剂 +9 闪电kkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
[考研] 电子信息346 +4 zuoshaodian 2026-04-08 4/200 2026-04-08 11:54 by zzucheup
信息提示
请填处理意见