计算化学公社

标题: 求助计算C3H4-3+反应路径的算法和基组的选择 [打印本页]

作者
Author: lyt0226 时间: 2021-6-4 20:45
标题: 求助计算C3H4-3+反应路径的算法和基组的选择
各位老师好，

如果想计算C3H4-3+碎裂成两、三个碎片这样的过程的反应路径应该选择什么算法和基组比较合适？

已发表的文献中用B3LYP/6-311G优化过渡态和中间态，我用同样的方法也成功找到了过渡态，但是试图换成用m062X/def2tzvpp去优化，却优化不到这个过渡态，是这个算法和基组不合适吗？

作者
Author: wzkchem5 时间: 2021-6-4 21:59
你的方法比文献的方法合理，文献的结果不可信。基组必须带极化函数才能用。
但是因为这个体系很小，其实完全可以用MRCISD+Q做，至少也可以用CCSD(T)，比M062X靠谱得多。

作者
Author: 喵星大佬 时间: 2021-6-4 22:29

wzkchem5 发表于 2021-6-4 21:59
你的方法比文献的方法合理，文献的结果不可信。基组必须带极化函数才能用。
但是因为这个体系很小，其实完 ...

MRCISD+Q操作起来太麻烦了，在瓜丝里至少可以用双杂化优化CCSD(T)单点

作者
Author: wzkchem5 时间: 2021-6-5 00:34

喵星大佬发表于 2021-6-4 15:29
MRCISD+Q操作起来太麻烦了，在瓜丝里至少可以用双杂化优化CCSD(T)单点

我觉得CCSD(T)/cc-pVTZ优化都没问题，单点用CCSD(T) TZ/QZ外推

作者
Author: lyt0226 时间: 2021-6-5 11:00

wzkchem5 发表于 2021-6-5 00:34
我觉得CCSD(T)/cc-pVTZ优化都没问题，单点用CCSD(T) TZ/QZ外推

谢谢您！

您说的基组外推指的是sob老师发表的“http://sobereva.com/172”这个内容吗？

作者
Author: lyt0226 时间: 2021-6-5 11:02

喵星大佬发表于 2021-6-4 22:29
MRCISD+Q操作起来太麻烦了，在瓜丝里至少可以用双杂化优化CCSD(T)单点

谢谢您！

那我如果用B2LYP-D3（BJ）/def2TZVPP优化，然后用CCSD（T）/aug-cc-pvtz进行单点计算,您看合理吗？

作者
Author: hdhxx123 时间: 2021-6-5 11:04

lyt0226 发表于 2021-6-5 11:02
谢谢您！

那我如果用B2LYP-D3（BJ）/def2TZVPP优化，然后用CCSD（T）/aug-cc-pvtz进行单点计算,您看合 ...

你说的泛函是B2PLYP-D3（BJ）还是B3LYP-D3（BJ）。。

作者
Author: lyt0226 时间: 2021-6-5 11:07

hdhxx123 发表于 2021-6-5 11:04
你说的泛函是B2PLYP-D3（BJ）还是B3LYP-D3（BJ）。。

我是试图在sob老师给的算法里面找到相对合适的“http://bbs.keinsci.com/thread-536-1-1.html”，里面列出的是B2PLYP-D3（BJ）

作者
Author: wzkchem5 时间: 2021-6-5 14:33

lyt0226 发表于 2021-6-5 04:00
谢谢您！

您说的基组外推指的是sob老师发表的“http://sobereva.com/172”这个内容吗？

对，外推以后误差可以比单独cc-pVTZ计算出来结果小好几倍

作者
Author: biogon 时间: 2021-6-5 16:59

喵星大佬发表于 2021-6-4 22:29
MRCISD+Q操作起来太麻烦了，在瓜丝里至少可以用双杂化优化CCSD(T)单点

用MOKIT解愁

作者
Author: 北大-陶豫 时间: 2021-6-5 18:08
琢磨了半天没懂 C3H4-3+ 是个什么东西……

作者
Author: 喵星大佬 时间: 2021-6-5 18:14

北大-陶豫发表于 2021-6-5 18:08
琢磨了半天没懂 C3H4-3+ 是个什么东西……

估计是质谱里的碎片离子吧

作者
Author: lyt0226 时间: 2021-6-6 09:45
本帖最后由 lyt0226 于 2021-6-6 09:47 编辑

北大-陶豫发表于 2021-6-5 18:08
琢磨了半天没懂 C3H4-3+ 是个什么东西……

就是丙二烯三价离子，它被电子碰撞或离子碰撞，或者激光辐射后电离掉三个电子，之后会由于原子间库伦排斥势发生碎裂，碎裂成不同的碎片，就是想研究后面的这个碎裂过程，里面可能涉及一些次序解离之类的过程

作者
Author: lyt0226 时间: 2021-6-6 10:37

wzkchem5 发表于 2021-6-5 00:34
我觉得CCSD(T)/cc-pVTZ优化都没问题，单点用CCSD(T) TZ/QZ外推

您好，我用CCSD/cc-pVTZ优化的时候总是出现“No analytic 2nd derivatives for this method”的报错，是什么原因？

作者
Author: biogon 时间: 2021-6-6 14:37

lyt0226 发表于 2021-6-6 10:37
您好，我用CCSD/cc-pVTZ优化的时候总是出现“No analytic 2nd derivatives for this method”的报错，是 ...

用数值梯度优化

作者
Author: snljty 时间: 2021-6-6 15:05

lyt0226 发表于 2021-6-6 10:37
您好，我用CCSD/cc-pVTZ优化的时候总是出现“No analytic 2nd derivatives for this method”的报错，是 ...

只是优化反应物不会有这个报错。如果是优化过渡态，看卢老师博文《简谈Gaussian里找过渡态的关键词opt=TS和QST2、3》http://bbs.keinsci.com/thread-12093-1-1.html ，在里面搜“值得一提的是”

作者
Author: lyt0226 时间: 2021-6-6 15:17

snljty 发表于 2021-6-6 15:05
只是优化反应物不会有这个报错。如果是优化过渡态，看卢老师博文《简谈Gaussian里找过渡态的关键词opt=TS ...

好的，看到了，谢谢！

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