计算化学公社

标题: 求助使用sobMECP程序怎么搜寻开壳层单重态和三重态之间的MECP？ [打印本页]

作者
Author: 埋头库库炫饭 时间: 2024-10-3 10:34
标题: 求助使用sobMECP程序怎么搜寻开壳层单重态和三重态之间的MECP？
本帖最后由埋头库库炫饭于 2024-10-3 10:36 编辑

在使用Sob老师的sobMECP程序时我遇到了几个问题，求各位老师和大佬解答一下。1、SobMECP程序只能搜寻闭壳层单重态(CSS)和三重态(Tri)之间MECP吗？2、要搜寻开壳层单重态(OSS)和三重态(Tri)之间的MECP该怎么设置？3、我的这个势能面从int5(Tri)到int6(CSS)存在MECP过程，如果不能搜寻OSS和Tri之间的MECP，那该怎么使用CSS和Tri之间的MECP描述图中这个过程？
(, 下载次数 Times of downloads: 30)

作者
Author: sobereva 时间: 2024-10-3 15:44
这种情况适合用SF-TDDFT，但Gaussian没有。GAMESS-US、ORCA都支持。非要用Gaussian就只能让sobMECP用UKS对称破缺计算

作者
Author: 埋头库库炫饭 时间: 2024-10-3 17:13

sobereva 发表于 2024-10-3 15:44
这种情况适合用SF-TDDFT，但Gaussian没有。GAMESS-US、ORCA都支持。非要用Gaussian就只能让sobMECP用UKS对 ...

感谢老师的回复，由于使用ORCA有些基组不支持，所以还是准备用UKS对称破缺计算。请问后续是只要原输入文件中加上UKS关键词就行了吗？

作者
Author: sobereva 时间: 2024-10-4 00:29

埋头库库炫饭发表于 2024-10-3 17:13
感谢老师的回复，由于使用ORCA有些基组不支持，所以还是准备用UKS对称破缺计算。请问后续是只要原输入文 ...

没那么简单
仔细看下文搞清楚原理
谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82

作者
Author: zjxitcc 时间: 2024-10-4 16:58
不建议使用UKS和SF-TDDFT，赶紧使用GAMESS的MRSF-TDDFT方法是目前的最好办法。参考阅读《Spin-flip方法中RODFT难收敛解决办法》Spin-flip方法中RODFT难收敛解决办法 - 量子化学 (Quantum Chemistry) - 计算化学公社 (keinsci.com)

作者
Author: xxxhhh 时间: 2025-7-6 18:47
您好，我看到您这个也是搜寻涉及过渡态的MECP，想请教一下，我在搜寻三重态和五重态之间的MECP时，感觉结构与预期结构相差很大，并且不收敛，想请教一下看您是怎么做的

作者
Author: 埋头库库炫饭 时间: 2025-7-6 19:35

xxxhhh 发表于 2025-7-6 18:47
您好，我看到您这个也是搜寻涉及过渡态的MECP，想请教一下，我在搜寻三重态和五重态之间的MECP时，感觉结构 ...

我的方法和http://bbs.keinsci.com/forum.php ... p;highlight=sobmecp博文一样。你可以用三重态和五重态的过渡态结构做初猜，都试一遍

作者
Author: xxxhhh 时间: 2025-7-7 16:52
本帖最后由 xxxhhh 于 2025-7-7 16:55 编辑

埋头库库炫饭发表于 2025-7-6 19:35
我的方法和http://bbs.keinsci.com/forum.php?mod=viewthread&tid=865&highlight=sobmecp博文一样。你可 ...

我想找这个黑点的MECP，我试过用两个过渡态作为初猜，用的也是这篇博文的方法，但是最后都跑成了RC，想问问你有这种情况出现吗，而且跑的默认100步都没有收敛

作者
Author: 埋头库库炫饭 时间: 2025-7-8 09:12

xxxhhh 发表于 2025-7-7 16:52
我想找这个黑点的MECP，我试过用两个过渡态作为初猜，用的也是这篇博文的方法，但是最后都跑成了RC，想问 ...

我得两个结构相差不是很大，所以在默认设置下就收敛了。如果你这最终还是还不行，可以试试找3IM1和5IM1之间得MECP，这样反应路径就是5IM1-3IM1-3TS1-3RC

作者
Author: xxxhhh 时间: 2025-7-8 10:27

埋头库库炫饭发表于 2025-7-8 09:12
我得两个结构相差不是很大，所以在默认设置下就收敛了。如果你这最终还是还不行，可以试试找3IM1和5IM1之 ...

好的好的，我的这个就是一个二氧化碳结合的过渡态，我用过渡态跑八九步差不多就变成反应物RC了，那我接下来试着把头部文件改成中间体IM1试试

欢迎光临计算化学公社 (http://bbs.keinsci.com/)