Field キーワードは計算に有限フィールドを追加することを指定します。Gaussian では、場として電気多極子(16 極子まで)または フェルミ接触項を指定できます。フィールドには、次の 2 つの形式のいずれかでパラメーターが必要です。 M±N または F(M)N。ここで M は多極子を指定し、F(M) は原子 M のフェルミ接触摂動を指定します(入力ファイルの分子仕様セクションの順序に従います)。N*0.0001 は、最初の形式では場の大きさを原子単位で、2 番目の形式ではフェルミ接触摂動の大きさを指定します。
たとえば、Field=X+10 は X 方向に 0.001 a.u. の電気双極子場を加えます。Field=XXYZ-20 は、指定した 16 極子場を 0.0020 a.u. の大きさで、デフォルト(標準配向で決まる向き)とは逆向きに加えます。同様に、Field=F(3)27 は原子 3 に対してスピン密度の 0.0027 倍の摂動を与えます。
係数はデカルト演算子行列の係数であることに注意してください。結果を解釈するときは、符号規則の選択に注意する必要があります。
すべてのパラメーターは入力方向です。
フィールド指定パラメーターは、必要に応じて他のオプションの間に置けます。Field を指定すると、アーカイブエントリの生成は無効になります。
オプション
Read
入力ストリームから 34 個の電気多極子コンポーネントの係数を自由形式で読み取ります。
OldRead
入力ストリームから 35 個の電気多極子成分の係数を古い形式の形式 (単極子項を含む) で読み取ります: 形式 3D20.10 を使用します (最初の成分は電荷です)。
RWF
読み取り/書き込みファイルから 35 個の多極子コンポーネントを取得します。
ERWF
読み取り/書き込みファイルから 3 つの電気双極子場の成分のみを抽出します。
Checkpoint
チェックポイントファイルから 35 個の多極子成分を読み取ります。Chk は Checkpoint の同義語です。Geom=Check では Checkpoint がデフォルトです。
EChk
チェックポイント ファイルから 3 つの電気双極子場の成分のみを抽出します。
NoChk
Geom=AllCheck を使用する場合、チェックポイントファイルから外部フィールド係数は取得しません。
制限事項
GVB 計算で対称性を有効にしたまま、指定した場が分子対称性を破る場合、有限場から正しい数値導関数が得られないことがあります。安全のため、GVB と Field を併用するときは常に Guess=NoSymm を使ってください。
一般に、外場によって波動関数の対称性が元の分子と異なる場合、数値導関数は不正確になることがあります。そのため、Field を用いた数値微分では NoSymm を指定するのが適切です。
実例
電場中で構造最適化を行うには、Opt=Z-Matrix NoSymm を用い、入力構造を通常の Z 行列座標または記号付きデカルト座標で与える必要があります。まず Z 行列を使う例を示します。
# RHF/3-21G Field=x+60 Opt=Z-Matrix NoSymm Z-Matrix optimization 0 1 C H 1 B1 H 1 B2 2 A1 H 1 B3 2 A2 3 D1 H 1 B4 2 A3 3 D2 B1 1.070000 B2 1.070000 B3 1.070000 B4 1.070000 A1 109.471203 A2 109.471203 A3 109.471231 D1 120.000015 D2 -119.999993
次に、記号付きデカルト座標を用いる例を示します。
# HF/6-31G(d) Opt=Z-Matrix Field=z-50 NoSymm Symbolic Cartesian coordinates optimization 0 1 O 0 x1 y1 z1 H 0 x2 y2 z2 H 0 x3 y3 z3 x1=0.0 y1=0.0 z1=0.12 x2=0.0 y2=0.75 z2=-0.46 x3=0.0 y3=-0.75 z3=-0.46