计算化学公社

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

[波函数分析求助] 如何从Gaussian产生的fchk文件计算分子轨道cube文件中的格点数据

[复制链接 Copy URL]

490

帖子

2

威望

5126

eV
积分
5656

Level 6 (一方通行)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 让你变成回忆 于 2018-10-13 23:07 编辑

如题,最近想了解一下如何从fchk文件的GTF出发,计算得到分子轨道cub文件中的格点数据。在Multiwfn的手册的主功能200的子功能3中没有列出具体的公式,请教sob老师有推荐的资料吗?谢谢!最近自己想写程序实现这个功能,了解一下这个过程。

1043

帖子

0

威望

4174

eV
积分
5217

Level 6 (一方通行)

2#
发表于 Post on 2018-10-13 19:34:29 来自手机 | 只看该作者 Only view this author
同问

6万

帖子

99

威望

5万

eV
积分
120148

管理员

公社社长

3#
发表于 Post on 2018-10-14 00:26:32 | 只看该作者 Only view this author
把这两个看了应该就差不多了
高斯fch文件与wfn波函数文件的介绍及转换方法
http://sobereva.com/55
利用wfn文件计算电子密度的代码的编写方法
http://sobereva.com/182
北京科音自然科学研究中心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

威望

5126

eV
积分
5656

Level 6 (一方通行)

4#
 楼主 Author| 发表于 Post on 2018-10-14 12:45:45 | 只看该作者 Only view this author
sobereva 发表于 2018-10-14 00:26
把这两个看了应该就差不多了
高斯fch文件与wfn波函数文件的介绍及转换方法
http://sobereva.com/55

好的。我去仔细读读。
谢谢sob老师!

490

帖子

2

威望

5126

eV
积分
5656

Level 6 (一方通行)

5#
 楼主 Author| 发表于 Post on 2018-10-16 20:58:28 | 只看该作者 Only view this author
本帖最后由 让你变成回忆 于 2018-10-16 21:01 编辑
sobereva 发表于 2018-10-14 00:26
把这两个看了应该就差不多了
高斯fch文件与wfn波函数文件的介绍及转换方法
http://sobereva.com/55

sob老师您好,您所提到的这两个博文我基本上看完了。
我现在的理解是:
您在《利用wfn文件计算电子密度的代码的编写方法》中需要手动输入的x,y,z这个点,其实就是cub中预先根据格点质量划分好的一些点,然后计算这些点处轨道波函数的数值即可。
请问是这个过程吗?
PS:学习各种程序产生的各种文件的内容和格式,Multiwfn的源代码真的是好东西。这次总算是把fch文件和wfn文件搞清楚了,在fch2wfn转换的那几个循环绕了许久。debug自己的代码一个很好的方法就是修改Multiwfn的代码,让它打印出自己需要的,再和自己写的对比。

6万

帖子

99

威望

5万

eV
积分
120148

管理员

公社社长

6#
发表于 Post on 2018-10-17 04:22:13 | 只看该作者 Only view this author
让你变成回忆 发表于 2018-10-16 20:58
sob老师您好,您所提到的这两个博文我基本上看完了。
我现在的理解是:
您在《利用wfn文件计算电子密度 ...

就是循环cube文件定义的每个格点,对每个格点位置调用计算密度的函数来计算

cube格式介绍参见
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

490

帖子

2

威望

5126

eV
积分
5656

Level 6 (一方通行)

7#
 楼主 Author| 发表于 Post on 2018-10-17 08:08:01 | 只看该作者 Only view this author
sobereva 发表于 2018-10-17 04:22
就是循环cube文件定义的每个格点,对每个格点位置调用计算密度的函数来计算

cube格式介绍参见

嗯嗯。谢谢sob老师。

本版积分规则 Credits rule

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

GMT+8, 2025-8-16 02:11 , Processed in 0.142181 second(s), 20 queries , Gzip On.

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