CMD RDFs obtained with 16 molecules of each molecular compound with the extended version of the CL&Pol polarizable force field by Goloviznina et al. (named FF1 hereafter) agree in the hydrogen bond distances, but with discrepancies in the g(r) for ChCl:EG (in blue in Figure ). Nevertheless, while in their system the interactions were underestimated by the original FF, in our case the interaction with the cation (Cl -ChOH H-bond) is overestimated and those with the neutral molecule (Cl -LAOH and Cl -LACOOH H-bond) are underestimated. This is specially striking for Cl -LAOH RDF, where at the first peak, the g(r) value reaches 0.58, considerably lower than what a non-correlated pair of particles would exhibit. We tune the values of the parameters for these particular interparticle interactions taking as reference the RDFs obtained by AIMD. The modified force field that includes these tuned parameters is named FF2 hereafter. The original and modified values of the parameters are shown in Table 1 . The agreement between AIMD and FF2 (see RDFs in red in Figure ) is very satisfactory.