计算化学公社

标题: 求助过渡态和IRC分析是否可以做到CCSD（T）层次 [打印本页]

作者
Author: liurui 时间: 2023-2-6 21:22
标题: 求助过渡态和IRC分析是否可以做到CCSD（T）层次
老师们，想问一下过渡态和IRC分析是否可以做到CCSD（T）层次来进行呢？需要算隧穿概率等问题，如果是算TS用ccsd(t)，molpro.gaussian.Q-Chem等大部分软件得到的过渡态频率都是数值频率，过渡态结构的数值频率可以当做解析频率用么，亦或者用解析频率的CFOUR，但是未找到怎样做过渡态IRC这种分析，想问有什么好办法么？谢谢各位老师

作者
Author: chands 时间: 2023-2-6 22:24
不能用CCSD(T)算IRC，太耗时了，算不动。

作者
Author: wjc404 时间: 2023-2-6 22:30
本帖最后由 wjc404 于 2023-2-6 22:50 编辑

如果有Gaussian和CFOUR，可以试试这个链接上提供的方法http://bbs.keinsci.com/forum.php ... =%B8%DF%D0%A1%CB%C4。
需要注意的是，用CCSD(T)需要搭配的基组至少是cc-pVTZ起步（或者如果有F12方法则从cc-pVDZ-F12起步），不然基组不完备的误差远超理论方法本身的误差，就没有意义了。在这种前提下，普通工作站能够算的动解析频率的有机体系一般不超过7个第二周期原子。对大多数体系而言，CCSD(T)存在的意义是做单点能校正。

作者
Author: sobereva 时间: 2023-2-7 02:24
说数值频率还是解析频率没意义，没有什么任务只要求必须是解析频率。

找过渡态、产生IRC用到双杂化程度就已经很精确了，撑死了CCSD也够了，用CCSD(T)毫无必要，除非是个位数原子的体系

不应脱离实际研究场景、实际意义，去过度追究计算技术问题

作者
Author: wzkchem5 时间: 2023-2-7 02:34
可以算，但极慢。
一般算隧穿系数的方法，都支持用低精度方法（例如DFT）算IRC，再用高精度方法（例如CCSD(T)）在IRC上面算单点。个人觉得这样引入的误差很可能小于用IRC计算隧穿系数这种做法本身的误差了，本来用IRC计算隧穿系数就是一种近似，严格一些的做法需要用PIMD甚至波包动力学等。

作者
Author: liurui 时间: 2023-2-7 17:30

chands 发表于 2023-2-6 22:24
不能用CCSD(T)算IRC，太耗时了，算不动。

收到，谢谢老师！

作者
Author: liurui 时间: 2023-2-7 17:30

wjc404 发表于 2023-2-6 22:30
如果有Gaussian和CFOUR，可以试试这个链接上提供的方法http://bbs.keinsci.com/forum.php ... =%B8%DF%D0%A ...

收到，谢谢老师，俺去试一下

作者
Author: liurui 时间: 2023-2-7 17:31

sobereva 发表于 2023-2-7 02:24
说数值频率还是解析频率没意义，没有什么任务只要求必须是解析频率。

找过渡态、产生IRC用到双杂化程度 ...

收到，谢谢老师！

作者
Author: liurui 时间: 2023-2-7 17:31

wzkchem5 发表于 2023-2-7 02:34
可以算，但极慢。
一般算隧穿系数的方法，都支持用低精度方法（例如DFT）算IRC，再用高精度方法（例如CCSD ...

收到，谢谢老师！

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