10.13.2 Numerical Calculation of Static Polarizabilities

Where analytic gradients are not available, static polarizabilities (only) can be computed via finite-difference in the applied field, which is known as the finite field (FF) approach. Beginning with Q-Chem 5.1, a sophisticated “Romberg” approach to FF differentiation is available, which includes procedures for assessing the stability of the results with respect to the finite-difference step size. The Romberg approach is described in Section 10.13.3. This section describes Q-Chem’s older approach to FF calculations based on straightforward application of small electric fields along the appropriate Cartesian directions.

Dipole moments can be calculated numerically as the first derivative of the energy with respect to $\overrightarrow{F}$ by setting JOBTYPE = DIPOLE and IDERIV = 0. If IDERIV is not specified explicitly, the dipole moment will be calculated analytically, which for post-Hartree–Fock levels of theory invokes a gradient calculation in order to utilize the relaxed wavefunction.

Similarly, set JOBTYPE = POLARIZABILITY for numerical evaluation of the static polarizability tensor $\overleftrightarrow{\alpha }$. This is performed by either first-order finite difference, taking first-order field derivatives of analytic dipole moments, or by second-order finite difference of the energy. The latter is useful (indeed, required) for methods where analytic gradients are not available, such as CCSD(T) for example. Note, however, that the electron cloud is formally unbound in the presence of static electric fields and therefore a bound solution is a consequence of using a finite basis set. (With analytic derivative techniques the perturbing field is infinitesimal so this is not an issue.) This fact, along with the overall sensitivity of numerical derivatives to the finite-difference step size, means that care must be taken in choosing the strength of the applied finite field.

To control the order for numerical differentiation with respect to the applied electric field, use IDERIV in the same manner as for geometric derivatives, i.e., for polarizabilties use IDERIV = 0 for second-order finite-difference of the energy and IDERIV = 1 for first-order finite difference of gradients. In addition, for numerical polarizabilities at the Hartree-Fock or DFT level set RESPONSE_POLAR = $-1$ in order to disable the analytic polarizability code.

In finite-difference geometric derivatives the $rem variable FDIFF_STEPSIZE controls the size of the nuclear displacements (see Section 9.2.1), but here it controls the magnitude of the electric field perturbations: