计算化学公社

 找回密码 Forget password
 注册 Register
Views: 12014|回复 Reply: 12
打印 Print 上一主题 Last thread 下一主题 Next thread

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

[复制链接 Copy URL]

183

帖子

0

威望

795

eV
积分
978

Level 4 (黑子)

本帖最后由 DUNING 于 2019-1-9 14:38 编辑

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

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

2#
发表于 Post on 2019-1-9 13:13:50 | 只看该作者 Only view this author
高斯fch文件与wfn波函数文件的介绍及转换方法
http://sobereva.com/55

Gaussian都拼不对,我怀疑自行开发程序的可能性
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

490

帖子

2

威望

5306

eV
积分
5836

Level 6 (一方通行)

3#
发表于 Post on 2019-1-9 13:42:24 | 只看该作者 Only view this author
还不如直接说你要干什么,没准Multiwfn就能做

183

帖子

0

威望

795

eV
积分
978

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2019-1-9 14:31:35 | 只看该作者 Only view this author
本帖最后由 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老师能帮我看看,我应该怎么做?本人目前还是个小白,所以问题也很低级,谢谢老师

3097

帖子

29

威望

1万

eV
积分
17221

Level 6 (一方通行)

5#
发表于 Post on 2019-1-9 16:05:45 | 只看该作者 Only view this author
DUNING 发表于 2019-1-9 14:31
谢谢老师的指正,现在已经改了,确实是我太大意。其实我就是在一篇文献中看到了一个库伦积分的公式,里面 ...

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

183

帖子

0

威望

795

eV
积分
978

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2019-1-9 16:34:24 | 只看该作者 Only view this author
liyuanhe211 发表于 2019-1-9 16:05
不知道你是就想要这个数,还是为了做练习。如果就想要这个数,计算过程中的库伦积分高斯好像直接能输出

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

3097

帖子

29

威望

1万

eV
积分
17221

Level 6 (一方通行)

7#
发表于 Post on 2019-1-9 16:39:33 | 只看该作者 Only view this author
DUNING 发表于 2019-1-9 16:34
老师,你好,我是想要这个数



183

帖子

0

威望

795

eV
积分
978

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2019-1-9 17:18:01 | 只看该作者 Only view this author

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

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

9#
发表于 Post on 2019-1-9 22:24:25 | 只看该作者 Only view this author
想解析地计算那个积分很不容易,你得自己折腾双电子积分代码,这是编写量化程序最复杂的部分之一
如果通过立方格点以数值方法积分倒还比较简单。把fch文件载入弄明白,之后再把格点数据这块搞明白就差不多了,参考
Gaussian型cube文件简介及读、写方法和简单应用
http://sobereva.com/125
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

1043

帖子

0

威望

4188

eV
积分
5231

Level 6 (一方通行)

10#
发表于 Post on 2019-1-9 22:54:15 来自手机 | 只看该作者 Only view this author
DUNING 发表于 2019-1-9 17:18
谢谢老师,我刚刚查了资料,发现双电子积分中的库伦积分和我看到文献中的这个公式还是有不一样的地方,文 ...

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

183

帖子

0

威望

795

eV
积分
978

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2019-1-10 10:01:39 | 只看该作者 Only view this author
granvia 发表于 2019-1-9 22:54
你的哈密顿算符里本来就有了介电常函数项吗? 另外,四个态的积分很正常。双电子积分就是对四个基函数的 ...

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

1043

帖子

0

威望

4188

eV
积分
5231

Level 6 (一方通行)

12#
发表于 Post on 2019-1-10 12:53:33 来自手机 | 只看该作者 Only view this author
DUNING 发表于 2019-1-10 10:01
对的,本来就有,文献说的是这个体系得考虑介电屏蔽效应,所以有了这一项

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

183

帖子

0

威望

795

eV
积分
978

Level 4 (黑子)

13#
 楼主 Author| 发表于 Post on 2019-1-10 17:01:31 | 只看该作者 Only view this author
granvia 发表于 2019-1-10 12:53
没算过这类问题。是否可以用高斯的溶剂模型选项来把你的介电函数加进去呢?

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

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-24 10:15 , Processed in 0.177533 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list