如何得到静电势在原子核处扣除原子核贡献的具体数值？

dxyan

看到文献说供体化合物的EPH.nuc值与该化合物的氢键酸性有很好的关联性，作者进而比较了不同化合物的酚羟基中氢原子核处扣除原子核贡献的静电势（EPH.nuc），不知道这个数值从哪里得到啊？什么是氢原子核处扣除原子核贡献的静电势呢？
问题较多，请求达人赐教，谢谢？

sobereva

静电势是由电子和原子核电荷一起贡献的，直接计算原子核处的静电势值显然不可行，因为会得到无限大的结果。如果在氢原子核处计算扣除了这个原子核贡献的静电势值，就可以衡量这个原子核（即质子）在分子中感受到的静电作用的大小，显然静电相互作用越强质子越难电离，故与酸性有相关性。

用Multiwfn可方便地计算。比如你要计算体系中某个氢原子（假设为5号原子）的原子核处扣除了其核电荷贡献的静电势值，就启动Multiwfn，载入体系波函数文件（.wfn、.wfx、.fch、.molden皆可），然后输入
1000
2  //选择用户自定义函数
39  //第39号用户自定义函数，对应于扣除某原子核贡献的静电势值，详见手册3.7节第39个函数的说明
1000
3  //选择扣除哪个原子的原子核贡献
5  //5号原子
1  //计算某个点上各种实空间函数的数值
a5   //计算5号原子核处的函数数值
输出文件中的User defined real space function后面的值就是你要的量，单位为a.u.。这个数值越负，说明这个质子与体系静电相互作用越强，酸性通常越弱，质子越难电离走。

dxyan

sobereva 发表于 2016-2-23 22:31
静电势是由电子和原子核电荷一起贡献的，直接计算原子核处的静电势值显然不可行，因为会得到无限大的结果。 ...

谢谢大神回复，我先自己试试。

快要学会了吧

你好，请问有没有相关文献可以分享一下呢~