Looking at both MAPE and 𝑅 2 , we can observe how the general trends of Fig. are replicated for all the methods at different thresholds. However, it is noteworthy that the MAPE (𝑅 2 ) decreases (increases) when increasing the 𝛿 TPA Q-CC3 threshold. Indeed, looking at Q-CCSD (E-CCSD) results, we can observe that it goes from 25.9% (29.6%) for the 𝑡 1 interval, to 10.4% (16.9%) for the 𝑡 2 interval, and to 5.4% (8.1%) for the 𝑡 3 interval. Similar decreases are obtained for the other methods (see Tables S11-S13 in the SI). Analogous trends can be found for the correlation. For the states in 𝑡 1 , Q-CCSD and E-CCSD are the only ones showing 𝑅 2 larger than 0.5 (0.84 and 0.82 respectively), while in 𝑡 2 all methods have a correlation larger than 0.4 (in which Q-CCSD excell with an 𝑅 2 of 0.85). On the other hand in the 𝑡 3 interval all methods show large correlation (𝑅 2 > 0.91). It is therefore a good news that as the molecules get more interesting for practical applications all methods become more trustworthy.