静电作用能，波利排斥能和轨道相互作用用能的计算

哈哈宝

文献中写： “Eelstat   electrostatic interaction energy between the fragments that are calculated with the frozen electron density distribution of A and B in the geometry of the complex AB”
“ EPauli is calculated by enforcing the Kohn-Sham determinant of AB, which results from superimposing fragments A and B, to obey the Pauli principle  by  antisymmetrization  and  renormalization. ” "Eorbital, is calculated in the final step of the ETS analysis when the Kohn-Sham orbitals relax  to  their  optimal  form"
似乎是看懂字面的意思了，但是具体到计算还是一锅浆糊，文献中给出了具体的数值（二茂铁）我想用高斯软件，同样方法和基组重复一下这个计算，试着得到静电能，排斥能和轨道相互作用能的数值，但是究竟如何下手呢？
1. 怎么A和B部分的冻结电子密度？
2. 如何强制A和B遵循电子的反对称和归一化在高斯计算中？


求高手指点！

brothers

Gaussian并不能做能量分解分析（EDA）。可以使用ADF（申请一个30天的试用）或者Gamess（免费）。其他软件用的相对少一些。
完全不必自己一步步计算各部分的能量。与常规的单点计算相比，做能量分解只需额外指定：
（1）能量分解模型（目前有十几种分解方案）；
（2）A与B fragments各自所包含的原子数目；
（3）A与B的电荷（原子归属和电荷的设定实际上反映了EDA分析的“随意性”）
然后就可以run了（好久没做了，印象中是不缺啥了）。论坛中有Gamess做能量分解的帖子，都搜出来看一遍，相信你三天之内就能跑成功。
但是对结果的分析仍然要求具有一定的理论基础，需要多看模型和文献。
可以参考Frenking以及苏培峰的文章。google学术一下energy decomposition analysis，最前几篇的引用率都很高，也都有分析案例。下下来仔细看看就对了。
祝你好运。

liyuanhe211

brothers

对了，不要用32位windows版本的Gamess-US，当年用它算出来一片一片的NaN，用了64位的就正常了。

sobereva

前面已经回复得比较全面了，补充一点Pauli不叫波利而叫泡利

哈哈宝

非常感谢楼上各位的详细指导！