计算化学公社

 找回密码 Forget password
 注册 Register
Views: 521|回复 Reply: 6
打印 Print 上一主题 Last thread 下一主题 Next thread

[PSI4] PSI4能量分解错误,计算的过程中莫名其妙地停下来了

[复制链接 Copy URL]

121

帖子

0

威望

911

eV
积分
1032

Level 4 (黑子)

各位老师,我在用PSI4对(硝仿-和一个阳离子)组成的离子化合物进行SAPT计算,采用的方法来自Multiwfn的金标准(sapt2+(3)dmp2/aug-cc-pVTZ),计算的过程中莫名其妙地停下来了,全程没搜到一个error,最后几行如图所示。
请各位老师指教一下这是怎么回事,谢谢各位老师!(PS:我算了另外一个离子化合物体系也是这样的结果)

202309161449102933..png (62.07 KB, 下载次数 Times of downloads: 27)

202309161449102933..png

121

帖子

0

威望

911

eV
积分
1032

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2023-9-18 21:10:07 | 只看该作者 Only view this author
zjxitcc 发表于 2023-9-18 14:21
内存不足,且体系超出理论方法应用上限。

好的,谢谢老师解惑,学生非常感谢

3813

帖子

4

威望

8002

eV
积分
11895

Level 6 (一方通行)

MOKIT开发者

6#
发表于 Post on 2023-9-18 14:21:22 | 只看该作者 Only view this author
黑色桃花 发表于 2023-9-16 15:53
谢谢老师的解惑,我会修改的,但是想请问一下您,就是这个算不动指的是什么意思,是我的核数不够这些吗还 ...

内存不足,且体系超出理论方法应用上限。
自动做多参考态计算的程序MOKIT

121

帖子

0

威望

911

eV
积分
1032

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2023-9-16 16:08:31 | 只看该作者 Only view this author
sobereva 发表于 2023-9-16 15:56
sapt2+(3)dmp2/aug-cc-pVTZ的耗时和对资源要求极高,对稍微大一点的体系都完全算不动,而且耗内存还极高

...

好的卢老师,我再学学sobEDA,还有就是http://bbs.keinsci.com/thread-10096-1-1.html中对G4MP2-6X有些问题,您有空的话,请您解惑一下,万分感谢卢老师

5万

帖子

99

威望

5万

eV
积分
112482

管理员

公社社长

4#
发表于 Post on 2023-9-16 15:56:51 | 只看该作者 Only view this author
sapt2+(3)dmp2/aug-cc-pVTZ的耗时极高,对稍微大一点的体系都完全算不动,而且耗内存/硬盘还极高,内存/硬盘不够直接就停了。要对不同级别的SAPT的耗时和对资源要求有基本的概念再去用

强烈推荐改用sobEDAw,算几百原子都没压力,结果和高阶SAPT还能很好吻合
使用sobEDA和sobEDAw方法做非常准确、快速、方便、普适的能量分解分析
http://sobereva.com/685http://bbs.keinsci.com/thread-39446-1-1.html


北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

121

帖子

0

威望

911

eV
积分
1032

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2023-9-16 15:53:32 | 只看该作者 Only view this author
zjxitcc 发表于 2023-9-16 15:38
“采用的方法来自Multiwfn的金标准”这话不太妥当,金/银/铜标准是SAPT原文献里提出/建议的,Multiwfn是接 ...

谢谢老师的解惑,我会修改的,但是想请问一下您,就是这个算不动指的是什么意思,是我的核数不够这些吗还是其他原因,请老师再给详细说下,我想了解一下,真的跪谢老师

3813

帖子

4

威望

8002

eV
积分
11895

Level 6 (一方通行)

MOKIT开发者

2#
发表于 Post on 2023-9-16 15:38:44 | 只看该作者 Only view this author
“采用的方法来自Multiwfn的金标准”这话不太妥当,金/银/铜标准是SAPT原文献里提出/建议的,Multiwfn是接受/沿用了这个规则。从目前提供的信息来看是这个体系对于金标准来说太大了,算不动。为验证这点,你可以算个小体系(例如H2O...NH3复合物),若小体系可以正常计算,而当前体系不行,基本上就是了。

解决办法:降低SAPT级别,例如银或铜级别;更换能量分解方法,例如GKS-EDA, sobEDA等等。
自动做多参考态计算的程序MOKIT

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-26 17:33 , Processed in 0.864379 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list