このキーワードを指定すると,核に働く力(つまり,エネルギーのグラジエント)を計算します。双極子モーメントも(MP2, CC, QCI, CIに対するエネルギーの適切な解析的導関数として)計算されます [202,447]。
エネルギーを数値的に微分して力を計算します。これは解析的グラジェントが利用できない方法に対してはデフォルトです。この手法では数値的不安定性を示すので,最適なステップ幅を各ケースに対して指定しなければならないことに注意してください。
力の数値的評価を再スタートします。
数値微分で用いるステップ幅を0.0001*Nにします。その単位は,Units=Bohrが指定されていない限り,デフォルトではオングストロームです。デフォルトのステップ幅は0.01 Åです。StepSizeはEnOnlyと組み合わせるときのみ有効です。
解析的グラジェントは全てのSCF波動関数,全DFT法,CIS, MP2, MP3, MP4(SDQ), CID, CISD, CCD, CCSD, QCISD, CASSCF, SAC-CI, 全半経験的方法に対して利用可能です。その他の手法では,力は数値的微分で決定されます。
次の出力では核に働く力が出力されています(この例は水の計算です):
***** AXES RESTORED TO ORIGINAL SET *****
-------------------------------------------------------------------
Center Atomic Forces (Hartrees/Bohr)
Number Number X Y Z
-------------------------------------------------------------------
1 8 -.049849321 .000000000 -.028780519
2 1 .046711997 .000000000 -.023346514
3 1 .003137324 .000000000 .052127033
-------------------------------------------------------------------
MAX .052127033 RMS .031211490
-------------------------------------------------------------------
Internal Coordinate Forces (Hartree/Bohr or radian)
Cent Atom N1 Length/X N2 Alpha/Y N3 Beta/Z J
-------------------------------------------------------------------
1 O
2 H 1 -.023347( 1)
3 H 1 -.023347( 2) 2 -.088273( 3)
-------------------------------------------------------------------
MAX .088272874 RMS .054412682
力は標準配向で決定されますが,出力される前にオリジナル(Z-matrix)の軸に直されます(出力中にそのように記述されます)。つぎに,内部座標を用いているときには,内部座標に関する対応する導関数(長さや角度はZ-matrixのものが使われます)が出力されます。力は,それらの後に各ケースで最大値と2乗平均平方根(RMS)の値で出力されます。
平日9:30~17:30 (土曜日、日曜日、祝祭日、年末年始、夏期休暇は、休日とさせていただきます。)