请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 3342|回复: 25

[新手求助] 能否详细讲解下CCSD(t)优化的使用注意事项以及技巧?

[复制链接]

65

帖子

0

威望

265

eV
积分
330

Level 3 能力者

发表于 2016-3-10 16:47:26 | 显示全部楼层 |阅读模式
本帖最后由 mmjjmmjj 于 2016-3-10 19:36 编辑

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

1万

帖子

25

威望

2万

eV
积分
43663

管理员

公社社长

发表于 2016-3-10 17:08:27 | 显示全部楼层
要写成CCSD(T),T大写。
甭用它来优化,CCSD对于优化目的已经足够,CCSD(T)主要是用来算单点的。绝大多数程序里CCSD(T)都没有解析梯度,故只能用有限差分计算梯度,再加上其单点计算量很大,所以也就能优化几个原子的体系。而且Gaussian里用它优化你还得把被优化的几何变量用变量方式表达,麻烦。另外,几何结构对于计算级别敏感度较低,从优化结果上一般很难看出CCSD和CCSD(T)的差异。
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

65

帖子

0

威望

265

eV
积分
330

Level 3 能力者

 楼主| 发表于 2016-3-10 19:25:09 | 显示全部楼层
本帖最后由 mmjjmmjj 于 2016-3-10 20:00 编辑

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

谢谢大神!

65

帖子

0

威望

265

eV
积分
330

Level 3 能力者

 楼主| 发表于 2016-3-10 19:28:06 | 显示全部楼层
本帖最后由 mmjjmmjj 于 2016-3-10 19:39 编辑
sobereva 发表于 2016-3-10 17:08
要写成CCSD(T),T大写。
甭用它来优化,CCSD对于优化目的已经足够,CCSD(T)主要是用来算单点的。绝大多数 ...

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

65

帖子

0

威望

265

eV
积分
330

Level 3 能力者

 楼主| 发表于 2016-3-10 19:35:46 | 显示全部楼层
sobereva 发表于 2016-3-10 17:08
要写成CCSD(T),T大写。
甭用它来优化,CCSD对于优化目的已经足够,CCSD(T)主要是用来算单点的。绝大多数 ...

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

1万

帖子

25

威望

2万

eV
积分
43663

管理员

公社社长

发表于 2016-3-11 06:14:05 | 显示全部楼层
mmjjmmjj 发表于 2016-3-10 19:25
哦。您这个太有用了!就是说1. 比较优化后的能量时用CCSD,如果是比较单点再用CCSD(T)处理?
2. 高斯里面 ...

1 没有所谓的“比较优化后的能量”。你要精确比较两个结构的能量,就先CCSD优化结构,然后CCSD(T)算单点比较能量。
2 等于
3
无标题.png
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

1万

帖子

25

威望

2万

eV
积分
43663

管理员

公社社长

发表于 2016-3-11 06:15:24 | 显示全部楼层
mmjjmmjj 发表于 2016-3-10 19:35
另,您能否在思想家公社的门口撰写一点关于CCSD系列的帖子?相信会有很多人收益的!

这没有什么好写的,基础理论在所有外文经典量化书里都有,你要想了解更多可以看
无标题.png
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

65

帖子

0

威望

265

eV
积分
330

Level 3 能力者

 楼主| 发表于 2016-3-11 16:12:28 | 显示全部楼层
sobereva 发表于 2016-3-11 06:14
1 没有所谓的“比较优化后的能量”。你要精确比较两个结构的能量,就先CCSD优化结构,然后CCSD(T)算单点 ...

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

1万

帖子

25

威望

2万

eV
积分
43663

管理员

公社社长

发表于 2016-3-11 16:41:26 | 显示全部楼层
mmjjmmjj 发表于 2016-3-11 16:12
比较能量时,如果已经用CCSD算出了结果,为什么还要用CCSD(T)再算一遍单点比较呢?

CCSD(T)的能量明显比CCSD准啊
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

65

帖子

0

威望

343

eV
积分
408

Level 3 能力者

发表于 2016-3-27 17:17:32 | 显示全部楼层
sobereva 发表于 2016-3-11 16:41
CCSD(T)的能量明显比CCSD准啊

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

1万

帖子

25

威望

2万

eV
积分
43663

管理员

公社社长

发表于 2016-4-3 13:48:02 | 显示全部楼层
413 发表于 2016-3-27 17:17
请教一下UCCSD(T)计算的结果已经经过零点能矫正了吗?为什么在单点计算的输出结果里看不到零点能矫正的信 ...

直接计算的是单点能,没有ZPE。你得做freq才能得到ZPE。
单个原子频率计算和分子体系一样,但你没必要做什么频率确认,因为单个原子没有振动,只有平动
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

65

帖子

0

威望

343

eV
积分
408

Level 3 能力者

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

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

1万

帖子

25

威望

2万

eV
积分
43663

管理员

公社社长

发表于 2016-4-3 14:35:47 | 显示全部楼层
413 发表于 2016-4-3 14:23
请教您一下,fchk文件中Total SCF Density指的是什么?下行输出会有三个相同的数值分别有什么意义呢?

对于CCSD(T),这就是HF的密度矩阵
不知道你说的哪三个数值
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

65

帖子

0

威望

343

eV
积分
408

Level 3 能力者

发表于 2016-4-3 14:51:43 | 显示全部楼层
sobereva 发表于 2016-4-3 14:35
对于CCSD(T),这就是HF的密度矩阵
不知道你说的哪三个数值

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

1万

帖子

25

威望

2万

eV
积分
43663

管理员

公社社长

发表于 2016-4-3 14:58:52 | 显示全部楼层
413 发表于 2016-4-3 14:51
抱歉, 我没表述清楚。
请教该矩阵的矩阵元个数与gaussian type primitive orbital有关系吗?
或该矩阵 ...

密度矩阵的唯独是nbasis*nbasis,这里nbasis是指基函数数目。
由于是对称矩阵,所以fch里只记录下三角部分。
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2019-7-21 07:07 , Processed in 0.156543 second(s), 26 queries .

快速回复 返回顶部 返回列表