Abstract:Objective: Commonly used cable equation-based approaches for determining the effects of electromagnetic fields on excitable cells make several simplifying assumptions that could limit their predictive power. Bidomain or "whole" finite element methods have been developed to fully couple cells and electric fields for more realistic neuron modeling. Here, we introduce a novel bidomain integral equation designed for determining the full electromagnetic coupling between stimulation devices and the intracellular, membrane, and extracellular regions of neurons. Methods: Our proposed boundary element formulation offers a solution to an integral equation that connects the device, tissue inhomogeneity, and cell membrane-induced E-fields. We solve this integral equation using first-order nodal elements and an unconditionally stable Crank-Nicholson time-stepping scheme. To validate and demonstrate our approach, we simulated cylindrical Hodgkin-Huxley axons and spherical cells in multiple brain stimulation scenarios. Main results: Comparison studies show that a boundary element approach produces accurate results for both electric and magnetic stimulation. Unlike bidomain finite element methods, the bidomain boundary element method does not require volume meshes containing features at multiple scales. As a result, modeling cells, or tightly packed populations of cells, with microscale features embedded in a macroscale head model, is made computationally tractable, and the relative placement of devices and cells can be varied without the need to generate a new mesh. Significance: Device-induced electromagnetic fields are commonly used to modulate brain activity for research and therapeutic applications. Bidomain solvers allow for the full incorporation of realistic cell geometries, device E-fields, and neuron populations. Thus, multi-cell studies of advanced neuronal mechanisms would greatly benefit from the development of fast-bidomain solvers to ensure scalability and the practical execution of neural network simulations with realistic neuron morphologies.

Efficient solution of bidomain equations in simulation of cardiac excitation anisotropic propagation

Parallel Solution In Simulation Of Cardiac Excitation Anisotropic Propagation

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

Progress in ion-channel-based myocardium modeling and simulation

Towards the Cardiac Equivalent Source Models in Electrocardiography and Magnetocardiography: A Simulation Study

Simulation of Biatrial Conduction Via Different Pathways During Sinus Rhythm with a Detailed Human Atrial Model

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

One-dimensional Simulation of Transmural Heterogeneity of Cardiac Cellular Electromechanics

Preconditioning the bidomain model with almost linear complexity

High-order parallel-in-time method for the monodomain equation in cardiac electrophysiology

Simulation of Active Cardiac Electromechanical Dynamics

An Effective Finite Difference Method for Simulation of Bioheat Transfer in Irregular Tissues

A Highly Efficient and Accurate Algorithm for Solving the Par Tial Differential Equation in Cardiac Tissue Models

Numerical study of POD-Galerkin-DEIM reduced order modeling of cardiac monodomain formulation

Integration of Different Cardiac Electrophysiological Models into a Single Simulation Pipeline

A multi-order smoothed particle hydrodynamics method for cardiac electromechanics with the Purkinje network

An explicit local space-time adaptive framework for monodomain models in cardiac electrophysiology

Algorithm Study on the Three-dimensional Cardiac Tissue Based on the Model of Ventricular Action Potential

An Efficient Parallel Numerical Modeling of Bioheat Transfer in Realistic Tissue Structure

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

Numerical Study on the Mathematical Model of the Cardiac Electrical Conduction