[スポンサーリンク]

一般的な話題

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

[スポンサーリンク]

 

今回は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 水素を一つフッ素に置換したもの

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

 

裏話

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

 

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

 

なみぃ

投稿者の記事一覧

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

関連記事

  1. ベンジル位アセタールを選択的に酸素酸化する不均一系触媒
  2. スルホキシイミンを用いた一級アミン合成法
  3. 研究室でDIY! ~明るい棚を作ろう~
  4. アメリカ大学院留学:卒業後の進路とインダストリー就活(1)
  5. 異なる“かたち”が共存するキメラ型超分子コポリマーを造る
  6. 研究室での英語【Part 3】
  7. 研究室での英語【Part 2】
  8. 糖のC-2位アリール化は甘くない

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

注目情報

ピックアップ記事

  1. 求電子的トリフルオロメチル化 Electrophilic Trifluoromethylation
  2. スポットライトリサーチムービー:動画であなたの研究を紹介します
  3. 「ヨーロッパで修士号と博士号を取得する」 ―ETH Zürichより―
  4. ケムステ版・ノーベル化学賞候補者リスト【2020年版】
  5. 総合化学大手5社の前期、4社が経常減益
  6. 立体選択的なスピロ環の合成
  7. シンクロトロンで実験してきました【アメリカで Ph.D. を取る: 研究の非日常の巻】
  8. 吉田善一 Zen-ichi Yoshida
  9. トリス(ペンタフルオロフェニル)ボラン : Tris(pentafluorophenyl)borane
  10. ケミストリ・ソングス【Part 2】

関連商品

ケムステYoutube

ケムステSlack

月別アーカイブ

2022年11月
 123456
78910111213
14151617181920
21222324252627
282930  

注目情報

最新記事

塩基が肝!シクロヘキセンのcis-1,3-カルボホウ素化反応

ニッケル触媒を用いたシクロヘキセンの位置および立体選択的なカルボホウ素化反応が開発された。用いる塩基…

中国へ行ってきました 西安・上海・北京編①

2015年(もう8年前ですね)、中国に講演旅行に行った際に記事を書きました(実は途中で断念し最後まで…

アゾ重合開始剤の特徴と選び方

ラジカル重合はビニルモノマーなどの重合に用いられる方法で、開始反応、成長反応、停止反応を素反応とする…

先端事例から深掘りする、マテリアルズ・インフォマティクスと計算科学の融合

開催日:2023/12/20 申込みはこちら■開催概要近年、少子高齢化、働き手の不足の影…

最新の電子顕微鏡法によりポリエチレン分子鎖の向きを可視化することに成功

第583回のスポットライトリサーチは、東北大学大学院 工学研究科 応用化学専攻 陣内研究室の狩野見 …

\脱炭素・サーキュラーエコノミーの実現/  マイクロ波を用いたケミカルリサイクル・金属製錬プロセスのご紹介

※本セミナーは、技術者および事業担当者向けです。脱炭素化と省エネに貢献するモノづくり技術の一つと…

【書籍】女性が科学の扉を開くとき:偏見と差別に対峙した六〇年 NSF(米国国立科学財団)長官を務めた科学者が語る

概要米国の女性科学者たちは科学界のジェンダーギャップにどのように向き合い,変えてきたのか ……

【太陽ホールディングス】新卒採用情報(2025卒)

■■求める人物像■■「大きな志と好奇心を持ちまだ見ぬ価値造像のために前進できる人…

細胞代謝学術セミナー全3回 主催:同仁化学研究所

細胞代謝研究をテーマに第一線でご活躍されている先生方をお招きし、同仁化学研究所主催の学術セミナーを全…

マテリアルズ・インフォマティクスにおける回帰手法の基礎

開催日:2023/12/06 申込みはこちら■開催概要マテリアルズ・インフォマティクスを…

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

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

PAGE TOP