计算化学公社

 找回密码 Forget password
 注册 Register
Views: 27201|回复 Reply: 27
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] 如何用高斯做CCSD(T)单点能计算?

[复制链接 Copy URL]

286

帖子

0

威望

1142

eV
积分
1428

Level 4 (黑子)

最近有个体系,不算大,不到60个原子,元素有C,N,O,H,结合了几个水分子。用M06-2X/6-31++G(d,p)优化读的能量,审稿人让做一下CC计算,( I highlyrecommend to run CC calculations (as SP on the DFT structures) to have muchmore definitive conclusions.),还提到相同的体系别人已经用CCSD(T)计算过了。我打算用CCSD(T)做一下单点计算,对比一下能量。请问关键字应该怎么写?还有:
1.基组怎么选?我看别人用的是aug-cc-pVDZ
2.CCSD和CCSD(T)的区别?T是指激发三重态吗?我的体系都是基态,没有激发态。
3.内存应该设置多大?
如果计算单点能的话,直接CCSD(T)/aug-cc-pVDZ 就可以吗?

期待大神回复~~~

评分 Rate

参与人数
Participants 2
eV +5 收起 理由
Reason
大喵说多谢 + 2 谢谢
iU_65018 + 3 赞!

查看全部评分 View all ratings

6万

帖子

99

威望

6万

eV
积分
125190

管理员

公社社长

2#
发表于 Post on 2019-5-2 05:10:15 | 只看该作者 Only view this author
接近六十个原子做CCSD(T)就甭想了,直接告诉审稿人算不动。用像样的基组,标准方式的CCSD(T)至多只能算二十个原子出头的无对称性的有机体系。这审稿人缺乏常识。

顶多也就是用ORCA做DLPNO-CCSD(T)计算,这样的体系才算得动
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

286

帖子

0

威望

1142

eV
积分
1428

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2019-5-2 05:52:33 | 只看该作者 Only view this author
sobereva 发表于 2019-5-2 05:10
接近六十个原子做CCSD(T)就甭想了,直接告诉审稿人算不动。用像样的基组,标准方式的CCSD(T)至多只能算二十 ...

请问 单点也算不动吗?
还有,那个ORCA是免费的吗?我到文献里看到DLPNO-CCSD(T)了,这是个方法吗?容易上手吗?
谢谢!

6万

帖子

99

威望

6万

eV
积分
125190

管理员

公社社长

4#
发表于 Post on 2019-5-2 06:21:05 | 只看该作者 Only view this author
lei234 发表于 2019-5-2 05:52
请问 单点也算不动吗?
还有,那个ORCA是免费的吗?我到文献里看到DLPNO-CCSD(T)了,这是个方法吗?容易 ...

我说的就是单点

是。还算容易上手

aug-cc-pVDZ这种2-zeta档次基组跟CCSD(T)完全不搭,用CCSD(T)必须从3-zeta起步

我不知道你想具体算什么问题,直接从M06-2X/6-31++G(d,p)这样很low的级别(基组太渣)提升到CCSD(T)完全没必要,用个双杂化泛函足够说事了

基组的问题看
谈谈量子化学中基组的选择
http://sobereva.com/336http://bbs.keinsci.com/thread-3545-1-1.html
下文专门提了双杂化和DLPNO-CCSD(T)的事
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272http://bbs.keinsci.com/thread-536-1-1.html

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

286

帖子

0

威望

1142

eV
积分
1428

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2019-5-2 06:52:04 | 只看该作者 Only view this author
sobereva 发表于 2019-5-2 06:21
我说的就是单点

是。还算容易上手

我主要对比结构的能量,涉及到一点过渡态。但不想都算,找两个结构象征性的对比一下能量想着。审稿人建议做CCSD(T),您说的双杂化泛函可以用高斯算吗?我看您博文说的大部分也得用ORCA呢。那我这种体系您建议选哪种双杂化泛函,如果不想学ORCA,只用高斯的话?非常感激!!!

6万

帖子

99

威望

6万

eV
积分
125190

管理员

公社社长

6#
发表于 Post on 2019-5-2 13:15:42 | 只看该作者 Only view this author
lei234 发表于 2019-5-2 06:52
我主要对比结构的能量,涉及到一点过渡态。但不想都算,找两个结构象征性的对比一下能量想着。审稿人建议 ...

非要用Gaussian就用B2PLYP-D3(BJ),关键词写成B2PLYPD3(从g09 D.01开始支持)
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

286

帖子

0

威望

1142

eV
积分
1428

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2019-5-4 01:13:25 | 只看该作者 Only view this author
本帖最后由 lei234 于 2019-5-4 01:21 编辑
sobereva 发表于 2019-5-2 13:15
非要用Gaussian就用B2PLYP-D3(BJ),关键词写成B2PLYPD3(从g09 D.01开始支持)

好的,还有就是基组的问题,对于您说的这种双杂化泛函,应该用什么基组合适呢?我看见您博文里说双杂化泛函的基组要求和MP2相同,至少从6-311G(2d,p)起步,我如果把适合双杂化泛函的基组从低到高排列一下,顺序是不是6-311G(2d,p)<cc-pVTZ,顺便问下,用cc-pVTZ的话,直接B2PLYPD3/cc-pVTZ 就可以吗?感谢!我体系有明显的氢键作用,我用了6-311++g(2d,p),另外加了溶剂化。您看合适吗?

6万

帖子

99

威望

6万

eV
积分
125190

管理员

公社社长

8#
发表于 Post on 2019-5-4 06:23:35 | 只看该作者 Only view this author
lei234 发表于 2019-5-4 01:13
好的,还有就是基组的问题,对于您说的这种双杂化泛函,应该用什么基组合适呢?我看见您博文里说双杂化泛 ...

至少用def2-TZVP,要更好就用def2-TZVPP或cc-pVTZ,再好就用cc-pVQZ或def2-QZVPP

没事甭给氢加弥散,我这里强调了
谈谈弥散函数和“月份”基组
http://sobereva.com/119
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

286

帖子

0

威望

1142

eV
积分
1428

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2019-5-7 05:29:05 | 只看该作者 Only view this author
sobereva 发表于 2019-5-4 06:23
至少用def2-TZVP,要更好就用def2-TZVPP或cc-pVTZ,再好就用cc-pVQZ或def2-QZVPP

没事甭给氢加弥散, ...

我用#p B2PLYPD3/Def2TZVP scrf=(smd,solvent=toluene)算单点,提示以下错误,不知道是怎么回事?还请赐教。

6万

帖子

99

威望

6万

eV
积分
125190

管理员

公社社长

10#
发表于 Post on 2019-5-7 08:13:05 | 只看该作者 Only view this author
lei234 发表于 2019-5-7 05:29
我用#p B2PLYPD3/Def2TZVP scrf=(smd,solvent=toluene)算单点,提示以下错误,不知道是怎么回事?还请赐 ...

硬盘不够
用maxdisk关键词事先告诉Gaussian硬盘剩余空间还多少,省得把硬盘写满了
要么就用ORCA(强烈推荐),比Gaussian算双杂化甚至能快两、三个数量级,而且占硬盘很少,此文靠后部分说了
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272http://bbs.keinsci.com/thread-536-1-1.html
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

286

帖子

0

威望

1142

eV
积分
1428

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2019-5-7 09:52:18 | 只看该作者 Only view this author
sobereva 发表于 2019-5-7 08:13
硬盘不够
用maxdisk关键词事先告诉Gaussian硬盘剩余空间还多少,省得把硬盘写满了
要么就用ORCA(强烈 ...

老师,非常感谢您的回复!
说明书里只讲了不同计算中maxdisk的作用,没有说怎么用....请问具体是怎么用的呀? 是写在关键字那里就可以了吗?还是maxdisk=多少 这样?它的意思是磁盘空间剩多少用多少吗?还是?
很多感谢!

6万

帖子

99

威望

6万

eV
积分
125190

管理员

公社社长

12#
发表于 Post on 2019-5-8 09:55:08 | 只看该作者 Only view this author
lei234 发表于 2019-5-7 09:52
老师,非常感谢您的回复!
说明书里只讲了不同计算中maxdisk的作用,没有说怎么用....请问具体是怎么用 ...

硬盘剩多少就maxdisk写多少,这样Gaussian如果事先知道硬盘空间不够大,要么提前报错,要么自动切换为占硬盘少的算法让计算得以进行
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

11

帖子

0

威望

58

eV
积分
69

Level 2 能力者

13#
发表于 Post on 2020-3-23 18:48:42 | 只看该作者 Only view this author
sobereva 发表于 2019-5-2 05:10
接近六十个原子做CCSD(T)就甭想了,直接告诉审稿人算不动。用像样的基组,标准方式的CCSD(T)至多只能算二十 ...

您好,在上面提到的“标准方式的CCSD(T)至多只能算二十个原子出头的无对称性的有机体系”有相关的文献中说明吗?谢谢您打扰了。

6万

帖子

99

威望

6万

eV
积分
125190

管理员

公社社长

14#
发表于 Post on 2020-3-25 15:43:07 | 只看该作者 Only view this author
bd1225an 发表于 2020-3-23 18:48
您好,在上面提到的“标准方式的CCSD(T)至多只能算二十个原子出头的无对称性的有机体系”有相关的文献中 ...

这是实际计算经验,没有什么文献
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

11

帖子

0

威望

58

eV
积分
69

Level 2 能力者

15#
发表于 Post on 2020-3-25 17:06:26 | 只看该作者 Only view this author
sobereva 发表于 2020-3-25 15:43
这是实际计算经验,没有什么文献

好的
谢谢老师

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-26 13:04 , Processed in 0.196345 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list