计算化学公社

标题: 一个Huckel分子轨道理论的计算软件（程序+源代码） [打印本页]

作者
Author: coolrainbow 时间: 2017-8-2 09:46
标题: 一个Huckel分子轨道理论的计算软件（程序+源代码）
程序HMO，给出分子的坐标和电荷后，就可以对分子进行HMO计算，输出电荷密度，键级，和能量，还可以给出相应的分子数据图像，十分方便。特别适合结构化学教学。
二进制程序：Windows和Linux的程序可以从此下载：http://www.zhjun-sci.com/software-others.php
源代码：此程序全部源代码公开，可以在Github下载：https://github.com/ZhangJunQCC/H ... Orbital-Calculation

如对于二苯乙烯，输入文件如下styrene.inp：

1.8 0
C                -3.33505854 -4.05509747 0.03172776
C                -1.93989854 -4.05509747 0.03172776
C                -1.24236054 -2.84734647 0.03172776
C                -1.94001454 -1.63883747 0.03052876
C                -3.33483954 -1.63891547 0.03004976
C                -4.03244054 -2.84712147 0.03104576
C                -1.16951047 -0.30544951 0.03044614
C                -1.84201131 0.87111682 0.02928423
C                -1.06360352 2.19990338 0.03203690
C                -1.75320092 3.41233984 0.03084857
C                0.33149892 2.19089168 0.03626174
C                -1.04780409 4.61603387 0.03357108
C                1.03679478 3.39425768 0.03850471
C                0.34687625 4.60719876 0.03679910

第一行的1.8，表示距离小于1.8 的两个原子都算成键；
第一行的0，表示原子的电荷为0.当然也可以是+1，-2 之类。
敲入一个简单的命令；

hmo styrene.inp > styrene.out

就可以得到一个out文件，里面可以看到轨道、能级、键级等：

MO Indices:             1             2             3
  Occupation:             2             2             2
Eigenvalues:    3.21050865    3.06407743    2.50466736
Energies:    a+2.21051b    a+2.06408b    a+1.50467b
         1    0.18305145    -0.25721905    0.36874498
         2    0.20231841    -0.26546002    0.27741927
         3    0.26417514    -0.29071098    0.04867874
         4    0.38164303    -0.33458996    -0.20417396
...
Bond orders:
         #          Bond order
      1-1          1.00000000
      1-2          0.65658339
      1-3          0.00000000
      1-4       -0.29424047
      1-5          0.00000000
      1-6          0.65658339
...
Atomic Properties:
         #    Electron Density          Pi Valence
         1          1.00000000          1.31316679
         2          1.00000000          1.33775485
         3          1.00000000          1.28257491
         4          1.00000000          1.63390015
         5          1.00000000          1.28257491

最后还得到一个文件叫做styrene.dot，它可以用Graphviz渲染成可视化的图像。只要安装了Graphviz，就可以调用如下的命令：

sfdp styrene.dot -Tpdf -o styrene.pdf

就可以得到如下的图像：

(, 下载次数 Times of downloads: 22)

作者
Author: Warm_Cloud 时间: 2017-8-3 11:21
张老师出品

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