The enthalpic contribution is primarily dictated by intermolecular interactions, including hydrogen bonds and π-π interactions. In this sense, it is expected that the dispersioncorrected DFT calculations in the Q/s approach will outperform the force field calculations applied in C/d. Indeed, the inaccuracies in the FF could partly explain some of the deficiencies in C/d, such as its apparent underestimation of tetramer populations. Therefore, enthalpy differences cannot explain the better performance of the C/d strategy. The key to understand this difference would necessarily lie in the entropic term, which is linked to the vast number of configurations arising from the significant molecular flexibility in loosely bound aggregates. This implies that many configurations are necessary to achieve population distributions achieved by the experimental system, which becomes unmanageable for the automated conformational search of the Q/s approach. Conversely, this conformational flexibility is, by construction, effectively captured by the MD simulation adopted in the C/d strategy.