[スポンサーリンク]

一般的な話題

ヒュッケル法(前編)~手計算で分子軌道を求めてみた~

[スポンサーリンク]

bergです。今回からは2回に分けて、量子化学の黎明期に有機化学の分子軌道論との橋渡しとしての役割を果たしたヒュッケル法(Hückel法)について詳しく見ていきます。

ヒュッケル法とは量子化学の初期に提案されたもっとも単純な分子軌道法のひとつで、発見者であるドイツの物理化学者、エーリヒ・ヒュッケル(Erich Armand Arthur Joseph Hückel)教授にちなんで名づけられたエポニムです。多彩な研究を精力的に行ったことでも知られ、ヒュッケル法、関連するヒュッケル則(Hückel則)にとどまらず、電解質溶液中のイオンの相互作用について統計力学を駆使して記述したデバイ-ヒュッケルの式(Debye–Hückel equation)にもその名を残しています。

1938年のヒュッケル教授(画像:Wikipedia

ヒュッケル法(、および、それを応用してヘテロ原子を含む分子の議論を可能とした拡張ヒュッケル法)では必ずしも実態に即した結果が得られるわけではありませんが、分子軌道計算用のプログラムを使用せずに簡単な手計算のみで示唆に富んだ結果が得られることから今日でも大学の量子化学・物理化学でなじみ深い題材とされています。後年ウッドワード(Robert Barns Woodward)教授や福井謙一教授によって開発され、「化学反応に関する我々人類の理解を深めてゆくうえでの一里塚」としてノーベル化学賞に輝くこととなったフロンティア軌道論にも大きく寄与しています[1]。

量子化学の黎明

現代物理学の根幹として量子化学の基礎ともなっている量子力学は、オーストリアの物理学者、エルヴィン・シュレディンガー教授が1926年に導出したシュレディンガー方程式(Schrödinger equation)に端を発しています。

1933年のシュレディンガー教授(画像:Wikipedia

シュレディンガー方程式では系の状態(量子状態)が時間とともに逐次変化していく(シュレディンガー描像)ものとみなし、その解は波動関数を与えます。ただし量子化学においては、その対象となる原子や分子について電子を束縛するエネルギー場とみなし、電子そのものの運動ではなく空間ごとの存在確率(電子密度)に注目した方が便利な場合があります。そこで、シュレディンガー方程式を時間と空間について変数分離し、以下の時間に依存しないシュレディンガー方程式として取り扱うことが一般的です。

ここではハミルトニアン演算子、Ψは空間座標を変数とする波動関数、Eはエネルギー固有値を表します。すなわちこれは固有値方程式であり、波動関数は原子/分子軌道の定常状態、エネルギー固有値はその軌道準位のエネルギー値を表します。それぞれの固有状態は直交します。

量子力学の分野において、以上のシュレディンガー方程式はまさに原点にして頂点とも呼ぶべき極めて重要な式です。しかし、そんな完璧なシュレディンガー方程式にも、量子化学に応用する上で致命的な問題点を抱えていました。

…それは、解析的に解くことができない点です。シュレディンガー方程式は電子が1つしかない水素原子については厳密解を求めることができ、それゆえにs軌道からf軌道まで様々な軌道の形状が調べられています。しかしながら、原子核と電子の数が3つ(三体問題)以上の多体問題においては正確に解くことは不可能です。

解法がない以上、様々な仮定や近似を重ねてなるべく実態の分子に近い計算結果を出すことが求められます。量子化学はこうして様々な近似を導入することで発展してきました。

ヒュッケル法と実例

(注:以下、わかりやすくするため導出を省いている部分があります)

ヒュッケルは二重結合と単結合が交互に繰り返し現れるような構造を持ったπ共役系の有機分子について、次のような仮定をして近似計算を試みました。

①分子軌道を構成原子の原子軌道の重ね合わせと仮定(分子軌道近似

②(エネルギー値の近い)p軌道以外の寄与を無視(π電子近似

③分子軌道を原子軌道の線形結合として記述(LCAO近似 ; Linear Combination of Atomic Orbital)

さて、さきほどの時間依存しないシュレディンガー方程式の両辺に左側から波動関数Ψの複素共役関数であるΨ*をかけると、

これを全範囲で積分すると、

よって、エネルギー固有値Eは以下のように書き下されます。

以下ではこの式を前提として議論していきます。

・エチレン(CH2=CH2

π結合をもつもっとも単純な有機分子として、まずエチレンを例に考えていきましょう(共役系とはいわない気がしますが)。

エチレン(画像:Wikipediaより改変)

エチレンは群論でいうところのD2h対称な分子であるため炭素原子に区別はないのですが、便宜的に1番、2番と連番しました。それぞれの炭素原子の原子軌道(p軌道)の波動関数をφ1、φ2とすると、分子軌道はその線形結合として近似するため

となります。これをさきほどのエネルギー固有値の式に代入して、

これをそのまま展開すると煩雑となるため、以下の文字を導入します。

・クーロン積分:α

・共鳴積分:β

・重なり積分

共鳴積分については、ほかの化合物の事例でも結合の有無によって値が規定されていると考えていただいて差し支えありません。

これらを用いて簡潔に展開すると、

軌道のエネルギーを最小とするc1、c2を求めるためにE=N/Dとおいて偏微分すると、

であるから、

同様に、

この非自明解は、

の行列式(永年方程式)を解くことで与えられます。

よって、エネルギー固有値の解は以下の2つとなります。

規格化条件の制約から、波動関数は

と表すことができます。

この結果の解釈としては、以下の図のように考えることができます。

エチレンの分子軌道(画像:Wikipediaより改変)

同位相の原子軌道を重ねたHOMOと逆位相のLUMOがあり、そのエネルギーギャップが2βです。

エチレンをはじめアルケンは求電子試薬の攻撃を受けて付加反応を起こしやすいことで知られていますが、求電子試薬はエチレンのHOMOの電子密度が高くなる方向、すなわちπ平面に対して鉛直方向から接近し、反応します。なお、各原子上の電子密度は原子軌道にかかる係数の2乗に比例します。LUMOでは炭素原子の中間に節面があり、位相が逆転します。

また、π結合次数は隣接する原子の原子軌道にかかる係数の積に電子数をかけたものとして表現できるので、

となります。実際にはσ結合もあるため結合次数は2です。

全電子エネルギーEtotalも単純に求められ、2(α+β)となります。

なお、途中の解き方として、

と定義すると、任意のn次の永年方程式はxのn次方程式となるので簡単になります。この性質を用いて、以下ではもっと複雑な有機分子に着目していきましょう。

・アリル共役系

3炭素からなるアリル共役系(アリルアニオン、アリルラジカル、アリルカチオン)の場合、永年方程式は3次になります。

サラスの公式を用いて解くと

となります。さきほどと同様に規格化条件を用いると

として波動関数が求まります。

アリル共役系の分子軌道(筆者作成)

・シクロプロペニル共役系

一方、環の閉じたシクロプロパン骨格を有するシクロプロペニル共役系(シクロプロペニルアニオン、シクロプロペニルラジカル、シクロプロペニルカチオン)ではどうなるでしょうか。

永年方程式は以下のようになります。

これを解くと、

が得られます。

ここで重解が得られたことから、シクロプロペニル共役系の場合には縮退した2つの軌道があります。規格化条件のみではこれらの解が一意に定まらず、Ψ2とΨ3直交することを利用してΨ2においてc3=0とする必要があります。これより、

が求められます。

シクロプロペニル共役系の分子軌道(筆者作成)

まったく異なる形状の軌道が同じエネルギーというのもなんとも不思議ですよね。

この全電子エネルギーをアリルと比較すると、アニオンでは不安定化している一方でカチオン・ラジカルでは安定化されていることがわかります。

特にシクロプロペン誘導体から調製されるシクロプロペニルカチオンは極めて安定で、最小の芳香族化学種としても知られています。(→過去記事

このように三員環構造の方が安定になるケースは水素原子のみで構成されたH3+(プロトン化水素分子)でもみられ、三中心二電子結合をもつ星間物質として比較的豊富に存在します。

・ブタジエン

さて、最小のポリエンともいえるブタジエンではどうなるでしょうか。

ブタジエン(画像:Wikipedia

永年方程式は

ですので、これを地道に次数下げして解くと

となります。これより波動関数は

と求まります。

ブタジエンの分子軌道(画像:Wikipediaより改変)

ブタジエンではπ共役のため、結合次数が構造式のような単結合や二重結合とはなりません。エチレンと同様に計算すると、π結合次数は

となり、π電子が非局在化していることがうかがえます。一方でπ電子密度は

ですべて同一です。

ブタジエンの結果をエチレンと比較すると、いくつか興味深い考察が得られます。

まず、ブタジエンがエチレン2分子よりも安定であることです。全電子エネルギーの差し引きにより、ブタジエンの方が

だけ安定化されていることがわかります。これが共鳴エネルギー(非局在化エネルギー)です。

また、HOMO-LUMOギャップエチレンの2βからブタジエンでは1.236βに減少していることも見逃せません。この傾向はさらにポリエンの鎖長(共役系)が長くなっても続き、UV-Visの吸収波長は長波長側へレッドシフトします(導出は割愛します)。

・シクロブタジエン

さて、ブタジエンの両端を環化させたシクロブタジエンは最小の反芳香族分子として知られています。

永年方程式は

となり、この解は

となります。重解に対しては直交条件を用いることで波動関数は

と求められます。

シクロブタジエンの分子軌道(筆者作成)

E=αの軌道が縮退しているため、基底状態が三重項となります。

ここで全電子エネルギーに注目すると、エチレン二分子の値よりも

も不安定化する結果となっています。そのため、シクロブタジエンは安定に存在できない反芳香族性を有する分子であり、依然としていくつかの誘導体が単離されるのみにとどまっています。

さて、シクロプロペニルカチオンやシクロブタジエンなどの環状共役系分子(平面・単環の場合)においては、ヒュッケル法を用いて求めたn員環のエネルギー固有値の解がE=αに中心を持つ半径2βの円(フロスト円)に内接する正n角形を用いて表されることが知られています。

シクロプロペニルカチオンのFrost円(画像:Wikipedia

Frost円を用いると、環状共役系が安定化される(芳香族性をもつ)のは、nを整数としてそのπ電子数が(4n+2)を満たすときであることが示されています。これをヒュッケル則と呼びます。逆に、π電子数が4n個のときは不安定化されて反芳香族性を持ちます。

非常に簡明なヒュッケル則ですが、欠点が皆無というわけではありません。たとえばシクロオクタテトラエン(π電子数8)は反芳香族性を持ちそうですが、不安定化を回避するために非平面の構造をとり、共役しない構造をとっています(非芳香族)。このように配座などの要因も絡むためにヒュッケル則は中員環・大員環のアヌレンにおいては必ずしも当てはまりません。

また、物理的にねじれた大員環ではp軌道の位相が逆転するため、行列式の最も右上と左下の項が-1となります。これを解くとヒュッケル則とは逆に、π電子数が4n個の時に安定化します。これをメビウス芳香族性と呼びます。(→過去記事

・ベンゼン(C6H6

ベンゼン(画像:Wikipedia

芳香族分子の代表格、ベンゼン。π電子数は6で、先のヒュッケル則を満たします。永年方程式は6次となりますが、対称性が高いためきれいに展開できます。

これを解いて

また、

ベンゼンの分子軌道(画像:Wikipediaより改変)

さすがにこれ以上複雑な分子について永年方程式を手計算で解いていくのはかなり無理がありそうですね。そこで、次回は特殊な計算科学ソフトウェアを使わずにMicrosoft Excel®のソルバー機能を用いて簡単に解を求める方法をご紹介します。得られた結果からナフタレン(C10H8)の反応性(なぜ低温ではα位上での置換が優先する?)や、フラーレン(C60)の電子構造にもせまっていきますのでお楽しみに!

関連書籍

参考文献

[1] 稲垣 都士, 池田 博隆, 化学と教育, 2019, 67(1), p.28-31,

https://doi.org/10.20665/kakyoshi.67.1_28

…フロンティア軌道論についての総説

berg

berg

投稿者の記事一覧

化学メーカー勤務。学生時代は有機をかじってました⌬
電気化学、表面処理、エレクトロニクスなど、勉強しながら執筆していく予定です

関連記事

  1. 交互に配列制御された高分子合成法の開発と機能開拓
  2. ルィセンコ騒動のはなし(後編)
  3. V字型分子が実現した固体状態の優れた光物性
  4. 複数のイオン電流を示す人工イオンチャネルの開発
  5. プラスマイナスエーテル!?
  6. 2018年1月20日:ケムステ主催「化学業界 企業研究セミナー」…
  7. 活性が大幅に向上したアンモニア合成触媒について
  8. 含ケイ素四員環-その2-

コメント、感想はこちらへ

注目情報

ピックアップ記事

  1. 実験を加速する最新機器たち|第9回「有機合成実験テクニック」(リケラボコラボレーション)
  2. 芳香族化合物のC–Hシリル化反応:第三の手法
  3. 保護により不斉を創る
  4. 有機合成テクニック集[ケムステ版]
  5. ナノ粒子の機能と応用 ?コロイダルシリカを中心に?【終了】
  6. 積水化学、工業用接着剤で米最大手と提携
  7. 大阪大学インタラクティブ合宿セミナーに参加しました
  8. 新形式の芳香族化合物を目指して~反芳香族シクロファンにおける三次元芳香族性の発現~
  9. デスソース
  10. 液体キセノン検出器

関連商品

ケムステYoutube

ケムステSlack

注目情報

注目情報

最新記事

広瀬すずさんがTikTok動画に初挑戦!「#AGCチャレンジ」を開始

TikTok For BusinessとAGC株式会社は、AGCをより多くの人に知っていただくことを…

新規性喪失の例外規定とは?

bergです。今回は論文投稿・学会発表と特許出願を同時に行うための新規性喪失の例外規定の適用手続きに…

新車の香りは「発がん性物質」の香り、1日20分嗅ぐだけで発がんリスクが高まる可能性

「新車の香り」には、がんや生殖障害、子どもの先天性欠損症などを引き起こす可能性があるベンゼンやホルム…

溶液を流すだけで誰でも簡単に高分子を合成できるリサイクル可能な不均一系ラジカル発生剤の開発

第 295 回のスポットライトリサーチは東京大学豊田研究室の博士課程 1 年 岡美奈実さんと修士課程…

Carl Boschの人生 その9

Tshozoです。書いてると色々膨らんで収集がつかなくなりますね。ということで前回の続き。W…

創薬・医療系ベンチャー支援プログラム”BlockbusterTOKYO” ビジネスプラン発表会を開催!

東京都が主催し、Beyond Next Ventures株式会社が運営するBlockbuster T…

酸化反応を駆使した(-)-deoxoapodineの世界最短合成

第294回のスポットライトリサーチは、吉田慶 博士にお願いしました。今回取り上げる研究は有機…

特許取得のための手続き

bergです。本記事では特許出願に必要な手続きについてかいつまんでご紹介します。皆さんの研究もひょっ…

Chem-Station Twitter

PAGE TOP