[スポンサーリンク]

一般的な話題

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

[スポンサーリンク]

 

今回は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. 胃薬のラニチジンに発がん性物質混入のおそれ ~簡易まとめ
  2. Dead Endを回避せよ!「全合成・極限からの一手」⑦(解答編…
  3. 流れる電子ッ!壊れるピリジンッ!含窒素多環式骨格構築!
  4. 可視光レドックス触媒を用いた芳香環へのC-Hアミノ化反応
  5. 原子半径・電気陰性度・中間体の安定性に起因する課題を打破〜担持N…
  6. 第99回日本化学会年会 付設展示会ケムステキャンペーン Part…
  7. 有機反応を俯瞰する ーシグマトロピー転位
  8. アステラス病態代謝研究会 2019年度助成募集

注目情報

ピックアップ記事

  1. グルコース (glucose)
  2. ケムステイブニングミキサー2024に参加しよう!
  3. 有機スペクトル解析ワークブック
  4. オペレーションはイノベーションの夢を見るか? その3+まとめ
  5. 第36回 生体を模倣する化学― Simon Webb教授
  6. 中学生の研究が米国の一流論文誌に掲載された
  7. クロム(η6-アレーン)カルボニル錯体 Cr(η6-arene)(CO)3 Complex
  8. 第99回日本化学会年会 付設展示会ケムステキャンペーン Part III
  9. 印象に残った天然物合成 2
  10. アカデミックから民間企業への転職について考えてみる

関連商品

ケムステYoutube

ケムステSlack

月別アーカイブ

2022年11月
 123456
78910111213
14151617181920
21222324252627
282930  

注目情報

最新記事

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

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…

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

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

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

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