Hierarchical Electrochemical Modeling and Simulation of Bio-Hybrid Interfaces

Emanuela Abbate,Matteo Porro,Thierry Nieus,Riccardo Sacco
DOI: https://doi.org/10.1016/j.cma.2015.11.024
2015-01-19
Abstract:In this article we propose and investigate a hierarchy of mathematical models based on partial differential equations (PDE) and ordinary differential equations (ODE) for the simulation of the biophysical phenomena occurring in the electrolyte fluid that connects a biological component (a single cell or a system of cells) and a solid-state device (a single silicon transistor or an array of transistors). The three members of the hierarchy, ordered by decreasing complexity, are: (i) a 3D Poisson-Nernst-Planck (PNP) PDE system for ion concentrations and electric potential; (ii) a 2D reduced PNP system for the same dependent variables as in (i); (iii) a 2D area-contact PDE system for electric potential coupled with a system of ODEs for ion concentrations. The backward Euler method is adopted for temporal semi-discretization and a fixed-point iteration based on Gummel's map is used to decouple system equations. Spatial discretization is performed using piecewise linear triangular finite elements stabilized via edge-based exponential fitting. Extensively conducted simulation results are in excellent agreement with existing analytical solutions of the PNP problem in radial coordinates and experimental and simulated data using simplified lumped parameter models.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to accurately model and simulate electrochemical phenomena in bio - solid hybrid interfaces in the field of bioelectronics. Specifically, the research focuses on the mathematical modeling of biophysical phenomena occurring in electrolyte fluids (connecting biological components such as single cells or cell systems and solid - state devices such as single silicon transistors or transistor arrays). ### Research Background and Problem Description Bio - electronic interfaces are bio - hybrid structures consisting of living cells attached to an electronic substrate and surrounded by an electrolyte bath. These interfaces have important applications in devices such as neurochips, where ionic currents caused by cell activity are conducted through the cell membrane and generate voltages outside the cell, which in turn affect the charge flow on the substrate. Therefore, understanding and accurately simulating the electrochemical processes at these interfaces are crucial for the development of new bio - electronic devices. ### Main Problems 1. **Multiscale Modeling**: Since a planar electrical conductor is formed at the cell - chip junction, the main physical phenomena (such as ionic electro - diffusion and gate voltage modulation) occur within this three - dimensional region. To accurately describe these phenomena, a multiscale perspective is required. 2. **Balance between Complexity and Computational Efficiency**: Although the 3D Poisson - Nernst - Planck (PNP) model is the most accurate, it has a high computational cost. Therefore, simplified models need to be developed to improve computational efficiency while ensuring accuracy. ### Proposed Solutions The authors propose a hierarchical model system based on partial differential equations (PDEs) and ordinary differential equations (ODEs), including: 1. **3D PNP System**: \[ \begin{aligned} &\frac{\partial c_i}{\partial t}+\nabla\cdot f_i(c_i,\phi) = 0, \\ &f_i(c_i,\phi)= - D_i\nabla c_i+\mu_i z_i|z_i|c_i E, \\ &D_i=\mu_i V_{th}|z_i|, \\ &\nabla\cdot E=\frac{1}{\epsilon}\rho, \\ &E = -\nabla\phi, \\ &\rho = q\sum_{i = 1}^M z_i c_i. \end{aligned} \] 2. **2D Reduced PNP System**: By averaging the 3D PNP equations in the direction perpendicular to the electrolyte gap, a 2D PNP system is obtained. 3. **2D Surface Contact Model**: Further simplification, ignoring the spatial dependence of ion concentrations in the electrolyte gap, results in a time - dependent 2D Poisson equation coupled with an ODE system describing ion dynamics. Through these hierarchical models, the authors are able to find an appropriate balance between different levels of complexity and computational cost, thereby more efficiently simulating the electrochemical phenomena in bio - solid hybrid interfaces.