计算化学公社

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

[综合交流] 求交流:构建过渡态的结构初猜有何经验、技巧

  [复制链接 Copy URL]

3097

帖子

29

威望

1万

eV
积分
17221

Level 6 (一方通行)

本帖最后由 liyuanhe211 于 2015-11-30 00:24 编辑

尤指一些经验性的结论和奇技淫巧。

抛个砖,有其他组的师兄曾说,一般有机反应涉及C-C键断裂的,可以将断裂键键长设置为ground state键长的1.4倍左右,有时效果较好。

292

帖子

8

威望

1696

eV
积分
2148

Level 5 (御坂)

2#
发表于 Post on 2015-11-30 01:35:54 | 只看该作者 Only view this author
可以参考这篇博文
<谈谈搜索反应过渡态的实用计策>
http://www.bccg.org:8888/wordpress/2015/06/18/reactiontips/

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
sobereva + 2

查看全部评分 View all ratings

52

帖子

0

威望

315

eV
积分
367

Level 3 能力者

3#
发表于 Post on 2015-11-30 01:42:33 | 只看该作者 Only view this author
听老师说的,放热反应过渡态结构与反应物更相似,可以在反应物基础的把将要断裂的键的键长延长个20~30%,放热反应则相反。

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
sobereva + 2

查看全部评分 View all ratings

127

帖子

1

威望

1231

eV
积分
1378

Level 4 (黑子)

4#
发表于 Post on 2015-11-30 01:53:03 | 只看该作者 Only view this author
有一个名叫Hammond's postulate的观点,也就是3楼提到的方法。 https://en.wikipedia.org/wiki/Hammond%27s_postulate
一般就是反应物结构稍微改一下作为初猜。
找比较难找的 TS 的话,先是随便把要断掉的键拉长,频率计算, 找到对应TS过程的频率,然后IOP(1/130)选择对应的那个频率 (因为有些不好找的过渡态 可能 初猜的 frequency 的最小mod对应的不是过渡态,所以需要iop一下)。

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
sobereva + 2

查看全部评分 View all ratings

292

帖子

8

威望

1696

eV
积分
2148

Level 5 (御坂)

5#
发表于 Post on 2015-11-30 01:53:30 | 只看该作者 Only view this author
minagami 发表于 2015-11-30 01:42
听老师说的,放热反应过渡态结构与反应物更相似,可以在反应物基础的把将要断裂的键的键长延长个20~30%,放 ...

反应过渡态结构与反应物或生成物相似的说法源于
Hammond假说 :如果化学反应中过渡态和反应物(或生成物)能量近似则结构近似。
但实际计算中,我感觉几乎没用。

55

帖子

0

威望

1201

eV
积分
1256

Level 4 (黑子)

6#
发表于 Post on 2015-11-30 05:38:20 | 只看该作者 Only view this author
我常用的是从初始反应物或者中间体沿着反应坐标扫描,一般会得到很好的过渡态初始结构,优化很容易达到收敛标准。

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
一颗赛艇 + 2 我很赞同

查看全部评分 View all ratings

407

帖子

0

威望

856

eV
积分
1263

Level 4 (黑子)

7#
发表于 Post on 2015-11-30 08:51:47 | 只看该作者 Only view this author
shx 发表于 2015-11-30 05:38
我常用的是从初始反应物或者中间体沿着反应坐标扫描,一般会得到很好的过渡态初始结构,优化很容易达到收敛 ...

请问您说的反应坐标是如何确定的呢?

292

帖子

8

威望

1696

eV
积分
2148

Level 5 (御坂)

8#
发表于 Post on 2015-11-30 10:27:39 | 只看该作者 Only view this author
小书童 发表于 2015-11-30 08:51
请问您说的反应坐标是如何确定的呢?

反应坐标可以理解为是反应过程中每一步里每个原子三个分量上笛卡尔坐标的偏移形成的一个新的坐标轴。其原点是过渡态的结构。楼上说的新键和旧键的断裂实际上是对反应坐标贡献最大的,因为那个时候其他原子几乎不动,主要是参与成键和断键的原子在动。
准确的反应坐标的单位是amu^1/2*Bohr,从量纲上看是距离乘以了质量的平方,意味着反应坐标进行了质量加权。假如我们对反应坐标求一阶导数,并做归一化,这个时候获得量就是反应路径方向向量。

评分 Rate

参与人数
Participants 2
eV +5 收起 理由
Reason
小书童 + 3
卡开发发 + 2

查看全部评分 View all ratings

3097

帖子

29

威望

1万

eV
积分
17221

Level 6 (一方通行)

9#
 楼主 Author| 发表于 Post on 2015-11-30 11:04:51 | 只看该作者 Only view this author
smutao 发表于 2015-11-30 01:53
反应过渡态结构与反应物或生成物相似的说法源于
Hammond假说 :如果化学反应中过渡态和反应物(或生成物 ...

感觉跑IRC出来经常违背Hammond~~

106

帖子

0

威望

4110

eV
积分
4216

Level 6 (一方通行)

10#
发表于 Post on 2015-11-30 22:14:52 | 只看该作者 Only view this author
请问怎么快速寻找二面角旋转而产生的过渡态呢?

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

11#
发表于 Post on 2015-11-30 22:59:36 | 只看该作者 Only view this author
幻七熏 发表于 2015-11-30 22:14
请问怎么快速寻找二面角旋转而产生的过渡态呢?


这个凭直觉一般是比较好猜的。也可以用半经验快速做个柔性扫描,取最高点当初猜。

评分 Rate

参与人数
Participants 2
eV +4 收起 理由
Reason
幻七熏 + 2 谢谢
chrischen1128 + 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

54

帖子

0

威望

591

eV
积分
645

Level 4 (黑子)

12#
发表于 Post on 2015-12-1 16:53:13 | 只看该作者 Only view this author
寻找过渡态就好像在做艺术工作。

164

帖子

0

威望

1501

eV
积分
1665

Level 5 (御坂)

Miss X

13#
发表于 Post on 2015-12-1 17:13:40 | 只看该作者 Only view this author
为嘛小的感觉找过渡态都是在耗费生命
真诚相待,淡淡相处,方天久地长^^

198

帖子

0

威望

2242

eV
积分
2440

Level 5 (御坂)

14#
发表于 Post on 2015-12-1 22:15:18 | 只看该作者 Only view this author
本帖最后由 kevin 于 2015-12-3 21:14 编辑
Shannon 发表于 2015-11-30 01:53
有一个名叫Hammond's postulate的观点,也就是3楼提到的方法。 https://en.wikipedia.org/wiki/Hammond%27s ...

请问什么叫“对应的那个虚频”这种不稳定点计算频率可能会得到很多虚频,怎么判断那个虚频是需要的,根据数值大小么?找到了这个虚频之后怎么优化过渡态呢?您说的整的iop(1/30)意思是保持这个虚频不变然后在这基础上对初猜的结构进行opt+freq么?谢谢啦
%mem=400MW
%nprocshared=4
#p opt=modredundant HF/STO-3G

Title Card Required

0 1
C                 -0.49913700    1.59776300   -0.17442300
C                 -1.95159000   -0.27732100   -0.24923300
C                 -0.88680000   -1.18545000    0.38371500
C                  0.52844700   -0.66681700    0.09069500
C                  0.64434600    0.80883300    0.47024700
H                 -0.40908300    1.55992900   -1.27533000
H                 -1.02881800   -1.18672900    1.47583100
H                  0.71844900   -0.77176000   -0.98841000
H                 -1.84073600   -0.31394700   -1.34646200
O                 -1.74849900    1.06386300    0.21331400
O                 -1.04805300   -2.49716300   -0.13659300
H                 -0.29937200   -3.02289400    0.18294200
O                  1.90402400    1.32140000    0.00608300
H                  1.94606800    2.26440700    0.22471500
C                 -3.37299700   -0.68985100    0.12383500
H                 -3.45939900   -1.77815900    0.02844400
H                 -3.55475600   -0.41498600    1.16722100
O                 -0.39874400    2.91349500    0.28975500
H                 -1.12052300    3.43698000   -0.08619900
O                 -4.36593900   -0.01591000   -0.64339200
H                 -4.44210200   -0.44438200   -1.50536900
H                  0.57336800    0.91984400    1.55935900
O                  1.40804000   -1.50944600    0.82820700
H                  2.30900600   -1.45103000    0.47415200
Cl                 4.37812300   -0.41062800   -0.46888000
H                  3.39955700    0.45055100   -0.26794700

5 13 1.43 S 20 0.05
13 26 1.75 S 20 -0.05
24 25 2.5 S 20 -0.05

老师,我想要得到的是,5 13 断键,24 25 成键, 25 26 断键

QQ截图20151203210834.png (41.34 KB, 下载次数 Times of downloads: 102)

QQ截图20151203210834.png
站在宇宙中心呼唤爱

3097

帖子

29

威望

1万

eV
积分
17221

Level 6 (一方通行)

15#
 楼主 Author| 发表于 Post on 2015-12-1 23:11:24 | 只看该作者 Only view this author
planet5460 发表于 2015-12-1 16:53
寻找过渡态就好像在做艺术工作。

不不不,寻找MECP才是

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
卡开发发 + 2 我很赞同

查看全部评分 View all ratings

本版积分规则 Credits rule

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

GMT+8, 2026-2-22 06:37 , Processed in 0.226274 second(s), 32 queries , Gzip On.

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