计算化学公社

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

[Multiwfn使用咨询] 轨道定域化后出现香蕉键怎么处理

[复制链接 Copy URL]

35

帖子

0

威望

345

eV
积分
380

Level 3 能力者

本帖最后由 rorrim 于 2026-6-22 13:44 编辑

笔者多次尝试后发现,在默认设置下,Multiwfn的自动检测π轨道的功能基本能够正确判断π轨道,但有时候会多判断π轨道有时候会少判断π轨道,导致LOL-π上一些位点等值面明显有问题,所以笔者会自己筛查一遍LMO看自动判断的π轨道是否正确,

一般根据等值面形状就能够判断是σ轨道还是π轨道,但笔者无法确定图中的氯富勒烯上的两个LMO是σ轨道还是π轨道,这两个LMO跟富勒烯上通常的σ轨道不太一样,如果调低等值面则与π轨道相似。
理论上这两个位置应该有π电子,如果不选择这两个LMO,LOL-π上这两个位置就没有等值面了


笔者之后意识到这可能是香蕉键,找到了同一个碳碳键上对应的两个香蕉键,请问在选择π轨道时该如何处理香蕉键,因为之后要计算LOL-π的等值面
以及为什么其余轨道都正常对应σ轨道与π轨道,这个区域确却是香蕉键



000118.png (596.77 KB, 下载次数 Times of downloads: 0)

000118.png

000226.png (620.29 KB, 下载次数 Times of downloads: 0)

000226.png

6万

帖子

99

威望

6万

eV
积分
128183

管理员

公社社长

2#
发表于 Post on 2026-6-21 02:24:11 | 只看该作者 Only view this author
直接贴图,不要把图像文件放压缩包里,已经把你的违规附件删除
认真看置顶的社员必读贴了解怎么正确贴图
北京科音自然科学研究中心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

4365

帖子

4

威望

1万

eV
积分
14490

Level 6 (一方通行)

MOKIT开发者

3#
发表于 Post on 2026-6-22 20:24:55 | 只看该作者 Only view this author
说轨道局域化,一定要同时交待局域化方法细节,例如常见三种局域化方法Boys/PM/ER你用的是哪一种。参与轨道局域化的轨道范围是什么,例如RKS双占据轨道第30-50号,UKS alpha自旋占据轨道第30-50号,RKS空轨道第51-70号,等等。
自动做多参考态计算的程序MOKIT

35

帖子

0

威望

345

eV
积分
380

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2026-6-22 21:46:48 | 只看该作者 Only view this author
zjxitcc 发表于 2026-6-22 20:24
说轨道局域化,一定要同时交待局域化方法细节,例如常见三种局域化方法Boys/PM/ER你用的是哪一种。参与轨道 ...

用的是基于Mulliken布居的Pipek-Mezey定域化,仅定域化了占据轨道,这个分子是闭壳层单重态且波函数稳定。除了我图上的那两个具有香蕉键的CC键以及与之对称的那个键的LMO,其余的LMO都是σ轨道或π轨道的形状

6万

帖子

99

威望

6万

eV
积分
128183

管理员

公社社长

5#
发表于 Post on 2026-6-22 21:52:58 | 只看该作者 Only view this author
用Multiwfn里的Pipek-Mezey结合Becke布居做轨道定域化再试,很可能就得到sigma+pi分离的LMO了
北京科音自然科学研究中心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

1

帖子

0

威望

579

eV
积分
580

Level 4 (黑子)

6#
发表于 Post on 2026-6-22 23:14:23 | 只看该作者 Only view this author
本帖最后由 131A 于 2026-6-22 23:24 编辑

曾经在尝试绘制Möbius型[8]Annulene的pi ELF时遇到过类似的问题。当对所有占据轨道进行定域化时,用Multiwfn默认的基于Mulliken布居的Pipek-Mezey定域化,得到反式双键处的pi轨道具有较强的香蕉键特征,无法被自动识别为pi-LMO。而当改为基于Lowdin布居或Becke布居的Pipek-Mezey定域化,该反式双键处则可以得到pi特征明显的、可以被Multiwfn自动识别出的LMO。

Pi_LMO.png (545.41 KB, 下载次数 Times of downloads: 0)

Pi_LMO.png

35

帖子

0

威望

345

eV
积分
380

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2026-6-23 17:15:50 | 只看该作者 Only view this author
sobereva 发表于 2026-6-22 21:52
用Multiwfn里的Pipek-Mezey结合Becke布居做轨道定域化再试,很可能就得到sigma+pi分离的LMO了

非常感谢,我这个分子的原子数有些多,用Becke布居做轨道定域化在集群上算了一个多小时都没算完,中途网络断开直接白算了,但是改成Lowdin布居后原本有香蕉键的位置就不会出现香蕉键了

35

帖子

0

威望

345

eV
积分
380

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2026-6-23 17:19:47 | 只看该作者 Only view this author
131A 发表于 2026-6-22 23:14
曾经在尝试绘制Möbius型[8]Annulene的pi ELF时遇到过类似的问题。当对所有占据轨道进行定域化时,用Mu ...

非常感谢,改成Lowdin布居后原本有香蕉键的位置就不会出现香蕉键了,但是Multiwfn的自动识别π轨道依旧有问题,还是会多识别并少识别一些π轨道,需要手动判断π轨道,可能是我算的这个分子有些复杂

本版积分规则 Credits rule

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

GMT+8, 2026-7-10 15:36 , Processed in 0.724996 second(s), 24 queries , Gzip On.

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