X

Search Results

Searching....

4.7 Dual-Basis Self-Consistent Field Calculations

4.7.1 Introduction

(July 14, 2022)

The dual-basis approximation 724 Liang W. Z., Head-Gordon M.
J. Phys. Chem. A
(2004), 108, pp. 3206.
Link
, 1138 Steele R. P. et al.
J. Chem. Phys.
(2006), 125, pp. 074108.
Link
, 1141 Steele R. P. et al.
J. Phys. Chem. A
(2006), 110, pp. 13915.
Link
, 288 DiStasio, Jr. R. A., Steele R. P., Head-Gordon M.
Mol. Phys.
(2007), 105, pp. 2731.
Link
, 1140 Steele R. P., Head-Gordon M.
Mol. Phys.
(2007), 105, pp. 2455.
Link
, 1137 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.58)

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 256 Curtiss L. A. et al.
J. Chem. Phys.
(1991), 94, pp. 7221.
Link
, 255 Curtiss L. A. et al.
J. Chem. Phys.
(1998), 109, pp. 7764.
Link
, 254 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.