[スポンサーリンク]

一般的な話題

ゼロから始める!量子化学計算~遷移状態を求める~

[スポンサーリンク]

 

今回はGaussian及び、GaussViewを用いて、マロンアルデヒドの分子内水素移動反応における遷移状態を求める方法を紹介します。

なお本記事は、大学の先生方の力を借りて、初めて量子化学計算を行った学部生の備忘録的な記事です。

 

こんな人に読んでほしい

・Gaussianのライセンスを持っていて自由に使えるが、何をしたら良いかわからない人

・とりあえず量子化学計算をやってみたい人

計算の流れ

今回遷移状態を求める反応はこちらです。

まず初めに、GaussViewでマロンアルデヒドの構造を作成します。このとき、ベンゼン環から構造を作ると比較的簡単に構造を作ることができます。

また、分子内で移動する水素は遷移状態と思われる位置に後から置きます。このとき実際の遷移状態に近ければ近いほど計算の精度が上がります。具体的には、環のやや内側寄りに置くと良いと思います。O-H結合は後から加えます。

次に計算のインプットファイルを作成します。Job TypeはOptimization+Frequencyにします。ここでOptimizationは構造最適化計算、Frequencyは振動数計算のことです。「Optimize to a」ではTS(Berny)を選択し、遷移状態を求めるよう指示します。「Calculate Force Constants」ではOnceを選択し、計算方法ではHartree-Fockを選択します。以上の設定で計算を行い、遷移状態を求めます。

次に先程の計算結果(logファイル)でresultのvibrationsを見て、虚の振動数(Freqの値が負)があることを確認します。虚の振動数があることは、得られた構造が局所的な安定構造 つまり遷移状態である可能性を意味します。虚の振動数を選択した状態でさらに計算を行います。

次の計算ではJob TypeでIRCを選択します。これは遷移状態が正しい反応物と生成物(ここでは水素が移動する前後のことを指す)までつながっているかを確認する計算のことです。また、Compute more pointsを選択すると計算する範囲を変えることができます。初期設定は10になっています。今回は30にして計算を行いましょう。次にGuessとSolvationをDefaultに、Additional Keywordsに何か入っている場合があるのでそこを空欄にして計算を行います。この計算結果から以下のグラフが得られます。

このグラフは分子全体のエネルギーの大きさを表したもので、大きいほどその構造が不安定、小さいほどその構造が安定していることを示しています。極大値が遷移状態です。また、遷移状態のマロンアルデヒドの構造をみると、水素が2つの酸素から等しい距離に位置しており、直観にもあう遷移状態になっています。

このグラフの下にもう一つグラフが出てきます。これは簡単に言うと一つ目のグラフを微分したもの(変数をそれぞれ微分したものの平均)であり、縦軸が0の点は一つ目のグラフにおいて極値であることを表します。

応用編

マロンアルデヒド分子内水素移動の反応式は対照的であるため、グラフも線対称なものが得られます。炭素鎖についている水素を別の原子に置換して計算しても面白いかもしれません。

画像 1 水素を一つフッ素に置換したもの

片側のホルミル基の水素をフッ素に置換したことで、グラフの概形が対称ではなくなっています。

 

裏話

私はこの計算をするために何度も大学の先生に質問をしました。普通の電卓と違って一筋縄ではいきませんね…。備忘録的な感覚で記事を書かせてもらいました。計算をやってみたいけど、何から手を付ければよいのかわからない、と困っている人の役に立ったら幸いです。

 

量子化学計算についてもっと詳しく知りたい人におすすめの本

[amazonjs asin=”4061543881″ locale=”JP” tmpl=”Small” title=”新版 すぐできる 量子化学計算ビギナーズマニュアル (KS化学専門書)”]

 

Avatar photo

なみぃ

投稿者の記事一覧

学部3年、化学系の勉強をしています。

関連記事

  1. 3回の分子内共役付加が導くブラシリカルジンの網羅的全合成
  2. AIと融合するバイオテクノロジー|越境と共創がもたらす革新的シン…
  3. Micro Flow Reactorで瞬間的変換を達成する
  4. DNAを切らずにゲノム編集-一塩基変換法の開発
  5. ChemDrawの使い方【作図編②:触媒サイクル】
  6. ドーパミンで音楽にシビれる
  7. 配位子だけじゃない!触媒になるホスフィン
  8. 英会話とプログラミングの話

注目情報

ピックアップ記事

  1. BASF International Summer Courses 2017  BASFワークショップ2017
  2. 人物でよみとく化学
  3. 【日産化学 23卒/Zoomウェビナー配信!】START your chemi-story あなたの化学をさがす 研究職限定 キャリアマッチングLIVE
  4. アニリン類のC–N結合に不斉炭素を挿入する
  5. ハロルド・クロトー Harold Walter Kroto
  6. 海外学会出張でeSIMを使ってみました
  7. 小山 靖人 Yasuhito Koyama
  8. リケラボとコラボして特集記事を配信します
  9. 就職か進学かの分かれ道
  10. edXで京都大学の無料講義配信が始まる!

関連商品

ケムステYoutube

ケムステSlack

月別アーカイブ

2022年11月
 123456
78910111213
14151617181920
21222324252627
282930  

注目情報

最新記事

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

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

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

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

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

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

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

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

OIST Science Challenge 2025 に参加しました

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

ペーパークラフトで MOFをつくる

第650回のスポットライトリサーチには、化学コミュニケーション賞2024を受賞された、岡山理科大学 …

月岡温泉で硫黄泉の pH の影響について考えてみた 【化学者が行く温泉巡りの旅】

臭い温泉に入りたい! というわけで、硫黄系温泉を巡る旅の後編です。前回の記事では群馬県草津温泉をご紹…

二酸化マンガンの極小ナノサイズ化で次世代電池や触媒の性能を底上げ!

第649回のスポットライトリサーチは、東北大学大学院環境科学研究科(本間研究室)博士課程後期2年の飯…

日本薬学会第145年会 に参加しよう!

3月27日~29日、福岡国際会議場にて 「日本薬学会第145年会」 が開催されま…

TLC分析がもっと楽に、正確に! ~TLC分析がアナログからデジタルに

薄層クロマトグラフィーは分離手法の一つとして、お金をかけず、安価な方法として現在…

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

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