cscf - クラウドでオンライン

プログラム：

NAME

cscf - ハートリー・フォック方程式を解きます

DESCRIPTION

プログラム cscf 反復手順を実行してハートリー・フォック方程式を解きます。

このプログラムは、D2h 対称性とそのサブグループおよび軌道占有に限定されています。
は整数である必要があります。 したがって、特定の純粋な角運動量状態は、
縮退軌道の部分的な占有は、現在のコードでは取得できません。 のために
たとえば、最低エネルギーから導出された線形ONOの2PIu（ダブレットPI u）状態
線形 (pi u)1 構成は 2B2u (ダブレット B2u) または 2B3u としてのみ計算できます。
3PIu (ダブレット PI u) 状態の (ダブレット B 2u) 成分、およびその結果の空間
波動関数はPI対称性を持ちません。 しかし、それはある意味では望ましいことであり、
エネルギーは曲げ角度の連続関数となるためです。 エネルギーの計算
2B2u (ダブレット B 2u) または 2B3u (ダブレット B 3u) としてベント構成を使用し、純粋な 2PIu を実行します。
線形ジオメトリでの (ダブレット PI u) 状態は、顕著な不連続性をもたらします。

ほとんどの場合、三重項状態は二重に縮退したものの二重占有によって生じる
(e')3 または (e")2 から生じる 2A2 (三重項 A 2) 状態などの軌道
D3h 対称の構成、または (pi g)3 または (pi の 2SIGMAg (三重項 SIGMA g) 状態)
Dinfh (D infinity h) 対称の u)2 構成は、適切な空間対称性を持ちます。
これらの同じ電子配置から生じる一重項状態は、本質的に
複数構成であるため、単一構成では十分に表現されません
波動関数。

参考文献

PKファイル方式:

1. RC ラフェネッティ、Chem. 物理学。 レット。 20（1973）335。

分子対称性と閉殻 HF 計算:

1. M.Dupuis、HFKing、Int. Ｊ．Ｑｕａｎｔ． 化学。 11 (1977) 613.

クローズドシェルの DIIS:

1. P. Pulay、Chem. 物理学。 レット。 73（1980）393。

2. P. プーレー、J. Comp。 化学。 3 (1982) 556。

開殻の結合係数 (アルファとベータ):

1. CCJ Roothaan、修道院長。 物理学。 32（1960）179。

ダンピング：

1. ハートリー博士、「原子構造の計算」（ワイリー：ニューヨーク）1957 年。

2. MC Zerner および M. Hehenberger、Chem. 物理学。 レット。 62 (1979) 550。

レベルシフト:

1. VR サンダースと IH ヒリアー、インターナショナルＪ．Ｑｕａｎｔ． 化学。 7 (1973) 699。

収束中 CSCF

困難なオープンシェルの場合には、適切なクローズドシェルを使用することをお勧めします。
最初に計算を実行し (余分な電子を追加または削除し)、次にこの SCF ベクトルを実行します。
目的の開殻波動関数の推定として使用できます。 TCSCF の場合は、常に次のようになります。
最初に閉じたシェル (またはおそらく適切なトリプレット) SCF を実行してから、使用するのが賢明です。
これは TCSCF の推測です。

オープン シェル システムの場合、0.5 ～ 3.0 のレベル シフト値が推奨されます。 から始めてください
最初の SCF 計算では高い値 (2.0 ～ 3.0) を使用し、その後、それを (0.5 ～ 1.0 に) 下げます。
後続の実行では、収束した SCF ベクトルを開始点として使用します。

このバージョンのコードではサポートが終了していることに注意することが非常に重要です。
オープンタイプ。 1 しなければなりません つかいます 　 NEW キーワード REFERENCE および マルチプ 〜へ 指定する 　 type of SCF
必要です。

入力 FORMAT

XNUMXμmの波長を持つ cscf プログラムは、デフォルトのキーワード パス (最初のキーワード パス) を介して検索します。 SCF その後 DEFAULT)
次のキーワードの場合:

LABEL = string
出力に含める文字列です。 この文字列は使用されません
プログラムによって。 デフォルトはありません。

WFN = string
これは、最終的に望ましいタイプの波動関数です。 デフォルトは SCF.

オープンタイプ is いいえ より長いです 。

REFERENCE = string
これは、実行したい SCF 計算のタイプを指定します。 それは次のいずれかです RHF
(閉殻シングレットの場合)、 ROHF (制限されたオープンシェル計算の場合)、 UHF
(無制限のオープンシェル計算の場合)、 ツーコン (XNUMX つの構成の場合
一重項)、または 特殊。 場合 特殊 が与えられた場合、アルファとベータの結合
係数は次のように指定する必要があります。 ALPHA および BETA キーワード。 デフォルトは RHF.

MULTP= 整数
分子の多重度を指定します。 デフォルトはシングレットです。

料金= 整数
分子の電荷を指定します。 デフォルトは 0 です。

文書 = 整数ベクトル
これにより、各既約範囲内の二重占有軌道の数が得られます。
表現。 デフォルトはありません。 これが指定されていない場合、CSCF は次のことを試みます。
コアハミルトニアンを使用して職業を推測します。

SOCC = 整数ベクトル
これにより、各既約範囲内の単一占有軌道の数が得られます。
表現。 デフォルトはありません。

ダータイプ = string
これは、最終的に実行される導関数の順序を指定します。 使用されています
を通じて、タンピングされたコーヒーベッドの上から均一にフィルターバスケットの内の粉に浸透していきます。 SCF 特定のファイルを書き込むかどうかを決定するプログラムであり、また、
波動関数のデフォルトの収束を決定するために使用されます。 デフォルトは
最初.

マキシター = 整数
これにより、最大反復回数が得られます。 デフォルトは 40 です。

収束 = 整数
これは、波動関数がどの程度厳密に収束するかを指定します。 収束とは、
密度行列の RMS 変化 (「デルタ P」) を、
与えられた値。 収束基準は 10**(-整数）。 次の場合、デフォルトは 7 です。
両言語で ダータイプ = NONEを および WFN = SCF が与えられ、それ以外の場合は 10 になります。

レベルシフト = リアル
レベルシフトを指定します。 デフォルトは 1 です。

直接 = boolean
SCF 計算を積分直接法で行うかどうかを指定します。 の
デフォルトはfalseです。

PRINT_MOS = boolean
分子軌道を印刷するかどうかを指定します。 デフォルトは false です。

あまり使用されない入力パラメータも多数あります。 もししないなら
次のオプションの意味を理解してから、それらのオプションが
入力。 ほとんどの場合、デフォルトのままで機能します。 これらは指定されています
次のキーワードで:

DELETE_INTS = boolean
次の場合、Integrals ファイルは消去されます。 WFN = SCF および ダータイプ = 最初 or ダータイプ = NONEを.
積分ファイルを保持したい場合は、次のように設定します DELETE_INTS = 偽。 デフォルトは
真。

再注文 = string
このパラメーターは、分子軌道の並べ替えを制御します。 に設定されている場合 BEFORE そうして
チェックポイント ファイルからの推定軌道が並べ替えられます。 に設定されている場合 AFTER、収束した
軌道はチェックポイント ファイルに書き込まれる前に並べ替えられます。 どちらでも
場合 モーオーダー 並べ替えマップを指定するには、パラメーターを指定する必要があります。 デフォルトは
軌道を並べ替えないこと。

モーオーダー = 整数ベクトル
これは分子軌道並べ替えベクトルを指定します。 以下の場合にのみ使用されます。
再注文 が設定されています。 このベクトルは、すべての軌道をその新しいインデックスにマッピングします。 モーオーダー = (0
2 1) 順序を変更した後、軌道 1 と 2 が交換されることを指定します。 のランク
このベクトルは MO の数と同じです。 指数はピッツァー順です
(対称性によって並べられ、次に各対称ブロック内のエネルギーによって並べられます)、基数 0。 CSCF
指定された MOORDER が異なる irrep からの軌道を混合する場合、失敗する可能性があります。 そこには
はデフォルトではありません。

ALPHA = 実数ベクトル
If オープンタイプ = 特殊の場合、このパラメータはアルファ結合係数を与えます。
このベクトルの要素の数は MM(MM+1)/2 です。MM は要素の数です。
単一占有分子軌道を含む既約表現。 そこには
はデフォルトではありません。

BETA = 実数ベクトル
If オープンタイプ = 特殊の場合、このパラメータはベータ結合係数を与えます。
このベクトルの要素の数は MM(MM+1)/2 です。MM は要素の数です。
単一占有分子軌道を含む既約表現。 そこには
はデフォルトではありません。

ゲス = string
このオプションは、CSCF が使用する固有ベクトルの初期推定のタイプを決定します。
現時点で有効なオプションは次のとおりです: (1) ゲス = CORE、それが使用されるようになります
計算を開始するためのコア ハミルトニアン固有ベクトル。 (2) ゲス = AUTO which
その結果、チェックポイント ファイル内の MO ベクトルを使用しようとするか、または次のような手段に頼ることになります。
そのファイルに固有ベクトルがない場合、コアは推測します。 デフォルトの場合 AUTO.

イプリント = 整数
これは印刷オプションです。 デフォルトは 0 です。

MO_OUT = boolean
対称性と職業を持つ軌道を最後に出力します。
計算。 デフォルトは true です。

ROTATE = boolean
これが false の場合、分子軌道は回転しません。 回転のみ
開殻系の仮想軌道に影響を与えます。 このパラメータは true である必要があります
相関勾配の場合は false であり、XNUMX 次導関数以降は false でなければなりません。
次の場合、デフォルトは false です。 WFN = SCF それ以外の場合は真です。

CHECK_ROT = boolean
上で説明した分子軌道回転をチェックして、列が存在しないことを確認します。
SCF 固有ベクトル行列は回転によって交換されます。 以下の場合は効果がありません ROTATE =
false。 デフォルトは true です。

CHECK_MO_ORTHOGONALITY = boolean
分子軌道が正規直交しているかどうかを確認します。 デバッグの場合にのみ役立ちます。 の
デフォルトはfalseです。

DIIS = boolean
これにより、dis が使用されるかどうかが決まります。 デフォルトは true です。

DISSTART = 整数
これにより、DIIS が使用される最初の反復が得られます。 デフォルトは 0 です。

NDIIS = 整数
これにより、diis プロシージャで使用する誤差行列の数が得られます。 デフォルト
クローズド シェルの場合は 6、オープン シェルの場合は 4、tcscf の場合は 3 です。

DISDAMP = リアル
これにより、dis 手順の減衰係数が得られます。 デフォルトは 0.0 です。
クローズド シェル、オープン シェルの場合は 0.02、tcscf の場合は 0.01。

増分 = リアル
これは、ci 係数が再計算される頻度を決定するために tcscf で使用されます。
小さな数値 (~0.25) を指定すると、ほぼ scf ごとに再計算されます。
反復。 デフォルトは 0.25 です。

DYN_ACC = boolean
直接 scf を実行する場合、動的積分精度カットオフを行うかどうかを指定します。
使用されます。 デフォルトは true (動的カットオフを使用) です。 最初の反復は次のとおりです
XNUMX 桁までの精度の積分で実行されます。 密度が収束した後
10^-5 または 30 回の反復が完了し、完全な積分精度が使用されます。 scfの場合
これを設定して動的カットオフを無効にすると、収束の問題が発生します。
変数を false にすると役立つかもしれません。

オルソグ_オンリー = boolean
CASSCF または軌道最適化のための他の非 HF/KS スキームでは、
現在のジオメトリに対して他のジオメトリの MO を再直交化するのに役立ちます。
これらは、新しい MO の初期推定として使用できます。 これは次の方法で実行できます。
CSCF を実行する オルソグ_オンリー = true。 軌道が直交化されると、
プログラムは SCF 計算を実行せずに終了します。 このキーワードは、
チェックポイント ファイルに以前の軌道がない場合は無視されます。 デフォルトは true
if WFN = デトカス.

30 5月、1991 cscf(1)

onworks.net サービスを使用してオンラインで cscf を使用する