This Appendix is adopted from the Ph.D. Thesis by Alston Misquitta [66].

The following procedure can be applied to construct an auxiliary basis set for each atom in the dimer under consideration. Denoted by is the decontracted basis set used in calculations for given dimer. The auxiliary basis, denoted by , is developed as follows:

- Construct an auxiliary basis as
the tensor product of with itself. That is,
(15) If G

^{}(l_{i},α_{i}) and G^{}(l_{j},α_{j}) are two basis functions of centered at the same point, where l_{i}and l_{j}are the angular quantum numbers and α_{i}and α_{j}are the exponents, then the product is a basis function belonging to centered at the same point and given by G^{}(l_{k},α_{k}) where l_{k}= l_{i}+ l_{j}and α_{k}= α_{i}+ α_{j}. The resulting basis is a (large) basis including high symmetry functions compared to the original basis. All the products involving basis functions from different centers are rejected. - Within each angular symmetry of
, a reduction is performed in the number of basis
functions as follows. Given an ϵ, if
there are n basis functions for which log
α
_{k1}, log α_{k2},...,log α_{kn}are in an ϵ neighborhood, then these n functions are replaced by one function with the exponent β = (α_{k1}α_{k2}...α_{kn})^{1∕n}. Perform this reduction for the whole basis set. - If necessary, the resulting basis can be reduced even further with the previous step repeated on the reduced basis with a different value of ϵ.
- Further pruning of the reduced
basis can be done as follows:
- Reject all functions of g-symmetry and higher. This is necessary if cadpac is used to obtain integrals.
- The ‘large’ exponents, particularly the ones of high symmetries, can generally be removed. The criterion for this removal is best found by trial and error and by monitoring the constraint conditions.

We have found the optimal values of ϵ used in the pruning process to be between 0.3 and 0.5. This procedure typically results in an auxiliary basis that is 2 to 3 times larger than the basis used to obtain molecular orbitals and eigenvalues. While in general the auxiliary basis for a molecule should be centered also on sites between pairs of atoms in addition to atomic sites, it is our experience that the use of only atomic centers is adequate. This is also true for the optimized auxiliary basis sets [14].