We present a perspective on the use of diffuse basis functions for electronic structure calculations by density functional theory and wave function theory. We especially emphasize minimally augmented basis sets and calendar basis sets. We base our conclusions on our previous experience with commonly computed quantities, such as bond energies, barrier heights, electron affinities, noncovalent (van der Waals and hydrogen bond) interaction energies, and ionization potentials, on Stephens et al.'s results for optical rotation and on our own new calculations (presented here) of polarizabilities and of potential energy curves of van der Waals complexes. We emphasize the benefits of partial augmentation of the higher-zeta basis sets in preference to full augmentation at a lower ζ level. Benefits and limitations of the use of fully, partially, and minimally augmented basis sets are reviewed for different electronic structure methods and molecular properties. We have found that minimal augmentation is almost always enough for density functional theory (DFT) when applied to ionization potentials, electron affinities, atomization energies, barrier heights, and hydrogen-bond energies. For electric dipole polarizabilities, we find that augmentation beyond minimal has an average effect of 8% at the polarized triple-ζ level and 5% at the polarized quadruple-ζ level. The effects are larger for potential energy curves of van der Waals complexes. The effects are also larger for wave function theory (WFT). Even for WFT though, full augmentation is not needed for most purposes, and a level of augmentation between minimal and full is optimal for most problems. The calendar basis sets named after the months provide a convergent sequence of partially augmented basis sets that can be used for such calculations. The jun-cc-pV(T+d)Z basis set is very useful for MP2-F12 calculations of barrier heights and hydrogen bond strengths.

Many methods for correcting harmonic partition functions for the presence of torsional motions employ some form of one-dimensional torsional treatment to replace the harmonic contribution of a specific normal mode. However, torsions are often strongly coupled to other degrees of freedom, especially other torsions and low-frequency bending motions, and this coupling can make assigning torsions to specific normal modes problematic. Here, we present a new class of methods, called multi-structural (MS) methods, that circumvents the need for such assignments by instead adjusting the harmonic results by torsional correction factors that are determined using internal coordinates. We present three versions of the MS method: (i) MS-AS based on including all structures (AS), i.e., all conformers generated by internal rotations; (ii) MS-ASCB based on all structures augmented with explicit conformational barrier (CB) information, i.e., including explicit calculations of all barrier heights for internal-rotation barriers between the conformers; and (iii) MS-RS based on including all conformers generated from a reference structure (RS) by independent torsions. In the MS-AS scheme, one has two options for obtaining the local periodicity parameters, one based on consideration of the nearly separable limit and one based on strongly coupled torsions. The latter involves assigning the local periodicities on the basis of Voronoi volumes. The methods are illustrated with calculations for ethanol, 1-butanol, and 1-pentyl radical as well as two one-dimensional torsional potentials. The MS-AS method is particularly interesting because it does not require any information about conformational barriers or about the paths that connect the various structures.

We combine the diffuse basis functions from the 6-31+G basis set of Pople and co-workers with the correlation-consistent basis sets of Dunning and co-workers. In both wave function and density functional calculations, the resulting basis sets reduce the basis set superposition error almost as much as the augmented correlation-consistent basis sets, although they are much smaller. In addition, in density functional calculations the new basis sets, called cc-pVxZ+ where x = D, T, Q, ..., or x = D+d, T+d, Q+d, ..., give very similar energetic predictions to the much larger aug-cc-pVxZ basis sets. However, energetics calculated from correlated wave function calculations are more slowly convergent with respect to the addition of diffuse functions. We also examined basis sets with the same number and type of functions as the cc-pVxZ+ sets but using the diffuse exponents of the aug-cc-pVxZ basis sets and found very similar performance to cc-pVxZ+; these basis sets are called minimally augmented cc-pVxZ, which we abbreviate as maug-cc-pVxZ.

We present sets of convergent, partially augmented basis set levels corresponding to subsets of the augmented "aug-cc-pV(n+d)Z" basis sets of Dunning and co-workers. We show that for many molecular properties a basis set fully augmented with diffuse functions is computationally expensive and almost always unnecessary. On the other hand, unaugmented cc-pV(n+d)Z basis sets are insufficient for many properties that require diffuse functions. Therefore, we propose using intermediate basis sets. We developed an efficient strategy for partial augmentation, and in this article, we test it and validate it. Sequentially deleting diffuse basis functions from the "aug" basis sets yields the "jul", "jun", "may", "apr", etc. basis sets. Tests of these basis sets for Møller-Plesset second-order perturbation theory (MP2) show the advantages of using these partially augmented basis sets and allow us to recommend which basis sets offer the best accuracy for a given number of basis functions for calculations on large systems. Similar truncations in the diffuse space can be performed for the aug-cc-pVxZ, aug-cc-pCVxZ, etc. basis sets.

Eliminating all but the s and p diffuse functions on the non-hydrogenic atoms and all diffuse functions on the hydrogen atoms from the aug-cc-pV(x+d)Z basis sets of Dunning and co-workers, where x = D, T, Q, ..., yields the previously proposed "minimally augmented" basis sets, called maug-cc-pV(x+d)Z. Here, we present extensive and systematic tests of these basis sets for density functional calculations of chemical reaction barrier heights, hydrogen bond energies, electron affinities, ionization potentials, and atomization energies. The tests show that the maug-cc-pV(x+d)Z basis sets are as accurate as the aug-cc-pV(x+d)Z ones for density functional calculations, but the computational cost savings are a factor of about two to seven.