计算化学公社

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

[CASTEP/Dmol3/MS] MS搜索过渡态自动配对出问题

[复制链接 Copy URL]

34

帖子

0

威望

87

eV
积分
121

Level 2 能力者

我是用MS中的Dmol3模块进行过渡态搜索的,在TS Search前需要对原子进行配对,但是利用其中的auto find进行配对之后,发现许多原子配对都十分混乱。由于我的模型是高岭土表面,因此原子数很多,auto find之后发现许多配对原子都不对,但是要是一个个检查的话又很麻烦。请问有没有什么方法可以让他配对的成功率高一些呢?我认为我的表面前后变化并没有非常大啊

3753

帖子

3

威望

1万

eV
积分
19669

Level 6 (一方通行)

围观吃瓜群众

21#
发表于 Post on 2021-9-9 18:11:33 | 只看该作者 Only view this author
ysmmsy8 发表于 2021-9-9 09:34
下图1是震动频率和强度

下图2是-10.6的频率对应的模型,与产物结构类似,只是结构角度的微调

有个原因,可能是晶面太小。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

8

帖子

0

威望

69

eV
积分
77

Level 2 能力者

20#
发表于 Post on 2021-9-9 16:31:50 | 只看该作者 Only view this author
wzkchem5 发表于 2021-9-9 15:32
第一确保初猜结构和过渡态足够近,第二把优化步长调得比平时小,再试

好滴,谢谢老师的解答~

1万

帖子

0

威望

8991

eV
积分
20739

Level 6 (一方通行)

19#
发表于 Post on 2021-9-9 15:32:49 | 只看该作者 Only view this author
ysmmsy8 发表于 2021-9-9 02:36
我根据这个过渡态搜索得到的轨迹图做了好几个中间体的初猜模型,结构优化后,不是反应物的结构,就是产物 ...

第一确保初猜结构和过渡态足够近,第二把优化步长调得比平时小,再试
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
本团队长期招收研究生,有意者可私信联系

8

帖子

0

威望

69

eV
积分
77

Level 2 能力者

18#
发表于 Post on 2021-9-9 09:36:23 | 只看该作者 Only view this author
wzkchem5 发表于 2021-9-8 21:48
建议从这个过渡态往产物的反方向优化,尝试找到中间体,再用TS search找这个中间体和反应物之间的过渡态

我根据这个过渡态搜索得到的轨迹图做了好几个中间体的初猜模型,结构优化后,不是反应物的结构,就是产物的结构。。。。。。

8

帖子

0

威望

69

eV
积分
77

Level 2 能力者

17#
发表于 Post on 2021-9-9 09:34:37 | 只看该作者 Only view this author
卡开发发 发表于 2021-9-9 00:07
可以,具体看看图和详细信息吧,比如振动模式等。


下图1是震动频率和强度

下图2是-10.6的频率对应的模型,与产物结构类似,只是结构角度的微调
下图3从左往右分别是反应物、过渡态、产物

下图4是过渡态搜索得到的能量

3753

帖子

3

威望

1万

eV
积分
19669

Level 6 (一方通行)

围观吃瓜群众

16#
发表于 Post on 2021-9-9 00:07:41 | 只看该作者 Only view this author
ysmmsy8 发表于 2021-9-8 19:46
产生的过渡态的频率计算有两个虚频,一个频率是-140,对应的强度为186.57 km/mol;一个频率是-10.6,对应 ...

可以,具体看看图和详细信息吧,比如振动模式等。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

1万

帖子

0

威望

8991

eV
积分
20739

Level 6 (一方通行)

15#
发表于 Post on 2021-9-8 21:48:38 | 只看该作者 Only view this author
ysmmsy8 发表于 2021-9-8 14:21
还没有进行TS confirm,只进行了TS search,得到的能量图如下图1所示,过渡态的能量基本介于反应物与产物 ...

建议从这个过渡态往产物的反方向优化,尝试找到中间体,再用TS search找这个中间体和反应物之间的过渡态
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
本团队长期招收研究生,有意者可私信联系

8

帖子

0

威望

69

eV
积分
77

Level 2 能力者

14#
发表于 Post on 2021-9-8 21:21:31 | 只看该作者 Only view this author
wzkchem5 发表于 2021-9-8 20:13
只是虚频振动方向不合理,还是TS confirm结果不合理?
当虚频很小的时候,虚频振动方向容易和直觉预期不 ...

还没有进行TS confirm,只进行了TS search,得到的能量图如下图1所示,过渡态的能量基本介于反应物与产物的能量之间,而并不是要比后两者高,由此做出判断这个过渡态搜索出来的结构不是过渡态,并且搜索出的过渡态与产物的结构很相似,图二所示(从左到右分别是反应物、过渡态、产物)

202109082121003479..png (115.94 KB, 下载次数 Times of downloads: 79)

202109082121003479..png

202109082121273234..png (48.33 KB, 下载次数 Times of downloads: 72)

202109082121273234..png

1万

帖子

0

威望

8991

eV
积分
20739

Level 6 (一方通行)

13#
发表于 Post on 2021-9-8 20:13:31 | 只看该作者 Only view this author
ysmmsy8 发表于 2021-9-8 12:46
产生的过渡态的频率计算有两个虚频,一个频率是-140,对应的强度为186.57 km/mol;一个频率是-10.6,对应 ...

只是虚频振动方向不合理,还是TS confirm结果不合理?
当虚频很小的时候,虚频振动方向容易和直觉预期不一致,但这不一定代表过渡态找错了。比如一根键断裂,可能快要全断的时候有一个二面角转了一下,然后接着才全断,这时虚频可能是二面角转动的方向,容易误导你认为过渡态找错了,但是沿着两边扫IRC或做结构优化就会发现这个过渡态和你的目标反应物和产物确实是相连的。
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
本团队长期招收研究生,有意者可私信联系

8

帖子

0

威望

69

eV
积分
77

Level 2 能力者

12#
发表于 Post on 2021-9-8 19:46:12 | 只看该作者 Only view this author
卡开发发 发表于 2021-9-8 13:20
最好先通过频率计算确定一下过渡态是不是过渡态,是的话再用TS confirm看看当中是否还存在其他中间体。

产生的过渡态的频率计算有两个虚频,一个频率是-140,对应的强度为186.57 km/mol;一个频率是-10.6,对应的强度为3480.91km/mol。但两个虚频对应的原子运动轨迹都不合理,且过渡态的能量基本和反应物与产物的能量持平。(老师,方便加您好友发详细的图片吗?)

3753

帖子

3

威望

1万

eV
积分
19669

Level 6 (一方通行)

围观吃瓜群众

11#
发表于 Post on 2021-9-8 13:20:20 | 只看该作者 Only view this author
ysmmsy8 发表于 2021-9-7 15:38
老师,如果产生的新的过渡态结构和产物的结构类似,考虑反应物与过渡态之间是否还存在其他的过渡态中间体 ...

最好先通过频率计算确定一下过渡态是不是过渡态,是的话再用TS confirm看看当中是否还存在其他中间体。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

8

帖子

0

威望

69

eV
积分
77

Level 2 能力者

10#
发表于 Post on 2021-9-7 15:38:06 | 只看该作者 Only view this author
卡开发发 发表于 2020-5-28 14:41
看新的结构与原来反应物产物是否结构相近。如果相近说明其中有个结构没优化好,取而代之来和另一个建立路 ...

老师,如果产生的新的过渡态结构和产物的结构类似,考虑反应物与过渡态之间是否还存在其他的过渡态中间体,因此把过渡态的模型拿出来进行结构优化后,还是与原来的产物结构超级相似,遇到这种情况改如何调整呢?

34

帖子

0

威望

87

eV
积分
121

Level 2 能力者

9#
 楼主 Author| 发表于 Post on 2020-5-28 15:29:29 | 只看该作者 Only view this author
卡开发发 发表于 2020-5-28 14:41
看新的结构与原来反应物产物是否结构相近。如果相近说明其中有个结构没优化好,取而代之来和另一个建立路 ...

哦哦,我知道了,太感谢了♪(・ω・)ノ

3753

帖子

3

威望

1万

eV
积分
19669

Level 6 (一方通行)

围观吃瓜群众

8#
发表于 Post on 2020-5-28 14:41:14 | 只看该作者 Only view this author
youtuBo 发表于 2020-5-28 11:33
那请问是进行构型优化之后,对得到的构型进行重新配对搜索过渡态吗?是将哪两个进行配对啊

看新的结构与原来反应物产物是否结构相近。如果相近说明其中有个结构没优化好,取而代之来和另一个建立路径;如果不同于另外两个,说明原来的反应可能并非基元反应,分别在新的结构和原来两个结构建立路径分别搜索过渡态。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

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

GMT+8, 2025-8-14 14:54 , Processed in 1.265283 second(s), 24 queries , Gzip On.

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