用G4MP2算有机分子的热力学常数总出错，应如何解决？

逗比哒蓝蓝 · 发表于 Post on 2019-9-13 04:31:59

本帖最后由逗比哒蓝蓝于 2019-9-13 04:44 编辑

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

LT-III-054.gjf (2.58 KB, 下载次数 Times of downloads: 44)

LT-III-054.zip (379.81 KB, 下载次数 Times of downloads: 7)

逗比哒蓝蓝 · 发表于 Post on 2019-9-13 04:35:00

本帖最后由逗比哒蓝蓝于 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

谢谢！

sobereva · 发表于 Post on 2019-9-13 07:00:00

36个原子体系用G4(MP2)根本算不动，别算了

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

逗比哒蓝蓝 · 发表于 Post on 2019-9-13 08:25:34

谢谢老师的回复。另外想问一下，我这个体系是不是不适合用热力学组合方法，只能用低级别算热校正量，高级别算单点能，最后两者相加的方法来算热力学量？

卡开发发 · 发表于 Post on 2019-9-13 10:48:28

本帖最后由卡开发发于 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）

逗比哒蓝蓝 · 发表于 Post on 2019-9-13 11:07:23

明白了，多谢各位老师。

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

[新手求助] 用G4MP2算有机分子的热力学常数总出错，应如何解决？

浏览过的版块