Force

このキーワードを指定すると，核に働く力（つまり，エネルギーのグラジエント）を計算します。双極子モーメントも（MP2, CC, QCI, CIに対するエネルギーの適切な解析的導関数として）計算されます [202,447]。

EnOnly
エネルギーを数値的に微分して力を計算します。これは解析的グラジェントが利用できない方法に対してはデフォルトです。この手法では数値的不安定性を示すので，最適なステップ幅を各ケースに対して指定しなければならないことに注意してください。

Restart
力の数値的評価を再スタートします。

StepSize=N
数値微分で用いるステップ幅を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)の値で出力されます。