GaussViewを用いた計算ファイルの作成③
量子化学計算ソフト Gaussian を利用する際、最初のハードルとなるのが計算ファイル（Input File）の作成です。Gaussian の入力ファイルはテキスト形式で作成できますが、計算条件や分子構造を一から記述するのは初心者にとって少し難しく感じられるかもしれません。
そこで活躍するのが GaussView です。GaussView を利用すると、分子構造の作成から計算条件の設定までをGUI上で直感的に操作でき、Gaussian の入力ファイルを簡単に作成できます。
本記事では、分子構造の作成から 計算のキーワードの設定までを解説します。これから Gaussian を使い始める方や、入力ファイルの作成方法を基礎から学びたい方はぜひ参考にしてください。
振動計算の設定
GaussViewを用いた計算ファイルの作成②で作成したアセトアルデヒドの構造最適化計算が終了したら、続いて振動計算を行います。
構造最適化によって得られた構造が真の安定構造（極小構造）であることを確認するためには、振動数計算が必要です。安定構造であれば虚振動（負の振動数）は現れず、すべての振動数が正の値となります。
今回は構造最適化が完了した構造を用いるため、構造最適化は行わず、振動計算（Freq）のみを設定します。
まず、GaussViewで構造最適化後のファイルを開き、[Calculate] → [Gaussian Calculation Setup...] をクリックします。
次に、[Job Type] タブを開き、[Optimization] を [Frequencies] に変更します。
なお、Gaussianでは構造最適化と振動数計算を [Opt+Freq] として一度に実行することも可能です。しかし、構造最適化が正常に終了している場合は、振動数計算のみを実行した方が計算時間を短縮できます。また、構造最適化が収束しなかった場合でも、振動数計算だけを再実行できるため、個別に計算を設定する方法を覚えておくと便利です。
計算手法（Method）や基底関数（Basis Set）は、構造最適化と同じ条件を使用してください。異なる条件を指定すると、最適化した構造と振動計算の計算レベルが一致しなくなるため、通常は同じ条件で計算を行います。
設定が完了したら、[Link 0] タブでメモリ容量やCPUコア数を確認し、必要に応じて変更します。最後に [Submit] をクリックすると、振動計算が開始されます。
GaussViewで振動モードとIRスペクトルを確認する
振動数計算が終了したら、GaussViewを用いて各振動モードやIRスペクトルを確認できます。
まず、振動数計算を実行したログファイル（.log または .out）をGaussViewで開きます。
続いて、メニューから [Results] → [Vibrations] をクリックします。
すると、振動数計算で得られた各振動モードの一覧が表示されます。
各振動モードを選択して [Start Animation] をクリックすると、原子の振動の様子をアニメーションで確認できます。例えば、C-H伸縮振動やC=O伸縮振動など、どの原子がどのように動いているかを視覚的に把握できます。
さらに、同じ画面から[Spectra]を選択するとIRスペクトルを表示することもできます。
IRスペクトルでは、横軸に波数（cm⁻¹）、縦軸にIR強度が表示されます。ピークを選択すると対応する振動モードが表示されるため、どの官能基の振動に由来するピークなのかを確認できます。
実験で測定したIRスペクトルと比較することで、分子構造の妥当性や官能基の帰属を検討することができます。
例えば、アセトアルデヒドでは以下のような振動が観測されます。
-約1700 cm⁻¹付近：C=O伸縮振動
-2800～3100 cm⁻¹付近：C-H伸縮振動
-1000～1500 cm⁻¹付近：変角振動や骨格振動
ただし、計算で得られる振動数は理論値であり、実験値とは若干のずれが生じることがあります。IRスペクトルは理論値として出力されるため、そのまま実験値と完全一致するわけではありません。特にDFT計算では、一般に振動数が系統的に高めに出る傾向があります。
そのため、実務的にはスケーリングファクター（scaling factor）を適用して補正した上で比較を行います。例えばB3LYP/6-31G(d)レベルでは、0.96前後の係数が用いられることが多く、これを掛けることで実験スペクトルとの整合性が改善されます。
また、ピーク強度についても注意が必要で、計算上のIR強度は相対的な指標であり、実験の吸光度と直接一致するものではありません。そのため、ピーク位置の帰属が主な目的となります。
アセトアルデヒドの場合であれば、C=O伸縮振動が最も強いピークとして現れることが多く、この特徴は構造解析において非常に重要な指標となります。
振動計算で分かること
ここまでで振動モードとIRスペクトルの確認方法について説明しましたが、振動数計算はそれ以外にも重要な情報を与えてくれます。
主な用途は以下のとおりです。
〇安定構造・遷移状態の判定
・虚振動が0個：安定構造（極小構造）
・虚振動が1個：遷移状態の候補
・虚振動が2個以上：構造が収束していない、または高次の鞍点である可能性
〇IRスペクトルの予測
計算から得られた振動数とIR強度を用いて、理論IRスペクトルを作成できます。実験スペクトルとの比較やピーク帰属に利用されます。
〇熱力学量の計算
振動数計算からは、ゼロ点振動エネルギー（ZPE）やエンタルピー、エントロピー、ギブズ自由エネルギーなども算出されます。反応エネルギーや反応自由エネルギーの評価には欠かせない計算です。
〇後続計算への利用
振動数計算の結果は、遷移状態探索（TS）やIRC計算の妥当性確認にも利用されます。特に遷移状態では、虚振動の方向が目的とする反応経路に対応していることを確認することが重要です。
関連記事
Gaussianのトラブルシューティングシリーズ
・Gaussianのトラブルシューティング①
・Gaussianのトラブルシューティング②
・Gaussianのトラブルシューティング③
・Gaussianのトラブルシューティング④
GaussViewを用いた計算ファイルの作成シリーズ
・GaussViewを用いた計算ファイルの作成①
・GaussViewを用いた計算ファイルの作成②
お問い合わせ
量子化学計算、分子シミュレーション、計算科学に関する技術的なご相談や受託解析のご依頼がございましたら、お気軽にお問い合わせください。
計算条件の検討、結果の解釈、反応機構解析、材料開発・創薬研究への活用など、幅広くサポートいたします。
