Python (Numpy) による Hückel 分子軌道法計算

Python はプログラミング言語ですが、豊富なライブラリを持ち、スプリプト言語のような感覚で高度な科学技術計算を簡単に行うことができます。

※　PythonによるHückel分子軌道法計算プログラムについては、簡単なヒュッケル分子軌道法（HMO法）計算プログラム のページを参照してください。

1. Numpyモジュールを読み込みます。

※　IPython(x,y)などのように、起動時にNumpyやScipyが読み込まれる設定になっている環境ではimportの必要はありません。

Python処理系を単独でインストールしている場合は、Numpyモジュールを追加インストールする必要があるかもしれません。

実行例　エチレン　　　C1＝C2

=  α + β       + )
=  α - β       - )

と求められました。

実行例　ブタジエン　　　C1＝C2-C3＝C4

一般には、eig()で求められた固有値は大きさの順に並んでいないため、降順に並べ替える必要があります。

　　 =  α + 1.618β   　
　　 =  α + 0.618β     
　　 =  α  - 0.618β   
　　 =  α - 1.618β    　　

　と求められました。

実行例3　ホルムアルデヒド　　　 　 　C1＝O2

=  α + 1.618 β     ψ₁ = 0.526φ₁ + 0.850φ₂
=  α - 0.618 β     ψ₂ = 0.850φ₁ - 0.526φ₂

と求められました。

---