请问各位，我这样做是否合理？

yjcmwgk

现在我们面对的是一个这样的体系。两个负2价的酞菁环，夹着一个正3价的Tb原子。
我想考察配位之后，两个酞菁环的静电势场是啥样。
我打算做法如下：
首先，优化整个配合物分子，拿到它的波函数文件。这个分子必然是开壳层分子，所以我全部采用开壳层的UDFT方法来计算。
第二步，将这个分子视为两部分，分别命名为FragmentPc（两个酞菁环）和FragmentTb（就是Tb这一个原子）
第三步，衡量每一个占据轨道的分布，也就是每一个占据轨道分别在FragmentPc和FragmentTb上有多少布居百分比。
第四步，如果某一个占据轨道在FragmentPc上的布居为n（必然的是 0 < n < 1），那么，这个分子轨道的占据数我就直接从1改为n。
第五步，删除FragmentTb所属的全部基函数。
第六步，在改了轨道占据数，也改了基函之后，重新得到新的波函数文件。然后用这个文件去做静电势计算。
请问这样做是否可行？

sobereva

Multiwfn里直接删除Tb的高斯函数就行了（主功能-4），没必要再调节占据数，否则最终效果就double-counting了。
另外，之后在Multiwfn里计算静电势之前，注意还得把Tb的原子核对静电势的贡献扣除。需要用自定义函数39，并在菜单里设定扣除哪个原子核的贡献，手册4.1.2的例子用到了这一点可参考。

yjcmwgk

sobereva 发表于 2017-4-12 11:11
Multiwfn里直接删除Tb的高斯函数就行了（主功能-4），没必要再调节占据数，否则最终效果就double-counting ...

非常感谢！原来sob君已经考虑了占据数问题了！再次感谢！

sobereva

yjcmwgk 发表于 2017-4-12 11:37
非常感谢！原来sob君已经考虑了占据数问题了！再次感谢！

其实不是明确地去考虑了这个。
删除高斯函数就像切蛋糕，切掉1/4，剩下的自然而然就是3/4了。