説明

Gaussian 16 は、さまざまな密度汎関数理論 (DFT) を提供します [ Hohenberg64 P. ホーエンバーグと W. コーン、「不均一電子ガス」 物理学。牧師, 136 (1964) B864-B71。DOI: PhysRev.136.B864 , Kohn65 W. コーンと L. J. シャム、「交換効果と相関効果を含む自己無矛盾方程式」 物理学。牧師, 140 (1965) A1133-A38。DOI: PhysRev.140.A1133 , Parr89 R. G. Parr and W. Yang, Density-functional theory of atoms and molecules (Oxford Univ. Press, Oxford, 1989)., Salahub89 d電子とf電子の挑戦、エド。 D. R. サラブと M. C. ゼルナー (ACS、ワシントン D.C.、1989 年)。DOI: 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 および E. Wimmer、「分子幾何学、振動、および反応エネルギーに対する密度汎関数Gaussian型軌道アプローチ」 J.Chem.物理学。, 96 (1992) 1280-303。DOI: 1.462165 , Becke92 A. D. Becke、「密度関数熱化学。I. 交換のみの勾配補正の効果」 J.Chem.物理学。, 96 (1992) 2155-60。DOI: 1.462066 , Gill92 P.M.W.Gill、B.G.Johnson、J.A.Pople、M.J.Frisch、「さまざまな基底関数を使用した Becke-Lee-Yang-Parr (B-LYP) 密度汎関数理論のパフォーマンス」 化学。物理学。レット。, 197 (1992) 499-505。DOI: 0009-2614(92)85807-M , Perdew92 J. P. Perdew、J. A. Chevary、S. H. Vosko、K. A. Jackson、M. R. Pederson、D. J. Singh、C. Fiolhais、「原子、分子、固体、表面: 交換と相関のための一般化勾配近似の応用」 物理学。 Rev.B, 46 (1992) 6671-87。DOI: PhysRevB.46.6671 , Scuseria92 G. E. Scuseria、「ハートリー・フォックと密度汎関数理論のハイブリッドによる結合クラスター結果の比較」 J.Chem.物理学。, 97 (1992) 7528-30。DOI: 1.463977 , Becke92a A. D. Becke、「密度関数熱化学。II. Perdew-Wang 一般化勾配相関補正の効果」 J.Chem.物理学。, 97 (1992) 9173-77。DOI: 1.463343 , Perdew92a J. P. Perdew および Y. Wang、「電子ガス相関エネルギーの正確かつシンプルな分析表現」 物理学。 Rev.B, 45 (1992) 13244-49。DOI: PhysRevB.45.13244 , Perdew93a J. P. Perdew、J. A. Chevary、S. H. Vosko、K. A. Jackson、M. R. Pederson、D. J. Singh、C. Fiolhais、「正誤表: 原子、分子、固体、表面 - 交換と相関のための一般化勾配近似の応用」 物理学。 Rev.B, 48 (1993) 4978.DOI: PhysRevB.48.4978.2 , Sosa93a C. Sosa およ​​び C. Lee、「非局所補正を使用した遷移構造の密度関数的記述: 水素分子へのシリレン挿入反応」 J.Chem.物理学。, 98 (1993) 8004-11。DOI: 1.464554 , Stephens94 P. J. Stephens、F. J. Devlin、M. J. Frisch、および C. F. Chabalowski、「密度関数力場を使用した振動吸収および円二色性スペクトルの非経験的計算」 J.Phys.化学。, 98 (1994) 11623-27。DOI: j100096a001 , Stephens94a P. J. Stephens、F. J. Devlin、C. S. Ashvar、C. F. Chabalowski、および M. J. Frisch、「振動円二色性スペクトルの理論的計算」 ファラデーについて議論します。, 99 (1994) 103-19。DOI: FD9949900103 , Ricca95 A. リッカとC. W. バウシュリッヒャー ジュニア、「歴代のH」2Fe(H の O 結合エネルギー2O)N+」 J.Phys.化学。, 99 (1995) 9003-07。DOI: j100022a010 DFT の手法とアプリケーションについては、上記の文献を参照してください。エネルギー [ Pople92 J. A. Pople、P. M. W. Gill、B. G. Johnson、「有限基底関数系内のコーン・シャム密度汎関数理論」 化学。物理学。レット。, 199 (1992) 557-60。DOI: 0009-2614(92)85009-Y ]、解析勾配、および真の解析周波数 [ Johnson93a B. G. Johnson および M. J. Frisch、「勾配補正された密度関数エネルギーの解析的二次導関数: 求積重み導関数の効果」 化学。物理学。レット。, 216 (1993) 133-40。DOI: 0009-2614(93)E1238-C , Johnson94 B. G. Johnson および M. J. Frisch、「勾配補正された密度関数エネルギーの解析二次導関数の実装」 J.Chem.物理学。, 100 (1994) 7429-42。DOI: 1.466887 , Stratmann97 R. E. Stratmann、J. C. Burant、G. E. Scuseria、および M. J. Frisch、「密度汎関数理論における調和振動周波数計算の改善」 J.Chem.物理学。, 106 (1997) 10175-83。DOI: 1.474047 これらはすべて、DFT モデルとして利用できます。

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

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

次の式で与えられる同じ最適メモリ サイズ freqmem DFT 周波数計算には推奨されます。

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

Note: 二重ハイブリッド汎関数については、 MP2 計算コストが似ているため、キーワードを使用します。

Accuracy Considerations

DFT 計算では、Hartree-Fock 計算の各主要フェーズに追加のステップが追加されます。このステップは、関数 (または関数のさまざまな導関数) の数値積分です。したがって、Hartree-Fock 計算における数値誤差の原因 (積分精度、SCF 収束、CPHF 収束) に加えて、DFT 計算の精度は数値積分で使用される点の数にも依存します。

UltraFine 統合グリッド ( Integral=UltraFine) は Gaussian 16 のデフォルトです。このグリッドは、合理的な追加コストで計算精度を大幅に向上させます。本番環境の DFT 計算ではこれより小さいグリッドを使用することはお勧めしません。を使用することが重要であることにも注意してください。 same エネルギーを比較するすべての計算用のグリッド (エネルギーの差や生成熱などの計算)。

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

基礎的概念

Background

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

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

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

V The nuclear repulsion energy.
P The density matrix.
⟨hP⟩ The one-electron (kinetic plus potential) energy.
1/2⟨PJ(P)⟩ The classical coulomb repulsion of the electrons.
-1/2⟨PK(P)⟩ The exchange energy resulting from the quantum (fermion) nature of electrons.

密度汎関数理論のコーン・シャム定式化では [ Kohn65 W. コーンと L. J. シャム、「交換効果と相関効果を含む自己無矛盾方程式」 物理学。牧師, 140 (1965) A1133-A38。DOI: PhysRev.140.A1133 ] では、単一の行列式の正確な交換 (HF) は、より一般的な式である交換相関汎関数に置き換えられます。これには、交換エネルギーと電子相関エネルギーの両方を説明する項が含まれますが、後者はハートリー-フォック理論には存在しません。

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

ここで EX[P] は交換関数、E は交換関数です。C[P] は相関関数です。

Kohn-Sham の定式化内では、Hartree-Fock 理論は密度汎関数理論の特殊なケースとみなすことができます。X[P] は交換積分 -1/2 と E で与えられます。C=0。密度汎関数理論で通常使用される汎関数は、密度および場合によっては密度勾配の関数の積分です。

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

ここで、メソッドは E にどの関数 f が使用されるかが異なります。X および (存在する場合) どの f が E に使用されるかC。純粋な DFT 法に加えて、Gaussian は、交換汎関数がハートリー-フォック交換と上記形式の関数積分との線形結合であるハイブリッド法をサポートします。提案された汎関数は、閉じた形式では評価できず、数値求積法によって解決される積分につながります。

Hybrid関数キーワード

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

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

Becke Three-Parameter Hybrid Functionals

これらの汎関数は、1993 年に Becke によって考案された形式を持っています [ Becke93a A. D. Becke、「密度汎関数熱化学。III. 正確な交換の役割」 J.Chem.物理学。, 98 (1993) 5648-52。DOI: 1.464913 ]:

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

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

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

B3LYP LYP 式によって提供される非局所相関と、局所相関 (関数 V ではない) に VWN 関数 III を使用します。 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. コーエンと N. C. ハンディ、「動的相関」 モル。物理学。, 99 (2001) 607-15。DOI: 00268970010023435 ].

Functionals Including Dispersion

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

Long-Range-Corrected Functionals

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

さらに、プレフィックス LC- 平尾と同僚の長い修正を適用するために、ほとんどの純粋汎関数に追加することができます。 Iikura01 飯倉宏、常田隆、柳井隆、平尾和久、「一般化勾配近似交換汎関数の長距離補正スキーム」 J.Chem.物理学。, 115 (2001) 3540-44。DOI: 1.1383587 ]: 例: LC-BLYP.

Other Hybrid Functionals

Functionals from the Truhlar Group

Functionals Employing PBE Correlation

Becke One-Parameter Hybrid Functionals

B1B95 キーワードは元の論文で定義されている Becke の 1 パラメーター ハイブリッド汎関数を指定するために使用されます [ Becke96 A. D. Becke、「密度関数熱化学。IV. 新しい動的相関関数と正確交換混合への影響」 J.Chem.物理学。, 104 (1996) 1040-46。DOI: 1.470829 ].

このプログラムは、Adamo と Barone によって実装された他の同様の 1 パラメーター ハイブリッド関数も提供します。 Adamo97 C. Adamo および V. Barone、「調整可能なパラメータを必要としない信頼性の高い断熱接続モデルに向けて」 化学。物理学。レット。, 274 (1997) 242-50。DOI: S0009-2614(97)00651-9 ]。あるバリエーションでは、 B1LYP、LYP 相関関数が使用されます (上記の B3LYP で説明したように)。別のバージョンでは、 mPW1PW91、Adamo と Barone によって修正された Perdew-Wang 交換を PW91 相関と組み合わせて使用​​します [ Adamo98 C. Adamo および V. Barone、「調整可能なパラメータを必要としない、改善された長距離動作と断熱接続方法を備えた交換機能: mPW および mPW1PW モデル」 J.Chem.物理学。, 108 (1998) 664-75。DOI: 1.475428 ];の mPW1LYP, mPW1PBEmPW3PBE バリエーションも用意されています。

Revisions to B97

Functionals with τ-Dependent Gradient-Corrected Correlation

Older Functionals

Half-and-Half Functionals

以下の機能は、下位互換性のためのみに含まれています。これらは次のとおりであることに注意してください。 not Becke が提案した「ハーフアンドハーフ」汎関数と同じ [ Becke93 A. D. Becke、「ハートリー・フォック理論と局所密度汎関数理論の新しい混合」 J.Chem.物理学。, 98 (1993) 1372-77。DOI: 1.464304 ].

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

User-Defined Hybrid Models

Gaussian 16 は、 一般的な形式:

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

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

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

  • IOp(3/76=mmmmmnnnnn) sets P1 to mmmmm/10000 と P2 to nnnnn/10000. P1 is usually set to either 1.0 または 0.0, depending on whether an exchange functional is desired または not, と any scaling is accomplished using P3 と P4.
  • IOp(3/77=mmmmmnnnnn) sets P3 to mmmmm/10000 と P4 to nnnnn/10000.
  • IOp(3/78=mmmmmnnnnn) sets P5 to mmmmm/10000 と P6 to nnnnn/10000.

例えば、 IOp(3/76=1000005000) P を設定します1 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 の値は P です。2、ScaDFX に指定された値のシーケンスは P4、P3、P6、P5.

Pure関数キーワード

キーワード: 純粋関数

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

Exchange Functionals

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

Correlation Functionals

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

相関関数の変動。 次の相関関数は、さまざまな相関関数のローカル項と非ローカル項を組み合わせたものです。

  • VP86: VWN5 local と P86 non-local correlation functional.
  • V5LYP: VWN5 local と LYP non-local correlation functional.

Standalone Pure Functionals

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

経験的分散力

Dispersion

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

PFD

APFD 関数から Petersson-Frisch 分散モデルを追加 [ Austin12 A. Austin、G. Petersson、M. J. Frisch、F. J. Dobek、G. Scalmani、および K. Throssell、「球面原子分散項を使用した密度汎関数」 J.Chem.理論と計算. 8 (2012) 4989.DOI: ct300778e ].

GD2

Grimme の分散の D2 バージョンを追加 [ Grimme06 S. Grimme、「長距離分散補正を使用して構築された半経験的 GGA タイプの密度汎関数」 J.コンプ。化学。, 27 (2006) 1787-99。DOI: 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

オリジナルの D3 ダンピング機能を備えた Grimme の分散の D3 バージョンを追加 [ Grimme10 S. Grimme、J. Antony、S. Ehrlich、H. Krieg、「94 要素 H-Pu の密度汎関数分散補正 (DFT-D) の一貫した正確な非経験パラメータ化」 J.Chem.物理学., 132 (2010) 154104.DOI: 1.3382344 ]。以下の表は、GD3 パラメーターが定義されている Gaussian 16 の関数のリストを示しています。これらの汎関数の場合、分散は次のように要求されます。 EmpiricalDispersion=GD3.

Functional S6 SR6 S8
B2PLYPD3 [ Goerigk11 L. Goerigk と S. Grimme, “Efficient と Accurate Double-Hybrid-Meta-GGA Density Functionals—Evaluation with the Extended GMTKN30 Database for General Main Group Thermochemistry, Kinetics, と Noncovalent Interactions,” J. Chem. Theory Comput., 7 (2011) 291-309. DOI: 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、L. Goerigk、「分散補正密度汎関数理論における減衰関数の効果」 J. Comp.化学. 32 (2011) 1456-65。DOI: jcc.21759 ]。以下の表は、GD3BJ パラメーターが定義されている Gaussian 16 の関数のリストを示しています。太字で強調表示された関数には、示されたキーワードが指定された場合にデフォルトでこの分散モデルが含まれます (例: B2PLYPD3)。残りの関数については、分散が要求されます。 EmpiricalDispersion=GD3BJ.

Functional S6 S8 ABJ1 ABJ2
B2PLYPD3 [ Goerigk11 L. Goerigk と S. Grimme, “Efficient と Accurate Double-Hybrid-Meta-GGA Density Functionals—Evaluation with the Extended GMTKN30 Database for General Main Group Thermochemistry, Kinetics, と Noncovalent Interactions,” J. Chem. Theory Comput., 7 (2011) 291-309. DOI: 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 である必要があります)。

適用範囲

Availability

エネルギー、分析勾配、および分析周波数。 ADMP calculations.

超分極率やラマン強度などの 3 次の特性は、3 次導関数が実装されていない関数 (交換関数) では利用できません。 G96, P86, PKZB, wPBEhPBEh;相関関数 PKZB;ハイブリッド汎関数 OHSE1PBEOHSE2PBE.

関連キーワード

関連キーワード

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

実例

実例

エネルギーは、Hartree-Fock 計算と同様の形式で DFT 計算で報告されます。から出力されるエネルギーは次のとおりです。 B3LYP calculation:

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