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 個の多極子成分を読み取ります。ChkCheckpoint の同義語です。Geom=Check では Checkpoint がデフォルトです。

EChk

チェックポイント ファイルから 3 つの電気双極子場の成分のみを抽出します。

NoChk

Geom=AllCheck を使用する場合、チェックポイントファイルから外部フィールド係数は取得しません。

適用範囲

利用可能性

単一点エネルギー、構造最適化、周波数、ForceScan 計算で利用できます。

制限事項

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