wzkchem5 发表于 2023-10-22 04:41 好的好的,谢谢老师的耐心指导,学生之后尝试用ORCA计算一下canonical CCSD(T)再进行对比,再次感谢您! |
zh7729 发表于 2023-10-21 15:46 你用orca算canonical CCSD(T)不是更好吗,不然要是审稿人比较刁钻,说你无法排除是DLPNO和canonical的差别虽然存在、但是和orca和高斯的差别抵消了,你怎么回答 |
wzkchem5 发表于 2023-10-21 20:26 谢谢老师的回复,学生明白了很多。还有一个问题就是我算的体系比较大,Gaussian的canonical CCSD(T)只有在def2-TZVPP下才能算得动,但因为这个基组比较小,所以我想通过ORCA的DLPNO-CCSD(T)把能量外推到CBS,首先我需要证明使用DLPNO在相同基组下计算出的能量与canonical CCSD(T)在精度上损失不大,那么可以直接把Gaussian算出的基于UHF的CCSD(T)能量和ORCA算出的基于ROHF计算出的DLPNO-CCSD(T)能量进行对比吗?谢谢您的耐心回复!麻烦您了 |
zh7729 发表于 2023-10-21 03:52 (1)对比是可以对比的,只不过对比结果意义有限,既不能用来说明ROHF reference和UHF reference的差异,也不能说明用不用DLPNO的差异,更不能说明orca和高斯的差异。一定要控制变量做对比的话,注意orca也可以做canonical CCSD(T),可以用来排除一些变量 (2)虽然可以强行从DLPNO-CCSD(T)能量里把UHF能量减掉,但是剩下的项的物理意义比较模糊,不能叫做相关能。UHF和UHF、ROHF和ROHF是可以在orca和高斯之间比较的,但为了严格可比,需要按http://sobereva.com/573的方法严格控制变量,才能严格比较 |
wzkchem5 发表于 2023-10-20 23:21 谢谢老师的回复,学生明白了。我想再问一下:(1)接下来我想比较一下ORCA的DLPNO-CCSD(T)和Gaussian的canonical CCSD(T)在def2-TZVPP基组下的能量精度差异,但Gaussian是以UHF为基础进行计算的,ORCA是以ROHF为基础计算的,这种情况下两者可以直接进行能量的对比吗?(2)另外如果我想分别把Gaussian和ORCA的CCSD(T)能量拆分成HF部分和电子相关部分后再进行比较,是不是Gaussian和ORCA在HF部分(分别为UHF和ROHF)和电子相关部分(分别以UHF和ROHF为基础)的能量没有可比性呢?如果没有可比性的话,我想把DLPNO-CCSD(T)的能量按照输出的结果拆分成UHF部分和电子相关能的部分,这样的操作是合理的吗?问题比较多,麻烦老师了! |
对,DLPNO-CCSD(T)可以理解为在QRO参考态的基础上做了一个ROHF的CCSD(T)计算,所以相关能是以ROHF计算为基础定义的,extrapolate HF能量的时候取的也是这一步ROHF计算的能量。 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-26 09:40 , Processed in 0.264716 second(s), 25 queries , Gzip On.