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

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

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对应上。老师您看有哪些地方可能会出现问题的地方吗，麻烦您了

描述得太不具体，没法说

开源免费的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