Juliaによる Hückel 分子軌道法計算

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

1. 線形代数の操作を行う前に、LinearAlgebraモジュールを読み込みます。
2. 固有値を求める関数はeigvals、固有ベクトルを求める関数はeigvecsです。eigenで、固有値と固有ベクトルを同時に求めることもできます。

実行例　エチレン　　　C1＝C2

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

と求められました。

降順にしたい場合は、reverse()を用います。

詳しくは、簡単なヒュッケル分子軌道法計算プログラムのダウンロードに掲載しているJuliaプログラムを見てください。

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

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

　と求められました。

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

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

と求められました。

---