计算化学公社

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

[综合交流] 请问这种大环结构可以用什么软件和方法找过渡态呢?

[复制链接 Copy URL]

516

帖子

1

威望

5052

eV
积分
5588

Level 6 (一方通行)

本帖最后由 niobium 于 2018-1-16 15:59 编辑

如图,请问各位,这两种结构之间的转变,可以用什么方法找到过渡态呢,或者确定一些能量比较高的构型?高斯的scan功能可以实现吗?但好像不太好确定扫描哪个变量。

捕获.PNG (74.78 KB, 下载次数 Times of downloads: 38)

捕获.PNG

2.jpg (538.36 KB, 下载次数 Times of downloads: 25)

2.jpg

1.jpg (559.22 KB, 下载次数 Times of downloads: 23)

1.jpg

DPB3-D3.gjf

7.4 KB, 下载次数 Times of downloads: 4

DPB3-C2.gjf

8.83 KB, 下载次数 Times of downloads: 1

2479

帖子

11

威望

6864

eV
积分
9563

Level 6 (一方通行)

2#
发表于 Post on 2018-1-16 16:03:43 | 只看该作者 Only view this author
大共轭的有机结构一般来说半经验足以应付。软件用XTB或者MOPAC。

如果追求精确结果,则DFT即可。软件的话高斯和ORCA是上佳选择。M06-2X是有机体系的首选泛函,基组6-311G*即可。

264

帖子

0

威望

2588

eV
积分
2852

Level 5 (御坂)

3#
发表于 Post on 2018-1-16 17:43:04 | 只看该作者 Only view this author
本帖最后由 Jasminer 于 2018-1-16 17:45 编辑

找过渡态,这不一定是基元反应吧?找“构型”用构象搜索,方法、软件不少。



516

帖子

1

威望

5052

eV
积分
5588

Level 6 (一方通行)

4#
 楼主 Author| 发表于 Post on 2018-1-16 17:46:37 | 只看该作者 Only view this author
Jasminer 发表于 2018-1-16 17:43
找过渡态,这不一定是基元反应吧?找“构型”用构象搜索,方法、软件不少。

确实应该不像是个基元反应。其实我也不一定是想找个过渡态,而是想说明从一个结构变到另一个结构很困难,不知道可以怎么用计算的证据来支撑。

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

5#
发表于 Post on 2018-1-16 17:48:06 | 只看该作者 Only view this author
不会是基元反应,一步过不去。
你要找能量比较高的构象,可以尝试做高温下的动力学(比如用xtb程序跑MD,或者高斯里的PM6D3结合BOMD),如果能有互变现象,取其中一些帧,分别优化。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

145

帖子

0

威望

3930

eV
积分
4075

Level 6 (一方通行)

6#
发表于 Post on 2018-1-16 18:34:05 | 只看该作者 Only view this author
niobium 发表于 2018-1-16 17:46
确实应该不像是个基元反应。其实我也不一定是想找个过渡态,而是想说明从一个结构变到另一个结构很困难, ...

既然一个左扭,一个右扭,要想扭过去得从不扭的状态过去啊,证明行不容易,想说明困难那就算一下中间态,看一下势垒呗。如果scan计算量会很大吧

1093

帖子

6

威望

6269

eV
积分
7482

Level 6 (一方通行)

7#
发表于 Post on 2018-1-16 19:19:18 | 只看该作者 Only view this author
用targeted MD跑一下动力学,然后找轨迹中间的结构

903

帖子

37

威望

5324

eV
积分
6967

Level 6 (一方通行)

8#
发表于 Post on 2018-1-18 16:57:02 | 只看该作者 Only view this author
本帖最后由 ggdh 于 2018-1-18 17:12 编辑


强行扫描二面角的结果,因为优化的时候最上面的苯苯二面角的扭转垒(74-80对应的二面角,从60到0到-60)跨不过,所以导致了上面的情况

输入文件是修改D3的那个文件如下:
pm6 opt(modredundant,loose)
.....
.....
.....
D 208 209 74 80 S 10 18.000000

有兴趣 其实还是可以继续做的。
手动帮扫描过程中的点越过能垒,或者是对扫描中的每个点做个固定二面角的动力学,找到对应二面角的constrained global min


本版积分规则 Credits rule

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

GMT+8, 2024-11-24 16:46 , Processed in 0.175050 second(s), 24 queries , Gzip On.

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