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

プログラム：

NAME

detci - 行列式構成相互作用プログラム

DESCRIPTION

プログラム デッチ で非相対論的電子ハミルトニアン演算子を対角化します。
スレーター決定基の基礎。 使用される決定要因のセット (CI 空間) は、
さまざまな方法。 プログラムは、
制限付きアクティブ スペース CI。 これには、フル CI までの CISD、CISDT、CISDTQ などが含まれます。
参照がすべての決定要因として選択される複数参照 CI と同様に、
いくつかのMO活性空間で最大n個の電子が励起されます。 これには CISD[T] が含まれます。
CISD[TQ]、および二次 CI (SOCI)。

参考文献

制限付きアクティブ スペース CI:

1. 完全および制限付きの決定基ベースの構成相互作用アルゴリズム
Configuration Interaction Spaces、J. Olsen、BO Roos、P. Jorgensen、および HJ Aa。
ジェンセン、Ｊ．Ｃｈｅｍ． 物理学。 89、2185 (1988)。

2. 完全構成相互作用 (FCI) 計算で XNUMX 億の限界を超える
Ｊ．Ｏｌｓｅｎ、Ｐ．Ｊｏｒｇｅｎｓｅｎ、およびＪ．Ｓｉｍｏｎｓ、Ｃｈｅｍ． 物理。 レット。 169, 463 (1990)。

三次仮想部分空間 (RAS IV):

1. 制限された XNUMX 重および XNUMX 重を含むコンパクトな変分波動関数
Ｅｘｃｉｔａｔｉｏｎｓ、ＣＤ ＳｈｅｒｒｉｌｌおよびＨＦ Ｓｃｈａｅｆｅｒ、Ｊ．Ｐｈｙｓ． 化学。 100、6069-6075
とします。

DETCI プログラム:

1. CD Sherrill、大規模な完全参照およびマルチ参照のための計算アルゴリズム
構成相互作用波動関数、博士論文、ジョージア大学、アテネ、
ジョージア州、1996 年。

ファイル REQUIRED

input.dat - 入力ファイル
file71 - 変換された一電子積分
file72 - 変換された XNUMX 電子積分

一時的 ファイル 中古

file50 - ハミルトニアンの対角線
file51 - CI ベクトル
file52 - シグマ ベクトル
file53 - D ファイル (修正ベクトル)

ファイル UPDATED

Output.dat - 出力ファイル

入力 FORMAT

次のコマンドライン引数が使用可能です。

-静かな これは次と同じ結果になります 印刷=0.

-o fname
出力ファイルのファイル名を指定します。 デフォルトはoutput.datです。

-e このオプションにより、合計 CI エネルギーまたはエネルギーが、という名前のファイルに書き込まれます。
detci_energys.dat.

-c 値
DETCAS 計算に役立つ、CI ベクトルのより緩やかな収束を提供します。 の
値は実数であり、次のような整数ではありません 収束. 使用される収束
の緩い方になります 値 と 収束.

このプログラムの追加入力はファイルから読み取られます。次のキーワードは
有効:

収束 = 整数
CI ベクトルに必要な収束。 の RMS が次の場合に収束が達成されます。
CI ベクトルの誤差は 10**(-n) 未満です。 デフォルトは、エネルギーの場合は 4、エネルギーの場合は 7 です。
グラデーション用。 これは、に見られるものと同じ CI ベクトル収束基準ではありません。
グガチ。

文書 = 整数配列
このベクトルは、各 irrep 内の二重占有軌道の数を示します。 がある
デフォルトはありません。

SOCC = 整数配列
このベクトルは、各 irrep 内の単一占有軌道の数を示します。 がある
デフォルトはありません。

DIAG_METHOD = string
これは、ハミルトニアンの対角化に使用する方法を指定します。 の
有効なオプションは次のとおりです。 RSP、H行列全体を形成し、libciomrを使用して対角化します
すべての固有値を取得するには (nb には巨大なメモリが必要です)。 オルセン、オルセンの
前処理付き逆部分空間法 (1990); ミトルシェンコフ、2x2 を使用するには
オルセン/デビッドソン法; と DAVIDSON （または SEM) リューの同時拡張を使用する
根が XNUMX つしか見つからない場合、Davidson 法と同じです。
SEMデバッグモードもあり、 セムテストを選択します。 SEM 方法が一番
堅牢ですが、ディスク上に 2(N*M)+1 個の CI ベクトルが必要です。ここで、N は最大です。
反復回数、M は根の数です。

前提条件者 = string
これは、選択した対角化で使用する前提条件のタイプを指定します
方法。 有効なオプションは次のとおりです。 DAVIDSON ハミルトニアン行列を近似する
対角要素によって; H0BLOCK_INV の正確なハミルトニアンを使用します
H0_BLOCKSIZE を明示的に反転します。 GEN_DAVIDSON これはスペクトルを行います
H0BLOCK の分解; ITER_INV 反復アプローチを使用して
H0BLOCK の補正ベクトル。 H0BLOCK_INV、GEN_DAVIDSON、および ITER_INV
アプローチはすべて形式的に同等ですが、ITER_INV は計算量が少なくなります。
高い。 デフォルトは DAVIDSON.

REFERENCE = string
参照関数のタイプを指定します。 これは RHF または ROHF です。 UHFと
TWOCON はサポートされていません。 ROHF の場合、多重度 1 は開殻を意味します。
シングレット。 プログラムは開殻一重項に対して実行されますが、まだ実行されていません
この場合、正しい XNUMX つの行列式参照を使用するように適切に調整されているため、
完全な CI を除いて、オープンシェルのシングレット参照で実行することはお勧めしません。

UPDATE = string
DAVIDSON 標準の DAVIDSON 更新または修正ベクトル式を採用し、
オルセン OLSEN 補正ベクトルを使用します。 デフォルトは DAVIDSON.

HD_OTF = boolean
TRUE の場合、ハミルトニアン行列の対角要素がオンザフライで計算されます。
それ以外の場合は、対角要素ベクトルがディスク上の別のファイルに書き込まれます。 デフォルト
は本当です。

HD_AVE = string
HD_EXACT 正確な対角エネルギーを使用して、展開ベクトルを生成します。
スピン対称性を破る。 HD_KAVE スピン結合の対角エネルギーを平均化します
スピン純粋展開ベクトルを生成するように設定します。 ORB_ENER 軌道の和を採用
スピン純粋な展開ベクトルを与えるエネルギー近似ですが、通常は
ダビッドソンの反復回数。 伝道者 軌道の和と差を使う
エネルギーを SCF 参照エネルギーと組み合わせて、スピンの純粋な膨張ベクトルを生成します。
ライニンガー から一電子の寄与を差し引く近似
軌道エネルギーを 0.5 倍し、XNUMX 電子の寄与を足し戻します。
で、スピンの純粋な拡張ベクトルを生成し、真にあなたによって開発され、次のように機能します
だけでなく 伝道者.

ノドファイル = boolean
NUM_ROOTS = 1 の場合のみ可能です。BVEC ファイルの最後のベクトル空間を使用して
別の DVEC ファイルを使用するのではなく、スクラッチ DVEC を書き込みます。

エネルギー_コンバージェンス = 整数
CI エネルギーに必要な収束。 単一点エネルギーのデフォルトは 6 です。
グラデーションまたは CASSCF の場合は 8 です。

EX_LVL = 整数
仮想軌道への励起の励起レベル (デフォルト 2、つまり CISD)。

VAL_EX_LVL = 整数
RAS II の軌道のリファレンスの励起レベル。 デフォルトはゼロです。

FCI = boolean
このフラグが TRUEの場合、文字列の格納は Full 用に単純化されます
CI と計算に必要なオーバーヘッドが少なくなります。 ただし、最終結果は
ときのものと同一であること FCI = 間違った情報. 次の場合、予期しない結果が生じる可能性があります FCI =
TRUE 焙煎が極度に未発達や過発達のコーヒーにて、クロロゲン酸の味わいへの影響は強くなり、金属を思わせる味わいと乾いたマウスフィールを感じさせます。 EX_LVL フル CI と一貫性がありません。

FROZEN_DOCC = 整数配列
各既約範囲内の最低エネルギー二重占有軌道の数
そこからの興奮がない表現。 の綿の注文
不可解な表現が使用されています。 デフォルトはゼロ ベクトルです。

FROZEN_UOCC = 整数ベクトル
各既約表現における最高エネルギーの空軌道の数
そこには興奮はありません。 デフォルトはゼロ ベクトルです。

フリーズコア = boolean
このオプションは、凍結したコア軌道を含めるかどうかを決定します
暗黙的に (true) または明示的に (false)。 前者の場合、エネルギッシュな
凍結されたコア軌道からの寄与は、XNUMX電子に折りたたまれます
変換プログラムによって計算された「凍結コアエネルギー」に変換されます。
デフォルトはtrueです。

EXPORT_VECTOR = boolean
これは、実行の最後に収束したベクトルを保存するかどうかを指定します。 の
ベクトルは、他のプログラムが使用できるように透過的な形式で保存されます
簡単に。 フォーマットは src/lib/libqt/slatrdset.h。 デフォルトは
偽。

NUM_EXPORT = 整数
If EXPORT_VECTOR が true に設定されている場合、これによりベクトルの数が決まります。
実行の最後にエクスポートする必要があります。 デフォルトは 1 です。

GUESS_VECTOR = string
これは、CI 反復で使用する推定ベクトルのタイプを指定します。 現在
SEM 反復法でのみ使用されます。 受け入れられる値は UNIT 単位ベクトル
推測してみて （NUM_ROOTS と NUM_INIT_VECS 両方とも 1 でなければなりません); H0_BLOCK 固有ベクトルを使用するには
H0 BLOCK サブマトリックスから (デフォルト); Dファイル 使用する NUM_ROOTS 以前に収束した
D ファイルのベクトル。 と MP2 MP2波動関数を推測として使用する（機能しない
現時点では）。

H0_BLOCKSIZE = 整数
このパラメーターは、ハミルトニアンの "H0" ブロックのサイズを指定します。
きっちり解決。 SCF エネルギーが最も低い n 個の行列式が選択され、
ハミルトニアンの部分行列は、これらの行列式を使用して形成されます。 この部分行列は
CI 反復の収束を加速するために使用されます。 オルセン と ミトルシェンコフ
反復スキーム、また、 SEM 場合のメソッド
GUESS_VECTOR = H0_BLOCK. デフォルトは 40 です。プログラムによって変更される場合があることに注意してください。
Ms=0 の場合の与えられたサイズ (Ms0 = TRUE) H0 ブロックに含まれていると判断した場合
時間反転対称性によって関連付けられた行列式のペアの XNUMX つのメンバーのみ。 為に
ブロック サイズが非常に小さい場合、これにより H0 ブロック全体が削除される可能性があります。 の
これが発生した場合、プログラムは警告を出力する必要があります。

H0_BLOCK_COUPLING_SIZE = 整数
一般化された範囲内のカップリング ブロックのサイズを指定するパラメーター
デビッドソンプレコンディショナー. デフォルト値は 1000 です。

MAX_DET = 整数
決定要因の最大数を設定します。 CI スペースがこれよりも大きい場合、
プログラムは中止します。 このオプションは、非常に大きな計算が行われないようにするために存在します。
偶然走る。 現在の開発段階では、デフォルトは 10000 ですが、
やがて引き上げられる。

マキシター = 整数
ハミルトニアンを対角化する最大反復回数。 デフォルトは 12 です。

Ms0 = boolean
If TRUE、状態の Ms=0 コンポーネントを使用します。 デフォルトは TRUE クローズドシェルの場合
〜へ 間違った情報 それ以外は。 関連する S パラメータに一致する最初のデバイスのリモートコントロール URL を返します。

Nプリント = 整数
この値は、出力される決定要素の数を指定します。
それらの係数、最終 CI の最も重要な決定要因のリスト
ベクター。 デフォルトは 20 です。

NUM_ROOTS = 整数
この値は、世俗から取得されるルートの数を示します。
方程式。 デフォルトは XNUMX です。 複数のルートが必要な場合は、設定します DIAG_METHOD
〜へ SEM (または、非常に小さなケースでは、 RSP or セムテスト).

NUM_INIT_VECS = 整数
CI 反復手順で使用する初期ベクトルの数。 デフォルトは
根の数。

OPDM = boolean
If TRUE 一粒子密度行列を計算し、 OPDM_書き込み デフォルトは
TRUE。 のデフォルト値 OPDM is 間違った情報.

OPDM_FILE = 整数
一粒子密度行列を書き込むファイル(ユニット番号) OPDM_書き込み =
TRUE. 現在のデフォルト値は 73 です。

OPDM_書き込み = boolean
XNUMX 粒子密度マトリックスをディスクに書き込むかどうかのフラグ。

OPDM_PRINT = boolean
XNUMX 粒子密度マトリックスを出力するかどうかのフラグ。

OPDM_DIAG = boolean
XNUMX 粒子密度行列を対角化するかどうかのフラグ。

ウォートノス = boolean
CI 自然軌道を PSIF_CHKPT に書き込むかどうかのフラグ。

ORBSファイル = 整数
各種 CI 自然軌道を記述するためのファイル (ユニット番号)。 デフォルト値は
76.

OPDM_AVE = boolean
を取得するために、複数のルートにわたって OPDM を平均化するかどうかのフラグ。
状態平均一粒子密度行列。 この密度行列は対角化できます
CI 自然軌道を取得します。

ORBS_ROOT = 整数
CI 自然軌道を書き込むルート番号を設定するためのフラグ
PSIF_CHKPT。 デフォルト値は 1 (最低ルート) です。

PRINT = 整数
このオプションは、出力の冗長性を決定します。 1 または 2 の値は、
最小限の印刷。値 3 は詳細な印刷を指定します。 値 4 または 5 は、
デバッグに使用されます。 テストケースが非常に小さい場合を除き、レベル 5 を使用しないでください (例:
STO H2O CISD)。

ROOT = 整数
XNUMX 粒子密度行列を書き出すルート (XNUMX 粒子密度
行列はすべてのルートに対して書き込まれます)。 州固有の CASSCF または CI に役立ちます
励起状態での最適化。

S = 整数
スピン量子数 S の値は、このオプションによって与えられます。 デフォルトは 0 です
(一重項)。 これが実際に使用される唯一のことは、
Ms=0 コンポーネントが使用されている場合の CI ベクトルの冗長な半分 (つまり、 Ms0 =
TRUE）。 S が整数でない場合、このパラメーターを入力する必要はありません。
そのような状態は Ms=0 コンポーネントを持つことができないためです。

TPDM = boolean
If TRUE 二粒子密度行列を計算し、 TPDM_書き込み デフォルトは
TRUE。 のデフォルト値 TPDM is 間違った情報.

TPDM_ファイル = 整数
次の場合に XNUMX 粒子密度行列を書き込むためのファイル (ユニット番号) TPDM_書き込み =
TRUE. 現在のデフォルト値は 74 です。

TPDM_書き込み = boolean
XNUMX 粒子密度マトリックスをディスクに書き込むかどうかのフラグ。

TPDM_PRINT = boolean
XNUMX 粒子密度マトリックスを出力するかどうかのフラグ。 通常、非常に
小さなケースのデバッグを除いて、悪い考えです。

また、PSI の初心者が使用する必要のない、あまり一般的ではない入力もあります。
使用しています。

ベンダゾリ = boolean
Bendazzoli らの論文に基づいて、いくつかのルーチンを使用してシグマを計算します。
遅くて価値がないようです。 最終的に消えるかもしれません。 のみ動作します
完全なCIであり、彼らの巧妙なスキームがどのように見えるかを覚えていません
一般に RAS に拡張されます。

CALC_SSQ = boolean
TRUE の場合、最終 CI の S^2 演算子の期待値を計算します
各ルートの波動関数。 原則として、DETCI は S^2 固有関数を生成する必要があります。
デフォルトは FALSE です。

COLLAPSE_SIZE 整数
Davidson 部分空間が折りたたまれたときに保持するベクトルの数を指定します (
マックスネクト 下）。 XNUMX より大きい場合、折りたたまれた部分空間が最良の状態を保持します。
前の n 回の反復に対する CI ベクトルの推定。 デフォルトは 1 です。

FIRST_TMP_UNIT = 整数
DETCI によって使用される最初のスクラッチ ファイルに関連付けられたファイル (ユニット) 番号を示します。
他のスクラッチ ファイルには、この時点から順に番号が付けられます。
H(diag)、C、S、D。これらの論理ファイルはそれぞれ、多数の物理ファイルを占有します。
さらにあいまいな入力パラメーターによって指定された NUM_HD_TMP_UNITS、
NUM_C_TMP_UNITS、 NUM_S_TMP_UNITS、 NUM_D_TMP_UNITS。 ユーザーが指定することもできます
パラメータを使用して、これらのセットごとに異なる開始点
FIRST_HD_TMP_UNIT など。 ファイルを複数のユニットに分割すると役立つ場合があります
に存在する大きなファイルに対処する際の整数のサイズの問題を回避します。
DETCI および PSI I/O ライブラリ内。 しかし、繰り返しますが、私はそれをテストして何を確認していません
起こります。 各セクションの最初のユニットは、FILES の見出しの下に出力されます。
DETCI 実行を開始するパラメーター出力。

FZC = boolean
凍結したコア軌道を真に凍結したものとして扱うかどうかを決定します (つまり、
計算から完全に欠落し、 FZC = TRUE) またはそれらが存在するかどうか
二重占有に制限されています（FZC = 間違った情報）。 GUGA CI プログラムでは、これは
FZC軌道とCOR軌道と呼ばれるものの違い。 一般に、
凍結コア軌道の積分は DETCI では必要ありませんが、必要になる場合があります
MCSCF またはグラデーション用。

アイコア = 整数
CI ベクトルのバッファリングの処理方法を指定します。 値が 0 の場合、プログラムは
一度に 1 つの RAS サブブロックで I/O を実行します。 2 は一度に CI ベクトル全体を使用します。 そしてXNUMX
一度に 0 つのirrep ブロックを使用します。 2 または XNUMX の値は、
I/O (反復で H を構築するときに C ベクトルの複数の読み取りが必要
DIAG_METHOD = SEM の場合はサブスペース)、必要なコア メモリは少なくて済みます。

私は止める = boolean
If TRUE 文字列情報が形成された後、その前に DETCI が停止します。
積分が読み取られます。 ユーザーができるように、最終的に整数に変更される場合があります
複数の停留所から選択します。

マックスネクト = 整数
ディスク上に保持できる Davidson サブスペース ベクトルの最大数を指定します。
CI 係数とシグマ ベクトル。 (XNUMX つの H(diag) ベクトルと数
D ベクトルの数は根の数に等しい)。 ディスク上のベクトルの数
の値に達する マックスネクト、Davidson部分空間はに折りたたまれます
COLLAPSE_SIZE 各ルートのベクトル。 これは、ディスク容量を節約するのに非常に役立ちます。
デフォルトは マキシター * NUM_ROOTS + NUM_INIT_VECS.

混合 = boolean
これは、"混合" RAS II/RAS III 励起が CI に許可されるかどうかを決定します。
スペース。 これは、RAS CI に追加の制約を設定する場合に便利です。

ミックス4 = boolean
これは MIXED キーワードに似ていますが、RAS IV への励起を指します。

ヌニッツ = 整数
C ベクトルの格納に使用されるスクラッチ ファイルの数 (および sigma にも使用)
ベクトル）。

OEI_ERASE = boolean
これは、プログラムがその後に一電子積分ファイルを消去するかどうかを決定します
読んだ。 デフォルトは最終的に true になりますが、開発中は
デフォルトはfalseです。

OEI_FILE = 整数
このキーワードにより、ユーザーは変換された XNUMX 電子積分ファイルを指定できます。
デフォルトは71です。

PRINT_CIBLKS = boolean
プログラムがすべてのブロックの概要を印刷するかどうかを指定します。
CI ベクトル (行列形式にキャストできます。参考文献を参照してください。)

R4S = boolean
RAS IV 文字列を最小限のセットに制限して、メモリを節約します。 あなたがいる場合
このオプションについて懸念がある場合は、David に連絡してアドバイスを求める必要があります。
DETCIの専門家。

REF_SYM = 整数
このオプションを使用すると、ユーザーは、
参照。 これはおそらく完全な CI の場合にのみ意味があり、おそらく意味がない
単位ベクトルの推測で動作します。 番号付けは、完全にゼロから始まります。
対称的な修復。

REPL_OTF = boolean
その場で文字列置換を行うかどうかを DETCI に指示します。 保存できます
膨大な量のメモリ (特に切り捨てられた CI の場合) ですが、やや不安定で、
しばらくテストされていません。 私が思い出したように、それは特定のクラスでのみ機能します
RAS計算。 デビッドに連絡して支援を受けてください。 最終的には、オンザフライ
代わりのものは、時間がかからないように、よりスマートな方法でやり直す必要があります
膨大な CPU 時間。 これらの線に沿った作業は、一度開始され、完了する可能性があります
いずれはね。

再起動 = boolean
このオプションにより、ユーザーは終了した DETCI 反復を再開できます。
時期尚早。 CI およびシグマ ベクトルがディスク上にあることを前提としています。 の数
によって指定されるベクトル RESTART_VECS ルートごとに XNUMX つのベクトルに折りたたまれます。

RESTART_VECS = 整数
If 再起動 = TRUE これは、読み取る CI (およびシグマ) ベクトルの数を指定します。
ディスク。 通常、これは、正常に完了した反復の数です。
以前の実行は、その実行の根の数に掛けます。

TEI_ERASE = boolean
これは、プログラムがその後に XNUMX 電子積分ファイルを消去するかどうかを決定します。
読んだ。 デフォルトは最終的に true になりますが、開発中は
デフォルトはfalseです。

TEI_FILE = 整数
このキーワードにより、ユーザーは変換された XNUMX 電子積分ファイルを指定できます。
デフォルトは72です。

MPN = boolean
このオプションが TRUE の場合、DETCI は MPn シリーズを k 次まで計算します。
maxnvect によって決定されます。 オープンシェルシステム (REF=ROHF、WFN = ZAPTN) の場合、DETCI
ZAPTn シリーズを計算します。 GUESS_VECTOR を TRUE に設定する必要があります。 HD_OTF は
TRUE に設定します。 HD_AVE は orb_ener に設定する必要があります。

保存_MPN2 = 整数
MPN が TRUE で WIGNER が TRUE の場合、このオプションは有効になります。 1 に設定した場合
MP(2n-1) エネルギーが節約されます。 2 に設定すると、MP(2n-2) エネルギーが節約されます。 もしあれば
その他の値 MPn エネルギーが節約されます。 デフォルトは 0 です。

9 2月、1996 デッチ(1)

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