Force field molecular dynamics simulations of the same system were employed to preequilibrate the system and generate initial configurations for the AIMD simulations in different ion pair geometries. After minimization, equilibration in the NPT ensemble was done for 0.5 ns to determine the average size of the simulation cell, followed by a 20 ns simulation in the NPT ensemble. We extracted from there 20 representative initial configurations of the ion pair.