24小时热门版块排行榜    

查看: 2533  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 350求调剂 +6 weudhdk 2026-03-19 6/300 2026-03-23 15:47 by tangyuan0840221
[考研] 263求调剂 +6 yqdszhdap- 2026-03-22 9/450 2026-03-23 12:57 by yqdszhdap-
[考研] 一志愿070300浙大化学358分,求调剂! +4 酥酥鱼.. 2026-03-21 4/200 2026-03-23 08:12 by Iveryant
[考研] 307求调剂 +11 冷笙123 2026-03-17 11/550 2026-03-22 20:16 by edmund7
[考研] 287求调剂 +8 晨昏线与星海 2026-03-19 9/450 2026-03-22 17:01 by i_cooler
[考研] 一志愿北京化工大学070300 学硕336求调剂 +5 vv迷 2026-03-21 8/400 2026-03-22 14:20 by ColorlessPI
[基金申请] 山东省面上项目限额评审 +4 石瑞0426 2026-03-19 4/200 2026-03-22 08:50 by Wei_ren
[考研] 278求调剂 +9 烟火先于春 2026-03-17 9/450 2026-03-21 17:47 by 学员8dgXkO
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 265求调剂 +3 Jack?k?y 2026-03-17 3/150 2026-03-21 03:17 by JourneyLucky
[考研] 301求调剂 +10 yy要上岸呀 2026-03-17 10/500 2026-03-21 03:14 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 321求调剂 +9 何润采123 2026-03-18 11/550 2026-03-20 23:19 by JourneyLucky
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +22 rare12345 2026-03-18 22/1100 2026-03-20 20:39 by zhukairuo
[考研] 085410人工智能专硕317求调剂(0854都可以) +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by 不168
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
信息提示
请填处理意见