新手求助orca DLPNO-CCSD(T)关键词书写

sakura1109

! DLPNO-CCSD(T) def2-TZVPP RIJK def2/JK def2-TZVPP/C normalPNO tightSCF CPCM KDIIS
%maxcore 20000
%pal nprocs 8 end

%CPCM
smd true
SMDsolvent "chlorobenzene"
End


* xyz 0 4坐标
学生刚接触orca使用DLPNO-CCSD(T)算单点能 ，想请大佬们看看以上关键词有没有明显错误？以及有以下几个疑问:
1.在帖子http://bbs.keinsci.com/thread-7062-1-1.html里看到sob老师提过“CCSD(T)别写RI-，速度不会有提升”，请问我需要把 RIJK def2/JK删掉吗？还是我的理解有误(最开始写的def2-TZVPP/JK，报错，在帖子http://bbs.keinsci.com/thread-26235-1-1.html里看到zjxitcc老师说“是def2/JK”，后来又看到http://bbs.keinsci.com/thread-31252-1-1.html冰释之川老师提到“在ORCA里def2-TZVPPD没有配套的辅助基组用于RI-J加速”，那def2-TZVPP也是吗，这个def2是指哪个基组呀？)
2.帖子http://bbs.keinsci.com/thread-20143-1-1.html里 wzkchem5老师说过“分子间的色散用DLPNO-CCSD(T)描述”，因为之前写上D3，会提示WARNING: Grimme dispersion and correlated wavefunction methods
         beyond MP2 are not to be used at the same time!
         Please delete either one from your input.
因此删除了D3关键词，所以其实用DLPNO-CCSD(T)本身就可以描述色散作用吗？
3.帖子http://bbs.keinsci.com/thread-20748-1-1.html里zjxitcc老师说“DLPNO计算应该加TightPNO，默认的有点渣。”DLPNO-CCSD(T)必须要用到tightPNO吗,但是体系有80个原子担心算不动，所以可以保持normalPNO吗？
4.KDIIS是提示不收敛就加上的，似乎它一般要搭配KDIIS+SOSCF使用，目前只加了KDIIS能跑，需要加上SOSCF吗？
体系有80个原子，根据推文https://mp.weixin.qq.com/s/V6vQeROVDulhxdLGFN68mQ预估可能至少一周才能出来结果，先不论结果能用否，目前希望至少关键词没大问题，谢谢各位老师！
还有一个问题5.之前不收敛写Grid5 FinalGrid6提示报错，后来发现帖子http://bbs.keinsci.com/thread-19550-1-1.html中sob老师说“当前写Grid5 FinalGrid6毫无意义，又不是DFT计算”，现在反应过来这个是后HF计算还是MP2呀？

wzkchem5

1. 所谓的不建议写RI只适用于CCSD(T)，不适用于DLPNO-CCSD(T)；DLPNO-CCSD(T)必须用RI。def2/JK适用于所有的不加弥散的def2开头的基组，对于加弥散的基组可能适用也可能不适用，保险起见可以用autoaux。
2. 是的。建议复习一下量化教材，虽然一般教材不讲DLPNO-CCSD(T)，但是会讲CCSD(T)。学懂了原理、知道为什么CCSD(T)自带色散、它和D3的色散有何不同，再去做计算。
3. 当且仅当算得动的时候，用tightPNO。normalPNO精度固然低一些，但是一般仍然比用DLPNO-CCSD(T)以外的方法要准。
4. KDIIS和SOSCF没有任何需要搭配的理由。可能你见过的有限的几个例子里，恰好要么两个都用，两个都不用，但这只是碰巧，或者个别几个人的习惯，并没有任何深层次理由。
5. 自己看量化教材耦合簇计算相关章节，学透了自己就知道了。不然换个稍微不一样的方法，你又不知道是属于什么方法了

sakura1109

wzkchem5 发表于 2023-3-20 18:52
1. 所谓的不建议写RI只适用于CCSD(T)，不适用于DLPNO-CCSD(T)；DLPNO-CCSD(T)必须用RI。def2/JK适用于所有 ...

谢谢老师的回复！确实在原理部分欠缺许多，由于任务繁重所以选择先跑起来再说。老师您能看看我的关键词有没有什么大问题吗？

wzkchem5

sakura1109 发表于 2023-3-20 12:09
谢谢老师的回复！确实在原理部分欠缺许多，由于任务繁重所以选择先跑起来再说。老师您能看看我的关键词有 ...

第一行关键词不用写CPCM（写了是不是错的我不知道，但不写一定是对的），其他没问题

sakura1109

wzkchem5 发表于 2023-3-20 20:01
第一行关键词不用写CPCM（写了是不是错的我不知道，但不写一定是对的），其他没问题

好的好的谢谢老师的回复！我以为写上模块同时也要写上关键词才能调用，目前确实没有任何报错罢了，在别处看到“ORCA5在读取关键词设置的时候 是按书写顺序。以 后出现的设置为准。所以如果开头 !行中书写了关键词，在后面模块部分也书写了，那如果有不一样的设置就按照后面模块中的设置为准。”
这么看确实把关键词中的CPCM删了保险

sakura1109

wzkchem5 发表于 2023-3-20 20:01
第一行关键词不用写CPCM（写了是不是错的我不知道，但不写一定是对的），其他没问题

再多问您一下，目前8核，每个核20000MB，这种情况出现ORCA finished by error termination in SCF应该是SCF不收敛而不是内存问题吧？

hebrewsnabla

sakura1109 发表于 2023-3-20 20:37
再多问您一下，目前8核，每个核20000MB，这种情况出现ORCA finished by error termination in SCF应该是S ...

SCF有没有收敛，可以查看SCF步骤的输出，不是靠猜。

sakura1109

hebrewsnabla 发表于 2023-3-20 20:44
SCF有没有收敛，可以查看SCF步骤的输出，不是靠猜。

我觉得您说的很有道理，我猛然发现自己完全不知道怎么分析orca的输出文件

wzkchem5

sakura1109 发表于 2023-3-20 13:51
我觉得您说的很有道理，我猛然发现自己完全不知道怎么分析orca的输出文件

但凡懂英语，都能看懂。
很多人会觉得，需要先学怎么分析输出文件，再去分析输出文件。其实很多时候并非如此，直接把输出文件当英文阅读理解题来做就行。

mfdsrax2

wzkchem5 发表于 2023-3-20 18:52
1. 所谓的不建议写RI只适用于CCSD(T)，不适用于DLPNO-CCSD(T)；DLPNO-CCSD(T)必须用RI。def2/JK适用于所有 ...

请问一下DLPNO-CCSD(T)为什么必须用RI，我之前的都没用也能算出来，需要重算吗？

zjxitcc

mfdsrax2 发表于 2023-3-21 09:57
请问一下DLPNO-CCSD(T)为什么必须用RI，我之前的都没用也能算出来，需要重算吗？

不用。输出文件里有RI信息，你可以自行查看。参考态HF步骤可以有RI，也可以没有；DLPNO-CC步骤必有RI。

wzkchem5

mfdsrax2 发表于 2023-3-21 02:57
请问一下DLPNO-CCSD(T)为什么必须用RI，我之前的都没用也能算出来，需要重算吗？

DLPNO-CCSD(T)计算即使不写RI，也默认打开RI

sakura1109

wzkchem5 发表于 2023-3-20 21:04
但凡懂英语，都能看懂。
很多人会觉得，需要先学怎么分析输出文件，再去分析输出文件。其实很多时候并非 ...

老师您好，我自己分析了一下输出文件，判断是SCF不收敛，正在尝试不同关键词解决，您能帮我看一下确认一下吗？有怀疑过内存，但是确实没看见关于memory的提示，谢谢

wzkchem5

sakura1109 发表于 2023-3-21 11:03
老师您好，我自己分析了一下输出文件，判断是SCF不收敛，正在尝试不同关键词解决，您能帮我看一下确认一 ...

stderr有输出吗？

sakura1109

wzkchem5 发表于 2023-3-21 18:18
stderr有输出吗？

老师，stderr指的是这个吗