求助：基于键临界点的电子密度估计内氢键键能

yangna11

各位老师好，我打算采用基于临界点的电子密度估计单体A和单体B之间的氢键键能，复合物的结构优化后计算单点，再进行AIM拓扑分析，我想请教的是，做AIM拓扑分析所用的.fch文件是算单点能产生的，可以加隐式溶剂模型吗，因为我是研究复合物在水溶液中的相互作用。
此外，我还用sobEDAw对复合物的能量进行了分解，这里能量的分解是针对于真空环境而言的，（几何结构在溶剂模型中优化，单点任务在真空中进行），那么我上面的AIM拓扑分析用的.fch文件是否也需要在真空环境下进行单点能量的计算，才可以与sobEDAw进行能量分解对应起来。如果需要的话是否还要加conterpoise=2计算单点能。

sobereva

对于使用下文介绍的JCC文章里我提出的经验公式，原文在什么条件下计算（真空）你就必须在什么条件下计算，否则没法套用那公式
透彻认识氢键本质、简单可靠地估计氢键强度：一篇2019年JCC上的重要研究文章介绍
http://sobereva.com/513（http://bbs.keinsci.com/thread-14600-1-1.html）

溶剂效应对氢键强度的影响可以自己手动计算，参考下文
谈谈分子间结合能的构成以及分解分析思想
http://sobereva.com/733（http://bbs.keinsci.com/thread-50890-1-1.html）

是

用“conterpoise=2”的意义不明，绝对不要把使用counterpoise当成计算相互作用能的必须方式。更何况当基组足够大的时候BSSE都小到可以忽略

yangna11

sobereva 发表于 2026-2-26 17:07
对于使用下文介绍的JCC文章里我提出的经验公式，原文在什么条件下计算（真空）你就必须在什么条件下计算， ...

好的，谢谢老师，还想请教一下老师，JCC文章中结构优化的用到的泛函和基组是B3LYP-D3(BJ)/def-TZVP；单点任务用到的泛函和基组B3LYP-D3(BJ)/ma-TZVPP，请问是这样的吗，如果是的话，那么我用sobEDAw进行能量分解分析的时候也要用同样的条件对结构进行优化和单点能的计算。因为我之前结构优化的条件是b3lyp/6-31g(d,p)，单点任务是b3lyp/6-311+g(d,p)。是否需要重新进行sobEDAw进行能量分解分析

sobereva

yangna11 发表于 2026-2-26 17:32
好的，谢谢老师，还想请教一下老师，JCC文章中结构优化的用到的泛函和基组是B3LYP-D3(BJ)/def-TZVP；单点 ...

不需要，不是同一件事