计算化学公社

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

[Gaussian/gview] 寻找Mo配合物过渡态的时候,为啥经常会出现有一个水分子跑去第二壳层的情况

[复制链接 Copy URL]

25

帖子

0

威望

109

eV
积分
134

Level 2 能力者

各位老师好!我想请问一下,我在寻找过渡态的时候,为啥经常会出现有一个水分子跑去第二壳层的情况?我该怎么去调整初猜结构呢?谢谢各位老师的指导。我的目的是左边的羟基进攻右边的Mo原子,羟基上的氢脱去与右边的羟基结合形成水分子的过渡态


1685975664617.png (77.38 KB, 下载次数 Times of downloads: 13)

输出的log文件

输出的log文件

MoO2(OH)2(H2O)2-TS-25-20230605.gjf

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

25

帖子

0

威望

109

eV
积分
134

Level 2 能力者

24#
 楼主 Author| 发表于 Post on 2023-6-27 21:57:09 | 只看该作者 Only view this author
wzkchem5 发表于 2023-6-27 21:28
取决于你认为你优化TS的成功率有多高。一次全去掉的话,如果收敛了就省很多时间,但是不收敛概率大;一个 ...

收到谢谢老师

1万

帖子

0

威望

7402

eV
积分
18171

Level 6 (一方通行)

23#
发表于 Post on 2023-6-27 21:28:45 | 只看该作者 Only view this author
abcd401 发表于 2023-6-27 13:39
(1)我限制的就是键长    opt=modredundant    H原子与其他两个原子之间的键长   然后用的是笛卡尔坐标 ...

取决于你认为你优化TS的成功率有多高。一次全去掉的话,如果收敛了就省很多时间,但是不收敛概率大;一个一个去掉的话,即使能收敛也会比较费时间(多做了几次结构优化),但是不收敛概率小
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

25

帖子

0

威望

109

eV
积分
134

Level 2 能力者

22#
 楼主 Author| 发表于 Post on 2023-6-27 20:39:50 | 只看该作者 Only view this author
wzkchem5 发表于 2023-6-27 17:54
(1)单限制一个原子肯定是不行的,等于没限制,因为分子的其他部分可以整体平移。至少得限制一个键长
...

(1)我限制的就是键长    opt=modredundant    H原子与其他两个原子之间的键长   然后用的是笛卡尔坐标优化
(2)老师,抱歉 ! 我重新发布一个新的帖子,回复好像贴不了图片。要是我同时限制了很多个条件,那我后面取消限制用TS优化时,是依次按顺序取消呢,还是所有限制条件同时取消    用TS优化??这个有点疑惑

1万

帖子

0

威望

7402

eV
积分
18171

Level 6 (一方通行)

21#
发表于 Post on 2023-6-27 17:54:59 | 只看该作者 Only view this author
abcd401 发表于 2023-6-27 10:38
老师,我还想请教一下   
①如果限制性优化这个过渡态   是只限制H原子还是要把H和两个O同时限制住??
...

(1)单限制一个原子肯定是不行的,等于没限制,因为分子的其他部分可以整体平移。至少得限制一个键长
(2)图片没贴成功。如果你确定你不希望那个H靠近,就通过限制某个或某几个键长让它不要靠近
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

25

帖子

0

威望

109

eV
积分
134

Level 2 能力者

20#
 楼主 Author| 发表于 Post on 2023-6-27 17:38:09 | 只看该作者 Only view this author
wzkchem5 发表于 2023-6-13 01:07
假如TS存在的话,这个做法确实有可能比柔性扫描更快地找到TS;但假如TS不存在的话,当且仅当初始键长足够 ...

老师,我还想请教一下   
①如果限制性优化这个过渡态   是只限制H原子还是要把H和两个O同时限制住??
②我在把H和O同时限制的时候发现一个问题  旁边水分子上的H会慢慢靠近中间的O    这个问题要怎么解决?C:\Users\Lenovo\Desktop

25

帖子

0

威望

109

eV
积分
134

Level 2 能力者

19#
 楼主 Author| 发表于 Post on 2023-6-14 15:58:25 | 只看该作者 Only view this author
wzkchem5 发表于 2023-6-13 01:09
稳定、乃至结构优化收敛且没有虚频,都不代表是全局能量极小值点,只能代表是局部能量极小值点。构型搜索 ...

好的  非常感谢老师的指导!

1万

帖子

0

威望

7402

eV
积分
18171

Level 6 (一方通行)

18#
发表于 Post on 2023-6-13 01:09:30 | 只看该作者 Only view this author
abcd401 发表于 2023-6-12 14:20
我现在要找MoO2(OH)2(H2O)2聚合成(H2O)2(OH)O2Mo-O-MoO2(OH)(H2O)2的过渡态   MoO2(OH)2(H2O)2单体我是从 ...

稳定、乃至结构优化收敛且没有虚频,都不代表是全局能量极小值点,只能代表是局部能量极小值点。构型搜索的目的在于尽量确保你得到的结构是全局能量最低的结构
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

1万

帖子

0

威望

7402

eV
积分
18171

Level 6 (一方通行)

17#
发表于 Post on 2023-6-13 01:07:27 | 只看该作者 Only view this author
abcd401 发表于 2023-6-12 13:48
好的谢谢!我试着做一下柔性扫描   我前几天跟我导师讨论了一下   他让我 “对这个氢与两个氧之间的键长  ...

假如TS存在的话,这个做法确实有可能比柔性扫描更快地找到TS;但假如TS不存在的话,当且仅当初始键长足够长且优化步长足够小的时候,可以用这个方法证明没有TS,但是还是不如柔性扫描严格。
另外注意我说的柔性扫描是用来解决“证明A到B没有过渡态”的问题的,你关心的应该是“已知A到B有过渡态C,想严格证明不存在另一个(或者一类)过渡态D”,这个正如我所说,无法严格证明,且即使做不严格的证明,也很难是单个的柔性扫描就能搞定的
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

25

帖子

0

威望

109

eV
积分
134

Level 2 能力者

16#
 楼主 Author| 发表于 Post on 2023-6-12 21:20:14 | 只看该作者 Only view this author
wzkchem5 发表于 2023-6-12 15:36
这种反应不能用手摆结构然后直接算,必须做构型搜索,否则结果不能发表。参见http://bbs.keinsci.com/thr ...

我现在要找MoO2(OH)2(H2O)2聚合成(H2O)2(OH)O2Mo-O-MoO2(OH)(H2O)2的过渡态   MoO2(OH)2(H2O)2单体我是从文献中提取坐标出来 自己优化    实际上这个单体是稳定存在的    过渡态是我把两个单体按文献中其他结构一样拼凑在一起再优化的

25

帖子

0

威望

109

eV
积分
134

Level 2 能力者

15#
 楼主 Author| 发表于 Post on 2023-6-12 20:48:33 | 只看该作者 Only view this author
wzkchem5 发表于 2023-6-12 15:35
举例:欲证明自由基A和自由基B复合形成单键是无垒反应,可以把A和B放在足够远的地方,然后做柔性扫描,逐 ...

好的谢谢!我试着做一下柔性扫描   我前几天跟我导师讨论了一下   他让我 “对这个氢与两个氧之间的键长  先做一个opt限制性优化  然后把约束去掉后先做TS优化看会不会靠近羟基   形成水分子”     不知道有没有这个必要?   

1万

帖子

0

威望

7402

eV
积分
18171

Level 6 (一方通行)

14#
发表于 Post on 2023-6-12 15:36:30 | 只看该作者 Only view this author
abcd401 发表于 2023-6-8 09:02
您好,我还想补充一下   这个反应物和产物都是我手动绘制出来的  假想!   但是真正怎样反应的不知道。而 ...

这种反应不能用手摆结构然后直接算,必须做构型搜索,否则结果不能发表。参见http://bbs.keinsci.com/thread-577-1-1.html
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

1万

帖子

0

威望

7402

eV
积分
18171

Level 6 (一方通行)

13#
发表于 Post on 2023-6-12 15:35:07 | 只看该作者 Only view this author
abcd401 发表于 2023-6-8 08:38
首先,谢谢您耐心的解答。我方便问一下  您说的“由A到B的反应,如果可以不经过过渡态直接反应过去(也就 ...

举例:欲证明自由基A和自由基B复合形成单键是无垒反应,可以把A和B放在足够远的地方,然后做柔性扫描,逐渐缩短A-B键长,如果一直扫到A-B完全成键(靠键长判断也行,靠键级判断也行)能量都是单调下降的,那么就是无垒反应。但注意两个前提条件:(1)初始结构A和B必须放得足够远,至于多远算远没有普适标准,我个人的标准是以A和B的范德华表面刚好不接触为限;(2)扫描步长必须足够小,同样地,多小算小没有普适标准,我一般取0.1埃
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

25

帖子

0

威望

109

eV
积分
134

Level 2 能力者

12#
 楼主 Author| 发表于 Post on 2023-6-8 16:02:51 | 只看该作者 Only view this author
wzkchem5 发表于 2023-6-6 23:55
对,所以可能客观上你的体系就是不存在这个水分子不跑到第二壳层的过渡态。
由A到B的反应,如果可以不经 ...

您好,我还想补充一下   这个反应物和产物都是我手动绘制出来的  假想!   但是真正怎样反应的不知道。而且我没有优化反应物     比如单体A   两个A反应得到产物B   我只优化了一个A      没有优化2个A在一起的时候    这个过渡态是我把优化后的A拼凑在一起形成的。  (可能我表达上没有那么的清晰    我想表达的意思是从反应物到产物做柔性扫描,我是对哪个坐标进行柔性扫描呢?

25

帖子

0

威望

109

eV
积分
134

Level 2 能力者

11#
 楼主 Author| 发表于 Post on 2023-6-8 15:38:22 | 只看该作者 Only view this author
wzkchem5 发表于 2023-6-6 23:55
对,所以可能客观上你的体系就是不存在这个水分子不跑到第二壳层的过渡态。
由A到B的反应,如果可以不经 ...

首先,谢谢您耐心的解答。我方便问一下  您说的“由A到B的反应,如果可以不经过过渡态直接反应过去(也就是所谓的无垒反应),一般比较好证明,一般来说适当选取扫描坐标做柔性扫描就可以说明”这个要怎么证明呢?对什么坐标做柔性扫描?

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

GMT+8, 2024-11-27 07:35 , Processed in 0.186388 second(s), 25 queries , Gzip On.

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