X

Search Results

Searching....

4.7 Dual-Basis Self-Consistent Field Calculations

4.7.1 Introduction

(November 19, 2024)

The dual-basis approximation 779 Liang W. Z., Head-Gordon M.
J. Phys. Chem. A
(2004), 108, pp. 3206.
Link
, 1213 Steele R. P. et al.
J. Chem. Phys.
(2006), 125, pp. 074108.
Link
, 1216 Steele R. P. et al.
J. Phys. Chem. A
(2006), 110, pp. 13915.
Link
, 313 DiStasio, Jr. R. A., Steele R. P., Head-Gordon M.
Mol. Phys.
(2007), 105, pp. 2731.
Link
, 1215 Steele R. P., Head-Gordon M.
Mol. Phys.
(2007), 105, pp. 2455.
Link
, 1212 Steele R. P., DiStasio, Jr. R. A., Head-Gordon M.
J. Chem. Theory Comput.
(2009), 5, pp. 1560.
Link
to self-consistent field (HF or DFT) energies provides an efficient means for obtaining large basis set effects at vastly less cost than a full SCF calculation in a large basis set. First, a full SCF calculation is performed in a chosen small basis (specified by BASIS2). Second, a single SCF-like step in the larger, target basis (specified, as usual, by BASIS) is used to perturbatively approximate the large basis energy. This correction amounts to a first-order approximation in the change in density matrix, after the single large-basis step:

Etotal=Esmall basis+tr[(Δ𝐏)𝐅]large basis. (4.59)

Here 𝐅 (in the large basis) is built from the converged (small basis) density matrix. Thus, only a single Fock build is required in the large basis set. Currently, HF and DFT energies (SP) as well as analytic first derivatives (FORCE or OPT) are available.

Note:  As of version 4.0, first derivatives of unrestricted dual-basis DFT energies—though correct—require a code-efficiency fix. We do not recommend use of these derivatives until this improvement has been made.

Across the G3 set 279 Curtiss L. A. et al.
J. Chem. Phys.
(1991), 94, pp. 7221.
Link
, 278 Curtiss L. A. et al.
J. Chem. Phys.
(1998), 109, pp. 7764.
Link
, 277 Curtiss L. A. et al.
J. Chem. Phys.
(2000), 112, pp. 7374.
Link
of 223 molecules, using cc-pVQZ, dual-basis errors for B3LYP are 0.04 kcal/mol (energy) and 0.03 kcal/mol (atomization energy per bond) and are at least an order of magnitude less than using a smaller basis set alone. These errors are obtained at roughly an order of magnitude savings in cost, relative to the full, target-basis calculation.