计算化学公社

 找回密码 Forget password
 注册 Register
Views: 15232|回复 Reply: 12
打印 Print 上一主题 Last thread 下一主题 Next thread

[Multiwfn功能更新] Multiwfn已支持计算任意轨道的pi成份

[复制链接 Copy URL]

6万

帖子

99

威望

5万

eV
积分
120192

管理员

公社社长

2020-Jan-15注:此贴介绍的方法笔者已经发表于Theor. Chem. Acc., 139, 25 (2020) DOI: 10.1007/s00214-019-2541-z,其中专门对pi电子结构分析的算法和应用做了详细、深入的介绍和讨论,给了许多应用例子。使用本帖提到的Multiwfn的功能的读者,除了引用Multiwfn的原文以外,也请同时引用上述文章


之前笔者在下文中专门介绍过Multiwfn的自动识别pi轨道的功能,对于研究pi电子结构非常有用
在Multiwfn中单独考察pi电子结构特征
http://sobereva.com/432http://bbs.keinsci.com/thread-10610-1-1.html

今日在官网http://sobereva.com/multiwfn更新的Multiwfn 3.6(dev)中对这个功能进行了扩展,现在已经可以自动计算任意轨道的pi成份。具体来说,对于一个给定的体系,先用Multiwfn的轨道定域化功能产生定域化轨道(LMO),然后进入pi轨道识别功能判断出pi型LMO,然后用户提供一个含有当前体系的其它轨道的文件(比如含有MO、NTO、自然轨道等都可以),这个文件里的轨道将会向被指认为pi型的LMO进行投影,那些pi型LMO对某个轨道的贡献就是这个轨道的pi成分了。这思路很简单,经过实测发现效果非常理想,下面是一个例子,cycloheptatriene:


可见轨道图形和pi成份对应很好。比如MO23,在一些C-C键轴上,以及好几个氢上都有显著等值面,和pi特征无关,因此MO23的pi成份不算特别大。而MO25肉眼一看就知道几乎都是pi特征的等值面,这是为什么算出来的它的pi成份很大。


此方法的原理在最新版Multiwfn手册3.100.22节有详细介绍,在手册4.100.22节给了计算出上图结果的具体例子。

这个功能用处很大。比如之前笔者介绍过AICD这个方法,见
使用AICD 2.0绘制磁感应电流图
http://sobereva.com/294http://bbs.keinsci.com/thread-1270-1-1.html
对于pi芳香性体系,我们考察AICD的时候往往都是想考察pi轨道产生的环电流,这需要在Gaussian输入文件里给定pi分子轨道编号,然而当体系比较大时,一方面一个一个看图挑pi轨道费劲,另一方面,当体系扭曲比较显著的时候,往往存在显著的sigma与pi的混合,有时候没有个客观的选取pi轨道的标准。而利用Multiwfn的这个功能,就可以定义一个阈值,比如80%,只把pi成份高于这个阈值的分子轨道列出来,然后用于AICD的pi环电流的计算,此时就没有什么含糊性了。

再比如,我们讨论电子激发问题的时候,用Multiwfn可以做NTO分析,使得电子激发往往能够通过一对NTO轨道的跃迁来描述,见
使用Multiwfn做自然跃迁轨道(NTO)分析
http://sobereva.com/377http://bbs.keinsci.com/thread-5980-1-1.html
用上面介绍的功能考察NTO的pi成分的话,我们就可以在文章中给出跃迁涉及的NTO轨道的pi成份,从而讨论跃迁的pi特征。


顺带一提,Multiwfn也可以计算轨道的孤对电子成份,需要先做轨道定域化找到对应于孤对电子的轨道,然后再用Multiwfn的轨道对应关系定量分析功能去考察分子轨道中这个孤对电子轨道的成份,例子看手册4.200.6.2节。

评分 Rate

参与人数
Participants 9
eV +43 收起 理由
Reason
汪文焕 + 5 赞!
lonemen + 3 谢谢分享
ezez + 5 赞!
Novice + 5 赞!
柒月小鱼 + 5 疯狂打call!!!
ntrip + 5 好物!
ggdh + 5 赞!
zsu007 + 5 赞!
三石草祭 + 5 233333

查看全部评分 View all ratings

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

1043

帖子

0

威望

4174

eV
积分
5217

Level 6 (一方通行)

2#
发表于 Post on 2019-5-5 17:07:55 来自手机 | 只看该作者 Only view this author
支持! 不知道能否处理过渡金属与配体之间的pi成键,比如CO的pi*对金属d轨道的反馈pi键?

908

帖子

37

威望

5435

eV
积分
7083

Level 6 (一方通行)

3#
发表于 Post on 2019-5-5 18:49:04 | 只看该作者 Only view this author
请问 这个方法是不是默认被挑选出来的LMO轨道是100%的Π轨道
会不会有LMO自己的的Π含量不纯的情况?

908

帖子

37

威望

5435

eV
积分
7083

Level 6 (一方通行)

4#
发表于 Post on 2019-5-5 18:51:55 | 只看该作者 Only view this author
granvia 发表于 2019-5-5 17:07
支持! 不知道能否处理过渡金属与配体之间的pi成键,比如CO的pi*对金属d轨道的反馈pi键?

你这个问题直接用LMO就可以观察,好像用不到这个功能。。

6万

帖子

99

威望

5万

eV
积分
120192

管理员

公社社长

5#
 楼主 Author| 发表于 Post on 2019-5-6 09:15:28 | 只看该作者 Only view this author
granvia 发表于 2019-5-5 17:07
支持! 不知道能否处理过渡金属与配体之间的pi成键,比如CO的pi*对金属d轨道的反馈pi键?

这是CDA范畴的
使用Multiwfn做电荷分解分析(CDA)、绘制轨道相互作用图
http://sobereva.com/166
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6万

帖子

99

威望

5万

eV
积分
120192

管理员

公社社长

6#
 楼主 Author| 发表于 Post on 2019-5-6 09:16:16 | 只看该作者 Only view this author
ggdh 发表于 2019-5-5 18:49
请问 这个方法是不是默认被挑选出来的LMO轨道是100%的Π轨道
会不会有LMO自己的的Π含量不纯的情况?

扭曲的体系,LMO也不是100%纯的pi,但一般来说可以认为近乎100%(不放心可以肉眼再看一下等值面),所以就把它当做纯的pi来判断其它类型轨道了
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

1043

帖子

0

威望

4174

eV
积分
5217

Level 6 (一方通行)

7#
发表于 Post on 2019-5-6 18:22:17 来自手机 | 只看该作者 Only view this author
sobereva 发表于 2019-5-6 09:15
这是CDA范畴的
使用Multiwfn做电荷分解分析(CDA)、绘制轨道相互作用图
http://sobereva.com/166

不是这个意思。我的意思是,金属配合物的分子轨道中能否实现自动sigma, pi(甚至delta)对称性的分离

6万

帖子

99

威望

5万

eV
积分
120192

管理员

公社社长

8#
 楼主 Author| 发表于 Post on 2019-5-7 08:14:53 | 只看该作者 Only view this author
granvia 发表于 2019-5-6 18:22
不是这个意思。我的意思是,金属配合物的分子轨道中能否实现自动sigma, pi(甚至delta)对称性的分离

delta不行
这个算法自动识别配合物的pi分子轨道应该没问题(我没测试过)
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

2

帖子

0

威望

1361

eV
积分
1363

Level 4 (黑子)

9#
发表于 Post on 2019-5-14 11:25:20 | 只看该作者 Only view this author
Sobereva 老师,请问在获取体系的LMO后,如何得到pi成分在MO中的百分比。
我目前在手册4.100.22中没有看到这一部分的描述。

6万

帖子

99

威望

5万

eV
积分
120192

管理员

公社社长

10#
 楼主 Author| 发表于 Post on 2019-5-14 20:01:03 | 只看该作者 Only view this author
张含宇 发表于 2019-5-14 11:25
Sobereva 老师,请问在获取体系的LMO后,如何得到pi成分在MO中的百分比。
我目前在手册4.100.22中没有看到 ...

手册太老,看最新的。怎么计算手册里写明了。注意看清楚程序和手册的发布日期
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

2

帖子

0

威望

1361

eV
积分
1363

Level 4 (黑子)

11#
发表于 Post on 2019-5-15 20:57:29 | 只看该作者 Only view this author
sobereva 发表于 2019-5-14 20:01
手册太老,看最新的。怎么计算手册里写明了。注意看清楚程序和手册的发布日期

找到了最新的手册,谢谢您

78

帖子

0

威望

1446

eV
积分
1524

Level 5 (御坂)

12#
发表于 Post on 2023-5-22 19:21:32 | 只看该作者 Only view this author
这个是非平面研究pi轨道成分,那对平面分子也一样可以这样分析?还有pi轨道成分得到,那剩下的是sigma成分吗?

6万

帖子

99

威望

5万

eV
积分
120192

管理员

公社社长

13#
 楼主 Author| 发表于 Post on 2023-5-23 02:04:03 | 只看该作者 Only view this author
wzy 发表于 2023-5-22 19:21
这个是非平面研究pi轨道成分,那对平面分子也一样可以这样分析?还有pi轨道成分得到,那剩下的是sigma成分 ...

在Multiwfn中单独考察pi电子结构特征
http://sobereva.com/432http://bbs.keinsci.com/thread-10610-1-1.html

里面不就有平面体系的分析例子
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-18 01:52 , Processed in 0.227091 second(s), 30 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list