State of hydrophobic and hydrophilic ionic liquids in aqueous solutions: are the ions fully dissociated?

Molecular dynamics simulations were performed for aqueous solutions of five ionic liquids (ILs): 1-ethyl-3-methylimidazolium ([C2mim]) bis(trifluoromethanesulfonyl) imide ([NTf2]), 1-n-butyl-3-methylimidazolium ([C4mim]) [NTf2], 1-n-hexyl-3-methylimidazolium ([C6mim]) [NTf2], [C2mim] ethylsulfate ([C2H5SO4]), and [C2mim] chloride (Cl) to determine whether the ions of these ILs are associated at relatively high dilutions and whether the association is governed by hydrophobicity/hydrophilicity of the ILs. The adaptive biasing force technique was applied to calculate the potential of mean force (PMF) for each IL ion pair. For all of the ILs, the PMF is characterized by two distinct contact minima in which the ions have different relative conformations. The hydrophobic ILs bearing the anion [NTf2](-) exist predominantly in the associative state; the strength of the association of these ILs increases with increase in the alkyl chain length. The most hydrophilic IL [C2mim] Cl was determined to be almost fully dissociated at the concentration examined in the study. [C2mim] [C2H5SO4] showed hydration behavior that was intermediate between that exhibited by the ILs in which the anion is substituted with either Cl(-) or [NTf2](-) paired with [C2mim](+). Association constants for these ILs were also computed. Radial distribution functions calculated by constraining the ions at the contact minima showed that hydration of the anion plays the dominant role in determining the microscopic behavior of these ILs in aqueous solutions.