计算化学公社

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

[综合交流] 从三重态到单重态之间的过渡态要怎么找

[复制链接 Copy URL]

410

帖子

0

威望

2519

eV
积分
2929

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
有一个A-->B的反应,A是三重态,B是单重态,从A到B的过渡态找了很久都没有找到,希望大家能出出招,具体情况如下

如果把过渡态设成三重态,找到的过渡态的IRC起点是A,但终点不是B;
如果把过渡态设成单重态,找到的过渡态的IRC终点是B,但起点不是A。

除此之外,也试过把A设成单重态,或着把B设成三重态,优化后都会散架
把过渡态设成单重态,然后用非限制性方法优化也没有成功

928

帖子

1

威望

8273

eV
积分
9221

Level 6 (一方通行)

2#
发表于 Post on 2023-1-29 18:24:03 | 只看该作者 Only view this author
本帖最后由 hebrewsnabla 于 2023-1-29 20:49 编辑

IRC是同一个态的势能面上的曲线,它不可能连接两个不同的态的结构。不能寻找一个不存在的东西。你得考虑交叉点。

1万

帖子

0

威望

9902

eV
积分
22161

Level 6 (一方通行)

3#
发表于 Post on 2023-1-29 18:25:56 | 只看该作者 Only view this author
这种反应要找MECP
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?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

6万

帖子

99

威望

6万

eV
积分
125245

管理员

公社社长

4#
发表于 Post on 2023-1-29 21:47:26 | 只看该作者 Only view this author
使用sobMECP程序结合Gaussian程序搜索极小能量交叉点
http://sobereva.com/286http://bbs.keinsci.com/thread-865-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

410

帖子

0

威望

2519

eV
积分
2929

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2023-2-11 17:37:25 | 只看该作者 Only view this author
本帖最后由 mfdsrax2 于 2023-2-11 17:39 编辑
wzkchem5 发表于 2023-1-29 18:25
这种反应要找MECP

我看了一下介绍,是把我现有的过渡态作为初始结构,用sobMECP去计算,找到MECP,作为这个反应的“过渡态”的意思吗?

1万

帖子

0

威望

9902

eV
积分
22161

Level 6 (一方通行)

6#
发表于 Post on 2023-2-11 17:47:30 | 只看该作者 Only view this author
mfdsrax2 发表于 2023-2-11 10:37
我看了一下介绍,是把过渡态作为初始结构,用sobMECP去计算,找到MECP,作为这个反应的“过渡态”的意思 ...

初始结构可以选单重态或三重态过渡态,但不是必须这样,如果有更好的MECP初始猜测(手画的也行,扫描出来的也行),也可以用。
可以用sobMECP,但如果你用的量化软件本身就支持找MECP,那也可以用量化软件自带的找MECP的功能。
当且仅当找出来的MECP在整个反应中是能量最高点时,应当用MECP能量来计算能垒(不过此时计算自由能垒可能有一些困难,因为MECP处的自由能没有明确定义,不能直接读频率计算的结果来获得;这个还请其他人补充,有没有比较好的办法能计算MECP是能量最高点的反应的自由能垒)。但是MECP未必是反应途径的能量最高点。例如可能A经过三重态过渡态,一路下坡,在下坡的路上经过MECP变成单重态。那么此时过渡态应该取作三重态过渡态,IRC应该这样做:(1)从过渡态跑反向IRC,跑到A附近;(2)从过渡态跑正向IRC,证明IRC经过MECP附近(未必严格经过MECP,这个是正常的);(3)从MECP跑downhill IRC跑到B。另外也可能A在上坡过程中经过MECP变成单重态,走单重态过渡态得到B,此时的处理方法可以自行类推。
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?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

410

帖子

0

威望

2519

eV
积分
2929

Level 5 (御坂)

7#
 楼主 Author| 发表于 Post on 2023-2-13 10:54:34 | 只看该作者 Only view this author
wzkchem5 发表于 2023-2-11 17:47
初始结构可以选单重态或三重态过渡态,但不是必须这样,如果有更好的MECP初始猜测(手画的也行,扫描出来 ...

非常感谢,高斯能够找MECP吗,或者ORCA?

178

帖子

3

威望

1922

eV
积分
2160

Level 5 (御坂)

8#
发表于 Post on 2023-2-13 11:08:31 | 只看该作者 Only view this author
mfdsrax2 发表于 2023-2-13 10:54
非常感谢,高斯能够找MECP吗,或者ORCA?

Gaussian(在DFT级别)不行,可以用sobMECP调Gaussian,ORCA可以

本版积分规则 Credits rule

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

GMT+8, 2026-2-28 19:20 , Processed in 0.161876 second(s), 20 queries , Gzip On.

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