Gaussian Keyword Program Development | コンフレックス株式会社

このセクションでは、Gaussian 16 の拡張やインターフェースを行う開発者に役立つキーワードとオプションについて説明します。また、非標準のルートと標準の方向の決定についても説明します。

ここで説明するキーワードとオプションは、新しいメソッドの開発やその他のデバッグ目的には役立ちますが、運用レベルの計算には推奨されません。

一般ジョブ再開

Restart

ここでは、デバッグ向けに設計された Restart の一般的な使い方を説明します。実運用での使用方法は Restart の説明を参照してください。このキーワードは読み取り/書き込みファイルを再利用して計算を再開します。Restart=L1 形式では、読み取り/書き込みファイルを再利用しつつ新しいルートを生成します。

元のルートを使用した再起動では、次の構文を使用して、特定のリンクの発生と、オーバーレイファイルとリンク揮発性ファイルをクリーンアップするか保持するかを指定できます。

 #P Restart [Ln[(m)]] [Clean|KeepOverlay|KeepAll]

すべてのパラメーターを指定すると、ジョブはリンク n の m 回目の出現位置から再開されます。Clean は、Link1 によってすべてのルーチンおよびオーバーレイの揮発性ファイルを削除します。KeepOverlay はオーバーレイ揮発性ファイルのみ保持し、リンク揮発性ファイルは保持しません。KeepAll はすべて保持します。既定値は、リンク内再起動用に読み取り/書き込みファイルが設定されている場合は KeepAll、それ以外は Clean です。

IOp 設定キーワード

IOp1= keyword

このキーワードはオペレーティングシステムインターフェイスのさまざまな詳細を制御します。オプションは標準ですが、すべてのバージョンですべてが実装されているわけではありません (または関連している場合もあります!)。

FileIODump	Dump FileIO tables at the end of each link. FDump is a synonym for this キーワード。
TimeStamp	Turn on time stamping. TStamp is a synonym for this キーワード。
FileIOPrint	Turn on additional debug print in FileIO.
Synch	Currently a no-op (appears in a few test jobs).
NoDFTJ	Turn off use of the pure Coulomb term for non-hybrid DFT (seldom useful in production jobs).
AbelianOnly	Force the use of only abelian symmetry (seldom useful in production jobs).
NoPackSort	Turn off packing addresses into 32-bits during sorts, even if the address space being sorted is < 2 ³¹ (seldom useful in production jobs).

IOp2

このオプションは、動的に割り当てられるメモリの最大量を設定します。 MDV と Core の同義語です IOp2 .

IOp33

これにより、標準のデバッグ印刷オプションが指定どおりに設定されます。たとえば、次のように設定されます。 IOp(33) オーバーレイ 2 のすべての呼び出しで 3 に、そして IOp(33) オーバーレイ 7 のすべての呼び出しで 1 に変わります。

 IOp33(2=3,7=1)

Gaussian 16 IOps リファレンスまた、すべての内部オプション (IOps) についても説明します。

CPHF

デバッグには次のオプションが使用されます。

KeepMicro	Keep all EE centers in CPHF, even for Opt=CalcFC または Opt=CalcAll with non-quadratic microiterations, where atoms that are not used in internal coordinates need not be included in the CPHF.
NoReuse	Do not reuse the electric field CPHF solution in the 2 ^nd (nuclear) CPHF during frequency calculations. The default is ReUse .
XY	Treat real と imaginary perturbations together. The opposite is NoXY , which does them separately. The default is to treat them separately if nuclear perturbations are also being done, but to treat them together if there are only electromagnetic perturbations.
ZVector	Use the Z-Vector method [ Diercksen81 G. H. F. Diercksen, B. O. Roos, and A. J. Sadlej, “Legitimate calculation of 1st-order molecular-properties in the case of limited CI functions – dipole-moments,” Chem. Phys., 59 (1981) 29-39. DOI: , Diercksen81a G. H. F. Diercksen and A. J. Sadlej, “Perturbation-theory of the electron correlation-effects for atomic and molecular-properties – 2nd-order and 3rd-order correlation corrections to molecular dipole-moments and polarizabilities,” J. Chem. Phys., 75 (1981) 1253-66. DOI: , Handy84 N. C. Handy and H. F. Schaefer III, “On the evaluation of analytic energy derivatives for correlated wave-functions,” J. Chem. Phys., 81 (1984) 5031-33. DOI: ] for post-SCF gradients. Allowed と the default if Hartree-Fock 2nd derivatives are not also requested. The NoZVector keyword says to use the full 3 × N _Atoms CPHF for post-SCF gradients.

FMM

デバッグには次のオプションが使用できます。

LMax= N	Specifies the maximum order multipole. The default is 25.
Levels= N	Specifies the number of levels to use in the FMM. The default is 8 for molecules と is adjusted dynamically for PBC.
Tolerance= N	Specifies the accuracy level as 10 ^-N . The default values for N are 11 except for pass 0 of the SCF where it is 7.
JBoxLen= N	Sets the minimum box length (size) to N /1000 Bohrs when doing J. By default, N is 2.5. The maximum of JBoxLen と KBoxLen is used if J と K are done at the same time. BoxLen is a synonym for JBoxLen .
KBoxLen= N	Sets the minimum box length (size) to N /1000 Bohrs when doing K. By default, N is 0.75. The maximum of KBoxLen と JBoxLen is used if K と J are done at the same time.
AllNearField	Turn on all near-field in FMM.
NoParallelCPHF	Forbid parallel execution in FMM during the CPHF phase. NoParCPHF is a synonym for this option.

Integral

デバッグには次のオプションが使用されます。

CNDO	Do calculation in main code using CNDO/2 integrals.
INDO	Do calculation in main code using INDO/2 integrals.
ZIndo1	Do calculation in main code using ZIndo/1 integrals.
ZIndoS	Do calculation in main code using ZIndo/S integrals.
DPRISM	Use the PRISM algorithm [ Gill94 P. M. W. Gill, in Advances in Quantum Chemistry, Vol. 25 (Academic Press, San Diego, CA, 1994) 141-205. DOI: ] for spdf integral derivatives. This is the default.
Rys1E	Evaluate one-electron integrals using the Rys method [ Dupuis76 M. Dupuis, J. Rys, and H. F. King, “Evaluation of molecular integrals over Gaussian basis functions,” J. Chem. Phys., 65 (1976) 111-16. DOI: , King76 H. F. King and M. Dupuis, “Numerical Integration Using Rys Polynomials,” J. Comp. Phys., 21 (1976) 144-65. DOI: , Rys83 J. Rys, M. Dupuis, and H. F. King, “Computation of electron repulsion integrals using the Rys quadrature method,” J. Comp. Chem., 4 (1983) 154-57. DOI: ], instead of the default method. This is necessary on machines with very limited memory.
Rys2E	If writing two-electron integrals, use Rys method (L314) [ Dupuis76 M. Dupuis, J. Rys, and H. F. King, “Evaluation of molecular integrals over Gaussian basis functions,” J. Chem. Phys., 65 (1976) 111-16. DOI: , King76 H. F. King and M. Dupuis, “Numerical Integration Using Rys Polynomials,” J. Comp. Phys., 21 (1976) 144-65. DOI: , Rys83 J. Rys, M. Dupuis, and H. F. King, “Computation of electron repulsion integrals using the Rys quadrature method,” J. Comp. Chem., 4 (1983) 154-57. DOI: , Schlegel84 H. B. Schlegel, J. S. Binkley, and J. A. Pople, “First and Second Derivatives of Two Electron Integrals over Cartesian Gaussians using Rys Polynomials,” J. Chem. Phys., 80 (1984) 1976-81. DOI: ]. This is slower than the default method, but may be needed for small memory machines と is chosen by default if regular (non-Raffenetti) integrals are requested (by the NoRaff option).
DSRys	Use scalar Rys integral derivative code. Can combine with Berny for df only using Rys.
Berny	Use Berny sp integral derivative と second derivative code (L702).
Pass	Pass specifies that the integrals be stored in memory via disk, と NoPass disables this. Synonymous with SCF=[No]Pass , which is the recommended usage.
NoJEngine	Forbid use of special Coulomb code.
NoSP	Do not use the special sp integral program (L311) when writing integrals to disk.
RevDagSam	Reverse choice of diagonal sampling in Prism.
NoSchwartz	Do not use Schwartz integral estimates (only use the heuristic set). Schwartz says to use the Schwartz integral estimates in addition to the heuristic set. The default is to use both.
RevRepFock	Reverse choice of Scat20 vs. replicated Fock matrices.
NoDFTCut	Turn off extra DFT cutoffs.
SplitSP	Split AO S=P shells into separate S と P shells. NoSplitSP is the default.
SplitSPDF	Split AO S=P=D と S=P=D=F shells into S=P, D, と F. NoSplitSPDF is the default.
SplitDBFSPDF	Split density S=P=D と S=P=D=F into S=P, D, と F. NoSplitDBFSPDF is the default.
NoGather	Forbid use of gather/scatter digestion, even when processing small numbers of density matrices. Splatter is a synonym for this option.
ForceNuc	Do nuclear-electron Coulomb with electron-electron.
SepJK	Do J と K in HF/hybrid DFT separately for testing.
Seq2E	Set up for parallel 2 electron integral evaluation but then do not run in parallel (for debugging).
SeqXC	Set up for parallel 2 electron integral evaluation but then do not run in parallel (for debugging).
SeqLinda	Cause Linda workers to run sequentially. Currently just makes the Linda workers other than the master run simultaneously but before the master.
BigAtoms	Make all atom sizes large in XC quadrature.
BigShells	Make all shell sizes large in XC quadrature.
NoSymAtGrid	Do not use (Abelian) symmetry to reduce grid points on symmetry-unique atoms.
LinMIO	Convert to linear storage in FoFCou for testing.
RevDistanceMatrix	Reverse choice of whether to precompute distance matrix during numerical quadrature. The default is to precompute for molecules but not for PBC.
NoDynParallel	Turn off dynamic work allocation.

Sparse

デバッグには次のオプションが使用されます。

Loose	Sets the cutoff to 5 * 10 ^-5 .
Medium	Sets the cutoff to 5 * 10 ^-7 . This is the default for semi-empirical methods.
Tight	Sets the cutoff to 1 * 10 ^-10 . This is the default for DFT methods.
N	Sets the cutoff to 1 * 10 ^-N .

リンク呼び出しと順序の変更

ExtraLinks

これは、追加のリンクの実行を要求します。これらは、通常のリンクの後のオーバーレイのすべてのインスタンスに追加されます。たとえば、ExtraLinks=L9997 を指定すると、オーバーレイ 99 の各インスタンスにリンク 9999 (デフォルト) と 9997 がこの順序で含まれます。

ExtraOverlays

このコマンドは、追加のオーバーレイカードを非標準ルート形式で読み取り、標準ルートの最終（オーバーレイ 99）カードの直前に挿入するように要求します。

Skip

ルート内の最初のオーバーレイカードをスキップします。 Skip=OvNNN 最初にオーバーレイ NNN が出現するまでスキップします。 Skip=M 最初の M カードをスキップします。

Use=L nnn

これにより、プログラム内の代替ルートが指定されます。次のオプションが利用可能です。

L123	Use L123 instead of L115 for IRC. This is the default for IRC, except for IRCMax jobs.
L402	Use old link 402 code for semi-empirical.
L503	Use link 503 for SCF.
L506	Use link 506 for ROHF.

Gaussian の標準方向

計算を実行する前に、分子の対称性を利用して、分子を標準配向と呼ばれる別の座標系に再配向することができます。構造の最適化では、すべてのステップで再配向が発生します。次にプログラムは、最適化中に分子の標準方向が 180 度反転したかどうかをチェックし、反転を回避します。これにより、GaussView で最適化や IRC などをアニメーション化する際のジャンプが回避され、SCFの収束が向上します。

このセクションでは、目標、考慮すべき要素、および分子の標準的な配向の軸を配置するためのさまざまなルールについて説明します。

選択目標

標準指向の規則を選択する際の目標は次のとおりです。

分子を再配向して 3×3 変換行列を簡潔にすること。
内部座標値は異なっていても整数情報が同一の 2 つの Z 行列（幾何最適化の連続点など）は、同じ標準配向を与えること。
同一分子を表す 2 つの異なる Z 行列は、原子番号付けの違いを除いて同じ座標を与えること。
対称性により 0 となる分子軌道係数の数を最大化すること。

一般的な考慮事項

標準的な向きについて考慮する必要がある要素は次のとおりです。

計算全体で右手系座標系を使用します。
分子の電荷中心が原点に来るように並進します。
入力時の順序に対して原子の並び替えは行いません。
デカルト軸の優先順位は X < Y < Z の順とみなします。

軸配置の規則

軸の回転と位置合わせの基準を以下に示します。これらの基準のいずれかを満たすために回転が必要な場合は、次のように定義される X、Y、または Z 軸を中心とした 180 度の回転である必要があります。

X		Rotate about Y
Y		Rotate about Z
Z		Rotate about X

回転軸または電荷主軸は、明確な結果が得られるまで次のテストを連続的に適用することに応じて、平行または逆平行の 2 つの方法のいずれかでデカルト軸と位置合わせできます。

軸上にある最大原子番号の原子群の座標和は正でなければなりません。
電荷の 3 次モーメントは正でなければなりません。
基準軸への原子座標射影の総和は正でなければなりません。
基準軸への射影が非ゼロとなる最初の原子は、その軸方向の射影が正でなければなりません。

電荷主軸配置の規則

他にルールがない場合、最大の主充電モーメントに対応する主軸は、利用可能な最も優先度の高いデカルト軸と一致する必要があります。個々の点グループには、次のような特別な考慮事項があります。

C _s		The molecular plane must be made coincident with the XY plane. Note that although this convention conflicts with Mulliken's suggestion, it is consistent with the character tables of Cotton と Herzberg. The molecule is then rotated about the Z axis according to the rules given below for C _n molecules.
C _2v		The molecular plane is placed in the YZ plane, following Mulliken’s recommendation for planar C _2v molecules. The following tests are successively applied for non-planar molecules: (1) The mirror plane with the most atoms is put in the YZ plane; (2) The mirror plane with the most non-hydrogen atoms is put in the YZ plane; (3) The mirror plane with the lowest numbered atom is made coincident with YZ. Finally, the axes of charge rules are applied (as described above).
Planar , D _2h		Following Mulliken's recommendation, the molecular plane is placed in the YZ plane. The molecule is rotated about the X axis so that the Z axis can pass through either the greater number of atoms, or, if this is not decisive, the greater number of bonds.
C _n		Follow the rules for general symmetric top molecules.
C _i		Translate but do not reorient.
C ₁		Translate but do not reorient.

対称コマ分子の特別規則

対称的な上部分子は、3 つの慣性モーメントのうち 2 つが等しいことで区別されます。したがって、3 次モーメントを基準軸として一意に特定することができ、原子の円集合を考慮することによって点グループが分析されます。

対称トップ分子には次のルールが適用されます。

一意軸は Z 軸に一致させます。
A circular-set of atoms is composed of atoms lying in a plane which have the same atomic number, と are equidistant from a reference axis perpendicular to the plane. Atoms on the reference axis are not included in any circular-set. A circular-set of atoms is generated by a proper rotation axis.
The key atom in a symmetric top molecule is the atom with the lowest number in the key circular-set. The following tests are carried out successively to find the key circular-set:
- どの集合が XY 平面に最も近いか。
- どの集合が Z 軸方向に正の射影を持つか。
- どの集合が Z 軸に最も近いか。
- どの集合が最小原子番号の原子から構成されるか。
orientation is then chosen for the specific point group:
- C _n , C _nh , S _n : The molecule is rotated about the Z axis to maximize the number of pairs of heavy atoms parallel to the Y axis. If no such arrangement is satisfactory, then the key atom is placed in the YZ plane to give it a positive Y coordinate.
- D _n , D _nh : One of the C ₂ axes is made coincident with the Y Cartesian axis. The tests described below are used to decide which C ₂ axis is so positioned.
- D _nd , C _nv : One of the vertical planes is made coincident with the YZ Cartesian plane. The tests below are used to decide which plane is so positioned.
The following are tests for selecting among axes for the D _n , D _nh , D _nd , と C _nv molecules:
- キー原子の Y 軸方向射影を最大化します。
- 2 つの配向で Y 軸方向射影が同じ最大値の場合、X 軸方向射影が最大の方を選択します。
- 分子が XY 平面内にある場合、配向指定を完了するために標準軸配向規則（上記）を X 軸に適用します。

球対称コマ分子の特別規則

球状上部分子は、等しい慣性モーメントを持つことで区別され、球状の原子セットを識別することで特徴付けることができます。

球状の原子セットは、原点から等距離にあり、同じ原子番号を持つ原子で構成されます。球状セットは、原点からの距離の増加と、任意の 1 つの距離での原子番号の増加に基づいて順序付けする必要があります。キー原子は、最初の球状セット内の最も小さい番号の原子です。

一般的にはそうではありませんが、適切な幾何学的制約があれば、D を持たせることは可能です。 _2d、D _2h、またはD ₂ 対称な頂点を持つ分子。このような分子は、上で与えられた規則に従って X、Y、Z 軸と一致する 3 つの垂直な 2 重軸を持っています。

非標準ルートの指定

オプションやリンクの組み合わせが必要な場合、それは通常のオプションとは大幅に異なります。標準ルートにアクセスすると、オーバーレイと関連オプションを含むリンクの完全なシーケンスを読み込むことができます。ジョブタイプ入力セクションは次の行で始まります。

    # 
    
    NonStd

これに、実行順序で目的のオーバーレイごとに 1 行が続き、オーバーレイ番号、スラッシュ、目的のオプション、もう 1 つのスラッシュ、オプションのリストが示されます。実行するリンク、最後にセミコロン:

Ov/Opt=val,Opt=val,…/リンク,リンク,…;

例えば：

7/5=3,7=4/2,3,16;

各リンクでオプション 5 を 3 に、オプション 7 を 4 に設定して、リンク 702、703、および 716 (この順序で) を経由する実行を指定します。すべてのオプションにデフォルト値がある場合、行は次のようになります。

7//2,3,16;

ルート仕様のさらなる特徴は、ジャンプ番号。これは、リンクリストの最後、セミコロンの直前に括弧内に示されています。現在のオーバーレイの完了後にどのオーバーレイ行が実行されるかを示します。省略した場合のデフォルト値は、 +0、プログラムがリスト内の次の行に進むことを示します (行をスキップしません)。ジャンプ番号がに設定されている場合 -4 一方、次のように

7//2,3,16(-4);

その後、オーバーレイで指定されたルートの 4 行前まで実行が続行されます (現在の行はカウントしません)。

この機能により、ルートにループを組み込むことができ、最適化の実行に役立ちます。プログラムチェーンルーチンへの引数はジャンプをオーバーライドできます。これは、構造の最適化中に、最適化が完了するまで一連のオーバーレイラインをループするために使用され、最適化が完了した時点で、ループの終わりに続く行が実行されます。

非標準ルートは通常、最初から作成されるのではなく、目的のものに最も類似した標準ルートによって生成されたシーケンスを出力および変更することによって構築されることに注意してください。これは、次の方法で最も簡単に実現できます。 testrt utility.

単純なルートの例。標準ルート:

# RHF/STO-3G

次の非標準ルートが生成されます。

 1/38=1/1;
 2/12=2,17=6,18=5,40=1/2;
 3/6=3,11=1,16=1,25=1,30=1,116=1/1,2,3;
 4//1;
 5/5=2,38=5/2;
 6/7=2,8=2,9=2,10=2,28=1/1;
 99/5=1,9=1/99;

結果として得られるプログラムのシーケンスを以下に示します。

A Simple Route Sequence

プログラム実行の基本的なシーケンスは、他のプログラムの実行シーケンスと同じです。アブ・イニシオただし、リンク 1 (ルートセクションの読み取りと解釈) が実際の計算に先行し、リンク 9999 (チェックポイントファイルへの書き込み) がその後に続く点が異なります。同様に、MP4 シングルポイントには、母集団分析 (リンク 601) とリンク 9999 の前に、積分変換 (リンク 801 および 804) と MP 計算 (リンク 913) が挿入されています。リンク 9999 は、完了時にジョブステップを自動的に終了します。

ループを含むルート。標準ルート:

# RHF/STO-3G Opt

次の標準ルートを生成します。

1	1/18=20,19=15,38=1/1,3;
2	2/9=110,12=2,17=6,18=5,40=1/2
3	3/6=3,11=1,16=1,25=1,30=1,71=1,116=1/1,2,3;
4	4//1;
5	5/5=2,38=5/2;
6	6/7=2,8=2,9=2,10=2,28=1/1;
7	7//1,2,3,16;
8	1/18=20,19=15/3(2);
9	2/9=110/2;
10	99//99;
11	2/9=110/2;
12	3/6=3,11=1,16=1,25=1,30=1,71=1,116=1/1,2,3;
13	4/5=5,16=3/1;
14	5/5=2,38=5/2;
15	7//1,2,3,16;
16	1/18=20,19=15/3(-5);
17	2/9=110/2;
18	6/7=2,8=2,9=2,10=2,19=2,28=1/1;
19	99/9=1/99;

プログラム実行の結果のシーケンスを以下に示します。

ループを含むルート

このルートはいくつかの考慮事項によって複雑になります。

最適化の最初の点は、後続ステップとは分けて扱う必要があります。これは、一度だけ実行すべき処理（初期分子仕様の読み込みや初期軌道の生成など）があるためです。
構造に対するループ処理が必要であり、最適化プログラム（この場合は Berny optimizer, Link 103）が、次の構造計算が必要か、すでに最適化済みかを判定します。
収束済み構造が与えられた場合、プログラムは勾配を1回だけ計算し、構造が最適化済みであることを認識して終了します。
集団解析と軌道出力は、既定では最初と最後の点でのみ実行され、中間構造（通常は情報量が少ない）では実行されません。

最初の点は、ルートに積分、推定、SCF、および積分導関数の 2 つの基本シーケンスを含めることによって処理されています。最初のシーケンスには、リンク 101 (初期構造を読み取るため)、リンク 103 (独自の初期化を行う) が含まれており、リンク 401 に初期推定を生成するように指示するオプションが設定されています。 2 番目のシーケンスは、最適化の過程でリンク 103 で生成された構造を使用し、次の初期推定として前の構造から波動関数を取得するようにリンク 401 に指示するオプションを設定しています。

8 行目の前方ジャンプには、リンク 103 が (特別なアクションを実行せずに) 正常に終了した場合、後続の行 (リンク 202 および 9999 を呼び出す) がスキップされるという効果があります。通常、リンク 103 の 2 回目の呼び出しでは、初期勾配が検査され、新しい構造が選択されます。次に実行されるリンクはリンク 202 で、新しい構造を処理し、その後にメインの最適化ループを構成する 2 番目のエネルギー + 勾配シーケンスの残りが続きます。リンク 103 の 2 回目の呼び出しで構造が収束していることが判明すると、リンク 103 はジャンプを抑制するフラグを付けて終了し、次の行によってリンク 202、601、および 9999 が呼び出され、ジョブが完了します。

行 11 ～ 16 はメインの最適化ループを形成します。これにより、最適化の 2 番目以降のポイントの積分、波動関数、および勾配が評価されます。これはリンク 103 で終了します。構造がまだ収束していない場合、リンク 103 は新しい構造を選択して通常どおり終了し、16 行目で後方ジャンプが実行され、次の行が 11 行目に処理されて新しいサイクルが始まります。リンク 103 は、構造が収束したことを検出すると、終了してジャンプを抑制し、最終行 (17 ～ 19) が処理されるようにします。

リンク 601 の最後のインスタンスは、要求に応じて、最終的な多極子モーメント、軌道および個体数解析を出力します。最後に、リンク 9999 はアーカイブエントリを生成し、ジョブステップを終了します。

MP および CI の最適化には、変換オーバーレイと相関オーバーレイ (8 と 9)、およびオーバーレイ 7 の前に挿入されたポスト SCF 勾配オーバーレイ (11 と 10) があります。同じ 2 相ルート構造が数値微分に使用され、周波数または分極率が生成されます。

のルート Opt=Restart 基本的には、最初のステップの特別な行が省略された、元の最適化のメインループにすぎません。リンク 103 の 2 回目の呼び出しは保持され、実際の再起動が行われます。

読み書き可能なファイル番号

以下はそのリストです読み書き可能なファイル。チェックポイントファイルに永続的に存在するものには、文字が付いています。 P、チェックポイントファイル上に一時的に存在するものには、文字がマークされます。 T . T ファイルは、最適化または数値頻度実行の再開に使用するために保存されますが、ジョブステップが正常に完了すると削除されます。

Type	RWF	説明
P	501	Gen array.
P	502	/LABEL/—Title と atomic orbital labels.
	503	Connectivity information (MxBond,0),NBond(NAtoms),IBond(MxBond,NAtoms),RBond(MxBond,NAtoms), where arrays are rounded to a multiple of IntPWP.
	504	Dipole derivative matrices (NTT,3,NAt3).
P	505	Array of copies of /Gen/ from potential surface scan.
P	506	Saved basis set information before massage, uncontraction, etc.
P	507	ZMAT/ と /ZSUBST/.
P	508	/IBF/ Integral Bugger Format.
	509	Incomplete integral buffer.
T	510	/FPINFO/ Fletcher-Powell optimization program data.
P	511	/GRDNT/ energy, First と second derivatives over variables, NVAR.
P	512	Pseudo-potential information.
P	513	/DIBF/ integral derivative buffer format.
	514	Overlap matrix, optionally followed by absolute overlap と absolute overlap over primitives.
	515	Core-Hamiltonian. There are four matrices here: H(α), the α core Hamiltonian; H(β), the β core Hamiltonian; G'(α), the α G' contribution to Fock matrix; G'(β), the β G' contribution to Fock matrix. H(α) と H(β) differ only if Fermi contact integrals have been added. The G' matrices are for perturbations which are really quadratic in the density (and hence have a factor of 1/2 in their contribution to the energy as compared to the true one-electron terms) but which are computed externally to the SCF.
	516	Kinetic energy と modifications to the α と β core Hamiltonian. These include ECP terms, Douglas-Kroll-Hess corrections, multipole perturbations と Fermi contact perturbations. The latter are used for calculations in which the nuclear と electronic Coulomb terms are computed together, such as the Harris functional と PBC calculations. For semi-empirical, holds the core Hamiltonian without nuclear attraction terms for use in the initial guess.
	517	Fermi contact integrals.
	518	Multipole integrals, in the order X,Y,Z,XX,YY,ZZ,XY,XZ,YZ,XXX,YYY,ZZZ,XYY,XXY,XXZ,XZZ,YZZ,YYZ,XYZ,XXXX,YYYY,ZZZZ, XXXY,XXXZ,YYYX,YYYZ,ZZZX,ZZZY,XXYY,XXZZ,YYZZ,XXYZ,YYXZ,ZZXY.
T	519	Common /OptEn/—optimization control for link 109.
T	520	Electronic state: count と packed string (1+9 integers).
P	521	Electronic state: count と packed string (1+9 integers).
P	522	Eigenvalues, alpha と if necessary, beta.
	523	Symmetry assignments.
P	524	MO coefficients, real alpha.
P	525	(no longer used)
P	526	MO coefficients, real beta.
P	527	(no longer used)
T	528	SCF density matrix, real alpha.
T	529	(no longer used)
T	530	SCF density matrix, real beta.
T	531	(no longer used)
T	532	SCF density matrix, real total.
T	533	(no longer used)
T	534	SCF density matrix, real spin.
	535	(no longer used)
	536	Fock matrix, real alpha.
	537	Fock matrix, imaginary alpha.
	538	Fock matrix, real beta.
	539	Fock matrix, imaginary beta.
	540	Molecular alpha-beta overlap (U), real.
	541	Molecular alpha-beta overlap (U), imaginary.
T	542	Pseudo-potential information.
T	543	Pseudo-potential information.
T	544	Pseudo-potential information.
P	545	/ORB/ – window information.
	546	Bucket entry points.
	547	Eigenvalues (double precision with window: always alpha と beta, even in RHF case).
P	548	MO coefficients (double precision with window, alpha と if necessary beta). Complex if necessary.
	549	Molecular orbital alpha-beta overlap, double precision with window.
T	550	Potential surface scan common block.
T	551	Symmetry operaiton info (permutations, transformation matrices, etc.)
P	552	Character strings containing the stoichiometric formula と framework group designation.
T	553	Temporary storage of common/gen/ during FP optimizations.
T	554	Alternate starting MO coefficients, from L918 to L503, real alpha. Also MO coefficients in S ^-1/2 basis for L509 と rotation angles from L914 to L508.
	555	Alternate starting MO coefficients, from L918 to L503, imaginary alpha.
T	556	Alternate starting MO coefficients, from L918 to L503, real beta. Also MO coefficients in S ^-1/2 basis for L509 と rotation angles from L914 to L508.
	557	Alternate starting MO coefficients, from L918 to L503, imaginary beta.
	558	Saved HF 2nd derivative information for G1, G2, etc.
	559	Common /MAP/.
	560	Core-Hamiltonian (a. o. basis) with 2 j – k part of deleted orbitals added in. (i.e. frozen core).
P	561	External point charges または SCIPCM informations.
P	562	Symmetry operations と character table in full point group.
T	563	Integer symmetry assignments (α).
T	564	Integer symmetry assignments (β).
T	565	Lists of symmetry equivqlent shells と basis functions.
T	566	Unused in G16.
T	567	GVB pair information (currently dimensioned for 100 paired orbitals).
P	568	Saved hamiltonian information from L504 と L506.
P	569	Saved read-in window.
P	570	Saved amplitudes (IAS1,IAS2,IAD1,IAD2,IAD3; only IAS1 と IAD2 for closed-shell).
	571	Energy weighted density matrix.
	572	Dipole-velocity integrals <Phi\|Del\|Phi'>, X, Y, と Z, followed by R × Del integrals (R × X, R × Y, R × Z).
	573	More SCIPCM information.
T	574	/MSINFO/ Murtaugh-Sargent program data.
T	575	/OPTGRD/ Gradient optimization program data for L103, L115, と L509.
T	576	/TESTS/ Control constants in L105.
T	577	Symmetry adapted basis function data.
T	578	A logical vector indicating which MO's are occupied.
T	579	NEQATM (NATOMS*NOP2) for symmetry.
T	580	NEQBAS (NBASISNOP2+NBas6DNOp2) for symmetry.
T	581	NSABF (NBASIS*NOP2) for symmetry. Followed by matching integer character table, always (8,8).
T	582	MAPROT (3*NBASIS) for symmetry.
T	583	MAPPER (NATOMS) for symmetry.
P	584	FXYZ (3*NATOMS) cartesian forces. During PSCF gradient runs, there will be two arrays here: first the PSCF gradient, then the HF only component (needed for PSCF with HF 2nd deriv).
P	585	FFXYZ (NAT3TT) cartesian force constants (lower triangle).
T	586	Info for L106, L110, と L111.
T	587	L107 (LST) data.
	588	Sx over cartesians in the ao basis.
	589	Hx over cartesians in the ao basis.
	590	F(x) over cartesians in the ao basis (all α, followed by all β for UHF) (without CPHF terms).
	591	U1(A,I) — MO coefficient derivatives with respect to electric field と nuclear coordinates.
	592	Electric field と nuclear P1 (AO basis).
	593	Electric field と nuclear W1 (AO basis).
	594	Electric field と nuclear S1 (MO basis).
	595	Magnetic field U1(A,I) — Del(X,Y,Z) then R × (X,Y,Z), 6 α followed by 6 β.
	596	Full MO Fock derivatives in the MO basis, including CPHF terms.
P	597	Configuration changes for Guess=Alter.
	598	User Name.
	599	Density basis set info: NDBFn, NVar, U0, DenBfn(4,NDBfn), ITypDB(NDBfn), Var(NVar), IJAnDB(NDBfn), IVar(4,NDBfn).
	600	Saved data for intra-link restart.
P	601	Saved structures, と possibly forces と force constants along reaction path. All structures, then all forces, then all force constants.
	602	Post-SCF two-particle density matrix.
P	603	Density Matrices at various levels of theory.
T	604	common /drt1/ from drt program … misc integer ci stuff, followed by variable dimension drt arrays.
P	605	Atomic charges from Mulliken Populations, ESP fits, etc. Bitmap followed by 0 または more NAtoms arrays. Bits 0/1/2/3/4 Mulliken/ESP-fit/Bader/NPA/APT.
	606	SCF orbital symmetries in Abelian point group. Alpha and, if necessary, beta, full set followed by windowed set.
	607	Window'd orbital symmetries like rw 606 (always alpha と beta).
	608	IBF for sorted integrals (normally on SAO unit).
	609	Bit map for sorted integrals (normally on SAO unit).
	610	Sorted AO integrals (normally on SAO unit).
	611	NTT maps for sorted integrals (normally on SAO unit).
	612	Some 1E generators for direct CI matrix element generation.
	613	Some more 1E generators for direct CI matrix element generation.
	614	Configuration information for CAS-MP2.
	615-616	Used for CAS-MP2.
	617	Spin-orbit integrals.
P	618	Nuclear coordinate third derivatives.
P	619	Electric field derivatives: 1 WP word bit map, dipole, dipole derivative, polarizability, dipole 2nd derivatives, polarizability derivatives, hyperpolarizability.
	620	Magnetic field derivatives for GIAOs.
	621	Susceptiblity と chemical shift tensors.
	622	Partial overlap derivatives (<Mu\|dNu/da>, NBasisNBasisNAt3).
P	623	Born-Oppenheimer 波動関数 derivatives (<Phi\|d2Phi/dadb> for electronic Phi と a,b nuclear, NAt3TT).
	624	Unused in G16.
	625	Expansion vectors と AY products from CPHF, in the order Y α, AY α, Y β, AY β.
	626	MCSCF MO 1PDM (NTT).
	627	MCSCF MO Lagrangian (NTT).
	628	MCSCF MO 2PDM (NTT,NTT) または NVTTTT.
	629	AO 2PDM (shell order).
T	630	MCSCF information.
	631	Post-SCF Lagrangian (TA, then TB if UHF).
	632	OV3NAtoms, followed by OV*NVar d2E/d(V,O)d(XYZ,Atom).
P	633	Excited-state CI densities.
T	634	SCF Restart information (alpha, then possibly beta MOs).
P	635	CIS と CASSCF CI coefficients と restart information.
	636	NBO analysis information.
	637	Natural orbitals generated by link 601.
	640	MCSCF data または CIS AO Tx's for 2nd derivatives.
	641	MCSCF data for 2nd derivatives.
	642	MCSCF data for 2nd derivatives.
	643	MCSCF data for 2nd derivatives.
	644	MCSCF data for 2nd derivatives.
	645	MCSCF data for 2nd derivatives.
	646	MCSCF data for 2nd derivatives.
	647	MCSCF data for 2nd derivatives.
	648	MCSCF data for 2nd derivatives.
	649	Eigenvalue derivatives (non-canonical form even if done canonically).
	650	2PDM derivatives, (LenTQ,NDeriv,ShellQuartet) order.
	651	Full U's, canonical または non-canonical as requested.
	652	Generalized density derivatives for the current method (NTT,NDeriv,IOpCl+1).
	653	Lagrangian derivatives for the current method (NTT,NDeriv,IOpCl+1).
	654	Gx(Gamma).
	655	G(Gamma).
	656	Non-symmetric S1 と S2 parts of Lagrangian for MP2 または CIS second derivatives.
	657	tIx と tIx/D matrices from L811 for L1112.
	658	L(x) from L1111.
	659	MO correlated W for correlated frequencies.
	660	2nd order CPHF results: Pia,xy, Sxy, Fxy (complete) all in MO basis, PSF α then PSF β if UHF.
	661	Computed electric field from L602.
	662	Points for electrostatic evaluation.
T	663	Saved information for L117 と L124.
	664	Spin projection data.
P	665	Redundant coordinate information.
	666	(no longer used)
	667	CIS AO Fock matrix.
	668	CIS Gx(T) matrices.
	669	Saved /ZMat/ と /ZSubst/ during redundant optimzations.
P	670	New format basis set data (compressed /B/).
P	671	New optimization (L103/L104) data.
P	672	Unused in G16.
	673	Global optimization data.
	674	ONIOM internal data.
	675	Saved files for LS during ONIOM.
	676	Saved files for MS during ONIOM.
	677	Saved files for LM during ONIOM.
	678	Saved files for HS during ONIOM.
	679	Saved files for MM during ONIOM.
	680	Saved files for LL during ONIOM.
	681	Saved files for HM during ONIOM.
	682	Saved files for ML during ONIOM.
	683	Saved files for HL during ONIOM.
	684	SABF information for DBFS: equivalent to files 577 と 581 for AOs.
	685	Cholesky U, または transformation to surviving basis functions.
	686	Cholesky U ^-1 .
	687	Molecular mechanics parameters.
	688	Density in orthogonal basis (α spin) for ADMP または sparse SCF.
	691	Saved initial files during ONIOM (gridpoint 17, hence 674+17=691).
	694	Permutation applied to MOs for post-SCF symmetry.
	695	Magnetic properties.
	696	Saved magnetic field density derivatives.
	698	Saved initial structure during geometry optimization, in standard orientation, also used for constraints with the force constants following the structure.
	699	Density in orthogonal basis (β spin) for ADMP または sparse SCF.
	700	Saved /Mol/ for ONIOM.
	701	Saved Trajectory/IRC/Optimization history.
	702	Fit density for Coulomb.
	703	Fit density for Coulomb.
	704	Saved XC contribution to electric field F(xa) for polar derivatives.
P	710	Basic PCM information.
P	711	Other PCM data.
P	712	Non equilibrium data for PCM.
	713	Saved information for RFO with ONIOM microiterations.
	714	Saved model system information for ONIOM microiterations.
	715	Saved rigid fragment information for ONIOM microiterations.
T	716	Saved copy of basis set data for counterpoise.
T	717	Saved copy of ECP data for counterpoise.
T	718	Saved copy of fitting basis for counterpoise.
	719	Saved DiNa information.
P	720	Saved DiNa information.
	721	Frequency-dependent properties.
	722	Derivatives of frequency-dependent properties.
	723	Density fitting matrices (metrics).
	724	Density fitting basis (same format as /B/).
	725	DBF symmetry information (NEqDBF(NDBF,NOp2),NEqDB6(NDBF6D,NOp2)).
	726	DBF shell symmetry information (NEqDBS(NDBShl,NOpAll)).
	727	F(x)(P-Pfit) for density fitting second derivatives.
	728	PBC cell replication information.
	729	Alternate new guess during optimizations.
	730	Counterpoise input specification.
	731	Counterpoise intermediate data.
	732	Basis set for finite nuclei.
	733	PBC Cell scalars と integer cell indices.
	734	State-specific input parameters for SAC-CI.
	735	Excitation lables of SAC と SAC-CI.
	736	Eigenvalues と eigenvectors of SAC と SAC-CI.
	737	H matrices と their indices of non-zero elements used for SAC/SAC-CI.
	738	Saved atomic parameters for DFTB/EHTSC.
	739	Temporary storage for imaginary core Hamiltonian perturbations.
	740	Orbital information for SAC/SAC-CI gradients と PES by GSUM.
	741	MOD Orbital information for SAC gradients.
	742	Saved quadrature grid.
	743	Alpha Fock matrices in orthonormal basis for ADMP, also alpha HF Fock matrix for non-HF post-SCF.
	744	Beta Fock matrices in orthonormal basis for ADMP, also beta HF Fock matrix for non-HF post-SCF.
	745	K-integration mesh information.
	746	Eigenvalues と orbitals at all k-points.
	747	Information for external low-level calculations for ONIOM.
	748	TS vector information for ONIOM TS optimizations.
	749	Conical intersection information for ONIOM.
	750	Not used in G16.
	751	Temporary storage for SO ECP integrals.
	752	Pseudo-canonical MO Fock matrix for ROMP と ROCC.
	753	Data for FD polar derivatives.
	754	Saved PCM charge derivatives.
	755	PCM inverse matrices.
	756	Charge information for ONIOM.
	757	MO:MO embedding charge data for L924.
	758	Derivatives of embedding charges, when computed explicitly.
	759	Basis set info for density embedding.
P	760	Full set of pseudocanonical orbitals for RO.
	761	Charges from external PCM iterations (both L117 と L124).
	762	Saved weights for non-symmetric Mulliken analysis.
	763	File for FC/HT integrals.
	764	File for FC/HT integrals.
P	765	Saved normal modes.
	766	Saved QuadMac vectors (temporary).
	767	CIS coefficients reordered by symmetry.
	768	Semi-empirical parameters.
	769	Saved MOs during numerical differentiation.
P	770	Saved ground-to-excited state energies と transition moments.
	771	EOM iteration information.
	772	Symmetry operations と character table in Abelian point group.
	989	Multi-step job information (1000 reals と 2000 integers).
	990	KJob info in some implementations.
	991	Holds file names, ID's と save flags.
	992	Used for link substitution information in some implementations.
	993	COMMON /INFO/
	994	COMMON /PHYCON/
	995	COMMON /MUNIT/
	996	COMMON /IOP/
P	997	COMMON /MOL/
P	998	COMMON /ILSW/
	999	Overlay data.