24小时热门版块排行榜    

查看: 2534  |  回复: 4

coolrainbow

木虫 (著名写手)

未来国家冻凉

[交流] 一个Huckel分子轨道理论的计算软件(程序+源代码) 已有4人参与

程序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

就可以得到如下的图像:

一个Huckel分子轨道理论的计算软件(程序+源代码)
回复此楼
技术博客:http://hi.baidu.com/coolrainbow/blog
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
2楼2017-08-02 10:56:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tcclab

银虫 (正式写手)

不错!顶!
3楼2017-08-03 15:28:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
谢谢分享
4楼2017-08-04 20:51:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wywzjzz

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
赞!不知道遇到相同本征值的时候,您怎么让本征向量看着更舒服一点的~
5楼2019-08-13 18:23:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 coolrainbow 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 接收2026硕士调剂(学硕+专硕) +4 allen-yin 2026-03-23 6/300 2026-03-23 15:04 by 汪!?!
[考研] 291求调剂 +5 孅華 2026-03-22 5/250 2026-03-23 09:20 by haoshis
[考研] 289材料与化工(085600)B区求调剂 +3 这么名字咋样 2026-03-22 4/200 2026-03-22 17:56 by 云民大李老师
[考研] 求调剂 +5 Zhangbod 2026-03-21 7/350 2026-03-22 13:13 by Zhangbod
[考研] 303求调剂 +5 安忆灵 2026-03-22 6/300 2026-03-22 12:46 by 素颜倾城1988
[考研] 考研调剂 +4 来好运来来来 2026-03-21 4/200 2026-03-22 12:15 by 星空星月
[基金申请] 山东省面上项目限额评审 +4 石瑞0426 2026-03-19 4/200 2026-03-22 08:50 by Wei_ren
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 求调剂 +3 .m.. 2026-03-21 4/200 2026-03-21 16:25 by barlinike
[考研] 265求调剂 +12 梁梁校校 2026-03-19 14/700 2026-03-21 13:38 by lature00
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 265求调剂 +3 Jack?k?y 2026-03-17 3/150 2026-03-21 03:17 by JourneyLucky
[考研] 304求调剂 +6 曼殊2266 2026-03-18 6/300 2026-03-21 00:32 by JourneyLucky
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 一志愿苏州大学材料求调剂,总分315(英一) +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] 材料学求调剂 +4 Stella_Yao 2026-03-20 4/200 2026-03-20 20:28 by ms629
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
信息提示
请填处理意见