説明

Gaussian 16 は、さまざまな密度汎関数理論 (DFT) を提供します [ Hohenberg64 P. Hohenberg and W. Kohn, “Inhomogeneous Electron Gas,” Phys. Rev., 136 (1964) B864-B71. DOI: 10.1103/PhysRev.136.B864 , Kohn65 W. Kohn and L. J. Sham, “Self-Consistent Equations Including Exchange and Correlation Effects,” Phys. Rev., 140 (1965) A1133-A38. DOI: 10.1103/PhysRev.140.A1133 , Parr89 R. G. Parr and W. Yang, Density-functional theory of atoms and molecules (Oxford Univ. Press, Oxford, 1989). , Salahub89 The Challenge of d and f Electrons, Ed. D. R. Salahub and M. C. Zerner (ACS, Washington, D.C., 1989). DOI: 10.1021/bk-1989-0394 ] モデル ([ Labanowski91 Density Functional Methods in Chemistry, Ed. J. K. Labanowski and J. W. Andzelm (Springer-Verlag, New York, 1991). , Andzelm92 J. Andzelm and E. Wimmer, “Density functional Gaussian-type-orbital approach to molecular geometries, vibrations, and reaction energies,” J. Chem. Phys., 96 (1992) 1280-303. DOI: 10.1063/1.462165 , Becke92 A. D. Becke, “Density-functional thermochemistry. I. The effect of the exchange-only gradient correction,” J. Chem. Phys., 96 (1992) 2155-60. DOI: 10.1063/1.462066 , Gill92 P. M. W. Gill, B. G. Johnson, J. A. Pople, and M. J. Frisch, “The performance of the Becke-Lee-Yang-Parr (B-LYP) density functional theory with various basis sets,” Chem. Phys. Lett., 197 (1992) 499-505. DOI: 10.1016/0009-2614(92)85807-M , Perdew92 J. P. Perdew, J. A. Chevary, S. H. Vosko, K. A. Jackson, M. R. Pederson, D. J. Singh, and C. Fiolhais, “Atoms, molecules, solids, and surfaces: Applications of the generalized gradient approximation for exchange and correlation,” Phys. Rev. B, 46 (1992) 6671-87. DOI: 10.1103/PhysRevB.46.6671 , Scuseria92 G. E. Scuseria, “Comparison of coupled-cluster results with a hybrid of Hartree-Fock and density functional theory,” J. Chem. Phys., 97 (1992) 7528-30. DOI: 10.1063/1.463977 , Becke92a A. D. Becke, “Density-functional thermochemistry. II. The effect of the Perdew-Wang generalized-gradient correlation correction,” J. Chem. Phys., 97 (1992) 9173-77. DOI: 10.1063/1.463343 , Perdew92a J. P. Perdew and Y. Wang, “Accurate and Simple Analytic Representation of the Electron Gas Correlation Energy,” Phys. Rev. B, 45 (1992) 13244-49. DOI: 10.1103/PhysRevB.45.13244 , Perdew93a J. P. Perdew, J. A. Chevary, S. H. Vosko, K. A. Jackson, M. R. Pederson, D. J. Singh, and C. Fiolhais, “Erratum: Atoms, molecules, solids, and surfaces - Applications of the generalized gradient approximation for exchange and correlation,” Phys. Rev. B, 48 (1993) 4978. DOI: 10.1103/PhysRevB.48.4978.2 , Sosa93a C. Sosa and C. Lee, “Density-functional description of transition structures using nonlocal corrections: Silylene insertion reactions into the hydrogen molecule,” J. Chem. Phys., 98 (1993) 8004-11. DOI: 10.1063/1.464554 , Stephens94 P. J. Stephens, F. J. Devlin, M. J. Frisch, and C. F. Chabalowski, “Ab initio Calculation of Vibrational Absorption and Circular Dichroism Spectra Using Density Functional Force Fields,” J. Phys. Chem., 98 (1994) 11623-27. DOI: 10.1021/j100096a001 , Stephens94a P. J. Stephens, F. J. Devlin, C. S. Ashvar, C. F. Chabalowski, and M. J. Frisch, “Theoretical Calculation of Vibrational Circular Dichroism Spectra,” Faraday Discuss., 99 (1994) 103-19. DOI: 10.1039/FD9949900103 , Ricca95 A. Ricca and C. W. Bauschlicher Jr., “Successive H2O binding energies for Fe(H2O)N+,” J. Phys. Chem., 99 (1995) 9003-07. DOI: 10.1021/j100022a010 DFT の手法とアプリケーションについては、上記の文献を参照してください。エネルギー [ Pople92 J. A. Pople, P. M. W. Gill, and B. G. Johnson, “Kohn-Sham density-functional theory within a finite basis set,” Chem. Phys. Lett., 199 (1992) 557-60. DOI: 10.1016/0009-2614(92)85009-Y ]、解析勾配、および真の解析周波数 [ Johnson93a B. G. Johnson and M. J. Frisch, “Analytic second derivatives of the gradient-corrected density functional energy: Effect of quadrature weight derivatives,” Chem. Phys. Lett., 216 (1993) 133-40. DOI: 10.1016/0009-2614(93)E1238-C , Johnson94 B. G. Johnson and M. J. Frisch, “An implementation of analytic second derivatives of the gradient-corrected density functional energy,” J. Chem. Phys., 100 (1994) 7429-42. DOI: 10.1063/1.466887 , Stratmann97 R. E. Stratmann, J. C. Burant, G. E. Scuseria, and M. J. Frisch, “Improving harmonic vibrational frequencies calculations in density functional theory,” J. Chem. Phys., 106 (1997) 10175-83. DOI: 10.1063/1.474047 これらはすべて、DFT モデルとして利用できます。

自己矛盾のない反応場 (SCRF) を DFT エネルギー、最適化、および周波数計算とともに使用して、溶液内のシステムをモデル化できます。

純粋な DFT 計算では、密度フィッティングを利用することがよくあります。詳細は 基底関数系 の説明を参照してください。

freqmem によって提示されるものと同じ最適メモリサイズが、DFT の振動解析(周波数計算)に対しても推奨されます。

分極率導関数(ラマン強度)と超分極率は、DFT の振動解析(周波数計算)ではデフォルトでは計算されません。これらを求めるには Freq=Raman を使用します。Polar 計算では、これらの量が計算されます。

Note: 二重ハイブリッド汎関数については、計算コストが MP2 と同程度であるため、MP2 キーワードを使用してください。

精度に関する考慮事項

DFT 計算では、Hartree-Fock 計算の各主要段階に追加のステップが加わります。このステップは、汎関数(またはその各種導関数)の数値積分です。したがって、Hartree-Fock 計算における数値誤差の原因(積分精度、SCF 収束、CPHF 収束)に加えて、DFT 計算の精度は数値積分で使用される格子点の数にも依存します。

UltraFine 積分グリッド(Integral=UltraFine)は Gaussian 16 のデフォルトです。このグリッドは、妥当な追加コストで計算精度を大幅に向上させます。本番用の DFT 計算では、これより粗いグリッドを使用することは推奨されません。また、エネルギー差や生成熱など、エネルギーを比較するすべての計算で、同じグリッドを使用することが重要です。

必要に応じて、より大きなグリッドを使用できます (例: 特定の種類のシステムの厳密な構造の最適化)。代替グリッドは、Integral=Grid ルートセクションのオプションで指定します。

基礎的概念

Background

ハートリー-フォック理論では、エネルギーは次のような形式になります。

EHF = V + ⟨hP⟩ + 1/2⟨PJ(P)⟩ – 1/2⟨PK(P)⟩

ここで、用語は次の意味を持ちます。

V 核反発エネルギー。
P 密度行列。
⟨hP⟩ 一電子エネルギー(運動エネルギー + ポテンシャルエネルギー)。
1/2⟨PJ(P)⟩ 電子間の古典的クーロン反発。
-1/2⟨PK(P)⟩ 電子の量子的性質(フェルミオン性)に由来する交換エネルギー。

密度汎関数理論のコーン・シャム定式化では [ Kohn65 W. Kohn and L. J. Sham, “Self-Consistent Equations Including Exchange and Correlation Effects,” Phys. Rev., 140 (1965) A1133-A38. DOI: 10.1103/PhysRev.140.A1133 ] では、単一の行列式の正確な交換 (HF) は、より一般的な式である交換相関汎関数に置き換えられます。これには、交換エネルギーと電子相関エネルギーの両方を説明する項が含まれますが、後者はハートリー-フォック理論には存在しません。

EKS = V + ⟨hP⟩ + 1/2⟨PJ(P)⟩ + EX[P] + EC[P]

ここで EX[P] は交換汎関数、EC[P] は相関汎関数です。

Kohn-Sham の定式化では、Hartree-Fock 理論は密度汎関数理論の特殊な場合とみなすことができます。この場合、EX[P] は交換積分 -1/2⟨PK(P)⟩ で与えられ、EC=0 です。密度汎関数理論で通常使用される汎関数は、密度、および場合によっては密度勾配の関数を積分したものです。

EX[P] = ∫f(ρα(r),ρβ(r),∇ρα(r),∇ρβ(r))dr

各手法の違いは、EX にどの関数 f を用いるか、また(存在する場合は)EC にどの関数 f を用いるかにあります。純粋 DFT 法に加えて、Gaussian は、交換汎関数が Hartree-Fock 交換と上記形式の関数積分との線形結合であるハイブリッド法をサポートします。提案されている汎関数は閉じた形では評価できず、数値求積によって評価する積分を含みます。

Hybrid関数キーワード

キーワード: ハイブリッド汎関数

いくつかの Hartree-Fock 交換と DFT 交換相関の混合を含むハイブリッド汎関数は、キーワードを介して利用できます。

Becke Three-Parameter Hybrid Functionals

これらの汎関数は、1993 年に Becke によって考案された形式を持っています [ Becke93a A. D. Becke, “Density-functional thermochemistry. III. The role of exact exchange,” J. Chem. Phys., 98 (1993) 5648-52. DOI: 10.1063/1.464913 ]:

A*EXSlater+(1-A)*EXHF+B*ΔEXBecke+ECVWN+C*ΔECnon-local

ここで、A, B, C は、G1 分子セットへのフィッティングを通じて Becke によって決定された定数です。

このハイブリッド汎関数にはいくつかのバリエーションがあります。

B3LYP は、LYP 式による非局所相関と、局所相関には VWN の汎関数 III(汎関数 V ではありません)を使用します。LYP には局所項と非局所項の両方が含まれるため、実際に使用される相関汎関数は次のようになります。

C*ECLYP+(1-C)*ECVWN

言い換えれば、LYP には VWN と本質的に同等のローカル項が含まれるため、VWN は必要な過剰なローカル相関を提供するために使用されます。

B3P86 は、Perdew 86 による非局所相関を用いる同じ形式の汎関数を指定します。B3PW91 は、Perdew/Wang 91 による非局所相関を用いるこの汎関数を指定します。

O3LYP は、B3LYP と同様の 3 パラメーター汎関数です。

A*EXLSD+(1-A)*EXHF+B*ΔEXOPTX+C*ΔECLYP+(1-C)ECVWN

ここで、A、B、および C は、参照 [ Cohen01 A. J. Cohen and N. C. Handy, “Dynamic correlation,” Mol. Phys., 99 (2001) 607-15. DOI: 10.1080/00268970010023435 ].

分散を含む汎関数

スタンドアロンのキーワード EmpiricalDispersion を用いて、さまざまな汎関数に分散補正スキームを指定することもできます。

長距離補正汎関数

交換汎関数の非クーロン部分は通常、長距離で急速に減衰し、非常に不正確になるため、高い軌道への電子励起などの過程をモデル化するには適していません。このような場合に対処するため、さまざまなスキームが考案されています。Gaussian 16 では、長距離補正を含む次の汎関数を提供しています。

さらに、接頭辞 LC- をほとんどの純粋汎関数に付けることで、Hirao らの長距離補正を適用できます [ Iikura01 H. Iikura, T. Tsuneda, T. Yanai, and K. Hirao, “Long-range correction scheme for generalized-gradient-approximation exchange functionals,” J. Chem. Phys., 115 (2001) 3540-44. DOI: 10.1063/1.1383587 ]。例: LC-BLYP.

その他のハイブリッド汎関数

Truhlar グループの汎関数

PBE 相関を用いる汎関数

Becke の 1 パラメーターハイブリッド汎関数

B1B95 キーワードは元の論文で定義されている Becke の 1 パラメーター ハイブリッド汎関数を指定するために使用されます [ Becke96 A. D. Becke, “Density-functional thermochemistry. IV. A new dynamical correlation functional and implications for exact-exchange mixing,” J. Chem. Phys., 104 (1996) 1040-46. DOI: 10.1063/1.470829 ].

このプログラムは、Adamo と Barone によって実装された、他の同様の 1 パラメーターハイブリッド汎関数も提供します。 Adamo97 C. Adamo and V. Barone, “Toward reliable adiabatic connection models free from adjustable parameters,” Chem. Phys. Lett., 274 (1997) 242-50. DOI: 10.1016/S0009-2614(97)00651-9 ]。あるバリエーションでは、B1LYP として、上記の B3LYP で説明した LYP 相関汎関数を使用します。別のバージョン mPW1PW91 では、Adamo と Barone によって修正された Perdew-Wang 交換を PW91 相関と組み合わせて使用します [ Adamo98 C. Adamo and V. Barone, “Exchange functionals with improved long-range behavior and adiabatic connection methods without adjustable parameters: The mPW and mPW1PW models,” J. Chem. Phys., 108 (1998) 664-75. DOI: 10.1063/1.475428 ]; mPW1LYP, mPW1PBE, mPW3PBE の各バリエーションも用意されています。

B97 の改訂版

Functionals with τ-Dependent Gradient-Corrected Correlation

旧来の汎関数

Half-and-Half 汎関数

以下の汎関数は、下位互換性のためだけに含まれています。これらは、Becke が提案した「half-and-half」汎関数 [ Becke93 A. D. Becke, “A new mixing of Hartree-Fock and local density-functional theories,” J. Chem. Phys., 98 (1993) 1372-77. DOI: 10.1063/1.464304 ] と同じではないことに注意してください。

  • BHandH: 0.5*EXHF + 0.5*EXLSDA + ECLYP
  • BHandHLYP: 0.5*EXHF + 0.5*EXLSDA + 0.5*ΔEXBecke88 + ECLYP

ユーザー定義ハイブリッドモデル

Gaussian 16 では、次の一般形式でユーザー定義のハイブリッド汎関数を指定できます。

P2EXHF + P1(P4EXSlater + P3ΔExnon-local) + P6EClocal + P5ΔECnon-local

利用可能な唯一の局所交換法は Slater(S)であり、局所交換だけが必要な場合に使用します。前述のとおり、任意の組み合わせ可能な非局所交換汎関数および相関汎関数を使用できます。

6 つのパラメーターの値は、プログラムのさまざまな非標準オプションを使用して指定されます。

  • IOp(3/76=mmmmmnnnnn) は、P1mmmmm/10000、P2nnnnn/10000 に設定します。P1 は通常、交換汎関数を使用するかどうかに応じて 1.0 または 0.0 に設定し、スケーリングは P3 と P4 で行います。
  • IOp(3/77=mmmmmnnnnn) は、P3mmmmm/10000、P4nnnnn/10000 に設定します。
  • IOp(3/78=mmmmmnnnnn) は、P5mmmmm/10000、P6nnnnn/10000 に設定します。

例えば、IOp(3/76=1000005000) は P1 を 1.0、P2 を 0.5 に設定します。すべての値は、必要な先頭のゼロを付けた 5 桁で表す必要があります。

以下は、B3LYP キーワードに対応する汎関数を指定するルートセクションです。

#P BLYP IOp(3/76=1000002000) IOp(3/77=0720008000) IOp(3/78=0810010000)

出力ファイルには、使用されている値が表示されます。

 IExCor=  402 DFT=T Ex=B+HF Corr=LYP ExCW=0 ScaHFX=  0.200000
                    ScaDFX=  0.800000  0.720000  1.000000  0.810000
                  

ここで、ScaHFX の値は P2 であり、ScaDFX に指定された値の並びは P4, P3, P6, P5 に対応します。

Pure汎関数キーワード

キーワード: Pure汎関数

さまざまな Pure DFT モデルの名前は、交換汎関数と相関汎関数の名前を組み合わせて付けられます。場合によっては、その分野で使われる標準的な同義語もキーワードとして使用できます。Pure 汎関数を指定するには、交換汎関数コンポーネントのキーワードと目的の相関汎関数のキーワードを組み合わせます。たとえば、Becke 交換汎関数(B)と LYP 相関汎関数は、BLYP キーワードで指定します。同様に、SVWN は Slater 交換(S)と VWN 相関汎関数を指定し、文献では LSDA(局所スピン密度近似)という同義語で知られています。LSDASVWN の同義語です。DFT 機能を持つ一部の他のソフトウェアパッケージでは、「LSDA」が要求された場合に SVWN5 相当のものを使用します。比較する場合は、各パッケージのドキュメントを注意深く確認してください。

交換汎関数

次の 交換汎関数は Gaussian 16 で利用できます。特に指定がない限り、使用可能な方法を生成するには、これらの交換汎関数を相関汎関数と組み合わせる必要があります。

相関汎関数

次の 相関汎関数が利用可能で、対応するキーワードコンポーネントごとにリストされています。これらはすべて、目的の交換汎関数のキーワードと組み合わせる必要があります。

相関汎関数のバリエーション。 次の相関汎関数は、さまざまな相関汎関数の局所項と非局所項を組み合わせたものです。

  • VP86: VWN5 の局所項と P86 の非局所相関汎関数を組み合わせたものです。
  • V5LYP: VWN5 の局所項と LYP の非局所相関汎関数を組み合わせたものです。

単独で使用する Pure 汎関数

次の Pure 汎関数は自己完結型であり、他の汎関数キーワードコンポーネントと組み合わせることはありません。

経験的分散力

Dispersion

EmpiricalDispersion キーワードにより経験的な分散が可能になります。次のオプションを取ります。

PFD

APFD 汎関数から Petersson-Frisch 分散モデルを追加します [ Austin12 A. Austin, G. Petersson, M. J. Frisch, F. J. Dobek, G. Scalmani, and K. Throssell, “A density functional with spherical atom dispersion terms,” J. Chem. Theory and Comput. 8 (2012) 4989. DOI: 10.1021/ct300778e ].

GD2

Grimme の分散の D2 バージョンを追加 [ Grimme06 S. Grimme, “Semiempirical GGA-type density functional constructed with a long-range dispersion correction,” J. Comp. Chem., 27 (2006) 1787-99. DOI: 10.1002/jcc.20495 ]。以下の表は、GD2 パラメーターが定義されている Gaussian 16 の汎関数のリストを示しています。太字で強調表示された汎関数には、示されたキーワードが指定された場合にデフォルトでこの分散モデルが含まれます(例: B2PLYPD)。それ以外の汎関数では、EmpiricalDispersion=GD2 によって分散補正を要求します。

Functional S6 SR6
B97D 1.2500 1.1000
B2PLYPD 0.5500 1.1000
mPW2PLYPD 0.4000 1.1000
PBEPBE 0.7500 1.1000
BLYP 1.2000 1.1000
B3LYP 1.0500 1.1000
BP86 1.0500 1.1000
TPSSTPSS 1.0000 1.1000

このモデルで使用される減衰関数には、固定値 6.0 の D6 パラメーターも含まれています。

この経験的分散法は、IOps(3/174,176) を介して他の汎関数で使用できます(SR6 は 1.1 である必要があります)。

wB97XD を独立したキーワードとして指定した場合、この分散モデルのうち S6 と SR6 の値がそれぞれ 1.0 と 1.1 であるバージョンを使用します。この汎関数は、GD3 モデルで使用される減衰関数と同様の減衰関数を使用し、D6 と IA6 の固定値はそれぞれ 6.0 と 12 です。

GD3

Grimme の分散補正の D3 版(オリジナルの D3 減衰関数を使用)を追加します [ Grimme10 S. Grimme, J. Antony, S. Ehrlich and H. Krieg, “A consistent and accurate ab initio parameterization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu,” J. Chem. Phys., 132 (2010) 154104. DOI: 10.1063/1.3382344 ]。以下の表は、GD3 パラメーターが定義されている Gaussian 16 の汎関数のリストを示しています。これらの汎関数では、EmpiricalDispersion=GD3 によって分散補正を要求します。

Functional S6 SR6 S8
B2PLYPD3 [ Goerigk11 L. Goerigk and S. Grimme, “Efficient and Accurate Double-Hybrid-Meta-GGA Density Functionals—Evaluation with the Extended GMTKN30 Database for General Main Group Thermochemistry, Kinetics, and Noncovalent Interactions,” J. Chem. Theory Comput., 7 (2011) 291-309. DOI: 10.1021/ct10046k ] 0.6400 1.4270 1.0220
B97D3 1.0000 0.8920 0.9090
B3LYP 1.0000 1.2610 1.7030
BLYP 1.0000 1.0940 1.6820
PBE1PBE 1.0000 1.2870 0.9280
TPSSTPSS 1.0000 1.1660 1.1050
PBEPBE 1.0000 1.2170 0.7220
BP86 1.0000 1.1390 1.6830
BPBE 1.0000 1.0870 2.0330
B3PW91 1.0000 1.1760 1.7750
BMK 1.0000 1.9310 2.1680
CAM–B3LYP 1.0000 1.3780 1.2170
LC-wPBE 1.0000 1.3550 1.2790
M05 1.0000 1.3730 0.5950
M052X 1.0000 1.4170 0.0000
M06L 1.0000 1.5810 0.0000
M06 1.0000 1.3250 0.0000
M062X 1.0000 1.6190 0.0000
M06HF 1.0000 1.4460 0.0000
PW6B95D3 1.0000 1.532 0.862

このモデルでは、固定値 1.0 の SR8 パラメータも使用します。このモデルで使用される減衰関数には、それぞれ 6.0、14、6.0、および 16 の固定値を持つ D6、IA6、D8、および IA8 パラメーターも含まれています。

この経験的分散法は、IOps(3/174-176) を介して他の汎関数で使用できます(S6 は 1.0 である必要があります)。

GD3BJ

Becke-Johnson ダンピングを使用した Grimme の分散の D3 バージョンを追加 [ Grimme11 S. Grimme, S. Ehrlich and L. Goerigk, “Effect of the damping function in dispersion corrected density functional theory,” J. Comp. Chem. 32 (2011) 1456-65. DOI: 10.1002/jcc.21759 ]。以下の表は、GD3BJ パラメーターが定義されている Gaussian 16 の汎関数のリストを示しています。太字で強調表示された汎関数には、示されたキーワードが指定された場合にデフォルトでこの分散モデルが含まれます(例: B2PLYPD3)。それ以外の汎関数では、EmpiricalDispersion=GD3BJ によって分散補正を要求します。

Functional S6 S8 ABJ1 ABJ2
B2PLYPD3 [ Goerigk11 L. Goerigk and S. Grimme, “Efficient and Accurate Double-Hybrid-Meta-GGA Density Functionals—Evaluation with the Extended GMTKN30 Database for General Main Group Thermochemistry, Kinetics, and Noncovalent Interactions,” J. Chem. Theory Comput., 7 (2011) 291-309. DOI: 10.1021/ct10046k ] 0.6400 0.9147 0.3065 5.0570
B97D3 1.0000 2.2609 0.5545 3.2297
PW6B95D3 1.0000 0.7257 0.2076 6.3750
B3LYP 1.0000 1.9889 0.3981 4.4211
BLYP 1.0000 2.6996 0.4298 4.2359
PBE1PBE 1.0000 1.2177 0.4145 4.8593
TPSSTPSS 1.0000 1.9435 0.4535 4.4752
PBEPBE 1.0000 0.7875 0.4289 4.4407
BP86 1.0000 3.2822 0.3946 4.8516
BPBE 1.0000 4.0728 0.4567 4.3908
B3PW91 1.0000 2.8524 0.4312 4.4693
BMK 1.0000 2.0860 0.1940 5.9197
CAM–B3LYP 1.0000 2.0674 0.3708 5.4743
LC-wPBE 1.0000 1.8541 0.3919 5.0897

この経験的分散法は、IOps(3/174-178) を介して他の汎関数で使用できます(S6 は 1.0 である必要があります)。

適用範囲

利用可能な計算

エネルギー、解析勾配、解析周波数、および ADMP 計算で利用できます。

超分極率やラマン強度などの三次の性質は、三次導関数が実装されていない汎関数では利用できません。該当する交換汎関数は G96, P86, PKZB, wPBEh, PBEh、相関汎関数は PKZB、ハイブリッド汎関数は OHSE1PBEOHSE2PBE です。

関連キーワード

関連キーワード

IOp, Int=Grid , Stable , TD , DenFit , B2PLYP , mPW2LYP

実例

実例

DFT 計算のエネルギーは、Hartree-Fock 計算と同様の形式で出力されます。B3LYP 計算から出力されるエネルギーの例を以下に示します。

 SCF Done:  E(RB3LYP) =  -75.3197099428     A.U. after    5 cycles