求助：用MP2方法算单点，从高斯输出文件中如何查看并计算

哇卡卡卡卡卡

各位老师好！审稿人让我用post HF方法去计算一下单点，因为文章中使用DFT方法计算的两个过渡态的能量差只有0.6 kcal/mol，所以我尝试了MP2方法、MP4方法、CBSQB3方法，还有ORCA的DLPNO-CCST(T)方法，目前存在一些疑惑，请各位老师指导！

①用MP2方法，能够正常结束，但是我不太清楚这两个过渡态的能量差是多少，我不知道如何去统计，在这个帖子中提到了查看用MP2方法计算的电子能量谈谈该从Gaussian输出文件中的什么地方读电子能量 - 思想家公社的门口：量子化学·分子模拟·二次元 (sobereva.com)，但是具体还是不太清楚应该怎么算。关于两个过渡态的out文件已经随附件上传。
②MP4方法，出现的错误是g_write，硬盘内存的原因，一开始遇到这个问题的时候我将高斯计算的临时文件/Scratch夹改成了我的主目录，结果还是会出现这种报错，这种情况应该怎么解决？
③用ORCA的DLPNO-CCSD(T)方法计算两个过渡态的单点，虽然结尾都是ORCA TERMINATED NORMALLY，但是我查看单点能量，波函数是不完全收敛的（FINAL SINGLE POINT ENERGY     -3088.079387890857   (Wavefunction not fully converged!)），我看网上都是直接出现SCF不收敛，查阅ORCA手册，尝试加了slowconv，以及降低收敛限，最终的结果都是一样的（正常结束，但是波函数不完全收敛），不知道这种情况应该怎么办？关于用ORCA计算的输出文件已随附件上传

希望各位老师方便的时候可以指点我一下，十分感谢！

sobereva

直接上DLPNO-CCSD(T)就完了，根本没必要考虑其它的。MP2的结果通常还不如DFT（泛函选择恰当时）

DLPNO-CCSD(T)计算前先手动做个HF计算，并且用论坛、手册里提到的各种方法把SCF弄收敛了。必要时也可以读取Gaussian产生的波函数，看下文。HF顺利做完了之后，DLPNO-CCSD(T)任务再去读取HF的gbw文件里的波函数当初猜
将Gaussian等程序收敛的波函数作为ORCA的初猜波函数的方法
http://sobereva.com/517（http://bbs.keinsci.com/thread-14905-1-1.html）

哇卡卡卡卡卡

sobereva 发表于 2023-5-15 19:47
直接上DLPNO-CCSD(T)就完了，根本没必要考虑其它的。MP2的结果通常还不如DFT（泛函选择恰当时）

DLPNO-C ...

谢谢sobereva老师，那MP2方法计算得到的输出文件，具体应该怎么看呢，我想也统计看一下能量差

zjxitcc

哇卡卡卡卡卡 发表于 2023-5-15 20:47
谢谢sobereva老师，那MP2方法计算得到的输出文件，具体应该怎么看呢，我想也统计看一下能量差

《谈谈该从Gaussian输出文件中的什么地方读电子能量》http://sobereva.com/488

哇卡卡卡卡卡

zjxitcc 发表于 2023-5-16 10:32
《谈谈该从Gaussian输出文件中的什么地方读电子能量》http://sobereva.com/488

老师，我根据帖子还有高斯手册，查到了对应的E值
老师我还有一个疑惑：我应该如何得到G值呢，第一次接触post HF方法，前面用DFT方法计算时，算freq时会得到gibbs校正值，那么在post HF方法中，我觉得就不能用前面的这个gibbs校正值了吧，是需要重新用MP2方法算一个freq还是用HF方法去算freq，或者都不是，希望老师方便的时候指导我一下

zjxitcc

哇卡卡卡卡卡 发表于 2023-5-16 15:45
老师，我根据帖子还有高斯手册，查到了对应的E值
老师我还有一个疑惑：我应该如何得到G值呢，第一次接触 ...

在此文中搜“粗人”和“讲究的人”
《谈谈隐式溶剂模型下溶解自由能和体系自由能的计算》http://sobereva.com/327

哇卡卡卡卡卡

zjxitcc 发表于 2023-5-16 15:48
在此文中搜“粗人”和“讲究的人”
《谈谈隐式溶剂模型下溶解自由能和体系自由能的计算》http://soberev ...

谢谢老师的指导！