计算化学公社

标题: [1,5]氢转移过渡态计算问题 [打印本页]

作者
Author: 李海恩 时间: 2023-1-12 14:14
标题: [1,5]氢转移过渡态计算问题
最近在尝试计算脂肪醛光激发 Norrish反应中的[1,5]氢转移过渡态，找了一些参考文献（例如DOI：10.26434/chemrxiv.12941702.v1）

请教各位大佬：我从SI复制了文献给的TS结构坐标，想用G16自己再算一次（泛函基组与文献一致），但最终得到的结构和文献中差距很大

请问这是什么原因导致呢？是使用软件不同？或输入文件我写得不对？

作者
Author: zjxitcc 时间: 2023-1-12 14:40
文献里有说这个过渡态是在S1还是T1下找的么？

作者
Author: 李海恩 时间: 2023-1-12 16:02

zjxitcc 发表于 2023-1-12 14:40
文献里有说这个过渡态是在S1还是T1下找的么？

S1和T1都给了！我这个是S1的坐标

作者
Author: fatpig 时间: 2023-1-12 18:53
估计是别人用了open-shell singlet，而你用了closed-shell singlet...

作者
Author: liyuanhe211 时间: 2023-1-12 20:50
从一般规律来说这个过渡态应该是存在的，但是你自己的结构明显是找过渡态失败回到底物了（就是不依赖文献，不管哪个版本的高斯，只要计算级别不太离谱，应该都能找到类似长相的过渡态），所以你自己的计算应该是有一定问题。

作者
Author: zjxitcc 时间: 2023-1-12 21:36

李海恩发表于 2023-1-12 16:02
S1和T1都给了！我这个是S1的坐标

有两种可能性：
（1）该体系在当前计算级别（CAM-B3LYP/def2TZVP）下基态是闭壳层单重态，则方法不要写UCAM，写CAM或RCAM。由于你写了UCAM，关键词singlets就失效了，算的未必是S1，而且还多费了一些机时。
（2）该体系在当前计算级别下基态是开壳层单重态，不能用这一行关键词企图一步计算达到目的，先检验一下基态波函数稳定性，排除或确定此种可能，再谈后续计算。

当前关键词既不符合（1），也不符合（2），不合理。

作者
Author: 李海恩 时间: 2023-1-13 12:35

zjxitcc 发表于 2023-1-12 21:36
有两种可能性：
（1）该体系在当前计算级别（CAM-B3LYP/def2TZVP）下基态是闭壳层单重态，则方法不要写U ...

好的我试试看，谢谢大佬！

作者
Author: 李海恩 时间: 2023-1-13 14:07

liyuanhe211 发表于 2023-1-12 20:50
从一般规律来说这个过渡态应该是存在的，但是你自己的结构明显是找过渡态失败回到底物了（就是不依赖文献， ...

谢谢大佬！已经修改关键词再算一次了T T

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