Multiphysics Modeling of Electrochemical Impedance Spectroscopy Responses of SAM-Modified Screen-Printed Electrodes

In this work, we present a multiphysics modeling approach capable of simulating electrochemical impedance spectroscopy (EIS) responses of screen-printed electrodes (SPEs) modified with self-assembled monolayers of 11-Mercaptoundecanoic acid (MUA). Commercially available gold SPEs are electrochemically characterized through experimental cyclic voltammetry and EIS measurements with 10 mM [Fe(CN)6]3-/4- redox couple in phosphate buffered saline before and after the surface immobilization of MUA at different concentrations. We design the multiphysics model through COMSOL Multiphysics (R) based on the 3D geometry of the devices under test. The model includes four different physics considering the metal/solution interface electrochemical phenomena, the ion and electron potentials and currents, and the measurement set-up. The model is calibrated through a set of experimental measurements, allowing the tuning of the parameters used by the model. We use the calibrated model to simulate the EIS response of MUA-modified SPEs, comparing the results with experimental data. The simulations fit the experimental curves well, following the variation of MUA concentration on the surface from 1 mu M to 100 mu M. The EIS parameters, retrieved through a CPE-modified Randles' circuit, confirm the consistency with the experimental data. Notably, the simulated surface coverage estimates and the variation of charge transfer resistance due to MUA-immobilization are well matched with their experimental counterparts, reporting only a 2% difference and being consistent with the experimental electrochemical behavior of the SPEs.