[スポンサーリンク]

一般的な話題

計算化学:DFT計算って何?Part II

 

近年、DFT計算により実験結果をサポートする論文が増えて来ています。実験屋さんでも計算化学についての知識を持つことは必須となりつつあります。

しかし、DFTって何なのでしょうか?DFT計算が一般的になりつつある昨今、「DFT計算って何?」という基本的な質問をしにくい雰囲気があると思います。

周りの人に思い切って聞くことができず、夜中にこっそりWikipediaで調べてみても「密度汎関数法」と書いてあるだけで「???」な感じです。そう、計算化学を分かりやすく説明しているウェブサイト、教科書はあまりないのです。計算化学が一部の人にしか普及しないのは、難しいとかではなく、バイブル的な教科書が存在しないためだと個人的には思っています。

 

前回の「汎関数って何?」という記事のpart IIに当たる今回の記事では、量子化学計算のかんたんな歴史から説明し、次の記事(PartIII)と併せて、初心者の人でもDFT計算が何かを理解できるようにしようと思います。

今回の記事でも細かいところは省略していくので、詳細を知りたい勉強熱心な読者の方は一番下の参考図書を読んで下さいね!

 

量子力学の発展

量子化学歴史

前回の記事に書いたように、「物質の性質、反応性は、その物質の電子状態が分かれば予測できる」という考えに基づき、有機化学者は量子化学計算を行なっています。電子状態を知る=シュレーディンガー方程式を解くということになっています。

 

このシュレーディンガー方程式ですが、1926年に発表されています。さらにシュレーディンガー方程式に相対論を盛り込んだディラック方程式は1928年に発表されています。

この時期は量子力学が一気に進展したときで、ハートリー法(1928年Hartree)、スレーター行列式(1929年Slater)、ハートリーフォック法(1930年Fock, Slater)、分子軌道法(1926年Hund、1927年Mulliken)、メラープリセット摂動法(1934年Moller, Plesset)、TDHF法(1930年Dirac)LDA交換汎関数(1930年Dirac)など、まだまだ書ききれないくらい多くの理論が生み出されました。この時期のキーワードは、「実際の分子中のシュレーディンガー方程式をどのように解いたら良いのだろうか?」でした。

 

余談ですが、量子力学の歴史を勉強していると、さまざまな理論の提唱、開発で大きな貢献をしているDirac、、、天才過ぎだろって思ってしまいます。Diracは、ノーベル賞の受賞を「有名になるのが嫌だ」という理由で辞退しようとして、周りの人達に全力で止められたという少し変わった人なのですが、、、すごい人ですね。

 

さて、量子力学の歴史の説明に戻ります。この後、「分子の中の電子運動の波動関数をどのように解釈すればよいか?」ということをキーワードにし、混成軌道モデル(1928年Pauling)、遷移状態理論(1935年Eyring)、LCMO近似(1929年Lennard-Jones、1938年Coulson)、化学反応原理(1936年Bell、1938年Evans, Polanyi)などが発表されました。

 

ここまで、いろいろな理論を羅列してきましたが、要約すると

「シュレーディンガー方程式を厳密に解くのは不可能→近似的に解く理論の開発が1930年前後に盛んに行なわれた」

です。

 

DFT (密度汎関数法)の歴史

さて、ここまで紹介した理論は「シュレーディンガー方程式を解くことは難しいので近似的に解いてしまおう!」という感じでした。それに対し、DFTは「裏技的にハミルトニアンを求めてしまおう!」みたいなイメージです。

DFT

DFTの基礎理論となるトーマス・フェルミ理論は1927年に発表されました。簡単に言うと電子密度だけでハミルトニアン演算子を表わすことが出来ますよ!という理論です。

 

しかし、この理論には解の一意性や汎関数の存在を保証する物理的な裏付けは何も無く、また化学結合すら全く再現できないため、1960年まで忘れ去られていました。現在では皆が使っているDFTも、当時は「使えないな、コレ」みたいな感じだったんですね。

 

1964年、トーマス・フェルミ理論のコンセプトの正しさを物理的に裏付ける定理が提案されました。ホーエンベルグ・コーン定理です。これは下記の2つの定理からなっています。

1.外場ポテンシャルは電子密度で決定される。

2.あらゆる電子密度について、常にエネルギーの変分原理が成り立つ。

 

しかし、第一定理の証明で電子密度が波動関数と一対一対応すると仮定していることが若干問題であり、これはV表現可能性問題と呼ばれています。1979年レヴィは制限つき探索法を提案してこの問題を解決しました。また、第二定理の証明においてもN表現可能性問題という問題がありました。

 

コーン・シャム方程式

前回の記事ではDFT計算=コーン・シャム方程式のような書き方をしたため、「トーマス・フェルミ理論?」「ホーエンベルグ・コーン定理?」と頭の中が混乱しているかもしれませんが、ここからいよいよコーン・シャム方程式が登場します。

Kohn-Sham2

ホーエンベルグ・コーン定理によりトーマス・フェルミ理論の正しさが立証されましたが、実際の電子状態の計算はまだ出来ませんでした。この定理に基づく計算法は、翌年1965年に発表されたコーン・シャム方程式を用いることによりやっと可能になりました。

コーン・シャム方程式では、運動エネルギーの計算に電子密度の汎関数ではなく、ハートリーフォック法と同様の独立近似の定式を利用しています。このことにより化学・固体物性の定量的な計算が可能となり、DFTの急速な拡大につながりました。

しかし、このことにより、トーマス・フェルミ理論で当初提唱された純粋なDFTと現在使われているDFT計算はイコールではなくなりました。つまり、多くの人が知らないでいるのですが、「DFT理論=コーン・シャム方程式」は厳密には違います。

ここまでを簡単にまとめると、

「1927年にDFT理論が提唱される→実用性なし→1964年にホーエンベルグ・コーン定理により証明される→1965年にコーン・シャム方程式が発表され実際の系に使えるようになった」

です。

ざっと1930年から1980年まで説明しましたが、長くなりましたので、続きは次回の記事で書きたいと思います。

 

参考図書

 

The following two tabs change content below.
ゼロ

ゼロ

女の子。研究所勤務。趣味は読書とハイキング ♪ ハンドルネームは村上龍の「愛と幻想のファシズム」の登場人物にちなんでま〜す。5 分後の世界、ヒュウガ・ウイルスも好き!

関連記事

  1. ちょっとした悩み
  2. Reaxys Prize 2011発表!
  3. イオンペアによるラジカルアニオン種の認識と立体制御法
  4. 向かう所敵なし?オレフィンメタセシス
  5. オキソニウムカチオンを飼いならす
  6. 複雑化合物合成にも適用可能なC-H酸化反応
  7. ヒト胚研究、ついに未知領域へ
  8. 混合原子価による芳香族性

コメント

  1. この記事へのコメントはありません。

  1. この記事へのトラックバックはありません。

注目情報

ピックアップ記事

  1. ついったー化学部
  2. “マイクロプラスチック”が海をただよう その1
  3. ナノってなんて素敵ナノ
  4. Keith Fagnou Organic Chemistry Symposium
  5. トムソン:2007年ノーベル賞の有力候補者を発表
  6. イヴ・ショーヴァン Yves Chauvin
  7. 細野 秀雄 Hideo Hosono
  8. 中国へ講演旅行へいってきました②
  9. 学振申請書を磨き上げるポイント ~自己評価欄 編(前編)~
  10. 過ぎ去りし器具への鎮魂歌

注目記事

関連商品

注目情報

試薬検索:東京化成工業



最新記事

ニホニウムグッズをAmazonでゲットだぜ!

まだまだ続くよニホニウムフィーバー!ということで、色々な書籍などは巷にあふれてますね。先日当…

高機能・高性能シリコーン材料創製の鍵となるシロキサン結合のワンポット形成

第92回のスポットライトリサーチは、産業技術総合研究所触媒化学融合研究センターケイ素化学チーム(島田…

極性表面積 polar surface area

極性表面積(polar surface area, PSA)とは、分子表面のうち極性を帯びている部分…

Late-Stage C(sp3)-H活性化法でステープルペプチドを作る

バルセロナ大学・Fernando Albericioらは、パラジウム触媒によるLate-Stage …

企業研究者たちの感動の瞬間: モノづくりに賭ける夢と情熱

概要有機合成化学およびプロセス化学に関係した研究の現場での,企業化や商品化に賭けるそれぞれの…

銅触媒と可視光が促進させる不斉四置換炭素構築型C-Nカップリング反応

カリフォルニア工科大学・Jonas C. PetersおよびGregory C. Fuの共同研究グル…

Chem-Station Twitter

PAGE TOP