Starting from the over 115,000 molecules in the FORMED database and excluding the molecules we had already screened in our previous study, we evaluate 4070 molecules at the CC2/aug-cc-pVDZ//GFN2-xTB level based on their size and their E(S 1 -T 1 ) computed at the ωB97X-D/6-31G(d)//GFN2-xTB using the improved workflow. All optimizations were carried out using Gaussian16, while CC2 computations used Turbomole 7.1. We select molecules with E(S 1 -T x ) < 0.01 eV for indepth exploration (T x being the triplet state with same electron configuration as S 1 ). Finally, we compute the excitation energies for leading candidates at the equations-of-motion coupled-cluster level (EOM-CCSD) as implemented in the QChem 5.1 package. We used CC2 as the workhorse in the screening process as it is more computationally efficient than EOM-CCSD. Based on the data in this and our previous studies, CC2 provides an excellent compromise between accuracy and efficiency for the initial screening steps because it systematically underestimates the singlet-triplet gaps (Figure ). This way, we avoid false negatives while the false positives will be removed at the final step that is carried out at the EOM-CCSD level (Figure ).