计算化学公社

标题: 请教CCSD(T)搭配基组以及适用分子大小的问题 [打印本页]

作者
Author: bd1225an 时间: 2020-3-7 12:00
标题: 请教CCSD(T)搭配基组以及适用分子大小的问题

我最近有篇投稿论文，是计算一个有机反应，使用的是m062x/6-311++G(2df,2p)计算单点能。最大的结构有41个原子，其中非氢原子是24个。

审稿意见里有一条是审稿人认为我们原来使用的方法精确度不够，提出让用CCSD(T)/6-311++g(2df,2p)进行计算。我印象中CCSD(T)方法计算不了太大的体系，尝试了一下，也确实是出错了，提示是

Leave Link 801 at Wed Mar 4 22:30:33 2020, MaxMem= 262144000 cpu: 34.2

(Enter /THFS/home/g-usx-wc/g09/l804.exe)

Closed-shell transformation, MDV= 262144000 ITran=4 ISComp=1.

Semi-Direct transformation.

ModeAB= 4 MOrb= 66 LenV= 240976548

LASXX= 29099767686 LTotXX= 29099767686 LenRXX= 58289304876

LTotAB= 29189537190 MaxLAS= 27824646960 LenRXY= 0

NonZer= 87389072562 LenScr=175463674368 LnRSAI= 27824646960

LnScr1= 55867586048 LExtra= 0 Total= 317445212252

MaxDsk= -1 SrtSym= T ITran= 4

GetIJB would need an additional 271334551 words of memory to use all 24 processors.

DoSDTr: NPSUse= 11

JobTyp=0 Pass 1: I= 1 to 66.

Complete sort for first half transformation.

Erroneous write. Write -1 instead of 1062160288.

fd = 4

orig len = 1922291360 left = 1062160288

g_write

根据以往的经验似乎是因为磁盘空间不够？但是其实我们这个分子大小我本来也怀疑这种方法算不动。

另外文献（J. Phys. Chem. A. 2020, 124, 90-100.）中提到用CCSD(T)/aug-cc-pVDZ-DK 这种方法能够计算数百个原子，但是aug-cc-pVDZ-DK这个基组我不太了解

我有两个问题向各位大神请教：

1. 如果非要用CCSD(T)计算的话，是否基组降低一些能够满足“41个原子，其中非氢原子24个”这样大小的体系的运算？降低到什么程度比较合适呢？不知道是否有人有这方面的经验。

2. aug-cc-pVDZ-DK这种基组应该如何实现呢？

作者
Author: snljty 时间: 2020-3-7 12:59
本帖最后由 snljty 于 2020-3-7 13:02 编辑

写了一大堆回复给删了，还是看卢老师的博文吧。
这几篇花点时间看一下，补充一下基础知识，上面涉及的一大堆问题就明白了。
方法方面：
简谈量子化学计算中DFT泛函的选择
http://bbs.keinsci.com/thread-536-1-1.html
各种后HF方法精度简单横测
http://bbs.keinsci.com/thread-6005-1-1.html

关于其中提到的ORCA和DLPNO-CCSD(T)/cc-pVTZ等：
详谈Multiwfn产生ORCA量子化学程序的输入文件的功能
http://bbs.keinsci.com/thread-13560-1-1.html
量子化学程序ORCA的安装方法
http://bbs.keinsci.com/thread-11697-1-1.html

关于基组的问题
谈谈量子化学中基组的选择
http://bbs.keinsci.com/thread-3545-1-1.html
谈谈赝势基组的选用
http://bbs.keinsci.com/thread-5625-1-1.html
在线基组和赝势数据库一览
http://bbs.keinsci.com/thread-2081-1-1.html

作者
Author: sobereva 时间: 2020-3-7 14:24
拿6-311++g(2df,2p)搭配CCSD(T)是只有外行人才想得出来的主意，还让你用于41个原子体系，估计这审稿人水平不高
CCSD(T)是不可能结合aug-cc-pVDZ-DK算几百原子的。只有用DLPNO、CIM之类做法才有可能。

-DK版本什么时候才应当用仔细看
在赝势下做波函数分析的一些说明
http://sobereva.com/156

其它该说的楼上都说了

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