- Search
- Download PDF

(May 16, 2021)

While atomic orbital standard basis sets are used to expand one-electron functions such
as molecular orbitals, *auxiliary basis sets* are also used in many Q-Chem jobs to
efficiently approximate products of one-electron functions, such as arise in electron
correlation methods.

For a molecule of fixed size, increasing the number of basis functions
*per atom*, $n$, leads to $\mathcal{O}({n}^{4})$ growth in the number of significant
four-center two-electron integrals, since the number of non-negligible
product charge distributions, $|\mu \nu \u27e9$, grows as $\mathcal{O}({n}^{2})$. As a result,
the use of large (high-quality) basis expansions is computationally costly.
Perhaps the most practical way around this “basis set quality” bottleneck is
the use of auxiliary basis expansions.
^{
311
}
Chem. Phys. Lett.

(1993),
208,
pp. 359.
Link
^{,}
^{
287
}
Phys. Chem. Chem. Phys.

(2000),
2,
pp. 2113.
Link
^{,}
^{
524
}
Proc. Natl. Acad. Sci. USA

(2005),
102,
pp. 6692.
Link
The ability to use auxiliary
basis sets to accelerate a variety of electron correlation methods, including
both energies and analytical gradients, is a major feature of
Q-Chem.

The auxiliary basis $\{|K\u27e9\}$ is used to approximate products of Gaussian basis functions:

$$|\mu \nu \u27e9\approx |\stackrel{~}{\mu \nu}\u27e9=\sum _{K}|K\u27e9{C}_{\mu \nu}^{K}$$ | (8.2) |

Auxiliary basis expansions were introduced long ago, and are now widely recognized as an effective and powerful approach, which is sometimes synonymously called resolution of the identity (RI) or density fitting (DF). When using auxiliary basis expansions, the rate of growth of computational cost of large-scale electronic structure calculations with $n$ is reduced to approximately ${n}^{3}$.

If $n$ is fixed and molecule size increases, auxiliary basis expansions reduce the pre-factor associated with the computation, while not altering the scaling. The important point is that the pre-factor can be reduced by 5 or 10 times or more. Such large speedups are possible because the number of auxiliary functions required to obtain reasonable accuracy, $X$, has been shown to be only about 3 or 4 times larger than $N$.

The auxiliary basis expansion coefficients, $\mathbf{C}$, are determined by minimizing the deviation between the fitted distribution and the actual distribution, $\u27e8\mu \nu -\stackrel{~}{\mu \nu}|\mu \nu -\stackrel{~}{\mu \nu}\u27e9$, which leads to the following set of linear equations:

$$\sum _{L}\u27e8K|L\u27e9{C}_{\mu \nu}^{L}=\u27e8K|\mu \nu \u27e9$$ | (8.3) |

Evidently solution of the fit equations requires only two- and three-center integrals, and as a result the (four-center) two-electron integrals can be approximated as the following optimal expression for a given choice of auxiliary basis set:

$$\u27e8\mu \nu |\lambda \sigma \u27e9\approx \u27e8\stackrel{~}{\mu \nu}|\stackrel{~}{\lambda \sigma}\u27e9=\sum _{K,L}{C}_{\mu \nu}^{L}\u27e8L|K\u27e9{C}_{\lambda \sigma}^{K}$$ | (8.4) |

In the limit where the auxiliary basis is complete (*i.e.* all products of AOs
are included), the fitting procedure described above will be exact. However,
the auxiliary basis is invariably incomplete (as mentioned above,
$X\approx 3N)$ because this is essential for obtaining increased computational
efficiency.

More details on Q-Chem’s use of RI methods is given in Section 6.6 on RI-MP2 and related methods, Section 6.18 on pairing methods, Section 6.10.6 on coupled cluster methods, Section 4.6.6 on DFT methods, and Section 7.12 on restricted active space methods. In the remainder of this section we focus on documenting the input associated with the auxiliary basis itself.

Q-Chem contains a variety of built-in auxiliary basis sets, that can be
specified by the *$rem* keyword AUX_BASIS.

AUX_BASIS

Sets the auxiliary basis set to be used

TYPE:

STRING

DEFAULT:

No default auxiliary basis set

OPTIONS:

General, Gen
User-defined. As for BASIS
Symbol
Use standard auxiliary basis sets as in the table below
Mixed
Use a combination of different basis sets

RECOMMENDATION:

Consult literature and Basis Set Exchange to aid your selection.

AUX_BASIS_J

Sets the auxiliary basis set for RI-J to be used or invokes RI-J

TYPE:

STRING

DEFAULT:

No default auxiliary basis set

OPTIONS:

General, Gen
User-defined. As for BASIS
Symbol
Use standard auxiliary basis sets as in the table below
Mixed
Use a combination of different basis sets

RECOMMENDATION:

Consult literature and Basis Set Exchange to aid your selection.

AUX_BASIS_K

Sets the auxiliary basis set for RI-K or occ-RI-K to be used or invokes occ-RI-K

TYPE:

STRING

DEFAULT:

No default auxiliary basis set

OPTIONS:

General, Gen
User-defined. As for BASIS
Symbol
Use standard auxiliary basis sets as in the table below
Mixed
Use a combination of different basis sets

RECOMMENDATION:

Consult literature and Basis Set Exchange to aid your selection.

AUX_BASIS_CORR

Sets the auxiliary basis set for RI-MP2 to be used or invokes RI-MP2 in case of double-hybrid DFT or MP2

TYPE:

STRING

DEFAULT:

No default auxiliary basis set

OPTIONS:

General, Gen
User-defined. As for BASIS
Symbol
Use standard auxiliary basis sets as in the table below
Mixed
Use a combination of different basis sets

RECOMMENDATION:

Consult literature and Basis Set Exchange to aid your selection.

Symbolic Name | Atoms Supported |
---|---|

RIJ-def2-XX${}^{a}$ | H, He, Li $\to $ Ne, Na $\to $ Ar, K $\to $ Br, Rb $\to $ Xe, Cs $\to $ La, Hf $\to $ Rn |

RIJK-def2-XX${}^{a}$ | H, He, Li $\to $ Ne, Na $\to $ Ar, K $\to $ Br, Rb $\to $ Xe, Cs $\to $ La, Hf $\to $ Rn |

RIMP2-def2-XX${}^{a}$ | H, He, Li $\to $ Ne, Na $\to $ Ar, K $\to $ Br, Rb $\to $ Xe, Cs $\to $ La, Hf $\to $ Rn |

RIMP2-VDZ${}^{b}$ | H, He, Li $\to $ F, Na $\to $ Cl, K $\to $ Br |

RIMP2-TZVPP
^{
1155
}
Chem. Phys. Lett.(1998), 294, pp. 143. Link |
H, He, Li $\to $ Ne, Na $\to $ Ar, Ga $\to $ Kr |

RIMP2-cc-pVDZ | H, He, Li $\to $ Ne, Na $\to $ Ar, Ga $\to $ Kr |

RIMP2-cc-pVTZ | H, He, Li $\to $ Ne, Na $\to $ Ar, Ga $\to $ Kr |

RIMP2-cc-pVQZ | H, He, Li $\to $ Ne, Na $\to $ Ar, Ga $\to $ Kr |

RIMP2-aug-cc-pVDZ | H, He, B $\to $ Ne, Al $\to $ Ar, Ga $\to $ Kr |

RIMP2-aug-cc-pVTZ | H, He, B $\to $ Ne, Al $\to $ Ar, Ga $\to $ Kr |

RIMP2-aug-cc-pVQZ | H, He, B $\to $ Ne, Al $\to $ Ar, Ga $\to $ Kr |

${}^{a}$: XX includes SV(P), SVP, SVPD (RIMP2 only), TZVP, TZVPD (RIMP2 only), TZVPP, TZVPPD (RIMP2 only), QZVP, QZVPP, QZVPPD (RIMP2 only)

${}^{b}$: This is the original SVP fitting set from Ref. 1155.

In addition to built-in auxiliary basis sets, it is also possible to enter user-defined
auxiliary basis sets using an *$aux_basis* (or *$aux_basis_j*, *$aux_basis_k*, *$aux_basis_corr*)
input section, whose syntax generally follows
the *$basis* input section documented above in Section 8.4.