计算化学公社

标题: 如何找到gaussian输出的波函数具体形式 [打印本页]

作者
Author: DUNING 时间: 2019-1-9 11:15
标题: 如何找到gaussian输出的波函数具体形式
本帖最后由 DUNING 于 2019-1-9 14:38 编辑

在得到gaussian的输出结果后，还想编写辅助程序进行后续分析，其中需要某些能态上的波函数的具体形式，请问这时应该在哪儿去找？我看了wfn文件中记录的数据，其中给出了所有原始高斯函数的指数，我怎么能把它对应到我要的几个能态上呢？

作者
Author: sobereva 时间: 2019-1-9 13:13
高斯fch文件与wfn波函数文件的介绍及转换方法
http://sobereva.com/55

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

作者
Author: 让你变成回忆 时间: 2019-1-9 13:42
还不如直接说你要干什么，没准Multiwfn就能做

作者
Author: DUNING 时间: 2019-1-9 14:31
本帖最后由 DUNING 于 2019-1-9 14:37 编辑

sobereva 发表于 2019-1-9 13:13
高斯fch文件与wfn波函数文件的介绍及转换方法
http://sobereva.com/55

谢谢老师的指正，现在已经改了，确实是我太大意。其实我就是在一篇文献中看到了一个库伦积分的公式，里面需要用到波函数，公式不是很难，就想着自己试试把这个积分算出来。C:\Users\Administrator\Desktop\1547015267(1).png这个公式贴在下面了，
希望Sob老师能帮我看看，我应该怎么做？本人目前还是个小白，所以问题也很低级，谢谢老师 (, 下载次数 Times of downloads: 78)

作者
Author: liyuanhe211 时间: 2019-1-9 16:05

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

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

作者
Author: DUNING 时间: 2019-1-9 16:34

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

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

作者
Author: liyuanhe211 时间: 2019-1-9 16:39

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

(, 下载次数 Times of downloads: 96)

作者
Author: DUNING 时间: 2019-1-9 17:18

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

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

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

作者
Author: granvia 时间: 2019-1-9 22:54

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

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

作者
Author: DUNING 时间: 2019-1-10 10:01

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

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

作者
Author: granvia 时间: 2019-1-10 12:53

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

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

作者
Author: DUNING 时间: 2019-1-10 17:01

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)