[スポンサーリンク]

一般的な話題

計算化学:汎関数って何?

先日の学会ポスター発表会場にて

 

発表者:「・・・というわけで、DFT計算によりこの反応経路を証明しました。」

質問者:「素晴らしい発表でした。しかし、それDFT計算じゃないですよね?」

発表者:「えっ!DFT計算ですけど。。。」

質問者:「DFTってB3LYPとかですよね?M06ってDFTじゃないですよね?」

発表者:「M06はDFT計算で使う汎関数のひとつですが、、、」

質問者:「、、、失礼しました。」

 

 

学会等で計算化学についての発表を見ていると「全くDFT計算分かっていないな」と感じる事が多々あります。そんな勘違いをなくすために、今回は汎関数について非常に簡単に説明したいと思います。

 

 

量子計算化学

量子化学計算では、ある分子の波動関数Ψを求める事が出来れば、その分子に関する様々な情報が得られるという考えに基づいています。それは例えば、分子の最安定構造、反応性、物性の予測などなどです。しかし、Schrödinger方程式を厳密に解く事は、実際にはほぼ不可能です。そこで、計算化学では厳密ではないにせよ、Schrödinger方程式の解を求める方法を用います。

ひとつ目は、ハミルトニアンを計算しやすい形に変える事です。これは、計算化学においては、どのような計算手法を選択するかということになります。ハミルトニアンの数式を解けるものに変えてしまうものの例として、HF法DFT法等があり、解ける部分を先に解き、解けない部分の効果を後から近似的に足す例として、MP2法などがあります。

ふたつ目は、波動関数の自由度を制限する事です。これは、計算化学においては、どのような基底関数を選択するかという事になります。これについては、今回の記事では触れず次回説明します。

量子化学計算では、inputファイルで「計算手法/基底関数」といった形で指定します。

 

DFT計算って何?

HF法では、ハミルトニアンの中の一部を近似し波動関数を直接求めようとするのに対し、DFT計算では、波動関数の代わりに電子密度を求めようとします。

DFT計算では、波動関数に対するSchrödinger方程式と等価な電子密度に対する方程式としてKohn-Sham方程式を用います。

Schrödinger 方程式

Kohn-Sham方程式は交換相関汎関数Veffを用います。この汎関数に何を用いるかによって計算の精度が決まります。
Kohn-Sham

汎関数

交換相関汎関数は、交換汎関数相関汎関数からなり、それぞれ電子のスピンに由来する効果電子相関に由来する効果を記述しています。

汎関数の種類

それぞれいろいろな関数が知られていますが、相性が良いものを組み合わせなければなりません。使いやすいように、これらがセットとなっている交換相関汎関数も下図に示すように知られています。

交換相関汎関数

 

これらの中で最も有名なものがB3LYPです。つまり、冒頭の会話で質問者は、「B3LYP=DFT計算」と考えていましたが、DFT計算で用いられる様々な汎関数のうちのひとつにB3LYPがあるという方が正しいのです。

B3LYPは最も良く用いられている汎関数ですが、B3LYPにも欠点はあります。水素結合、van der Waals力を見積もれなかったり、遷移金属の電子状態を正確に計算できなかったりという点です。

万能な汎関数というものは存在せず、系によって最適な汎関数を用いなければなりません。欠点を修正するために長距離補正や分散力補正を組み込んだものや、計算時間を短縮するために半経験的汎関数などもあります。この辺については、各自でお調べください。

次回は、基底関数のことについて簡単に触れたいと思います。

 

関連書籍

The following two tabs change content below.
ゼロ

ゼロ

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

関連記事

  1. ネオ元素周期表
  2. 多彩な蛍光を発する単一分子有機化合物をつくる
  3. 文具に凝るといふことを化学者もしてみむとてするなり⑪:どっちもク…
  4. 2004年ノーベル化学賞『ユビキチン―プロテアソーム系の発見』
  5. ReaxysPrize2015ファイナリスト発表!
  6. 【追悼企画】生命現象の鍵を追い求めてー坂神洋次教授
  7. 実験ノートの字について
  8. 光分解性シアニン色素をADCのリンカーに組み込む

コメント

  1. うーむ

  2. pure DFT以外は認めない原理主義者の質問かと思ってたら、hybridのB3LYPとか例に上げやがってみたいな話かと思ったら全然違ってた。

  3. 「DFT計算では、波動関数に対するSchrödinger方程式と等価な電子密度に対する方程式としてKohn-Sham方程式を用います。Kohn-Sham方程式は交換相関汎関数Veffを用います。」

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

注目情報

ピックアップ記事

  1. 手術中にガン組織を見分ける標識試薬
  2. 塩化インジウム(III):Indium(III) Chloride
  3. 構造式の効果
  4. コルベ電解 Kolbe Electrolysis
  5. Essential Reagents for Organic Synthesis
  6. 抗がん剤大量生産に期待 山大農学部豊増助教授 有機化合物生成の遺伝子発見
  7. 超難関天然物 Palau’amine・ついに陥落
  8. 小林 修 Shu Kobayashi
  9. C70の中に水分子を閉じ込める
  10. 芳香族性に関する新概念と近赤外吸収制御への応用

注目記事

関連商品

注目情報

試薬検索:東京化成工業



注目情報

最新記事

芳香族ニトロ化合物のクロスカップリング反応

第116回のスポットライトリサーチは、京都大学大学院工学研究科 材料化学専攻 有機材料化学講座(中尾…

バイエルスドルフという会社 ~NIVEA、8×4の生みの親~

Tshozoです。女装とかそういう趣味は無いのですが嫁さん(実在)に付き合って化粧品コーナを回ること…

化学系プレプリントサーバー「ChemRxiv」のβ版が運用開始

2017年8月14日、米国化学会(ACS)は、化学分野のプレプリントサーバー“ChemRxiv”のベ…

光触媒で人工光合成!二酸化炭素を効率的に資源化できる新触媒の開発

第115回のスポットライトリサーチは、東京工業大学 理学院 化学系 博士後期課程2年の栗木 亮さんに…

誰も教えてくれなかった 実験ノートの書き方 (研究を成功させるための秘訣)

概要悪い例とよい例を比較しながら,実験ノートを具体的にどう書けばよいのかを懇切丁寧に説明する…

神経変性疾患関連凝集タンパク質分解誘導剤の開発

第114回のスポットライトリサーチは、東京大学大学院薬学系研究科博士後期課程2年の山下 博子(やまし…

Chem-Station Twitter

PAGE TOP