Choosing the right electrode representation for modeling real bioelectronic interfaces: a comprehensive guide

Abstract

Objective. Producing realistic numerical models of neurostimulation electrodes in contact with the electrolyte and tissue, for use in time-domain finite element method simulations while maintaining a reasonable computational burden remains a challenge. We aim to provide a straightforward experimental-theoretical hybrid approach for common electrode materials (Ti, TiN, ITO, Au, Pt, IrOx) that are relevant to the research field of bioelectronics, along with all the information necessary to replicate our approach in arbitrary geometry for real-life experimental applications. Approach. We used electrochemical impedance spectroscopy (EIS) to extract the electrode parameters in the AC regime under different DC biases. The pulsed electrode response was obtained by fast amperometry (FA) to optimize and verify the previously obtained electrode parameters in a COMSOL Multiphysics model. For optimization of the electrode parameters a constant phase element (CPE) needed to be implemented in time-domain. Main results. We find that the parameters obtained by EIS can be used to accurately simulate pulsed response only close to the electrode open circuit potential, while at other potentials we give corrections to the obtained parameters, based on FA measurements. We also find that for many electrodes (Au, TiN, Pt, and IrOx), it is important to implement a distributed CPE rather than an ideal capacitor for estimating the electrode double-layer capacitance. We outline and provide examples for the novel time-domain implementation of the CPE for finite element method simulations in COMSOL Multiphysics. Significance. An overview of electrode parameters for some common electrode materials can be a valuable and useful tool in numerical bioelectronics models. A provided FEM implementation model can be readily adapted to arbitrary electrode geometries and used for various applications. Finally, the presented methodology for parametrization of electrode materials can be used for any materials of interest which were not covered by this work.Objective. Producing realistic numerical models of neurostimulation electrodes in contact with the electrolyte and tissue, for use in time-domain finite element method simulations while maintaining a reasonable computational burden remains a challenge. We aim to provide a straightforward experimental-theoretical hybrid approach for common electrode materials (Ti, TiN, ITO, Au, Pt, IrOx) that are relevant to the research field of bioelectronics, along with all the information necessary to replicate our approach in arbitrary geometry for real-life experimental applications. Approach. We used electrochemical impedance spectroscopy (EIS) to extract the electrode parameters in the AC regime under different DC biases. The pulsed electrode response was obtained by fast amperometry (FA) to optimize and verify the previously obtained electrode parameters in a COMSOL Multiphysics model. For optimization of the electrode parameters a constant phase element (CPE) needed to be implemented in time-domain. Main results. We find that the parameters obtained by EIS can be used to accurately simulate pulsed response only close to the electrode open circuit potential, while at other potentials we give corrections to the obtained parameters, based on FA measurements. We also find that for many electrodes (Au, TiN, Pt, and IrOx), it is important to implement a distributed CPE rather than an ideal capacitor for estimating the electrode double-layer capacitance. We outline and provide examples for the novel time-domain implementation of the CPE for finite element method simulations in COMSOL Multiphysics. Significance. An overview of electrode parameters for some common electrode materials can be a valuable and useful tool in numerical bioelectronics models. A provided FEM implementation model can be readily adapted to arbitrary electrode geometries and used for various applications. Finally, the presented methodology for parametrization of electrode materials can be used for any materials of interest which were not covered by this work.