使用PSI4做对称匹配微扰理论(SAPT)能量分解计算

dgwfly · 发表于 Post on 2022-1-23 04:33:45

在计算离子SCN-和分子C54H18间SAPT能量分解出现了一些问题，在短距离（小于0.3nm）时，用bronze standard 计算出现的induction energy很奇怪，甚至出现了正值，可能是 scale的问题，对E_exch-ind 的计算不是很好。有的文献也指出SAPT对含离子的体系，又是强结合的情况计算的不是很好，更高级的SAPT2+(3)可能能得更好的效果。但我的体系太大，不能用更高级别的SAPT去做测试。大家对计算含离子短距离的体系有什么建议吗？
附件是我的输入输出文件

sobereva · 发表于 Post on 2022-1-23 13:25:32

dgwfly 发表于 2022-1-23 04:33
在计算离子SCN-和分子C54H18间SAPT能量分解出现了一些问题，在短距离（小于0.3nm）时，用bronze standard ...

尝试用SAPT(DFT)，这种大小做得动

Lym · 发表于 Post on 2023-1-12 16:27:59

您好，请问我想利用SAPT讨论溶剂对探针对碱土金属离子结合的贡献，因此计算探针离子(-1;1）和剩余部分(碱土金属+溶剂: +2;1）间SAPT能量分解。在自己笔记本的虚拟机上进行计算的，计算了一半屏幕出现提示出错，想知道如何解决这个问题，麻烦老师指点一下

Ca3BMeCN.inp (2.76 KB, 下载次数 Times of downloads: 2)

Ca3BMeCN.out (141.09 KB, 下载次数 Times of downloads: 1) 屏幕提示:

sobereva · 发表于 Post on 2023-1-13 06:24:43

Lym 发表于 2023-1-12 16:27
您好，请问我想利用SAPT讨论溶剂对探针对碱土金属离子结合的贡献，因此计算探针离子(-1;1）和剩余部分(碱土 ...

可能临时文件把硬盘满了

另外，SAPT0结合def2-TZVPP是很不好的选择，更何况单体还有阴离子的。仔细看下文了解该用什么
使用PSI4做对称匹配微扰理论(SAPT)能量分解计算
http://sobereva.com/526（http://bbs.keinsci.com/thread-15902-1-1.html）

七尺贱 · 发表于 Post on 2023-4-26 13:58:04

sobereva 发表于 2023-1-13 06:24
可能临时文件把硬盘满了

另外，SAPT0结合def2-TZVPP是很不好的选择，更何况单体还有阴离子的。仔细看 ...

想问一下卢老师，SAPT0做能量分解，SAPT0算不算DFT范畴

sobereva · 发表于 Post on 2023-4-27 02:24:00

七尺贱发表于 2023-4-26 13:58
想问一下卢老师，SAPT0做能量分解，SAPT0算不算DFT范畴

明显不算，毫无关系

SAPT(DFT)或DFT-SAPT才跟DFT有关系

七尺贱 · 发表于 Post on 2023-4-27 12:28:11

sobereva 发表于 2023-4-27 02:24
明显不算，毫无关系

SAPT(DFT)或DFT-SAPT才跟DFT有关系

谢谢老师

自由风 · 发表于 Post on 2023-6-29 20:00:59

我刚安装的Psi4计算时报错（见下图），是我未安装好的问题吗？我用sob老师给的算例是可以完成计算的。

sobereva · 发表于 Post on 2023-6-29 23:32:09

自由风发表于 2023-6-29 20:00
我刚安装的Psi4计算时报错（见下图），是我未安装好的问题吗？我用sob老师给的算例是可以完成计算的。

是
写明了没装libecpint模块

自由风 · 发表于 Post on 2023-6-30 09:13:37

sobereva 发表于 2023-6-29 23:32
是
写明了没装libecpint模块

谢谢回复！我是按照您帖子（http://sobereva.com/526）用installer包安装的，没有单独安装libecpint模块，像这种问题还有补救的措施吗？

自由风 · 发表于 Post on 2023-6-30 13:41:38

sobereva 发表于 2023-6-29 23:32
是
写明了没装libecpint模块

我用您网盘提供的1.5版本，安装之后就可以了

自由风 · 发表于 Post on 2023-7-1 09:37:41

我用Psi4做能量分解分析报错，找了很久也没找到类似的情况

sobereva · 发表于 Post on 2023-7-2 11:16:47

自由风发表于 2023-7-1 09:37
我用Psi4做能量分解分析报错，找了很久也没找到类似的情况

甭用M06-2X这种meta系列的泛函

DLJS · 发表于 Post on 2023-8-3 11:38:34

sob老师，您好，我按照老师的方法做能量分解总是中断，没有明显报错，体系七百多个原子，给了180个G是不是太大了算不了导致的呢？

sobereva · 发表于 Post on 2023-8-3 17:41:10

DLJS 发表于 2023-8-3 11:38
sob老师，您好，我按照老师的方法做能量分解总是中断，没有明显报错，体系七百多个原子，给了180个G是不是 ...

强烈建议用sobEDA代替SAPT做大体系能量分解，速度又快（比DFT算单点耗时只高一倍多）又不耗内存和硬盘，结果还非常合理，看原文以及里面的教程：
Tian Lu, Qinxue Chen, A simple, efficient and universal energy decomposition analysis method based on dispersion-corrected density functional theory
ChemRxiv (2023) DOI: 10.26434/chemrxiv-2023-n79rz

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

[PSI4] 使用PSI4做对称匹配微扰理论(SAPT)能量分解计算