Electrochemical Impedance Spectroscopy (EIS). The metal-ligand coordinated PDMS samples are sandwiched between two pieces of gold coated SiO 2 substrates with 100 µm Teflon sheet spacers placed on both sides. The gold electrode had an area of 1cm Force field parameters and initial configurations for the PDMS and ions are generated with the Enhanced Monte Carlo software 73 using the PCFF framework. The point charge assigned to the ion pairs is rescaled by 0.7. The charge downscale procedure is commonly used for classical force field where a static point-charge model can lead to an exaggeration of ionic interaction strength. The velocity-Verlet integrator with a time step of 0.5 f s is used for motion evolution; this value is selected as the largest one ensuring the reasonable drift of atoms within a single time step. The relaxation process consists of steepest descent and conjugate gradient energy minimization, followed by equilibration in the canonical (nVT) and isothermal-isobaric (nPT) statistical ensembles at the target temperature and pressure of 303 K and 1 bar, for a total duration of 5 ns. The Nose-Hoover thermostats and barostats are used with damping parameters set to 50 and 500 f s, respectively. The production run is conducted in the nVT ensemble at 303 K for 50 ns. For non-bonded interactions, long-range electrostatic interactions are calculated by using the particle-particle particlemesh solver with a 1.2 nm cutoff distance while Lennard-Jones interactions are truncated at 1.2 nm with the long-range Van der Waals tail corrections included for energy and pressure modification.