Similarly, anions were tethered to both surfaces with a grafting density equal to that of the polymer chains. Each solvent molecule was represented as an oligomer with 5 monomer beads. To account for the interaction of ions with the oligomeric solvent and polymer chains, we employ a polymeric Stockmayer model in which a freely rotating dipole moment of constant strength µ is embedded into a finite-size particle, while other bonded and nonbonded interactions are treated using the Kremer-Grest polymer model. The outputs of the MD simulations are used to characterize three features: (i) the effective pair potential U (r) between functionalized spherical nanoparticles, (ii) local cation diffusion coefficients in the direction parallel to the flat surface (D ), and (iii) the spatial distribution of cations in the region between the nanoparticles.