ヒュッケル法のソースを表示

{{出典の明記|date=2018年8月}}
'''ヒュッケル法'''（ヒュッケルほう、{{lang-en-short|H&uuml;ckel method}}）は、[[ドイツ]]の[[エーリヒ・ヒュッケル]]によって提案された[[分子軌道法]]である。[[エチレン]]や1,3-[[ブタジエン]]、[[ベンゼン]]等の[[π電子]][[共役系|共役]]分子に適用する例が入門用としてよく用いられる。

ヒュッケル法では電子に関する[[積分]]に対して以下のような近似を導入する。
*[[重なり積分]]の値は同じ原子軌道同士では1、異なる原子軌道の間では0
*クーロン積分（ハミルトニアン行列の対角要素）の値は同じ種類の原子では等しいものとし、パラメータ<math>\alpha</math>によって表す
*共鳴積分（ハミルトニアン行列の非対角要素）の間に結合を持つ原子間でのみ0でない値をもつとし、パラメータ<math>\beta</math>によって表す。

しばしば紹介されるπ共役系に対する計算では、その大胆な仮定にもかかわらず定性的に正しい結果を与え成功をおさめた。しかしながら結合を持つ原子間と持たない原子間とで積分の扱いを変えるため、構造を特定できていない分子に対しては適用できない。その欠点を改善するために、導入する近似を少なくし計算する積分の量を少し増やした[[拡張ヒュッケル法]]と呼ばれる方法があり、そちらと対比して単純ヒュッケル法とも呼ぶ。現在では計算機の性能の向上などにより単純ヒュッケル法を用いる必要性はほとんどなく、ただヒュッケル法、と言った場合には[[拡張ヒュッケル法]]を指すことが多い。

==ヒュッケル法の仮定==
:<math>S_{ij}
= \int \phi_i^* \phi_j d\tau
= \begin{cases} 1 & i= j \\ 0 & i \ne j \end{cases}</math>

[[重なり積分]]の値 <math>S_{ij}</math> は同じ原子軌道同士では1であり、異なる原子軌道間では0である。

:<math>H_{ij}
= \int \phi_i^* \hat{H} \phi_j d\tau
=  \begin{cases} \alpha & i= j \\ \beta & i \ne j \end{cases}</math>

<math>\alpha</math>をクーロン積分と呼び、<math>\beta</math>を共鳴積分と呼ぶ。結合を持たない原子間では<math>\beta</math>の値は0。

==1,3-ブタジエンの例==
[[ファイル:2p-1-3-butadiene.png|right|]]1,3-ブタジエンの[[分子軌道]] <math>\psi</math> を4つのπ軌道の線形結合で表す。<math>c_1 \dots c_4</math>は各π軌道の分子軌道への寄与の大きさである。
{{Indent|<math>\psi = c_1\phi_1 + c_2\phi_2 + c_3\phi_3 + c_4\phi_4 </math>}}

変分法によってエネルギー<math>E</math>を求めると、

{{Indent|<math>E = \frac{\int \psi^* \hat H \psi d \tau}{\int \psi^* \psi d \tau}</math>}}

この右辺を極小にする<math>c_1 \dots c_4</math>を定めることになる。実際に<math>\psi</math>を代入すると、

{{Indent|<math>E = \frac{\sum c_i c_j H_{ij}}{\sum c_i c_j S_{ij}}</math>}}

と表せる。これを極小とするような<math>c_1 \dots c_4</math>を得る条件は、

{{Indent|<math>\frac{\partial E}{\partial c_k}=0</math>　　　<math>\left( k=1,2,3,4 \right)</math>}}

である。これらの式が<math>c_1=c_2=\dots=0</math>以外の解を持つには、

:<math>
  \begin{vmatrix}
    H_{11} - ES_{11} & H_{12} - ES_{12} & \dots & H_{14} - ES_{14} \\
    H_{21} - ES_{21} & H_{22} - ES_{22} & \dots & H_{24} - ES_{24} \\
    \dots & \dots & \ddots & \vdots & \\
    H_{41} - ES_{41} & H_{42} - ES_{42} & \dots & H_{44} - ES_{44} \\
  \end{vmatrix}= 0
</math>

である必要がある。これを永年方程式と呼ぶ。
ヒュッケル法の仮定より、

{{Indent|<math>
    \begin{vmatrix}
         \alpha-E & \beta & 0 & 0 \\
         \beta & \alpha-E & \beta & 0 \\
         0 & \beta & \alpha-E & \beta & \\
         0 & 0 & \beta & \alpha-E \\
    \end{vmatrix}= 0
</math>}}

であるから、左辺を展開して、

{{Indent|<math>
\left( \alpha-E \right)
    \begin{vmatrix}
         \alpha-E & \beta & 0 \\
         \beta & \alpha-E & \beta & \\
         \ 0 & \beta & \alpha-E \\
    \end{vmatrix}
-\beta
    \begin{vmatrix}
         \beta & \beta & 0 \\
         0 & \alpha-E & \beta & \\
         0 & \beta & \alpha-E \\
    \end{vmatrix}
</math><br />
<math>
    =\left( \alpha-E \right)^4 -3\left( \alpha-E \right)^2 \beta^2 + \beta^4 
    =0
</math>}}

<math>x = \left( \alpha-E \right)^2 / \beta^2</math>　と置くと、行列式は

{{Indent|<math>
    x^2 - 3x + 1 = 0
</math>}}
となり、その根は<math>x = 2.62 , 0.38</math>である。したがって、ブタジエンの4個の分子軌道のエネルギーは、

{{Indent|<math>E = \alpha \pm 1.62\beta, \alpha \pm 0.62\beta</math>}}

である。

{{DEFAULTSORT:ひゆつけるほう}}
[[Category:半経験的分子軌道法]]
[[Category:量子力学]]
[[Category:量子化学]]
[[Category:計算化学]]
[[Category:分子物理学]]
[[Category:電子軌道]]
[[Category:エーリヒ・ヒュッケル]]
[[Category:化学のエポニム]]