| 查看: 7101 | 回复: 6 | ||||
[求助]
高斯中的分子轨道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 |
» 猜你喜欢
北京师范大学胡德平课题组以及袁萌伟课题组考研招生
已经有0人回复
-
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有132人回复
-
祈福论文小修
已经有21人回复
-
香港理工大学计算材料与物理课题组招聘博后
已经有15人回复
-
香港理工大学计算材料与物理课题组招聘博后
已经有15人回复
-
西班牙Institut Català de Nanociència i Nanotecnologia, CSIC-BIST招收公派博士生
已经有0人回复
-
请问四氢呋喃溶解的聚合物用甲醇沉淀时，如何使沉淀过程加速？
已经有2人回复
-
英国利物浦大学屠昕和雍雪课题组招收CSC博士生
已经有26人回复
-
北京交通大学王熙教授团队招收2026年入学催化方向博士研究生（还有1个名额）
已经有0人回复
-
华北电力大学（北京）（第一性原理计算）博士招生——学博，专博各1人
已经有45人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
怎么解读homo和lomo分子轨道图？
已经有11人回复- 【求助】NBO轨道指的是nbo电荷的分布情况吗？
已经有3人回复
- Gaussian计算得到的前线分子轨道图能说明什么？
已经有7人回复
- 高斯计算中离子电荷的处理
已经有15人回复
- 大家做NBO分析都用其中的什么内容呢？总觉得自己分析不够透彻...
已经有21人回复
- 新手，请高手帮忙，计算小分子中某个碳原子上的电子云密度
已经有15人回复
- 【求助】如何利用高斯计算电子云密度
已经有7人回复
- 【求助】密立根电荷和NBO电荷的差异
已经有4人回复
- 【小卒个人文集】写给初学者:小议形式电荷对量化初学者的误导
已经有27人回复
sobereva
至尊木虫 (著名写手)
本人已永久离开小木虫
- QC强帖: 17
- 应助: 226 (大学生)
- 贵宾: 5.4
- 金币: 2180.5
- 红花: 217
- 帖子: 1241
- 在线: 366.8小时
- 虫号: 281271
- 注册: 2006-09-23
- 专业: 宇宙学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
gkf
gkf高: 金币+5, 这么晚，讲解这么多，不容易啊！！ 2013-06-26 02:20:05
zhou2009: 编辑内容 2013-06-26 06:52
zheyi-1984: 金币+20, ★★★★★最佳答案 2013-07-05 20:04:21
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
发现sobereva老师一直很勤恳，大半夜的应助，非常感激。3楼2013-06-26 08:53:27
4楼2013-06-26 15:44:36
5楼2013-06-26 15:46:06
6楼2013-06-26 15:55:07
sobereva
至尊木虫 (著名写手)
本人已永久离开小木虫
- QC强帖: 17
- 应助: 226 (大学生)
- 贵宾: 5.4
- 金币: 2180.5
- 红花: 217
- 帖子: 1241
- 在线: 366.8小时
- 虫号: 281271
- 注册: 2006-09-23
- 专业: 宇宙学
7楼2013-06-26 16:04:47