granvia 发表于 2019-1-10 12:53
没算过这类问题。是否可以用高斯的溶剂模型选项来把你的介电函数加进去呢？

它是一个函数，随着r变化，不是常量，如果考虑用溶剂化效应加进去，不知道应该怎么加呢？

DUNING 发表于 2019-1-10 10:01
对的，本来就有，文献说的是这个体系得考虑介电屏蔽效应，所以有了这一项

没算过这类问题。是否可以用高斯的溶剂模型选项来把你的介电函数加进去呢？

granvia 发表于 2019-1-9 22:54
你的哈密顿算符里本来就有了介电常函数项吗？ 另外，四个态的积分很正常。双电子积分就是对四个基函数的 ...

对的，本来就有，文献说的是这个体系得考虑介电屏蔽效应，所以有了这一项

DUNING 发表于 2019-1-9 17:18
谢谢老师，我刚刚查了资料，发现双电子积分中的库伦积分和我看到文献中的这个公式还是有不一样的地方，文 ...

你的哈密顿算符里本来就有了介电常函数项吗？ 另外，四个态的积分很正常。双电子积分就是对四个基函数的积分，然后做个线性代数就得到你想要的四个态的积分。

想解析地计算那个积分很不容易，你得自己折腾双电子积分代码，这是编写量化程序最复杂的部分之一
如果通过立方格点以数值方法积分倒还比较简单。把fch文件载入弄明白，之后再把格点数据这块搞明白就差不多了，参考
Gaussian型cube文件简介及读、写方法和简单应用
http://sobereva.com/125

liyuanhe211 发表于 2019-1-9 16:39

谢谢老师，我刚刚查了资料，发现双电子积分中的库伦积分和我看到文献中的这个公式还是有不一样的地方，文献中的公式在分母中有介电函数项，并且文献中的公式好像涉及到四个态（该计算主要针对俄歇过程）。不知道我的理解对不对？

DUNING 发表于 2019-1-9 16:34
老师，你好，我是想要这个数

liyuanhe211 发表于 2019-1-9 16:05
不知道你是就想要这个数，还是为了做练习。如果就想要这个数，计算过程中的库伦积分高斯好像直接能输出

老师，你好，我是想要这个数

DUNING 发表于 2019-1-9 14:31
谢谢老师的指正，现在已经改了，确实是我太大意。其实我就是在一篇文献中看到了一个库伦积分的公式，里面 ...

不知道你是就想要这个数，还是为了做练习。如果就想要这个数，计算过程中的库伦积分高斯好像直接能输出

还不如直接说你要干什么，没准Multiwfn就能做

高斯fch文件与wfn波函数文件的介绍及转换方法
http://sobereva.com/55

Gaussian都拼不对，我怀疑自行开发程序的可能性