读了http://sobereva.com/435，里面说：“注意事项2：cubegen计算静电势是基于.fch文件里的密度矩阵实现的。有的时候.fch文件里有多种密度矩阵，默认情况下使用SCF密度矩阵。如果你做的是后HF、TDDFT等计算，为了计算后HF波函数或激发态波函数的静电势，你需要修改settings.ini里的cubegendenstype参数成为对应的密度矩阵标识。比如，你用# MP2/cc-pVTZ density关键词产生了.fch文件，那么里面既有SCF密度矩阵也有MP2密度矩阵，如果你不改cubegendenstype参数，那么利用cubegen算的静电势将是HF级别的；如果你把cubegendenstype参数改为MP2，则cubegen算的静电势将是MP2级别的。更多关于.fch文件里密度矩阵的信息见《在Multiwfn中基于fch产生自然轨道的方法与激发态波函数、自旋自然轨道分析实例》（http://sobereva.com/403）文中的说明。”
也许我需要把settings.ini文件里面的cubegendenstype= SCF改一下，我在00802.fchk里面搜索到了包含SCF density, CI density、CI Rho(1) density。但还不知道对于我的CAM-B3LYP方法下基态和激发态都改成什么。

liaorongbao 发表于 2021-8-1 19:41
下边有个分子结构图，在这个位置重新做了个案例分子，分子式CH3(C2H2)3Cl计算方法基组和溶剂#P CAM-B3LYP ...

MPI是分子静电势对分子表面积分得到，坐标为主要变量。“当原子核位置固定后每一种电子云分布都得到一个特定的偶极矩，每一种电子云分布又同时对应一个特定的 MPI”后半句的这个理解夸大了电子云分布对MPI数值的影响程度，MPI的定义与计算在上面的博文以及Multiwfn手册3.15都有详述，以及需要注意的是，MPI是基于分子表面（而可能不是原子极近距离如1、2埃内）的计算，我个人理解为电子云分布变化影响在你目前所选分子表面上微乎其微，可以尝试着改变分子表面的选择（但我只用过范德华半径下的分子表面，别的情况合理性未知）。

liaorongbao 发表于 2021-8-1 12:41
下边有个分子结构图，在这个位置重新做了个案例分子，分子式CH3(C2H2)3Cl计算方法基组和溶剂#P CAM-B3LYP ...

是不是用错了波函数，即使对激发态计算也用了基态波函数算的MPI？

hdhxx123 发表于 2021-7-31 17:07
http://sobereva.com/518看文中MPI的定义，是基于静电势的计算，体现的是体系表面电荷分布不均匀性。因此与 ...

谢谢！感觉应该是您说的这个思路

http://sobereva.com/518看文中MPI的定义，是基于静电势的计算，体现的是体系表面电荷分布不均匀性。因此与结构有关而与偶极矩关系不大。上面这篇文章也评价了偶极矩的理论上表现的物理意义与其使用范围，建议仔细研读。