Abstract:Objective.Increasing complexity in extracellular stimulation experiments and neural implant design also requires realistic computer simulations capable of modeling the neural activity of nerve cells under the influence of an electrical stimulus. Classical model approaches are often based on simplifications, are not able to correctly calculate the electric field generated by complex electrode designs, and do not consider electrical effects of the cell on its surrounding. A more accurate approach is the finite element method (FEM), which provides necessary techniques to solve the Poisson equation for complex geometries under consideration of electrical tissue properties. Especially in situations where neurons experience large and non-symmetric extracellular potential gradients, a FEM solution that implements the cell membrane model can improve the computer simulation results. To investigate the response of neurons in an electric field generated by complex electrode designs, a FEM framework for extracellular stimulation was developed in COMSOL.Approach.Methods to implement morphologically- and biophysically-detailed neurons including active Hodgkin-Huxley (HH) cell membrane dynamics as well as the stimulation setup are described in detail. Covered methods are (a) development of cell and electrode geometries including meshing strategies, (b) assignment of physics for the conducting spaces and the realization of active electrodes, (c) implementation of the HH model, and (d) coupling of the physics to get a fully described model.Main results.Several implementation examples are briefly presented: (a) a full FEM implementation of a HH model cell stimulated with a honeycomb electrode, (b) the electric field of a cochlear electrode placed inside the cochlea, and (c) a proof of concept implementation of a detailed double-cable cell membrane model for myelinated nerve fibers.Significance.The presented concepts and methods provide basic and advanced techniques to realize a full FEM framework for innovative studies of neural excitation in response to extracellular stimulation.

Cell Electropermeabilization Modeling via Multiple Traces Formulation and Time Semi-Implicit Coupling

Application of Efficient Numerical Methods in Solution of Ordinary Differential Equations for Modeling Electrical Activity in Cardiac Cells

Boundary integral formulation of the cell-by-cell model of cardiac electrophysiology

A splitting, discontinuous Galerkin solver for the cell-by-cell electroneutral Nernst-Planck framework

Scalable approximation and solvers for ionic electrodiffusion in cellular geometries

Accurate numerical simulation of electrodiffusion and water movement in brain tissue

Efficient Numerical Methods for Solving the Membrane Equations of Human Ventricular Cell Model

Cut finite element discretizations of cell-by-cell EMI electrophysiology models

Finite Element Simulation of Ionic Electrodiffusion in Cellular Geometries

Cell membrane electroporation modeling: A multiphysics approach

Hierarchical Electrochemical Modeling and Simulation of Bio-Hybrid Interfaces

Modeling Excitable Cells with the EMI Equations: Spectral Analysis and Iterative Solution Strategy

A Discontinuous Galerkin Method for the Extracellular Membrane Intracellular Model

Finite element simulation of ionicelectrodiffusion in cellular geometries

A Boundary Element Method of Bidomain Modeling for Predicting Cellular Responses to Electromagnetic Fields

Dense cell-by-cell systems of PDEs: approximation, spectral analysis, and preconditioning

A finite element method framework to model extracellular neural stimulation

Numerical study on the effect of capacitively coupled electrical stimulation on biological cells considering model uncertainties