求助：利用氢键的键临界点(BCP)位置的属性，批量计算分子间氢键键能的方法

吃西瓜的佩奇

卢老师您好！我在利用帖子http://sobereva.com/522中基于键临界点的电子密度估计内氢键键能的方法，计算水分子和另一个分子形成的二聚体中O-H···O氢键的键能，由于需要计算的模型有上千个，所以我需要写脚本并调用Multiwfn程序批量计算。但在识别对应内氢键的BCP的序号时每一个研究模型的序号并不统一。所以想向您请教怎么才能解决这个问题？我目前想到的是识别形成分子间氢键的O-H···O原子的序号，指定原子进行计算，但是不知如何实现，特来请教您是否可行，或者您有没有更好的办法可以帮助我解决这一难题，谢谢您！

sobereva

我不知道你的这些体系的原子有什么规律没有。如果每个体系的形成你要考察的那个氢键的原子序号都是统一的，可以在搜索时用主功能2里的1 Search CPs from given starting points，然后选2 Using an atom or midpoint of two atoms as starting point从而使用指定的两个原子的中点作为搜索临界点的初猜位置，这样就能搜索到对应那个氢键的BCP，而且这是1号CP。

例如