计算化学公社

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

[新手求助] 搜索弱相互作用过渡态的问题

[复制链接 Copy URL]

21

帖子

0

威望

85

eV
积分
106

Level 2 能力者

各位好,想请教过渡态相关的问题。构型一的弱相互作用类型是氢键,构型二弱相互作用类型是碳键。想找弱相互作用的转化过程中经历了怎样的过渡态,反应机理。初猜总是不合理,QSTN的方法也有尝试,也不理想,求指教。构型三是自己想出的中间体,优化后结构稳定,可以参考。

构型一.png (24.17 KB, 下载次数 Times of downloads: 24)

弱相互作用类型氢键

弱相互作用类型氢键

构型二.png (18.04 KB, 下载次数 Times of downloads: 34)

弱相互作用类型碳键

弱相互作用类型碳键

构型三.png (21.15 KB, 下载次数 Times of downloads: 21)

中间体

中间体

6万

帖子

99

威望

6万

eV
积分
125140

管理员

公社社长

2#
发表于 Post on 2017-8-14 16:30:59 | 只看该作者 Only view this author
复合物构型肯定不止这三种,再尝试看能否找到其它极小点构型(可以借助molclus),以其为参考再找过渡态。
注意帖子标题要反映帖子内容,给你改了
北京科音自然科学研究中心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

21

帖子

0

威望

85

eV
积分
106

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2017-8-14 19:46:00 | 只看该作者 Only view this author
sobereva 发表于 2017-8-14 16:30
复合物构型肯定不止这三种,再尝试看能否找到其它极小点构型(可以借助molclus),以其为参考再找过渡态。
...

好的,谢谢sob老师,我用这个程序试一试,下次发帖子注意。

22

帖子

0

威望

444

eV
积分
466

Level 3 能力者

4#
发表于 Post on 2017-8-16 08:57:24 | 只看该作者 Only view this author
请问一下,不同弱相互作用构型之间也有过渡态吗?若存在的话,过渡态是是不是有虚频呢?

6万

帖子

99

威望

6万

eV
积分
125140

管理员

公社社长

5#
发表于 Post on 2017-8-16 12:49:26 | 只看该作者 Only view this author
小巨人 发表于 2017-8-16 08:57
请问一下,不同弱相互作用构型之间也有过渡态吗?若存在的话,过渡态是是不是有虚频呢?

有,有
北京科音自然科学研究中心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

21

帖子

0

威望

85

eV
积分
106

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2017-8-17 09:51:54 | 只看该作者 Only view this author
在研究,开始不明白构想搜索,找不到合理的构型,sob老师说的先搜索合理的构象,然后在着基础上寻找过渡态,很好的方法,受益。谢谢

21

帖子

0

威望

85

eV
积分
106

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2017-9-5 15:53:15 | 只看该作者 Only view this author
Sob老师,您好。我最近研究了gentor程序,由构型一到构型二,在设置gentor里setting中,怎么选择旋转才合理呢。

构型一.png (24.44 KB, 下载次数 Times of downloads: 22)

构型一.png

构型二.png (18.31 KB, 下载次数 Times of downloads: 19)

构型二.png

6万

帖子

99

威望

6万

eV
积分
125140

管理员

公社社长

8#
发表于 Post on 2017-9-6 05:42:13 | 只看该作者 Only view this author
茶新味 发表于 2017-9-5 15:53
Sob老师,您好。我最近研究了gentor程序,由构型一到构型二,在设置gentor里setting中,怎么选择旋转才合理 ...

二聚体适合genmer而不是gentor
北京科音自然科学研究中心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

21

帖子

0

威望

85

eV
积分
106

Level 2 能力者

9#
 楼主 Author| 发表于 Post on 2017-9-6 21:15:30 | 只看该作者 Only view this author
Sob老师,您好,我使用genmer寻找二聚体的分子构象。最开始设置生成20个团簇数目,出现isostat-1,然后我把生成的cluster保存成traj.xyz,用molclus调用GAUSS产生isostat2,产生五个构型,如图三,我想找构型一二聚体之间形成的氢键弱相互作用到构型5形成的碳键弱相互作用,寻找不同相互作用之间过渡态。从genmer得到的结果来看,这二者存在过渡态吗?这种想法可行吗?

isostat-1.png (39.72 KB, 下载次数 Times of downloads: 28)

isostat-1.png

isostat-2.png (22.63 KB, 下载次数 Times of downloads: 21)

isostat-2.png

图三.png (91.03 KB, 下载次数 Times of downloads: 19)

图三.png

6万

帖子

99

威望

6万

eV
积分
125140

管理员

公社社长

10#
发表于 Post on 2017-9-7 02:04:52 | 只看该作者 Only view this author
茶新味 发表于 2017-9-6 21:15
Sob老师,您好,我使用genmer寻找二聚体的分子构象。最开始设置生成20个团簇数目,出现isostat-1,然后 ...


1,2,3本质都一样,只需要考虑一个。
虽然当前已经考虑了20种初始构型,但不一定已经全部都极小点都已经找到了,不排除还有其它极小点的可能。

4-5之间的过渡态根据当前两个结构的信息,应该容易能给出合适的初猜。但4与1/2/3之间结构差异还比较大,恐怕只经历一个过渡态还过不去,你可以先猜一些过渡态算算试试。

另外,一定注意合理选择计算级别。
北京科音自然科学研究中心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

21

帖子

0

威望

85

eV
积分
106

Level 2 能力者

11#
 楼主 Author| 发表于 Post on 2017-9-11 09:55:26 | 只看该作者 Only view this author
本帖最后由 茶新味 于 2017-9-11 10:35 编辑

Sob,老师,您好。我又从新考虑50个构型,最后算出来还是类似的以上三种构型。应该是全部找到,但是我在手动猜测过渡态的过程,总会出现负十几或者几十的虚频,然后走了IRC,有的时候只走几步,有的时候曲线震荡强烈。如图1和2和3。图三主要是沿着CH3F的三个H原子转动。老师,我计算的这种不同弱相互作用之间会有合理的过渡态吗?做了好几个月,越做越没自信。

图一.png (41.2 KB, 下载次数 Times of downloads: 20)

图一.png

图二.png (54.85 KB, 下载次数 Times of downloads: 14)

图二.png

图三.png (70.66 KB, 下载次数 Times of downloads: 28)

图三.png

6万

帖子

99

威望

6万

eV
积分
125140

管理员

公社社长

12#
发表于 Post on 2017-9-11 17:21:40 | 只看该作者 Only view this author
茶新味 发表于 2017-9-11 09:55
Sob,老师,您好。我又从新考虑50个构型,最后算出来还是类似的以上三种构型。应该是全部找到,但是我在手 ...


弱相互作用过渡态附近势能面相当平滑,处理起来并不容易,必须在数值精度层面上达到较高要求
诸如找过渡态、走IRC的时候都用int=ultrafine结合calcall,并且找过渡态时候用tight收敛限,可有效避免数值噪音。但没有万能的办法,免不了折腾一番。
北京科音自然科学研究中心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

21

帖子

0

威望

85

eV
积分
106

Level 2 能力者

13#
 楼主 Author| 发表于 Post on 2017-9-14 14:11:25 | 只看该作者 Only view this author
Sob老师,您好。弱相互作用之前过渡态虚频负十几或者几十可以吗?这种过渡态好像并不像有机反应机理一样有很大的虚频,走IRC如果可以走通,可以认为就是合理过渡态吗?

6万

帖子

99

威望

6万

eV
积分
125140

管理员

公社社长

14#
发表于 Post on 2017-9-14 14:38:22 | 只看该作者 Only view this author
茶新味 发表于 2017-9-14 14:11
Sob老师,您好。弱相互作用之前过渡态虚频负十几或者几十可以吗?这种过渡态好像并不像有机反应机理一样有 ...

弱相互作用体系势能面平缓,对应的构型转变的过渡态在这种虚频大小是很正常的,IRC走通了就算成功
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2026-2-20 18:31 , Processed in 0.173620 second(s), 23 queries , Gzip On.

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