[スポンサーリンク]

一般的な話題

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

[スポンサーリンク]

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

 

発表者:「・・・というわけで、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力を見積もれなかったり、遷移金属の電子状態を正確に計算できなかったりという点です。

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

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

 

関連書籍

ゼロ

ゼロ

投稿者の記事一覧

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

関連記事

  1. Carl Boschの人生 その2
  2. 一重項励起子開裂を利用した世界初の有機EL素子
  3. 高分子を”見る” その2
  4. 誤解してない? 電子の軌道は”軌道”では…
  5. Carl Boschの人生 その7
  6. Ru触媒で異なるアルキン同士をantiで付加させる
  7. NMRの測定がうまくいかないとき(2)
  8. とにかく見やすい!論文チェックアプリの新定番『Researche…

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

注目情報

ピックアップ記事

  1. 指向性進化法 Directed Evolution
  2. グリチルリチン酸 (glycyrrhizic acid)
  3. それは夢から始まったーベンゼンの構造提唱から150年
  4. 口頭発表での緊張しない6つのヒント
  5. トヨタ、世界初「省ネオジム耐熱磁石」開発
  6. 2007年度ノーベル化学賞を予想!(3)
  7. natureasia.com & Natureダイジェスト オンラインセミナー開催
  8. クリック反応の反応機構が覆される
  9. NHC (Bode触媒2) と酸共触媒を用いるα,β-不飽和アルデヒドとカルコンの[4+2]環化反応
  10. コランニュレンの安定結合を切る

関連商品

ケムステYoutube

ケムステSlack

注目情報

注目情報

最新記事

第19回ケムステVシンポ「化学者だって起業するっつーの」を開催します!

少し前に化学者のつぶやきからこのような記事が出ました:【ケムステSlackに訊いて見た④】化学系学生…

10種類のスパチュラを試してみた

大好評、「試してみた」シリーズの第6弾。今回は試薬の秤量にか欠かせない、…

第48回「分子の光応答に基づく新現象・新機能の創出」森本 正和 教授

久々の研究者へのインタビューです。第48回は、立教大学の森本正和先生にお願いいたしました。第17回ケ…

畠山琢次 Takuji Hatakeyama

畠山琢次 (はたけやま たくじ)は、日本の化学者である。専門は有機合成化学,材料化学。2021年現在…

DNA origami入門 ―基礎から学ぶDNAナノ構造体の設計技法―

(さらに…)…

NBSでのブロモ化に、酢酸アンモニウムをひとつまみ

芳香環のブロモ化といえば、構造活性相関の取得はもちろんの事、カップリング反応の足場と…

森本 正和 Masakazu Morimoto

森本 正和(もりもと まさかず、MORIMOTO Masakazu)は、日本の化学者である。専門は有…

「リジェネロン国際学生科学技術フェア(ISEF)」をご存じですか?

近年、中高生向けの科学プログラムやコンテストがいっそうの充実を見せています。未来の化学者育成に少なか…

Chem-Station Twitter

PAGE TOP