Excited states exhibiting double excitation character are notoriously di cult to model using conventional singlereference methods, such as adiabatic time-dependent density-functional theory (TD-DFT) or equation-of-motion coupled cluster (EOM-CC). In addition, these states are typical experimentally "dark" making their detection in photo-absorption spectra very challenging. Nonetheless, they play a key role in the faithful description of many physical, chemical, and biological processes. In the present work, we provide accurate reference excitation energies for transitions involving a substantial amount of double excitation using a series of increasingly large di use-containing atomic basis sets. Our set gathers 20 vertical transitions from 14 small-and medium-size molecules (acrolein, benzene, beryllium atom, butadiene, carbon dimer and trimer, ethylene, formaldehyde, glyoxal, hexatriene, nitrosomethane, nitroxyl, pyrazine, and tetrazine). Depending on the size of the molecule, selected con guration interaction (sCI) and/or multicon gurational (CASSCF, CASPT2, (X)MS-CASPT2 and NEVPT2) calculations are performed in order to obtain reliable estimates of the vertical transition energies. In addition, coupled cluster approaches including at least contributions from iterative triples (such as CC3, CCSDT, CCSDTQ, and CCSDTQP) are assessed. Our results clearly evidence that the error in CC methods is intimately related to the amount of double excitation character of the transition. For "pure" double excitations (i.e. for transitions which do not mix with single excitations), the error in CC3 can easily reach 1 eV, while it goes down to few tenths of an eV for more common transitions (like in trans-butadiene) involving a signi cant amount of singles. As expected, CC approaches including quadruples yield highly accurate results for any type of transitions. e quality of the excitation energies obtained with multicon gurational methods is harder to predict. We have found that the overall accuracy of these methods is highly dependent of both the system and the selected active space. e inclusion of the σ and σ orbitals in the active space, even for transitions involving mostly π and π orbitals, is mandatory in order to reach high accuracy. A theoretical best estimate (TBE) is reported for each transition. We believe that these reference data will be valuable for future methodological developments aiming at accurately describing double excitations. TOC graphical abstract arXiv:1811.12861v1 [physics.chem-ph] 30 Nov 2018