Distinct dissociation kinetics between ion pairs: Solvent-coordinate free-energy landscape analysis

Different ion pairs exhibit different dissociation kinetics; however, while the nature of this process is vital for understanding various molecular systems, the underlying mechanism remains unclear. In this study, to examine the origin of different kinetic rate constants for this process, molecular dynamics simulations were conducted for LiCl, NaCl, KCl, and CsCl in water. The results showed substantial differences in dissociation rate constant, following the trend kLiCl < kNaCl < kKCl < kCsCl. Analysis of the free-energy landscape with a solvent reaction coordinate and subsequent rate component analysis showed that the differences in these rate constants arose predominantly from the variation in solvent-state distribution between the ion pairs. The formation of a water-bridging configuration, in which the water molecule binds to an anion and a cation simultaneously, was identified as a key step in this process: water-bridge formation lowers the related dissociation free-energy barrier, thereby increasing the probability of ion-pair dissociation. Consequently, a higher probability of water-bridge formation leads to a higher ion-pair dissociation rate.