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

Area of Science:

• Computational Neuroscience
• Biophysics
• Electrophysiology

Background:

• Traditional cable equation models simplify electromagnetic field effects on excitable cells, limiting predictive accuracy.
• Bidomain finite element methods offer more realistic neuron modeling by coupling cells and electric fields.
• Accurate modeling of electromagnetic field interactions is crucial for brain stimulation research and therapies.

Approach:

• Introduced a novel bidomain integral equation formulation for comprehensive electromagnetic coupling analysis.
• Utilized first-order nodal elements and a Crank-Nicholson time-stepping scheme for solving the integral equation.
• Validated the approach through simulations of Hodgkin-Huxley axons and spherical cells in various brain stimulation scenarios.

Key Points:

• The boundary element method accurately predicts electric and magnetic stimulation effects.
• Unlike finite element methods, it avoids multi-scale volume meshing, simplifying complex models.
• Enables computationally tractable modeling of cells with microscale features in macroscale models.

Conclusions:

• The bidomain boundary element method provides a computationally efficient and accurate solution for neuron modeling.
• Facilitates realistic neural network simulations with complex neuron morphologies for advanced research.
• Advances the development of fast bidomain solvers for scalable neural simulations in brain modulation applications.