X

Search Results

Searching....

4.7 Dual-Basis Self-Consistent Field Calculations

4.7.1 Introduction

(April 13, 2024)

The dual-basis approximation 746 Liang W. Z., Head-Gordon M.
J. Phys. Chem. A
(2004), 108, pp. 3206.
Link
, 1170 Steele R. P. et al.
J. Chem. Phys.
(2006), 125, pp. 074108.
Link
, 1173 Steele R. P. et al.
J. Phys. Chem. A
(2006), 110, pp. 13915.
Link
, 298 DiStasio, Jr. R. A., Steele R. P., Head-Gordon M.
Mol. Phys.
(2007), 105, pp. 2731.
Link
, 1172 Steele R. P., Head-Gordon M.
Mol. Phys.
(2007), 105, pp. 2455.
Link
, 1169 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 265 Curtiss L. A. et al.
J. Chem. Phys.
(1991), 94, pp. 7221.
Link
, 264 Curtiss L. A. et al.
J. Chem. Phys.
(1998), 109, pp. 7764.
Link
, 263 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.