A boundary element method of bidomain modeling for predicting cellular responses to electromagnetic fields

Area of Science:

• Computational neuroscience
• Biophysics
• Electromagnetic modeling

Background:

• Existing cable equation models for electromagnetic field effects on neurons have limitations due to simplifying assumptions.
• Bidomain finite element methods offer more realistic neuron modeling by fully coupling cells and electric fields.

Purpose of the Study:

• To introduce a novel bidomain integral equation formulation for simulating electromagnetic coupling between stimulation devices and neuronal regions.
• To develop a boundary element method for accurate neuron modeling in brain stimulation scenarios.

Main Methods:

• Developed a boundary element formulation to solve an integral equation for electric and magnetic fields.
• Employed first-order nodal elements and a Crank-Nicholson time-stepping scheme.
• Validated the method using simulations of Hodgkin-Huxley axons and spherical cells.

Main Results:

• The boundary element method accurately simulates both electric and magnetic stimulation effects.
• This method eliminates the need for multi-scale volume meshes required by finite element methods.
• Simplifies modeling of complex cell populations and allows flexible device placement without re-meshing.

Conclusions:

• The bidomain boundary element method provides an accurate and computationally efficient alternative for simulating electromagnetic interactions with neurons.
• This approach facilitates multi-cell studies and neural network simulations with realistic neuron morphologies for research and therapeutic applications.