Vibrational spectroscopy is a fundamental tool to investigate local atomic arrangements and the effect of the environment, provided that the spectral features can be correctly assigned. This can be challenging in experiments and simulations when double peaks are present because they can have different origins. Fermi dyads are a common class of such doublets, stemming from the resonance of the fundamental excitation of a mode with the overtone of another. We present a new, efficient approach to unambiguously characterize Fermi resonances in density functional theory (DFT) based simulations of condensed phase systems. With it, the spectral features can be assigned and the two resonating modes identified. We also show how data from DFT simulations employing classical nuclear dynamics can be post-processed and combined with a perturbative quantum treatment at a finite temperature to include analytically thermal quantum nuclear effects. The inclusion of these effects is crucial to correct some of the qualitative failures of the Newtonian dynamics simulations at a low temperature such as, in particular, the behavior of the frequency splitting of the Fermi dyad. We show, by comparing with experimental data for the paradigmatic case of supercritical CO2, that these thermal quantum effects can be substantial even at ambient conditions and that our scheme provides an accurate and computationally convenient approach to account for them. Published by AIP Publishing.