AMOEBA-IL follows the AMOEBA procedure to parametrize the bonded terms. The non-bonded terms are initially parametrized by comparing individual inter-molecular interaction contributions of gas-phase dimers to QM energy decomposition analysis (QM-EDA) reference values as described in. In the current work, we have used DFT-based Symmetry Adapted Perturbation Theory (SAPT-DFT) to obtain the reference QM inter-molecular interactions. Once the initial non-bonded terms are parametrized, bulk simulations are carried out to calculate densities and heats of vaporization to determine the accuracy of the parameters, and if needed, adjust the Halgren parameters to optimize agreement for both gas-phase and bulk properties. As mentioned above, the distributed multipoles are obtained by using GEM-DM. The procedure to obtain GEM-DM has been discussed in detail previously in Refs. Briefly, GEM-DM multipoles are extracted from fitted GEM densities directly. The GEM densities are fitted by minimizing the error of the inter-molecular Coulomb interaction for a series of dimers, here composed of the initial [Im + ]-based cation structures and an already parametrized water molecule 34 as a function of the distance between the monomers (cation position fixed) on the same axis with increments of 0.5 Å from one dimer to the next one. The inter-molecular Coulomb interaction reference values were obtained using the symmetry-adapted perturbation theory (SAPT) method implemented in Psi4. Validation for the calculated total energies of the parametrized fragments was done by comparing counterpoise corrected total inter-molecular interactions computed with Gaussian16 package at the MP2(full)/aug-cc-PVTZ level of theory. Additionally, the total inter-molecular interaction at the second-order perturbation was calculated with SAPT as implemented in Psi4. On the other hand, the Coulomb and total inter-molecular energies reproduced by the developed parameters were calculated with the analyze program of TINKER. Comparison between structure minimization with GEM-DM and