close

Gaussianよくあるご質問

SCFやDFTのエネルギーをすべての構成要素に分解するにはどうしたらよいですか?

エネルギーの完全な内訳は、通常の出力の一部ではありませんが、入力に2つの変更を加えることで取得できます。 まず、計算のルート・セクションにExtraLinks=L608を追加してください。 次に、形状指定の後にinputを追加して、どのような計算を行うかをLink 608に伝えます。 必要な情報は、オンラインのG16 IOPs リファレンスにあるIOp(3/74)のリストに載っています。

この例では、追加された情報が別の色で強調されています。

#P b3lyp/6-31g* ExtraLinks=L608
Gaussian Test Job 135: Fe=O perpendicular to ethene, in triplet state.

この場合,B3LYP の計算であるため,値 -5 が使用され,B3LYP に対応する IOp(3/74) の値は -5 です。 他のメソッドでのこのパラメータの適切な値は、IOpsリファレンス、 またはルートカードの#Pで実行した以前のDFT計算の出力から見つけることができます( 以下の部分出力でハイライトされています)。

Hartree-Fockのエネルギーの場合、値は100ですが、これはデフォルト値であるため、出力されません。

リンク608からの出力はどういう意味ですか?

先の入力ファイルからの出力である以下の出力では、重要な部分をハイライトしています。 各用語の説明は、出力例の後に記載しています。

下記は、エネルギーを構成するさまざまな要素です。

  • ET:運動エネルギー
  • EV:電子-原子核引力
  • EJ:クーロン項
  • EK:DFT密度を用いたハートリー・フォック交換法
  • ENuc:原子核-原子核反発
  • ScaHFX:HF交換のスケールファクタ
  • ScaDFX:DFT交換のスケールファクタ
  • ENTVJ:ENuc + ET + EV + EJ
  • Ex:交換項 = EK * ScaHFX + DFT-Exchange * ScaDFX
  • ETot:トータルエネルギー(SCFエネルギー)=ET+EV+EJ+Ex+Ec+ENuc

*注意:DFT交換の値は、EK、ScaHFX、ScaDFX、および上記のEx方程式を使用して取得できます。