As an example of a "real-world" application where the HF+D approach might be preferred to DFT+D, we consider a data set of ion pairs corresponding to common constituents of room-temperature ionic liquids. For such systems, Grimme et al. report that DFT calculations sometimes face convergence issues related to delocaliza- tion error and unphysical charge transfer, at least when semilocal functionals are employed. The HF+D(EFP) method provides a potentially useful alternative in such cases, and indeed could be used as a general way to probe whether delocalization error is problematic in a particular case, without sacrificing the quantitatively important dispersion interaction. The data set assembled for this work is a subset of the IL-2013 database and consists of 195 unique ion pairs. We shift each dimer along the vector connecting the cation and anion centers of mass, in order to generate eight different geometries for each dimer corresponding to intermolecular separations ranging from 0.9R e to 2.0R e , where R e is the equilibrium separation from the IL-2013 database. We call this data set IL195×8, and it consists of 1,560 individual structures whose interaction energies we have computed at the level of coupled-cluster theory with single, double, and perturbative triple excitations [CCSD(T)], using the domain-localized pair natural orbital (DLPNO) implementation of CCSD(T) in the ORCA program. Tight cutoffs are used because the accuracy of the DLPNO approximation is known to be quite sensitive to these thresholds for noncovalent interaction energies. DLPNO-CCSD(T) interaction energies were extrapolated to the complete basis-set (CBS) limit using a two-point formula with aug-cc-pVTZ and aug-cc-pVQZ. Additional details are provided in the Supplementary Material.