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

  1. 行列要素を記述し、関数 matrix で行列に変換します。
  2. 関数 eigenvectors で固有値と固有ベクトルを求め、関数 float で実数表示します。

※ 固有値が大きさの順に出力されない点と、固有ベクトルが正規化されていないことに注意してください。

Ethylene

エチレン

E_1 =  α + β      \[CurlyPhi]_1 = c_1 (φ_1 + φ_2)
E_2 =  α - β      \[CurlyPhi]_1 = c_2 (φ_1 - φ_2)

と求められました。

Allyl

アリル

 

E_1 =  α + 1.414 β      \[CurlyPhi]_1 = c_1 (φ_1 + 1.414 φ_2 + φ_3)
E_2 =  α                     \[CurlyPhi]_2 = c_2 (φ_1 - φ_2)
E_3 =  α - 1.414β       \[CurlyPhi]_3 = c_1 (φ_1 - 1.414 φ_2 + φ_3)

と求められました。

Butadiene

ブタジエン

 固有ベクトルを正規化すると、

  E_1 =  α + 1.618β    
  E_2 =  α + 0.618β      
  E_3 =  α  - 0.618β    
   =  α - 1.618β    

 と求められました。


メニューに戻る