计算化学公社

标题: 静电作用能,波利排斥能和轨道相互作用用能的计算 [打印本页]

作者
Author:
哈哈宝    时间: 2015-11-17 17:04
标题: 静电作用能,波利排斥能和轨道相互作用用能的计算
文献中写: “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遵循电子的反对称和归一化在高斯计算中?


求高手指点!

作者
Author:
brothers    时间: 2015-11-17 20:14
本帖最后由 brothers 于 2015-11-17 20:26 编辑

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

作者
Author:
liyuanhe211    时间: 2015-11-17 20:29
(, 下载次数 Times of downloads: 30)



作者
Author:
brothers    时间: 2015-11-17 20:39
对了,不要用32位windows版本的Gamess-US,当年用它算出来一片一片的NaN,用了64位的就正常了。
作者
Author:
sobereva    时间: 2015-11-18 02:33
前面已经回复得比较全面了,补充一点Pauli不叫波利而叫泡利
作者
Author:
哈哈宝    时间: 2015-11-18 08:06
非常感谢楼上各位的详细指导!




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3