24小时热门版块排行榜    

查看: 2834  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 5/250 2026-08-18 11:25 by 無關想念
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +43 医学老男孩 2026-08-13 97/4850 2026-08-18 11:07 by 医学老男孩
[基金申请] 时间戳变了,能看出什么问题? +5 基诺咪客 2026-08-17 6/300 2026-08-18 11:06 by hunter无悔
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +4 Tide man 2026-08-14 4/200 2026-08-18 10:57 by marinewhale
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +12 Ldrop2023 2026-08-13 15/750 2026-08-17 15:02 by 小豌豆_发芽
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
信息提示
请填处理意见