计算化学公社

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

[新手求助] 用高斯计算时如何判断一个多重态过渡态是自由基偶联还是电子对进攻的极性机理

[复制链接 Copy URL]

15

帖子

1

威望

447

eV
积分
482

Level 3 能力者

本帖最后由 czhang8 于 2025-4-30 17:30 编辑

如题,本菜鸟现在的课题中有一步C-C键的生成,已知能量最低的过渡态是三重态。而反应物可能为:
1. 三重态亲电试剂与单重态亲核试剂 电子对进攻

2. 1中反应物经单电子转移后生成的 四重态 与 二重态 自由基偶联 实际有四个单电子 但是有两个自旋相反 即将成对

请问各位大佬有什么办法知道这个三重态是哪种构成的吗?看TS的 Mulliken spin population?根据我的理解,似乎在过渡态中,即使对二重态反应物,单电子也并不会分布在除成键位置以外的地方,可能并不能分辨。
抑或是这两种情况的过渡态 我们认为是同一个?

多谢解答!

1万

帖子

0

威望

8942

eV
积分
20684

Level 6 (一方通行)

2#
发表于 Post on 2025-5-1 11:37:27 | 只看该作者 Only view this author
看population可以大体看出来,但不要用Mulliken population,Mulliken对基组太敏感了,当基组比较大时没有物理意义。可以用Hirshfeld或NPA population

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
czhang8 + 5 赞!

查看全部评分 View all ratings

Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

15

帖子

1

威望

447

eV
积分
482

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2025-5-1 11:51:27 | 只看该作者 Only view this author
本帖最后由 czhang8 于 2025-4-30 19:55 编辑
wzkchem5 发表于 2025-4-30 19:37
看population可以大体看出来,但不要用Mulliken population,Mulliken对基组太敏感了,当基组比较大时没有 ...

感谢大佬!请问population去分辨这两种情况可以关注哪些区别呢?成键原子吗?

另外想请问一下 高斯opt任务结果中的mulliken spin population可以作为一个非常初级的参考 来看OSS / CSS吗?还是说这个也完全不可靠呢?

1万

帖子

0

威望

8942

eV
积分
20684

Level 6 (一方通行)

4#
发表于 Post on 2025-5-1 12:04:31 | 只看该作者 Only view this author
czhang8 发表于 2025-5-1 11:51
感谢大佬!请问population去分辨这两种情况可以关注哪些区别呢?成键原子吗?

另外想请问一下 高斯opt ...

最简单的就是看亲核试剂的自旋密度,如果远小于1,就是机理(1);如果接近1,就是机理(2);如果既不是很接近0又不是很接近1,就看IRC上的其他结构的自旋密度

评分 Rate

参与人数
Participants 1
eV +3 收起 理由
Reason
czhang8 + 3 谢谢

查看全部评分 View all ratings

Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

15

帖子

1

威望

447

eV
积分
482

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2025-5-1 12:19:33 | 只看该作者 Only view this author
wzkchem5 发表于 2025-4-30 20:04
最简单的就是看亲核试剂的自旋密度,如果远小于1,就是机理(1);如果接近1,就是机理(2);如果既不是 ...

谢谢您!您说mulliken spin 特敏感不推荐,那opt任务给出来的mulliken spin可以做oss/css的初步判断吗?

6万

帖子

99

威望

5万

eV
积分
120060

管理员

公社社长

6#
发表于 Post on 2025-5-1 12:52:58 | 只看该作者 Only view this author
czhang8 发表于 2025-5-1 12:19
谢谢您!您说mulliken spin 特敏感不推荐,那opt任务给出来的mulliken spin可以做oss/css的初步判断吗?

没弥散函数时Mulliken自旋布居还能凑合用,做个定性判断;而带弥散函数时绝对不行。此时建议用Multiwfn计算Hirshfeld自旋布居,参考下文
谈谈自旋密度、自旋布居以及在Multiwfn中的绘制和计算
http://sobereva.com/353http://bbs.keinsci.com/thread-4422-1-1.html

Mulliken布居分析对基组的极度敏感性在下文的图5展现得十分充分
《原子电荷计算方法的对比》
http://www.whxb.pku.edu.cn/CN/abstract/abstract27818.shtml
下文的我的原子电荷/布居分析综述文章也提到了这点
一篇深入浅出、完整全面介绍原子电荷的综述文章已发表!
http://sobereva.com/714http://bbs.keinsci.com/thread-46067-1-1.html

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
czhang8 + 2 赞!

查看全部评分 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

15

帖子

1

威望

447

eV
积分
482

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2025-7-8 12:59:05 | 只看该作者 Only view this author
wzkchem5 发表于 2025-4-30 19:37
看population可以大体看出来,但不要用Mulliken population,Mulliken对基组太敏感了,当基组比较大时没有 ...

大佬您好!我遵循您的指导进行了尝试,排除了单电子氧化还原的可能。现在我需要单独比较一个偶联过程是两边各出一个电子的自由基偶联,还是极性机理的亲电试剂-亲核试剂的偶联。请问还可以通过Hirshfeld自旋布局判断吗?我发现成键位点上自旋布居数有,但比较小,反而是成键位点alpha位置的C(原本和反应位点形成双键)自旋密度比较显著(但是也不到0.5)。这种可以判定为单电子偶联吗?
我就是很困惑如果是单电子机理,为什么过渡态结构中成键位点自旋密度只有0.0x?我原本理解会有接近1;而如果是极性机理,是不是不应该有自旋密度在alphaC上?

1万

帖子

0

威望

8942

eV
积分
20684

Level 6 (一方通行)

8#
发表于 Post on 2025-7-8 13:27:57 | 只看该作者 Only view this author
czhang8 发表于 2025-7-8 12:59
大佬您好!我遵循您的指导进行了尝试,排除了单电子氧化还原的可能。现在我需要单独比较一个偶联过程是两 ...

可以。在IRC上各个点处算Hirshfeld布居,可以知道电子是否有转移,以及何时转移。最好不要只看一个结构(如过渡态)
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

15

帖子

1

威望

447

eV
积分
482

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2025-7-9 05:33:36 | 只看该作者 Only view this author
wzkchem5 发表于 2025-7-7 21:27
可以。在IRC上各个点处算Hirshfeld布居,可以知道电子是否有转移,以及何时转移。最好不要只看一个结构( ...

谢谢Wang老师!就是取IRC的结构 算一个单点能生成波函数文件 然后分析对吗?生成IRC结构的波函数文件时需不需要读取什么别的波函数文件(比如该过渡态的chk文件)呢?(我担心只从结构出发去算波函数可能会让单电子分布等性质差异比较大)

1万

帖子

0

威望

8942

eV
积分
20684

Level 6 (一方通行)

10#
发表于 Post on 2025-7-9 12:40:59 | 只看该作者 Only view this author
czhang8 发表于 2025-7-9 05:33
谢谢Wang老师!就是取IRC的结构 算一个单点能生成波函数文件 然后分析对吗?生成IRC结构的波函数文件时需 ...

对,对于IRC上的每个结构都算波函数。
IRC计算过程中本来就会生成每个结构的波函数,只是有的软件不保存下来而已。如果要重新生成波函数,最好读过渡态的波函数做初猜

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
czhang8 + 5 谢谢

查看全部评分 View all ratings

Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

本版积分规则 Credits rule

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

GMT+8, 2025-8-12 15:30 , Processed in 0.368827 second(s), 21 queries , Gzip On.

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