求助分子轨道绘制问题

小白努力学习 · 发表于 Post on 2024-10-26 20:17:02

本帖最后由小白努力学习于 2024-10-26 20:18 编辑

求救各位老师一个问题：我想自己绘制分子轨道，目前利用Gaussian转换的fch文件读取参数进行组合，原始高斯基函数用的是下面图片的公式。但是绘制的和GaussView不同，大概是什么原因造成的？

sobereva · 发表于 Post on 2024-10-26 23:00:48

描述得太不具体，没法说

开源免费的Multiwfn是看分子轨道非常好用的程序
使用Multiwfn观看分子轨道
http://sobereva.com/269（http://bbs.keinsci.com/thread-462-1-1.html）

可以看Multiwfn的function.f90里的function fmo函数了解轨道波函数是怎么计算的。Multiwfn主页的Resources页面里有Modifying source code of Multiwfn幻灯片可以帮助了解Multiwfn的数据结构。载入fch获取波函数信息的方式可参看Multiwfn的fileIO.f90里的subroutine readfch以及下文
高斯fch文件与wfn波函数文件的介绍及转换方法
http://sobereva.com/55

小白努力学习 · 发表于 Post on 2024-10-27 10:24:29

sobereva 发表于 2024-10-26 23:00
描述得太不具体，没法说

开源免费的Multiwfn是看分子轨道非常好用的程序

老师有需求需要自己写一个可视化程序。我是读取fch的Primitive exponents ，Contraction coefficients ，P(S=P) Contraction coefficients ，Alpha MO coefficients 。
按照fch文件的数据顺序进行原子轨道组合，组合出的原子轨道数目和fch文件的也对应上了。首先归一化球谐函数的高斯函数，后带入X,Y,Z格点的数据，计算R,THETA,PHI，带入归一化后的高斯函数进行计算，组合原子轨道，在带入分子轨道系数组合分子轨道。但是P轨道的可视化一直是错误的，是一个半球形，无法和GaussView对应上。老师您看有哪些地方可能会出现问题的地方吗，麻烦您了

wzkchem5 · 发表于 Post on 2024-10-27 14:45:30

小白努力学习发表于 2024-10-27 03:24
老师有需求需要自己写一个可视化程序。我是读取fch的Primitive exponents ，Contraction coefficients ， ...

是不是只作了轨道值>0的等值面，没有作轨道值<0的等值面？
我以前写过一个开源的轨道可视化程序http://bbs.keinsci.com/forum.php ... 1&fromuid=21811，可以参考一下

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[建模与可视化] 求助分子轨道绘制问题