计算化学公社

标题: 高斯计算CCSD（T）计算单点能和频率的续算问题 [打印本页]

作者
Author: danr 时间: 2023-9-21 10:11
标题: 高斯计算CCSD（T）计算单点能和频率的续算问题
用高斯的ccsd(t)做单点能和频率计算的时候，计算了很多天了，然后断电了，怎么进行续算

作者
Author: niobium 时间: 2023-9-21 10:15
多大的体系？单点建议转投ORCA，用DLPNO-CCSD(T)，频率不知道你做的目的是什么？

作者
Author: danr 时间: 2023-9-21 13:09

niobium 发表于 2023-9-21 10:15
多大的体系？单点建议转投ORCA，用DLPNO-CCSD(T)，频率不知道你做的目的是什么？

体系不大，体系一共八九个原子，但是计算资源不是很多。也算了一周多了。频率计算是想要体系的自由能。就想知道怎么续算，学习一下，以后也好用。

作者
Author: zjxitcc 时间: 2023-9-21 13:35

danr 发表于 2023-9-21 13:09
体系不大，体系一共八九个原子，但是计算资源不是很多。也算了一周多了。频率计算是想要体系的自由能。就 ...

问一下，为啥没提到结构优化，你结构优化收敛了，已经进入了freq阶段？那具体是CCSD振幅迭代阶段，还是(T)能量计算阶段，抑或是计算数值导数中某一阶段？

作者
Author: sobereva 时间: 2023-9-21 14:12
opt freq撑死了用CCSD也就够了，用解析梯度都没有的CCSD(T)做等于自虐，而且连红外强度都得不到，而且结果比CCSD也几乎不会有可察觉的进步。尤其是算自由能的目的，热力学校正量带来的误差相对于电子能量的误差来说只是个零头，用CCSD(T)完全是费力不讨好
用CCSD(T)做振动分析一定要真正懂得自己在做什么，不是真懂就绝对不要做

作者
Author: niobium 时间: 2023-9-21 20:15

zjxitcc 发表于 2023-9-21 13:35
问一下，为啥没提到结构优化，你结构优化收敛了，已经进入了freq阶段？那具体是CCSD振幅迭代阶段，还是(T ...

目测他不知道opt和freq要在同一级别。

作者
Author: danr 时间: 2023-9-22 10:31

sobereva 发表于 2023-9-21 14:12
opt freq撑死了用CCSD也就够了，用解析梯度都没有的CCSD(T)做等于自虐，而且连红外强度都得不到，而且结果 ...

谢谢sob老师

作者
Author: danr 时间: 2023-9-22 10:34

niobium 发表于 2023-9-21 20:15
目测他不知道opt和freq要在同一级别。

所以就是说freq的话，必须和opt的方法一致是吗。不能就是说我用DFT优化的结构，然后用CC来计算freq得到自由能，这样是不准确的。但是DFT的结构做cc的单点是可以的吧

作者
Author: danr 时间: 2023-9-22 10:35

zjxitcc 发表于 2023-9-21 13:35
问一下，为啥没提到结构优化，你结构优化收敛了，已经进入了freq阶段？那具体是CCSD振幅迭代阶段，还是(T ...

额，想着用DFT的优化的结构来着

作者
Author: zjxitcc 时间: 2023-9-22 10:47
本帖最后由 zjxitcc 于 2023-9-22 10:49 编辑

danr 发表于 2023-9-22 10:34
所以就是说freq的话，必须和opt的方法一致是吗。不能就是说我用DFT优化的结构，然后用CC来计算freq得到自 ...

单点是可以。直接freq不行，因为DFT优化的几何结构在CCSD(T)下力不为零，未经CCSD(T)的opt、直接算freq得到热力学量无意义。

作者
Author: danr 时间: 2023-11-25 11:22

zjxitcc 发表于 2023-9-22 10:47
单点是可以。直接freq不行，因为DFT优化的几何结构在CCSD(T)下力不为零，未经CCSD(T)的opt、直接算freq得 ...

我想问一下，我用DFT优化，用CCSD（T）计算单点能，那么CCSD（T）的单点能的零点能修正是加上DFT的零点能进行修正吗

作者
Author: zjxitcc 时间: 2023-11-25 13:24

danr 发表于 2023-11-25 11:22
我想问一下，我用DFT优化，用CCSD（T）计算单点能，那么CCSD（T）的单点能的零点能修正是加上DFT的零点能 ...

是的

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