CCSD(T)算不动较大体系有没有什么解决办法

mfdsrax2 · 发表于 Post on 2021-11-8 20:08:38

我需要计算反应能垒，有一个中间体特别大，有6个硝基，42个非氢原子（只含CNO），其他体系是用CCSD(T)/def2TZVPP//M062X/6-311+g(2d,p)算的，这个有点算不动，有没有什么解决办法

wzkchem5 · 发表于 Post on 2021-11-8 20:25:28

用orca的DLPNO-CCSD(T)算，计算时间大概也就是DFT的10倍的量级。这么大的体系，高斯之类的软件不可能算得动。
因为DLPNO引入了一些近似，绝对能量和普通的CCSD(T)不可比，所以其他的CCSD(T)单点也要用orca重算才能比较。不过opt+freq不需要重算

zjxitcc · 发表于 Post on 2021-11-8 20:26:40

这都2021年快要结束了，DLPNO-CCSD(T)已经很普及了，下载个ORCA，用上DLPNO-CCSD(T)，岂非爽哉？
《用ORCA做DLPNO-CCSD(T)计算》，《ORCA 5.0安装及运行》

zjxitcc · 发表于 Post on 2021-11-8 20:28:51

wzkchem5 发表于 2021-11-8 20:25
用orca的DLPNO-CCSD(T)算，计算时间大概也就是DFT的10倍的量级。这么大的体系，高斯之类的软件不可能算得动 ...

闭壳层DLPNO-CCSD(T1) VeryTightSCF TightPNO搭配同样的基组，再开RI，也应该是可以比的吧，都超级接近了。

红米饭1234 · 发表于 Post on 2021-11-8 21:05:03

反应位垒没必要用后HF做单点，基组改大一点就行了。推荐 M062X//6-311+g(2d,p)//M062X/6-311g*

sobereva · 发表于 Post on 2021-11-8 21:10:40

用DLPNO-CCSD(T)/cc-pVTZ，看下文相应部分
详谈Multiwfn产生ORCA量子化学程序的输入文件的功能
http://sobereva.com/490（http://bbs.keinsci.com/thread-13560-1-1.html）

如果算着还吃力，在ORCA里用PWPB95-D3(BJ)/def2-QZVPP，看上文，用Multiwfn可以直接产生输入文件，非常方便

PS：如果对势垒精度要求较高，M062X/6-311+g(2d,p)级别不够用。

wzkchem5 · 发表于 Post on 2021-11-8 22:08:29

zjxitcc 发表于 2021-11-8 13:28
闭壳层DLPNO-CCSD(T1) VeryTightSCF TightPNO搭配同样的基组，再开RI，也应该是可以比的吧，都超级接近了 ...

反正如果是我的话，除非针对我的体系自己做benchmark确认一下，否则不太敢这么减。。。
而且这样做的话写文章也不好写，写我只用高斯或者只用orca都不反映实际情况，如果写有的用高斯有的用orca又会被审稿人喷得体无完肤。。。

mfdsrax2 · 发表于 Post on 2021-11-8 23:00:32

sobereva 发表于 2021-11-8 21:10
用DLPNO-CCSD(T)/cc-pVTZ，看下文相应部分
详谈Multiwfn产生ORCA量子化学程序的输入文件的功能
http://so ...

PS：如果对势垒精度要求较高，M062X/6-311+g(2d,p)级别不够用。

有哪些推荐的优化级别？

sobereva · 发表于 Post on 2021-11-8 23:01:36

mfdsrax2 发表于 2021-11-8 23:00
PS：如果对势垒精度要求较高，M062X/6-311+g(2d,p)级别不够用。

有哪些推荐的优化级别？

我针对的是5L说的，这个级别对于算能量不够用，优化完全足够。实际上加弥散都完全多余

mfdsrax2 · 发表于 Post on 2021-12-4 20:56:14

zjxitcc 发表于 2021-11-8 20:28
闭壳层DLPNO-CCSD(T1) VeryTightSCF TightPNO搭配同样的基组，再开RI，也应该是可以比的吧，都超级接近了 ...

你好，我试了
! DLPNO-CCSD(T1) tightPNO RIJK def2-TZVPP def2-TZVPP/JK def2-TZVPP/C verytightSCF noautostart miniprint nopop

但是提示下面的问题，def2-TZVPP不需要JK吗？
         !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
                              INPUT ERROR
         UNRECOGNIZED OR DUPLICATED KEYWORD(S) IN SIMPLE INPUT LINE
                     DEF2-TZVPP/JK
         !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

biogon · 发表于 Post on 2021-12-4 22:23:32

mfdsrax2 发表于 2021-12-4 20:56
你好，我试了
! DLPNO-CCSD(T1) tightPNO RIJK def2-TZVPP def2-TZVPP/JK def2-TZVPP/C verytightSCF no ...

自动生成关键词一下先确认版本

zjxitcc · 发表于 Post on 2021-12-4 22:24:52

mfdsrax2 发表于 2021-12-4 20:56
你好，我试了
! DLPNO-CCSD(T1) tightPNO RIJK def2-TZVPP def2-TZVPP/JK def2-TZVPP/C verytightSCF no ...

是def2/JK

mfdsrax2 · 发表于 Post on 2021-12-4 23:11:08

zjxitcc 发表于 2021-12-4 22:24
是def2/JK

好的可以了，程序运行时有两个警告，这个有影响吗？

================================================================================
                                    WARNINGS
                     Please study these warnings very carefully!
================================================================================

WARNING: MDCI localization with Augmented Hessian Foster-Boys
  ===> : Switching off randomization!

WARNING: Post HF methods need fully converged wavefunctions
  ===> : Setting SCFConvForced true
      You can overwrite this default with %scf ConvForced false

INFO : the flag for use of LIBINT has been found!

mfdsrax2 · 发表于 Post on 2021-12-4 23:13:20

zjxitcc 发表于 2021-12-4 22:24
是def2/JK

还有一个问题，RIJK def2/JK和def2-TZVPP/C都是加速方法吗，如果不写是否可以提高精度？

zjxitcc · 发表于 Post on 2021-12-4 23:43:41

mfdsrax2 发表于 2021-12-4 23:13
还有一个问题，RIJK def2/JK和def2-TZVPP/C都是加速方法吗，如果不写是否可以提高精度？

看如下几篇，了解相关的基本知识
《在ORCA的SCF计算中使用RI近似》
《RI-MP2简介及用ORCA做RI-MP2计算》
《用ORCA做DLPNO-CCSD(T)计算》

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] CCSD(T)算不动较大体系有没有什么解决办法