24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 2567  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 271求调剂 +30 2261744733 2026-04-11 34/1700 2026-04-13 15:24 by Kilig0317
[考研] 化工学硕 285求调剂 +29 Wisjxn 2026-04-07 29/1450 2026-04-13 13:47 by floriea
[教师之家] 山东双非院校考核超级无底线,领导幸灾乐祸,教师遭殃恐 +3 qut2026 2026-04-11 7/350 2026-04-12 20:24 by qut2026
[考研] 电子信息270求调剂 +18 terminal469 2026-04-07 18/900 2026-04-12 16:23 by ajpv风雷
[考研] 085410 273求调剂 +10 X1999 2026-04-09 10/500 2026-04-12 09:24 by 逆水乘风
[考研] 22专硕求调剂 +6 haoyun上岸 2026-04-11 8/400 2026-04-11 23:21 by labixiaoqiao
[考研] 求调剂 +11 翩翩一书生 2026-04-09 11/550 2026-04-11 19:57 by 逆水乘风
[考研] 298求调剂 +9 钉叮咚冬瓜 2026-04-07 11/550 2026-04-11 09:35 by zhq0425
[考研] 0854调剂 +4 长弓傲 2026-04-09 4/200 2026-04-11 09:18 by 猪会飞
[考研] 284求调剂 +12 archer.. 2026-04-10 13/650 2026-04-11 08:44 by zhq0425
[考研] 化学工程与技术324调剂 +23 孙常华 2026-04-09 25/1250 2026-04-11 00:07 by 骑牛渡寒江
[考研] 368求调剂 +3 17385968172 2026-04-10 3/150 2026-04-10 20:12 by 电子奥德彪
[考研] 085601初试330分找调剂 +10 流心奶黄包l 2026-04-09 10/500 2026-04-10 08:14 by Sammy2
[考研] 材料专硕调剂 +16 哈哈哈吼吼吼哈 2026-04-07 17/850 2026-04-09 21:16 by wutongshun
[考研] 070300化学 求调剂 +13 73372112 2026-04-08 13/650 2026-04-09 20:22 by maddjdld
[考研] 调剂 +19 2261744733 2026-04-08 19/950 2026-04-09 19:11 by vgtyfty
[考研] 生物学328分求调剂 +9 闪电kkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
[考研] 264求调剂 +11 麦小叮当 2026-04-07 11/550 2026-04-08 16:05 by 一只好果子?
[考研] 388求调剂 +6 四川王涛 2026-04-07 8/400 2026-04-08 00:17 by JourneyLucky
[考研] 22408 318分求调剂 +4 勤奋的小笼包 2026-04-06 6/300 2026-04-07 15:05 by 纸鹤555
信息提示
请填处理意见