请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 161|回复: 12

[Gaussian/gview] 如何找到gaussian输出的波函数具体形式

[复制链接]

84

帖子

0

威望

240

eV
积分
324

Level 3 能力者

发表于 2019-1-9 11:15:27 | 显示全部楼层 |阅读模式
本帖最后由 DUNING 于 2019-1-9 14:38 编辑

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

1万

帖子

25

威望

2万

eV
积分
37489

管理员

公社社长

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

Gaussian都拼不对,我怀疑自行开发程序的可能性
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群1号:18616395,2号:466017436。超过4000人,用于交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

291

帖子

2

威望

1492

eV
积分
1823

Level 5 (御坂)

发表于 2019-1-9 13:42:24 | 显示全部楼层
还不如直接说你要干什么,没准Multiwfn就能做

84

帖子

0

威望

240

eV
积分
324

Level 3 能力者

 楼主| 发表于 2019-1-9 14:31:35 | 显示全部楼层
本帖最后由 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老师能帮我看看,我应该怎么做?本人目前还是个小白,所以问题也很低级,谢谢老师 1547015267(1).png

2376

帖子

23

威望

7950

eV
积分
10786

Level 6 (一方通行)

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

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

84

帖子

0

威望

240

eV
积分
324

Level 3 能力者

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

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

2376

帖子

23

威望

7950

eV
积分
10786

Level 6 (一方通行)

发表于 2019-1-9 16:39:33 | 显示全部楼层
DUNING 发表于 2019-1-9 16:34
老师,你好,我是想要这个数

snap0000782.jpg


84

帖子

0

威望

240

eV
积分
324

Level 3 能力者

 楼主| 发表于 2019-1-9 17:18:01 | 显示全部楼层

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

1万

帖子

25

威望

2万

eV
积分
37489

管理员

公社社长

发表于 2019-1-9 22:24:25 | 显示全部楼层
想解析地计算那个积分很不容易,你得自己折腾双电子积分代码,这是编写量化程序最复杂的部分之一
如果通过立方格点以数值方法积分倒还比较简单。把fch文件载入弄明白,之后再把格点数据这块搞明白就差不多了,参考
Gaussian型cube文件简介及读、写方法和简单应用
http://sobereva.com/125
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群1号:18616395,2号:466017436。超过4000人,用于交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

203

帖子

0

威望

422

eV
积分
625

Level 4 (黑子)

发表于 2019-1-9 22:54:15 来自手机 | 显示全部楼层
DUNING 发表于 2019-1-9 17:18
谢谢老师,我刚刚查了资料,发现双电子积分中的库伦积分和我看到文献中的这个公式还是有不一样的地方,文 ...

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

84

帖子

0

威望

240

eV
积分
324

Level 3 能力者

 楼主| 发表于 2019-1-10 10:01:39 | 显示全部楼层
granvia 发表于 2019-1-9 22:54
你的哈密顿算符里本来就有了介电常函数项吗? 另外,四个态的积分很正常。双电子积分就是对四个基函数的 ...

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

203

帖子

0

威望

422

eV
积分
625

Level 4 (黑子)

发表于 2019-1-10 12:53:33 来自手机 | 显示全部楼层
DUNING 发表于 2019-1-10 10:01
对的,本来就有,文献说的是这个体系得考虑介电屏蔽效应,所以有了这一项

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

84

帖子

0

威望

240

eV
积分
324

Level 3 能力者

 楼主| 发表于 2019-1-10 17:01:31 | 显示全部楼层
granvia 发表于 2019-1-10 12:53
没算过这类问题。是否可以用高斯的溶剂模型选项来把你的介电函数加进去呢?

它是一个函数,随着r变化,不是常量,如果考虑用溶剂化效应加进去,不知道应该怎么加呢?
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2019-1-20 18:00 , Processed in 0.119002 second(s), 26 queries .

快速回复 返回顶部 返回列表