In this work, we have deployed two alternative strategies to account for the populations of different phenol aggregates, i.e., the Q/s and C/d methodologies described in Sections 2.2 and 2.3. In the following sections, we first introduce and compare the populations of aggregates provided by each of them and eventually use them to compute the average IR spectra, focusing on the O H stretching region.