Searching....

# 4.7.3 Dual-Basis Dynamics

(December 20, 2021)

The ability to compute SCF and MP2 energies and forces at reduced cost makes dual-basis calculations attractive for ab initio molecular dynamics simulations, which are described in Section 9.9. Dual-basis BOMD has demonstrated 1073 Steele R. P., Head-Gordon M., Tully J. C.
J. Phys. Chem. A
(2010), 114, pp. 11853.
savings of 58%, even relative to state-of-the-art, Fock-extrapolated BOMD. Savings are further increased to 71% for dual-basis RI-MP2 dynamics. Notably, these timings outperform estimates of extended Lagrangian (“Car-Parrinello”) dynamics, without detrimental energy conservation artifacts that are sometimes observed in the latter. 454 Herbert J. M., Head-Gordon M.
J. Chem. Phys.
(2004), 121, pp. 11542.

Two algorithm improvements make modest but worthwhile improvements to dual-basis dynamics. First, the iterative, small-basis calculation can benefit from Fock matrix extrapolation. 454 Herbert J. M., Head-Gordon M.
J. Chem. Phys.
(2004), 121, pp. 11542.
Second, extrapolation of the response equations (“$Z$-vector” equations) for nuclear forces further increases efficiency. 1076 Steele R. P., Tully J. C.
Chem. Phys. Lett.
(2010), 500, pp. 167.
(See Section 9.9.) Q-Chem automatically adjusts to extrapolate in the proper basis set when DUAL_BASIS_ENERGY is activated.