计算化学公社

标题: 如何从Gaussian产生的fchk文件计算分子轨道cube文件中的格点数据 [打印本页]

作者
Author: 让你变成回忆 时间: 2018-10-13 18:17
标题: 如何从Gaussian产生的fchk文件计算分子轨道cube文件中的格点数据
本帖最后由让你变成回忆于 2018-10-13 23:07 编辑

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

作者
Author: granvia 时间: 2018-10-13 19:34
同问

作者
Author: sobereva 时间: 2018-10-14 00:26
把这两个看了应该就差不多了
高斯fch文件与wfn波函数文件的介绍及转换方法
http://sobereva.com/55
利用wfn文件计算电子密度的代码的编写方法
http://sobereva.com/182

作者
Author: 让你变成回忆 时间: 2018-10-14 12:45

sobereva 发表于 2018-10-14 00:26
把这两个看了应该就差不多了
高斯fch文件与wfn波函数文件的介绍及转换方法
http://sobereva.com/55

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

作者
Author: 让你变成回忆 时间: 2018-10-16 20:58
本帖最后由让你变成回忆于 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的代码，让它打印出自己需要的，再和自己写的对比。

作者
Author: sobereva 时间: 2018-10-17 04:22

让你变成回忆发表于 2018-10-16 20:58
sob老师您好，您所提到的这两个博文我基本上看完了。
我现在的理解是：
您在《利用wfn文件计算电子密度 ...

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

cube格式介绍参见
Gaussian型cube文件简介及读、写方法和简单应用
http://sobereva.com/125

作者
Author: 让你变成回忆 时间: 2018-10-17 08:08

sobereva 发表于 2018-10-17 04:22
就是循环cube文件定义的每个格点，对每个格点位置调用计算密度的函数来计算

cube格式介绍参见

嗯嗯。谢谢sob老师。

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