[スポンサーリンク]

一般的な話題

計算化学:DFTって何? PartIII

[スポンサーリンク]

前回の記事では、1930年前後から1980年頃までの量子化学の発展について簡単に説明しました。

今回は、DFT計算で使われる汎関数の種類、gaussianでの計算の流れなどを説明します。

前回と今回の記事を読むことにより、DFT計算の歴史的背景、理論、問題点、汎関数の種類、gaussianでの計算の流れなどが理解できます。これで、「DFTって何?」と聞かれても答えられるようになると思います。

今回の記事でも細かいところは省略していくので、興味のある方は一番下にある参考図書を読んで下さいね!

1930年前後から1980年頃までの歴史において、前回触れなかったけれど、重要なものとして、基底関数の開発(1950年Boys)、ローターン方程式(1950年Roothaan, Hall)などがあります。また、現在の計算化学に必須と言っても過言ではないGaussian(1970年、Pople)、Gamess(1982年)などのプログラムの開発もこの頃です。

 

様々な汎関数

以前の記事でお話ししたように、交換相関汎関数はコーン・シャム方程式の中で唯一近似されている部分です。よって、DFT計算の信頼性は、どの汎関数を選ぶかで決まります。これまで、星の数ほどの汎関数が開発されてきましたが、それらは次のように分類できます。

 

局所密度近似(LDA)汎関数
一般化勾配近似(GGA)汎関数
メタGGA汎関数
混成汎関数(ハイブリッド汎関数)
半経験的汎関数
プログレッシブ汎関数

 

シュレーディンガー方程式とは違い、交換相関汎関数の理想的な式の形は未知ですので、それぞれの研究者が、それぞれ好きなパラメーターで汎関数を作りました。
ある人は、交換相関汎関数は電子密度に比例すると言い、またある人は電子密度の微分に比例すると言った、、、というようにどの変数を使うか、ということで分類されています。さらには、幾つかの汎関数を組み合わせたものが良い(ハイブリッド汎関数)という人まで出てきて、カオスな感じです。
個々の汎関数の詳しい内容については、こちらの参考書を読むことをオススメします。式の成り立ちを化学者にもわかりやすく書いてあります。
今回は、有機合成化学者がよく使う代表的な汎関数を 2 つ紹介します。

まずは、B3LYPですっ!

B3LYP

B3LYP は、最初の混成汎関数であり、有機化学の世界でもっともよく使われている汎関数です(でした)。1993年に Becke により発表されました。

この汎関数にはハートリーフォック交換積分の混合率だけでなく、B88交換汎関数、LYP相関汎関数、そしてLDA交換・相関汎関数からの差のそれぞれの混合率の計三個のパラメータが含まれています。

B3LYP

半経験的パラメータは、a1 = 0.2、a2 = 0.72、a3 = 0.81 の 3 つです。これらのパラメーターをどのようにして決めたかというと、ベンチマークと言われる指標で、実験値と計算値が近くなるように設定されました。なので、数学的な意味はありません。こんなことを書くと、「なんだよ!DFT 計算って不正確じゃん。騙したナァァァ!」って思うかもしれませんが、そんなことはありません。むしろ、シュレーディンガー方程式を近似的に解く MP2 法などよりも計算結果と実験結果が一致します。

そもそも、

ヒトはなぜ計算化学を利用するか?

という、ある種、人生の根源的な問いについて考えてみると、計算結果と実験結果が一致するからに他ならないという解が得られるはずです。数学的にいくら正しくても、実験結果との誤差が大きければ誰も使いません。同様のことが、酵素や遺伝子配列などの相同性検索にも言えると思います。

機能が近ければ配列も近いはず!

という仮定は、進化的にも直感的にも正しいと感じるかもしれませんが、厳密に数学的に証明されたわけではありません。なぜ、研究者が相同性検索を使うかといえば、実験結果と一致するからです!!DFT 計算も同様です。

M06系

M062-X_Donald

また、最近よく使われ始めているのは、ドナルド らが開発した M06 系だと思います。こちらは半経験的汎関数に分類されており、38個の半経験的パラメータが入っています。しかし、M06が今後ずっと主流になるとは考えられてはおらず、このような半経験的汎関数は洗練された汎関数が開発されるまでの過渡的な汎関数と考えられています

遷移金属の入っている系には  M06、遷移金属の入っていない系には  M06-2x  が使われているような印象を受けます。また、弱い相互作用を見積もることが可能な汎関数ですので、ペリ環状反応などでは、B3LYP で構造最適化したのちに M062X  で一点計算している例をよく見かけます。

 

実際の計算の流れ

最近の計算プログラムは非常に使いやすくなっているため、まったく背景知識が無い人でも、化合物のエネルギーを計算したり、構造最適化を行なったり、スペクトルの計算が出来てしまいます。例えるならば、テレビゲーム

 

とりあえずインプットファイルを作り計算させれば、結果は出てきます。しかし、その結果が正しいのか?どのような意味があるか?を判断するには有機化学の知識が必要となります。また、エラーが起きた時に正しく対処するためには、背景にある理論や計算プログラムを良く理解しておく必要があります。

計算化学って言うけど、コンピューターが実際にどのような計算を行なってくれるのか答えられる人は少ないです。
例えば、分子構造最適化の計算手順は、大まかに以下のようになっています。

1.座標の読み込み
2.エネルギーの計算(Roothaan 方程式 SCF 計算)
3.力の計算
4.構造最適化

計算終了後のlog ファイルを見ると、どの部分が律速か、どのパラメーターを変えれば早くなるかなどが分かるかと思いますが、長くなってしまったので、詳しいことは次回、 計算化学:DFTって何? PartIV で書きます。

ここまで読んでいただいて、ありがとうございました。

Avatar photo

ゼロ

投稿者の記事一覧

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

関連記事

  1. あなたの天秤、正確ですか?
  2. CSJジャーナルフォーラム「ジャーナルの将来像を考える」
  3. Nature Chemistry:Research Highli…
  4. 製薬会社のテレビCMがステキです
  5. CO酸化触媒として機能する、“無保護”合金型ナノ粒子を担持した基…
  6. 【技術系スタートアップ合同フォーラムのお知らせ】 ディープテック…
  7. 触媒と光で脳内のアミロイドβを酸素化
  8. メチオニン選択的なタンパク質修飾反応

注目情報

ピックアップ記事

  1. 山口マクロラクトン化 Yamaguchi Macrolactonizaion
  2. 不斉をあざ(Aza)やかに(Ni)制御!Aza-Heck環化/還元的カップリング
  3. Independence Day
  4. 【マイクロ波化学(株) 石油化学/プラスチック業界向けウェビナー】 マイクロ波による新事業 石油化学・プラスチック業界のための脱炭素・電化ソリューション
  5. ジョン・ケンドリュー John C. Kendrew
  6. Nazarov環化を利用した全合成研究
  7. 再生医療ーChemical Times特集より
  8. 健康食品 高まる開発熱 新素材も続々
  9. 植物毒素の全合成と細胞死におけるオルガネラの現象発見
  10. ベシャンプ還元 Bechamp Reduction

関連商品

ケムステYoutube

ケムステSlack

月別アーカイブ

2016年7月
 123
45678910
11121314151617
18192021222324
25262728293031

注目情報

最新記事

創薬懇話会2025 in 大津

日時2025年6月19日(木)~6月20日(金)宿泊型セミナー会場ホテル…

理研の研究者が考える未来のバイオ技術とは?

bergです。昨今、環境問題や資源問題の関心の高まりから人工酵素や微生物を利用した化学合成やバイオテ…

水を含み湿度に応答するラメラ構造ポリマー材料の開発

第651回のスポットライトリサーチは、京都大学大学院工学研究科(大内研究室)の堀池優貴 さんにお願い…

第57回有機金属若手の会 夏の学校

案内:今年度も、有機金属若手の会夏の学校を2泊3日の合宿形式で開催します。有機金…

高用量ビタミンB12がALSに治療効果を発揮する。しかし流通問題も。

2024年11月20日、エーザイ株式会社は、筋萎縮性側索硬化症用剤「ロゼバラミン…

第23回次世代を担う有機化学シンポジウム

「若手研究者が口頭発表する機会や自由闊達にディスカッションする場を増やし、若手の研究活動をエンカレッ…

ペロブスカイト太陽電池開発におけるマテリアルズ・インフォマティクスの活用

持続可能な社会の実現に向けて、太陽電池は太陽光発電における中心的な要素として注目…

有機合成化学協会誌2025年3月号:チェーンウォーキング・カルコゲン結合・有機電解反応・ロタキサン・配位重合

有機合成化学協会が発行する有機合成化学協会誌、2025年3月号がオンラインで公開されています!…

CIPイノベーション共創プログラム「未来の医療を支えるバイオベンチャーの新たな戦略」

日本化学会第105春季年会(2025)で開催されるシンポジウムの一つに、CIPセッション「未来の医療…

OIST Science Challenge 2025 に参加しました

2025年3月15日から22日にかけて沖縄科学技術大学院大学 (OIST) にて開催された Scie…

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

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