Concerning the remaining hydrogen bonds in the system, a summary is provided in Figure . Overall, FF1 tends to underestimate the strength of the hydrogen bond network in comparison to DFT, with a unique exception (Figure ). Regarding the influence of the modified force fields in the H-bond network, two distinct scenarios emerge: those in which LA acts as a hydrogen bond donor, interacting with either choline or another LA ((a) and (b) in Figure ). For these specific H-bonds, the corrections to the original force field slightly increase the deviation from AIMD calculations, despite the three force fields exhibiting similar results. This behavior is likely attributed to the fact that the parameter correction strengthens both Cl -LAOH and Cl -LACOOH interactions, resulting in fewer available LAs to directly coodinate with other LAs and cholines.