计算化学公社

标题: 能否详细讲解下CCSD(t)优化的使用注意事项以及技巧？ [打印本页]

作者
Author: mmjjmmjj 时间: 2016-3-10 16:47
标题: 能否详细讲解下CCSD(t)优化的使用注意事项以及技巧？
本帖最后由 mmjjmmjj 于 2016-3-10 19:36 编辑

哪位大侠能否详细讲解下CCSD(T)优化的使用注意事项以及技巧呢？会用，但是感觉使用并不是很得当。

作者
Author: sobereva 时间: 2016-3-10 17:08
要写成CCSD(T)，T大写。
甭用它来优化，CCSD对于优化目的已经足够，CCSD(T)主要是用来算单点的。绝大多数程序里CCSD(T)都没有解析梯度，故只能用有限差分计算梯度，再加上其单点计算量很大，所以也就能优化几个原子的体系。而且Gaussian里用它优化你还得把被优化的几何变量用变量方式表达，麻烦。另外，几何结构对于计算级别敏感度较低，从优化结果上一般很难看出CCSD和CCSD(T)的差异。

作者
Author: mmjjmmjj 时间: 2016-3-10 19:25
本帖最后由 mmjjmmjj 于 2016-3-10 20:00 编辑

哦。您这个太有用了！就是说1. 比较优化后的能量时用CCSD，如果是比较单点再用CCSD(T)处理？
2. 高斯里面算CCSD(T)中给出的CCSD的结果是否等于用CCSD算的结果？3. 用CCSD(T)优化，如何表明优化哪些量？如果没有标明，会有什么后果呢？

谢谢大神！

作者
Author: mmjjmmjj 时间: 2016-3-10 19:28
本帖最后由 mmjjmmjj 于 2016-3-10 19:39 编辑

sobereva 发表于 2016-3-10 17:08
要写成CCSD(T)，T大写。
甭用它来优化，CCSD对于优化目的已经足够，CCSD(T)主要是用来算单点的。绝大多数 ...

另外，最近发表拙作一篇，在文章中应用了Multiwfn。确实非常好用实在的软件！

作者
Author: mmjjmmjj 时间: 2016-3-10 19:35

sobereva 发表于 2016-3-10 17:08
要写成CCSD(T)，T大写。
甭用它来优化，CCSD对于优化目的已经足够，CCSD(T)主要是用来算单点的。绝大多数 ...

另，您能否在思想家公社的门口撰写一点关于CCSD系列的帖子？相信会有很多人收益的！

作者
Author: sobereva 时间: 2016-3-11 06:14

mmjjmmjj 发表于 2016-3-10 19:25
哦。您这个太有用了！就是说1. 比较优化后的能量时用CCSD，如果是比较单点再用CCSD(T)处理？
2. 高斯里面 ...

1 没有所谓的“比较优化后的能量”。你要精确比较两个结构的能量，就先CCSD优化结构，然后CCSD(T)算单点比较能量。
2 等于
3
(, 下载次数 Times of downloads: 35)

作者
Author: sobereva 时间: 2016-3-11 06:15

mmjjmmjj 发表于 2016-3-10 19:35
另，您能否在思想家公社的门口撰写一点关于CCSD系列的帖子？相信会有很多人收益的！

这没有什么好写的，基础理论在所有外文经典量化书里都有，你要想了解更多可以看
(, 下载次数 Times of downloads: 46)

作者
Author: mmjjmmjj 时间: 2016-3-11 16:12

sobereva 发表于 2016-3-11 06:14
1 没有所谓的“比较优化后的能量”。你要精确比较两个结构的能量，就先CCSD优化结构，然后CCSD(T)算单点 ...

比较能量时，如果已经用CCSD算出了结果，为什么还要用CCSD(T)再算一遍单点比较呢？

作者
Author: sobereva 时间: 2016-3-11 16:41

mmjjmmjj 发表于 2016-3-11 16:12
比较能量时，如果已经用CCSD算出了结果，为什么还要用CCSD(T)再算一遍单点比较呢？

CCSD(T)的能量明显比CCSD准啊

作者
Author: 413 时间: 2016-3-27 17:17

sobereva 发表于 2016-3-11 16:41
CCSD(T)的能量明显比CCSD准啊

请教一下UCCSD(T)计算的结果已经经过零点能矫正了吗？为什么在单点计算的输出结果里看不到零点能矫正的信息？
另请教个位数原子用ccsd优化时，频率确认怎么做呢？

作者
Author: sobereva 时间: 2016-4-3 13:48

413 发表于 2016-3-27 17:17
请教一下UCCSD(T)计算的结果已经经过零点能矫正了吗？为什么在单点计算的输出结果里看不到零点能矫正的信 ...

直接计算的是单点能，没有ZPE。你得做freq才能得到ZPE。
单个原子频率计算和分子体系一样，但你没必要做什么频率确认，因为单个原子没有振动，只有平动

作者
Author: 413 时间: 2016-4-3 14:23

sobereva 发表于 2016-4-3 13:48
直接计算的是单点能，没有ZPE。你得做freq才能得到ZPE。
单个原子频率计算和分子体系一样，但你没必要做 ...

请教您一下，fchk文件中Total SCF Density指的是什么？下行输出会有三个相同的数值分别有什么意义呢？

作者
Author: sobereva 时间: 2016-4-3 14:35

413 发表于 2016-4-3 14:23
请教您一下，fchk文件中Total SCF Density指的是什么？下行输出会有三个相同的数值分别有什么意义呢？

对于CCSD(T)，这就是HF的密度矩阵
不知道你说的哪三个数值

作者
Author: 413 时间: 2016-4-3 14:51

sobereva 发表于 2016-4-3 14:35
对于CCSD(T)，这就是HF的密度矩阵
不知道你说的哪三个数值

抱歉，我没表述清楚。
请教该矩阵的矩阵元个数与gaussian type primitive orbital有关系吗？
或该矩阵元个数怎么得来的呢？

作者
Author: sobereva 时间: 2016-4-3 14:58

413 发表于 2016-4-3 14:51
抱歉，我没表述清楚。
请教该矩阵的矩阵元个数与gaussian type primitive orbital有关系吗？
或该矩阵 ...

密度矩阵的唯独是nbasis*nbasis，这里nbasis是指基函数数目。
由于是对称矩阵，所以fch里只记录下三角部分。

作者
Author: 413 时间: 2016-4-3 16:25

sobereva 发表于 2016-4-3 14:58
密度矩阵的唯独是nbasis*nbasis，这里nbasis是指基函数数目。
由于是对称矩阵，所以fch里只记录下三角部 ...

请教总电子密度可否由fchk中的MO coeffs平方算得？

作者
Author: sobereva 时间: 2016-4-4 01:00

413 发表于 2016-4-3 16:25
请教总电子密度可否由fchk中的MO coeffs平方算得？

不是这么算的，这么得到的什么都不是。获得密度你还得考虑基函数的数学表达式

作者
Author: 413 时间: 2016-4-5 09:48

sobereva 发表于 2016-4-4 01:00
不是这么算的，这么得到的什么都不是。获得密度你还得考虑基函数的数学表达式

通过计算density matrix与basis function相乘呢？

作者
Author: sobereva 时间: 2016-4-5 12:02

413 发表于 2016-4-5 09:48
通过计算density matrix与basis function相乘呢？

密度矩阵每个矩阵元和相应的两个基函数的乘积相乘，然后全部加和，则基函数的值算的是r处的，结果就是r处的电子密度。

作者
Author: 413 时间: 2016-4-5 14:57

sobereva 发表于 2016-4-5 12:02
密度矩阵每个矩阵元和相应的两个基函数的乘积相乘，然后全部加和，则基函数的值算的是r处的，结果就是r处 ...

就是核间距为r处的分子的总density？

作者
Author: sobereva 时间: 2016-4-5 15:17

413 发表于 2016-4-5 14:57
就是核间距为r处的分子的总density？

r是指三维空间中的一个点，是矢量

作者
Author: wave0403 时间: 2016-4-9 04:58
提示: 作者被禁止或删除内容自动屏蔽

作者
Author: liuxm 时间: 2017-3-5 12:15
sob老师，真的是受益匪浅。

作者
Author: sxh526 时间: 2019-5-15 09:43

sobereva 发表于 2016-3-10 17:08
要写成CCSD(T)，T大写。
甭用它来优化，CCSD对于优化目的已经足够，CCSD(T)主要是用来算单点的。绝大多数 ...

老师请问审稿人让用CCSD(T)优化几何结构，是否有相关文献可以反驳这一意见呢，因为最小体系都有8个重原子，高斯而且算到一半莫名其妙暂停，实在是算不动，这个小问题一直卡着

作者
Author: sobereva 时间: 2019-5-15 14:28

sxh526 发表于 2019-5-15 09:43
老师请问审稿人让用CCSD(T)优化几何结构，是否有相关文献可以反驳这一意见呢，因为最小体系都有8个重原子 ...

直接说算不动
要么你用PSI4、CFOUR、较新版Molpro之类有CCSD(T)解析梯度的程序去优化

作者
Author: sxh526 时间: 2019-5-15 15:17

sobereva 发表于 2019-5-15 14:28
直接说算不动
要么你用PSI4、CFOUR、较新版Molpro之类有CCSD(T)解析梯度的程序去优化

好的，谢谢老师

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