4.6.3 Linear Scaling Exchange (LinK) Matrix Evaluation

Hartree-Fock calculations and the popular hybrid density functionals such as B3LYP also require two-electron integrals to evaluate the exchange energy associated with a single determinant. There is no useful multipole expansion for the exchange energy, because the bra and ket of the two-electron integral are coupled by the density matrix, which carries the effect of exchange. Fortunately, density matrix elements decay exponentially with distance for systems that have a HOMO/LUMO gap. ¹⁰⁹⁹ Schwegler E., Challacombe M., Head-Gordon M.
J. Chem. Phys.
(1997), 106, pp. 9708. Link The better the insulator, the more localized the electronic structure, and the faster the rate of exponential decay. Therefore, for insulators, there are only a linear number of numerically significant contributions to the exchange energy. With intelligent numerical thresholding, it is possible to rigorously evaluate the exchange matrix in linear scaling effort. For this purpose, Q-Chem contains the linear scaling K (LinK) method ⁹⁰⁶ Ochsenfeld C., White C. A., Head-Gordon M.
J. Chem. Phys.
(1998), 109, pp. 1663. Link to evaluate both exchange energies and their gradients ⁹⁰⁷ Ochsenfeld C.
Chem. Phys. Lett.
(2000), 327, pp. 216. Link in linear scaling effort (provided the density matrix is highly sparse). The LinK method essentially reduces to the conventional direct SCF method for exchange in the small molecule limit (by adding no significant overhead), while yielding large speedups for (very) large systems where the density matrix is indeed highly sparse. For full details, we refer the reader to the original papers. ⁹⁰⁶ Ochsenfeld C., White C. A., Head-Gordon M.
J. Chem. Phys.
(1998), 109, pp. 1663. Link ^{, 907} Ochsenfeld C.
Chem. Phys. Lett.
(2000), 327, pp. 216. Link LinK can be explicitly requested by the following option (although Q-Chem automatically switches it on when the program believes it is the preferable algorithm).