计算化学公社

 找回密码 Forget password
 注册 Register

PSI4 SAPT2+运算报错

查看数: 2477 | 评论数: 14 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2022-7-24 21:30

正文摘要:

我用PSI4计算的时候出现了以下的报错内容,想请教一下原因是什么呢?我该如何改正? 输入文件内容如下: memory 10 GB molecule dimer { 0 1 C      1.40580462   -2.299252 ...

回复 Reply

fy_111 发表于 Post on 2025-7-17 21:00:50
zjxitcc 发表于 2025-7-17 19:59
不合理。这个问题不在于使用哪种EDA方法,而在于你这么划分是否合理,是否能研究你的问题。打个不太恰当 ...

谢谢老师的解答!老师的意思是对于这种多环芳烃的结构难以直接划分出分子片段,不适合通过能量分解去量化分子片段的相互作用,不知道我的理解是否有误?
那我是否可以采用其他方法间接量化作用力的强弱,比如AIM临界点电子密度。
注:这两个片段间同时存在静电引力(羰基和胺)和立体排斥(原子距离小于范德华半径和)相互作用
zjxitcc 发表于 Post on 2025-7-17 19:59:22
fy_111 发表于 2025-7-17 18:58
老师,那采用“使用Multiwfn做基于分子力场的能量分解分析 http://sobereva.com/442 ”这个方法是不是更 ...

不合理。这个问题不在于使用哪种EDA方法,而在于你这么划分是否合理,是否能研究你的问题。打个不太恰当的比方,你要用电锯割小麦,我说用电锯不合适,你反过来问我你该用哪个牌子的电锯。
fy_111 发表于 Post on 2025-7-17 18:58:37
zjxitcc 发表于 2025-7-17 18:39
这个划分片段的方式看起来有点粗暴,连共轭环都切了,建议到PSI4官方论坛提问https://forum.psicode.org/ ...

老师,那采用“使用Multiwfn做基于分子力场的能量分解分析 http://sobereva.com/442 ”这个方法是不是更合理呢?基于分子已经优化好的电子结构去分析,不存在额外的片段计算。
zjxitcc 发表于 Post on 2025-7-17 18:39:00
fy_111 发表于 2025-7-17 16:36
老师您好,看到您关于片段划分的回复,我意识到我最近在使用PSI的F/I-SAPT时一直报错的原因可能和片段划 ...

这个划分片段的方式看起来有点粗暴,连共轭环都切了,建议到PSI4官方论坛提问https://forum.psicode.org/。如果有其他手段可以研究你的问题,应当优先选择其他更恰当的手段,否则结果变成“你敢算、审稿人不敢信”,那最后也是白忙活。
fy_111 发表于 Post on 2025-7-17 16:36:16
zjxitcc 发表于 2022-7-26 10:59
少氢的方式太离谱了,片段2含奇数个电子,不可能是单重态。能量分解有很多种方法,但无论做哪种,大前提 ...

老师您好,看到您关于片段划分的回复,我意识到我最近在使用PSI的F/I-SAPT时一直报错的原因可能和片段划分不合理有关,特此想请老师帮忙看一下我重新划分的片段(如图:黄、绿、灰三个片段,分别在红色虚线处断键)及电荷及自旋多重度的设置(三个片段分别为:0 1,0 1,0 1)是否合理,谢谢老师!


202507171635571371..png (117.3 KB, 下载次数 Times of downloads: 3)

片段划分

片段划分
zjxitcc 发表于 Post on 2022-9-25 18:33:38
星斗如盘 发表于 2022-9-25 18:04
谢谢,大体系只能用超算了,请问100个原子级别,使用PSI4进行EDA时,一般推荐什么级别进行计算呢?(如果 ...
《使用PSI4做对称匹配微扰理论(SAPT)能量分解计算》http://sobereva.com/526
星斗如盘 发表于 Post on 2022-9-25 18:04:27
zjxitcc 发表于 2022-9-25 18:01
所有量化程序都遵从memory内存,disk硬盘

谢谢,大体系只能用超算了,请问100个原子级别,使用PSI4进行EDA时,一般推荐什么级别进行计算呢?(如果超算性能够的话)
zjxitcc 发表于 Post on 2022-9-25 18:01:55
星斗如盘 发表于 2022-9-25 17:59
请教一下,PSI4的memory指的是运行内存,还是硬盘容量?

所有量化程序都遵从memory内存,disk硬盘
星斗如盘 发表于 Post on 2022-9-25 17:59:41
zjxitcc 发表于 2022-7-24 22:06
(1)报错内容没有(正确地)展示出来;(2)10GB内存太少了,少说100GB。

请教一下,PSI4的memory指的是运行内存,还是硬盘容量?
zjxitcc 发表于 Post on 2022-7-26 10:59:01
Luckynan 发表于 2022-7-26 10:50
好的,因为这两个片段都是在一个化合物上截取出来的,因此看起来少氢。我是想分析一个化合物上的两个结构 ...

少氢的方式太离谱了,片段2含奇数个电子,不可能是单重态。能量分解有很多种方法,但无论做哪种,大前提都是片段的几何结构、电荷和自旋多重度要合理设置,不能随意切键、指认,即使程序不报错,算出来也未必合理。

你当前的目的适合F-SAPT或I-SAPT,见https://psicode.org/psi4manual/master/fisapt.html
又或GKS-EDA(划分成3个片段),见《GKS-EDA计算简介》。
Luckynan 发表于 Post on 2022-7-26 10:50:37
zjxitcc 发表于 2022-7-25 17:21
(1)想回复他人,应该在对应楼层点击“回复”,否则别人不会收到回复消息;
(2)你的两个片段结构有问 ...

好的,因为这两个片段都是在一个化合物上截取出来的,因此看起来少氢。我是想分析一个化合物上的两个结构之间的SAPT分析,该如何修改输入文件呢?
zjxitcc 发表于 Post on 2022-7-25 17:21:16
Luckynan 发表于 2022-7-25 16:11
非常抱歉,下面是报错内容,非常感谢
--> dimer = geometry("""
    0 1

(1)想回复他人,应该在对应楼层点击“回复”,否则别人不会收到回复消息;
(2)你的两个片段结构有问题,结构似乎都少了(>=1个)氢原子。
zjxitcc 发表于 Post on 2022-7-24 22:06:49
(1)报错内容没有(正确地)展示出来;(2)10GB内存太少了,少说100GB。

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

GMT+8, 2025-8-13 03:08 , Processed in 0.156188 second(s), 25 queries , Gzip On.

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