The dual-basis
approximation
724
J. Phys. Chem. A
(2004),
108,
pp. 3206.
Link
,
1138
J. Chem. Phys.
(2006),
125,
pp. 074108.
Link
,
1141
J. Phys. Chem. A
(2006),
110,
pp. 13915.
Link
,
288
Mol. Phys.
(2007),
105,
pp. 2731.
Link
,
1140
Mol. Phys.
(2007),
105,
pp. 2455.
Link
,
1137
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:
(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
J. Chem. Phys.
(1991),
94,
pp. 7221.
Link
,
255
J. Chem. Phys.
(1998),
109,
pp. 7764.
Link
,
254
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.