计算化学公社

标题: 用G4MP2算有机分子的热力学常数总出错，应如何解决？ [打印本页]

作者
Author: 逗比哒蓝蓝 时间: 2019-9-13 04:31
标题: 用G4MP2算有机分子的热力学常数总出错，应如何解决？
本帖最后由逗比哒蓝蓝于 2019-9-13 04:44 编辑

各位老师好，我正在用G4(MP2)算质子化三苯基氧膦的热力学常数，尝试过几次都没有成功。out文件里显示结构优化已经完成了，但是之后马上就停止了。这里上传了gjf和out文件，请各位老师帮忙看一下问题在哪里。谢谢！
(, 下载次数 Times of downloads: 44) (, 下载次数 Times of downloads: 7)

作者
Author: 逗比哒蓝蓝 时间: 2019-9-13 04:35
本帖最后由逗比哒蓝蓝于 2019-9-13 04:44 编辑

out文件的末尾附在这里。

ANorm= 0.1392992745D+01
E2 = -0.2557223766D+01 EUMP2 = -0.11089293691523D+04
Keep R2 and R3 ints in memory in canonical form, NReq=3224595433.
CIDS: In Core Option IDoMem=          1.
Iterations=  50 Convergence= 0.100D-06
Produce multiple copies of IABC intergrals
Iteration Nr. 1
**********************
DD1Dir will call FoFMem 1 times, MxPair=    2450
NAB=  1225 NAA=    0 NBB=    0.
MP4(R+Q)=  0.68881760D-01
201665.50user 513.66system 2:22:30elapsed 2364%CPU (0avgtext+0avgdata 61438844maxresident)k
63539248inputs+250256712outputs (978major+101799482minor)pagefaults 0swaps

谢谢！

作者
Author: sobereva 时间: 2019-9-13 07:00
36个原子体系用G4(MP2)根本算不动，别算了

另外，那叫热力学量，不叫热力学常数

作者
Author: 逗比哒蓝蓝 时间: 2019-9-13 08:25
谢谢老师的回复。另外想问一下，我这个体系是不是不适合用热力学组合方法，只能用低级别算热校正量，高级别算单点能，最后两者相加的方法来算热力学量？

作者
Author: 卡开发发 时间: 2019-9-13 10:48
本帖最后由卡开发发于 2019-9-13 10:49 编辑

逗比哒蓝蓝发表于 2019-9-13 08:25
谢谢老师的回复。另外想问一下，我这个体系是不是不适合用热力学组合方法，只能用低级别算热校正量，高级别 ...

低级别opt freq（比如M062X/6-31G*）+scale factor来处理热力学校正量，然后ORCA用双杂化（RI-PWPB95-D3/def2-QZVPP）来计算单点，30几原子原则上没问题，结果应该也不错。（参考http://sobereva.com/214）

作者
Author: 逗比哒蓝蓝 时间: 2019-9-13 11:07
明白了，多谢各位老师。

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