Gaussianのトラブルシューティング④
Gaussianを用いた量子化学計算では、計算条件の設定ミスや構造データの不備などにより、さまざまなエラーが発生することがあります。そこで本記事では、Gaussianを使用する際によく遭遇するエラーを例に、その原因と対処法について解説します。
これからGaussianを使い始める方はもちろん、日常的に計算を行っている方のトラブルシューティングの参考になれば幸いです。
※計算エラーの原因は系や計算条件によって異なります。本記事で紹介する内容はあくまで一例であり、必ずしも全てのケースで解決できるものではありません。
エラーメッセージ
Convergence failure.
Error termination via Lnk1e in /usr/local/intel/g16/l508.exe at Wed May 13 17:45:27 2026.
原因
このエラーは、SCF（Self-Consistent Field）計算が収束しなかった場合に発生します。
Gaussianでは電子密度が一定の収束条件を満たすまで反復計算を行いますが、何らかの理由で収束しない場合、
Convergence failure.
というエラーが出力され、計算が終了します。
特に以下のような系で発生しやすい傾向があります。
・遷移金属錯体
・開殻系（ラジカル）
・電荷を持つ分子
・金属クラスター
・初期構造が不適切な場合
・HOMO-LUMOギャップが小さい系
また、構造最適化の途中で分子構造が大きく変化した場合にも発生することがあります。
SCF収束エラーは、Gaussianで遭遇する機会の多いエラーの一つです。
対処法
SCF収束エラーが発生した場合は、まずアウトプットファイルを確認し、SCF計算がどのような挙動を示しているかを把握してください。大きく分けると、以下のようなパターンがよく見られます。
〇パターン①：SCFが徐々に収束しているが収束しきれていない
この場合は、SCF計算自体は順調に進行しているものの、
・最大反復回数（MaxCycle）に達した
・収束判定まであと少しだった
などの可能性があります。
対処法としては、
scf=(maxcycle=512)
のようにSCFの最大反復回数を増やすことで、正常に収束する場合があります。
〇パターン②：SCFエネルギーが大きく振動している
この場合は、電子状態が安定せず、SCF計算が振動している可能性があります。
特に、
・遷移金属錯体
・開殻系（ラジカル）
・金属クラスター
・HOMO-LUMOギャップが小さい系
で発生しやすい傾向があります。
このような場合は、
scf=(xqc,maxcycle=512)
や
scf=(fermi,xqc,maxcycle=512)
を指定することで改善することがあります。
〇パターン③：何度試しても収束しない
SCFの設定を変更しても改善しない場合は、入力構造や計算条件そのものに原因がある可能性があります。
特に、
・不自然な結合長や原子配置
・電荷（Charge）の設定ミス
・スピン多重度（Multiplicity）の設定ミス
・初期構造が平衡構造から大きく離れている
といった問題がないか確認してください。
また、一度計算レベルを変えて構造最適化を行ってから再計算することで、SCFが収束するケースもあります。
よく使うキーワード
scf=(maxcycle=N)
SCF計算の最大反復回数を N 回に設定します。収束まであと少しという場合に有効です。
scf=xqc
通常のSCFアルゴリズムで収束しない場合に、より強力な収束アルゴリズムを使用します。
scf=fermi
電子占有数を平滑化してSCF収束を助けます。特に遷移金属錯体や金属クラスターで有効な場合があります。
guess=mix
初期軌道の生成時にα軌道とβ軌道を混合します。
主に
・開殻系
・ラジカル
・遷移金属錯体
で有効です。
int=ultrafine
数値積分グリッドを細かくするオプションです。特にM06系汎関数や遷移金属を含む系では、SCFや構造最適化の安定化に役立つことがあります。
関連記事
Gaussianのトラブルシューティングシリーズ
・Gaussianのトラブルシューティング①
・Gaussianのトラブルシューティング②
・Gaussianのトラブルシューティング③
お問い合わせ
量子化学計算、分子シミュレーション、計算科学に関する技術的なご相談や受託解析のご依頼がございましたら、お気軽にお問い合わせください。
計算条件の検討、結果の解釈、反応機構解析、材料開発・創薬研究への活用など、幅広くサポートいたします。
