In this study, we carried out molecular dynamics simulations to examine the molecular mechanism for ionic liquid pair association in aqueous solutions. We chose the commonly studied imidazolium-based ionic liquid pairs. We computed potentials of mean force (PMF) for four systems: 1,3-dimethlylimidazoliumchloride, 1,3-dimethlylimidazolium iodide, 1-methly-3-octylimidazolium chloride, and 1-methly-3-octylimidazolium iodide. Our PMF studies show a stronger interaction for the ion pairs of systems involving dimethlylimidazolium as the cation species compared with that of the systems containing octylimidazolium. This result indicates a decrease in ion-pair association as the cation alkyl tail length increases. We also studied the kinetics of ion-pair dissociation using different rate theories such as the Grote-Hynes and Kramer's theories. As expected, the computed rate results significantly deviated from results obtained from transition state theory because it does not account for dynamical solvent effects. Dissociative barrier curvatures are found to be very small for the systems investigated because the transmission coefficients computed using Grote-Hynes theory and Kramer's theory are approximately equal. Our analysis of the rotational dynamics of cations revealed that the time scales for molecular reorientation are longer for cations with longer alkyl tails.