[スポンサーリンク]

一般的な話題

ヒュッケル法(後編)~Excelでフラーレンの電子構造を予測してみた!~

[スポンサーリンク]

bergです。今回は量子化学の黎明期に有機化学の分子軌道論との橋渡しとしての役割を果たしたヒュッケル法(Hückel法)の2回目です。前回はシュレディンガー方程式からヒュッケル法の概要、エチレン、ブタジエン、ベンゼンなどの簡単な分子での計算を追っていき、芳香族性の判定法(ヒュッケル法、フロスト円)までたどり着きましたね。ただし、ヒュッケル法ではn個のπ電子を有する共役系について知見を得るにはn次の行列式(永年方程式)を解く必要があり、複雑な分子では手計算で求めるのが困難でした。そこで、2回目の今回はMicrosoft Excel®のソルバー機能を用いて簡単に解を求める方法をご紹介してまいります。

前回のおさらい

・量子力学の基盤となるシュレディンガー方程式は、厳密には解けない(水素原子除く)。

・ヒュッケル法はπ共役系分子について、分子軌道近似π電子近似LCAO近似を施すことで単純化して解く方法。

・シュレディンガー方程式の左側からΨ*をかけて積分したものをエネルギー固有値Eについて解き、これが最安定となる係数を求める。

・上記はEの各係数での偏微分が0となることを意味し、これをまとめると永年方程式(行列式)が0となる条件を求めることと同値。

・永年方程式の解から各準位のエネルギー固有値と分子軌道のかたちが分かる。

・ポリエン・環状共役系ではエチレンより安定化するケース→共鳴エネルギー(非局在化エネルギー)、芳香族性

・逆にシクロブタジエンでは不安定化→反芳香族性

・芳香族性(4n+2)/反芳香族性(4n)の判定はヒュッケル則に従う。

・ねじれた環では安定性が逆転する(メビウス芳香族性

ざっとこんなところでしょうか。それではいよいよ今回の内容に進んでいきましょう。

ナフタレンの反応性を考察してみた!

まず手始めに、ナフタレン(C10H8)を例にMicrosoft Excel®のソルバー機能を用いて簡単に解を求める方法をご紹介します。

ナフタレン(画像:Wikipedia

スルホン化が代表例ですが、ナフタレンへの芳香族求電子置換(SEAr)反応が一般的に低温ではα位に、高温ではβ位に起こりやすいことは高校化学でもよく題材に取り上げられ、学部生向けの教科書にはそれぞれ速度論支配熱力学支配であることが説明されています。混み合っていないβ置換体が熱力学的に安定なのはともかくとして、それではなぜα位の方が早く反応が進行するのでしょうか?これは従来の有機電子論では説明がつかず、実はフロンティア軌道論の揺籃の地として決定的な役割を果たした難問でもあります[1]。

——————————————————-

以下の計算環境は

ソフトウェア:Microsoft Excel 2019 64bit版

CPU:Intel Core i9-9980HK 2.40GHz

仕様可能メモリ:25.6GB

で実行しています。

——————————————————-

Excelのソルバー機能とは、関数の値が指定値となるような変数の値を探索する機能です。今回の場合、行列式(関数)の値が0となるようなx(変数)を全て求めることになります。

まずは行列を作成します。変数x以外の定数項は手動で入力、xの項はQ5セルから呼び出しています。永年方程式の解(行列式)はMDETERM関数で求めています。

ここまで来たらあとはソルバーを回すだけです。「データ」タブの右側、「分析」のカテゴリのソルバーアイコンを選択し、目的セル(行列式の値)を指定値(0)にするために変数(x)を変化させます。x>3の領域には解がなさそうなので念のため確認してみます。

案の定解なしでした。これでx<3の範囲で探索すればよいことになります。

(解があると予想される領域で見つからない場合にはオプションから「制約条件の制度」を緩めます)

見つかりました。最大の解はx=2.303くらいでした。回答をレポートとして別タブに出力しておくと後々便利です。

このように計算にかかった時間までご丁寧に出力されます。

解が得られるたびに制約条件を変えて次の解を探索していくと、ナフタレンについては、

と定まります。どれも0.1秒足らずで計算できました。10個のエネルギー固有値が重複なく得られたので、軌道の縮退がないことが示されました。

案外重解が一つもなかったので、HOMOのエネルギーがα+0.618β、LUMOがα-0.618βとわかります。

さて、問題はHOMO(x=-0.618)のローブの形です。得られたxの値を永年方程式の導出過程に現れる行列のα-Eに代入し、行列と係数ベクトルの積(MMULT関数で計算可能です)が0となり、2乗和が1となるような10個の係数群を求めます。条件を増やさないと一意には定まらないため対称性に着目して同様にソルバーで解くと、

α位の係数が0.4253、β位が0.2629となります。ゆえにHOMOの電子密度は、

となり、α位がβ位の2.7倍近いという結果が得られ、実験結果の説明がつきました。福井謙一教授らはこのHOMO(フロンティア軌道)の電子密度の差異に着目し、これが求電子試薬との反応性を規定していると考えてフロンティア軌道論を提唱しました。このような計算で分子の反応性を予測できるというのは非常に興味深いですね。

フラーレンの永年方程式を解いてみた!

このたび、π共役系分子の金字塔ともいえるフラーレン(C60過去記事)についても永年方程式が解けるのか試してみました!

フラーレン(画像:Wikipedia

永年方程式は60次になりますが、画像の連番を参考にひたすら作成していきます。

構造式とにらみ合って打ち込ましたが、かなりしんどいです(笑)これでもほとんどの項は0です。なんだか模様に規則性がありそうですね。

ソルバーでしらみつぶしに解を求めていくと

13の解が見つかりました。一部の解は20秒近く計算にかかりました(笑)

答え合わせをすると、フラーレンのエネルギー準位は15個ある[2]ようなのですが、いろいろと条件を変えて試行してもどうしても2つの軌道を発見できませんでした。60次関数では勾配が極めて急峻になるため、見逃しが発生しやすくなるのかもしれません。また、π平面の歪みによって縮退が解ける可能性も考えられそうです。随分苦行した割には寂しい結果となりましたが、特別なソフトを使わずにここまでの結果が得られるとわかったことは大きな収穫でした。学生時代以来ひさびさに量子化学に触れて、良い頭の体操になったと思います。

フラーレンはその高度に対照的な構造のため、高次に縮退した電子構造を持ちます。とりわけLUMOが三重縮退していることから電子受容能に優れ、最大6電子還元まで受けることが知られています。

このような特性を活かした[2+3]双極子付加反応をはじめとする化学修飾が容易であり、有機半導体などの機能性材料をはじめとする応用が嘱望されています。

・・・

ここまでご覧のように、ヒュッケル法は大胆な近似を用いた極めてシンプルな計算手法でありながら、それなりに化合物の物性を説明できる結論を導ける優れた手法です。とはいえ、手計算やExcelでの求解には限界もあります。現在はGAMESSGaussianなどの優れた計算科学ソフトウェアがあり、ハートリー=フォック法(HF法)密度汎関数法などのヒュッケル法よりはるかに高度な計算(非経験的分子軌道法)を誰でも手軽に行うことがでるようになりました。特にGAMESSは無償版も提供されており、非常にありがたい時代になったといえますね。

関連書籍

現代量子化学の基礎

現代量子化学の基礎

中島 威, 藤村 勇一
¥3,740(as of 06/13 18:27)
Amazon product information
「量子化学」のことが一冊でまるごとわかる

「量子化学」のことが一冊でまるごとわかる

齋藤 勝裕
¥1,870(as of 06/13 21:52)
Amazon product information

参考文献

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

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

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

[2] M. S. Golden et al, Journal of Physics: Condensed Matter, 1995, 7,8219-8247.

https://doi.org/10.1088/0953-8984%2F7%2F43%2F004

…フラーレンの電子構造について

gaming voltammetry

berg

投稿者の記事一覧

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

関連記事

  1. 4つの異なる配位結合を持つ不斉金属原子でキラル錯体を組み上げる!…
  2. タングトリンの触媒的不斉全合成
  3. 「フント則を破る」励起一重項と三重項のエネルギーが逆転した発光材…
  4. 有機合成化学協会誌2022年10月号:トリフルオロメチル基・気体…
  5. ホウ素と窒素で何を運ぶ?
  6. 個性あるTOC
  7. sp2-カルボカチオンを用いた炭化水素アリール化
  8. 「ソーシャルメディアを活用したスタートアップの価値向上」 Blo…

注目情報

ピックアップ記事

  1. 革新的医薬品の科学 薬理・薬物動態・代謝・安全性から合成まで
  2. ダイヤモンド構造と芳香族分子を結合させ新たな機能性分子を創製
  3. カリウム Potassium 細胞内に多量に含まれる元素
  4. アラスカのカブトムシは「分子の防寒コート」で身を守る
  5. エイモス・B・スミス III Amos B. Smith III
  6. ヒドロメタル化 Hydrometalation
  7. スタチンのふるさとを訪ねて
  8. 森本 正和 Masakazu Morimoto
  9. エチオ・リザード Ezio Rizzardo
  10. 東大薬小林教授がアメリカ化学会賞を受賞

関連商品

ケムステYoutube

ケムステSlack

月別アーカイブ

2021年2月
1234567
891011121314
15161718192021
22232425262728

注目情報

最新記事

水-有機溶媒の二液相間電子伝達により進行する人工光合成反応

第662回のスポットライトリサーチは、京都大学 大学院工学研究科 物質エネルギー化学専攻 阿部竜研究…

ケムステイブニングミキサー 2025 報告

3月26日から29日の日本化学会第105春季年会に参加されたみなさま、おつかれさまでした!運営に…

【テーマ別ショートウェビナー】今こそ変革の時!マイクロ波が拓く脱炭素時代のプロセス革新

■ウェビナー概要プロセスの脱炭素化及び効率化のキーテクノロジーである”マイクロ波…

予期せぬパラジウム移動を経る環化反応でベンゾヘテロールを作る

1,2-Pd移動を含む予期せぬ連続反応として進行することがわかり、高収率で生成物が得られた。 合…

【27卒】太陽HD研究開発 1day仕事体験

太陽HDでの研究開発職を体感してみませんか?私たちの研究活動についてより近くで体験していただく場…

熱がダメなら光当てれば?Lugdunomycinの全合成

光化学を駆使した、天然物Lugdunomycinの全合成が報告された。紫外光照射による異性化でイソベ…

第59回有機反応若手の会

開催概要有機反応若手の会は、全国の有機化学を研究する大学院生を中心とした若手研究…

多環式分子を一挙に合成!新たなo-キノジメタン生成法の開発

第661回のスポットライトリサーチは、早稲田大学大学院先進理工学研究科(山口潤一郎研究室)博士課程1…

可視光でスイッチON!C(sp3)–Hにヨウ素をシャトル!

不活性なC(sp3)–H結合のヨウ素化反応が報告された。シャトル触媒と光励起Pdの概念を融合させ、ヨ…

化学研究者がAIを味方につける時代―専門性を武器にキャリアを広げる方法―

化学の専門性を活かしながら、これからの時代に求められるスキルを身につけたい——。…

実験器具・用品を試してみたシリーズ

スポットライトリサーチムービー

PAGE TOP