[スポンサーリンク]

一般的な話題

計算化学: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 で書きます。

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

ゼロ

ゼロ

投稿者の記事一覧

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

関連記事

  1. 3回の分子内共役付加が導くブラシリカルジンの網羅的全合成
  2. Carl Boschの人生 その4
  3. おっさんマウスが小学生マウスを襲う?待ったの決め手はフェロモンに…
  4. アルミニウムで水素分子を活性化する
  5. 力をかけると塩酸が放出される高分子材料
  6. 三次元アクアナノシートの創製! 〜ジャイロイド構造が生み出す高速…
  7. 非天然アミノ酸合成に有用な不斉ロジウム触媒の反応機構解明
  8. 金属内包フラーレンを使った分子レーダーの創製

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

注目情報

ピックアップ記事

  1. 当量と容器サイズでヒドロアミノアルキル化反応を制御する
  2. 三菱化学グループも石化製品を値上げ、原油高で価格転嫁
  3. タンニンでさび防ぐ効果 八王子の会社
  4. コーリー・バクシ・柴田還元 Corey-Bakshi-Shibata (CBS) Reduction
  5. 第117回―「感染症治療を志向したケミカルバイオロジー研究」Erin Carlson准教授
  6. 2017年ケムステ人気記事ランキング
  7. ウーリンス試薬 Woollins’ Reagent
  8. PdとTiがVECsの反応性をひっくり返す?!
  9. 投票!2016年ノーベル化学賞は誰の手に??
  10. 「遠隔位のC-H結合を触媒的に酸化する」―イリノイ大学アーバナ・シャンペーン校・M.C.White研より

関連商品

ケムステYoutube

ケムステSlack

月別アーカイブ

2016年7月
« 6月   8月 »
 123
45678910
11121314151617
18192021222324
25262728293031

注目情報

注目情報

最新記事

ポンコツ博士の海外奮闘録 外伝② 〜J-1 VISA取得編〜

ポンコツシリーズ番外編 その2 J-1 VISA取得までの余談と最近日本で問題になった事件を経験した…

結合をアリーヴェデルチ! Agarozizanol Bの全合成

セスキテルペンAgarozizanol Bの全合成が初めて達成された。光照射下で進行するカスケード反…

有機合成化学協会誌2022年1月号:無保護ケチミン・高周期典型金属・フラビン触媒・機能性ペプチド・人工核酸・脂質様材料

有機合成化学協会が発行する有機合成化学協会誌、2022年1月号がオンライン公開されました。本…

第167回―「バイオ原料の活用を目指した重合法の開発」John Spevacek博士

第167回の海外化学者インタビューは、ジョン・スペヴァセック博士です。Aspen Research社…

繊維強化プラスチックの耐衝撃性を凌ぐゴム材料を開発

名古屋大学大学院工学研究科有機・高分子化学専攻の 野呂 篤史講師らの研究グループは、日本ゼオンと共同…

反応化学の活躍できる場を広げたい!【ケムステ×Hey!Labo 糖化学ノックインインタビュー②】

2021年度科学研究費助成事業 学術変革領域研究(B)に採択された『糖鎖ケミカルノックインが拓く膜動…

UiO-66: 堅牢なジルコニウムクラスターの面心立方格子

UiO-66 は六核ジルコニウムオキシクラスターを SBU に持ち、高い熱安定性 · 化学安定性を示…

危ない試薬・面倒な試薬の便利な代替品

実験室レベルでは、未だに危険な試薬を扱わざるを得ない場合も多いかと思います。tert…

赤外線の化学利用:近赤外からテラヘルツまで

(さらに…)…

【誤解してない?】4s軌道はいつも3d軌道より低いわけではない

3d 遷移金属は、多くが (3d)n(4s)2 という中途半端に 3d 軌道が埋まったまま 4s 軌…

Chem-Station Twitter

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

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

PAGE TOP