The effects of polarization on the rotational diffusion of ions in organic ionic plastic crystals

Organic ionic plastic crystals (OIPCs), which consist of organic molecular ions, are considered excellent candidates for solid electrolytes due to their high ionic conductivity in solid phases. Molecular ions undergo either rotational or conformational relaxation at certain temperatures in OIPCs. There have been molecular simulations to understand the rotational motion. The polarizability of ions was, however, often ignored in simulations due to the high computational cost. Since the polarizability may affect the translational diffusion, the ionic conductivity, and the phase transition of ionic liquids, it should be of interest to investigate how the polarizability would affect the rotational diffusion of ions in solid phases. In this work, we perform extensive atomistic molecular dynamics simulations for two different kinds of OIPCs, 1-methyl-3-methylimidazolium hexafluorophosphate ([MMIM][PF₆]) and 1-butyl-3-methylimidazolium hexafluorophosphate ([BMIM][PF₆]). We employ various simulation models for ions by turning on and off the polarization in their interaction potentials. We find that the polarizability hardly affects the density, the crystalline structure, and the phase transition of both OIPCs. However, a certain rotational motion, especially the rotational diffusion of PF₆ ^- in [MMIM][PF₆] OIPCs, is enhanced by a factor of up to four when the polarizability is turned on. The PF₆ ^- in [MMIM][PF₆] OIPCs undergoes rotational hopping motions more significantly due to polarizability. We find that the rotational diffusion of a certain ion can be enhanced only when the polarization results in a significant change in the dipole moment of the neighboring ions around the ion.