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)₆]^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^® 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 µM to 100 µ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.