[スポンサーリンク]

一般的な話題

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

[スポンサーリンク]

 

今回は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. まっすぐなペプチドがつまらないなら「さあ輪になって踊ろ!」
  3. 渡辺化学工業ってどんな会社?
  4. 中性ケイ素触媒でヒドロシリル化
  5. MOF 内の水分子吸着過程の解析とそれに基づく水蒸気捕集技術の向…
  6. 来年の応募に向けて!:SciFinder Future Lead…
  7. 音声読み上げソフトで書類チェック
  8. 第三回 ケムステVシンポ「若手化学者、海外経験を語る」を開催しま…

注目情報

ピックアップ記事

  1. 【スポットライトリサーチ】汎用金属粉を使ってアンモニアが合成できたはなし
  2. edXで京都大学の無料講義配信が始まる!
  3. 自己会合・解離機構に基づく蛍光応答性プローブを用いたエクソソーム高感度検出
  4. 鈴木 章 Akira Suzuki
  5. 第二回触媒科学国際シンポジウム
  6. ゲルのやわらかさの秘密:「負のエネルギー弾性」を発見
  7. タンパク質機能をチロシン選択的な修飾で可逆的に制御する
  8. 超原子価臭素試薬を用いた脂肪族C-Hアミノ化反応
  9. 計算化学記事まとめ
  10. ネイチャー論文で絶対立体配置の”誤審”

関連商品

ケムステYoutube

ケムステSlack

月別アーカイブ

2022年11月
 123456
78910111213
14151617181920
21222324252627
282930  

注目情報

最新記事

MDSのはなし 骨髄異形成症候群とそのお薬の開発状況 その1

Tshozoです。今回はかなり限定した疾患とそれに対するお薬の開発の中身をまとめておこうと思いま…

第42回メディシナルケミストリーシンポジウム

テーマAI×創薬 無限能可能性!? ノーベル賞研究が拓く創薬の未来を探る…

山口 潤一郎 Junichiro Yamaguchi

山口潤一郎(やまぐちじゅんいちろう、1979年1月4日–)は日本の有機化学者である。早稲田大学教授 …

ナノグラフェンの高速水素化に成功!メカノケミカル法を用いた芳香環の水素化

第660回のスポットライトリサーチは、名古屋大学大学院理学研究科(有機化学研究室)博士後期課程3年の…

第32回光学活性化合物シンポジウム

第32回光学活性化合物シンポジウムのご案内光学活性化合物の合成および機能創出に関する研究で顕著な…

位置・立体選択的に糖を重水素化するフロー合成法を確立 ― Ru/C触媒カートリッジで150時間以上の連続運転を実証 ―

第 659回のスポットライトリサーチは、岐阜薬科大学大学院 アドバンストケミストリー…

【JAICI Science Dictionary Pro (JSD Pro)】CAS SciFinder®と一緒に活用したいサイエンス辞書サービス

ケムステ読者の皆様には、CAS が提供する科学情報検索ツール CAS SciFind…

有機合成化学協会誌2025年5月号:特集号 有機合成化学の力量を活かした構造有機化学のフロンティア

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

ジョセップ・コルネラ Josep Cornella

ジョセップ・コルネラ(Josep Cornella、1985年2月2日–)はスペイン出身の有機・無機…

電気化学と数理モデルを活用して、複雑な酵素反応の解析に成功

第658回のスポットライトリサーチは、京都大学大学院 農学研究科(生体機能化学研究室)修士2年の市川…

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

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

PAGE TOP